东亚地区云微物理量分布特征的CloudSat卫星观测研究

本文利用2007~2010年整四年最新可利用的CloudSat卫星资料, 对东亚地区(15°~60°N, 70°~150°E)云的微物理量包括冰/液态水含量、冰/液态水路径、云滴数浓度和有效半径等的分布特征和季节变化进行了分析。本文将整个东亚地区划分为北方、南方、西北、青藏高原地区和东部海域五个子区域进行研究, 结果显示:东亚地区冰水路径值的范围基本在700 g m-2以下, 高值区分布在北纬40度以南区域, 在南方地区夏季的平均值最大, 为394.3 g m-2, 而在西北地区冬季的平均值最小, 为78.5 g m-2;而液态水路径的范围基本在600 g m-2以下, 冬季在东部海域的值最大, 达到300.8 g m-2, 夏季最大值为281.5 g m-2, 分布在南方地区上空。冰水含量的最高值为170 mg m-3, 发生在8 km附近, 南方地区夏季的值达到最大, 青藏高原地区的季节差异最大;而液态水含量在东亚地区的范围小于360 mg m-3, 垂直廓线从10 km向下基本呈现逐渐增大的趋势, 峰值位于1~2 km高度上。冰云云滴数浓度在东亚地区的范围在150 L-1以下, 水云云滴数浓度的值小于80 cm-3, 垂直廓线的峰值均在夏季最大。冰云有效半径在东亚地区的最大值为90 μm, 发生在5 km左右;水云有效半径在东亚地区的值分布在10 km以下, 最大值为10~12 μm, 基本位于1~2 km高度上。从概率分布函数来看, 东亚地区冰/水云云滴数浓度的分布呈现明显的双峰型, 其他量基本为单峰型。本文的结果可以为全球和区域气候模式在东亚地区对以上云微物理量的模拟提供一定的观测参考依据。     

Quasi‐static analysis of elastic behavior for some systems having higher fracture densities

Elastic behavior of geomechanical systems with interacting (but not intersecting) fractures is treated using generalizations of the Backus and the Schoenberg–Muir methods for analyzing layered systems whose layers are intrinsically anisotropic due to locally aligned fractures. By permitting the axis of symmetry of the locally anisotropic compliance matrix for individual layers to differ from that of the layering direction, we derive analytical formulas for interacting fractured regions with arbitrary orientations to each other. This procedure provides a systematic tool for studying how contiguous, but not yet intersecting, fractured domains interact, and provides a direct (though approximate) means of predicting when and how such interactions lead to more dramatic weakening effects and ultimately to failure of these complicated systems. The method permits decomposition of the system elastic behavior into specific eigenmodes that can all be analyzed, and provides a better understanding about which of these specific modes are expected to be most important to the evolving failure process. Copyright © 2009 John Wiley & Sons, Ltd.     

Effective stress concept in unsaturated soils: Clarification and validation of a unified framework

The effective stress principle, conventionally applied in saturated soils, is reviewed for constitutive modelling purposes. The assumptions for the applicability of Terzaghi's single effective stress are recalled and its advantages are inventoried. The possible stress frameworks applicable to unsaturated soil modelling are reassessed in a comparative manner, specifically the Bishop's single effective stress, the independent stress variables approach and the generalized stress framework. The latter considerations lead to the definition of a unified stress context, suitable for modelling soils under different saturation states. In order to qualify the implications brought by the proposed stress framework, several experimental data sets are re‐examined in the light of the generalized effective stress. The critical state lines (CSLs) at different saturation states tend to converge remarkably towards a unique saturated line in the deviatoric stress versus mean effective stress plane. The effective stress interpretation is also applied to isotropic paths and compared with conventional net stress conception. The accent is finally laid on a second key feature for constitutive frameworks based on a unified stress, namely the sufficiency of a unique mechanical yield surface besides the unique CSL. Copyright © 2007 John Wiley & Sons, Ltd.     

Matching hydrologic response to measured effective hydraulic conductivity

The objective of this study was to test the practicability of defining hydrologic response units as combinations of soil, land use and topography for modelling infiltration at the hillslope and catchment scales. In an experimental catchment in the East African Highlands (Kwalei, Tanzania), three methods of measuring infiltration were compared for their ability to capture the spatial variability of effective hydraulic conductivity: the constant head (CH) method; the tension infiltration (TI) method; and the mini‐rainfall simulation (RS) method. The three methods yielded different probability distributions of effective hydraulic conductivity and suggested different types of hydrologic response units. Independently from these measurements, the occurrence of infiltration‐excess overland flow was monitored over an area of 6 ha by means of overland flow detectors. The observed pattern of overland flow occurrence did not match any of the patterns suggested by the infiltration measurements. Instead, clusters of spots with overland flow were practically independent from field borders. Geostatistical analysis of the overland flow confirmed the absence of spatial correlation for distances over 40 m. The RS method yielded the pattern closest to the observations, probably because the method simulated better the processes that trigger infiltration‐excess overland flow, i.e. soil sealing and infiltration through macroporosity. The RS hydrologic response unit correlated significantly with observed overland flow frequency. However, the location of clusters and ‘hot spots’ of overland flow remained largely unexplained by land use, soil and topographic variables. It is concluded that using such landscape variables to define hydrologic units may create artificial boundaries that do no correspond to physical realities, especially if the stochastic component within hydrologic units is neglected. Copyright © 2005 John Wiley & Sons, Ltd.     

