Evaluating the structure and magnitude of the ash plume during the initial phase of the 2010 Eyjafjallajökull eruption using lidar observations and NAME simulations

The Eyjafjallajökull volcano in Iceland erupted explosively on 14 April 2010, emitting a plume of ash into the atmosphere. The ash was transported from Iceland toward Europe where mostly cloud-free skies allowed ground-based lidars at Chilbolton in England and Leipzig in Germany to estimate the mass concentration in the ash cloud as it passed overhead. The UK Met Office's Numerical Atmospheric-dispersion Modeling Environment (NAME) has been used to simulate the evolution of the ash cloud from the Eyjafjallajökull volcano during the initial phase of the ash emissions, 14–16 April 2010. NAME captures the timing and sloped structure of the ash layer observed over Leipzig, close to the central axis of the ash cloud. Relatively small errors in the ash cloud position, probably caused by the cumulative effect of errors in the driving meteorology en route, result in a timing error at distances far from the central axis of the ash cloud. Taking the timing error into account, NAME is able to capture the sloped ash layer over the UK. Comparison of the lidar observations and NAME simulations has allowed an estimation of the plume height time series to be made. It is necessary to include in the model input the large variations in plume height in order to accurately predict the ash cloud structure at long range. Quantitative comparison with the mass concentrations at Leipzig and Chilbolton suggest that around 3% of the total emitted mass is transported as far as these sites by small (<100 μm diameter) ash particles.

On the entrainment assumption in Schatzmann's integral plume model

The behaviour of stationary, non-passive plumes can be simulated in a reasonably simple and accurate way by integral models. One of the key requirements of these models, but also one of their less well-founded aspects, is the entrainment assumption, which parameterizes turbulent mixing between the plume and the environment. The entrainment assumption developed by Schatzmann and adjusted to a set of experimental results requires four constants and an ad hoc hypothesis to eliminate undesirable terms. With this assumption, Schatzmann’s model exhibits numerical instability for certain cases of plumes with small velocity excesses, due to very fast radius growth. The purpose of this paper is to present an alternative entrainment assumption based on a first-order turbulence closure, which only requires two adjustable constants and seems to solve this problem. The asymptotic behaviour of the new formulation is studied and compared to previous ones. The validation tests presented by Schatzmann are repeated and it is found that the new formulation not only eliminates numerical instability but also predicts more plausible growth rates for jets in co-flowing streams.

Ice Shelf Water plume flow beneath Filchner-Ronne Ice Shelf, Antarctica

[1] A two-dimensional plume model is used to study the interaction between Filchner-Ronne Ice Shelf, Antarctica and its underlying ocean cavity. Ice Shelf Water (ISW) plumes are initiated by the freshwater released from a melting ice shelf and, if they rise, may become supercooled and deposit marine ice due to the pressure increase in the in situ freezing temperature. The aim of this modeling study is to determine the origin of the thick accretions of marine ice at the base of Filchner-Ronne Ice Shelf and thus improve our understanding of ISW flow paths. The model domain is defined from measurements of ice shelf draft, and from this ISW the model is able to predict plumes that exit the cavity in the correct locations. The modeled plumes also produce basal freezing rates that account for measured marine ice thicknesses in the western part of Ronne Ice Shelf. We find that the freezing rate and plume properties are significantly influenced by the confluence of plumes from different meltwater sources. We are less successful in matching observations of marine ice under the rest of Filchner-Ronne Ice Shelf, which we attribute primarily to this model’s neglect of circulations in the ocean outside the plume.

The effect of a new drag-law parameterization on ice shelf water plume dynamics

A drag law accounting for Ekman rotation adjacent to a flat, horizontal bou ndary is proposed for use in a plume model that is written in terms of the depth-mean velocity. The drag l aw contains a variable turning angle between the mean velocity and the drag imposed by the turbulent bound ary layer. The effect of the variable turning angle in the drag law is studied for a plume of ice shelf wat er (ISW) ascending and turning beneath an Antarctic ice shelf with draft decreasing away from the groundi ng line. As the ISW plume ascends the sloping ice shelf–ocean boundary, it can melt the ice shelf, wh ich alters the buoyancy forcing driving the plume motion. Under these conditions, the typical turning ang le is of order 10° over most of the plume area for a range of drag coefficients (the minus sign arises for th e Southern Hemisphere). The rotation of the drag with respect to the mean velocity is found to be signifi cant if the drag coefficient exceeds 0.003; in this case the plume body propagates farther along and across the b ase of the ice shelf than a plume with the standard quadratic drag law with no turning angle.

