A soil water and heat transfer model including changes in soil frost and thaw fronts

Freeze-thaw processes in soils,including changes in frost and thaw fronts(FTFs),are important physical processes.The movement of FTFs affects soil hydrothermal characteristics,as well as energy and water exchanges between the land surface and the atmosphere and hydrothermal processes in the land surface.This paper reduces the issue of soil freezing and thawing to a multiple moving-boundary problem and develops a soil water and heat transfer model which considers the effects of FTF on soil hydrothermal processes.A local adaptive variable-grid method is used to discretize the model.Sensitivity tests based on the hierarchical structure of the Community Land Model(CLM)show that multiple FTFs can be continuously tracked,which overcomes the difficulties of isotherms that cannot simultaneously simulate multiple FTFs in the same soil layer.The local adaptive variable-grid method is stable and offers computational efficiency several times greater than the high-resolution case.The simulated FTF depths,soil temperatures,and soil moisture values fit well with the observed data,which further demonstrates the potential application of this simulation to the land-surface process model.     

Role of heat and water transport in frost heaving of fine-grained porous media under negligible overburden pressure

An equation accurately describing the rate of frost heave is derived by using the mixture theory of continuum mechanics. It is shown that the rate of frost heave is determined mainly by the rate of heat removal and the rate of water intake. When the phase equilibrium holds in the system, the relation between the rate of heat removal and the rate of water intake is shown to depend mainly on the phase composition data of a given medium.By studying reported experimental data, it is found that the phase equilibrium may hold until th rate of heat removal reaches a certain critical value. When the rate of heat removal exceeds this critical value, the phase equilibrium may possibly be disrupted for some media.     

Coupling of a simultaneous heat and water model with a distributed hydrological model and evaluation of the combined model in a cold region watershed

Snow and frozen soil prevail in cold regions worldwide, and the integration of these processes is crucial in hydrological models. In this study, a combined model was developed by fully coupling a simultaneous heat and water model with a geomorphologically based distributed hydrological model. The combined model simulates vertical and lateral water transfer as well as vertical heat fluxes and is capable of representing the effects of frozen soil and snowmelt on hydrological processes in cold regions. This model was evaluated by using in situ observations in the Binggou watershed, an experimental watershed for cold region hydrology of the Watershed Allied Telemetry Experimental Research Project. Results showed that the model was able to predict soil freezing and thawing, unfrozen soil water content, and snow depth reasonably well. The simulated hydrograph was in good agreement with the in situ observation. The Nash–Sutcliffe coefficient of daily discharge was 0.744 for the entire simulation period, 0.472 from April to June, and 0.711 from June to November. This model can improve our understanding of hydrological processes in cold regions and assess the impacts of global warming on hydrological cycles and water resources. Copyright © 2012 John Wiley & Sons, Ltd.     

Transport of water in frozen soil IV. Analysis of experimental results on the effects of ice content

Effects of ice content on the transport of water in frozen soil are studied experimentally and theoretically under isothermal conditions. A physical law, that the flux of water in unsaturated frozen soil is proportional to the gradient of total water content is proposed. Theoretical justification is made by the use of the two-phase flow theory. The experimental results are shown to support the proposed physical law. The results of this study are presented in two parts and this is the second paper describing the theoretical aspects of the study.     

An improved process-based representation of stream solute transport in the soil and water assessment tools