Gypsies in the palace: experimentalist's view on the use of 3‐D physics‐based simulation of hillslope hydrological response

As a fundamental unit of the landscape, hillslopes are studied for their retention and release of water and nutrients across a wide range of ecosystems. The understanding of these near‐surface processes is relevant to issues of runoff generation, groundwater–surface water interactions, catchment export of nutrients, dissolved organic carbon, contaminants (e.g. mercury) and ultimately surface water health. We develop a 3‐D physics‐based representation of the Panola Mountain Research Watershed experimental hillslope using the TOUGH2 sub‐surface flow and transport simulator. A recent investigation of sub‐surface flow within this experimental hillslope has generated important knowledge of threshold rainfall‐runoff response and its relation to patterns of transient water table development. This work has identified components of the 3‐D sub‐surface, such as bedrock topography, that contribute to changing connectivity in saturated zones and the generation of sub‐surface stormflow. Here, we test the ability of a 3‐D hillslope model (both calibrated and uncalibrated) to simulate forested hillslope rainfall‐runoff response and internal transient sub‐surface stormflow dynamics. We also provide a transparent illustration of physics‐based model development, issues of parameterization, examples of model rejection and usefulness of data types (e.g. runoff, mean soil moisture and transient water table depth) to the model enterprise. Our simulations show the inability of an uncalibrated model based on laboratory and field characterization of soil properties and topography to successfully simulate the integrated hydrological response or the distributed water table within the soil profile. Although not an uncommon result, the failure of the field‐based characterized model to represent system behaviour is an important challenge that continues to vex scientists at many scales. We focus our attention particularly on examining the influence of bedrock permeability, soil anisotropy and drainable porosity on the development of patterns of transient groundwater and sub‐surface flow. Internal dynamics of transient water table development prove to be essential in determining appropriate model parameterization. Copyright © 2010 John Wiley & Sons, Ltd.     

Finite element formulation for poroelastic problem with zero effective stress boundary condition

Zero effective stress boundary condition along with constant fluid flux is commonly encountered in geotechnical applications such as uncased borehole stability, fluid injection and production at an uncased borehole, hydraulic fracturing and sand production. This complex boundary condition introduces high nonlinearity in the numerical simulation. Conventional iterative methods such as Newton–Raphson method are required to solve this nonlinear problem iteratively, which involve huge computing time and also pose numerical difficulties on the convergence. To overcome this numerical difficulty and hence reduce the computing time, a novel numerical technique is proposed in this paper. Its performance is evaluated using a numerical example simulating fluid injection around an uncased borehole. Copyright © 2009 John Wiley & Sons, Ltd.     

Functional behaviour of pedotransfer functions in soil water flow simulation

Soil hydraulic parameter values for large‐scale modelling cannot be obtained by direct methods. Pedotransfer functions (PTFs) that relate soil hydraulic properties (SHPs) to generally available soil texture data may provide an alternative. A considerable number of PTF models has been developed, the application of three recent PTFs is evaluated. As a first step sets of SHPs derived from the PTFs are compared with measured sets of SHPs for three sites. No good agreement was found statistically between measured and PTF results or between PTF results. As a second step and from a practical point of view results for three hydrologically functional variables were compared and evaluated. The three selected functional variables are saturated hydraulic conductivity, k₀, in relation to infiltration excess runoff, available soil water amounts for evapotranspiration and water table depth for a specified upward flux or capillary rise. Derived k₀ distributions from PTFs show substantially less variance than from the measured data at all three sites. This can have a considerable impact on infiltration excess runoff, depending on the actual rainfall regime. Simulated available soil water amounts for evapotranspiration for some combinations of PTFs and sites are close to those obtained for measured SHPs, however, no consistency in results can be detected. Water table depths for specified upward flux densities using PTF derived SHPs are generally deeper than those based on measured SHPs and means a potentially higher water availability. Overall, differences in capillary rise among the selected PTFs and between measured and PTF based results are again inconsistent and show no clear relationship with soil texture. Finally, as a third step, effective SHPs were calculated by using spatially averaged texture as PTF input representing areal average behaviour. For these effective SHPs the calculated effective values for the three selected functional variables appear to be close to the areally averaged values obtained with step 2. The selected functional variables thus appear to depend linearly on the PTFs over the range for which the data are representative. This suggests that for our specific PTFs areal mean or effective values for the three functional variables can be obtained fairly accurately from a single measurement of a bulk collection of soil samples as input. Copyright © 2003 John Wiley & Sons, Ltd.