地幔柱的发生机制研究——以Ferrar为例

利用三维数值模拟方法,根据地幔对流控制方程,以180 Ma喷发的Ferrar大火成岩省(LIP)为例,在模型中引入Pangea超大陆、大型横波低速带(LLSVPs)和Pangea超大陆边缘的俯冲带,模拟地幔对流过程,研究其对应地幔热柱从地球内部热边界层(例如核幔边界)的生成过程,并讨论导致该地幔柱产生的相关因素.结果表...     

Multiple volcanic episodes of flood basalts caused by thermochemical mantle plumes

The hypothesis that a single mushroom-like mantle plume head can generate a large igneous province within a few million years has been widely accepted. The Siberian Traps at the Permian-Triassic boundary and the Deccan Traps at the Cretaceous-Tertiary boundary were probably erupted within one million years. These large eruptions have been linked to mass extinctions. But recent geochronological data reveal more than one pulse of major eruptions with diverse magma flux within several flood basalts extending over tens of million years. This observation indicates that the processes leading to large igneous provinces are more complicated than the purely thermal, single-stage plume model suggests. Here we present numerical experiments to demonstrate that the entrainment of a dense eclogite-derived material at the base of the mantle by thermal plumes can develop secondary instabilities due to the interaction between thermal and compositional buoyancy forces. The characteristic timescales of the development of the secondary instabilities and the variation of the plume strength are compatible with the observations. Such a process may contribute to multiple episodes of large igneous provinces.     

High 3He/4He ratios in picritic basalts from Baffin Island and the role of a mixed reservoir in mantle plumes

The high 3He/4He ratio of volcanic rocks thought to be derived from mantle plumes is taken as evidence for the existence of a mantle reservoir that has remained largely undegassed since the Earth's accretion. The helium isotope composition of this reservoir places constraints on the origin of volatiles within the Earth and on the evolution and structure of the Earth's mantle. Here we show that olivine phenocrysts in picritic basalts presumably derived from the proto-Iceland plume at Baffin Island, Canada, have the highest magmatic 3He/4He ratios yet recorded. A strong correlation between 3He/4He and 87Sr/86Sr, 143Nd/144Nd and trace element ratios demonstrate that the 3He-rich end-member is present in basalts that are derived from large-volume melts of depleted upper-mantle rocks. This reservoir is consistent with the recharging of depleted upper-mantle rocks by small volumes of primordial volatile-rich lower-mantle material at a thermal boundary layer between convectively isolated reservoirs. The highest 3He/4He basalts from Hawaii and Iceland plot on the observed mixing trend. This indicates that a 3He-recharged depleted mantle (HRDM) reservoir may be the principal source of high 3He/4He in mantle plumes, and may explain why the helium concentration of the 'plume' component in ocean island basalts is lower than that predicted for a two-layer, steady-state model of mantle structure.     

Whole-mantle convection and the transition-zone water filter

Because of their distinct chemical signatures, ocean-island and mid-ocean-ridge basalts are traditionally inferred to arise from separate, isolated reservoirs in the Earth's mantle. Such mantle reservoir models, however, typically satisfy geochemical constraints, but not geophysical observations. Here we propose an alternative hypothesis that, rather than being divided into isolated reservoirs, the mantle is filtered at the 410-km-deep discontinuity. We propose that, as the ascending ambient mantle (forced up by the downward flux of subducting slabs) rises out of the high-water-solubility transition zone (between the 660 km and 410 km discontinuities) into the low-solubility upper mantle above 410 km, it undergoes dehydration-induced partial melting that filters out incompatible elements. The filtered, dry and depleted solid phase continues to rise to become the source material for mid-ocean-ridge basalts. The wet, enriched melt residue may be denser than the surrounding solid and accordingly trapped at the 410 km boundary until slab entrainment returns it to the deeper mantle. The filter could be suppressed for both mantle plumes (which therefore generate wetter and more enriched ocean-island basalts) as well as the hotter Archaean mantle (thereby allowing for early production of enriched continental crust). We propose that the transition-zone water-filter model can explain many geochemical observations while avoiding the major pitfalls of invoking isolated mantle reservoirs.     

Seismological constraints on a possible plume root at the core-mantle boundary

