Quantifying volcanic ash dispersal and impact of the Campanian Ignimbrite super-eruption

1] We apply a novel computational approach to assess, for the first time, volcanic ash dispersal during the Campanian Ignimbrite (Italy) super-eruption providing insights into eruption dynamics and the impact of this gigantic event. The method uses a 3D time-dependent computational ash dispersion model, a set of wind fields, and more than 100 thickness measurements of the CI tephra deposit. Results reveal that the CI eruption dispersed 250–300 km3 of ash over ∼3.7 million km2. The injection of such a large quantity of ash (and volatiles) into the atmosphere would have caused a volcanic winter during the Heinrich Event 4, the coldest and driest climatic episode of the Last Glacial period. Fluorine-bearing leachate from the volcanic ash and acid rain would have further affected food sources and severely impacted Late Middle-Early Upper Paleolithic groups in Southern and Eastern Europe

Triboelectric charging of volcanic ash from the 2011 Grímsvötn eruption

The plume from the 2011 eruption of Grímsvötn was highly electrically charged, as shown by the considerable lightning activity measured by the United Kingdom Met Office’s low-frequency lightning detection network. Previous measurements of volcanic plumes have shown that ash particles are electrically charged up to hundreds of kilometers away from the vent, which indicates that the ash continues to charge in the plume [R. G. Harrison, K. A. Nicoll, Z. Ulanowski, and T. A. Mather, Environ. Res. Lett. 5 024004 (2010); H. Hatakeyama J. Meteorol. Soc. Jpn. 27 372 (1949)]. In this Letter, we study triboelectric charging of different size fractions of a sample of volcanic ash experimentally. Consistently with previous work, we find that the particle size distribution is a determining factor in the charging. Specifically, our laboratory experiments demonstrate that the normalized span of the particle size distribution plays an important role in the magnitude of charging generated. The influence of the normalized span on plume charging suggests that all ash plumes are likely to be charged, with implications for remote sensing and plume lifetime through scavenging effects.

Horizontal and vertical structure of the Eyjafjallajökull ash cloud over the UK: a comparison of airborne lidar observations and simulations

During April and May 2010 the ash cloud from the eruption of the Icelandic volcano Eyjafjallajökull caused widespread disruption to aviation over northern Europe. The location and impact of the eruption led to a wealth of observations of the ash cloud were being obtained which can be used to assess modelling of the long range transport of ash in the troposphere. The UK FAAM (Facility for Airborne Atmospheric Measurements) BAe-146-301 research aircraft overflew the ash cloud on a number of days during May. The aircraft carries a downward looking lidar which detected the ash layer through the backscatter of the laser light. In this study ash concentrations derived from the lidar are compared with simulations of the ash cloud made with NAME (Numerical Atmospheric-dispersion Modelling Environment), a general purpose atmospheric transport and dispersion model. The simulated ash clouds are compared to the lidar data to determine how well NAME simulates the horizontal and vertical structure of the ash clouds. Comparison between the ash concentrations derived from the lidar and those from NAME is used to define the fraction of ash emitted in the eruption that is transported over long distances compared to the total emission of tephra. In making these comparisons possible position errors in the simulated ash clouds are identified and accounted for. The ash layers seen by the lidar considered in this study were thin, with typical depths of 550–750 m. The vertical structure of the ash cloud simulated by NAME was generally consistent with the observed ash layers, although the layers in the simulated ash clouds that are identified with observed ash layers are about twice the depth of the observed layers. The structure of the simulated ash clouds were sensitive to the profile of ash emissions that was assumed. In terms of horizontal and vertical structure the best results were obtained by assuming that the emission occurred at the top of the eruption plume, consistent with the observed structure of eruption plumes. However, early in the period when the intensity of the eruption was low, assuming that the emission of ash was uniform with height gives better guidance on the horizontal and vertical structure of the ash cloud. Comparison of the lidar concentrations with those from NAME show that 2–5% of the total mass erupted by the volcano remained in the ash cloud over the United Kingdom.

Aircraft observations and model simulations of concentration and particle size distribution in the Eyjafjallajökull volcanic ash cloud

