Using postglacial sea level,crustal velocities and gravity-rate-of-change to constrain the influence of thermal effects on mantle lateral heterogeneities

Lateral heterogeneities in the mantle can be caused by thermal, chemical and non-isotropic pre-stress effects. Here, we investigate the possibility of using observations of the glacial isostatic adjustment (GIA) process to constrain the thermal contribution to lateral variations in mantle viscosity. In particular, global historic relative sea level, GPS in Laurentide and Fennoscandia, altimetry together with tide-gauge data in the Great Lakes area, and GRACE data in Laurentide are used. The lateral viscosity perturbations are inferred from the seismic tomography model S20A by inserting the scaling factor β to determine the contribution of thermal effects versus compositional heterogeneity and non-isotropic pre-stress effects on lateral heterogeneity in mantle viscosity. When β = 1, lateral velocity variations are caused by thermal effects alone. With β < 1, the contribution of thermal effect decreases, so that for β = 0, there is no lateral viscosity variation and the Earth is laterally homogeneous. These lateral viscosity variations are superposed on four different reference models which differ significantly in the lower mantle viscosity. The Coupled Laplace Finite Element method is used to predict the GIA response on a spherical, self-gravitating, compressible, viscoelastic Earth with self-gravitating oceans, induced by the ICE-4G deglaciation model.Results show that the effect of β on uplift rates and gravity rate-of-change is not simple and involves the trade-off between the contribution of lateral viscosity variations in the transition zone and in the lower mantle. Models with small viscosity contrast in the lower mantle cannot explain the observed uplift rates in Laurentide and Fennoscandia. However, the RF3S20 model with a reference viscosity profile simplified from Peltier's VM2 with the value of β around 0.2–0.4 is found to explain most of the global RSL data, the uplift rates in Laurentide and Fennoscandia and the BIFROST horizontal velocity data. In addition, the changes in GIA signals caused by changes in the value of β are large enough to be detected by the data, although uncertainty in other parameters in the GIA models still exists. This may encourage us to further utilize GIA observations to constrain the thermal effect on mantle lateral heterogeneity as geodetic and satellite gravity measurements are improved.     

Effect of lateral viscosity variations in the core-mantle boundary region on predictions of the long-wavelength geoid

Seismic studies of the lowermost mantle suggest that the core-mantle boundary (CMB) region is strongly laterally heterogeneous over both local and global scales. These heterogeneities are likely to be associated with significant lateral viscosity variations that may influence the shape of the long-wavelength non-hydrostatic geoid. In the present paper we investigate the effect of these lateral viscosity variations on the solution of the inverse problem known as the inferences of viscosity from the geoid. We find that the presence of lateral viscosity variations in the CMB region can significantly improve the percentage fit of the predicted data with observations (from 42 to 70% in case of free-air gravity) while the basic characterisics of the mantle viscosity model, namely the viscosity increase with depth and the rate of layering, remain more or less the same as in the case of the best-fitting radially symmetric viscosity models. Assuming that viscosity is laterally dependent in the CMB region, and radially dependent elsewhere, we determine the largescale features of the viscosity structure in the lowermost mantle. The viscosity pattern found for the CMB region shows a high density of hotspots above the regions of higher-than-average viscosity. This result suggests an important role for petrological heterogeneities in the lowermost mantle, potentially associated with a post-perovskite phase transition. Another potential interpretation is that the lateral viscosity variations derived for the CMB region correspond in reality to lateral variations in the mechanical conditions at the CMB boundary or to large-scale undulations of a chemically distinct layer at the lowermost mantle.     

大地测量观测和相对海平面联合约束的冰川均衡调整模型

利用大地测量和历史相对海平面变化数据,结合地震剪切波层析模型,联合确定了新的末次冰期冰川均衡调整（GIA）模型,其中地幔黏滞度不仅沿径向而且沿横向变化.研究思路是,先尝试性地选择比例系数β,利用与地震剪切波速异常的线性关系,计算地幔黏滞度横向扰动,并与横向均匀的黏滞度参考模型叠加给出3D地幔黏滞度模型;再利用耦合拉普拉斯方程的有限元法算法进行GIA预测;然后,重复该过程,直到预测与观测之间的吻合满意为止.主要结论有：（1）给出了横向非均匀的地幔黏滞度模型（RF3L20(β=0.4)）,发现了黏滞度显著的横向非均匀性和其对GIA预测的显著影响,指出横向非均匀不完全是由热效应引起的,可能还与化学组分等其他因素有关,该模型可用于地幔动力学研究.（2）给出了全球现今多种GIA预测速率,可为板块运动、陆地水储量、海水质量变化和冰川冰雪质量非平衡监测提供重要的改正.     

