Thermal-image-derived dynamics of vertical ash plumes at Santiaguito volcano,Guatemala

Vertical ash plumes were imaged at Santiaguito (Guatemala) using a thermal camera to capture plume ascent dynamics. The plumes comprised a convecting plume front fed by a steady feeder plume. Of the 25 plumes imaged, 24 had a gas thrust region within which ascent velocities were 15–50 m s⁻¹. A transition to buoyant ascent occurred 20 to 50 m above the vent, where ascent velocities declined to 4–15 m s⁻¹. Plumes that attained greater heights had higher heat contents, wider feeder plumes and higher buoyant ascent velocities.     

Dynamics of co-ignimbrite plumes generated from pyroclastic flows of Mount St. Helens (7 August 1980)

 Four co-ignimbrite plumes were generated along the flow path of the pyroclastic flow of 7 August 1980 at Mount St. Helens. Three of the plumes were generated in discrete pulses which can be linked to changes in slope along the channel. One plume was generated at the mouth of the channel where the flow decelerated markedly as it moved onto the lower slopes of the pumice plain. Plume generation here may be triggered by enhanced mixing due to a hydraulic jump associated with an abrupt slope change. Measurements of plume ascent velocity and width show that the co-ignimbrite plumes increased in velocity with height. The plumes have initial velocities of 1–2 m/s. Two of the plumes reached a velocity maximum (4.6 and 8.8 m/s, respectively, at heights of 270 and 315 m above the flow) and thereafter decelerated. The other plumes reached velocities of 6.2 and 13 m/s. The four plumes become systematically less energetic downstream as measured by their ascent rates, which can be interpreted as a consequence of decreasing interaction of the pyroclastic flow front with the atmosphere. Theoretical models of both co-ignimbrite plumes and discrete co-ignimbrite clouds assume that there is no initial momentum, and both are able to predict the observed acceleration stage. The rising plumes mix with and heat air and sediment out particles causing their buoyancy to increase. Theoretical models agree well with observations and suggest that the initial motion of the ascending material is best described as a discrete thermal cloud which expands as it entrains air, whereas the subsequent motion of the head may become influenced by material supplied from the following plume. The models agree well with observations for an initial temperature of the ash and air mixture in the range of 500–600 K, which is in turn consistent with the measured initial ash temperature of around 920 K. Ash masses of 3.4×10⁵ to 1.8×10⁶ kg are estimated. Received: 11 January 1996 / Accepted: 7 October 1996     

Mantle plumes in the models of quasi-turbulent thermal convection

The process of multiple self-nucleation and ascent of mantle plumes is studied in the numerical models of thermal convection. The plumes are observed even in the simplest isoviscous models of thermal convection that leave aside the more complex rheology of the material, thermochemical effects, phase transformations, etc., which, although controlling the features of plumes, are not necessary for their formation. The origin of plumes is mainly due to the instability of the mantle flows at highly intense (low-viscous) thermal convection. At high viscosity, convective flows form regular cells. As viscosity decreases, the ascending and descending flows become narrower and unsteady. At a further decrease in viscosity, the ascending plumes assume a mushroom-like shape and occasionally change their position in the mantle. The lifetime of each flow can attain 100 Ma. Using markers allows visualizing the evolution of the shape of the mantle plumes.     

Volcanic eruption clouds and the thermal power output of explosive eruptions

The maximum height attained by a volcanic eruption cloud is principally determined by the convective buoyancy of the mixture of volcanic gas + entrained air + fine-sized pyroclasts within the cloud. The thermal energy supplied to convection processes within an eruption cloud is derived from the cooling of pyroclastic material and volcanic gases discharged by an explosive eruption. Observational data from six recent eruptions indicates that the maximum height attained by volcanic eruption clouds is positively correlated with the rate at which pyroclastic material is produced by an explosive eruption (correlation coefficient r = + 0.97). The ascent of industrial hot gas plumes is also governed by the thermal convection process. Empirical scaling relationships between plume height and thermal flux have been developed for industrial plumes. Applying these scaling relationships to volcanic eruption clouds suggests that the rate at which thermal energy is released into the atmosphere by an explosive eruption increases in an approximately linear manner as an eruption's pyroclastic production rate increases.     

On the rise of strongly tilted mantle plume tails

