一次西南低涡过程的地面场中尺度分析

本文运用带通滤波分析方法,对四川省1981年7月由西南低涡等天气系统引起的大暴雨过程作了地面气象要素场的中尺度诊断分析。结果表明,带通滤波后的场具有明显的中尺度特征,其中一些物理量场的适当配置与6小时后的暴雨中心有较好的对应关系。说明了此带通滤波方法在实际业务预报中具有一定的使用价值。     

长江中下游地区地下水卤素元素的形成及其分布规律

本文以丰富的实际资料，论证了地下水的卤素元素（F、Cl、Br、I）的形成、含量及其分布规律与含水介质成分、上覆岩土性质、地下水退流条件、氧化还原环境、地下水矿化度之间的关系。根据江汉平原东部区和鄱阳湖区地下水中Br、I元素的调查研究结果及它们形成的控制因素与分布规律，结合长江三角洲南部区水文地球化学环境条件分析对比，指出该区是一个形成Br、I矿泉水的有利地区。     

台湾海峡西岸重力场变化梯度场图像分析

讨论了应用曲面拟合技术和利用离散的重力观测点的重力变化对特定研究区的重力场变化进行数值模拟以获取该区重力变化空间梯度场图像的具体方法;并应用该方法获得了2004年3月福建安溪3.6级地震前的重力变化空间梯度场图像,发现空间梯度场图像能清晰地反映出地震前重力场群体性变化的差异性;利用此方法对2005年6月~2006年3月间台湾海峡西岸重力场变化进行数值模拟,获得其空间梯度场图像。     

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.     

Correcting systematic effects in a large set of photometric light curves

We suggest a new algorithm to remove systematic effects in a large set of light curves obtained by a photometric survey. The algorithm can remove systematic effects, such as those associated with atmospheric extinction, detector efficiency, or point spread function changes over the detector. The algorithm works without any prior knowledge of the effects, as long as they linearly appear in many stars of the sample. The approach, which was originally developed to remove atmospheric extinction effects, is based on a lower rank approximation of matrices, an approach which has already been suggested and used in chemometrics, for example. The proposed algorithm is especially useful in cases where the uncertainties of the measurements are unequal. For equal uncertainties, the algorithm reduces to the Principal Component Analysis (PCA) algorithm. We present a simulation to demonstrate the effectiveness of the proposed algorithm and we point out its potential, in the search for transit candidates in particular.     

The NAIAD experiment for WIMP searches at Boulby mine and recent results

The NAIAD experiment (NaI Advanced Detector) for weakly interacting massive particle (WIMP) dark matter searches at Boulby mine (UK) is described. The detector consists of an array of encapsulated and unencapsulated NaI(Tl) crystals with high light yield. Six crystals are collecting data at present. Data accumulated by four of them (10.6 kg × year exposure) have been used to set upper limits on the WIMP–nucleon spin-independent and WIMP–proton spin-dependent cross-sections. Pulse shape analysis has been applied to discriminate between nuclear recoils, as may be caused by WIMP interactions, and electron recoils due to gamma background. Various calibrations of crystals are presented.     

Bond rolling resistance and its effect on yielding of bonded granulates by DEM analyses

A discrete element modelling of bonded granulates and investigation on the bond effect on their behaviour are very important to geomechanics. This paper presents a two‐dimensional (2‐D) discrete element theory for bonded granulates with bond rolling resistance and provides a numerical investigation into the effect of bond rolling resistance on the yielding of bonded granulates. The model consists of mechanical contact models and equations governing the motion of bonded particles. The key point of the theory is that the assumption in the original bond contact model previously proposed by the authors (55th CSCE‐ASCE Conference, Hamilton, Ont., Canada, 2002; 313–320; J. Eng. Mech. (ASCE) 2005; 131 (11):1209–1213) that bonded particles are in contact at discrete points, is here replaced by a more reliable assumption that bonded particles are in contact over a width. By making the idealization that the bond contact width is continuously distributed with the normal/tangential basic elements (BE) (each BE is composed of spring, dashpot, bond, slider or divider), we establish a bond rolling contact model together with bond normal/tangential contact models, and also relate the governing equations to local equilibrium. Only one physical parameter β needs to be introduced in the theory in comparison to the original bond discrete element model. The model has been implemented into a 2‐D distinct element method code, NS2D. Using the NS2D, a total of 86 1‐D, constant stress ratio, and biaxial compressions tests have been carried out on the bonded granular samples of different densities, bonding strengths and rolling resistances. The numerical results show that: (i) the new theory predicts a larger internal friction angle, a larger yielding stress, more brittle behaviour and larger final broken contact ratio than the original bond model; (ii) the yielding stress increases nonlinearly with the increasing value of β, and (iii) the first‐yield curve (initiation of bond breakage), which define a zone of none bond breakage and which shape and size are affected by the material density, is amplified by the bond rolling resistance in analogous to that predicted by the original bond model. Copyright © 2006 John Wiley & Sons, Ltd.     

Virtual Huanghe River System： Framework and Technology

1 Introduction Huanghe (Yellow) River basin is located in 32°–42°N, 96°–119°E. The area of the catchment is more than 752,000km2. The river is 5464km long with a drop in elevation of 4830m. Among the whole area, the moun- tainous and stone area accounts for 29%, loess and hills area 46%, sandy area 11% and plain area 14%, respec- tively. Different natural landscapes exist in this area. The Huanghe River flows through the Loess Plateau, where the soil is eroded seriously (Wang, 2002;…