The Eyjafjallajökull volcano in Iceland emitted a cloud of ash into the atmosphere during April and May 2010. Over the UK the ash cloud was observed by the FAAM BAe-146 Atmospheric Research Aircraft which was equipped with in-situ probes measuring the concentration of volcanic ash carried by particles of varying sizes. The UK Met Office Numerical Atmospheric-dispersion Modelling Environment (NAME) has been used to simulate the evolution of the ash cloud emitted by the Eyjafjallajökull volcano during the period 4–18 May 2010. In the NAME simulations the processes controlling the evolution of the concentration and particle size distribution include sedimentation and deposition of particles, horizontal dispersion and vertical wind shear. For travel times between 24 and 72 h, a 1/t relationship describes the evolution of the concentration at the centre of the ash cloud and the particle size distribution remains fairly constant. Although NAME does not represent the effects of microphysical processes, it can capture the observed decrease in concentration with travel time in this period. This suggests that, for this eruption, microphysical processes play a small role in determining the evolution of the distal ash cloud. Quantitative comparison with observations shows that NAME can simulate the observed column-integrated mass if around 4% of the total emitted mass is assumed to be transported as far as the UK by small particles (< 30 μm diameter). NAME can also simulate the observed particle size distribution if a distal particle size distribution that contains a large fraction of < 10 μm diameter particles is used, consistent with the idea that phraetomagmatic volcanoes, such as Eyjafjallajökull, emit very fine particles.

The effect of the Madden-Julian Oscillation on station rainfall and river level in the Fly River system, Papua New Guinea

The Madden-Julian oscillation (MJO) is the dominant mode of intraseasonal variability in tropical rainfall on the large scale, but its signal is often obscured in individual station data, where effects are most directly felt at the local level. The Fly River system, Papua New Guinea, is one of the wettest regions on Earth and is at the heart of the MJO envelope. A 16 year time series of daily precipitation at 15 stations along the river system exhibits strong MJO modulation in rainfall. At each station, the difference in rainfall rate between active and suppressed MJO conditions is typically 40% of the station mean. The spread of rainfall between individual MJO events was small enough such that the rainfall distributions between wet and dry phases of the MJO were clearly separated at the catchment level. This implies that successful prediction of the large-scale MJO envelope will have a practical use for forecasting local rainfall. In the steep topography of the New Guinea Highlands, the mean and MJO signal in station precipitation is twice that in the satellite Tropical Rainfall Measuring Mission 3B42HQ product, emphasizing the need for ground-truthing satellite-based precipitation measurements. A clear MJO signal is also present in the river level, which peaks simultaneously with MJO precipitation input in its upper reaches but lags the precipitation by approximately 18 days on the flood plains.

Volcanic ash layer depth: processes and mechanisms

The long duration of the 2010 Eyjafjallajökull eruption provided a unique opportunity to measure a widely dispersed volcanic ash cloud. Layers of volcanic ash were observed by the European Aerosol Research Lidar Network with a mean depth of 1.2 km and standard deviation of 0.9 km. In this paper we evaluate the ability of the Met Office's Numerical Atmospheric-dispersion Modelling Environment (NAME) to simulate the observed ash layers and examine the processes controlling their depth. NAME simulates distal ash layer depths exceptionally well with a mean depth of 1.2 km and standard deviation of 0.7 km. The dominant process determining the depth of ash layers over Europe is the balance between the vertical wind shear (which acts to reduce the depth of the ash layers) and vertical turbulent mixing (which acts to deepen the layers). Interestingly, differential sedimentation of ash particles and the volcano vertical emission profile play relatively minor roles.

Spatial evaluation of volcanic ash forecasts using satellite observations

The decision to close airspace in the event of a volcanic eruption is based on hazard maps of predicted ash extent. These are produced using output from volcanic ash transport and dispersion (VATD)models. In this paper an objectivemetric to evaluate the spatial accuracy of VATD simulations relative to satellite retrievals of volcanic ash is presented. The 5 metric is based on the fractions skill score (FSS). Thismeasure of skill provides more information than traditional point-bypoint metrics, such as success index and Pearson correlation coefficient, as it takes into the account spatial scale overwhich skill is being assessed. The FSS determines the scale overwhich a simulation has skill and can differentiate between a "near miss" and a forecast that is badly misplaced. The 10 idealised scenarios presented show that even simulations with considerable displacement errors have useful skill when evaluated over neighbourhood scales of 200–700km2. This method could be used to compare forecasts produced by different VATDs or using different model parameters, assess the impact of assimilating satellite retrieved ash data and evaluate VATD forecasts over a long time period.

Using data insertion with the NAME model to simulate the 8 May 2010 Eyjafjallajökull volcanic ash cloud

A data insertion method, where a dispersion model is initialized from ash properties derived from a series of satellite observations, is used to model the 8 May 2010 Eyjafjallajökull volcanic ash cloud which extended from Iceland to northern Spain. We also briefly discuss the application of this method to the April 2010 phase of the Eyjafjallajökull eruption and the May 2011 Grímsvötn eruption. An advantage of this method is that very little knowledge about the eruption itself is required because some of the usual eruption source parameters are not used. The method may therefore be useful for remote volcanoes where good satellite observations of the erupted material are available, but little is known about the properties of the actual eruption. It does, however, have a number of limitations related to the quality and availability of the observations. We demonstrate that, using certain configurations, the data insertion method is able to capture the structure of a thin filament of ash extending over northern Spain that is not fully captured by other modeling methods. It also verifies well against the satellite observations according to the quantitative object-based quality metric, SAL—structure, amplitude, location, and the spatial coverage metric, Figure of Merit in Space.

