A New Method for the Determination of Area-Averaged Turbulent Surface Fluxes from Low-Level Flights Using Inverse Models

The low-level flight method (LLF) has been combined with linear inverse models (IM) resulting in an LLF+IM method for the determination of area-averaged turbulent surface fluxes. With this combination, the vertical divergences of the turbulent latent and sensible heat fluxes were calculated from horizontal flights. The statistical errors of the derived turbulent surface fluxes were significantly reduced. The LLF+IM method was tested both in numerical and field experiments. Large-eddy simulations (LES) were performed to compare ‘true’ flux profiles with ‘measurements’ of simulated flights in an idealised convective boundary layer. Small differences between the ‘true’ and the ‘measured’ fluxes were found, but the vertical flux divergences were correctly calculated by the LLF+IM method. The LLF+IM method was then applied to data collected during two flights with the Helipod, a turbulence probe carried by a helicopter, and with the research aircraft Do 128 in the LITFASS-98 field campaign. The derived surface fluxes were compared with results from eddy-covariance surface stations and with large-aperture scintillometer data. The comparison showed that the LLF+IM method worked well for the sensible heat flux at 77 and 200 m flight levels, and also for the latent heat flux at the lowest level. The model quality control indicated failures for the latent heat flux at the 200 m level (and higher), which were probably due to large moisture fluctuations that could not be modelled using linear assumptions. Finally the LLF+IM method was applied to more than twenty low-level flights from the LITFASS-2003 experiment. Comparison with aggregated surface flux data revealed good agreement for the sensible heat flux but larger discrepancies and a higher statistical uncertainty for the latent heat flux     

Temperature–Humidity Dissimilarity and Heat-to-water-vapour Transport Efficiency Above and Within a Pine Forest Canopy: the Role of the Bowen Ratio

Over the past 15 years atmospheric surface-layer experiments over heterogeneous canopies have shown that the vertical transfer of sensible heat and water vapour exhibit a strong dissimilarity. In particular, the sensible-heat-to-water-vapour transport efficiencies generally exceed unity. One of the main consequences is that evaporation (latent heat flux) computed by the flux-variance method is overestimated, as persistently demonstrated by comparisons with evaporation obtained with the eddy-correlation method. Various authors proposed to take into account the temperature–humidity dissimilarity to extend the applicability of the flux-variance method in order to compute evaporation from non-uniform surfaces. They attempted to connect the sensible-heat-to-water-vapour transport efficiency (λ) to the correlation coefficient between temperature and humidity turbulent fluctuations (R _Tq ). This approach was found to be successful over ‘wet’ surfaces for which λ can be approximated by R _Tq and ‘dry’ surfaces for which λ can be approximated by 1/R _Tq . However, no solution has been proposed until now for intermediate hydrological conditions. We investigated this question using eddy-correlation measurements above and inside a pine forest canopy. For both levels, our data present a strong likeness with previously published results over heterogeneous surfaces. In particular, they confirm that λ is R _Tq in wet conditions and 1/R _Tq in dry conditions. Moreover, we defined the range of the Bowen ratio (Bo) values for which those two approximations are valid (below 0.1 and greater than 1, respectively) and established a relationship between λ, R _Tq and Bo for the intermediate range of Bo. We are confident that this new parameterization will enlarge the applicability of the flux-variance method to all kinds of heterogeneous surfaces in various hydrological conditions     

Landscape Roughness Parameters for Sherwood Forest – Experimental Results

The aim of this work is to present experimentally evaluated effective roughnesses (zoe) of a partly forested landscape. Although the ratio of boundary-layer height to obstacle size was only of the order of 50, there still seemed to exist a height range of 75–200 m where surface-layer similarity was approximately valid. Attempts were made to use conventional wind profile analysis to evaluate zoe, but the small height range and the large number of variables initially led to unacceptable uncertainties. Fixing the displacement height zd, rather than fitting it, reduced the data scatter to an acceptable level. The profile-derived roughness lengths zop obtained in this way were in good agreement with previous work, and with an alternative roughness length estimate zof for which flux-derived profile parameters u* and * were used. This implies that the profile-derived roughnesses were consistent with the measured surface-layer momentum flux. Comparison of both roughness estimates also yielded an improved estimate of the displacement height. Besides this, the authors tested a landscape roughness evaluation method which makes use of the gustiness parameter Tu = u/U in the surface layer. The results obtained by this method were in fair agreement with the profile-derived data. In previous work, the gustiness method was advocated because it could be used at relatively low levels, perhaps even within the roughness sub-layer. At the present measuring site, this was not the case as the gustiness method was only valid in an approximate way, and for a limited height range.     

利用高阶交错网格有限差分法模拟地震波在非均匀孔隙介质中的传播

烃类储集层是一种复合多相介质，在固体颗粒的空隙中含有气体或液体. 研究弹性波在该类地层中的传播规律对于油气勘探开发，特别对于全波列声波测井有重要意义. 为了提高孔隙弹性介质数值模拟的计算效率，本文采用改进显式交错网格有限差分算法取代常用的空间域四阶和时间域二阶的速度 - 应力有限差分算法，算法的空间域为八阶、时间域为二阶. 虽然计算的时间步长略小于空间域四阶的情形，但高阶有限差分算法可以选择较粗糙的网格，因此补偿了计算的低效；同时高阶交错网格有限差分算法的空间频散性比低阶算法小. 利用该算法计算了一个两层模型的波场，同时还模拟了等效弹性和孔隙弹性模型中波的传播. 结果表明慢波及其影响明显，尽管慢波衰减很快，但被某一界面反射后，转换形成的P波和S波仍以正常的方式传播，且比慢波衰减小.     

