Mantle plumes link magnetic superchrons to phanerozoic mass depletion events

The four most recent large mass extinction events in the Phanerozoic – the Cretaceous–Tertiary (KT), the Triassic–Jurassic (TJ), and the Permo-Triassic (PT) and Guadalupian–Tatarian (GT) doublet – are associated with a major flood basalt eruption, with the timing of peak volcanic activity corresponding within measurement uncertainties to the extinction event. Three magnetic superchrons precede the four largest Phanerozoic extinctions. The Cretaceous Long Normal Superchron (duration  35 Myr) precedes the KT and the Permian Kiaman Long Reversed Superchron ( 50 Myr) precedes the PT–GT doublet. In addition, the newly recognized Ordovician Moyero Long Reversed Superchron ( 30 Myr) precedes the end-Ordovician extinction event. There is a 10–20 Myr delay between the end of each superchron and the subsequent mass depletion event, both of which represent distant outliers from their respective populations. We propose that deep mantle plumes link these seemingly unrelated phenomena. Long-term ( 200 Myr) variations in mantle convection possibly associated with the Wilson cycle induce temporal and spatial variations in heat flow at the core–mantle boundary. Polarity reversals are frequent when core heat flow is high and infrequent when it is low. Thermal instabilities in the D”-layer of the mantle increase core heat flow, end the magnetic superchron, and generate deep mantle plumes. The plumes ascend through the mantle on a 20 Myr time scale, producing continental flood basalt (trap) eruptions, rapid climatic change, and massive faunal depletions.     

The geomagnetic field at the Paleozoic/Mesozoic and Mesozoic/Cenozoic boundaries and lower mantle plumes

The data on the amplitude of variations in the direction and paleointensity of the geomagnetic field and the frequency of reversals throughout the last 50 Myr near the Paleozoic/Mesozoic and Mesozoic/Cenozoic boundaries, characterized by peaks of magmatic activity of Siberian and Deccan traps, and data on the amplitude of variations in the geomagnetic field direction relative to contemporary world magnetic anomalies are generalized. The boundaries of geological eras are not fixed in recorded paleointensity, polarity, reversal frequency, and variations in the geomagnetic field direction. Against the background of the “normal” field, nearly the same tendency of an increase in the amplitude of field direction variations is observed toward epicenters of contemporary lower mantle plumes; Greenland, Deccan, and Siberian superplumes; and world magnetic anomalies. This suggests a common origin of lower mantle plumes of various formation times, world magnetic anomalies, and the rise in the amplitude of geomagnetic field variations; i.e., all these phenomena are due to a local excitation in the upper part of the liquid core. Large plumes arise in intervals of the most significant changes in the paleointensity (drops or rises), while no correlation exists between the plume generation and the reversal frequency: times of plume formation correlate with the very diverse patterns of the frequency of reversals, from their total absence to maximum frequencies, implying that world magnetic anomalies, variations in the magnetic field direction and paleointensity, and plumes, on the one hand, and field reversals, on the other, have different sources. The time interval between magmatic activity of a plume at the Earth’s surface and its origination at the core-mantle boundary (the time of the plume rise toward the surface) amounts to 20–50 Myr in all cases considered. Different rise times are apparently associated with different paths of the plume rise, “delays” in the plume upward movement, and so on. The spread in “delay” times of each plume can be attributed to uncertainties in age determinations of paleomagnetic study objects and/or the natural remanent magnetization, but it is more probable that this is a result of the formation of a series of plumes (superplumes) in approximately the same region at the core-mantle boundary in the aforementioned time interval. Such an interpretation is supported by the existence of compact clusters of higher field direction amplitudes between 300 and 200 Ma that are possible regions of formation of world magnetic anomalies and plumes.     

On the advection of sharp material interfaces in geodynamic problems: entrainment of the D″ layer