Recent seismological discoveries have indicated that the Earth's core-mantle boundary is far more complex than a simple boundary between the molten outer core and the silicate mantle. Instead, its structural complexities probably rival those of the Earth's crust. Some regions of the lowermost mantle have been observed to have seismic wave speed reductions of at least 10 per cent, which appear not to be global in extent. Here we present robust evidence for an 8.5-km-thick and approximately 50-km-wide pocket of dense, partially molten material at the core-mantle boundary east of Australia. Array analyses of an anomalous precursor to the reflected seismic wave ScP reveal compressional and shear-wave velocity reductions of 8 and 25 per cent, respectively, and a 10 per cent increase in density of the partially molten aggregate. Seismological data are incompatible with a basal layer composed of pure melt, and thus require a mechanism to prevent downward percolation of dense melt within the layer. This may be possible by trapping of melt by cumulus crystal growth following melt drainage from an anomalously hot overlying region of the lowermost mantle. This magmatic evolution and the resulting cumulate structure seem to be associated with overlying thermal instabilities, and thus may mark a root zone of an upwelling plume.     

Preservation of ancient and fertile lithospheric mantle beneath the southwestern United States

Stable continental regions, free from tectonic activity, are generally found only within ancient cratons-the centres of continents which formed in the Archaean era, 4.0-2.5 Gyr ago. But in the Cordilleran mountain belt of western North America some younger (middle Proterozoic) regions have remained stable, whereas some older (late Archaean) regions have been tectonically disturbed, suggesting that age alone does not determine lithospheric strength and crustal stability. Here we report rhenium-osmium isotope and mineral compositions of peridotite xenoliths from two regions of the Cordilleran mountain belt. We found that the younger, undeformed Colorado plateau is underlain by lithospheric mantle that is 'depleted' (deficient in minerals extracted by partial melting of the rock), whereas the older (Archaean), yet deformed, southern Basin and Range province is underlain by 'fertile' lithospheric mantle (not depleted by melt extraction). We suggest that the apparent relationship between composition and lithospheric strength, inferred from different degrees of crustal deformation, occurs because depleted mantle is intrinsically less dense than fertile mantle (due to iron having been lost when melt was extracted from the rock). This allows the depleted mantle to form a thicker thermal boundary layer between the deep convecting mantle and the crust, thus reducing tectonic activity at the surface. The inference that not all Archaean crust developed a strong and thick thermal boundary layer leads to the possibility that such ancient crust may have been overlooked because of its intensive reworking or lost from the geological record owing to preferential recycling.     

Water-rich basalts at mid-ocean-ridge cold spots

Although water is only present in trace amounts in the suboceanic upper mantle, it is thought to play a significant role in affecting mantle viscosity, melting and the generation of crust at mid-ocean ridges. The concentration of water in oceanic basalts has been observed to stay below 0.2 wt%, except for water-rich basalts sampled near hotspots and generated by 'wet' mantle plumes. Here, however, we report unusually high water content in basaltic glasses from a cold region of the mid-ocean-ridge system in the equatorial Atlantic Ocean. These basalts are sodium-rich, having been generated by low degrees of melting of the mantle, and contain unusually high ratios of light versus heavy rare-earth elements, implying the presence of garnet in the melting region. We infer that water-rich basalts from such regions of thermal minima derive from low degrees of 'wet' melting greater than 60 km deep in the mantle, with minor dilution by melts produced by shallower 'dry' melting--a view supported by numerical modelling. We therefore conclude that oceanic basalts are water-rich not only near hotspots, but also at 'cold spots'.     

Transient high temperatures in mantle plume heads inferred from magnesian olivines in Phanerozoic picrites

Both scaled laboratory experiments and numerical models of terrestrial mantle plumes produce 'balloon-on-a-string' structures, with a bulbous head followed by a stem-like tail. Discussions have focused on whether their initial upwelling heads are hotter than the tails or cooler, as a result of entrainment of ambient mantle during ascent, and also on whether initial plume upwelling is a newtonian or non-newtonian process. The temperature of the mantle delivered to the base of the lithosphere is a critical parameter in such debates. Dry continental magmas can normally contribute little to this topic because their hottest (ultramafic) examples can be expected to be trapped, owing to their density, beneath the Moho. Here we report a rare case in which olivine (with 93.3% forsterite; Mg2SiO4) phenocrysts, precipitated from an unerupted komatiitic melt (approximately 24% MgO) of the Tristan mantle plume head 132 Myr ago, were carried to upper-crust levels in northwest Namibia by less Mg-rich (9.6-18.5% MgO) magmas. We infer that the hidden melt, generated when the plume impinged on the base of the lithosphere, originated in the mantle with a potential temperature of approximately 1,700 degrees C. This is approximately 400 degrees C above ambient and much hotter than the temperatures previously calculated for steady-state Phanerozoic mantle plumes. Published data show that the same conclusion can be reached for the initial Iceland and Galapagos plumes.     

