Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Effect of spatial variation characteristics on contouring of design storm depth

The ordinary kriging method, a geostatistical interpolation technique, was applied for developing contour maps of design storm depth in northern Taiwan using intensity–duration–frequency (IDF) data. Results of variogram modelling on design storm depths indicate that the design storms can be categorized into two distinct storm types: (i) storms of short duration and high spatial variation and (ii) storms of long duration and less spatial variation. For storms of the first category, the influence range of rainfall depth decreases when the recurrence interval increases, owing to the increasing degree of their spatial independence. However, for storms of the second category, the influence range of rainfall depth does not change significantly and has an average of approximately 72 km. For very extreme events, such as events of short duration and long recurrence interval, we do not recommend usage of the established design storm contours, because most of the interstation distances exceed the influence ranges. Our study concludes that the influence range of the design storm depth is dependent on the design duration and recurrence interval and is a key factor in developing design storm contours. Copyright © 2003 John Wiley & Sons, Ltd.     

Exposure dating and reinterpretation of coarse debris accumulations (‘rock glaciers’) in the Cairngorm Mountains,Scotland

Relict rock glaciers have considerable potential for contributing to palaeoclimatic reconstruction, but this potential is often undermined by lack of dating control and problems of interpretation. Here we reinvestigate and date four proposed ‘rock glaciers’ in the Cairngorm Mountains and show that the morphology of only one of these appears consistent with that of a true rock glacier produced by creep of underlying ice or ice‐rich sediment. All four features comprise rockslide or rock avalanche runout debris, and the possibility that all four represent unmodified runout accumulations cannot be discounted. Surface exposure dating of the four debris accumulations using cosmogenic ¹⁰Be produced uncertainty‐weighted mean ages of 15.4 ± 0.8 ka, 16.2 ± 1.0 ka, 12.1 ± 0.6 ka and 12.7 ± 0.8 ka. All four ages imply emplacement under cold stadial conditions, two prior to the Windermere Interstade of ca. 14.5–12.9 cal. ka BP and two during the Loch Lomond Stade of ca. 12.9–11.5 cal. ka BP. The above ages indicate that paraglacial rock‐slope failure on granite rockwalls occurred within a few millennia after deglaciation. The mean exposure ages obtained for runout debris at two sites – Strath Nethy (16.2 ± 1.0 ka) and Lairig Ghru (15.4 ± 0.8 ka) – are consistent with basal radiocarbon ages from Loch Etteridge, 22 km to the southwest (mean = 15.6 ± 0.3 cal. ka BP) and imply widespread deglaciation of the Cairngorms and adjacent valleys before 15 ka and possibly 16 ka. Copyright © 2008 John Wiley & Sons, Ltd.     

Uncoupling of coupled flows in soil—a finite element method

Coupled flow of water, chemicals, heat and electrical potential in soil are of significance in a variety of circumstances. The problem is characterized by the coupling between different flows, i.e. a flow of one type driven by gradients of other types, and by the dual nature of certain flows, i.e. combined convection and conduction. Effective numerical solutions to the problem are challenged due to the coupling and the dual nature. In this paper, we first present a general expression that can be used to represent various types of coupled flows in soil. A finite element method is then proposed to solve the generalized coupled flows of convection-conduction pattern. The unknown vector is first decomposed into two parts, a convective part forming a hyperbolic system and a conductive part forming a parabolic system. At each time step, the hyperbolic system is solved analytically to give an initial solution. To solve the multi-dimensional hyperbolic system, we assume that a common eigenspace exists for the coefficient matrices, so that the system can be uncoupled by transforming the unknown vector to the common eigenspace. The uncoupled system is solved by the method of characteristics. Using the solution of the hyperbolic system as the initial condition, we then solve the parabolic system by a Galerkin finite element method for space discretization and a finite difference scheme for time stepping. The proposed technique can be used for solving multi-dimensional, transient, coupled or simultaneous flows of convection-conduction type. Application to a flow example shows that the technique indeed exhibits optimality in convergence and in stability.     