The D″ layer is a dense and chemically distinct layer at the base of the convecting mantle. Numerical modeling of the entrainment of this layer by mantle convection requires the solution of the advective transport equation without introducing numerical diffusion across sharp material boundaries. We use our improved second moment numerical method to solve the equation. The method conserves the amount of material and the first and second moments of material distribution in each control volume. We first consider two examples of isothermal Rayleigh–Taylor instability to illustrate the performance of our method by comparing our results with those of a number of field, tracer and marker chain methods. We show that the performance of our method in minimizing the numerical diffusion is better than the field methods and comparable to the tracer and marker chain methods. We then study the instability of the dense D″ layer and its interaction with the overlying mantle. A range of density contrast between the D″ layer and the mantle, layer thickness, and the Rayleigh number, Ra, is examined. We show that for higher values of these parameters, the amount of entrainment decreases and the layer remains stable over longer periods of time. For very thick D″ layers and high Ra values, internal convection can take place within the layer.     

地球磁极倒转与白垩纪超静磁带的非线性分析

地球磁场多次发生南北(正负)磁极位置的变换和白垩纪超静磁带(CNS)的异常现象,这已为大家所公认.但造成这种异常现象的原因,则是迄今未能很好解答的一个难题. 应用非线性理论对地球磁极倒转和白垩纪超静磁带进行了分析, 认为超静磁带事件意味着地球核幔相互作用和外核流体运动可能处于能量最低的状态,地球磁场系统通过不断地与外界交换物质和能量,维持一种空间或时间的有序结构.在121～83Ma期间,无外星撞击地球引起地磁极性倒转,可能是白垩纪超静磁带出现的原因之一.地球磁场极性的随机倒转具有混沌运动的自逆转特性,混沌理论给地磁极性倒转提出了一个简明的动力机制解释.     

Paleointensity of the geomagnetic field in the Cretaceous (from Cretaceous rocks of Mongolia)

A representative collection of Cretaceous rocks of Mongolia is used for the study of the magnetic properties of the rocks and for determination of the paleodirections and paleointensities H _anc of the geomagnetic field. The characteristic NRM component in the samples is recognized in the temperature interval from 200 to 620–660°C. The values of H _anc are determined by the Thellier-Coe method with observance of all present-day requirements regarding the reliability of such kind of results. Comparison of data in the literature on paleointensity in the Cretaceous superchron and in the Miocene supports the hypothesis of the inverse correlation between the average intensity of the paleofield and the frequency of geomagnetic reversals. The increase in the average intensities is accompanied by an appreciable increase in the variance of the virtual dipole moment (VDM). We suggest that the visible increase in the average VDM value in the superchron is due to the greater variability of VDM in this period compared to the Miocene.     

Three distinct reversing modes in the geodynamo

The data that describe the long-term reversing behavior of the geodynamo show strong and sudden changes in magnetic reversal frequency. This concerns both the onset and the end of superchrons and most probably the occurrence of episodes characterized by extreme geomagnetic reversal frequency (>10–15 rev./Myr). To account for the complexity observed in geomagnetic reversal frequency evolution, we propose a simple scenario in which the geodynamo operates in three distinct reversing modes: i—a “normal” reversing mode generating geomagnetic polarity reversals according to a stationary random process, with on average a reversal rate of ～3 rev./Myr; ii—a non-reversing “superchron” mode characterizing long time intervals without reversal; iii—a hyper-active reversing mode characterized by an extreme geomagnetic reversal frequency. The transitions between the different reversing modes would be sudden, i.e., on the Myr time scale. Following previous studies, we suggest that in the past, the occurrence of these transitions has been modulated by thermal conditions at the core-mantle boundary governed by mantle dynamics. It might also be possible that they were more frequent during the Precambrian, before the nucleation of the inner core, because of a stronger influence on geodynamo activity of the thermal conditions at the core-mantle boundary.     

Isotope systematics of noble gases in the Earth's mantle: possible sources of primordial isotopes and implications for mantle structure