Hydrological models have long been used to study the interactions between land, surface and groundwater systems, and to predict and manage water quantity and quality. The soil and water assessment tool (SWAT), a widely used hydrological model, can simulate various ecohydrological processes on land and subsequently route the water quality constituents through surface and subsurface waters. So far, in-stream solute transport algorithms of the SWAT model have only been minimally revised, even though it has been acknowledged that an improvement of in-stream process representation can contribute to better model performance with respect to water quality. In this study, we aim to incorporate a new and improved solute transport model into the SWAT model framework. The new process-based model was developed using in-stream process equations from two well established models—the One-dimensional Transport with Inflow and Storage model and the Enhanced Stream Water Quality Model. The modified SWAT model (Mir-SWAT) was tested for water quality predictions in a study watershed in Germany. Compared to the standard SWAT model, Mir-SWAT improved dissolved oxygen (DO) predictions by removing extreme low values of DO (<6 mg/L) simulated by SWAT. Phosphate concentration peaks were reduced during high flows and a better match of daily predicted and measured values was attained using the Mir-SWAT model (R² = 0.17, NSE = −0.65, RSR = 1.29 with SWAT; R² = 0.28, NSE = −0.04, RSR = 1.02 with Mir-SWAT). In addition, Mir-SWAT performed better than the SWAT model in terms of Chlorophyll-a content particularly during winter months, improving the NSE and RSR for monthly average Chl-a by 74 and 42%, respectively. With the new model improvements, we aim to increase confidence in the stream solute transport component of the model, improve the understanding of nutrient dynamics in the stream, and to extend the applicability of SWAT for reach-scale analysis and management.     

Sediment transport capacity and erosion processes: model concepts and reality

Sediment transport capacity, T_c, defined as the maximum amount of sediment that a flow can carry, is the basic concept in determining detachment and deposition processes in current process-based erosion models. Although defined conceptually and used extensively in modelling erosion, T_c was rarely measured. Recently, a series of laboratory studies designed to quantify effects of surface hydrologic conditions on erosion processes produced data sets feasible to evaluate the concept of T_c. A dual-box system, consisting of 1·8 m long sediment feeder box and a 5 m long test box, was used. Depending on the relative magnitudes of sediment delivery from feeder and test boxes, five scenarios are proposed ranging from deposition-dominated to transport-dominated sediment regimes. Results showed that at 5 per cent slope under seepage or 10 per cent slope under drainage conditions, the runoff from the feeder box caused in the additional sediment transport in the test box, indicating a transport-dominated sediment regime. At 5 per cent slope under drainage conditions, deposition occurred at low rainfall intensities. Increases in slope steepness, rainfall intensity and soil erodibility shifted the dominant erosion process from deposition to transport. Erosion process concepts from the Meyer–Wishmeier, Foster–Meyer and Rose models were compared with the experimental data, and the Rose model was found to best describe processes occurring during rain. A process-based erosion model needs to have components that can represent surface conditions and physical processes and their dynamic interactions. Copyright © 1999 John Wiley & Sons, Ltd.     

平扭耦联隔震体系的简化模型及有关参数变化的计算分析

本文根据隔震结构的动力学特点，提出了空间平扭耦联隔震结构简化计算模型，该模型仅有四个自由度，做定性分析和近似计算十分方便，力学概念明确可靠，文中并对这种体系的动力特性和反应，采用反应谱方法进行分析计算。还对平扭耦联隔震结构的重要特征参数对结构的动力特性和动力反应的影响，进行了较为全面的分析研究。     

Interbasin freshwater,heat and salt transport through the boundaries of the East and South China Seas from a variable-grid global ocean circulation model

It has long been recognized that the circulation in the East China Sea (ECS) and Japan/East Sea (JES) is closely related with that in Pacific, especially with the Kuroshio (e.g., Nitani[1], Hi-daka[2]). Based on current measurements in the Taiwan Strait a…     

Simulating daily soil water under foothills fescue grazing with the soil and water assessment tool model (Alberta,Canada)

