Cumulate nodules as evidence for convective fractionation in a phonolite magma chamber

The physical mechanism by which chemical zonation develops in magma chambers has been controversial partly because unambiguous geological constraints have been lacking. The 11,000 years B.P. eruption of Laacher See Volcano produced a zoned tephra deposit and also ejected crystal-rich nodules which provide a snapshot of the materials crystallising at the magma chamber margins. New data on petrography and chemical compositions of nodules, their cumulate minerals and interstitial glasses are used to deduce the chemical evolution of the phonolite melt due to fractional crystallisation of the mineral assemblages. These data, together with those on the vertical zonation of the melt in the bulk of the chamber, are shown to be consistent with a model of stratification of the chamber by convective fractionation, in which a thin boundary layer of residual melt from fractional crystallisation ascends at the chamber side and accumulates at the roof. Crystallisation could have provided buoyancy to drive convection by enriching incompatible volatile components (mainly water) in the residual melt. Available fluid dynamic studies of single- and double-diffusive boundary layers are used to assess convection in the Laacher See chamber. The boundary layer is likely to have been: (1) laminar, which implies that the density gradient in the chamber steepened upwards; (2) in the counterflow regime, in which compositional and thermal layers flow in opposite directions; and (3) thin ( 10 cm), estimated from theory for a flat wall, suggesting that wall morphology could be important in determining boundary layer characteristics. Estimates of mass transfer rates due to this mechanism suggest that the chamber could have become stratified in a time of the order of 10³ years.     

Incipient magma chamber formation as a result of repetitive intrusions

An analytical solution for periodic magma intrusions in conduits was developed to study the onset of shallow magma chamber formation. The solution is based on determining when a repetitive series of intrusions can cause the wall rock of a conduit to reach its melt temperature. The results show that magma chamber formation in conduits is a strong function of the volume rate of intrusion and that magma chamber formation is likely when the intrusion rate exceeds 10^?3 km³/ yr. which agrees with observations by other investigators. Once this critical value of intrusion rate is reached, magma chambers are likely to begin forming after only a few intrusive pulses (less than ten). Results for both cylindrical conduits and dikes show cylindrical conduits are more favourable for the formation of shallow magma chambers.     

Volumetric characteristics of lava flows from interferometric radar and multispectral satellite data: the 1995 Fernandina and 1998 Cerro Azul eruptions in the western Galápagos

We have used a suite of remotely sensed data, numerical lava flow modeling, and field observations to determine quantitative characteristics of the 1995 Fernandina and 1998 Cerro Azul eruptions in the western Galápagos Islands. Flank lava flow areas, volumes, instantaneous effusion rates, and average effusion rates were all determined for these two eruptions, for which only limited syn-eruptive field observations are available. Using data from SPOT, TOPSAR, ERS-1, and ERS-2, we determined that the 1995 Fernandina flow covers a subaerial area of 6.5×10⁶ m² and has a subaerial dense rock equivalent (DRE) volume of 42×10⁶ m³. Field observations, ATSR satellite data, and the FLOWGO numerical model allow us to determine that the effusion rate declined exponentially from a high of ~60–200 m³ s^-1 during the first few hours to <5 m³ s^-1 prior to ceasing after 73 days, with a mean effusion rate of 4–16 m³ s^-1. Integrating the ATSR-derived, exponentially declining effusion rate over the eruption duration produces a total (subaerial + submarine) DRE volume of between 27 and 100×10⁶ m³, the range in values being due to differing assumptions about heat loss characteristics; only values in the higher part of this range are consistent with the independently derived subaerial volume. Using SPOT, TOPSAR, ERS-1, and ERS-2 data, we determine that the 1998 Cerro Azul flow is 16 km long, covers 16 km², and has a DRE volume of 54×10⁶ m³. FLOWGO produces at-vent velocity and effusion rate values of 11 m s^-1 and ~600 m³ s^-1, respectively. The velocity value agrees well with the 12 m s^-1 estimated in the field. The mean effusion rate (total DRE volume/duration) was 7–47 m³ s^-1. Dike dimensions, fissure lengths, and pressure gradients along the conduit based on magma chamber depth estimates of 3–5 km produce mean effusion rates for the two eruptions that range over nearly four orders of magnitude, the range being due to uncertainty in the magma viscosity, dike dimensions, and pressure gradient between magma chamber and vent. Although somewhat consistent with mean effusion rates from other techniques, their wide range makes them less useful. The exponentially declining effusion rates during both eruptions are consistent with release of elastic strain being the driving mechanism of the eruptions. Our results provide independent input parameters for previously published theoretical relationships between magma chamber pressurization and eruption rates that constrain chamber volumes and increases in volume prior to eruption, as well as time constants of exponential decay during the eruption. The results and theoretical relationships combine to indicate that at both volcanoes probably 25–30% of the volumetric increase in the magma chamber erupted as lava onto the surface. In both eruptions the lava flow volumes are less than 1% of the magma chamber volume.     