树木年轮重建阿勒泰西部1481-2004年6-9月降水量序列

根据采自阿勒泰西部地区5个采点的树木年轮样本, 建立了该地区的树轮年表. 通过相关普查发现, 其差值年表序列与该地区当年6-9月的降水量存在明显的负相关关系, 且具有明确的树木生理学意义. 用沙勒哈(t, t 1)两个树轮差值年表序列可较好地重建该地区1481-2004年524 a来的当年6-9月的降水量, 且经过交叉检验表明重建方程稳定可靠. 分析发现, 阿勒泰西部地区524 a来的重建降水量序列具有9个偏干阶段和9个偏湿阶段, 并且具有34.8 a、 5.0～5.1 a、 4.3～4.4 a、 3.9～4.0 a、 3.7 a和3.4 a的显著干湿变化准周期. 重建降水量序列在1518年、 1548年、 1634年、 1765年、 1856年前后发生过突变.     

海上矿石运输最优船队的组建

针对海上矿石运输，介绍了一种完成定量任务并能使投资费用最少为目标的最优船队组建方法，同时对易变参数进行了灵敏度分析，使决策更加科学化     

A new type of Shipboard Meteorological Satellite Receiving-processing System

INTRODUCTIONTheSMSRPSisstalledontheshipnavigatingovertheseaandreceivessatellitecloudmapsatanytimetoprovidereliablereal-timedataofmeteorologyandocean.Itisanimportantequipmentforsafeguardofshipnavigation.Chinahasthousandsofoceanicships.Butalmostallships,withonlyveryfewexception,arenotstalledbySMSRPS.Thecausesareasfollows:1.Thetechniqueisverycomplicated.2.Thedevelopmentcostisveryhigh.3.Ifweintroduceforeignequipment,thecostistooexpensive.4.Foreignequipmentanditssoftwarearecompletelyclo…     

利用卫星遥感资料对南海北部陆架海洋表层温度锋的分析

利用7年(1993～1999)月平均的SST卫星遥感资料，分析了南海北部陆架区域海洋表层温度锋在一年中的逐月变化特征，表明南海北部陆架海洋表层温度锋存在明显的季节内变化。结合风场的卫星遥感资料，分析了东北季风对南海北部陆架温度锋的影响，表明东北季风风速的增加有利于温度锋强度的增强。通过对黑潮南海流套入侵较强的1999年2月与流套入侵较弱的1998年2月的SST卫星遥感资料的对比分析，考察了黑潮南海流套的入侵对南海北部陆架温度锋的影响，结果表明黑潮流套的较强入侵能够增加陆架温度锋的强度，对温度锋的走向也会产生一定的影响。     

全站仪似水准法和对向法高程测量的比较研究

通过对似水准法和对向法高程测量原理、误差来源及精度的比较分析,我们发现随着距离和竖直角度的增大,对向法高程测量中误差的变化大于似水准法高程测量中误差的变化;当两观测点间的水平距离小于等于1km时,对向法高程测量精度一般高于似水准法高程测量精度,但是当两观测点间的水平距离大于1km时,似水准法高程测量精度一般高于对向法高程测量精度。     

高精度重磁测量在商都地区CO2气藏勘查中的应用

通过在商都地区开展高精度重磁测量,提高了对盆地基底结构、构造、断裂特征与断裂分布以及汉诺坝组玄武岩分布的认识.商都坳陷具有单断式的结构特征,即西断东超.坳陷走向为北西向,由西北、东南两个次级洼槽组成;区内分布有两条北西南东向展布的深大断裂,为幔源CO2的运移通道,控制着CO2气藏的分布;局部强磁性、高密度异常为汉诺坝玄武岩的反映,是小丹岱组CO2气藏的良好盖层,控制着小丹岱组CO2气藏的分布.