A Modified Soil Moisture Model for Two‐Layer Soil

There can be marked variations in soil hydraulic properties in a soil vertical profile from the soil surface to the base of the root zone. Many existing two‐layer soil moisture (TLSM) models cannot well describe typical stratified soil profiles. A modified TLSM model is presented in this study. The modified model results and those from two existing models are compared with field observations. The modified TLSM model had the best agreement with the field observations. In both the surface layer and the root zone layer, the root mean square errors of soil moisture estimated by the modified model were smaller than those for the other models. The parameters in the modified TLSM model are relatively easy to determine. The modified TLSM model offers clear advantages over current TLSM models.     

Implications of measurement metrics on soil freezing curves: A simulation of freeze–thaw hysteresis

Soil freeze–thaw events have important implications for water resources, flood risk, land productivity, and climate change. A property of these phenomena is the relationship between unfrozen water content and sub-freezing temperature, known as the soil freezing characteristic curve (SFC). It is documented that this relationship exhibits hysteretic behaviour when frozen soil thaws, leading to the definition of the soil thawing characteristic curve (STC). Although explanations have been given for SFC/STC hysteresis, the effect that ‘scale’ – particularly ‘measurement scale’ – may have on these curves has received little attention. The most commonly used measurement scale metric is the ‘support’, which is the spatial (or temporal) unit within which the measured variable is integrated or soil volume sampled. We show (a) measurement support can influence the range and shape of the SFC and (b) hysteresis can be attributed, in part, to the support and location of the measurements comprising the SFC/STC. We simulated lab measured temperature, volumetric water content (VWC), and permittivity from soil samples undergoing freeze–thaw transitions using Hydrus-1D and a modified Dobson permittivity model. To assess the effect of measurement support and location on SFC/STC, we masked the simulated temperature and VWC/permittivity extent to match the instrument's support and location. By creating a detailed simulation of the intra- and inter-support variability associated with the penetration of a freezing front, we demonstrate how measurement support and location can influence the temperature range over which water freezing events are captured. We show it is possible to simulate hysteresis in homogenous media with purely geometric considerations, suggesting that SFC/STC hysteresis may be more of an apparent phenomenon than mechanistically real. Lastly, we develop an understanding of how the location and support of soil temperature and VWC/permittivity measurements influence the temperature range over which water freezing events are captured.     

Upscaling of field-scale soil moisture measurements using distributed land surface modeling

Accurate coarse-scale soil moisture information is required for robust validation of current- and next-generation soil moisture products derived from spaceborne radiometers. Due to large amounts of land surface and rainfall heterogeneity, such information is difficult to obtain from existing ground-based networks of soil moisture sensors. Using ground-based field data collected during the Soil Moisture Experiment in 2002 (SMEX02), the potential for using distributed modeling predictions of the land surface as an upscaling tool for field-scale soil moisture observations is examined. Results demonstrate that distributed models are capable of accurately capturing a significant level of field-scale soil moisture heterogeneity observed during SMEX02. A simple soil moisture upscaling strategy based on the merger of ground-based observations with modeling predictions is developed and shown to be more robust during SMEX02 than upscaling approaches that utilize either field-scale ground observations or model predictions in isolation.     

Stress-dissipative heat geotemperature precursor mechanism of earthquakes

ＩｎｔｒｏｄｕｃｔｉｏｎＭａｎｙｅａｒｔｈｑｕａｋｅｓｈａｖｅｐｒｏｖｅｎｔｈａｔｔｈｅｒｅｅｘｉｓｔｅｄｓｈａｌｏｗｇｅｏｔｅｍｐｅｒａｔｕｒｅａｎｏｍａｌｉｅｓａｒｏｕｎｄｔｈｅｅａｒｔｈｑｕａｋｅｒｅｇｉｏｎｓｂｅｆｏｒｅｔｈｅｅｖｅｎｔｓｏｃ...     

Multifield coupling model and its applications for pile foundation in permafrost

In the construction of Qinghai-Tibet railway,to avoid diseases caused by frost heave and thaw col-lapse of frozen ground,besides the normal bridges over the rivers,a lot of dry bridge structures have been built to replace subgrade in the regions of high tem-perature and high ice content frozen soil.So,the problems on forming mechanism of bearing capacity of pile foundation in cold regions already become one of hot spot problems in frozen soil engineering.Freezing force and frost heave force ar…     

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.     

