上地幔变黏度小尺度对流的数值研究

基于二维模型，利用有限元方法，研究上地幔-岩石圈系统的变黏度小尺度对流. 考虑该系统的黏度随温度以指数形式变化，数值结果表明，当黏度随温度变化较剧烈时，由于低温高黏度，系统的最上部物质不参与对流，系统发育形成一个类似于岩石圈的静止盖层. 计算表面热流、地形起伏及重力异常与对流格局有较好的相关性，高热流、上升地形对应于对流的上升区，反之低热流、下降地形与对流的下降区对应.     

利用地震层析成像数据计算地幔对流新模型的探讨

假设地幔地震层析成像数据对应的地幔横向不均匀结构是地幔热对流的结果. 将地震层析成像数据转化为地幔温度(或密度)不均匀分布，考虑热流体动力学的三个基本方程，顾及热输运方程中的非线性项，直接将地震层析成像转化的地幔温度不均匀分布作为内部荷载引入基本方程, 反演计算地幔对流. 本文在利用地震层析成像数据计算地幔对流模型的新理论和方法的基础上，用SH12WM13地震层析成像模型数据，计算了全球地幔对流格局. 结果表明，对流格局不仅依赖地震层析成像数据，而且在很大程度上受地幔动力学框架、热动力参数和边界条件的所确定的系统响应函数的影响. 显示了地幔中复杂的对流格局，特别是区域性层状对流以及多层对流环可能在地幔中存在的现象.     

三维重力内波和对流云街的启动机制

在一个下部为弱不稳定层、上部为稳定层的二层模式中，求解线性化的大气动力—热力学方程组，表明一个孤立山丘所引发的重力内波，在向下游传播过程中与边界层中对流相互作用，形成大气低层的涡漩活动，其图象与观测到的对流云街结构相似。     

Numerical models of mantle convection with secular cooling

A 2-D time-dependent finite-difference numerical model is used to investigate the thermal character and evolution of a convecting layer which is cooling as it convects. Two basic cooling modes are considered: in the first, both upper and lower boundaries are cooled at the same rate, while maintaining the same temperature difference across the layer; in the second, the lower boundary temperature decreases with time while the upper boundary temperature is fixed at 0°C. The first cooling mode simulates the effects of internal heating while the second simulates planetary cooling as mantle convection extracts heat from, and thereby cools, the Earth's core. The mathematical analogue between the effects of cooling and internal heating is verified for finite-amplitude convection. It is found that after an initial transient period the central core of a steady but vigorous convection cell cools at a constant rate which is governed by the rate of cooling of the boundaries and the viscosity structure of the layer. For upper-mantle models the transient stage lasts for about 30 per cent of the age of the Earth, while for the whole mantle it lasts for longer than the age of the Earth. Consequently, in our models the bulk cooling of the mantle lags behind the cooling of the core-mantle boundary. Models with temperature-dependent viscosity are found to cool in the same manner as models with depth-dependent viscosity; the rate of cooling is controlled primarily by the horizontally averaged variation of viscosity with depth. If the Earth's mantle cools in a similar fashion, secular cooling of the planet may be insensitive to lateral variations of viscosity.     

三维溶质运移问题的分步广义迎风解法

对对流占优的三维溶质运移问题提出了分步广义迎风解法，首先利用Ｎ．Ｎ．Ｙａ－ｎｅｎｋｏ对水动力弥散方程分步求解的思想，将原来的一个定解问题分解为两个定解问题即对流定解问题和扩散定解问题，对对流定解问题采用广义迎风对偶单元均衡法求解，对扩散定解问题采用一般的Ｇａｌｅｒｋｉｎ有限元法求解，不仅避免了用一般有限元法和有限差分法求解对流占优的地下水水质数学模型时常出现数值弥散和过量问题，而且避免了求节点速度这一步，节省运算步骤，对井点的浓度变化给出了更合适的求解方法。     

岩浆活动旋回与地球多层对流系统

从地球历史的节律角度讨论岩浆活动旋回。岩浆活动旋回可按时间与空间分出相互对应的不同尺度，进而可与深度尺度不同的地球多层对流系统相互对应，从而给予了岩浆活动旋回的动力学意义。岩浆活动的时间周期愈长，空间尺度亦愈大，需要的维持岩浆源存在的热能、挥发分供给的数量与深度亦愈大。     

Crystal morphologies in pillow basalts: implications for mid-ocean ridge processes

We present the results of a detailed petrological study of a sparsely phyric basalt (MAPCO CH98-DR11) dredged along the Mid-Atlantic Ridge (30°41′N). The sample contains microphenocrysts of olivine that display four different rapid-growth morphologies. Comparison of these morphologies with those obtained in dynamic crystallization experiments allows us to constrain the thermal history of the sample. The dendritic morphology (swallowtail, chain and lattice olivine) is directly related to the final quenching during magma–seawater interaction. In contrast, the three other morphologies, namely the complex polyhedral crystal, the closed hopper and the complex swallowtail morphology result from several cycles of cooling–heating (corresponding to a maximum degree of undercooling of 20–25°C) during crystal growth. These thermal variations occurred before eruption and are interpreted to be the result of turbulent convection in a small magmatic body beneath the ridge. The results suggest that the Mid-Atlantic Ridge is underlain by a mush zone that releases batches of liquid during tectonic segregation. Aphyric basalts are emitted during eruptions controlled by the tectonic activity, whereas phyric basalts correspond to small fractions of magma from the mush zone mobilized by reinjections of primitive magmas.     

Magmatic effects of the lunar late heavy bombardment

Large volumes of mare basalts are present on the surface of the moon, located preferentially in large impact basins. Mechanisms relating impact basins and mare basalt eruptions have previously been suggested: lunar impacts removed low-density material that may have inhibited eruption, and created cracks for fluid flow [Icarus 139 (1999) 246], and lunar basins have long been described as catchments for magma (e.g., [Rev. Geophys. Space Phys. 18 (1980) 107] and references therein). We present a new model for melt creation under near side lunar basins that is triggered by the impacts themselves. Magma can be produced in two stages. First, crater excavation depressurizes underlying material such that it may melt in-situ. Second, the cratered lithosphere rises isostatically, warping isotherms at the lithosphere-asthenosphere boundary which may initiate convection, in which adiabatic melting can occur. The first stage produces by far the largest volume of melt, but convective melting can continue for up to 350 Ma. We propose that giant impacts account for a large portion of the volume and longevity of mare basalt volcanism, as well as for several compositional groups, including high alumina, high titanium, KREEP-rich, and picritic magmas.     

The influence of thermochemical convection on the fixity of mantle plumes

A general feature of both isochemical and thermochemical studies of mantle convection is that horizontal plume velocities tend to be smaller than typical convective velocities, however, it is not clear which system leads to a greater fixity of mantle plumes. We perform two- and three-dimensional numerical calculations and compare both thermochemical and isochemical cases with similar convective vigor to determine whether presence of a dense component in the mantle can lead to smaller ratios of horizontal plume velocity to surface velocity. We investigate different viscosity and density contrasts between chemical components in the thermochemical calculations, and we perform isochemical calculations with both free-slip and no-slip bottom boundary conditions. We then compare both visually and quantitatively the results of the thermochemical and isochemical calculations to determine which leads to greater plume fixity. We find that horizontal plume velocities for thermochemical calculations are similar to those from isochemical calculations with no-slip bottom boundary conditions. In addition, we find that plumes tend to be more fixed for isochemical cases with free-slip bottom boundary conditions for two-dimensional calculations, however, in three dimensions, we find that plume fixity is similar to that observed in thermochemical calculations.