The Earth's mantle contains a mixture of primordial noble gases, in particular solar-type helium and neon, and radiogenic rare gases from long-lived U, ²³²Th, ⁴⁰K and short-lived ¹²⁹I, ²⁴⁴Pu. Rocks derived from deep mantle plume magmatism like on Hawaii or Iceland contain a higher proportion of primordial nuclides than rocks from the shallow upper mantle, e.g. mid ocean ridge basalts (MORBs). This is widely regarded as the key evidence for survival of a less degassed and more “primitive” reservoir within the lower mantle. We present an evaluation of noble gas composition showing the shallow mantle to have about five times more radiogenic (relative to primordial) isotopes than Hawaii/Iceland-type plume reservoirs, no matter if short- or long-lived decay systems are considered. This fundamental property suggests that both MORB and plume-type noble gases are mixtures of: (1) a homogeneous radiogenic component present throughout most of the mantle and (2) a uniform primordial noble gas component with very minor radiogenic ingrowth. This conclusion depends crucially on the observed excess of radiogenic Xe in plume-derived rocks, and is only valid if this Xe excess is inherent to the plume sources.Possible sources of the primordial component of mantle plume reservoirs—and possibly also the MORB mantle—could be mantle reservoirs that remained relatively isolated over most of Earth's history (“blobs”, a deep abyssal layer, or the D” layer), but these need a considerable concentration of primordial gases to compensate U, Th, K decay over 4.5 Ga. Earth's core is evaluated as an alternative viable source feeding primordial nuclides into mantle reservoirs: even low metal-silicate partitioning coefficients allow sufficient primordial noble gases to be incorporated into the early forming core, as the undifferentiated proto-Earth was initially gas-rich. Massive mantle degassing soon after core formation then provides the opposite concentration gradient that allows primordial noble gases reentering the mantle at the core-mantle boundary, probably via partial mantle melts. Another possible source of primordial noble gases in Earth's mantle are subducted sediments containing extraterrestrial dust with solar He and Ne, but this supply mechanism crucially depends on largely unconstrained parameters. The latter two scenarios do not require the preservation of a “primitive” mantle reservoir over 4.5 Ga, and can potentially better reconcile increasing geochemical evidence of recycled lithospheric components in mantle plumes and seismic evidence for whole mantle convection.     

Structure of hydrothermal plumes at the Logatchev vent field, 14°45′N, Mid-Atlantic Ridge: evidence from geochemical and geophysical data

In the Seventh cruise of R/V “Professor Logatchev” anomalies of natural electric field (EF), Eh and pS were discovered using a towed instrument package (RIFT) at 14°45′N on the MAR (Logatchev hydrothermal field). The anomalous zone (AZ) is situated close (10–35 m) to two low-temperature venting areas of degrading sulphides and a black smoker (Irina-Microsmoke) forming a distinct buoyant plume. Over or close to the main area of high-temperature venting situated to the south-east from the AZ, no EF or Eh anomalies were observed. According to the results of Mir dives the highly mineralised solutions from smoking craters at the main mound mostly form non-buoyant plumes (reverse-plumes). The buoyant plume structure shows the differentiation of the electrical and Eh fields within the plume. Maxima of the EF, Eh and E_H2S anomalies were revealed in the lower part (15 m) of the plume. The negative redox potential plume coupled with a sulphide anomaly is more localized in comparison with the EF. This observation indicates a distinct change in the composition of buoyant plume water, which may be due to the formation and fallout of early formed Fe sulphide particles soon after venting.     

Estimates of heat and pyroclast discharge by volcanic eruptions based upon the eruption cloud and steady plume observations

Fumarolic steam plumes and eruption clouds rise like convetive turbulent columns into the atmosphere. Formulae are presented here for estimating the heat power of plumes, the production rate of juvenile pyroclasts ejected during eruptions and the heat output of fumaroles. Their accuracy is tested using the well-studied examples of eruptions of Kamchatkan volcanoes.The Briggs (1969) formula may be used in observing the ascending part of a plume in crosswinds. The best results have been obtained using the CONCAWE formula which permits estimation of the heat power in crosswinds based on the axis height of a horizontal part of a maintained plume. Three connected equations have been suggested for a stable atmosphere and calm weather conditions. The first one, which is applicable for heights ranging from 100 m to 1 km, is the formula proposed by Morton et al. (1956). This equation changes for higher layers of the troposphere (1–10 km) and stratosphere (10–55 km).A classification scale was constructed allowing us to compare volcanic eruptions and fumarolic activity in terms of the intensity of their plumes.The described method is useful for volcano surveillance; it helps in the study of the energetics and mechanics of volcanic and magmatic processes.     