Grazing is common in the foothills fescue grasslands and may influence the seasonal soil‐water patterns, which in turn determine range productivity. Hydrological modelling using the soil and water assessment tool (SWAT) is becoming widely adopted throughout North America especially for simulation of stream flow and runoff in small and large basins. Although applications of the SWAT model have been wide, little attention has been paid to the model's ability to simulate soil‐water patterns in small watersheds. Thus a daily profile of soil water was simulated with SWAT using data collected from the Stavely Range Sub‐station in the foothills of south‐western Alberta, Canada. Three small watersheds were established using a combination of natural and artificial barriers in 1996–97. The watersheds were subjected to no grazing (control), heavy grazing (2·4 animal unit months (AUM) per hectare) or very heavy grazing (4·8 AUM ha^?1). Soil‐water measurements were conducted at four slope positions within each watershed (upper, middle, lower and 5 m close to the collector drain), every 2 weeks annually from 1998 to 2000 using a downhole CPN 503 neutron moisture meter. Calibration of the model was conducted using 1998 soil‐water data and resulted in Nash–Sutcliffe coefficient (EF or R²) and regression coefficient of determination (r²) values of 0·77 and 0·85, respectively. Model graphical and statistical evaluation was conducted using the soil‐water data collected in 1999 and 2000. During the evaluation period, soil water was simulated reasonably with an overall EF of 0·70, r² of 0·72 and a root mean square error (RMSE) of 18·01. The model had a general tendency to overpredict soil water under relatively dry soil conditions, but to underpredict soil water under wet conditions. Sensitivity analysis indicated that absolute relative sensitivity indices of input parameters in soil‐water simulation were in the following order; available water capacity > bulk density > runoff curve number > fraction of field capacity (FFCB) > saturated hydraulic conductivity. Thus these data were critical inputs to ensure reasonable simulation of soil‐water patterns. Overall, the model performed satisfactorily in simulating soil‐water patterns in all three watersheds with a daily time‐step and indicates a great potential for monitoring soil‐water resources in small watersheds. Copyright © 2004 John Wiley & Sons, Ltd.     

Transport of water in frozen soil: V, Method for measuring the vapor diffusivity when ice is absent

A new experimental method is introduced for determining the relative magnitudes of liquid and vapor diffusion by using a small amount of soluble chemical as a tracer. The theoretical justification of the method is presented for the case where ice is absent. The feasibility of the method is demonstrated by an experiment using marine-deposited clay.     

Bulk density effects on soil hydrologic and thermal characteristics: A numerical investigation

Soil bulk density (ρ_b) is commonly treated as static in studies of land surface dynamics. Magnitudes of errors associated with this assumption are largely unknown. Our objectives were to (a) quantify ρ_b effects on soil hydrologic and thermal properties and (b) evaluate effects of ρ_b on surface energy balance and heat and water transfer. We evaluated 6 soil properties, volumetric heat capacity, thermal conductivity, soil thermal diffusivity, water retention characteristics, hydraulic conductivity, and vapour diffusivity, over a range of ρ_b, using a combination of 6 models. Thermal conductivity, water retention, hydraulic conductivity, and vapour diffusivity were most sensitive to ρ_b, each changing by fractions greater than the associated fractional changes in ρ_b. A 10% change in ρ_b led to 10–11% change in thermal conductivity, 6–11% change in saturated and residual water content, 49–54% change in saturated hydraulic conductivity, and 80% change in vapour diffusivity. Subsequently, 3 field seasons were simulated with a numerical model (HYDRUS‐1D) for a range of ρ_b values. When ρ_b increased 25% (from 1.2 to 1.5 Mg m^?3), soil temperature variation decreased by 2.1 °C in shallow layers and increased by 1 °C in subsurface layers. Surface water content differed by 0.02 m³ m^?3 for various ρ_b values during drying events but differences mostly disappeared in the subsurface. Matric potential varied by >100 m of water. Surface energy balance showed clear trends with ρ_b. Latent heat flux decreased 6%, sensible heat flux increased 9%, and magnitude of ground heat flux varied by 18% (with a 25% ρ_b increase). Transient ρ_b impacted surface conditions and fluxes, and clearly, it warrants consideration in field and modelling investigations.     

Effects of antecedent moisture and macroporosity on infiltration and water flow in frozen soil