GIA-induced secular variations in the Earth's long wavelength gravity field: Influence of 3-D viscosity variations

Predictions of present day secular variations in the Earth's long wavelength geopotential driven by glacial isostatic adjustment (GIA) have previously been analyzed to infer the radial profile of mantle viscosity and to constrain ongoing cryospheric mass balance. These predictions have been based on spherically symmetric Earth models. We explore the impact of lateral variations in mantle viscosity using a new finite-volume formulation for computing the response of 3-D Maxwell viscoelastic Earth models. The geometry of the viscosity field is constrained from seismic-to-mographic images of mantle structure, while the amplitude of the lateral viscosity variations is tuned by a free parameter in the modeling. We focus on the zonal J˙_? harmonics for degrees ? = 2,…,8 and demonstrate that large-scale lateral viscosity variations of two to three orders of magnitude have a modest, 5-10%, impact on predictions of J˙₂. In contrast, predictions of higher degree harmonics show a much greater sensitivity to lateral variation in viscosity structure. We conclude that future analyses of secular trends (for degree ? > 2) estimated from ongoing (GRACE, CHAMP) satellite missions must incorporate GIA predictions based on 3-D viscoelastic Earth models.     

Possible effects of lateral viscosity variations induced by plate-tectonic mechanism on geoid inferred from numerical models of mantle convection

In a traditional analytical method, the convective features of Earth’s mantle have been inferred from surface signatures obtained by the geodynamic model only with depth-dependent viscosity structure. The moving and subducting plates, however, bring lateral viscosity variations in the mantle. To clarify the effects of lateral viscosity variations caused by the plate-tectonic mechanism, I have first studied systematically instantaneous dynamic flow calculations using new density-viscosity models only with vertical viscosity variations in a three-dimensional spherical shell. I find that the geoid high arises over subduction zones only when the vertical viscosity contrast between the upper mantle and the lower mantle is O(10³) to O(10⁴), which seems to be much larger than the viscosity contrast suggested by other studies. I next show that this discrepancy may be removed when I consider the lateral viscosity variation caused by the plate-tectonic mechanism using two-dimensional numerical models of mantle convection with self-consistently moving and subducting plates, and suggest that the observed geoid anomaly on the Earth’s surface is significantly affected by plate-tectonic mechanism as a first-order effect.     

Using postglacial sea level, crustal velocities and gravity-rate-of-change to constrain the influence of thermal effects on mantle lateral heterogeneities

Lateral heterogeneities in the mantle can be caused by thermal, chemical and non-isotropic pre-stress effects. Here, we investigate the possibility of using observations of the glacial isostatic adjustment (GIA) process to constrain the thermal contribution to lateral variations in mantle viscosity. In particular, global historic relative sea level, GPS in Laurentide and Fennoscandia, altimetry together with tide-gauge data in the Great Lakes area, and GRACE data in Laurentide are used. The lateral viscosity perturbations are inferred from the seismic tomography model S20A by inserting the scaling factor β to determine the contribution of thermal effects versus compositional heterogeneity and non-isotropic pre-stress effects on lateral heterogeneity in mantle viscosity. When β = 1, lateral velocity variations are caused by thermal effects alone. With β < 1, the contribution of thermal effect decreases, so that for β = 0, there is no lateral viscosity variation and the Earth is laterally homogeneous. These lateral viscosity variations are superposed on four different reference models which differ significantly in the lower mantle viscosity. The Coupled Laplace Finite Element method is used to predict the GIA response on a spherical, self-gravitating, compressible, viscoelastic Earth with self-gravitating oceans, induced by the ICE-4G deglaciation model.Results show that the effect of β on uplift rates and gravity rate-of-change is not simple and involves the trade-off between the contribution of lateral viscosity variations in the transition zone and in the lower mantle. Models with small viscosity contrast in the lower mantle cannot explain the observed uplift rates in Laurentide and Fennoscandia. However, the RF3S20 model with a reference viscosity profile simplified from Peltier's VM2 with the value of β around 0.2–0.4 is found to explain most of the global RSL data, the uplift rates in Laurentide and Fennoscandia and the BIFROST horizontal velocity data. In addition, the changes in GIA signals caused by changes in the value of β are large enough to be detected by the data, although uncertainty in other parameters in the GIA models still exists. This may encourage us to further utilize GIA observations to constrain the thermal effect on mantle lateral heterogeneity as geodetic and satellite gravity measurements are improved.     

