行星磁场与卫星引力及其相对运动的统计关系

从知道某些天体具有磁场起，人们就对其磁场的起源提出种种解释，例如有电池说，转子说，化石说，发电机说等等。但由于这些学说都分别与某些观测事实相抵触而未被公认，因此星球磁场的起源一直是未能解决的命题。余先河先生提出，星球的磁场起源可能与所受的引力有关，行星的磁场强度正比于其所受卫星的引力；正比于卫星与行星的引力连线转动的相对角速度。本文对这两方面的命题分别进行了相关分析，得到其相关系数分别为：ｒ＝０．８４８１和ｒ＝０．８４２５，它们都在ａ＝０．０１的信度水平上相关。结果表明余新河关于行星磁场起源的设想是有基础的。本文还对其统计结果和可能机制进行了讨论。     

洛坝铅锌矿床的构造演化分析

洛坝铅锌矿床位于黄诸关韧脆性变形带所波及的范围内,断裂、褶皱、劈理、节理等形态构造类型均很发育,经对小型构造的统计及应力场分析证明,矿区变形主要是在海西—印支期长期存在的南北向古应力场持续作用下形成的。     

Modelling Gas And Dust Release From Comet Hale–Bopp

Numerical simulations of the evolving activity of comet Hale-Bopp are presented, assuming a porous, spherical nucleus, 20 km in radius, made of dust and gas-laden amorphous ice. The main effects included are: crystallization of amorphous ice and release of occluded gas, condensation, sublimation and flow of gases through the pores, changing pore sizes, and flow of dust grains. The model parameters, such as initial pore size and porosity, emissivity, dust grain size, are varied in order to match the observed activity. In all cases, a sharp rise in the activity of the nucleus occurs at a large heliocentric distance pre-perihelion, marked by a few orders of magnitude increase in the CO and the CO₂ fluxes and in the rate of dust emission. This is due to the onset of crystallization, advancing down to a few meters below the surface, accompanied by release of the trapped gases. A period of sustained, but variable, activity ensues. The emission of water molecules is found to surpass that of CO at a heliocentric distance of 3 AU. Thereafter the activity is largely determined by the behaviour of the dust. If a dust mantle is allowed to build up, the water production rate does not increase dramatically towards perihelion; if most of the dust is ejected, the surface activity increases rapidly, producing a very bright comet. This revised version was published online in July 2006 with corrections to the Cover Date.     

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

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

Transient‐state groundwater flow in various geometries with wells at arbitrary positions: analytical methods

We have developed a method for analytically solving the porous medium flow equation in many different geometries for horizontal (two‐dimensional), homogeneous and isotropic aquifers containing impermeable boundaries and any number of pumping or injection wells located at arbitrary positions within the system. Solutions and results are presented for rectangular and circular aquifers but the method presented here is easily extendible to many geometries. Results are also presented for systems where constant head boundary conditions can be emulated internal to the aquifer boundary. Recommendations for extensions of the present work are briefly discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

The evolution of binary fractions in globular clusters

We study the evolution of binary stars in globular clusters using a new Monte Carlo approach combining a population synthesis code ( startrack ) and a simple treatment of dynamical interactions in the dense cluster core using a new tool for computing three- and four-body interactions ( fewbody ). We find that the combination of stellar evolution and dynamical interactions (binary–single and binary–binary) leads to a rapid depletion of the binary population in the cluster core. The maximum binary fraction today in the core of a typical dense cluster such as 47 Tuc, assuming an initial binary fraction of 100 per cent, is only ∼ 5–10 per cent. We show that this is in good agreement with recent Hubble Space Telescope observations of close binaries in the core of 47 Tuc, provided that a realistic distribution of binary periods is used to interpret the results. Our findings also have important consequences for the dynamical modelling of globular clusters, suggesting that 'realistic models' should incorporate much larger initial binary fractions than has usually been the case in the past.     

Atmospheric radiative transfer simulation for atmospheric correction of remote sensing data

IntroductionThe radiance leaving the earth-atmosphere sys-tem which can be sensed by a satellite borne ra-diometer is the sum of radiation emission fromtheearth surface and each atmospheric level that aretransmittedtothe top of the atmosphere.The radia-tion emissionfromthe earthsurface andthe radianceof each atmospheric level can be separated fromtheradiance at the top the atmospheric level a satellitemeasured.Thus,the earth surface parameters willbe retrieved from the surface radiance after a…     

利用Lageos 1 SLR资料解算1990-2001年的地球定向参数

利用 12年的Lageos 1卫星激光测距资料 (1990 - 2 0 0 1)解算得到了地球定向参数 (EOP) ,将该序列的结果与同期的EOP(IERS)C0 4进行比较 ,其外符精度为 :极移XP— 0 .4 0mas,YP— 0 .4 2mas ,日长变化Dr— 0 .0 35ms。