Utilization of water resources,ecological balance and land desertification in the Tarim Basin,Xinjiang

Because of the human exploitation and utilization of water resources in the Tarim Basin, the water resources consumption has changed from mainly natural ecosystem to artificial oasis ecosystem, and the environment has changed correspondingly. The basic changes are: desertification and oasis development coexist, both “the human being advance and the desert retreat” and “the desert advance and the human being retreat” coexist, but the latter is dominant. In the upper reaches, water volume drawing to irrigated agricultural areas has increased, artificial oases have been enlarging and moving from the deltas in the lower reaches of many rivers to the piedmont plains. In the middle and lower reaches of the Tarim River, the stream flow has decreased, old oases have declined, natural vegetations have been degenerating, desertification has been enlarging, and the environment has deteriorated. The transition regions, which consist of forestlands, grasslands and waters between the desert and the oases, have been decreasing continuously, their shelter function to the oases has been weakened, and the desert is threatening the oases seriously.     

多种驱动力作用下东亚大陆形变及应力场演化

将印度板块持续地向北推进、下伏地幔小尺度对流对增厚大陆岩石层的搬离作用以及剥蚀作用视为形成现今东亚大陆形变和应力场格局的主驱动力。在一梯形区域内，利用数值模拟的方法，研究了东亚大陆在不同的边界条件、不同的剥蚀率系数及不同的岩石力学参数条件下的形变及应力场格局。与现代空间大地测量技术 CGPS)以及利用地震观测得到的结果进行了对比。结果表明，本文模型预测的结果与上述的观测结果有较好的吻合，其西部地区比东部吻合得更好。说明控制东亚大陆西部形变和应力场基本格局的主驱动力，来源于印度板块对欧亚板块的碰撞、挤压，而对东部地区还应当考虑其与太平洋板块和菲律宾板块的相互作用。与此同时，下伏地幔小尺度对流对增厚大陆岩石层的搬离作用以及风化剥蚀对应力场的演化过程也不可忽视。      

2001年四川省雅江-康定6.0级地震烈度考察及发震构造背景雏议

2001年2月23日四川省发生的6.0级地震，其宏观震中位于雅江县与康定县之间的高山峡谷地带．极震区烈度可达Ⅷ 度，Ⅷ 度区、Ⅶ 度区和Ⅵ 度区面积分别为180km2，1472 km2和3998 km2，等震线分布总体形态呈椭圆形，长轴近南北向．震区建筑物普遍受损，极震区山地灾害较为严重．初步分析认为，此次地震的发生与理塘断裂和玉农希断裂的长期活动有关，主干断裂活动造成的地壳应力调整和集中，最终在块体内部的次级断裂上释放而产生地震．      

由GPS观测结果推导中国大陆现今水平应变场

以中国大陆及周边近400个GPS测站的水平运动速率为基础，给出了现今地壳水平应变场结果表明:①中国大陆水平应变为西强东弱，剪应变数值大于正应变数值(绝对值)，应变量级一般为10-8/a，局部区域达到10-7/a，但应变分布不均匀；②南北向应变最突出的部位为中国西南部西段的喜马拉雅条带、西北部的36N~42N段及柴达木断块的北缘；③东西向应变西边缘变化最大．此外，由西向东还具有正负交替的变化特征；④REN(东-北向剪切应变)与Rmax(最大剪切应变)数值较大的区域分别是喜马拉雅条带、西北部的36N～42N段、柴达木断块的西部、川滇菱形块体，以及阿拉善、祁连及塔里木断块的交界区；⑤青藏块体周边以面收缩为主，内部则以面膨胀为主．其以北的地区以面收缩为主．西界数值最大，东部数值最小(除燕山构造带外)；⑥西部西区主压应变为南北向，主张应变为东西向．西部东缘区主压应变为近东西向，主张应变为近南北向．川滇菱形块体主应变的方向发生了很大的变化，北部地区为东西压南北张，南部地区则恰好相反；⑦中国大陆的应变模式可能是断块模式与连续模式的组合．此外，小尺度优势应变可能是剪切应变．造成上述结果除与印度板块的碰撞及边界耦合有关外，还与深部物质的活动及地壳介质的物性有密切的关系．必须指出，由于GPS测站在空间上分布的不均匀性，那么，由此而来的应变场，其应变尺度也不一样．      

