胶州湾三维变动边界的潮流数值模拟

基于Blumberg等（1986）的河口、陆架和海洋模式，引入变边界处理技术，建立了胶州湾三维变动边界的潮流模型，模型以σ-正交曲线坐标下三维非线性潮波方程为基本方程，引入2.5阶瑞封闭方程组，采用分裂算子法数值求解方程组，利用湍封闭方程求解垂直紊动粘滞系数，采用干湿网格方法模拟潮流漫滩过程，三维变动边界潮流模型计算结果与实测值吻合良好。     

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.     

2.5D direct‐current resistivity forward modelling and inversion by finite‐element–infinite‐element coupled method

To reduce the numerical errors arising from the improper enforcement of the artificial boundary conditions on the distant surface that encloses the underground part of the subsurface, we present a finite‐element–infinite‐element coupled method to significantly reduce the computation time and memory cost in the 2.5D direct‐current resistivity inversion. We first present the boundary value problem of the secondary potential. Then, a new type of infinite element is analysed and applied to replace the conventionally used mixed boundary condition on the distant boundary. In the internal domain, a standard finite‐element method is used to derive the final system of linear equations. With a novel shape function for infinite elements at the subsurface boundary, the final system matrix is sparse, symmetric, and independent of source electrodes. Through lower upper decomposition, the multi‐pole potentials can be swiftly obtained by simple back‐substitutions. We embed the newly developed forward solution to the inversion procedure. To compute the sensitivity matrix, we adopt the efficient adjoint equation approach to further reduce the computation cost. Finally, several synthetic examples are tested to show the efficiency of inversion.     

保山井水位异常的数值模拟

利用1995年7月12日云南孟连Ms7.3地震和1996年2月3日云南丽江Ms7．0地震前保山井水位异常变化的原始资料,并根据介质体膨胀变化率和水位变化的耦合关系式,用有限差分法模拟了保山井水位的变化情况。通过拟合实际观测曲线,得出观测井附近含水层内平均应力变化率的演化情况。     

Centennial to millennial geomagnetic secular variation

A time-varying spherical harmonic model of the palaeomagnetic field for 0–7 ka is used to investigate large-scale global geomagnetic secular variation on centennial to millennial scales. We study dipole moment evolution over the past 7 kyr, and estimate its rate of change using the Gauss coefficients of degree 1 (dipole coefficients) from the CALS7K.2 field model and by two alternative methods that confirm the robustness of the predicted variations. All methods show substantial dipole moment variation on timescales ranging from centennial to millennial. The dipole moment from CALS7K.2 has the best resolution and is able to resolve the general decrease in dipole moment seen in historical observations since about 1830. The currently observed rate of dipole decay is underestimated by CALS7K.2, but is still not extraordinarily strong in comparison to the rates of change shown by the model over the whole 7 kyr interval. Truly continuous phases of dipole decrease or increase are decadal to centennial in length rather than longer-term features. The general large-scale secular variation shows substantial changes in power in higher spherical harmonic degrees on similar timescales to the dipole. Comparisons are made between statistical variations calculated directly from CALS7K.2 and longer-term palaeosecular variation models: CALS7K.2 has lower overall variance in the dipole and quadrupole terms, but exhibits an imbalance between dispersion in   g ¹₂  and   h ¹₂  , suggestive of long-term non-zonal structure in the secular variations.     

Discussion of ‘Modelling viscous damping in nonlinear response history analysis of buildings for earthquake excitation’ by Anil K. Chopra and Frank McKenna

This discussion deals with recommendations in the paper on appropriate damping formulations for use in nonlinear response history analysis of buildings. Concern over potentially excessive damping forces and moments should extend beyond the damping moments produced by the stiffness proportional part of Rayleigh damping that corresponds to rotational springs used to explicitly model plastic hinges. The key to an appropriate damping formulation for nonlinear analysis is a realistic mechanism that allows all damping forces and moments to be meaningfully assessed. Then features can be added to keep these forces and moments within reasonable bounds. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical simulation of sloshing waves in a 3D tank based on a finite element method

The sloshing waves in a three dimensional (3D) tank are analysed using a finite element method based on the fully non-linear wave potential theory. When the tank is undergoing two dimensional (2D) motion, the calculated results are found to be in very good agreement with other published data. Extensive calculation has been made for the tank in 3D motion. As in 2D motion, in addition to normal standing waves, travelling waves and bores are also observed. It is found that high pressures occur in various circumstances, which could have important implications for many engineering designs.     

东亚夏季风降水和环流数值模拟对不同对流参数化方案的敏感性

A 5-level spectral AGCM (ImPKU-5LAGCM) is used to examine the sensitivity of the simulated results of the summer monsoon rainfall and circulation in East Asia to different cumulus parameterization schemes in the climatological-mean case and in the cases of weak and strong Asian summer monsoons,respectively. The results simulated with the Arakawa-Schubert＇s(hereafter A-S＇s), Kuo＇s and Manabe＇s cumulus parameterization schemes show that these simulated distributions of the summer monsoon rainfall and circulation in East Asia depend strongly on the cumulus parameterization schemes either in the climatological-mean case or in the cases of weak and strong Asian summer monsoons. From the simulated results, it might be shown that the Kuo scheme appears to be more suitable for the simulation of the summer monsoon rainfall and circulation in East Asia than the A-S scheme or the Manabe scheme, although the A-S scheme is somewhat better in the simulations of the tropical rainfall. This might be due to that the Kuo＇s cumulus parameterization scheme is able to reflect well the characteristics of rainfall cloud system in the East Asian summer monsoon region, where the rainfall system used to be a mixing of cumulus and stratus.     

黄渤海带鱼(Trichiurus lepturus)渔业生物学参考点估计研究

It is important to find a reliable method to estimate maximum sustainable yield(MSY) or total allowable catch(TAC) for fishery management, especially when the data availability is limited which is a case in China. A recently developed method(CMSY) is a data-poor method, which requires only catch data, resilience and exploitation history at the first and final years of the catch data. CMSY was used in this study to estimate the biological reference points for Largehead hairtail(Trichiurus lepturus, Temminck and Schlegel) in the Yellow Sea and Bohai Sea, based on the fishery data from China Fishery Statistical Year Books during 1986 to 2012. Additionally,Bayesian state-space Schaefer surplus production model(BSM) and the classical surplus production models(Schaefer and Fox) performed by software CEDA and ASPIC, were also projected in this study to compare with the performance of CMSY. The estimated MSYs from all models are about 19.7×104–27.0×104 t, while CMSY and BSM yielded more reasonable population parameter estimates(the intrinsic population growth rate and the carrying capacity). The biological reference points of B/BMSY smaller than 1.0, while F/FMSY higher than 1.0 revealed an over-exploitation of the fishery, indicating that more conservative management strategies are required for Largehead hairtail fishery.     

水汽方程积分计算格式的数值试验

该文把Easter(1993)的水汽方程积分方案应用到球坐标的有限区模式（LAFS）中,并作了两组比较性试验。Easter的方案是把水汽的通量方程和连续方程联立求解。在计算中不会出现水汽负值。而且稳定性好。本文所做的两组试验表明。Easter的方案比LAFS的原方案所产生的计算性扩散小。因而降水预报更接近实况。但CPU时间增加了16%。