Voluminous lava-like precursor to a major ash-flow tuff: low-column pyroclastic eruption of the Pagosa Peak Dacite, San Juan volcanic field, Colorado

The Pagosa Peak Dacite is an unusual pyroclastic deposit that immediately predated eruption of the enormous Fish Canyon Tuff (5000 km³) from the La Garita caldera at 28 Ma. The Pagosa Peak Dacite is thick (to 1 km), voluminous (>200 km³), and has a high aspect ratio (1:50) similar to those of silicic lava flows. It contains a high proportion (40–60%) of juvenile clasts (to 3–4 m) emplaced as viscous magma that was less vesiculated than typical pumice. Accidental lithic fragments are absent above the basal 5–10% of the unit. Thick densely welded proximal deposits flowed rheomorphically due to gravitational spreading, despite the very high viscosity of the crystal-rich magma, resulting in a macroscopic appearance similar to flow-layered silicic lava. Although it is a separate depositional unit, the Pagosa Peak Dacite is indistinguishable from the overlying Fish Canyon Tuff in bulk-rock chemistry, phenocryst compositions, and ⁴⁰Ar/³⁹Ar age.The unusual characteristics of this deposit are interpreted as consequences of eruption by low-column pyroclastic fountaining and lateral transport as dense, poorly inflated pyroclastic flows. The inferred eruptive style may be in part related to synchronous disruption of the southern margin of the Fish Canyon magma chamber by block faulting. The Pagosa Peak eruptive sources are apparently buried in the southern La Garita caldera, where northerly extensions of observed syneruptive faults served as fissure vents. Cumulative vent cross-sections were large, leading to relatively low emission velocities for a given discharge rate. Many successive pyroclastic flows accumulated sufficiently rapidly to weld densely as a cooling unit up to 1000 m thick and to retain heat adequately to permit rheomorphic flow. Explosive potential of the magma may have been reduced by degassing during ascent through fissure conduits, leading to fracture-dominated magma fragmentation at low vesicularity. Subsequent collapse of the 75×35 km² La Garita caldera and eruption of the Fish Canyon Tuff were probably triggered by destabilization of the chamber roof as magma was withdrawn during the Pagosa Peak eruption.     

Magma contamination by direct wall rock interaction: constraints from xenoliths from the walls of a carbonate-hosted magma chamber (Vesuvius 1944 eruption)

During the 1944 eruption of Vesuvius different types of xenoliths were ejected. They represent fragments of the walls of a low volume (<0.5 km³) shallow (3–4 km depth) magma chamber. The study of these xenoliths enables us to estimate the amount of contamination occurring at the boundary of a high-T alkaline magma chamber hosted in carbonate rocks. The process of contamination of the magma by carbonates can be modelled, using isotopic and chemical data, as a mixing between magma and marbles. Mass exchanges occur at the boundary between the crystallizing magma and marble wall rocks, where endoskarn forms. The contamination of the solidification front of the chamber is very limited. The solidification front and the skarn shell effectively isolate the interior of the magma chamber from new inputs of contaminants from the carbonate wall rocks. Therefore, the main volume of magma, hosted in the magma chamber, did not undergo any significant mass exchange with the wall rocks.     