Heat and mass transfer in a thermochemical plume under an oceanic plate far from the mid-ocean ridge axis

Heat and mass transfer processes in the conduit of a thermochemical plume located beneath an oceanic plate far from a mid-ocean ridge (MOR) proceed under conditions of horizontal convective flows penetrating the plume conduit. In the region of a mantle flow approaching the plume conduit (in the frontal part of the conduit), the mantle material heats and melts. The melt moves through the plume conduit at the average velocity of flow v and is crystallized on the opposite side of the conduit. The heat and the chemical dope transferred by the conduit to the mantle flow are carried away by crystallized mantle material at the velocity v. The main equations of heat and mass transfer are obtained for a thermochemical plume interacting with a horizontal convective mantle flow. The joint multiparameter problem of heat and mass transfer is solved for a thermochemical plume located far from an MOR axis. The dope concentration at the base of the plume is found as a function of the Lewis number. The Lewis numbers and, accordingly, the diffusion coefficients of the chemical dope in the plume conduit far from the MOR axis are determined.     

The dynamical and thermal structure of deep mantle plumes

If the interpretation of the D″ layer at the base of the mantle as a thermal boundary layer, with a temperature increment in the order of 800 K, is correct, then the formation of deep-mantle plumes to vent material from it appears inevitable. We demonstrate quantitatively that the strong temperature dependence of viscosity guides the upward flow into long-lived chimneys that are ～ 20 km in diameter near the base of the mantle and decrease in width with progressive upward softening and partial melting of plume material. The speed of flow up the axis of the plume is correspondingly fast; 1.6 m y^?1 at the base and 4.8 m y^?1 at 670 km depth. Thermal diffusive spreading of a heated plume is compensated by a slow horizontal convergence of mantle material toward the chimney in response to the lower pressure there. This convergence, which contributes only a small increment to the flux of material up the plume, also serves to throttle the flow in the chimney. The global plume mass flux necessary to transport 1.6 × 10¹² W of core heat upward through the mantle is 1.8 × 10⁶ kg s^?1. At its base, plume material is probably still significantly below its solidus or eutectic temperature, but substantial partial melting is very likely as it rises. We speculate that a small fraction of this fluid component eventually emerges at the surface in “hot spots”, with the fate of the remainder being unknown. The behaviour and properties of D″ and of plumes are closely coupled. Not only are plumes a necessary consequence of a thermal boundary layer, but their existence is impossible without that layer.     

Crustal heat transfer in the Taupo Volcanic Zone (New Zealand): comparison with other volcanic arcs and explanatory heat source models