The rise of an initially horizontal, buoyant cylinder of fluid through a denser fluid at low Reynolds number is used to look at the ascent of strongly tilted mantle plumes through the mantle. Such ascents are characterized by (1) the growth of instabilities and (2) the development of a thermal wake downstream. Three-dimensional numerical experiments were carried out to examine these features. An hybrid particle-in-cell finite element method was used to look at the rise of non-diffusing cylinders and, a standard finite element method was used to look at the diffusing case. First the experiments show that the timescale of the fastest growing instability vary with the Rayleigh number and the viscosity ratio. In particular the growth rate decreases as the Rayleigh number decreases, in agreement with our analysis of the laboratory experiments of Kerr et al. (2008). Second the experiments show that the length of the thermal wake increases with the Rayleigh number but the change in viscosity has almost no influence on the wake length. Applied to strongly tilted mantle plumes we conclude that such plumes cannot be unstable given the plume timescales. We also discuss the application of this conclusion to weakly tilted plumes. Besides, this study allows to predict that mantle plumes are unlikely to have developed a significant thermal wake by the time they reach the surface. Finally, the resolution that is required to allow for the growth of mantle plume tails by combined diffusion and thermal entrainment is shown to represent a challenge for the large scale mantle convection simulations.     

Scoria cone formation through a violent Strombolian eruption: Irao Volcano, SW Japan

Scoria cones are common volcanic features and are thought to most commonly develop through the deposition of ballistics produced by gentle Strombolian eruptions and the outward sliding of talus. However, some historic scoria cones have been observed to form with phases of more energetic violent Strombolian eruptions (e.g., the 1943–1952 eruption of Parícutin, central Mexico; the 1975 eruption of Tolbachik, Kamchatka), maintaining volcanic plumes several kilometers in height, sometimes simultaneous with active effusive lava flows. Geologic evidence shows that violent Strombolian eruptions during cone formation may be more common than is generally perceived, and therefore it is important to obtain additional insights about such eruptions to better assess volcanic hazards. We studied Irao Volcano, the largest basaltic monogenetic volcano in the Abu Monogenetic Volcano Group, SW Japan. The geologic features of this volcano are consistent with a violent Strombolian eruption, including voluminous ash and fine lapilli beds (on order of 10^?1 km³ DRE) with simultaneous scoria cone formation and lava effusion from the base of the cone. The characteristics of the volcanic products suggest that the rate of magma ascent decreased gradually throughout the eruption and that less explosive Strombolian eruptions increased in frequency during the later stages of activity. During the eruption sequence, the chemical composition of the magma became more differentiated. A new K–Ar age determination for phlogopite crystallized within basalt dates the formation of Irao Volcano at 0.4?±?0.05 Ma.     

On the possibility of crater formation associated with an ascending plume

A hypothetical possibility of a qualitative explanation of large crater formation on the surfaces of the Moon and Mercury is discussed in terms of the concept of thermal mantle plumes. Prerequisites to this hypothesis are revealed under the assumption that the model equation of state of SiO₂ exhibiting an anomaly (a negative coefficient of thermal expansion) in the range of states approximately corresponding to average conditions typical of mantles of minor planets is applicable, in a first approximation, to mantle material. The anomaly reduces the buoyancy of hot plume material in such a way that, under conditions of moderate overheating, only relatively high columns comparable in size to the mantle are capable of ascending from the mantle bottom to the crust; allows cold peripheral material surrounding the hot column to be pushed away; causes compaction of the vertical zone of the contact of the column with the surrounding medium at the first stages after the plume ascent; and leads to compaction of the deep mantle due to the long-term heat supply. Such phenomena can lead to vertical craterlike deformations of the crust in areas of ascending large plumes whose presence can be supposed at early stages of the existence of minor planets. Significant implications of such an anomaly for geophysical processes can also be postulated.     

Determination of first-order degradation rate constants from monitoring networks