A model for the origin of large silicic magma chambers: precursors of caldera-forming eruptions

The relatively low rates of magma production in island arcs and continental extensional settings require that the volume of silicic magma involved in large catastrophic caldera-forming (CCF) eruptions must accumulate over periods of 10 ⁵ to 10 ⁶ years. We address the question of why buoyant and otherwise eruptible high-silica magma should accumulate for long times in shallow chambers rather than erupt more continuously as magma is supplied from greater depths. Our hypothesis is that the viscoelastic behavior of magma chamber wall rocks may prevent an accumulation of overpressure sufficient to generate rhyolite dikes that can propagate to the surface and cause an eruption. The critical overpressure required for eruption is based on the model of Rubin (1995a). An approximate analytical model is used to evaluate the controls on magma overpressure for a continuously or episodically replenished spherical magma chamber contained in wall rocks with a Maxwell viscoelastic rheology. The governing parameters are the long-term magma supply, the magma chamber volume, and the effective viscosity of the wall rocks. The long-term magma supply, a parameter that is not typically incorporated into dike formation models, can be constrained from observations and melt generation models. For effective wall-rock viscosities in the range 10 ¹⁸ to 10 ²⁰ Pa s ^–1, dynamical regimes are identified that lead to the suppression of dikes capable of propagating to the surface. Frequent small eruptions that relieve magma chamber overpressure are favored when the chamber volume is small relative to the magma supply and when the wall rocks are cool. Magma storage, leading to conditions suitable for a CCF eruption, is favored for larger magma chambers (>10 ² km ³) with warm wall rocks that have a low effective viscosity. Magma storage is further enhanced by regional tectonic extension, high magma crystal contents, and if the effective wall-rock viscosity is lowered by microfracturing, fluid infiltration, or metamorphic reactions. The long-term magma supply rate and chamber volume are important controls on eruption frequency for all magma chamber sizes. The model can explain certain aspects of the frequency, volume, and spatial distribution of small-volume silicic eruptions in caldera systems, and helps account for the large size of granitic plutons, their association with extensional settings and high thermal gradients, and the fact that they usually post-date associated volcanic deposits.     

Degassing of stratified magma by compositional convection

A new model is proposed for passive degassing from sub-volcanic magma chambers. The water content in stably stratified shallow magma chamber will be equated to its solubility at the upper boundary by convection. Water from a lower layer high in water content can enrich the contact zone of the upper layer and lead to further convective overturn of this boundary layer. A complete set of equations describing convection with bubble formation and dissolution is reduced to a simplified form by assuming a small bubble content. The development and pattern of flow driven by vesiculation is modeled numerically in a 2D magma chamber for relatively low Raleigh numbers (5×10⁵). Bubbles rising from the magma will collect near the roof in a layer of 8–10 vol% and then escape upward to fumaroles. The Stokes flux of bubbles escaping from an andesitic magma with viscosity 10⁴ P and a top surface of about 500×500 m corresponds with observed total magmatic water fluxes of 35 kg/s. Pressure within the chamber is buffered by elastic (and local visco-elastic) deformations in the solid rocks bounding the chamber to the range between ambient and close to lithostatic values. In a chamber closed to fresh magma inputs, the decrease in volume due to such gentle volatile escape lowers the reference pressure. Bubbles flux from the lower layer induced by variation of the saturation level around stratification boundary may be efficient mechanism for the water transport between layers.     

Application of <Superscript>137</Superscript>Cs fingerprinting technique to interpreting sediment production records from reservoir deposits in a small catchment of the Hilly Loess Plateau,China