Infiltration into frozen soil plays an important role in soil freeze–thaw and snowmelt-driven hydrological processes. To better understand the complex thermal energy and water transport mechanisms involved, the influence of antecedent moisture content and macroporosity on infiltration into frozen soil was investigated. Ponded infiltration experiments on frozen macroporous and non-macroporous soil columns revealed that dry macroporous soil produced infiltration rates reaching 10³ to 10⁴ mm day⁻¹, two to three orders of magnitude larger than dry non-macroporous soil. Results suggest that rapid infiltration and drainage were a result of preferential flow through initially air-filled macropores. Using recorded flow rates and measured macropore characteristics, calculations indicated that a combination of both saturated flow and unsaturated film flow likely occurred within macropores. Under wet conditions, regardless of the presence of macropores, infiltration was restricted by the slow thawing rate of pore ice, producing infiltration rates of 2.8 to 5.0 mm day⁻¹. Reduced preferential flow under wet conditions was attributed to a combination of soil swelling, due to smectite-rich clay (that reduced macropore volume), and pore ice blockage within macropores. In comparison, dry soil column experiments demonstrated that macropores provided conduits for water and thermal energy to bypass the frozen matrix during infiltration, reducing thaw rates compared with non-macroporous soils. Overall, results showed the dominant control of antecedent moisture content on the initiation, timing, and magnitude of infiltration and flow in frozen macroporous soils, as well as the important role of macropore connectivity. The study provides an important data set that can aid the development of hydrological models that consider the interacting effects of soil freeze–thaw and preferential flow on snowmelt partitioning in cold regions.     

基于双曲线模型的冻土中锚杆荷载传递特性研究

为研究不同冻结温度及含水率对冻结粉土中锚杆抗拔性能的影响,根据冻土与混凝土接触面的力学特性和变形规律,采用双曲线模型描述锚杆-冻土接触面的剪切特性.基于荷载传递法,建立考虑温度和含水率影响的锚杆荷载传递方程,采用有限差分法进行求解,得到锚固段的剪应力、轴力及承载力计算式;通过ABAQUS数值模拟验证荷载传递方程的合理性...     

The simulation of heat and water exchange in the boreal spruce forest by the land-surface model SWAP

All previous versions of a physically based land-surface model SWAP have assumed for simplicity that vegetation is fully covered by snow during the cold season. Such assumption is reasonable only for the regions dominated by short vegetation or for warm climates where snow processes are absent. The major goals of this paper are (i) modification of the latest version of SWAP by incorporation of tall vegetation into the cold-season parameterizations to make the model applicable for simulating heat and water transfer within a boreal forest biome and (ii) validation of the modified version using the data from a forested catchment located in the boreal zone. Modification of SWAP required to parameterize radiative and turbulent exchange between the forest crown and forest floor, partitioning snowfall between interception by the canopy (in doing so, snow interception differs from rain interception) and falling to the ground, formation of snow cover on the forest crown and forest floor including snow accumulation (both in solid and liquid fractions), snow evaporation, and snowmelt. The advanced model was validated using a set of hydrometeorological data measured during 18 years (1966–1983) at the Tayozhniy catchment (covered by boreal spruce forest), Valdai, Russia. Simulations of annual and monthly snow/rain interception, daily runoff at the catchment outlet, snow density, snow depth, snow water equivalent, soil water storage in three layers (0–20, 0–50 and 0–100 cm), and monthly evapotranspiration from the catchment were compared with observations. Analysis of the results of validation shows that the new version of the model SWAP reproduces the heat and water exchange processes occurring in mid-latitude boreal forest quite reasonable.     

A land surface soil moisture data assimilation framework in consideration of the model subgrid-scale heterogeneity and soil water thawing and freezing