In this article, different strategies for estimating first-order degradation rate constants from measured field data are compared by application to multiple, synthetic, contaminant plumes. The plumes were generated by numerical simulation of contaminant transport and degradation in virtual heterogeneous aquifers. These sites were then individually and independently investigated on the computer by installation of extensive networks of observation wells. From the data measured at the wells, that is, contaminant concentrations, hydraulic conductivities, and heads, first-order degradation rates were estimated by three 1D centerline methods, which use only measurements located on the plume axis, and a two-dimensional method, which uses all concentration measurements available downgradient from the contaminant source. Results for both strategies show that the true rate constant used for the numerical simulation of the plumes in general tends to be overestimated. Overestimation is stronger for narrow plumes from small source zones, with an average overestimation factor of about 5 and single values ranging from 0.5 to 20, decreasing for wider plumes, with an average overestimation factor of about 2 and similar spread. Reasons for this overestimation are identified in the velocity calculation, the dispersivity parameterization, and off-centerline measurements. For narrow plumes, the one- and the two-dimensional strategies show approximately the same amount of overestimation. For wider plumes, however, incorporation of all measurements in the two-dimensional approach reduces the estimation error. No significant relation between the number of observation wells in the monitoring network and the quality of the estimated rate constant is found for the two-dimensional approach.     

Age of the Hawaiian-Emperor bend

⁴⁰Ar/³⁹Ar age data on alkalic and tholeiitic basalts from Diakakuji and Kinmei Seamounts in the vicinity of the Hawaiian-Emperor bend indicate that these volcanoes are about 41 and 39 m.y. old, respectively. Combined with previously published age data on Yuryaku and Ko¯ko Seamounts, the new data indicate that the best age for the bend is 42.0 ± 1.4 m.y.Petrochemical data indicate that the volcanic rocks recovered from bend seamounts are indistinguishable from Hawaiian volcanic rocks, strengthening the hypothesis that the Hawaiian-Emperor bend is part of the Hawaiian volcanic chain.⁴⁰Ar/³⁹Ar total fusion ages on altered whole-rock basalt samples are consistent with feldspar ages and with⁴⁰Ar/³⁹Ar incremental heating data and appear to reflect the crystallization ages of the samples even though conventional K-Ar ages are significantly younger. The cause of this effect is not known but it may be due to low-temperature loss of³⁹Ar from nonretentive montmorillonite clays that have also lost⁴⁰Ar.     

Magma dynamics, crystallization, and chemical differentiation of the 1959 Kilauea Iki lava lake, Hawaii, revisited

Using constraints from an extensive database of geological and geochemical observations along with results from fluid mechanical studies of convection in magma chambers, we identify the main physical processes at work during the solidification of the 1959 Kilauea Iki lava lakes. In turn, we investigate their quantitative influence on the crystallization and chemical differentiation of the magma, and on the development of the internal structure of the lava lake. In contrast to previous studies, vigorous stirring in the magma, driven predominately by the descent of dense crystal-laden thermal plumes from the roof solidification front and the ascent of buoyant compositional plumes due to the in situ growth of olivine crystals at the floor, is predicted to have been an inevitable consequence of very strong cooling at the roof and floor. The flow is expected to have caused extensive but imperfect mixing over most of the cooling history of the magma, producing minor compositional stratification at the roof and thermal stratification at the floor. The efficient stirring of the large roof cooling is expected to have resulted in significant internal nucleation of olivine crystals, which ultimately settled to the floor. Additional forcing due to either crystal sedimentation or the ascent of gas bubbles is not expected to have increased significantly the amount of mixing. In addition to convection in the magma, circulation driven by the convection of buoyant interstitial melt in highly permeable crystal-melt mushes forming the roof and the floor of the lava lake is envisaged to have produced a net upward flow of evolved magma from the floor during solidification. In the floor zone, mush convection may have caused the formation of axisymmetric chimneys through which evolved magma drained from deep within the floor into the overlying magma and potentially the roof. We hypothesize that the highly evolved, pipe-like ‘vertical olivine-rich bodies’ (VORBs) [Bull. Volcanol. 43 (1980) 675] observed in the floor zone, of the lake are fossil chimneys. In the roof zone, buoyant residual liquid both produced at the roof solidification front and gained from the floor as a result of incomplete convective mixing is envisaged to have percolated or ‘leaked‘ into the overlying highly-permeable cumulate, displacing less buoyant interstitial melt downward. The results from Rayleigh fractionation-type models formulated using boundary conditions based on a quantitative understanding of the convection in the magma indicate that most of the incompatible element variation over the height of the lake can be explained as a consequence of a combination of crystal settling and the extensive but imperfect convective mixing of buoyant residual liquid released from the floor solidification front. The remaining chemical variation is understood in terms of the additional influences of mush convection in the roof and floor on the vertical distribution of incompatible elements. Although cooling was concentrated at the roof of the lake, the floor zone is found to be thicker than the roof zone, implying that it grew more quickly. The large growth rate of the floor is explained as a consequence of a combination of the substantial sedimentation of olivine crystals and more rapid in situ crystallization due to both a higher liquidus temperature and enhanced cooling resulting from imperfect thermal and chemical mixing.     