According to variations of 137Cs and clay contents, 44 flood couplets were identified in a profile of reservoir deposit with a vertical length of 28.12 m in the Yuntaishan Gully. Couplet 27 at the middle of the profile had the highest average 137Cs content of 12.65 Bq kg-1, which indicated the 1963s' deposits, then 137Cs content decreased both downward and upward in the profile. The second top and bottom couplets had average 137Cs contents of 2.15 Bq kg-1 and 0.92 Bq kg-1, respectively. By integrated analysis of reservoir construction and management history, variations of 137Cs contents over the profile, sediment yields of flood couplets and rainfall data during the period of 1958-1970, individual storms related to the flood couplets were identified. 44 floods with a total sediment yield of 2.36×104 m3 occurred and flood events in a year varied between 1 and 10 times during the period of 1960-1970. 7-10 flood events occurred during the wet period of 1961-1964 with very wet autumn, while only 1-2 events during the dry period of 1965-1969. Average annual specific sediment yield was 1.29×104 t km-2 a-1 for the Yuntaishan Gully during the period of 1960-1970, which was slightly higher than 1.11 ×104 t km-2 a-1 for the Upper Yanhe River Basin above the Ganguyi Hydrological Station and slightly lower than 1.40 ×104 t km-2 a-1 for the nearby Zhifang Gully during the same period. Annual specific sediment yields for the Yuntaishan Gully were correlated to the wet season's rainfalls well.     

Kinetic disequilibrium of C, He, Ar and carbon isotopes during degassing of mid-ocean ridge basalts

Vesicle characteristics (vesicularity, largest vesicle size, number of vesicles/cm²), CO₂-δ¹³C and CO₂-⁴He-⁴⁰Ar-⁴⁰Ar/³⁶Ar in vesicles and CO₂-δ¹³C in the glass have been measured in 19 tholeiitic basalt glasses from the Easter Microplate East Ridge (East Pacific Rise) collected at 3 different sites (26°S East Ridge, Pito Seamount and Pito Deep at 23°S).Carbon supersaturation values (C_melt/C_solubility) vary from 1.3 to 4.3. Carbon supersaturation values are strongly correlated with the number of vesicles/cm². There is also a correlation between number of vesicles/cm² and vesicle size. At the Pito Seamount site, there is a negative correlation between carbon supersaturation values and observed carbon isotope fractionation between CO₂ in vesicles and carbon dissolved in the glass (Δ¹³C_observed). High ⁴He/⁴⁰Ar* ratios in vesicles (from 49 to 190) are observed in both the most and least carbon supersaturated samples, while samples with intermediate carbon supersaturation have the lowest ⁴He/⁴⁰Ar* ratios (16±1). These correlations show that most quenched melts record different disequilibrium to equilibrium states during closed-system degassing.The samples showing the highest carbon supersaturation (4.3) have the highest ⁴He/⁴⁰Ar* (from 94 to 190). This observation shows for the first time that the ⁴He/⁴⁰Ar* ratio can be kinetically fractionated during incomplete degassing of magmas from the magma chamber to the seafloor. This result implies that high ⁴He/⁴⁰Ar* ratios are not a systematic indicator of open-system degassing (Rayleigh distillation) and that caution should be taken when using this ratio for any degassing correction.A two-stage degassing model, with the first stage being a closed-system degassing occurring between the source and the magma chamber, and the second stage of degassing (with a mode varying from open-system degassing to different degrees of kinetic closed-system degassing) taking place between the magma chamber and eruption on the seafloor, is the most appropriate to describe the degassing of MORB. Reconstructing initial carbon content of the magma prior to degassing and extrapolating the results to the entire ridge system results in a carbon flux of 1.6_-0.3^+0.6×10¹⁴ g/year. This value implies vigorous exchange of carbon between the mantle and the surface throughout geological times.     

Methane Oxidation in the Water Column of Xiangxi Bay,Three Gorges Reservoir