The Taupo Volcanic Zone (TVZ) is a 200-km-long volcanic arc segment which developed ≤2 Ma ago within the continental crust of the North Island of New Zealand and lies at the southern end of the much larger Tonga-Kermadec arc system. The total crustal heat transfer of the TVZ is at present c. 2600 MW/100 km, most of the heat being transferred by convective geothermal systems. The rate of transfer is anomalously high in comparison to that of other active arcs, and arguably the highest world wide for such a setting. Heat transfer of other active arcs appear to vary almost linearly with subduction speed (about 150 MW/100 km for 10 mm/yr). The mass rate of common type arc extrusions (basalts, andesites, dacites) also increases almost linearly with subduction speed. This allows separation of the TVZ heat transfer into a “normal” component, associated with extrusions and intrusions of andesites and dacites (about 600 MW/100 km), and an “anomalous” component of about 2000 MW/100 km, related to extrusions and intrusions of rhyolitic melts whose generation is not directly controlled by subduction processes.Rhyolitic melts in the TVZ are partial melts of dominantly crustal origin. Comparison with other arcs indicates that the long-term extrusion rate of TVZ rhyolites (about 400 kg/s per 100 km) is also the highest world wide for this setting. The occurrence of voluminous Quaternary rhyolitic pyroclastics is a rare phenomenon and appears to be associated with a few arc segments (TVZ, Sumatra, Kyushu) that undergo significant crustal deformation.Various models have been proposed to explain the phenomenon of the anomalously high heat transfer within the TVZ. Models which require only heat transfer from plumes and subcrustal melts, either ponded at the crust/mantle boundary or intruding a spreading crust, are not suitable because the associated heat transfer at the contact is too low by a factor 2 to explain the required transfer rate of about 0.8 W/m² representing the “anomalous” crustal heat component of the TVZ. Heat generation by focussed plastic deformation within the ductile lithosphere is an alternative mechanism to explain “endogenous crustal heating” which yields heating rates that are also too low by a factor of two, although important parameters (average yield strength of lithosphere and opening rate of the TVZ) are not well known. A further search for a suitable combination of heat source models is required.     

Electrical signatures due to thermal anomalies along mobile belts reactivated by the trail and outburst of mantle plume: Evidences from the Indian subcontinent

In this study the geodynamical scenario along with concepts of mantle plume and mobile belts is utilized to show that most of the existing and potential high thermal regions fall along the (mobile arms affected by the outburst and) traces of mantle plumes. Effects of channeling and partitioning of thermomagmatic flux (TMF) due to these mantle plumes along the mobile belts, particularly near the triple junctions, can be seen in the form of high heat flow and presence of hot springs. Triple junctions manifest over the Indian lithosphere are: Kutch-Cambay, Narmada Son-Godavari, Tapi-Mahanadi, Tapi-Damodar, Pondicherry region, Gulf of Mannar and SW corner of the subcontinent (off-shore), etc. Apart from mobile belts, the deltaic regions of Krishna, Godavari, Ganga, Cauvery, Narmada-Tapi and Indus, etc., are also posses higher level of thermal anomalies as these regions seem to have been substantially influenced by outbursts and traces of Reunion, Kerguelen, Marion and Crozet hotspots. This is reflected from the correlation between plume affected mobile belts and high heat flow regions, large number of hot springs, anomalous electrical conductivity and also deformation or seismicity. Such correlation can be seen along Cambay, west coast trend, Narmada-Son lineament zone, Godavari-Damodar grabens and Bengal basin. Himalayan belt being ongoing collision zone, thermal anomalies are identified in the form of hot springs along the Himalayan arc. At some locations, which might be junction of tectonic trends, there exist significantly large thermal outputs. Puga in Himalayan region is one such example, as seen from high heat flow (max. 468 mW m^− 2) and geothermal gradient (234 °C/km max.). Similarly, Tatapani in Narmada Son Lineament (NSL) region is another such example. The present study discusses the correlation between thermal reservoirs identified by magnetotelluric (MT) study results and plume activity and suggests the need for systematic and detailed MT investigations along plume activated mobile strips in other regions to search for geodynamical history and geothermal resources.     

Deconvolution of overcritical P-wave reflections from the top of D″

We have developed a method to retrieve overcritical PdP reflections from the top of D″ from the data. The differential travel-time, amplitude ratio and phase advance of the reflected PdP waveform relative to the direct P waveform are accurately obtained. Kirchhoff synthetics have been used to model the reflected waveforms for a set of different incidence angles i and velocity contrasts δα at the reflector which were then used as reference wavelets for the deconvolution.The results from a Gräfenberg data set differ markedly from the values predicted by the one-dimensional PWDK model derived from the same data. However, there are large uncertainties in i and δα. This could be caused by the presence of significant lateral heterogeneity in D″ or the need to model other phases which may exist such as the diving waves PDP and PDDP. The method is also applied to two other data examples. The trade-offs between different model parameters are examined.     