On the interaction between small- and large-scale convection and postglacial rebound flow in a power-law mantle

Some consequences arising from the superposition of flows of two different kinds or scales in a non-Newtonian mantle are discussed and applied to the cases mantle convection plus postglacial rebound flow as well as small- plus large-scale mantle convection. If the two flow types have similar magnitude, the apparent rheology of both flows becomes anisotropic and the apparent viscosity for one flow depends on the geometry of the other. If one flow has a magnitude significantly larger than the other, the apparent viscosity for the weak flow is linear but develops direction-dependent variations about a factorn (n being the power exponent of the rheology). For the rebound flow lateral variations of the apparent viscosity about at least 3 are predicted and changes in the flow geometry and relaxation time are possible. On the other hand, rebound flow may weaken the apparent viscosity for convection. Secondary convection under moving plates may be influenced by the apparent anisotropic rheology. Other mechanisms leading to viscous anisotropy during shearing may increase this effect. A linear stability analysis for the onset of convection with anisotropic linear rheology shows that the critical Rayleigh number decreases and the aspect ratio of the movement cells increases for decreasing horizontal shear viscosity (normal viscosity held constant). Applied to the mantle, this model weakens the preference of convection rolls along the direction of plate motion. Under slowly moving plates, rolls perpendicular to the plate motion seem to have a slight preference. These results could be useful for resolving the question of Newtonian versus non-Newtonian or isotropic versus anisotropic mantle rheology.     

New Perspectives on Mantle Dynamics from High-resolution Seismic Tomographic Model P1200

—Recently a high-resolution tomographic model, the P1200, based on P-wave travel times was developed, which allowed for detailed imaging of the top 1200 km of the mantle. This model was used in diverse ways to study mantle viscosity structure and geodynamical processes. In the spatial domain there are lateral variations in the transition zone, suggesting interaction between the lower-mantle plumes and the region from 600 km to 1000 km. Some examples shown here include the continental region underneath Manchuria, Ukraine and South Africa, where horizontal structures lie above or below the 660 km discontinuity. The blockage of upwelling is observed under central Africa and the interaction between the upwelling and the transition zone under the slow Icelandic region appears to be complex. An expansion of the aspherical seismic velocities has been taken out to spherical harmonics of degree 60. For degrees exceeding around 10, the spectra at various depths decay with a power-law like dependence on the degree, with the logarithmic slopes in the asymptotic portion of the spectra containing values between 2 and 2.6. These spectral results may suggest the time-dependent nature of mantle convection. Details of the viscosity structure in the top 1200 km of the mantle have been inferred both from global and regional geoid data and from the high-resolution tomographic model. We have considered only the intermediate degrees (l = 12–25) in the nonlinear inversion with a genetic algorithm approach. Several families of acceptable viscosity profiles are found for both oceanic and global data. The families of solutions for the two data sets have different characteristics. Most of the solutions asociated with the global geoid data show the presence of asthenosphere below the lithosphere. In other families a low viscosity zone between 400 and 600 km depth is found to lie atop a viscosity jump. Other families evidence a viscosity decrease across the 660 km discontinuity. Solutions from oceanic geoid show basically two low viscosity zones one lying right below the lithosphere; the other right under 660-km depth. All of these results bespeak clearly the plausible existence of strong vertical viscosity stratification in the top 1000 km of the mantle. The presence of the second asthenosphere may have important dynamical ramifications on issues pertaining to layered mantle convection. Numerical modelling of mantle convection with two phase transitions and a realistic temperature- and pressure-dependent viscosity demonstrates that a low viscosity region under the endothermic phase transition can indeed be generated self-consistently in time-dependent situations involving a partially layered configuration in an axisymmetric spherical-shell model.     

On a spectral method of solving the Stokes equation