A case for mantle plumes

The existence of at least several plumes in the Earth＇s mantle can be inferred with few assumptions from well-established observations. As well, thermal mantle plumes can be predicted from well-established and quantified fluid dynamics and a plausible assumption about the Earth＇s early thermal state. Some additional important observations, especially of flood basalts and rift-related magmatism, have been shown to be plausibly consistent with the physical theory. Recent claims to have detected plumes using seismic tomography may comprise the most direct evidence for plumes, but plume tails are likely to be difficult to resolve definitively and the claims need to be well tested. Although significant questions remain about its viability, the plume hypothesis thus seems to be well worth continued investigation. Nevertheless there are many non-plate-related magmatic phenomena whose association with plumes is unclear or unlikely. Compositional buoyancy has recently been shown potentially to substantially complicate the dynamics of plumes, and this may lead to explanations for a wider range of phenomena, including “headless” hotspot tracks, than purely thermal plumes.     

Insights into the dynamics of mantle plumes from uranium-series geochemistry

The long-standing paradigm that hotspot volcanoes such as Hawaii or Iceland represent the surface expression of mantle plumes--hot, buoyant upwelling regions beneath the Earth's lithosphere--has recently been the focus of controversy. Whether mantle plumes exist or not is pivotal for our understanding of the thermal, dynamic and compositional evolution of the Earth's mantle. Here we show that uranium-series disequilibria measured in hotspot lavas indicate that hotspots are indeed associated with hot and buoyant upwellings and that weaker (low buoyancy flux) hotspots such as Iceland and the Azores are characterized by lower excess temperatures than stronger hotspots such as Hawaii. This direct link between buoyancy flux and mantle temperature is evidence for the existence of mantle plumes.     

临界理查森数对不稳定状态下边界层模拟的影响

在边界层参数化方案中,临界理查森数可以用来判断湍流状态和计算边界层高度。大部分模式会将临界理查森数设置为常数,但实际上它会发生改变。本文基于WRF（Weather Research and Forecasting）中尺度数值预报模式和YSU（Yonsei University）边界层参数化方案（以下简称WRF/YSU）,以一次夏季的晴好天气过程为背景,研究了当临界理查森数发生变化,会对气象场和边界层湍流产生的影响。结果表明：（1）当临界理查森数增加时,大范围区域的边界层的高度会增加。在时空上能够体现出来。在空间上,边界层高度分布不均。从时间上看,在下午边界层高度受到的影响最大;（2）理查森数的临界值增加,湍流扩散系数也增大,水汽混合比增加,并且地面热通量总体增大,又会进一步使边界层高度增大;（3）YSU方案中湍流扩散方程由一个局部项、一个非局部项和一个夹卷项组成。对于热通量,夹卷项对总热通量的贡献与局部项和非局部项的贡献相当。对于水汽和动量通量,局地项及夹卷通量项在边界层中最大。     

Deep-mantle high-viscosity flow and thermochemical structure inferred from seismic and geodynamic data

Surface geophysical data that are related to the process of thermal convection in the Earth's mantle provide constraints on the rheological properties and density structure of the mantle. We show that these convection-related data imply the existence of a region of very high effective viscosity near 2,000 km depth. This inference is obtained using a viscous-flow model based on recent high-resolution seismic models of three-dimensional structure in the mantle. The high-viscosity layer near 2,000 km depth results in a re-organization of flow from short to long horizontal length scales, which agrees with seismic tomographic observations of very long wavelength structures in the deep mantle. The high-viscosity region also strongly suppresses flow-induced deformation and convective mixing in the deep mantle. Here we predict compositional and thermal heterogeneity in this region, using viscous-flow calculations based on the new viscosity profile, together with independent mineral physics data. These maps are consistent with the anti-correlation of anomalies in seismic shear and bulk sound velocity in the deep mantle. The maps also show that mega-plumes in the lower mantle below the central Pacific and Africa are, despite the presence of compositional heterogeneity, buoyant and actively upwelling structures.     

