一维扩散方程整体观测资料下的初值变分同化反演

本文考虑一维扩散方程的反问题,利用变分同化方法通过观测资料来确定方程中的未知初值,通过分析观测误差对于初值误差的影响,证明变分同化初值收敛于原问题的真实参数,并得到了参数的收敛精度。同时将得到的初值代入预报模式中,得到预报解,并分析了预报解的收敛性和预报误差。     

Mais comment s'écoule donc un glacier ? Aperçu historique

Ice and snow have often helped physicists understand the world. On the contrary it has taken them a very long time to understand the flow of the glaciers. Naturalists only began to take an interest in glaciers at the beginning of the 19th century during the last phase of glacier advances. When the glacier flow from the upslope direction became obvious, it was then necessary to understand how it flowed. It was only in 1840, the year of the Antarctica ice sheet discovery by Dumont d'Urville, that two books laid the basis for the future field of glaciology: one by Agassiz on the ice age and glaciers, the other one by canon Rendu on glacier theory. During the 19th century, ice flow theories, adopted by most of the leading scientists, were based on melting/refreezing processes. Even though the word ‘fluid’ was first used in 1773 to describe ice, more the 130 years would have to go by before the laws of fluid mechanics were applied to ice. Even now, the parameter of Glen's law, which is used by glaciologists to model ice deformation, can take a very wide range of values, so that no unique ice flow law has yet been defined. To cite this article: F. Rémy, L. Testut, C. R. Geoscience 338 (2006).     

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.     

剥蚀及地幔作用下青藏高原隆升过程的数值模拟

修改了England和Mckenzie的黏性薄层流变模型中控制大陆形变的连续性方程,将剥蚀作用对高原隆升演化的影响直接引入该方程,并考虑下伏地幔小尺度对流对增厚岩石层的搬离作用对高原隆升演化后期的影响,用有限差分法直接模拟青藏高原隆升过程. 数值模拟结果所显示的高原隆升演化过程与实际观测资料吻合较好,揭示了高原隆升演化过程的非平稳和多阶段的特性;同时还表明上地幔小尺度对流对岩石层底部的搬离作用可能是最近8Ma以来高原快速隆升的主导机制.     

带阻尼的一维地震反射波成像

本文是文献［８］的继续,文中提出了地震勘探的带阻尼波动方程反演的新方法,应用此反演算法可以通过微机识别地层的波阻抗,这里把地层分成１００层,地震数据含有１０％的误差,且反演计算时间仅仅是正演计算时间的三倍．文献［１］－［４］的计算方法是基于波动方程的混合问题,而文献［８］所提出的反演算法是基于波动方程的哥西问题,许多例子充分说明后者更有效．     

甘肃河西走廊沙尘暴与赤道中,东太平洋海温之间的遥相关分析

春、夏季甘肃河西走廊沙尘暴发生次数与赤道中、东太平洋海温之间遥相关分析结果表明：春、夏季河西走廊沙尘暴发生次数与前二年赤道中、东太平洋海温的负相关最好,时间过远或过近,相关性反而较差。利用前期的赤道中、东太平洋海温可以较好地预报出河西走廊春季沙尘暴发生的多寡。     

非线性洋面温度反演

本文从辐射方程出发，简述洋面温度与亮度温度之间的非线性关系，进行非线性洋面温度反演的方案设计，用正演算法 计算出非线性洋面温度反演公式的回归系数，对样进行检验。     

幂函数剪切模量成层土非线性地震反应的半解析算法

运用文献[1]所建议的动态应力一应变关系及其推广的Masing加卸载准则,考虑土料在地震等产生的不规则加载条件下的非线性滞回特征,将增量法与相应场地地震线性反应解析解[2]相结合,提出了该动力非线性方程的半解析时域算法,基于改进的一维剪切梁模型,对剪切模量是其深度的某一幂函数的成层非均质土层,建立了求解土体地震反应的非线性分析技术。针对文献[2]中的土层剖面,做了计算、分析和讨论。