Coupled conduit and atmospheric dispersal dynamics of the AD 79 Plinian eruption of Vesuvius

The AD 79 eruption of Vesuvius is certainly one of the most investigated explosive eruptions in the world. This makes it particularly suitable for the application of numerical models since we can be quite confident about input data, and the model predictions can be compared with field-based reconstruction of the eruption dynamics. Magma ascent along the volcanic conduit and the dispersal of pyroclasts in the atmosphere were simulated. The conduit and atmospheric domain were coupled through the flow conditions computed at the conduit exit. We simulated two different peak phases of the eruption which correspond to the emplacement of the white and gray magma types that produced Plinian fallout deposits with interlayered pyroclastic flow units during the gray phase. The input data, independently constrained and representative of each of the two eruptive phases, consist of liquid magma composition, crystal and water content, mass flow rate, and pressure–temperature–depth of the magma at the conduit entrance. A parametric study was performed on the less constrained variables such as microlite content of magma, pressure at the conduit entrance, and particle size representative of the eruptive mixture. Numerical results are substantially consistent with the reconstructed eruptive dynamics. In particular, the white eruption phase is found to lead to a fully buoyant eruption plume in all cases investigated, whereas the gray phase shows a more transitional character, i.e. the simultaneous production of a buoyant convective plume and pyroclastic surges, with a significant influence of the microlite content of magma in determining the partition of pyroclast mass between convective plumes and pyroclastic flows.     

Flow rates in the axial hot springs of the East Pacific Rise (21°N): Implications for the heat budget and the formation of massive sulfide deposits

The first in-situ measurements of flow rates in submarine hot springs (temperatures between 275° and 350°C) were made in the hot springs of the East Pacific Rise at 21°N during November, 1981. The flow rates ranged between 0.7 and 2.4 m/s. We estimate that the total rates at which heat and mass are transported from the three hot spring areas—OBS, National Geographic, and Southwest Vents—are ca. 2.2×10⁸ W and 150 kg of fluid/s respectively. The very high rate of heat loss precludes the existence of hot springs as steady state features of the East Pacific Rise at 21°N. We estimate that hot springs are active for a maximum of ca. 40,000 years. Chemical and/or mechanical clogging of the vents can drastically reduce the lifetime of the hot springs.Using a simple model of plume development [1], we calculated that less than 3% of the sulfide particles that are entrained in the hydrothermal plumes at 21°N settle from the plumes before dispersal by a lateral submarine current at a height of ca. 250 m above the seafloor [2]. The time that is required to form the sulfide mounds that surround the hot springs by the accumulation of sulfide particles that settle from the plumes is ca. 70–85 years. The current rate at which sulfide particles settle from the hydrothermal plumes and accumulate near the hot springs appears to be too small to lead to the formation of large massive sulfide deposits such as the Mavrouni deposit on Cyprus [3], within the time limits that are suggested by our thermal energy balance calculations.     

动态地幔柱尾管结构

动态地幔柱模式被广泛用于讨论地球科学中的一些重要课题,如巨大火成区的成因,冈瓦纳古陆解体的原因,板块内中小尺度动力过程的驱动因素等．但是这个基于实验研究而建立的模式中,忽略了地幔柱尾管特征及其作用．地幔柱尾管内温度和速度分布是研究地幔柱上升过程的必要条件．本文从控制尾管结构的基本方程出发,给出了一个定常轴对称地幔柱温度和速度分布的近似分析解．从而得到尾管结构的基本特征：影响尾管内温度分布的主导因素是地幔柱的热流通量,而尾管内上涌速度的大小则不仅取决于热流通量,主要是取决于地幔粘度随深度的变化方式．结果表明,对弱地幔柱,尾管的热损失可能是不可忽略的,而对强地幔柱,径向质量传递可能是不可忽略的.     

A simple method for calculating growth rates of petroleum hydrocarbon plumes