The Ensemble Kalman Filter (EnKF) is well known and widely used in land data assimilation for its high precision and simple operation. The land surface models used as the forecast operator in a land data assimilation system are usually designed to consider the model subgrid-heterogeneity and soil water thawing and freezing. To neglect their effects could lead to some errors in soil moisture assimilation. The dual EnKF method is employed in soil moisture data assimilation to build a soil moisture data as- similation framework based on the NCAR Community Land Model version 2.0 (CLM 2.0) in considera- tion of the effects of the model subgrid-heterogeneity and soil water thawing and freezing: Liquid volumetric soil moisture content in a given fraction is assimilated through the state filter process, while solid volumetric soil moisture content in the same fraction and solid/liquid volumetric soil moisture in the other fractions are optimized by the parameter filter. Preliminary experiments show that this dual EnKF-based assimilation framework can assimilate soil moisture more effectively and precisely than the usual EnKF-based assimilation framework without considering the model subgrid-scale heteroge- neity and soil water thawing and freezing. With the improvement of soil moisture simulation, the soil temperature-simulated precision can be also improved to some extent.     

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.     

Validation of a soil water balance model using soil water content and pressure head data

The validation of soil water balance models and the evaluation of the quality of the model predictions at field‐scale require time‐series of in situ measured model outputs. In our study, we have validated such a model using a 6‐year period with time‐series of automatically recorded, daily volumetric soil water contents measured with the time‐domain reflectometry with intelligent microelements (TRIME) method and daily pressure heads measured with tensiometers. The comparisons of simulated with measured soil water contents and pressure heads were analysed using the modelling efficiency index (IA) and the square root of the mean square error (RMSE) in order to evaluate the prediction quality of the model. In our study, IA and RMSE, obtained either from the comparison of simulated with measured soil water contents or the comparison of calculated with observed pressure heads, in some cases lead to different results regarding the evaluation of the simulation quality of the soil water balance model. For example, a good fit between simulated and observed soil water contents does not necessarily result in a comparably good fit between the corresponding calculated and measured pressure heads. Therefore, a combined use of both measurement techniques, which takes into account their respective advantages and disadvantages, gives a more complete overview on the simulation quality of the soil water balance model than the single use of one of those techniques. Copyright © 2004 John Wiley & Sons, Ltd.     

A model of the redistribution of disaggregated soil particles by rainsplash

A model of the dispersion of splash droplets from a single raindrop impact on a sloping soil surface is combined with a theory of the entrainment of mineral particles from a disaggregated mixture in splash droplets to obtain a model of the dispersion of such particles by a raindrop impact. Stochastic modelling techniques extend this further to a model of the spatial redistribution of soil on a plot after a period of rainfall. Since the model is probabilistic and physically based it enables the incorporation of further advances in the understanding of splash erosion at all stages and can simulate the effect of the stochastic nature of rainfall and soil properties on the process. Several different situations are simulated. These include the movement of marked soil particles from point sources and the spatial patterns of erosion on a sloping plot. The model can also simulate the differential erosion of different soil particle size fractions.     

Development and testing of the Revised Morgan-Morgan-Finney (RMMF) soil erosion model under different pedological datasets

The study aims to investigate the effect of soil properties delineation on erosion modelling. To that end, the soil attributes of the Venetikos River catchment, northwestern Greece, are described using two pedological datasets, i.e. field samples and classification maps. The goal is to select the most appropriate for the accurate estimation of erosion. The Revised Morgan-Morgan-Finney (RMMF) model is developed per base map (annual or multi-annual), keeping all other parameters unchanged. Modelled sediment yield (SY) values are validated against “observed” ones, calculated utilizing the sediment rating curve methodology. Overall, the classification maps approach (164.35 t km^-2 year^?1) performed better than the soil samples one (82.97 t km^-2 year^?1), displaying higher convergence to the synthetic SY (548.9 t km^-2 year^?1). The discrepancy among approaches is attributed to the different computation methodologies (thus pedological background) used. Both approximations successfully identified the high-risk erosion areas. The same conclusions arose from the multi-annual application of the model.