Estimation of potential evapotranspiration over the Yellow River basin: reference crop evaporation or Shuttleworth–Wallace?

Potential evapotranspiration (PET) is a key input to hydrological models. Its estimation has often been via the Penman–Monteith (P–M) equation, most recently in the form of an estimate of reference evapotranspiration (RET) as recommended by FAO‐56. In this paper the Shuttleworth–Wallace (S–W) model is implemented to estimate PET directly in a form that recognizes vegetation diversity and temporal change without reference to experimental measurements and without calibration. The threshold values of vegetation parameters are drawn from the literature based on the International Geosphere–Biosphere Programme land cover classification. The spatial and temporal variation of the LAI of vegetation is derived from the composite NOAA‐AVHRR normalized difference vegetation index (NDVI) using a method based on the SiB2 model, and the Climate Research Unit database is used to provide the required meteorological data. All these data inputs are publicly and globally available. Consequently, the implementation of the S–W model developed in this study is applicable at the global scale, an essential requirement if it is to be applied in data‐poor or ungauged large basins. A comparison is made between the FAO‐56 method and the S–W model when applied to the Yellow River basin for the whole of the last century. The resulting estimates of RET and PET and their association with vegetation types and leaf area index (LAI) are examined over the whole basin both annual and monthly and at six specific points. The effect of NDVI on the PET estimate is further evaluated by replacing the monthly NDVI product with the 10‐day product. Multiple regression relationships between monthly PET, RET, LAI, and climatic variables are explored for categories of vegetation types. The estimated RET is a good climatic index that adequately reflects the temporal change and spatial distribution of climate over the basin, but the PET estimated using the S–W model not only reflects the changes in climate, but also the vegetation distribution and the development of vegetation in response to climate. Although good statistical relationships can be established between PET, RET and/or climatic variables, applying these relationships likely will result in large errors because of the strong non‐linearity and scatter between the PET and the LAI of vegetation. It is concluded that use of the implementation of the S–W model described in this study results in a physically sound estimate of PET that accounts for changing land surface conditions. Copyright © 2006 John Wiley & Sons, Ltd.     

First‐year post‐fire erosion rates in Bitterroot National Forest,Montana

Accelerated runoff and erosion commonly occur following forest fires due to combustion of protective forest floor material, which results in bare soil being exposed to overland flow and raindrop impact, as well as water repellent soil conditions. After the 2000 Valley Complex Fires in the Bitterroot National Forest of west‐central Montana, four sets of six hillslope plots were established to measure first‐year post‐wildfire erosion rates on steep slopes (greater than 50%) that had burned with high severity. Silt fences were installed at the base of each plot to trap eroded sediment from a contributing area of 100 m². Rain gauges were installed to correlate rain event characteristics to the event sediment yield. After each sediment‐producing rain event, the collected sediment was removed from the silt fence and weighed on site, and a sub‐sample taken to determine dry weight, particle size distribution, organic matter content, and nutrient content of the eroded material. Rainfall intensity was the only significant factor in determining post‐fire erosion rates from individual storm events. Short duration, high intensity thunderstorms with a maximum 10‐min rainfall intensity of 75 mm h^?1 caused the highest erosion rates (greater than 20 t ha^?1). Long duration, low intensity rains produced little erosion (less than 0·01 t ha^?1). Total C and N in the collected sediment varied directly with the organic matter; because the collected sediment was mostly mineral soil, the C and N content was small. Minimal amounts of Mg, Ca, and K were detected in the eroded sediments. The mean annual erosion rate predicted by Disturbed WEPP (Water Erosion Prediction Project) was 15% less than the mean annual erosion rate measured, which is within the accuracy range of the model. Published in 2007 by John Wiley & Sons, Ltd.     

Non‐local damage modelling of concrete: a procedure for the determination of model parameters

This paper presents a procedure for the determination of parameters of non‐local damage models. This is to assure a consistent response of a non‐local damage model, as choice of the internal length and other parameters of the model are varied. Correlations between the internal length and other parameters governing the local constitutive behaviour of the model are addressed and exploited. Focus is put on the relationship between the internal length of the non‐local model and the width of the fracture process zone. Numerical examples are used to demonstrate the rigour of the proposed method. Copyright © 2006 John Wiley & Sons, Ltd.     

Identifying parameters controlling soil delayed behaviour from laboratory and in situ pressuremeter testing

The aim of this paper is to present a methodology for identifying the soil parameters controlling the delayed behaviour from laboratory and in situ pressuremeter tests by using an elasto‐viscoplastic model (EVP‐MCC) based on Perzyna's overstress theory and on the elasto‐plastic Modified Cam Clay model. The influence of both the model parameters and the soil permeability was studied under the loading condition of pressuremeter tests by coupling the proposed model equations with Biot's consolidation theory. On the basis of the parametric study, a methodology for identifying model parameters and soil permeability by inverse analysis from three levels of constant strain rate pressuremeter tests was then proposed and applied on tests performed on natural Saint‐Herblain clay. The methodology was validated by comparing the optimized values of soil parameters and the values of the same parameters obtained from laboratory test results, and also by using the identified parameters to simulate other tests on the same samples. The analysis of the drainage condition and the strain rate effect during a pressuremeter test demonstrated the coupled influence of consolidation and viscous effects on the test results. The numerical results also showed that the inverse analysis procedure could successfully determine the parameters controlling the time‐dependent soil behaviour. Copyright © 2008 John Wiley & Sons, Ltd.     