Spin crossover and iron-rich silicate melt in the Earth's deep mantle

A melt has greater volume than a silicate solid of the same composition. But this difference diminishes at high pressure, and the possibility that a melt sufficiently enriched in the heavy element iron might then become more dense than solids at the pressures in the interior of the Earth (and other terrestrial bodies) has long been a source of considerable speculation. The occurrence of such dense silicate melts in the Earth's lowermost mantle would carry important consequences for its physical and chemical evolution and could provide a unifying model for explaining a variety of observed features in the core-mantle boundary region. Recent theoretical calculations combined with estimates of iron partitioning between (Mg,Fe)SiO(3) perovskite and melt at shallower mantle conditions suggest that melt is more dense than solids at pressures in the Earth's deepest mantle, consistent with analysis of shockwave experiments. Here we extend measurements of iron partitioning over the entire mantle pressure range, and find a precipitous change at pressures greater than ～76?GPa, resulting in strong iron enrichment in melts. Additional X-ray emission spectroscopy measurements on (Mg(0.95)Fe(0.05))SiO(3) glass indicate a spin collapse around 70?GPa, suggesting that the observed change in iron partitioning could be explained by a spin crossover of iron (from high-spin to low-spin) in silicate melt. These results imply that (Mg,Fe)SiO(3) liquid becomes more dense than coexisting solid at ～1,800?km depth in the lower mantle. Soon after the Earth's formation, the heat dissipated by accretion and internal differentiation could have produced a dense melt layer up to ～1,000?km in thickness underneath the solid mantle. We also infer that (Mg,Fe)SiO(3) perovskite is on the liquidus at deep mantle conditions, and predict that fractional crystallization of dense magma would have evolved towards an iron-rich and silicon-poor composition, consistent with seismic inferences of structures in the core-mantle boundary region.     

Evidence against a chondritic Earth

The (142)Nd/(144)Nd ratio of the Earth is greater than the solar ratio as inferred from chondritic meteorites, which challenges a fundamental assumption of modern geochemistry--that the composition of the silicate Earth is 'chondritic', meaning that it has refractory element ratios identical to those found in chondrites. The popular explanation for this and other paradoxes of mantle geochemistry, a hidden layer deep in the mantle enriched in incompatible elements, is inconsistent with the heat flux carried by mantle plumes. Either the matter from which the Earth formed was not chondritic, or the Earth has lost matter by collisional erosion in the later stages of planet formation.     

Anisotropy of thermal diffusivity in the upper mantle.

Heat transfer in the mantle is a key process controlling the Earth's dynamics. Upper-mantle mineral phases, especially olivine, have been shown to display highly anisotropic thermal diffusivity at ambient conditions, and seismic anisotropy data show that preferred orientations of olivine induced by deformation are coherent at large scales (>50 km) in the upper mantle. Thus heat transport in the upper mantle should be anisotropic. But the thermal anisotropy of mantle minerals at high temperature and its relationship with deformation have not been well constrained. Here we present petrophysical modelling and laboratory measurements of thermal diffusivity in deformed mantle rocks between temperatures of 290 and 1,250 K that demonstrate that deformation may induce a significant anisotropy of thermal diffusivity in the uppermost mantle. We found that heat transport parallel to the flow direction is up to 30 per cent faster than that normal to the flow plane. Such a strain-induced thermal anisotropy implies that the upper-mantle temperature distribution, rheology and, consequently, its dynamics, will depend on deformation history. In oceans, resistive drag flow would result in lower vertical diffusivities in both the lithosphere and asthenosphere and hence in less effective heat transfer from the convective mantle. In continents, olivine orientations frozen in the lithosphere may induce anisotropic heating above mantle plumes, favouring the reactivation of pre-existing structures.     

A crystallizing dense magma ocean at the base of the Earth's mantle