Plume capture by divergent plate motions: implications for the distribution of hotspots,geochemistry of mid-ocean ridge basalts,and estimates of the heat flux at the core–mantle boundary

The coexistence of stationary mantle plumes with plate-scale flow is problematic in geodynamics. We present results from laboratory experiments aimed at understanding the effects of an imposed large-scale circulation on thermal convection at high Rayleigh number (10⁶≤Ra≤10⁹) in a fluid with a temperature-dependent viscosity. In a large tank, a layer of corn syrup is heated from below while being stirred by large-scale flow due to the opposing motions of a pair of conveyor belts immersed in the syrup at the top of the tank. Three regimes are observed, depending on the ratio V of the imposed horizontal flow velocity to the rise velocity of plumes ascending from the hot boundary, and on the ratio λ of the viscosity of the interior fluid to the viscosity of the hottest fluid in contact with the bottom boundary. When V≪1 and λ≥1, large-scale circulation has a negligible effect on convection and the heat flux is due to the formation and rise of randomly spaced plumes. When V>10 and λ>100, plume formation is suppressed entirely, and the heat flux is carried by a sheet-like upwelling located in the center of the tank. At intermediate V, and depending on λ, established plume conduits are advected along the bottom boundary and ascending plumes are focused towards the central upwelling. Heat transfer across the layer occurs through a combination of ascending plumes and large-scale flow. Scaling analyses show that the bottom boundary layer thickness and, in turn, the basal heat flux q depend on the Peclet number, Pe, and λ. When λ>10, q∝Pe^1/2 and when λ→1, q∝(Peλ)^1/3, consistent with classical scalings. When applied to the Earth, our results suggest that plate-driven mantle flow focuses ascending plumes towards upwellings in the central Pacific and Africa as well as into mid-ocean ridges. Furthermore, plumes may be captured by strong upwelling flow beneath fast-spreading ridges. This behavior may explain why hotspots are more abundant near slow-spreading ridges than fast-spreading ridges and may also explain some observed variations of mid-ocean ridge basalt (MORB) geochemistry with spreading rate. Moreover, our results suggest that a potentially significant fraction of the core heat flux is due to plumes that are drawn into upwelling flows beneath ridges and not observed as hotspots.     

Plume Mode of Thermal Convection in the Earth’s Mantle

In our previous works, based on numerical models, it was shown that under certain conditions a hot material can rise in portions in the tails of thermal mantle plumes. The spectrum of these pulsations can correspond to the observed spectra of catastrophic hotspot eruptions. Since most of the existing numerical models of thermal convection for the mantle of the present Earth do not reveal these pulsations, in this work, we analyze the physical cause and initiation conditions of pulsations of thermal plumes. The results of a numerical solution of the thermal convection equations for a material with varying parameters in the extended Boussinesq approximation are presented. It is shown how the structure of the convection is transformed with the increase of convection intensity. At the Rayleigh numbers Ra > 10⁶, convection becomes unsteady, and the configuration of the ascending and descending flows changes. The new flow emerging at the mantle bottom acquires a mushroom shape with a head and a tail. After the rise of the plume’s head to the surface, the tail remains in the mantle in the form of a quasi-stationary hot steam. It turns out that at Ra ~ 5 × 10⁷, the thermal mantle plume becomes pulsating and its tail is in fact a heated channel through which the hot material rises in successive portions. At the Rayleigh numbers Ra > 5 × 10⁸, the tail of the thermal plume breaks and the plume becomes a regular conveyor of separate ascending portions of the hot material, which are referred to as thermals. Thus, thermal convection with pulsating plumes takes place at the transitional stage from the regime of quasi-stationary plumes to the regime of thermals.     

Transitions in the style of mantle convection at high Rayleigh numbers