Consumption of aquifer Fe(III) during biodegradation of ground water contaminants may result in expansion of a contaminant plume, changing the outlook for monitored natural attenuation. Data from two research sites contaminated with petroleum hydrocarbons show that toluene and xylenes degrade under methanogenic conditions, but the benzene and ethylbenzene plumes grow as aquifer Fe(III) supplies are depleted. By considering a one-dimensional reaction front in a constant unidirectional flow field, it is possible to derive a simple expression for the growth rate of a benzene plume. The method balances the mass flux of benzene with the Fe(III) content of the aquifer, assuming that the biodegradation reaction is instantaneous. The resulting expression shows that the benzene front migration is retarded relative to the ground water velocity by a factor that depends on the concentrations of hydrocarbon and bioavailable Fe(III). The method provides good agreement with benzene plumes at a crude oil study site in Minnesota and a gasoline site in South Carolina. Compared to the South Carolina site, the Minnesota site has 25% higher benzene flux but eight times the Fe(III), leading to about one-sixth the expansion rate. Although it was developed for benzene, toluene, ethylbenzene, and xylenes, the growth-rate estimation method may have applications to contaminant plumes from other persistent contaminant sources.     

Pre-eruptive and syn-eruptive conditions in the Black Butte, California dacite: Insight into crystallization kinetics in a silicic magma system

A series of experiments and petrographic analyses have been run to determine the pre-eruption phase equilibria and ascent dynamics of dacitic lavas composing Black Butte, a dome complex on the flank of Mount Shasta, California. Major and trace element analyses indicate that the Black Butte magma shared a common parent with contemporaneously erupted magmas at the Shasta summit. The Black Butte lava phenocryst phase assemblage (20 v.%) consists of amphibole, plagioclase (core An_77.5), and Fe–Ti oxides in a fine-grained (< 0.5 mm) groundmass of plagioclase, pyroxene, Fe–Ti oxides, amphibole, and cristobalite. The phenocryst assemblage and crystal compositions are reproduced experimentally between 890 °C and 910 °C, ≥ 300 MPa, X_H2O = 1, and oxygen fugacity = NNO + 1. This study has quantified the extent of three crystallization processes occurring in the Black Butte dacite that can be used to discern ascent processes. Magma ascent rate was quantified using the widths of amphibole breakdown rims in natural samples, using an experimental calibration of rim development in a similar magma at relevant conditions. The majority of rims are 34 ± 10 μm thick, suggesting a time-integrated magma ascent rate of 0.004–0.006 m/s among all four dome lobes. This is comparable to values for effusive samples from the 1980 Mount St. Helens eruption and slightly faster than those estimated at Montserrat. A gap between the compositions of plagioclase phenocryst cores and microlites suggests that while phenocryst growth was continuous throughout ascent, microlite formation did not occur until significantly into ascent. The duration of crystallization is estimated using the magma reservoir depth and ascent rate, as determined from phase equilibria and amphibole rim widths, respectively. Textural analysis of the natural plagioclase crystals yields maximum growth rates of plagioclase phenocryst rims and groundmass microlites of 8.7 × 10^− 8 and 2.5 × 10^− 8 mm/s, respectively. These rates are comparable to values determined from time-sequenced samples of dacite erupted effusively from Mount St. Helens during 1980 and indicate that syneruptive crystallization processes were important during the Black Butte eruptive cycle.     

Complex dynamics of small-moderate volcanic events: the example of the 2011 rhyolitic Cordón Caulle eruption,Chile