中国与美国西部平均构造剪应力值的比较

提出了一种估计一个地区平均构造剪应力值的方法．用这种方法，选取哈佛大学公布的1977～1999年共15 993次地震的矩张量数据，估计了美国西部19个地区、中国和邻区43个地区(每个地区为1010范围)的平均环境剪应力值．结果表明，美国西部南加州西部海域和南加州应力最高，达13.7和12.0 MPa，然后向北、向南和向东逐渐递减，但最小也达8.7 MPa，是最高值的63%．中国新疆西北部地区和西藏察隅地区应力水平最高，达17.2和12.9 MPa，比美国高．中国的华北、云南、四川、台湾和美国南加州的应力水平差不多．中国南北地震带的应力水平为13 MPa左右，比南加州略高．两个重要地区的平均剪应力值分布图，提供了地学的基本数据．这些结果可为研究地震活动的大背景提供依据，对研究强地面运动参数(如峰值加速度及反应谱等)的衰减关系也是有用的．      

用应变积累释放模型研究中国大陆地块分区地震活动

依据活动地块假说及活动边界的划分，对中国大陆地震活动进行分区．在此基础上，应用应变积累释放模型，通过地震应变积累释放的计算，研究了中国大陆各活动地块分区的地震活动性，初步探讨了各地块的地震活动阶段划分．此外，还对应变积累释放模型的应用前提条件进行了必要的讨论，并对模型结果中可能出现的问题给予一定分析．      

Efficient methods for solving water flow in variably saturated soils under prescribed flux infiltration

The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.     

Screening model for volatile pollutants in dual porosity soils

This paper develops mass fraction models for transport and fate of agricultural pollutants in structured two-region soils. Mass fraction index models, based on a semi-infinite domain solution, are derived that describe leaching at depth, vapor losses through soil surface, absorption, and degradation in the dynamic- and stagnant-water soil regions. The models predict that leaching is the result of the combined effect of the upward vapor-phase transport relative to downward advection, residence time relative to half-life, dispersion, and lateral diffusive mass transfer. Simulations show that leached fraction of volatile compounds does not always decrease monotonically with increased residence time relative to the pollutant half-life, as a result of complex interactions among the different physical and biochemical processes. The results show that leaching, volatilization, and degradation losses can be affected significantly by lateral diffusive mass transfer into immobile-water regions and advection relative to dispersion (i.e. Peclet number) in the mobile-water regions. It is shown that solute diffusion into the immobile phase and subsequent biochemical decay reduces leaching and vapor losses through soil surface. Potential use of the modified leaching index for the screening of selected pesticides is illustrated for different soil textures and infiltration rates. The analysis may be useful to the management of pesticides and the design of landfills.     

Common-Reflection-Surface Stack for Converted Waves

The finite-offset (FO) common-reflection-surface (CRS) stack has been shown to be able to handle not only P-P or S-S but also arbitrarily converted reflections. It can provide different stack sections such as common-offset (CO), common-midpoint (CMP) and common-shot (CS) sections with significantly increased signal-to-noise ratio from the multi-coverage pre-stack seismic data in a data-driven way. It is our purpose in this paper to demonstrate the performance of the FO CRS stack on data involving converted waves in inhomogeneous layered media. In order to do this we apply the FO CRS stack for common-offset to a synthetic seismic data set involving P-P as well as P-S converted primary reflections. We show that the FO CRS stack yields convincing improvement of the image quality in the presence of noisy data and successfully extracts kinematic wavefield attributes useful for further analyses. The extracted emergence angle information is used to achieve a complete separation of the wavefield into its P-P and P-S wave components, given the FO CRS stacked horizontal and vertical component sections.     

Optimal placement of dampers and actuators based on stochastic approach

A general method is developed for optimal application of dampers and actuators by installing them at optimal location on seismic-resistant structures. The study includes development of a statistical criterion, formulation of a general optimization problem and establishment of a solution procedure. Numerical analysis of the seismic response in time-history of controlled structures is used to verify the proposed method for optimal device application and to demonstrate the effectiveness of seismic response control with optimal device location. This study shows that the proposed method for the optimal device application is simple and general, and that the optimally applied dampers and actuators are very efficient for seismic response reduction.