Spatial scale effect on seasonal streamflows in permafrost catchments on the Qinghai–Tibet Plateau

The scale issue is of central concern in hydrological processes to understand the potential upscaling or downscaling methodologies, and to develop models for scaling the dominant processes at different scales and in different environments. In this study, a typical permafrost watershed in the Qinghai‐Tibet Plateau was selected. Its hydrological processes were monitored for 4 years from 2004 to 2008, measuring the effects of freezing and thawing depth of active soil layers on runoff processes. To identify the nature and cause of variation in the runoff response in different size catchments, catchments ranging from 1·07 to 112 km² were identified in the watershed. The results indicated that the variation of runoff coefficients showed a ‘V’ shape with increasing catchment size during the spring and autumn seasons, when the active soil was subjected to thawing or freezing processes. A two‐stage method was proposed to create runoff scaling models to indicate the effects of scale on runoff processes. In summer, the scaling transition model followed an exponential function for mean daily discharge, whereas the scaling model for flood flow exhibited a linear function. In autumn, the runoff process transition across multiple scales followed an exponential function with air temperature as the driving factor. These scaling models demonstrate relatively high simulation efficiency and precision, and provide a practical way for upscaling or downscaling runoff processes in a medium‐size permafrost watershed. For permafrost catchments of this scale, the results show that the synergistic effect of scale and vegetation cover is an important driving factor in the runoff response. Copyright © 2011 John Wiley & Sons, Ltd.     

Muskingum equation based downstream sediment flow simulation models for a river system

Applications of sediment transport and water flow characteristics based sediment transport simulation models for a river system are presented in this study. An existing water–sediment model and a new sediment–water model are used to formulate the simulation models representing water and sediment movement in a river system. The sediment–water model parameters account for water flow characteristics embodying sediment transport properties of a section. The models are revised formulations of the multiple water inflows model describing water movement through a river system as given by the Muskingum principle. The models are applied to a river system in Mississippi River basin to estimate downstream sediment concentration, sediment discharge, and water discharge. River system and the river section parameters are estimated using a revised and the original multiple water inflows models by applying the genetic algorithm. The models estimate downstream sediment transport rates on the basis of upstream sediment/water flow rates to a system. Model performance is evaluated by using standard statistical criteria;downstream water discharge resulting from the original multiple water inflows model using the estimated river system parameters indicate that the revised models satisfactorily describe water movement through a river system. Results obtained in the study demonstrate the applicability of the sediment transport and water flow characteristics-based simulation models in predicting downstream sediment transport and water flow rates in a river system.     

季冻区砂性土堤防水分场温度场及水分迁移试验研究

以黑龙江干流堤防工程实际环境为研究基础,依托水分迁移试验装置,测试了干流堤防典型砂性土试样在冻融循环下的温度场、水分场、应力场的分布情况.结果表明:堤顶混凝土公路破坏与堤身不均匀沉降有关,温度变化引起堤基含水率出现梯度变化,从而出现应力场变化,且温度梯度含水率梯度呈线性关系.地基稳定冻结深度达到1.12 m.结合实测数...     

Grid-size dependence of Cauchy boundary conditions used to simulate stream–aquifer interactions

This work examines the simulation of stream–aquifer interactions as grids are refined vertically and horizontally and suggests that traditional methods for calculating conductance can produce inappropriate values when the grid size is changed. Instead, different grid resolutions require different estimated values. Grid refinement strategies considered include global refinement of the entire model and local refinement of part of the stream. Three methods of calculating the conductance of the Cauchy boundary conditions are investigated. Single- and multi-layer models with narrow and wide streams produced stream leakages that differ by as much as 122% as the grid is refined. Similar results occur for globally and locally refined grids, but the latter required as little as one-quarter the computer execution time and memory and thus are useful for addressing some scale issues of stream–aquifer interactions. Results suggest that existing grid-size criteria for simulating stream–aquifer interactions are useful for one-layer models, but inadequate for three-dimensional models. The grid dependence of the conductance terms suggests that values for refined models using, for example, finite difference or finite-element methods, cannot be determined from previous coarse-grid models or field measurements. Our examples demonstrate the need for a method of obtaining conductances that can be translated to different grid resolutions and provide definitive test cases for investigating alternative conductance formulations.     