Validity of the ice shelf water plume concept beneath Filchner-Ronne Ice Shelf

In winter, brine rejection from sea ice formation and export in the Weddell Sea, offshore of Filchner-Ronne Ice Shelf (FRIS), leads to the formation of High Salinity Shelf Water (HSSW). This dense water mass enters the cavity beneath FRIS by sinking southward down the sloping continental shelf towards the grounding line. Melting occurs when the HSSW encounters the ice shelf, and the meltwater released cools and freshens the HSSW to form a water mass known as Ice Shelf Water (ISW). If this ISW rises, the ‘ice pump’ is initiated (Lewis and Perkin, 1986), whereby the ascending ISW becomes supercooled and deposits marine ice at shallower locations due to the pressure increase in the in-situ freezing temperature. Sandh¨ager et al. (2004) were able to infer the thickness patterns of marine ice deposits at the base of FRIS (figure 1), so the primary aim of this work is to try to understand the ocean flows that determine these patterns. The plume model we use to investigate ISW flow is described fully by Holland and Feltham (accepted) so only a relatively brief outline is presented here. The plume is simulated by combining a parameterisation of ice shelf basal interaction and a multiplesize- class frazil dynamics model with an unsteady, depth-averaged reduced-gravity plume model. In the model an active region of ISW evolves above and within an expanse of stagnant ambient fluid, which is considered to be ice-free and has fixed profiles of temperature and salinity. The two main assumptions of the model are that there is a well-mixed layer underneath the ice shelf and that the ambient fluid outside the plume is stagnant with fixed properties. The topography of the ice shelf that the plume flows beneath is set to the FRIS ice shelf draft calculated by Sandh¨ager et al. (2004) masked with the grounding line from the Antarctic Digital Database (ADD Consortium, 2002). To initiate the plumes, we assume that the intrusion of dense HSSW initially causes melting at the points on the grounding line where the glaciological tributaries feeding FRIS go afloat.

Observations on the flight paths of the day-flying moth Virbia lamae during periods of mate location: do males have a strategy for contacting the pheromone plume?

1. To maximize the probability of rapid contact with a female’s pheromone plume, the trajectories of male foraging flights might be expected to be directed with respect to wind flow and also to be energetically efficient. 2. Flights directed either upwind, downwind, or crosswind have been proposed as optimal strategies for rapid and/or energetically efficient plume contact. Other possible strategies are random and Lévy walks, which have trajectories and turn frequencies that are not dictated by the direction of wind flow. 3. The planar flight paths of males of the day-active moth Virbia lamae were recorded during the customary time of its sexual activity. 4. We found no directional preference in these foraging flights with respect to the direction of contemporaneous wind flow, but, because crosswind encompasses twice the possible orientations of either upwind or downwind, a random orientation is in effect a de facto crosswind strategy. 5. A crosswind preference should be favoured when the plume extends farther downwind than crosswind, and this strategy is realized by V. lamae males by a random orientation of their trajectories with respect to current wind direction

The plume also rises: trajectories of pheromone plumes issuing from point sources in an orchard canopy at night

The trajectories of pheromone plumes in canopied habitats, such as orchards, have been little studied. We documented the capture of male navel orangeworm moths, Amyelois transitella, in female-baited traps positioned at 5 levels, from ground level to the canopy top, at approximately 6 m above ground, in almond orchards. Males were captured in similar proportions at all levels, suggesting that they do not favor a particular height during ranging flight. A 3-D sonic anemometer was used to establish patterns of wind flow and temperature at 6 heights from 2.08 to 6.65 m in an almond orchard with a 5 m high canopy, every 3 h over 72 h. The horizontal velocity of wind flow was highest above the canopy, where its directionality also was the most consistent. During the time of A. transitella mating (0300–0600), there was a net vertical displacement upward. Vertical buoyancy combined with only minor reductions in the distance that plumes will travel in the lower compared to the upper canopy suggest that the optimal height for release of pheromone from high-release-rate sources, such as aerosol dispensers (“puffers”), that are deployed at low densities (e.g., 3 per ha.) would be at mid or low in the canopy, thereby facilitating dispersion of disruptant throughout the canopy. Optimal placement of aerosol dispensers will vary with the behavioral ecology of the target pest; however, our results suggest that current protocols, which generally propose dispenser placement in the upper third of the canopy, should be reevaluated.

Transport of the smoke plume from Chiado’s fire in Lisbon (Portugal) sensed by atmospheric electric field measurements