The pattern and style of mantle convection govern the thermal evolution, internal dynamics, and large-scale surface deformation of the terrestrial planets. In order to characterize the nature of heat transport and convective behaviour at Rayleigh numbers, Ra, appropriate for planetary mantles (between 10⁴ and 10⁸), we perform a set of laboratory experiments. Convection is driven by a temperature gradient imposed between two rigid surfaces, and there is no internal heating. As the Rayleigh number is increased, two transitions in convective behaviour occur. First we observe a change from steady to time-dependent convection at Ra≈10⁵. A second transition occurs at higher Rayleigh numbers, Ra≈5×10⁶, with large-scale time-dependent flow being replaced by isolated rising and sinking plumes. Corresponding to the latter transition, the exponent β in the power law relating the Nusselt number Nu to the Rayleigh number (NuRa^β) is reduced. Both rising and sinking plumes always consist of plume heads followed by tails. There is no characteristic frequency for the formation of plumes.     

Long-range dependence in the Cenozoic reversal record

The frequency distribution of intervals between Cenozoic geomagnetic reversals approximates a power law, while their occurrence over time shows temporal clustering of short and long intervals. Application of the Aggregate Variance and Absolute Value methods suggest long-range dependence in this time series, a possible indication that the geodynamo operates as a self-organised complex system. This hypothesis may allow the Cretaceous superchron to be considered as an integral part of the ordinary reversal regime attested in the Cenozoic record.     

Cretaceous plutonism in Central Tibet: an example of post-collision magmatism?

The magmatic province of the northern Lhasa Terrane includes an Early Cretaceous (120–130 Ma) plutonic event, and a Late Cretaceous (80–110 Ma) volcanic event. The plutonic association constitutes an older suite of granodiorites, monzogranites and tonalites and a younger peraluminous leucogranite facies. Plutonism occurred about 20 Ma after obduction of the Banggong ophiolite, following closure between the Lhasa and Qiantang Terranes.The earlier suite is of broadly calc-alkaline in composition but differs from arc-related magmas in that only more evolved compositions are represented (SiO₂ > 58%) and Rb/Zr ratios are elevated relative to the Gangdese batholith to the south. Trace-element and isotopic constraints are consistent with derivation from a Late Proterozoic amphibole-bearing crustal source requiring temperatures > 950°C during anatexis. The leucogranites require a pelitic source which is tentatively identified as the Nyaingentanglha basement exposed south of the plutonic province. Unlike the High Himalaya leucogranites, trace elements and field relations require a high degree of melting at source (> 50%) suggesting fluid-absent melting at temperatures > 850°C. Such high crustal temperatures indicate convective heat transfer from the mantle.Thermal constraints together with a tectonic setting of post-emplacement uplift followed by a marine transgression in the northern Lhasa Terrane can not be reconciled with a model of tectonically thickened crust but are consistent with post-collision attenuation of the lithosphere.     

Heat transfer between a thermochemical plume channel and the surrounding mantle in the presence of horizontal mantle flow

Heat and mass transfer processes in the conduit of a thermochemical plume located beneath an oceanic plate far from a mid-ocean ridge (MOR) proceed under conditions of horizontal convective flows penetrating the plume conduit. In the region of a mantle flow approaching the plume conduit (in the frontal part of the conduit), the mantle material heats and melts. The melt moves through the plume conduit at the average velocity of flow v and is crystallized on the opposite side of the conduit (in the frontal part of the conduit). The heat and the chemical dope transferred by the conduit to the mantle flow are carried away by crystallized mantle material at the velocity v.The local coefficients of heat transfer at the boundary of the plume conduit are theoretically determined and the balance of heat fluxes through the side of the plume conduit per linear meter of the conduit height. The total heat generation rate, transmitted by the Hawaiian plume into the upper and lower mantle, is evaluated. With the use of regular patterns of heat transfer in the lower mantle, which is modeled on the horizontal layer, heated from below and cooled from above, the diameter of the plume source, the kinematic viscosity of the melt in the plume conduit, and the velocity of horizontal lower-mantle flows are evaluated and the dependences of the temperature drop, viscosity and Rayleigh number for the lower mantle on the diameter of the plume source are presented.