The distribution of geochemical species in the Earth's interior is largely controlled by fractional melting and crystallization processes that are intimately linked to the thermal state and evolution of the mantle. The existence of patches of dense partial melt at the base of the Earth's mantle, together with estimates of melting temperatures for deep mantle phases and the amount of cooling of the underlying core required to maintain a geodynamo throughout much of the Earth's history, suggest that more extensive deep melting occurred in the past. Here we show that a stable layer of dense melt formed at the base of the mantle early in the Earth's history would have undergone slow fractional crystallization, and would be an ideal candidate for an unsampled geochemical reservoir hosting a variety of incompatible species (most notably the missing budget of heat-producing elements) for an initial basal magma ocean thickness of about 1,000 km. Differences in 142Nd/144Nd ratios between chondrites and terrestrial rocks can be explained by fractional crystallization with a decay timescale of the order of 1 Gyr. These combined constraints yield thermal evolution models in which radiogenic heat production and latent heat exchange prevent early cooling of the core and possibly delay the onset of the geodynamo to 3.4-4 Gyr ago.     

A doubling of the post-perovskite phase boundary and structure of the Earth's lowermost mantle

The thermal structure of the Earth's lowermost mantle--the D' layer spanning depths of approximately 2,600-2,900 kilometres--is key to understanding the dynamical state and history of our planet. Earth's temperature profile (the geotherm) is mostly constrained by phase transitions, such as freezing at the inner-core boundary or changes in crystal structure within the solid mantle, that are detected as discontinuities in seismic wave speed and for which the pressure and temperature conditions can be constrained by experiment and theory. A recently discovered phase transition at pressures of the D' layer is ideally situated to reveal the thermal structure of the lowermost mantle, where no phase transitions were previously known to exist. Here we show that a pair of seismic discontinuities observed in some regions of D' can be explained by the same phase transition as the result of a double-crossing of the phase boundary by the geotherm at two different depths. This simple model can also explain why a seismic discontinuity is not observed in some other regions, and provides new constraints for the magnitude of temperature variations within D'.     

An early lunar core dynamo driven by thermochemical mantle convection

Although the Moon currently has no internally generated magnetic field, palaeomagnetic data, combined with radiometric ages of Apollo samples, provide evidence for such a magnetic field from approximately 3.9 to 3.6 billion years (Gyr) ago, possibly owing to an ancient lunar dynamo. But the presence of a lunar dynamo during this time period is difficult to explain, because thermal evolution models for the Moon yield insufficient core heat flux to power a dynamo after approximately 4.2 Gyr ago. Here we show that a transient increase in core heat flux after an overturn of an initially stratified lunar mantle might explain the existence and timing of an early lunar dynamo. Using a three-dimensional spherical convection model, we show that a dense layer, enriched in radioactive elements (a 'thermal blanket'), at the base of the lunar mantle can initially prevent core cooling, thereby inhibiting core convection and magnetic field generation. Subsequent radioactive heating progressively increases the buoyancy of the thermal blanket, ultimately causing it to rise back into the mantle. The removal of the thermal blanket, proposed to explain the eruption of thorium- and titanium-rich lunar mare basalts, plausibly results in a core heat flux sufficient to power a short-lived lunar dynamo.     

Mantle plumes： Why the current skepticism？

The present reappraisal of the mantle plume hypothesis is perhaps the most exciting current debate in Earth science. Nevertheless, the fundamental reasons for why it has arisen are often not well understood. They are that 1） many observations do not agree with the predictions of the original model, 2） it is possible that convection of the sort required to generate thermal plumes in the Earth＇s mantle does not occur, 3） so many variants of the original model have been invoked to accommodate conflicting data that the plume hypthesis is in practice no longer testable, and 4） alternative models are viable, though these have been largely neglected by researchers. Regardless of the final outcome, the present vigorous debate is to be welcomed since it is likely to stimulate new discoveries in a way that unquestioning acceptance of the conventional plume model will not.     

对流边界层顶部夹卷速度参数化的模拟研究

提出一个关于夹卷速度的新参数化方案，即混合层顶部夹卷速度W。可由对流边界层特征量W*、θ*、zi及覆盖逆温强度δθ^ /δz确定，对流槽实验模拟的测量结果证实了该参数化方案的合理性，并由实验数据拟合得出了参数化方案中的系数。数值模拟结果表明此方案可直接预报对流边界层高度zi，并能很好地预报出zi的变化过程，具有很强的实用性。