The influence of viscous and latent heating on crystal-rich magma flow in a conduit

The flow dynamics of crystal-rich high-viscosity magma is likely to be strongly influenced by viscous and latent heat release. Viscous heating is observed to play an important role in the dynamics of fluids with temperature-dependent viscosities. The growth of microlite crystals and the accompanying release of latent heat should play a similar role in raising fluid temperatures. Earlier models of viscous heating in magmas have shown the potential for unstable (thermal runaway) flow as described by a Gruntfest number, using an Arrhenius temperature dependence for the viscosity, but have not considered crystal growth or latent heating. We present a theoretical model for magma flow in an axisymmetric conduit and consider both heating effects using Finite Element Method techniques. We consider a constant mass flux in a 1-D infinitesimal conduit segment with isothermal and adiabatic boundary conditions and Newtonian and non-Newtonian magma flow properties. We find that the growth of crystals acts to stabilize the flow field and make the magma less likely to experience a thermal runaway. The additional heating influences crystal growth and can counteract supercooling from degassing-induced crystallization and drive the residual melt composition back towards the liquidus temperature. We illustrate the models with results generated using parameters appropriate for the andesite lava dome-forming eruption at Soufriere Hills Volcano, Montserrat. These results emphasize the radial variability of the magma. Both viscous and latent heating effects are shown to be capable of playing a significant role in the eruption dynamics of Soufriere Hills Volcano. Latent heating is a factor in the top two kilometres of the conduit and may be responsible for relatively short-term (days) transients. Viscous heating is less restricted spatially, but because thermal runaway requires periods of hundreds of days to be achieved, the process is likely to be interrupted. Our models show that thermal evolution of the conduit walls could lead to an increase in the effective diameter of flow and an increase in flux at constant magma pressure.

Pb-210-Ra-226 and Ra-228-Th-232 systematics in young arc lavas: implications for magma degassing and ascent rates

New data show that island arc rocks have (Pb-210/Ra-226)(o) ratios which range from as low as 0.24 up to 2.88. In contrast, (Ra-22S/Th-232) appears always within error of I suggesting that the large Ra-226-excesses observed in arc rocks were generated more than 30 years ago. This places a maximum estimate on melt ascent velocities of around 4000 m/year and provides further confidence that the Ra-226 excesses reflect deep (source) processes rather than shallow level alteration or seawater contamination. Conversely, partial melting must have occurred more than 30 years prior to eruption. The Pb-210 deficits are most readily explained by protracted magma degassing. Using published numerical models, the data suggest that degassing occurred continuously for periods up to several decades just prior to eruption but no link with eruption periodicity was found. Longer periods are required if degassing is discontinuous, less than 100% efficient or if magma is recharged or stored after degassing. The long durations suggest much of this degassing occurs at depth with implications for the formation of hydrothermal and copper-porphyry systems. A suite of lavas erupted in 1985-1986 from Sangeang Api volcano in the Sunda arc are characterised by deficits of Pb-210 relative to Ra-226 from which 6-8 years of continuous Rn-222 degassing would be inferred from recent numerical models. These data also form a linear (Pb-210)/Pb-(Ra-226)/Pb array which might be interpreted as a 71-year isochron. However, the array passes through the origin suggesting displacement downwards from the equiline in response to degassing and so the slope of the array is inferred not to have any age significance. Simple modelling shows that the range of (Ra-226)/Pb ratios requires thousands of years to develop consistent with differentiation occurring in response to cooling at the base of the crust. Thus, degassing post-dated, and was not responsible for magma differentiation. The formation, migration and extraction of gas bubbles must be extremely efficient in mafic magma whereas the higher viscosity of more siliceous magmas retards the process and can lead to Pb-210 excesses. A possible negative correlation between (Pb-210/Ra-226)(o) and SO2 emission rate requires further testing but may have implications for future eruptions. (C) 2004 Elsevier B.V. All rights reserved.

Tracing pre-eruptive magma degassing using (Pb-210/Ra-226) disequilibria in the volcanic deposits of the 1980-1986 eruption of Mount St. Helens