Plastic and evolutionary responses to heat stress in a temperate dung fly: negative correlation between basal and induced heat tolerance?

Extreme weather events such as heat waves are becoming more frequent and intense. Populations can cope with elevated heat stress by evolving higher basal heat tolerance (evolutionary response) and/or stronger induced heat tolerance (plastic response). However, there is ongoing debate about whether basal and induced heat tolerance are negatively correlated and whether adaptive potential in heat tolerance is sufficient under ongoing climate warming. To evaluate the evolutionary potential of basal and induced heat tolerance, we performed experimental evolution on a temperate source 4 population of the dung fly Sepsis punctum. Offspring of flies adapted to three thermal selection regimes (Hot, Cold and Reference) were subjected to acute heat stress after having been exposed to either a hot-acclimation or non-acclimation pretreatment. As different traits may respond differently to temperature stress, several physiological and life history traits were assessed. Condition dependence of the response was evaluated by exposing juveniles to different levels of developmental (food restriction/rearing density) stress. Heat knockdown times were highest, whereas acclimation effects were lowest in the Hot selection regime, indicating a negative association between basal and induced heat tolerance. However, survival, adult longevity, fecundity and fertility did not show such a pattern. Acclimation had positive effects in heat-shocked flies, but in the absence of heat stress hot-acclimated flies had reduced life spans relative to nonacclimated ones, thereby revealing a potential cost of acclimation. Moreover, body size positively affected heat tolerance and unstressed individuals were less prone to heat stress than stressed flies, offering support for energetic costs associated with heat tolerance. Overall, our results indicate that heat tolerance of temperate insects can evolve under rising temperatures, but this response could be limited by a negative relationship between basal and induced thermotolerance, and may involve some but not other fitness-related traits.

How accurate are volcanic ash simulations of the 2010 Eyjafjallajökull eruption?

In the event of a volcanic eruption the decision to close airspace is based on forecast ash maps, produced using volcanic ash transport and dispersion models. In this paper we quantitatively evaluate the spatial skill of volcanic ash simulations using satellite retrievals of ash from the Eyja allajökull eruption during the period from 7 to 16 May 2010. We find that at the start of this period, 7–10 May, the model (FLEXible PARTicle) has excellent skill and can predict the spatial distribution of the satellite-retrieved ash to within 0.5∘ × 0.5∘ latitude/longitude. However, on 10 May there is a decrease in the spatial accuracy of the model to 2.5∘× 2.5∘ latitude/longitude, and between 11 and 12 May the simulated ash location errors grow rapidly. On 11 May ash is located close to a bifurcation point in the atmosphere, resulting in a rapid divergence in the modeled and satellite ash locations. In general, the model skill reduces as the residence time of ash increases. However, the error growth is not always steady. Rapid increases in error growth are linked to key points in the ash trajectories. Ensemble modeling using perturbed meteorological data would help to represent this uncertainty, and assimilation of satellite ash data would help to reduce uncertainty in volcanic ash forecasts.

Pozzolanic behavior of bamboo leaf ash: Characterization and determination of the kinetic parameters