A method of solving the Stokes equation for a spherical mantle model by expansion in spherical harmonics was developed by Hager and O’Connell [1979]. However, this method is applicable only if the viscosity depends solely on depth. In this case, the Stokes equation reduces to a system of independent equations for each harmonic. Given lateral variations in viscosity, the Stokes equation contains terms in the form of products of harmonics, which invalidates all advantages of harmonic expansion. Zhang and Christensen [1993] developed a perturbation method for the case when terms containing products of lateral viscosity variations are small. These terms are first calculated from the preceding iteration and are then expanded in a series of harmonic functions. As a result, equations for harmonics remain independent. An evident advantage of the spectral method is the simplicity of the technique of incorporating the self-gravitation and compressibility effects. Moreover, this method partially eliminates difficulties related to the singularities at poles. As yet, it has not been applied in practice, possibly because the equations presented in [Zhang and Christensen, 1993] contain misprints that have not been elucidated in the literature. In the present work, a system of equations is derived for the spectral-iterative method of solving the Stokes equation and the errata present in formulas of Zhang and Christensen [1993] and significantly affecting results of calculations are analyzed.     

Viscosity distribution in the mantle convection models

Viscosity is a fundamental property of the mantle which determines the global geodynamical processes. According to the microscopic theory of defects and laboratory experiments, viscosity exponentially depends on temperature and pressure, with activation energy and activation volume being the parameters. The existing laboratory measurements are conducted with much higher strain rates than in the mantle and have significant uncertainty. The data on postglacial rebound only allow the depth distributions of viscosity to be reconstructed. Therefore, spatial distributions (along the depth and lateral) are as of now determined from the models of mantle convection which are calculated by the numerical solution of the convection equations, together with the viscosity dependences on pressure and temperature (PT-dependences). The PT-dependences of viscosity which are presently used in the numerical modeling of convection give a large scatter in the estimates for the lower mantle, which reaches several orders of magnitude. In this paper, it is shown that it is possible to achieve agreement between the calculated depth distributions of viscosity throughout the entire mantle and the postglacial rebound data. For this purpose, the values of the volume and energy of activation for the upper mantle can be taken from the laboratory experiments, and for the lower mantle, the activation volume should be reduced twice at the 660-km phase transition boundary. Next, the reduction in viscosity by an order of magnitude revealed at the depths below 2000 km by the postglacial rebound data can be accounted for by the presence of heavy hot material at the mantle bottom in the LLSVP zones. The models of viscosity spatial distribution throughout the entire mantle with the lithospheric plates are presented.     

A preliminary study of the effects of some flow parameters in the generation of poloidal and toroidal energies within a 3-D spherical thermal-convective system with variable viscosity

The effects of variable viscosity on flow dynamics within spherical shells are investigated using a finite-element thermal convection model, and preliminary result for cases with relatively low Rayleigh numbers and small viscosity contrasts are reported. These results demonstrate some general effects of viscosity variation on mantle dynamics, and, in particular, the generation of toroidal energy. Since lateral viscosity variations are necessary in the generation of toroidal motion in a thermally driven convective system, it is not surprising our results show that flows with greater viscosity contrasts produce greater amounts of toroidal energy. Our preliminary study further shows that solutions become more time-dependent as viscosity contrasts increase. Increasing the Rayleigh number is also found to increase the magnitude of toroidal energy. Internal heating, on the other hand, appears to lead to less toroidal energy compared wth bottom heating because it tends to produce a thermally more uniform interior and thus smaller viscosity variations.     

On asymptotics of low-frequency oscillations of the liquid core: 2. Frequency dependence of the inertial core-mantle coupling

New, unique information on the inertial and dissipative coupling of the liquid core and the mantle has been retrieved from modern high-precision (radiointerferometer and GPS) data on tidal variations in the rotation velocity and nutation of the Earth. Comparison of theoretical and observed data provided new estimates for the dynamic flattening of the outer liquid and the inner solid cores, mantle quality factor, viscosity of the liquid core, and electromagnetic coupling of the liquid core and the mantle [Molodensky, 2004, 2006]. As was shown in the first part of the paper [Molodensky, 2008] (further referred to as [I]), generation of eddy flows in Proudman-Taylor columns, whose orientation is controlled by the topography of the liquid core-mantle boundary, should be taken into account for correct estimation of the inertial coupling (see formulas (8) and (34) in [I]). The range of periods within which this effect plays a significant role is determined by the decay time of these flows. This time is estimated in the paper for the case where dissipation is related to viscous friction at the core-mantle boundary or with the electromagnetic coupling of the liquid core and the mantle. Because of significant uncertainties in modern data on the viscosity of the liquid core, the magnetic field intensity at the core-mantle boundary, and the electrical conductivity of the lower mantle, the dissipative coupling of the liquid core and the mantle cannot be calculated as yet. However, as shown in the paper, the decay time of eddy flows is connected with the attenuation time of subdiurnal free nutation and with the liquid core viscosity. This enables the estimation of the frequency dependence of the dissipative coupling in a fairly wide range. It is shown that the range of periods for which relations (8) and (34) in [I] are valid encompasses the best-studied length-of-day variations and, therefore, these relations are applicable to analysis of the majority of modern data.     