Disequilibria between Pb-210 and Ra-226 can be used to trace magma degassing, because the intermediate nuclides, particularly Rn-222, are volatile. Products of the 1980-1986 eruptions of Mount St. Helens have been analysed for (Pb-210/Ra-226). Both excesses and deficits of Pb-210 are encountered suggesting rapid gas transfer. The time scale of diffuse, non-eruptive gas escape prior to 1980 as documented by Pb-210 deficits is on the order of a decade using the model developed by Gauthier and Condomines (Earth Planet. Sci. Lett. 172 (1999) 111-126) for a non-renewed magma chamber and efficient Rn removal. The time required to build-up Pb-210 excess is much shorter (months) as can be observed from steady increases of (Pb-210/Ra-226) with time during 1980-1982. The formation of Pb-210 excess requires both rapid gas transport through the magma and periodic blocking of gas escape routes. Superposed on this time trend is the natural variability of (Pb-210/Ra-226) in a single eruption caused by tapping magma from various depths. The two time scales of gas transport, to create both Pb-210 deficits and Pb-210 excesses, cannot be reconciled in a single event. Rather Pb-210 deficits are associated with pre-eruptive diffuse degassing, while Pb-210 excesses document the more vigorous degassing associated with eruption and recharge of the system. (c) 2006 Elsevier B.V. All rights reserved.

The influence of different artificial soil types on the acute toxicity of carbendazim to the earthworm Eisenia fetida in laboratory toxicity tests

Field populations of earthworms have shown a varied response in mortality to the fungicide carbendazim, the toxic reference substance used in agrochemical field trials. The aim of this study was to determine the influence of soil conditions as a potential cause of this variation. Laboratory acute toxicity tests were conducted using a range of artificial soils with varying soil components (organic matter, clay, pH and moisture). Batch adsorption/desorption studies were run to determine the influence of the soil properties on carbendazim behaviour. Adsorption was shown to be correlated with organic matter content and pH and this in turn could be linked to Eisenia fetida mortality, with lower mortality occurring with increased adsorption. Overall while E.fetida mortality did vary significantly between several of the soils the calculated LC50 values in the different soils did not cover a wide range (6.04-16.00 mg kg(-1)), showing that under these laboratory conditions soil components did not greatly influence carbendazim toxicity to E.fietida. (c) 2007 Elsevier Masson SAS. All rights reserved.

Pb-210-Ra-226 and Ra-228-Th-232 systematics in young arc lavas: implications for magma degassing and ascent rates

New data show that island arc rocks have (Pb-210/Ra-226)(o) ratios which range from as low as 0.24 up to 2.88. In contrast, (Ra-22S/Th-232) appears always within error of I suggesting that the large Ra-226-excesses observed in arc rocks were generated more than 30 years ago. This places a maximum estimate on melt ascent velocities of around 4000 m/year and provides further confidence that the Ra-226 excesses reflect deep (source) processes rather than shallow level alteration or seawater contamination. Conversely, partial melting must have occurred more than 30 years prior to eruption. The Pb-210 deficits are most readily explained by protracted magma degassing. Using published numerical models, the data suggest that degassing occurred continuously for periods up to several decades just prior to eruption but no link with eruption periodicity was found. Longer periods are required if degassing is discontinuous, less than 100% efficient or if magma is recharged or stored after degassing. The long durations suggest much of this degassing occurs at depth with implications for the formation of hydrothermal and copper-porphyry systems. A suite of lavas erupted in 1985-1986 from Sangeang Api volcano in the Sunda arc are characterised by deficits of Pb-210 relative to Ra-226 from which 6-8 years of continuous Rn-222 degassing would be inferred from recent numerical models. These data also form a linear (Pb-210)/Pb-(Ra-226)/Pb array which might be interpreted as a 71-year isochron. However, the array passes through the origin suggesting displacement downwards from the equiline in response to degassing and so the slope of the array is inferred not to have any age significance. Simple modelling shows that the range of (Ra-226)/Pb ratios requires thousands of years to develop consistent with differentiation occurring in response to cooling at the base of the crust. Thus, degassing post-dated, and was not responsible for magma differentiation. The formation, migration and extraction of gas bubbles must be extremely efficient in mafic magma whereas the higher viscosity of more siliceous magmas retards the process and can lead to Pb-210 excesses. A possible negative correlation between (Pb-210/Ra-226)(o) and SO2 emission rate requires further testing but may have implications for future eruptions. (C) 2004 Elsevier B.V. All rights reserved.

The magma budget of Volcan Arenal, Costa Rica from 1980 to 2004

