从地震模拟看匀阻段与大震的关系

运用滑块模型建立了一个地震活动的计算机模拟系统。对地震活动计算机模型进行反复试验，认识到断层失稳过程中摩擦强度分布的重要性。把在断层滑动中摩擦阻力分布比较均匀的段称为匀阻段。对大地震而言，断层摩擦阻力的大或小都不是决定因素，而具有较长的匀阻段才是发生大震的首要条件     

涡旋诱发重联模型（Ⅰ）——动力学特性

本文首先从理论上系统地研究了涡旋诱发重联的基本动力学特性.包括:涡旋诱发重联的物理机制;Kelvin-Helmholtz（K-H）不稳定性和撕裂模不稳定性耦合过程;磁流体涡管的基本特性;不同尺度涡旋间的相互作用,动能和磁能间的转化及准稳态时流线、磁力线、等涡度线和等电流密度线的相似性等.用二维MHD数值模拟方法研究了上述一些特性,并与理论分析结论进行了比较.结果表明,模拟结果证实了理论推论.涡旋诱发重联是产生局部磁场重联的重要机制之一,该模型在磁顶区的应用将另文发表.     

断层地震孕育和发生的不稳定性模型

在以前的工作中,考虑直立走滑型断层地震,假设断层面微元破裂强度遵循Weibull概率分布,由细观力学方法推导出断层面的宏观本构关系是一个非线性函数,表现为弹性-软化塑性特征,在此基础上用稳定性理论研究了地震稳定性问题.而实际断层大多是倾斜的,为此,本文首先建立了由围岩和倾斜断层构成的平面地震力学模型,采用宏观的断层载荷-变形的全过程曲线,详细讨论了倾斜断层地震的不稳定性问题.结果表明,远场一旦施加位移,断层也同时错动,这可能与实际情况不符合.为了更好的模拟断层的初始能量累计过程,进一步对断层本构模型进行改进.考虑断层面破裂强度,采用Coulomb破裂准则,则断层表现为刚塑性本构关系,只有当断层面剪应力达到一个临界值时,断层才开始错动.研究表明,对于倾斜断层地震,与直立走滑型断层地震一样,系统刚度比β（围岩切线刚度与断层刚度最大值之比）是决定地震失稳的重要参数,只有当ββ≥1时,仅仅是断层无震滑动,不会发生地震.在远场应施以位移形式边界条件,以致地震失稳发生在平衡路径的位移转向点并伴有应力突跳.     

地球流体力学的研究与进展

简要介绍中国科学院大气物理研究所七十多年来在理论与计算地球流体力学方面的若干研究及其新的进展.在理论地球流体力学方面,介绍了长波动力学及线性稳定性问题、弱非线性理论及行星波动力学以及用Arnold方法(能量-Casimir方法)研究大气和海洋中各种流体运动的非线性稳定性问题的成果.此外,对扰动演变、扰动和基流相互作用及热带大气动力学中的第二类不稳定条件(CISK)也作了简要的介绍.在计算地球流体力学方面,主要内容包括:用物理观点和数学分析相结合的方法阐述了造成计算紊乱和计算不稳定的机理,论证计算稳定性、算     

On the physical model of earthquake precursor fields and the mechanism of precursors＇timespace distribution──origin and evidences of the strong body earthquake──generating model

Ｏｎｔｈｅｐｈｙｓｉｃａｌｍｏｄｅｌｏｆｅａｒｔｈｑｕａｋｅｐｒｅｃｕｒｓｏｒｆｉｅｌｄｓａｎｄｔｈｅｍｅｃｈａｎｉｓｍｏｆｐｒｅｃｕｒｓｏｒｓ＇ｔｉｍｅｓｐａｃｅｄｉｓｔｒｉｂｕｔｉｏｎ──ｏｒｉｇｉｎａｎｄｅｖｉｄｅｎｃｅｓｏｆ　ｔｈｅｓｔｒｏｎ...     

黑潮延伸体上游中尺度涡场的年代际振荡及其相关机制

黑潮延伸体上游区域的中尺度涡场的涡动能和涡特征尺度存在显著地年代际振荡,和黑潮延伸体路径的年代际变化有很好的相关性。当黑潮延伸体路径比较稳定时,其上游区域涡动能比较高,涡特征尺度比较大,反之相反。通过对黑潮延伸体上游区域的中尺度涡场进行集合分析发现:当黑潮延伸体处于稳定状态时,上游涡场几乎是各向均匀地,有轻微的径向伸长;而当黑潮延伸体处于不稳定状态时,上游的中尺度涡场有显著地纬向伸长。对与中尺度涡场的产生相关的线性斜压不稳定和正压不稳定进行了计算分析,结果显示,线性斜压不稳定不是控制中尺度涡场年代际变化的机制,而正压不稳定对中尺度涡场的年代际变化有积极的贡献。不稳定产生的中尺度涡之间存在非线性涡-涡相互作用。     

Was a dacite dome implicated in the 9,500 <Emphasis Type="SmallCaps">b.p.</Emphasis> collapse of Mt Ruapehu? A palaeomagnetic investigation