Mantle rheology: radial and lateral viscosity variations inferred from microphysical creep laws

Inferences on the rheology of the mantle based on theoretical and experimental rate equations for steady state creep are discussed and compared with results from geophysical models. The radial increase of viscosity by one to three orders of magnitude across the mantle, required by inversion of postglacial rebound and geodynamic data, is confirmed by microphysical models based on the estimation of continuous and discontinuous changes of creep parameters with depth. The upper mantle (viscosity 10²⁰–10²¹ Pa s) is likely to show non-Newtonian rheology (power-law creep) for average grain sizes larger than 0.1 mm as an order of magnitude. Given the variability of both grain size and stress conditions, local regions of linear rheology can be present. The rheology of transition zone and lower mantle (viscosity 10²²–10²⁴ Pa s) cannot be definitely resolved at present. Estimation of creep parameters leads to possible nonlinear or mixed rheology, if grain sizes are not lower than 0.1 mm and flow conditions can be approximated by a constant strain rate of about 10⁻¹⁵ s⁻¹. This conclusion can be modified by different flow conditions (e.g. a decrease in strain rate or constant viscous dissipation). Furthermore, experiments on fine-grained garnetites and perovskite analogues have shown that diffusion creep is predominant at laboratory conditions. However, the pressure dependence of creep in these phases is unknown, and therefore direct extrapolation to lower mantle conditions is necessarily speculative. Lateral variations of viscosity, largest in the upper and lowermost mantle (up to 2–4 orders of magnitude) are predicted by models based on lateral temperature anomalies derived from seismic tomographic models.     

Estimation of the deep density distribution in the lithosphere of central and Southern Asia using data about free mantle surface depth

The results of the selection of a model of the deep density distribution in the lithosphere of Central and Southern Asia, which explains the previously revealed dependence of the free mantle surface depth on the thickness of the crust [Artemjev, 1975], are described. It is shown that this dependence can be caused by variations in the mantle’s density with depth. Models of the continental and oceanic mantles with an increase in the linear density over depth are selected for the region of Asia. The level of the free surface depth in the oceanic mantle is higher than in the continental mantle. The observed dependence on the crustal thickness can also be used for determining nonlinear density variations with depth under the assumption that lateral density variations in this dependence are of a random character.     

Nonlinear free thermal convection in a spherical shell:A variable viscosity model

Introduction The velocity field of surface plate motion can be split into a poloidal and a toroidal parts.At the Earth′s surface,the toroidal component is manifested by the existence of transform faults,and the poloidal component by the presence of convergence and divergence,i.e.spreading and subduc-tion zones.They have coupled each other and completely depicted the characteristics of plate tec-tonic motions.The mechanism of poloidal field has been studied fairly clearly which is related to …     

Geoid and topography of Venus in various thermal convection models

Important though indirect information about the internal structure of Venus is provided by its topography and geoid. In the last decades this information has been used to constrain the Venus mantle viscosity structure and its dynamic regime. Recently, the geodynamic inversion of the Venus?? geoid and topography resulted in a group of best fitting viscosity profiles. We use these viscosity models here as an input to our mantle convection code. We carry out simulations of the Venus?? mantle evolution in a 3D spherical shell with depth dependent viscosity and check whether the character of the dynamic topography and the geoid represented by their power spectra fits the observed quantities. We compare the results with several other models obtained for different viscosity stratifications (constant, constant with highly viscous lithosphere, linear increase of viscosity). Further, we estimate the effect of other factors such as internal heating and varying Rayleigh number. We use a 2D spherical axisymmetric convection code to study the effect of lateral viscosity variations. In these 2D models we monitor the topography and the geoid developing above the axisymmetric plume and compare them with the observed elevations of Venus?? geoid and topography in several Regia. Though none of the models fits observed data perfectly, we can generally conclude, that the best fit between the observed and predicted quantities is reached for viscosity profiles with 200 km thick lithosphere followed by a gradual increase of viscosity with depth and with the upper mantle viscosity of 2 × 10 ²¹ Pa s. For all viscosity profiles the predicted geoid and topography spectra match the observed ones only up to the degree 40, thus indicating other than dynamic origin of these quantities for higher degrees.     