The Chiado’s fire that affected the city centre of Lisbon (Portugal) occurred on 25th August 1988 and had a significant human and environmental impact. This fire was considered the most significant hazard to have occurred in Lisbon city centre after the major earthquake of 1755. A clear signature of this fire is found in the atmospheric electric field data recorded at Portela meteorological station about 8 km NE from the site where the fire started at Chiado. Measurements were made using a Benndorf electrograph with a probe at 1 m height. The atmospheric electric field reached 510 V/m when the wind direction was coming from SW to NE, favourable to the transport of the smoke plume from Chiado to Portela. Such observations agree with predictions using Hysplit air mass trajectory modelling and have been used to estimate the smoke concentration to be ~0.4 mg/m3. It is demonstrated that atmospheric electric field measurements were therefore extremely sensitive to Chiado’s fire. This result is of particular current interest in using networks of atmospheric electric field sensors to complement existing optical and meteorological observations for fire monitoring.

The Coronation mantle and the Westminster sanctuary pavement

The Great Sanctuary pavement at Westminster Abbey, on which the kings of England were crowned, features a design of rotated squares, around a central roundel of unusual Egyptian alabaster, set within a circle of bright blue glass. The roundel’s visual resemblance to the world and the firmament is underlined by a surrounding inscription, identifying it as the world, the macrocosm. This paper proposes that the design and meaning of the Great Sanctuary pavement reflect the square macrocosmic mantle in which the English king was wrapped during the coronation ceremony. This paper underlines the special nature of the English coronation mantle within the context of other royal and imperial cosmic mantles. Its square form was unique, and it featured in the English coronation liturgies in a way that European mantles did not. The English liturgies stressed its role as a representation of the world or macrocosm

Upper-mantle seismic anisotropy from SKS splitting in the South American stable platform: A test of asthenospheric flow models beneath the lithosphere

Upper-mantle seismic anisotropy has been extensively used to infer both present and past deformation processes at lithospheric and asthenospheric depths. Analysis of shear-wave splitting (mainly from core-refracted SKS phases) provides information regarding upper-mantle anisotropy. We present average measurements of fast-polarization directions at 21 new sites in poorly sampled regions of intra-plate South America, such as northern and northeastern Brazil. Despite sparse data coverage for the South American stable platform, consistent orientations are observed over hundreds of kilometers. Over most of the continent, the fast-polarization direction tends to be close to the absolute plate motion direction given by the hotspot reference model HS3-NUVEL-1A. A previous global comparison of the SKS fast-polarization directions with flow models of the upper mantle showed relatively poor correlation on the continents, which was interpreted as evidence for a large contribution of ""frozen"" anisotropy in the lithosphere. For the South American plate, our data indicate that one of the reasons for the poor correlation may have been the relatively coarse model of lithospheric thicknesses. We suggest that improved models of upper-mantle flow that are based on more detailed lithospheric thicknesses in South America may help to explain most of the observed anisotropy patterns.

Petrogenesis of Early Cretaceous silicic volcanism in SE Uruguay: The role of mantle and crustal sources

Early Cretaceous (similar to 129 Ma) silicic rocks crop out in SE Uruguay between the Laguna Merin and Santa Lucia basins in the Lascano, Sierra Sao Miguel. Salamanca and Minas areas They are mostly rhyolites with minor quartz-trachytes and are nearly contemporaneous with the Parana-Etendeka igneous province and with the first stages of South Atlantic Ocean opening A strong geochemical variability (particularly evident from Rb/Nb, Nb/Y trace element ratios) and a wide range of Sr-Nd isotopic ratios ((143)Nd/(144)Nd((129)) = 0.51178-0.51209, (87)Sr/(86)Sr((129)) = 0.70840-0.72417) characterize these rocks Geochemistry allows to distiniguish two compositional groups, corresponding to the north-eastern (Lascano and Sierra Sao Miguel, emplaced on the Neo-Proterozoic southern sector of the Dom Feliciano mobile belt) and south-eastern localities (Salamanca, Minas, emplace on the much older (Archean) Nico Perez teriane or on the boundary between the Dom Feliciano and Nico Perez termites) These compositional differences between the two groups are explained by variable mantle source and crust contributions. The origin of the silicic magmas is best explained by complex processes involving assimilation and fractional crystallization and mixing of a basaltic magma with upper crustal lithologies, for Lascano and Sierra Sao Miguel rhyolites. In the Salamanea and Minas rocks genesis, a stronger contribution from lower crust is indicated.