Using topographic data collected by radar interferometry, stereo-photogrammetry, and field survey we have measured the changing surface of Volcan Arenal in Costa Rica over the period from 1980 to 2004. During this time this young volcano has mainly effused basaltic andesite lava, continuing the activity that began in 1968. Explosive products form only a few percent of the volumetric output. We have calculated digital elevation models for the years 1961, 1988 and 1997 and modified existing models for 2000 and 2004. From these we have estimated the volume of lava effused and coupled this with the data presented by an earlier study for 1968-1980. We find that a dense rock equivalent volume of 551 M m(3) was effused from 1968 to 2004. The dense rock equivalent effusion rate fell from about 2 m(3) s(-1) to about 0.1-0.2 m(3) s(-1) over the same period, with an average rate of about 0.5 m(3) s(-1). Between 1980 and 2004, the average effusion rate was 0.36 m(3) s(-1), a similar rate to that measured between 1974 and 1980. There have been two significant deviations from this long-term rate. The effusion rate increased from 1984 to 1991, at the same time as explosivity increased. After a period of moderate effusion rates in the 1990s, the rate fell to lower levels around 1999. (c) 2006 Elsevier B.V. All rights reserved.

Ground deformation at Soufriere Hills Volcano, Montserrat during 1998-2000 measured by radar interferometry and GPS

We examine the motion of the ground surface on the Soufriere Hills Volcano, Montserrat between 1998 and 2000 using radar interferometry (InSAR). To minimise the effects of variable atmospheric water vapour on the InSAR measurements we use independently-derived measurements of the radar path delay from six continuous GPS receivers. The surfaces providing a measurable inter-ferometric signal are those on pyroclastic flow deposits, mainly emplaced in 1997. Three types of surface motion can be discriminated. Firstly, the surfaces of thick, valley-filling deposits subsided at rates of 150-120 mm/year in the year after emplacement to 50-30 mm/year two years later. This must be due to contraction and settling effects during cooling. The second type is the near-field motion localised within about one kilometre of the dome. Both subsidence and uplift events are seen and though the former could be due to surface gravitational effects, the latter may reflect shallow (< 1 km) pressurisation effects within the conduit/dome. Far-field motions of the surface away from the deeply buried valleys are interpreted as crustal strains. Because the flux of magma to the surface stopped from March 1998 to November 1999 and then resumed from November 1999 through 2000, we use InSAR data from these two periods to test the crustal strain behaviour of three models of magma supply: open, depleting and unbalanced. The InSAR observations of strain gradients of 75-80 mm/year/krn uplift during the period of quiescence on the western side of the volcano are consistent with an unbalanced model in which magma supply into a crustal magma chamber continues during quiescence, raising chamber pressure that is then released upon resumption of effusion. GPS motion vectors agree qualitatively with the InSAR displacements but are of smaller magnitude. The discrepancy may be due to inaccurate compensation for atmospheric delays in the InSAR data. (c) 2005 Elsevier B.V. All rights reserved.

Factors affecting distribution of vegetation types on abandoned cropland in the hilly-gullied Loess Plateau region of China

A study was conducted in the forest-steppe region of the Loess Plateau to provide insight into the factors affecting the process of vegetation establishment, and to provide recommendations for the selection of indigenous species in order to speed up the succession process and to allow the establishment of vegetation more resistant to soil erosion. Four distinctive vegetation types were identified, and their distribution was affected not only by the time since abandonment but also by other environmental factors, mainly soil water and total P in the upper soil layers. One of the vegetation types, dominated by Artemisia scoparia, formed the early successional stage after abandonment while the other three types formed later successional stages with their distribution determined by the soil water content and total P. It can be concluded that the selection of appropriate species for introduction to accelerate succession should be determined by the local conditions and especially the total P concentration and soil water content.

Bayesian clustering of farm types using the mixtures model

A Bayesian method of classifying observations that are assumed to come from a number of distinct subpopulations is outlined. The method is illustrated with simulated data and applied to the classification of farms according to their level and variability of income. The resultant classification shows a greater diversity of technical charactersitics within farm types than is conventionally the case. The range of mean farm income between groups in the new classification is wider than that of the conventional method and the variability of income within groups is narrower. Results show that the highest income group in 2000 included large specialist dairy farmers and pig and poultry producers, whilst in 2001 it included large and small specialist dairy farms and large mixed dairy and arable farms. In both years the lowest income group is dominated by non-milk producing livestock farms.