Four field campaigns are carried out to quantify the methane (CH₄) oxidation rate in Xiangxi Bay (XXB) of the Three Gorges Reservoir (TGR), China. The water depth of the sampling site varied from 13 to 30 m resulting from the water level fluctuation of the TGR. The CH₄ oxidation rates are measured in situ as the decline of dissolved CH₄ concentration versus time in incubated, and those rates. The CH₄ oxidation rates range from 1.18 × 10^?3 to 3.69 × 10^?3 µmol L^?1 h^?1, with higher values and stronger variation during summer. A static floating chamber method is used to measure CH₄ emitted to the atmosphere resulting in an annual mean flux of 4.79 µmol m^?2 h^?1. The CH₄ emission rate is significantly negatively correlated with the water level. The results show that a large fraction of CH₄ is consumed in the water column with a range of 28.97–55.90 µmol m^?2 h^?1, accounting for ≈69–98% of the total CH₄ input into the water column, and more than 90% is consumed outside the summer, when the water level is lowest. Water depth, which is dominated by water level of the TGR, is a potentially important driver for CH₄ oxidation and atmospheric emission in the tributary bay.     

Magma convection and mixing dynamics as a source of Ultra-Long-Period oscillations

Many volcanic eruptions are shortly preceded by injection of new magma into a pre-existing, shallow (<10 km) magma chamber, causing convection and mixing between the incoming and resident magmas. These processes may trigger dyke propagation and further magma rise, inducing long-term (days to months) volcano deformation, seismic swarms, gravity anomalies, and changes in the composition of volcanic plumes and fumaroles, eventually culminating in an eruption. Although new magma injection into shallow magma chambers can lead to hazardous event, such injection is still not systematically detected and recognized. Here, we present the results of numerical simulations of magma convection and mixing in geometrically complex magmatic systems, and describe the multiparametric dynamics associated with buoyant magma injection. Our results reveal unexpected pressure trends and pressure oscillations in the Ultra-Long-Period (ULP) range of minutes, related to the generation of discrete plumes of rising magma. Very long pressure oscillation wavelengths translate into comparably ULP ground displacements with amplitudes of order 10⁻⁴–10⁻² m. Thus, new magma injection into magma chambers beneath volcanoes can be revealed by ULP ground displacement measured at the surface.     

The 1975 sub-terminal lavas, mount etna: a case history of the formation of a compound lava field

The 1975 sub-terminal activity was characterised by low effusion rates (0.3–0.5 m³ s⁻¹) and the formation of a compound lava field composed of many thousands of flow units. Several boccas were active simultaneously and effusion rates from individual boccas varied from about 10⁻⁴ to 0.25 m³s⁻¹. The morphology of lava flows was determined by effusion rate (E): aa flows with well-developed channels and levees formed when E > 2 × 10⁻³ m³ s⁻¹, small pahoehoe flows formed when 2 × 10⁻³ m³ s⁻¹ >E > 5 > 10⁻⁴ m³ s⁻¹ and pahoehoe toes formed when E < 5 × 10⁻⁴ m³ s⁻¹. There was very little variation with time in the effusion temperature, composition or phenocryst content of the lava.New boccas were commonly formed at the fronts of mature lava flows which had either ceased to flow or were moving slowly. These secondary boccas developed when fluid lava in the interior of mature aa flows either found a weakness in the flow front or was exposed by avalanching of the moving flow front. The resulting release of fluid lava was accompanied by either partial drainage of the mature flow or by the formation of a lava tube in the parent flow. The temperature of the lava forming the new bocca decreased with increasing distance from the source bocca (0.035°C m⁻¹). It is demonstrated from the rate of temperature decrease and from theoretical considerations that many of the Etna lavas still contained a substantial proportion of uncooled material in their interior as they came to rest. The formation of secondary boccas is postulated to be one reason why direct measurements of effusion rates tend, in general, to overestimate the total effusion rates of sub-terminal Etna lava fields.     

Hydrothermal vent waters at 13°N on the East Pacific Rise: isotopic composition and gas concentration