Age and geochemistry of mantle peridotites and diorite dykes from the Baldissero body: Insights into the Paleozoic-Mesozoic evolution of the Southern Alps

Trace element and isotopic data obtained for mantle spinel Iherzolites and diorite dykes from the Baldissero massif (Ivrea-Verbano Zone, Western Italy) provide new, valuable constraints on the petrologic and geodynamic evolution of the Southern Alps in Paleozoic to Mesozoic times. Whole rock and mineral chemistry indicates that Baldissero Iherzolites can be regarded as refractory mantle residues following limited melt extraction. In particular, the Light Rare Earth Elements (LREE)-depleted and fractionated compositions of whole rock and clinopyroxene closely match modelling results for refractory residues after low degrees (similar to 4-5%) of near-fractional melting of depleted mantle, possibly under garnet-facies conditions. Following this, the peridotite sequence experienced subsolidus re-equilibration at lithospheric spinel-facies conditions and intrusion of several generations of dykes. However, Iherzolites far from dykes show very modest metasomatic changes, as evidenced by the crystallisation of accessory titanian pargasite and the occurrence of very slight enrichments in highly incompatible trace elements (e.g. Nb). The Re-Os data for Iherzolites far from the dykes yield a 376 Ma (Upper Devonian) model age that is considered to record a partial melting event related to the Variscan orogenic cycle s.l. Dioritic dykes cutting the mantle sequence have whole rock, clinopyroxene and plagioclase characterised by high radiogenic Nd and low radiogenic Sr, which point to a depleted to slightly enriched mantle source. Whole rock and mafic phases of diorites have high Mg# values that positively correlate with the incompatible trace element concentrations. The peridotite at the dyke contact is enriched in orthopyroxene, iron and incompatible trace elements with respect to the Iherzolites far from dykes. Numerical simulations indicate that the geochemical characteristics of the diorites can be explained by flow of a hydrous, silica-saturated melt accompanied by reaction with the ambient peridotite and fractional crystallisation. The composition of the more primitive melts calculated in equilibrium with the diorite minerals show tholeiitic to transitional affinity. Internal Sm-Nd, three-point isochrons obtained for two dykes suggest an Upper Triassic-Lower Jurassic emplacement age (from 204 31 to 198 29 Ma). Mesozoic igneous events are unknown in the southern Ivrea-Verbano Zone (IVZ), but the intrusion of hydrous melts, mostly silica-saturated, have been well documented in the Finero region, i.e. the northernmost part of IVZ and Triassic magmatism with calc-alkaline to shoshonitic affinity is abundant throughout the Central-Eastern Alps. The geochemical and chronological features of the Baldissero diorites shed new light on the geodynamic evolution of the Southern Alps before the opening of the Jurassic Tethys. (C) 2010 Elsevier B.V. All rights reserved.

Zr- AND Ba-RICH MINERALS FROM THE PONTE NOVA ALKALINE MAFIC-ULTRAMAFIC MASSIF, SOUTHEASTERN BRAZIL: INDICATION OF AN ENRICHED MANTLE SOURCE

Zirconium- and Ba-rich minerals are found in gabbroic rocks from the Ponte Nova alkaline mafic-ultramafic massif in southeastern Brazil. The unusual mineralogical assemblage includes zirconolite, baddeleyite, Ba-rich alkali feldspar, and Ba- and Ti-rich biotite. Zirconolite of the Ponte Nova massif has higher levels of Zr (up to 1.172 apfu) than those registered in other terrestrial rocks and a prominent enrichment in the light rare-earth elements. Baddeleyite contains small quantities of Hf, Ti, and Fe. The Ba-rich alkali feldspar and Ba- and Ti-rich biotite contain up to 9.25 and 7.35 wt% BaO, respectively, and the biotite contains up to 12.01 wt% TiO(2). In the different intrusions of the Ponte Nova massif, such an unusual assemblage typifies the residual magma after the crystallization of clinopyroxene and olivine from previously enriched basanitic parental magma. The different trends of enrichments in REE and Th + U found for zirconolite of the intrusions of the Ponte Nova massif provide a better understanding of the variable degrees of enrichment of incompatible elements of the distinct gabbroic bodies. A lithospheric mantle source enriched in incompatible elements by the metasomatic action of volatile-rich fluids and with the presence of phlogopite or amphibole (or both) and other minor accessory phases could explain the presence of the Zr- and Ba-rich minerals in this gabbroic massif.