大别山及其周边现今垂直形变研究

基于1980年、2005年和2015年3期区域精密水准网观测资料,利用线性动态平差模型计算获取大别山地区水准网长期垂直运动速度场图像。研究发现,淮河平原地表下沉较为严重,大别山呈现弱隆升趋势运动,长江谷地边缘地区较大别山区呈现明显的隆升运动。跨郯庐断裂带水准剖面结果显示,垂直运动与地形呈负相关和弱相关。     

捷联式航空重力扰动矢量解算及误差补偿

介绍捷联式航空重力矢量测量的基本原理,利用波数相关滤波（WCF）、比力线性校正的方法对重力扰动的水平和垂向分量进行误差补偿,评定内符合精度。对国产捷联式航空重力矢量仪的某次试验数据进行处理,以验证和评估航空重力矢量仪的性能。结果表明,在半波长分辨率为7.5 km时,6条重复测线的重力扰动水平分量经波数相关滤波处理后,东、北向分量的平均内符合精度分别从9.77 mGal、9.18 mGal提高到5.95 mGal、3.83 mGal;对比力的垂向分量线性校正后,将解算的重力扰动垂向分量再用WCF方法处理,其平均内符合精度从1.27 mGal提高到0.59 mGal。     

Dominant climate factors influencing the Arctic runoff and association between the Arctic runoff and sea ice

By using the Arctic runoff data from R-ArcticNET V4.0 and ArcticRIMS, trends of four major rivers flowing into the Arctic Ocean, whose climate factor plays an important role in determining the variability of the Arctic runoff, are investigated. The results show that for the past 30 years, the trend of the Arctic runoff is seasonally dependent. There is a significant trend in spring and winter and a significant decreasing trend in summer, leading to the reduced seasonal cycle. In spring, surface air temperature is the dominant factor influencing the four rivers. In summer, precipitation is the most important factor for Lena and Mackenzie, while snow cover is the most important factor for Yenisei and Ob. For Mackenzie, atmospheric circulation does play an important role for all the seasons, which is not the case for the Eurasian rivers. The authors further discuss the relationships between the Arctic runoff and sea ice. Significant negative correlation is found at the mouth of the rivers into the Arctic Ocean in spring, while significant positive correlation is observed just at the north of the mouths of the rivers into the Arctic in summer. In addition, each river has different relationship with sea ice in the eastern Greenland Sea.     

双线偏振雷达判别降水粒子类型技术及其检验

在综合分析国内外降水粒子相态识别方法的基础上,结合我国双线偏振雷达的应用需求,建立了一套利用双线偏振雷达探测资料识别降水粒子类型的识别模式,并根据美国KOUN雷达观测资料对该模式进行了分析验证,同时检验了模式中的几个偏振参数:差分反射率因子ZDR、差传播相移KDP、水平偏振和垂直偏振回波功率零滞后互相关系数ρHV(0)对不同降水粒子的识别效果。通过分析认为:利用该识别模式得到的结果基本合理,可以反映降水区域内降水粒子的相态结构,但还需要对资料作进一步的研究。同时在识别各类降水粒子方面,ZDR的效果要高于KDP和ρHV(0),在缺乏相关资料的情况下可以考虑用ZH和ZDR配对识别的结果作为真实的结果。     

应用线性合成概率方法探讨青藏高原北部地区强震趋势

基于周期谱分析基础上的线性合成概率预测方法,对青藏高原北部地区Ms≥7.0、Ms≥6.5强震进行了试验性研究,并利用其结果作为预测依据,对青藏高原北部地区强震趋势做了初步估计.利用该方法,本文还对东昆仑断裂带1900年以来的Ms≥7.0强震进行了研究.     

Investigation of highly efficient algorithms for solving linear equations in the discontinuous deformation analysis method

Large‐scale engineering computing using the discontinuous deformation analysis (DDA) method is time‐consuming, which hinders the application of the DDA method. The simulation result of a typical numerical example indicates that the linear equation solver is a key factor that affects the efficiency of the DDA method. In this paper, highly efficient algorithms for solving linear equations are investigated, and two modifications of the DDA programme are presented. The first modification is a linear equation solver with high efficiency. The block Jacobi (BJ) iterative method and the block conjugate gradient with Jacobi pre‐processing (Jacobi‐PCG) iterative method are introduced, and the key operations are detailed, including the matrix‐vector product and the diagonal matrix inversion. Another modification consists of a parallel linear equation solver, which is separately constructed based on the multi‐thread and CPU‐GPU heterogeneous platforms with OpenMP and CUDA, respectively. The simulation results from several numerical examples using the modified DDA programme demonstrate that the Jacobi‐PCG is a better iterative method for large‐scale engineering computing and that adoptive parallel strategies can greatly enhance computational efficiency. Copyright © 2015 John Wiley & Sons, Ltd.     