Magnetic constraints on the spatial distribution of seismic activity

The lithospheric magnetic field (LMF) in China and its surrounding are calculated using the spherical harmonic coefficients given by the NGDC-720 model. The LMF comes from the magnetization of minerals in the crust and in the uppermost mantle. It may, therefore, provide unique insight into lithospheric tectonic processes and mechanisms. Here, we study the geomagnetic manifestation of active tectonic blocks, and find a close correlation between the LMF and seismicity. Many large faults are found to closely overlap with magnetic anomalies, or are distributed along the boundaries of magnetic anomalies. Earthquakes in these fault regions have occurred on the boundaries of magnetic anomalies, or in the transition zones between positive and negative anomalies. We analyze the components of the LMF, and the LMFs at different altitudes, finding that the vertical component, B _z at 200 km, is the most related to seismic activity. Relevant physical mechanisms are also discussed. We propose that the stress or viscosity differences caused by temperature variations, which manifest in the LMF, may be the predominant reason for the correlation between the LMF and seismic activity along large faults.     

Short scale heterogeneity in the lowermost mantle: insights from PcP-P and ScS-S data

We compare lateral variations at the base of the mantle as inferred from a global dataset of PcP-P travel time residuals, measured on broadband records, and existing P and S tomographic velocity models, as well as ScS-S travel time data in some selected regions. In many regions, the PcP-P dataset implies short scale lateral variations that are not resolved by global tomographic models, except under eastern Eurasia, where data and models describe a broad region of fast velocity anomalies across which variations appear to be of thermal origin. In other regions, such as central America and southeastern Africa, correlated short scale lateral variations (several hundred kilometers) are observed in PcP and ScS, implying large but not excessive values for the ratio R=∂ ln V_s/∂ ln V_p (∼2.5). On the other hand, in at least two instances, in the heart of the African Plume and on the edge of the Pacific Plume, variations in P and S velocities appear to be incompatible, implying strong lateral gradients across compositionally different domains, possibly also involving topography on the core-mantle boundary. One should be cautious in estimating R at the base of the mantle from global datasets, as different smoothing and sampling of P and S datasets may result in strong biases and meaningless results.     

Influence of variable uncertainties in seismic tomography models on constraining mantle viscosity from geoid observations

The radial viscosity structure of the Earth is explored on the basis of the geoid observations. The variations of uncertainty in seismic tomography models are accounted for when finding the radial viscosity structure. The new methodology we propose attempts to fit more closely those features of the geoid that are better constrained by tomography models and avoids to fit those features that are poorly constrained. This approach is particularly important because the error of geoid predictions caused by uncertainties in seismic tomography models is overwhelmingly larger than the noise in the geoid measurements. The synthetic tests indicate that the viscosity structures obtained by disregarding the uncertainty variations in seismic tomography models can be biased depending on the geoid spectral band and on the ‘input’ seismic tomography model. When the uncertainty variations in seismic models are considered in the inversion process, results do not indicate a viscosity in the transition zone lower than in the upper mantle. A robust feature found with the new method is a viscosity in the upper mantle two orders of magnitude smaller than in the lower mantle. The error covariance of seismic tomography models is critical for the method we suggest. A covariance matrix rigorously derived by seismologists should help to even more reliably infer the viscosity structure and relation between anomalies in density and seismic velocities from surface observations such as the geoid, and thus lead to a better knowledge of the Earth interior.     

地幔对流环型场的激发(Ⅰ)

讨论了地幔内部的粘滞度及施加在地表和CMB 的边界条件对地幔对流环型场的激发分析表明,当粘滞度侧向均匀时,环型场与极型场自然解耦,且环型场不影响重力位,当粘滞度侧向不均时,环型场与极型场耦合在一起.两者共同影响重力位.当引入板块运动速度时,边界条件非零,也能激发环型场;对侧向均匀粘滞度地幔,零边界条件不能激发环型场