After decades of repose, Puyehue-Cordón Caulle Volcano (Chile) erupted in June 2011 following a month of continuously increasing seismic activity. The eruption dispersed a large volume of rhyolitic tephra over a wide area and was characterized by complex dynamics. During the initial climactic phase of the eruption (24–30 h on 4–5 June), 11–14-km-high plumes dispersed most of the erupted tephra eastward towards Argentina, reaching as far as the Atlantic Ocean. This first eruptive phase was followed by activity of lower intensity, leading to the development of a complex stratigraphic sequence, mainly due to rapid shifts in wind direction and eruptive style. The resulting tephra deposits consist of 13 main layers grouped into four units. Each layer was characterized based on its dispersal direction, sedimentological features, and on the main characteristics of the juvenile fraction (texture, density, petrography, chemistry). The lowest part of the eruptive sequence (Unit I), corresponding to the tephra emitted between 4 and 5 June, is composed of alternating lapilli layers with a total estimated volume of ca. 0.75 km³; these layers record the highest intensity phase, during which a bent-over plume dispersed tephra towards the southeast-east, with negligible up-wind sedimentation. Products emitted during 5–6 June (Unit II) signaled an abrupt shift in wind direction towards the north, leading to the deposition of a coarse ash deposit in the northern sector (ca. 0.21 km³ in volume), followed by a resumption of easterly directed winds. A third phase (Unit III) began on 7 June and resulted in tephra deposits in the eastern sector and ballistic bombs around the vent area. A final phase (Unit IV) started after 15 June and was characterized by the emission of fine-grained white tephra from ash-charged plumes during low-level activity and the extrusion of a viscous lava flow. Timing and duration of the first eruptive phases were constrained based on comparison of the dispersal of the main tephra layers with satellite images, showing that most of the tephra was emitted during the first 72 h of the event. The analyzed juvenile material tightly clusters within the rhyolitic field, with negligible chemical variations through the eruptive sequence. Textural observations reveal that changes in eruption intensity (and consequently in magma ascent velocity within the conduit) and complex interactions between gas-rich and gas-depleted magma portions during ascent resulted in vesicular clasts with variable degrees of shear localization, and possibly in the large heterogeneity of the juvenile material.     

Review of Quantitative Surveys of the Length and Stability of MTBE,TBA, and Benzene Plumes in Groundwater at UST Sites

Quantitative information regarding the length and stability condition of groundwater plumes of benzene, methyl tert‐butyl ether (MTBE), and tert‐butyl alcohol (TBA) has been compiled from thousands of underground storage tank (UST) sites in the United States where gasoline fuel releases have occurred. This paper presents a review and summary of 13 published scientific surveys, of which 10 address benzene and/or MTBE plumes only, and 3 address benzene, MTBE, and TBA plumes. These data show the observed lengths of benzene and MTBE plumes to be relatively consistent among various regions and hydrogeologic settings, with median lengths at a delineation limit of 10 µg/L falling into relatively narrow ranges from 101 to 185 feet for benzene and 110 to 178 feet for MTBE. The observed statistical distributions of MTBE and benzene plumes show the two plume types to be of comparable lengths, with 90th percentile MTBE plume lengths moderately exceeding benzene plume lengths by 16% at a 10‐µg/L delineation limit (400 feet vs. 345 feet) and 25% at a 5‐µg/L delineation limit (530 feet vs. 425 feet). Stability analyses for benzene and MTBE plumes found 94 and 93% of these plumes, respectively, to be in a nonexpanding condition, and over 91% of individual monitoring wells to exhibit nonincreasing concentration trends. Three published studies addressing TBA found TBA plumes to be of comparable length to MTBE and benzene plumes, with 86% of wells in one study showing nonincreasing concentration trends.     

The transport of xenoliths in magmas

Field measurements of the rheology of Hawaiian and Etnean lavas have shown that, on eruption, they behave as Bingham liquids with yield strengths in the range 70–450 N/m². Ultramafic nodules entrained in a Bingham liquid cannot settle unless the stress they impose on the liquid exceeds approximately 5–7 times its yield strength. Consequently magmas with yield strengths of 10–1000 N/m² can transport xenoliths up to 30 cm diameter without settling occurring. The size of nodules commonly observed rarely exceeds 30 cm. A review of experimental data shows that, when conditions are appropriate for settling, the terminal velocities of nodules in magmas are substantially slower if a Bingham rather than Newtonian Model is assumed.The view that large nodules imply fast rates of magma ascent is rejected. A case is presented for slow rates of ascent being more suitable for nodule transport as there is more opportunity for cooling, crystallisation and hence development of a yield strength. The relative abundance of nodules in the alkaline suite may be a consequence of their slow rates of ascent, whereas their absence in tholeiitic melts may be a consequence of rapid ascent rates. This interpretation is compatible with deductions on their relative rates of ascent based on other geological evidence.     

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.     

Nature and origin of microphenocrysts in a basalt

In the upper 4 cm of a basaltic lava 2 kilometers from its source there is little or no significant change in the size, shape, number, and volume of plagioclase microphenocrysts. A maximum rate of clongation of 0.04 micrometers per second would be consistent with the data. These facts suggest that the microphenocrysts developed prior to extrusion, probably during ascent and efferveseence.