Evaluating the stress intensity factors of anisotropic bimaterials using boundary element method

This paper presents a boundary element method (BEM) procedure for a linear elastic fracture mechanics analysis in two‐dimensional anisotropic bimaterials. In this formulation, a displacement integral equation is only collocated on the uncracked boundary, and a traction integral equation is only collocated on one side of the crack surface. A fundamental solution (Green's function) for anisotropic bimaterials is also derived and implemented into the boundary integral formulation so that except for the interfacial crack part, the discretization along the interface can be avoided. A special crack‐tip element is introduced to capture the exact crack‐tip behavior. A computer program using FORTRAN has been developed to effectively calculate the stress intensity factors of an anisotropic bimaterial. This BEM program has been verified to have a good accuracy with previous studies. In addition, a central cracked bimaterial Brazilian specimen constituting cement and gypsum is prepared to conduct the Brazilian test under diametral loading. The result shows that the numerical analysis can predict relatively well the direction of crack initiation and the path of crack propagation. Copyright © 2007 John Wiley & Sons, Ltd.     

方格星虫多糖对小鼠免疫活性的影响

研究了方格星虫多糖对小鼠免疫功能的影响。结果表明:方格星虫粗多糖(CSG)能够对抗环磷酰胺引起的小鼠胸腺萎缩和脾萎缩,明显拮抗白细胞减少,对正常小鼠免疫器官没有明显的增重作用;对正常小鼠腹腔巨噬细胞吞噬中性红能力有显著的增强作用。方格星虫多糖粗品(CSG)和精制多糖(SG1和SG2)能够促进小鼠脾淋巴细胞增殖,且能与ConA协同作用。     

大姚、民乐-山丹地震前的中短期前兆

利用川滇及青藏块体北缘地区上世纪80年代至2003年的跨断层流动形变观测资料，借助应变强度比时序曲线，通过丽江Ms7．0、景泰Ms5．9等震例研究，获取与2003年7月大姚Ms6．2、Ms6．1地震和10月民乐-山丹Ms6．1地震孕育有关的中短期前兆异常。分析表明：①大姚、民乐-山丹地震前1年内异常场地数远超过场地总数的1／3，相对近源区异常在发震地区占一半；②大姚Ms6．2地震前半年内，盐源-丽江地块及其邻近区域以及北东向延伸的川滇交界地区异常明显增多，相对集中。2003年7～9月大幅减弱；而鲜水河断裂、楚雄-通海断裂南段有所增强；③民乐-山丹震前十余天至数月内相对近源区及其邻近构造区域呈现异常相对分散-相对集中且局部增强-一临震前减弱的时空演变特征。震后一个月内青藏块体北缘地区仍有近1／3的场地呈现异常(海原一六盘山断裂带尤为突出)。     

Homogenization‐based analysis of anisotropic damage in brittle materials with unilateral effect and interactions between microcracks

This paper is devoted to micromechanical modeling of induced anisotropic damage in brittle geomaterials. The formulation of the model is based on a proper homogenization procedure by taking into account unilateral effects and interactions between microcracks. The homogenization procedure is developed in the framework of Eshelby's inclusion solution and Ponte‐Castaneda and Willis (J. Mech. Phys. Solids 1995; 43 :1919–1951) estimate. The homogenization technique is combined with the thermodynamics framework at microscopic level for the determination of damage evolution law. A rigorous crack opening–closure transition condition is established and an energy‐release‐rate‐based damage criterion is proposed. Computational aspects on the implementation of micromechanical model are also discussed. The proposed model is evaluated by comparing numerical predictions with experimental data for various laboratory tests on concrete. Parametric studies on unilateral effects and influences of microcracks interactions are finally performed and analyzed. Copyright © 2008 John Wiley & Sons, Ltd.     

Discrete modelling of time‐dependent rockfill behaviour

This modelling study deals with the time‐dependent behaviour of rockfill media, which is of particular interest during the life of rockfill dams. Breakage of rock blocks and crack propagation are the main processes responsible for rockfill creep and collapse. The modelling procedure presented here is performed on two scales: on the rock block scale, where the grain is taken to be an assembly of rigid particles initially endowed with cohesive bonds, and on the rockfill scale, which is taken to involve a set of breakable grains interacting via contact and friction processes. The grain breakage process is described in term of a thermodynamically consistent damage interface model, where the damage is a gradual delayed process. This model was implemented in a non‐smooth contact dynamics code. The effects of the main parameters involved were analysed by performing numerical studies. The ability of the model to predict the creep behaviour of rockfill media is confirmed by presenting several simulations. Copyright © 2008 John Wiley & Sons, Ltd.     

地下水对混凝土的pHs分解性侵蚀判别式

针对地下水对混凝土分解性侵蚀指标pHs现行评价方法中存在水样间pHs横向没有可对比性和难以定量绘制pHs分解性侵蚀分区图等两大缺陷,推出了一个新的pHs分解性侵蚀判别式sk1,用编制的"地下水对混凝土分解性侵蚀判别式计算程序",完成了新定义的判别式sk1的计算,计算结果不但具备了水样间的每种分解性侵蚀横向对比的功能,而且还具备了定量绘制每种分解性侵蚀分区图和定量绘制3种分解性侵蚀综合评价迭加图的功能。通过实例应用,形象地印证了现行评价方法存在的缺陷及现推出的判别式、定量编绘侵蚀分区图的可行性。