MODELING SEDIMENT PROCESSES ON SMALL WATERSHEDS:A CONCEPTUAL FRAMEWORK Ⅰ. BROAD SHALLOW FLOW PROCESSE

A conceptual modeling framework for developing process-based mathematical models of sediment generation, transport, and deposition on broad shallow flow areas is presented. The governing equations relevant to process-based modeling of broad shallow flow sediment processes on small watersheds are presented and the effects of space and time averaging on the predictive equations are described. Starting from the most general one-dimensional, unsteady model of sediment processes, simpler model structures are obtained and the successive simplifications made on the governing equations in order to obtain simpler and less detailed formulations are described. Specific model formulations are given for illustrative purposes and applications of these models to erosion and sediment yield prediction from broad shallow flow areas are shown using sediment data from rainfall simulator plots. In spite of some progress made in the development of process-based erosion and sediment yield models from broad shallow flow areas, further developmental modeling efforts must be based on a clear separation between hydrologic and hydraulic processes, and the soil properties which are significant for each.     

A comparison of non-linear least square and GLUE for model calibration and uncertainty estimation for pesticide transport in soils

The problems of calibrating soil hydraulic and transport parameters are well documented, particularly when data are limited. Programs such as CXTFIT, UUCODE and PEST, based on well established principles of statistical inference, will often provide good fits to limited observations giving the impression that a useful model of a particular soil system has been obtained. This may be the case, but such an approach may grossly underestimate the uncertainties associated with future predictions of the system and resulting dependent variables. In this paper, this is illustrated by an application of CXTFIT within the generalised likelihood uncertainty estimation (GLUE) approach to model calibration which is based on a quite different philosophy. CXTFIT gives very good fits to the observed breakthrough curves for several different model formulations, resulting in very small parameter uncertainty estimates. The application of GLUE, however, shows that much wider ranges of parameter values can provide acceptable fits to the data. The wider range of potential outcomes should be more robust in model prediction, especially when used to constrain field scale models.     

Using targeted short‐term field investigations to calibrate and evaluate the structure of a hydrological model

This study combines the application of a hydrological model with the use of field data derived from short period measurement campaigns at two sites, one a low topography forested area and the other a steep grassland catchment. The main objective was to determine if the structure of the widely used Pitman model could be considered appropriate for simulating the field data. The model is typically applied at coarse spatial and temporal (1 month) scales, while the tests reported here use data from small catchments and are applied in a daily version of the model. The results demonstrate the importance of ensuring that field observations are measuring the same hydrological variables as the model simulations. At one study site, there was a mismatch in the soil moisture data that was corrected by incorporating a two‐layer soil algorithm into the model. The model results from both field sites identified the sensitivity of the model to assumptions about evaporative demands and indicate that the model structure is very sensitive to the potential evaporation inputs. The overall conclusion is that the model structure is generally appropriate for simulating the hydrological responses at the two sites, but that there remain some unresolved uncertainties about specific model components and the use of certain types of input data. The study lends support for the future development of a more complete daily version of this widely used model. Copyright © 2013 John Wiley & Sons, Ltd.     

Seismic response of earth dams: some recent developments

The paper focuses on theoretical methods for estimating the dynamic response of earth dams to earthquake ground excitation. Following an outline of the historical developments in this field, the basic concepts/models for response analysis are introduced and their salient features, advantages and limitations are elucidated. The major phenomena associated with, and factors influencing, the response are identified and studied. Particular emphasis is accorded to inhomogeneity due to dependence of soil stiffness on confining pressure, nonrectangular canyon geometry, and nonlinear-inelastic soil behaviour. Several new formulations that have evolved over the last five years are outlined and characteristic results provide considerable insight into the problem. The simplicity of some of these formulations is underlined and attempts are made to compare their predictions with measurements from full-scale, natural and man-made, forced vibration tests. The basic validity as well as the limitations of the proposed analysis methods is demonstrated and topics of needed further research are suggested.     

On the prediction of the Toce alpine basin floods with distributed hydrologic models