Elastic fields in two joined transversely isotropic media of infinite extent as a result of rectangular loading

This paper presents the closed‐form solutions for the elastic fields in two bonded rocks induced by rectangular loadings. Each of the two bonded rocks behaves as a transversely isotropic linear elastic solid of semi‐infinite extent. They are completely bonded together at a horizontal surface. The rectangular loadings are body forces along either vertical or horizontal directions and are uniformly applied on a rectangular area. The rectangular area is embedded in the two bonded rocks and is parallel to the horizontal interface. The classical integral transforms are used in the solution formulation, and the elastic solutions are expressed in the forms of elementary harmonic functions for the rectangular loadings. The stresses and displacements in the rocks induced by both the horizontal and vertical body forces are also presented. The numerical results illustrate the important effect of the anisotropic bimaterial properties on the stress and displacement fields. The solutions can be easily implemented for numerical calculations and applied to problems encountered in rock mechanics and engineering. Copyright © 2011 John Wiley & Sons, Ltd.     

Abrupt behaviours of streamflow and sediment load variations of the Yangtze River basin,China

Monthly sediment load and streamflow series spanning 1963–2004 from four hydrological stations situation in the main stem of the Yangtze River, China, are analysed using scanning t‐test and the simple two‐phase linear regression scheme. Results indicate significant changes in the sediment load and streamflow from the upper reach to the lower reach of the Yangtze River. Relatively consistent positive coherency relations can be detected between streamflow and sediment load in the upper reach and negative coherency in the middle and lower reaches. Interestingly, negative coherency is found mainly for larger time scales. Changes in sediment load are the result mainly of human influence; specifically, the construction of water reservoirs may be the major cause of negative coherency. Accentuating the human influence from the upper to the lower reach results in inconsistent correlations between sediment load and streamflow. Decreasing sediment load being observed in recent years has the potential to alter the topographical properties of the river channel and the consequent development and recession of the Yangtze Delta. Results of this study are of practical significance for river channel management and evaluation of the influence of human activities on the hydrological regimes of large rivers. Copyright © 2012 John Wiley & Sons, Ltd.     

Preconditioned IDR(s) iterative solver for non‐symmetric linear system associated with FEM analysis of shallow foundation

Non‐associated flow rule is essential when the popular Mohr–Coulomb model is used to model nonlinear behavior of soil. The global tangent stiffness matrix in nonlinear finite element analysis becomes non‐symmetric when this non‐associated flow rule is applied. Efficient solution of this large‐scale non‐symmetric linear system is of practical importance. The standard Krylov solver for a non‐symmetric solver is Bi‐CGSTAB. The Induced Dimension Reduction [IDR(s)] solver was proposed in the scientific computing literature relatively recently. Numerical studies of a drained strip footing problem on homogenous soil layer show that IDR(s = 6) is more efficient than Bi‐CGSTAB when the preconditioner is the incomplete factorization with zero fill‐in of global stiffness matrix K_ep (ILU(0)‐K_ep). Iteration time is reduced by 40% by using IDR(s = 6) with ILU(0)‐K_ep. To further reduce computational cost, the global stiffness matrix K_ep is divided into two parts. The first part is the linear elastic stiffness matrix K_e, which is formed only once at the beginning of solution step. The second part is a low‐rank matrix Δ, which is re‐formed at each Newton–Raphson iteration. Numerical studies show that IDR(s = 6) with this ILU(0)‐K_e preconditioner is more time effective than IDR(s = 6) with ILU(0)‐K_ep when the percentage of yielded Gauss points in the mesh is less than 15%. The total computation time is reduced by 60% when all the recommended optimizing methods are used. Copyright © 2013 John Wiley & Sons, Ltd.     

An application of Karmarkar's interior-point linear programming algorithm for multi-reservoir operations optimization

Optimization of multi-reservoir systems operations is typically a very large scale optimization problem. The following are the three types of optimization problems solved using linear programming (LP): (i) deterministic optimization for multiple periods involving fine stage intervals, for example, from an hour to a week (ii) implicit stochastic optimization using multiple years of inflow data, and (iii) explicit stochastic optimization using probability distributions of inflow data. Until recently, the revised simplex method has been the most efficient solution method available for solving large scale LP problems. In this paper, we show that an implementation of the Karmarkar's interior-point LP algorithm with a newly developed stopping criterion solves optimization problems of large multi-reservoir operations more efficiently than the simplex method. For example, using a Micro VAX II minicomputer, a 40 year, monthly stage, two-reservoir system optimization problem is solved 7.8 times faster than the advanced simplex code in MINOS 5.0. The advantage of this method is expected to be greater as the size of the problem grows from two reservoirs to multiples of reservoirs. This paper presents the details of the implementation and testing and in addition, some other features of the Karmarkar's algorithm which makes it a valuable optimization tool are illuminated.