Soil disturbed using a strip tillage implement on a range of soil types and the effects on sugar beet establishment

The area of soil disturbed using a single tine is well documented. However, modern strip tillage implements using a tine and disc design have not been assessed in the UK or in mainland Europe. Using a strip tillage implement has potential benefits for European agriculture where economic returns and sustainability are key issues. Using a strip tillage system a narrow zone is cultivated leaving most of the straw residue on the soil surface. Small field plot experiments were undertaken on three soil types and the operating parameters of forward speed, tine depth and tine design were investigated together with measurements of seedbed tilth and crop emergence. The type of tine used was found to be the primary factor in achieving the required volume of disturbance within a narrow zone whilst maintaining an area of undisturbed soil with straw residue on the surface. The winged tine produced greater disturbance at a given depth compared with the knife tine. Increasing forward speed did not consistently increase the volume of disturbance. In a sandy clay loam the tilth created and emergence of sugar beet by strip tillage and ploughing were similar but on a sandy loam the strip tillage treatments generally gave a finer tilth but poorer emergence particularly at greater working depth.

Comparison of the euryarchaeal microbial community in guts and food-soil of the soil-feeding termite Cubitermes fungifaber across different soil types

Termites are an important component of tropical soil communities and have a significant affect on the structure and nutrient content of soil. Digestion in termites is related to gut structure, gut physico-chemical conditions and gut symbiotic microbiota. Here we describe the use of 16S rRNA gene sequencing and Terminal-restriction Fragment Length Polymorphism (T-RFLP) analysis to examine methanogenic Archaea (MA) in the guts and food-soil of the soil-feeder Cubitermes fungifaber Sjostedt across a range of soil types. If they are strictly vertically inherited, then MA in guts should be the same in all individuals even if the soils differ across sites. In contrast, gut MA should reflect what is present in soil if populations are merely a reflection of what is ingested as the insects forage. We show clear differences between the euryarchaeal communities in termite guts and in food-soils from five different sites. Analysis of 16S rRNA gene clones indicated little overlap between the gut and soil communities. Gut clones were related to a termite-derived Methanomicrobiales cluster, to Methanobrevibacter and, surprisingly, to the haloalkaliphile Natronococcus. Soil clones clustered with Methanosarcina, Methanomicrococcus or Rice Cluster I. T-RFLP analysis indicated that the archaeal communities in the soil samples differed from site to site, whereas those in termite guts were similar between sites. There was some overlap between the gut and soil communities but these may represent transient populations in either guts or soil. Our data does not support the hypothesis that termite gut MA are derived from their food soil but also does not support a purely vertical transmission of gut microflora.

Comparison of Euryarchaea strains in the guts and food-soil of the soil-feeding termite Cubitermes fungifaber across different soil types

Termites are an important component of tropical soil communities and have a significant effect on the structure and nutrient content of soil. Digestion in termites is related to gut structure, gut physicochemical conditions, and gut symbiotic microbiota. Here we describe the use of 16S rRNA gene sequencing and terminal-restriction fragment length polymorphism (T-RFLP) analysis to examine methanogenic archaea (MA) in the guts and food-soil of the soil-feeder Cubitermes fungifaber Sjostedt across a range of soil types. If these MA are strictly vertically inherited, then the MA in guts should be the same in all individuals even if the soils differ across sites. In contrast, gut MA should reflect what is present in soil if populations are merely a reflection of what is ingested as the insects forage. We show clear differences between the euryarchaeal communities in termite guts and in food-soils from five different sites. Analysis of 16S rRNA gene clones indicated little overlap between the gut and soil communities. Gut clones were related to a termite-derived Methanomicrobiales cluster, to Methanobrevibacter and, surprisingly, to the haloalkaliphile Natronococcus. Soil clones clustered with Methanosarcina, Methanomicrococcus, or rice cluster I. T-RFLP analysis indicated that the archaeal communities in the soil samples differed from site to site, whereas those in termite guts were similar between sites. There was some overlap between the gut and soil communities, but these may represent transient populations in either guts or soil. Our data do not support the hypothesis that termite gut MA are derived from their food-soil but also do not support a purely vertical transmission of gut microflora.