A constitutive model for bonded geomaterials subject to mechanical and/or chemical degradation

The mechanical behaviour of bonded geomaterials is described by means of an elastoplastic strain‐hardening model. The internal variables, taking into account the ‘history’ of the material, depend on the plastic strains experienced and on a conveniently defined scalar measure of damage induced by weathering and/or chemical degradation. For the sake of simplicity, it is assumed that only internal variables are affected by mechanical and chemical history of the material. Despite this simplifying assumption, it can be shown that many interesting phenomena exhibited by weathered bonded geomaterials can be successfully described. For instance, (i) the transition from brittle to ductile behaviour with increasing pressure of a calcarenite with collapsing internal structure, (ii) the complex behaviour of chalk and other calcareous materials in oedometric tests, (iii) the chemically induced variation of the stress and strain state of such kind of materials, are all phenomena that can be qualitatively reproduced. Several comparisons with experimental data show that the model can capture the observed behaviour also quantitatively. Copyright © 2003 John Wiley & Sons, Ltd.     

Results from One-Year Continuous Operation of a Large Aperture Scintillometer over a Heterogeneous Land Surface

A large-aperture scintillometer (LAS) was operated continuouslyduring a period of more than one year over a heterogeneous land surface in Central Europeat the transition between marine and continental climates. The LAS measurements of the refractiveindex structure parameter, C_N ², were used to estimate the sensible heat flux. Thiswas possible for about 60to 80% of the time under daytime conditions during thesummer, with lower values obtained for the cold season (October to March). Using datafrom a three-week long field experiment, the LAS-based heat flux was compared with a weighedaverage of local heat flux measurements over the main land use classes (forest, agriculture,water) in the area, resulting in reasonable agreement. LAS-based heat fluxes were then used forcomparison with the heat flux values of a numerical weather prediction model. An over-predictionof the model heat flux was found in summer but the modelled values were lower than the LASderived data during the cold season.     

我国对流层二氧化碳非均匀动态分布特征及其成因

人类活动导致大气中温室气体浓度上升,是全球气候变暖的主要原因之一.本文针对已经连续运行13年的AIRS（Atmospheric Infrared Sounder）卫星反演的对流层中层CO₂浓度资料,利用地基观测结果对其进行验证,同时结合多种相关资料对我国区域CO₂浓度的时空分布及季节变化进行研究.结果显示：北半球30°N—60°N是CO₂浓度高值带,低值中心主要出现在15°S—15°N,140°W向东至100°E的低纬地区.地基观测与AIRS卫星反演结果基本一致,年增长率约为1.926 ppmv·a^-1.我国区域CO₂浓度空间分布上呈现北高南低的非均匀分布特征,4个高值中心分别位于东北地区西南部、内蒙古西部、新疆地区东部和西部,低值中心在云南和西藏地区.我国区域CO₂浓度有明显的季节变化特征,最高值出现在春季,冬夏季次之,秋季最低,其季节演变特点与风场的输送、降水量的清除和植被的吸收等密切相关.     

A depth‐integrated,coupled SPH model for flow‐like landslides and related phenomena

In the past decades, flow‐like catastrophic landslides caused many victims and important economic damage around the world. It is therefore important to predict their path, velocity and depth in order to provide adequate mitigation and protection measures. This paper presents a model that incorporates coupling between pore pressures and the solid skeleton inside the avalanching mass. A depth‐integrated, coupled, mathematical model is derived from the velocity–pressure version of the Biot–Zienkiewicz model, which is used in soil dynamics. The equations are complemented with simple rheological equations describing soil behaviour and are discretized using the SPH method. The accuracy of the model is assessed using a series of benchmarks, and then it is applied to back‐analyse the propagation stage of some catastrophic flow‐like slope movements for which field data are available. Copyright © 2008 John Wiley & Sons, Ltd.     

岩土材料塑性损伤模型及拱坝的变形局部化分析

常规的弹塑性模型难以模拟材料破坏时由于内部损伤累积引起的变形局部化带形成过程。将岩土材料抽象为含有孔洞的球形体胞单元的集合体,采用一个标量--孔洞体积百分比来刻划单元的损伤程度,基于细观损伤力学提出了一个适用于岩土材料的塑性损伤模型。将该模型通过用户子程序嵌入到大型有限元商业软件MRAC中,研究岩土材料的局部化剪切带形成过程。采用该模型,成功地对平面应力条件下岩土材料单轴受压和单轴受拉试件的局部化破坏分别进行了数值模拟。采用该模型,对拱坝在超水载作用下的变形局部化进行了研究,数值计算结果表明,上游坝踵处最先发生变形局部化现象。     

岩土类材料松弛模量的确定

岩土材料,即,土、冻土、岩石以及混凝土具有明显的徐变和应力松弛现象。通过徐变试验可容易确定徐变柔量;而应力松弛试验难度很大,故难以通过试验确定其松弛模量,以致岩土材料的松弛模量一直没有合适方法来确定。作者直接根据岩土材料的徐变柔量,通过解析计算给出了其松弛模量的显式算式,并且给出了三个解析算例。