Gas concentrations and isotopic compositions of water have been measured in hydrothermal waters from 13°N on the East Pacific Rise. In the most Mg-depleted samples ( 5 × 10⁻³ moles/kg) the gas concentrations are: 3–4.5 × 10⁻⁵ cm³ STP/kg helium, 0.62–1.24 cm³ STP/kg CH₄, 10.80–16.71 × 10⁻³ moles/kg CO₂. The samples contain large quantities (95–126 cm³/kg) of H₂ and some carbon monoxide (0.26–0.36 cm³/kg) which result from reaction with the titanium sampling bottles. δ¹³C in methane and CO₂ (−16.6 to −19.5 and −4.1 to −5.5 respectively) indicate temperatures between 475 and 550°C, whereas δ¹³C_CO is compatible with formation by reduction of CO₂ on Ti at 350°C close to the sampling temperature.³He/⁴He are very homogeneous at (7.5 ± 0.1)R_A(³He/⁴He = 1.0 × 10⁻⁵) and very similar to already published data as well as CH₄/³He ratios between 1.4 and 2.1 × 10⁶.¹⁸O and D in water show enrichments from 0.39 to 0.69‰ and from 0.62 to 1.49‰ respectively. These values correspond to W/R ratios of 0.4–7. The distinct¹⁸O enrichments indicate that the isotopic composition of the oceans is not completely buffered by the hydrothermal circulations. The³He-enthalpy relationship is discussed in terms of both hydrothermal heat flux and³He mantle flux.     

Urbanization effects on shoreline phytobenthos: a multiscale approach at lake extent

To understand how littoral biota respond to anthropogenic disturbances, limnologists seek to detect the scale at which patterns and processes occur. We conducted an extensive study on the shoreline phytobenthos of Lake Garda (Italy) with the following main objectives: (i) to examine the importance of urbanization for species distribution within a set of hierarchical spatial scales (10¹–10⁴ m), and then (ii) to test the spatio-temporal interactions on a reduced set of scales (10¹–10² m, and 10¹–10² days). Results showed that most of the variation in most abundant species and habitat characteristics occurred at the spatial scale of 10¹–10² m. Species richness was positively related with microheterogeneity, but the relationship occurred only at low urbanization and not at highly-urbanized sites where artificial shores were less heterogeneous. The similarity of species assemblages was regulated by two interacting processes, one operating at a fine spatial scale (10² m), reflecting the physical-habitat requirements of the species, and the other one operating at a broader scale (10⁴ m) in relation to the N–S nitrogen gradient. Overall, time explained 73 % of the total variation of species assemblages, space 7 %, and 20 % was explained by the interaction between space and time (the patch scale, 10s of m, and area scale, 100s of m, interacted with the finest temporal scale, 10s of days). This interaction might be explained by the process of species recruitment operating at different rates at the two spatial scales. Since the largest variation in species assemblages was at the temporal scale (due to the seasonal succession of phytobenthos), we recommend collecting at least one sample per season when monitoring littoral habitats.     

Eruptive versus non-eruptive behaviour of large calderas: the example of Campi Flegrei caldera (southern Italy)

Caldera eruptions are among the most hazardous of natural phenomena. Many calderas around the world are active and are characterised by recurrent uplift and subsidence periods due to the dynamics of their magma reservoirs. These periods of unrest are, in some cases, accompanied by eruptions. At Campi Flegrei caldera (CFc), which is an area characterised by very high volcanic risk, the recurrence of this behaviour has stimulated the study of the rock rheology around the magma chamber, in order to estimate the likelihood of an eruption. This study considers different scenarios of shallow crustal behaviour, taking into account the earlier models of CFc ground deformation and caldera eruptions, and including recent geophysical investigations of the area. A semi-quantitative evaluation of the different factors that lead to magma storage or to its eruption (such as magma chamber size, wall-rock viscosity, temperature, and regional tectonic strain rate) is reported here for elastic and viscoelastic conditions. Considering the large magmatic sources of the CFc ignimbrite eruptions (400–2,000 km³) and a wall-rock viscosity between 10¹⁸ and 10²⁰ Pa s, the conditions for eruptive failure are difficult to attain. Smaller source dimensions (a few cubic kilometres) promote the condition for fracture (eruption) rather than for the flow of wall rock. We also analyse the influence of the regional extensional stress regime on magma storage and eruptions, and the thermal stress as a possible source of caldera uplift. The present study also emphasises the difficulty of distinguishing eruption and non-eruption scenarios at CFc, since an unambiguous model that accounts for the rock rheology, magma-source dimensions and locations and regional stress field influences is still lacking.     

