长江流域蒸发皿蒸发量的区域变化特征及影响因素

基于长江流域148个气象站1980—2017年的蒸发皿观测数据,将旋转经验正交函数（REOF）和模糊C均值聚类（FCM）相结合,对流域蒸发皿蒸发量（PE）进行分区,然后运用Modified Mann-Kendall检验和多元逐步回归等方法,分析各子区域PE的变化特征并识别其主要影响因子。结果表明：1) REOF的前4个空间模态显示流域PE存在5个主要的异常敏感区,基于这4个空间模态,流域PE在空间上可划分为9个子区域。2）在年尺度上,各区PE呈不同程度的增加趋势,其中中部盆地区上升速率最大（111.28 mm/10 a）,西部高原区上升速率最小（12.5 mm/10 a）;而在季节尺度上,秋、冬季流域PE呈显著上升趋势,春、夏季PE的变化具有明显的区域差异性,部分地区PE为下降趋势。3）影响PE变化的主要因子因地而异,但大多子区域的PE变化与平均气温和饱和水汽压差的变化显著相关。     

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

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

1958-2017年太湖蒸发量年际变化趋势及主控因子

湖泊蒸发对气候变化非常敏感,是水文循环响应气候变化的指示因子,因此研究湖泊蒸发的控制因素,对于理解区域水文循环有重要意义.本文利用太湖中尺度涡度通量网避风港站观测数据校正JRA-55再分析资料,驱动CLM4.0-LISSS模型,并利用2012-2017年涡度相关通量数据和湖表面温度数据检验模型模拟蒸发结果,验证了该模型在太湖的适用性;估算了1958-2017年间太湖的湖面蒸发量,并利用Manner-Kendall趋势检验分析了湖面蒸发的变化趋势,寻找太湖实际蒸发的年际变化的主控因子.结果如下:校正后的JRA-55再分析资料模拟的太湖蒸发与观测值之间存在季节偏差,但是季节偏差在年尺度上相互抵消,再分析资料可用于年际尺度太湖蒸发变化的模拟;1958-2017年间太湖蒸发量以1977年为界,先下降(-3.6 mm/a),后增加(2.3 mm/a);多元逐步回归结果表明,向下的短波辐射是太湖1958-2017年间太湖蒸发变化的主控因子,向下的长波辐射、气温、比湿也对湖泊蒸发年际变化有一定影响,但是风速对蒸发量的年际变化影响不大.     

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.     

On the role of soil moisture in daytime evolution of temperatures

This paper explores the relationship between temperature, evaporation and soil moisture using a planetary boundary layer (PBL) model. It focuses on illustrating and quantifying the effect of soil moisture on the evolution of daytime temperatures. A simple convective PBL model coupled with the Penman–Monteith (PM) equation is used to estimate evapotranspiration. Following calibration and sensitivity analysis, the model was used to simulate the relative impact of dry and wet soil moisture conditions on daytime temperatures by changing the surface resistance parameter in the PM equation. It was found that the maximum temperature that can be reached during a day is constrained by the amount of soil moisture and the available net radiation, confirming previously published results. Higher temperatures can be reached with greater net radiation and dry soil moisture conditions. Copyright © 2012 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.     

A formulation for diagnostic anomaly models

It is shown that the anomaly model governing the difference between two integrations of the equations of motion can be written in a form that is formally linear in the anomalies, by choosing the basic state to be the mean of the results of the two integrations.