With the objective of improving flood predictions, in recent years sophisticated continuous hydrologic models that include complex land‐surface sub‐models have been developed. This has produced a significant increase in parameterization; consequently, applications of distributed models to ungauged basins lacking specific data from field campaigns may become redundant. The objective of this paper is to produce a parsimonious and robust distributed hydrologic model for flood predictions in Italian alpine basins. Application is made to the Toce basin (area 1534 km²). The Toce basin was a case study of the RAPHAEL European Union research project, during which a comprehensive set of hydrologic, meteorological and physiographic data were collected, including the hydrologic analysis of the 1996–1997 period. Two major floods occurred during this period. We compare the FEST04 event model (which computes rainfall abstraction and antecedent soil moisture conditions through the simple Soil Conservation Service curve number method) and two continuous hydrologic models, SDM and TDM (which differ in soil water balance scheme, and base flow and runoff generation computations). The simple FEST04 event model demonstrated good performance in the prediction of the 1997 flood, but shows limits in the prediction of the long and moderate 1996 flood. More robust predictions are obtained with the parsimonious SDM continuous hydrologic model, which uses a simple one‐layer soil water balance model and an infiltration excess mechanism for runoff generation, and demonstrates good performance in both long‐term runoff modelling and flood predictions. Instead, the use of a more sophisticated continuous hydrologic model, the TDM, that simulates soil moisture dynamics in two layers of soil, and computes runoff and base flow using some TOPMODEL concepts, does not seem to be advantageous for this alpine basin. Copyright © 2006 John Wiley & Sons, Ltd.     

Prediction of sand transport rates by waves and currents in the coastal zone

The predictions of a sand transport research model and Bijker's (J. Waterways, Harbours Coastal Eng. Div. ASCE 97 (WW4) (1971) 687) engineering model are compared with data obtained in wave-current conditions at three field sites. A key element in the present study is that the bed roughness at the three sites has been estimated from predictions of the sand ripple dimensions. The comparisons between suspended sand concentrations and transport rates show that a considerable amount of uncertainty (factor ±5 or more) arises when individual predictions are compared with the measurements. However, the overall bias in each set of comparisons is smaller than this, with overall agreement being within a factor of ±2 in most cases. While the results demonstrate that research models, adapted for field application, may be used to make practical sand transport computations with as much accuracy as engineering formulations, the true benefit of research models lies in the improved understanding of transport processes that they provide. This is illustrated with reference to the mechanism of grain size sorting caused by oblique incidence of waves on a current.     

宁波深厚软土动力特性研究及其地震响应分析

利用Davidenkov模型对宁波深厚软土地区土动力试验数据进行拟合分析,确定相关拟合参数,结果表明该模型能较好地描述软土的动剪切模量和动阻尼比随动剪应变的变化规律。然后,基于等价线性化方法以及ANSYS软件的APDL语言进行简单的二次开发,编写Davidenkov模型计算程序,并利用经典的SHAKE软件进行验证,结果表明该程序计算结果很好,能解决目前商业软件中缺乏合适的土体非线性动力本构模型的不足。最后,利用开发的模型计算分析宁波深厚软土的地震响应规律,并与规范推荐的简化计算公式进行比较分析,结果表明规范推荐的公式得到的土体响应偏大,对于地下结构的抗震设计偏于保守。研究成果对确定深厚软土地区土体动力特性及其地震响应研究提供了合理的计算方法。     

A modified ISBA surface scheme for modeling the hydrology of Athabasca River Basin with GCM-scale data

A soil–vegetation–atmosphere transfer model (SVAT), interactions between the soil–biosphere–atmosphere (ISBA) of Météo France, is modified and applied to the Athabasca River Basin (ARB) to model its water and energy fluxes. Two meteorological datasets are used: the archived forecasts from the Meteorological Survey of Canada’s Global Environmental Multiscale Model (GEM) and the European Centre for Mid-range Weather Forecasts global re-analysis (ERA-40), representing spatial scales typical of a weather forecasting model and a global circulation model (GCM), respectively. The original treatment of soil moisture and rainfall in ISBA (OISBA) is modified to statistically account for sub-grid heterogeneity of soil moisture and rainfall to produce new, highly non-linear formulations for surface and sub-surface runoff (MISBA). These new formulations can be readily applied to most existing SVATs. Stand alone mode simulations using the GEM data demonstrate that MISBA significantly improves streamflow predictions despite requiring two fewer parameters than OISBA. Simulations using the ERA-40 data show that it is possible to reproduce the annual variation in monthly, mean annual, and annual minimum flows at GCM scales without using downscaling techniques. Finally, simulations using a simple downscaling scheme show that the better performance of higher resolution datasets can be primarily attributed to improved representation of local variation of land cover, topography, and climate.