The paper presents a characterization and study of the pozzolanic behavior between calcium hydroxide (CH) and bamboo leaf ash (BLAsh), which was obtained by calcining bamboo leaves at 600 degrees C for 2 h in a laboratory electric furnace. To evaluate the pozzolanic behavior the conductometric method was used, which is based on the measurement of the electrical conductivity in a BLAsh/CH solution with the reaction time. Later, the kinetic parameters are quantified by applying a kinetic-diffusive model. The kinetic parameters that characterize the process (in particular, the reaction rate constant and free energy of activation) were determined with relative accuracy in the fitting process of the model. The pozzolanic activity is quantitatively evaluated according to the values obtained of the kinetic parameters. Other experimental techniques, such as X-ray diffraction (XRD) and scanning electron microscopy (SEM), were also employed. The results show that this kind of ash is formed by silica with a completely amorphous nature and a high pozzolanic activity. The correlation between the values of free energy of activation (Delta G(#)) and the reaction rate constants (K) are in correspondence with the theoretical studies about the rate processes reported in the literature. (C) 2010 Elsevier Ltd. All rights reserved.

A study of the chemical and physical properties of cashew nut shell ash for use in cement materials.

A study of the chemical and physical properties of cashew nut shell ash for use in cement materials. Ash occupies a prominent place among agro-industrial wastes, as it is derived from energy generation processes. Several types of ash have pozzolanic reactivity, and might be used as replacement material for cement, resulting in less energy waste and lower cost. This work aimed to investigate the physical and chemical properties of the cashew nut shell ash (CNSA), by performing the following measurement tests: chemical analysis, bulk density, specific mass, leaching and solubilization process, X-ray diffraction (XrD), specific surface area (BET) and pozzolanicity analysis with cement and lime. The results indicate a low reactivity of CNSA and the presence of heavy metals, alkalis and phenol.

Knockdown resistance in pyrethroid-resistant horn fly (Diptera: Muscidae) populations in Brazil

To investigate the kdr (knockdown resistance) resistance-associated gene mutation and determine its frequency in pyrethroid-resistant horn fly (Haematobia irritans) populations, a total of 1,804 horn flies of 37 different populations from all Brazilian regions (North, Northeast, Central-West, Southeast, and South) were molecular screened through polymerase chain reaction (PCR). The kdr gene was not detected in 87.08% of the flies. However, the gene was amplified in 12.92% of the flies, of which 11.70% were resistant heterozygous and 1.22% were resistant homozygous. Deviation from Hardy-Weinberg equilibrium (HWE) was found only in 1 ranch with an excess of heterozygous. When populations were grouped by region, three metapopulations showed significant deviations of HWE (Central-West population, South population and Southeast population). This indicates that populations are isolated one from another and kdr occurrence seems to be an independent effect probably reflecting the insecticide strategy used by each ranch. Although resistance to pyrethroids is disseminated throughout Brazil, only 48% of resistant populations had kdr flies, and the frequency of kdr individuals in each of these resistant populations was quite low. But this study shows that, with the apparent exception of the Northeast region, the kdr mechanism associated with pyrethroid resistance occurs all over Brazil.

Spatial and temporal patterns of occurrence of Lutzomyia sand fly species in an endemic area for cutaneous leishmaniasis in the Atlantic Forest region of northeast Brazil

Sand fly populations of different ecological niches in the Amaraji endemic American Cutaneous Leishmaniasis (ACL) focus of the Pernambuco Atlantic Forest region of northeastern Brazil were monitored spatiotemporally. Lutzomyia whitmani was dominant in all niches but occurred in smaller numbers in forested locations. L. whitmani was significantly less seasonal than the other species, being present throughout the year while other species were more abundant between February and April. These results suggest that L. whitmani may potentially be the principal vector of ACL in the region, even though the sand fly fauna was diverse: 88% were L.whitmani and 12% belonged to 11 other species. Two other species, L. complexa (1.3%) and L. migonei (0.8%), considered to be ACL vectors in other regions, were also present. This detailed picture of the sand fly population`s abundance and spatiotemporal distribution provides a basis for future modeling studies of forecasting sand fly activity patterns and ACL occurence.

30 million years of Permian volcanism recorded in the Choiyoi igneous province (W Argentina) and their source for younger ash fall deposits in the Parana Basin: SHRIMP U-Pb zircon geochronology evidence

We present four SHRIMP U-Pb zircon ages for the Choiyoi igneous province from the San Rafael Block, central-western Argentina. Dated samples come from the Yacimiento Los Reyunos Formation (281.4 +/- 2.5 Ma) of the Cochico Group (Lower Choiyoi section: andesitic breccias, dacitic to rhyolitic ignimbrites and continental conglomerates). Agua de los Burros Formation (264.8 +/- 2.3 Ma and 264.5 +/- 3.0 Ma) and Cerro Carrizalito Formation (251.9 +/- 2.7 Ma Upper Choiyoi section: rhyolitic ignimbrites and pyroclastic flows) spanning the entire Permian succession of the Choiyoi igneous province. A single ziron from the El Imperial Formation, that is overlain unconformably by the Choiyoi succession, yielded an early Permian age (297.2 +/- 5.3 Ma). while the main detrital zircon population indicated an Ordovician age (453.7 +/- 8.1 Ma). The new data establishes a more precise Permian age (Artinskian-Lopingian) for the section studied spanning 30 Ma of volcanic activity. Volcanological observations for the Choiyoi succession support the occurrence of explosive eruptions of plinian to ultraplinian magnitudes, capable of injecting enormous volumes of tephra in the troposphere-stratosphere. The new SHRIMP ages indicate contemporaneity between the Choyoi succession and the upper part of the Parana Basin late Paleozoic section, from the Irad up to the Rio do Rasto formations, encompassing about 24 Ma. Geochemical data show a general congruence in compositional and tectonic settings between the volcanics and Parana Basin Permian ash fall derived layers of bentonites. Thickness and granulometry of ash fall layers broadly fit into the depletion curve versus distance from the remote source vent of ultraplinian eruptions. Thus, we consider that the Choiyoi igneous province was the source of ash fall deposits in the upper Permian section of the Parana Basin. Data presented here allow a more consistent correlation between tectono-volcanic Permian events along the paleo-Pacific margin of southwestern Gondwana and the geological evolution of neighboring Paleozoic foreland basins in South America and Africa. (C) 2010 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.