The Avellino plinian eruption of Somma-Vesuvius (3760 y.B.P.): the progressive evolution from magmatic to hydromagmatic style

A detailed stratigraphic analysis of the Avellino plinian deposit of the Somma-Vesuvius volcano shows a complicated eruptive sequence controlled by a combination of magmatic and hydromagmatic processes. The role of external water on the eruptive dynamics was most relevant in the very early phase of the eruption when the groundwater explosively interacted with a rising, gas-exolving magma body creating the first conduit. This phase generated pyroclastic surge and phreatoplinian deposits followed by a rapidly increasing discharge of a gas-rich, pure magmatic phase which erupted as the most violent plinian episode. This continuing plinian phase tapped the magma chamber, generating about 2.9 km³ of reverse-graded fallout pumice, more differentiated at the base and more primitive at the top (white and gray pumice). A giant, plinian column, rapidly grew up reaching a maximum height of 36 km.The progressive magma evacuation at a maximum discharge rate of 10⁸ kg/s that accompanied a decrease of magmatic volatile content in the lower primitive magma allowed external water to enter the magma chamber, resulting in a drastic change in the eruptive style and deposit type. Early wet hydromagmatic events were followed by dry ones and only a few, subordinated magmatic phases. A thick, impressive sequence of pyroclastic surge bedsets of over 430 km² in area with a total volume of about 1 km³ is the visible result of this hydromagmatic phase.     

Nonlinear dynamics of boussinesq convection in a deep rotating spherical shell-i

Abstract

Convection in a rotating spherical shell has wide application for understanding the dynamics of the atmospheres and interiors of many celestial bodies. In this paper we review linear results for convection in a shell of finite depth at substantial but not asymptotically large Taylor numbers, present nonlinear multimode calculations for similar conditions, and discuss the model and results in the context of the problem of solar convection and differential rotation. Detailed nonlinear calculations are presented for Taylor number T = 10⁵, Prandtl number P = 1, and Rayleigh number R between 1 |MX 10⁴ and 4 |MX 10⁴ (which is between about 4 and 16 times critical) for a shell of depth 20% of the outer radius. Sixteen longitudinal wave numbers are usually included (all even wave numbers m between 0 and 30) the amplitudes of which are computed on a staggered grid in the meridian plane.

The kinetic energy spectrum shows a peak in the wave number range m = 12–18 at R = 10⁴, which straddles the critical wave number m = 14 predicted by linear theory. These are modes which peak near the equator. The spectrum shows a second strong peak at m = 0, which represents the differential rotation driven by the peak convective modes. As R is increased, the amplitude of low wave numbers increases relative to high wave numbers as convection fills in in high and middle latitudes, and as the longitudinal scale of equatorial convection grows. By R = 3 |MX 10⁴, m = 8 is the peak convective mode. There is a clear minimum in the total kinetic energy at middle latitudes relative to low and high, well into the nonlinear regime, representing the continued dominance of equatorial and polar modes found in the linear case. The kinetic energy spectrum for m > 0 is maintained primarily by buoyancy work in each mode, but with substantial nonlinear transfer of kinetic energy from the peak modes to both lower and higher wave numbers.