Constraining the process by which volcanoes become unstable is difficult. Several models have been proposed to explain the driving forces which cause volcanic edifices to catastrophically collapse. These include models for destabilisation of volcanic flanks by wedging due to dyke intrusion and the weakening of mechanical properties by pressurisation of pore fluids. It is not known which, if any, of the models are relevant to particular sector collapse events. Recent developments in the palaeomagnetic estimation of emplacement temperatures of volcaniclastic rocks have shown that even relatively low emplacement temperatures can be recorded by volcaniclastics with high fidelity. We have carried out a palaeomagnetic study of emplacement temperatures to investigate the role of igneous activity in the initiation of the 9,500 b.p. Murimotu sector collapse of Mt Ruapehu, New Zealand. This debris avalanche deposit has three fades which are stratigraphically superimposed, and the lowermost fades contains three lithological assemblages representing different segments of the edifice which were transported with little internal mixing within the flow. We have determined that some of the dacite-bearing assemblage 1, fades 1 was hot (∼350 °C) during transport and emplacement, whereas none of the other lithological assemblages of fades contained hot material. Our interpretation is that a dacite dome was active on the ancient Ruapehu edifice immediately prior to the Murimotu sector collapse. The partially cooled carapace of the dome and material shed from this part was incorporated into the avalanche deposit, along with cold lavas and volcaniclastics. We have not found evidence for incorporation of material at or close to magmatic temperatures, at least in the sampled locations. Our palaeomagnetic work allows us to develop a comprehensive, new palaeomagnetic classification of volcaniclastics. Published online: 25 January 2003 Editorial responsibility: D. Dingwell     

Baroclinic energetics and zonal plane distribution of monsoon disturbances

The adiabatic, quasi-geostraphic, 25-layer, numerical, linear model with Ekman boundary layer friction is utilised to perform the baroclinic stability analysis of the mean monsoon zonal wind profile. It is shown thec _i is a function of the resultant wavenumber alone. This relation is able to explain the effects of the lateral walls on the unstable waves.The energetics and zonal plane distribution of the short and long preferred viscous waves are computed. The upward motion of the short wave together with the warm (cold) core lies to the west of the surface trough position above (below) 850 mb. Further, it is shown that the main source of kinetic energy for the wave lies in the middle layer (850–700 mb) which is transported to the lower and upper layers. Computed is found to be in good agreement with observed values.     

Lagrangian meshfree particles method (SPH) for large deformation and failure flows of geomaterial using elastic–plastic soil constitutive model

Simulation of large deformation and post‐failure of geomaterial in the framework of smoothed particle hydrodynamics (SPH) are presented in this study. The Drucker–Prager model with associated and non‐associated plastic flow rules is implemented into the SPH code to describe elastic–plastic soil behavior. In contrast to previous work on SPH for solids, where the hydrostatic pressure is often estimated from density by an equation of state, this study proposes to calculate the hydrostatic pressure of soil directly from constitutive models. Results obtained in this paper show that the original SPH method, which has been successfully applied to a vast range of problems, is unable to directly solve elastic–plastic flows of soil because of the so‐called SPH tensile instability. This numerical instability may result in unrealistic fracture and particles clustering in SPH simulation. For non‐cohesive soil, the instability is not serious and can be completely removed by using a tension cracking treatment from soil constitutive model and thereby give realistic soil behavior. However, the serious tensile instability that is found in SPH application for cohesive soil requires a special treatment to overcome this problem. In this paper, an artificial stress method is applied to remove the SPH numerical instability in cohesive soil. A number of numerical tests are carried out to check the capability of SPH in the current application. Numerical results are then compared with experimental and finite element method solutions. The good agreement obtained from these comparisons suggests that SPH can be extended to general geotechnical problems. Copyright © 2008 John Wiley & Sons, Ltd.     

中、低纬度地区电离层不同层结间的耦合对夜晚F区不规则结构的影响

电离层不规则结构的形成和演化与电离层等离子体不稳定性密切相关 .长期以来 ,在中低纬区扩展 -F现象的研究中 ,没有考虑电离层上下层结之间的相互作用 .本文从理论上全面探讨了在中纬度地区 ,E区可变的Pedersen电导率和Hall电导率 ,与F区可变的Pedersen电导率的共同作用对F区梯度漂移不稳定性的影响 ,导出了存在这种耦合时电离层梯度漂移不稳定性的统一表达式 .这一耦合理论不但解释了实际观测中发现的 ,在某些地方电离层F区顶部的不稳定性发生率要高于F区底部这一现象 ;同时还表明 ,电离层E区与F区的耦合对F区夜晚梯度漂移不稳定性的形成不仅会有阻碍作用 ,同时还使得中、低纬度地区的扰动增长具有了方向的选择性 .该理论的一个重要结论是 ,F区中某一地区能否发展出梯度漂移不稳定性 ,并不完全由当地电离层F区的状态决定 ,同一根磁力线连接的、位于不同纬度地区的E区层结对其发展和演化也会有相当大的影响 .