For R = 1 to 2 |MX 10⁴, the differential rotation takes the form of an equatorial acceleration, with angular velocity generally decreasing with latitude away from the equator (as on the sun) and decreasing inwards. By R = 4 |MX 10⁴, this equatorial profile has completely reversed, with angular velocity increasing with depth and latitude. Also, a polar vortex which has positive rotation relative to the reference frame (no evidence of which has been seen on the sun) builds up as soon as polar modes become important. Meridional circulation is quite weak relative to differential rotation at R = 10⁴, but grows relative to it as R is increased. This circulation takes the farm of a single cell of large latitudinal extent in equatorial regions, with upward flow near the equator, together with a series of narrower cells in high latitudes. It is maintained primarily by axisymmetric buoyancy forces. The differential rotation is maintained at all R primarily by Reynolds stresses, rather than meridional circulation. Angular momentum transport toward the equator for R = 1–2 |MX 10⁴ maintains the equatorial acceleration while radially inward transport maintains the opposite profile at R = 4 |MX 10⁴.

The total heat flux out the top of the convective shell always shows two peaks for the range of R studied, one at the equator and the other near the poles (no significant variation with latitude is seen on the sun), while heat flux in at the bottom shows only a polar peak at large R. The meridional circulation and convective cells transport heat toward the equator to maintain this difference.

The helicity of the convection plus the differential rotation produced by it suggest the system may be capable of driving a field reversing dynamo, but the toroidal field may migrate with lime in each cycle toward the poles and equator, rather than just toward the equator as apparently occurs on the sun.

We finally outline additions to the physics of the model to make it more realistic for solar application.     

The intensity of plinian eruptions

Peak intensities (magma discharge rate) of 45 Pleistocene and Holocene plinian eruptions have been inferred from lithic dispersal patterns by using a theoretical model of pyroclast fallout from eruption columns. Values range over three orders of magnitude from 1.6 × 10⁶ to 1.1 × 10⁹ kg/s. Magnitudes (total erupted mass) also vary over about three orders of magnitude from 2.0 × 10¹¹ to 6.8 × 10¹⁴ kg and include several large ignimbrite-forming events with associated caldera formation. Intensity is found to be positively correlated with the magnitude when total erupted mass (tephra fall, surges and pyroclastic flows) is considered. Initial plinian fall phases with intensities in excess of 2.0 × 10⁸ kg/s typically herald the onset of major pyroclastic flow generation and subsequent caldera collapse. During eruptions of large magnitude, the transition to pyroclastic flows is likely to be the result of high intensity, whereas the generation of pyroclastic flows in small magnitude eruptions may occur more often by reduction of magmatic volatile content or some transient change in magma properties. The correlation between plinian fall intensity and total magnitude suggests that the rate of magma discharge is related to the size of the chamber being tapped. A simple model is presented to account for the variation in intensity by progressive enlargement of conduits and vents and excess pressure at the chamber roof caused by buoyant forces acting on the chamber as it resides in the crust. Both processes are fundamentally linked to the absolute size of the pre-eruption reservoir. The data suggest that sustained eruption column heights (i.e. magma discharge rates) are indicators of eventual eruption magnitude, and perhaps eruptive style, and thus are key parameters to monitor in order to assess the temporal evolution of plinian eruptions.     

Heterogeneous loss of OH on NaCl and NH4 NO 3 at tropospheric temperatures

The uptake coefficients () for OH radicals on some dry salts of tropospheric interest (NaCl and NH₄NO₃) have been investigated as a function of temperature using the flow tube technique combined with an EPR spectrometer as a detection method. The temperature dependence of -values measured over the temperature range 245–340 K can be expressed in Arrhenius form: ^OH_NaCl=(1.2±0.7)×10⁻⁵exp[(1750±200)/T] and ^OH_NH4NO3=(1.4±0.5)×10⁻⁴exp[(1000±100)/T]. These Arrhenius expressions lead to very similar -values (\approx4×10⁻³) for both salts studied at 300 K. It is shown that the heterogeneous OH sinks on solids aerosol play a very minor role in tropospheric chemistry in comparison with the homogeneous sinks.