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.

Role of coronal mass ejections in the heliospheric Hale cycle

The 11-year solar cycle variation in the heliospheric magnetic field strength can be explained by the temporary buildup of closed flux released by coronal mass ejections (CMEs). If this explanation is correct, and the total open magnetic flux is conserved, then the interplanetary-CME closed flux must eventually open via reconnection with open flux close to the Sun. In this case each CME will move the reconnected open flux by at least the CME footpoint separation distance. Since the polarity of CME footpoints tends to follow a pattern similar to the Hale cycle of sunspot polarity, repeated CME eruption and subsequent reconnection will naturally result in latitudinal transport of open solar flux. We demonstrate how this process can reverse the coronal and heliospheric fields, and we calculate that the amount of flux involved is sufficient to accomplish the reversal within the 11 years of the solar cycle.

Variation of the mixing state of Saharan dust particles with atmospheric transport

Mineral dust is an important aerosol species in the Earth’s atmosphere and has a major source within North Africa, of which the Sahara forms the major part. Aerosol Time of Flight Mass Spectrometry (ATOFMS) is first used to determine the mixing state of dust particles collected from the land surface in the Saharan region, showing low abundance of species such as nitrate and sulphate internally mixed with the dust mineral matrix. These data are then compared with the ATOFMS single particle mass spectra of Saharan dust particles detected in the marine atmosphere in the vicinity of the Cape Verde islands, which are further compared with those from particles with longer atmospheric residence sampled at a coastal station at Mace Head, Ireland. Saharan dust particles collected near the Cape Verde Islands showed increased internally mixed nitrate but no sulphate, whilst Saharan dust particles collected on the coast of Ireland showed a very high degree of internally mixed secondary species including nitrate, sulphate and methanesulphonate. This uptake of secondary species will change the pH and hygroscopic properties of the aerosol dust and thus can influence the budgets of other reactive gases, as well as influencing the radiative properties of the particles and the availability of metals for dissolution.

Influence of gut morphology on passive transport of freshwater bryozoans by waterfowl in Donana (southwestern Spain)

Waterbirds have been proposed as important vectors for the passive dispersal of those aquatic invertebrates and plants that lack a capacity for active dispersal between isolated water bodies. We analysed the frequency of internal transport of bryozoan propagules (statoblasts) by waterbirds in Donana, Spain, by examining their presence in the intestines and ceca of dead birds and analysing the role of different aspects of gut characteristics in explaining variation in the presence/absence and abundance of statoblasts. Of the 228 samples examined, 7.9% presented intact statoblasts of Plumatella fungosa (Pallas, 1768), Plumatella emarginata Allman, 1844, and two unidentified Plumatella species. For a given bird species, individuals with heavier gizzards and shorter ceca had a lower incidence and abundance of statoblasts in the lower gut. Grit mass and intestine length were unrelated to the presence or abundance of statoblasts. Our results suggest that waterbirds frequently transport bryozoans on a local scale, with lighter gizzards and longer ceca favouring such transport. Lighter gizzards are likely to destroy fewer propagules before they reach the lower gut. Species and individuals with longer ceca are particularly good candidates for long-distance dispersal of bryozoans, given the longer passage time of propagules that enter the ceca.

Estimating relationships between air mass origin and chemical composition

Observations of a chemical at a point in the atmosphere typically show sudden transitions between episodes of high and low concentration. Often these are associated with a rapid change in the origin of air arriving at the site. Lagrangian chemical models riding along trajectories can reproduce such transitions, but small timing errors from trajectory phase errors dramatically reduce the correlation between modeled concentrations and observations. Here the origin averaging technique is introduced to obtain maps of average concentration as a function of air mass origin for the East Atlantic Summer Experiment 1996 (EASE96, a ground-based chemistry campaign). These maps are used to construct origin averaged time series which enable comparison between a chemistry model and observations with phase errors factored out. The amount of the observed signal explained by trajectory changes can be quantified, as can the systematic model errors as a function of air mass origin. The Cambridge Tropospheric Trajectory model of Chemistry and Transport (CiTTyCAT) can account for over 70% of the observed ozone signal variance during EASE96 when phase errors are side-stepped by origin averaging. The dramatic increase in correlation (from 23% without averaging) cannot be achieved by time averaging. The success of the model is attributed to the strong relationship between changes in ozone along trajectories and their origin and its ability to simulate those changes. The model performs less well for longer-lived chemical constituents because the initial conditions 5 days before arrival are insufficiently well known.

A comparison of stratosphere-troposphere transport in convection-permitting and convection-parameterizing simulations of three mesoscale convective systems

The transport of stratospheric air into the troposphere within deep convection was investigated using the Met Office Unified Model version 6.1. Three cases were simulated in which convective systems formed over the UK in the summer of 2005. For each of these three cases, simulations were performed on a grid having 4 km horizontal grid spacing in which the convection was parameterized and on a grid having 1 km horizontal grid spacing, which permitted explicit representation of the largest energy-containing scales of deep convection. Cross-tropopause transport was diagnosed using passive tracers that were initialized above the dynamically defined tropopause (2 potential vorticity unit surface) with a mixing ratio of 1. Although the synoptic-scale environment and triggering mechanisms varied between the cases, the total simulated transport was similar in all three cases. The total stratosphere-to-troposphere transport over the lifetime of the convective systems ranged from 25 to 100 kg/m2 across the simulated convective systems and resolutions, which corresponds to ∼5–20% of the total mass located within a stratospheric column extending 2 km above the tropopause. In all simulations, the transport into the lower troposphere (defined as below 3.5 km elevation) accounted for ∼1% of the total transport across the tropopause. In the 4 km runs most of the transport was due to parameterized convection, whereas in the 1 km runs the transport was due to explicitly resolved convection. The largest difference between the simulations with different resolutions occurred in the one case of midlevel convection considered, in which the total transport in the 1 km grid spacing simulation with explicit convection was 4 times that in the 4 km grid spacing simulation with parameterized convection. Although the total cross-tropopause transport was similar, stratospheric tracer was deposited more deeply to near-surface elevations in the convection-parameterizing simulations than in convection-permitting simulations.

Ultra-distal tephra deposits from super-eruptions: examples from Toba, Indonesia and Taupo Volcanic Zone, New Zealand

Voluminous rhyolitic eruptions from Toba, Indonesia, and Taupo Volcanic Zone (TVZ), New Zealand, have dispersed volcanic ash over vast areas in the late Quaternary. The ~74 ka Youngest Toba Tuff (YTT) eruption deposited ash over the Bay of Bengal and the Indian subcontinent to the west. The ~340 ka Whakamaru eruption (TVZ) deposited the widespread Rangitawa Tephra, dominantly to the southeast (in addition to occurrences northwest of vent), extending across the landmass of New Zealand, and the South Pacific Ocean and Tasman Sea, with distal terrestrial exposures on the Chatham Islands. These super-eruptions involved ~2500 km^3 and ~1500 km3 of magma (dense-rock equivalent; DRE), respectively. Ultra-distal terrestrial exposures of YTT at two localities in India, Middle Son Valley, Madhya Pradesh, and Jurreru River Valley, Andhra Pradesh, at distances of >2000 km from the source caldera, show a basal ‘primary’ ashfall unit ~4 cm thick, although deposits containing reworked ash are up to ~3 m in total thickness. Exposures of Rangitawa Tephra on the Chatham Islands, >900 km from the source caldera, are ~15-30 cm thick. At more proximal localities (~200 km from source), Rangitawa Tephra is ~55-70 cm thick and characterized by a crystal-rich basal layer and normal grading. Both distal tephra deposits are characterized by very-fine ash (with high PM10 fractions) and are crystal-poor. Glass chemistry, stratigraphy and grain-size data for these distal tephra deposits are presented with comparisons of their correlation, dispersal and preservation. Using field observations, ash transport and deposition were modeled for both eruptions using a semi-analytical model (HAZMAP), with assumptions concerning average wind direction and strength during eruption, column shape and vent size. Model outputs provide new insights into eruption dynamics and better estimates of eruption volumes associ- ated with tephra fallout. Modeling based on observed YTT distal tephra thicknesses indicate a relatively low (<40 km high), very turbulent eruption column, consistent with deposition from a co-ignimbrite cloud extending over a broad region. Similarly, the Whakamaru eruption was modeled as producing a predominantly Plinian column (~45 km high), with dispersal to the southeast by strong prevailing winds. Significant ash fallout of the main dispersal direction, to the northwest of source, cannot be replicated in this modeling. The widespread dispersal of large volumes of fine ash from both eruptions may have had global environmental consequences, acutely affecting areas up to thousands of kilometers from vent.

Enhanced viability of corneal epithelial cells for efficient transport/storage using a structurally-modified calcium alginate hydrogel

Aims: Therapeutic limbal epithelial stem cells could be managed more efficiently if clinically validated batches were transported for ‘on-demand’ use. Materials & methods: In this study, corneal epithelial cell viability in calcium alginate hydrogels was examined under cell culture, ambient and chilled conditions for up to 7 days. Results: Cell viability improved as gel internal pore size increased, and was further enhanced with modification of the gel from a mass to a thin disc. Ambient storage conditions were optimal for supporting cell viability in gel discs. Cell viability in gel discs was significantly enhanced with increases in pore size mediated by hydroxyethyl cellulose. Conclusion: Our novel methodology of controlling alginate gel shape and pore size together provides a more practical and economical alternative to established corneal tissue/cell storage methods.

Estimating the volume of tephra deposits: a new simple strategy

Volume determination of tephra deposits is necessary for the assessment of the dynamics and hazards of explosive volcanoes. Several methods have been proposed during the past 40 years that include the analysis of crystal concentration of large pumices, integrations of various thinning relationships, and the inversion of field observations using analytical and computational models. Regardless of their strong dependence on tephra-deposit exposure and distribution of isomass/isopach contours, empirical integrations of deposit thinning trends still represent the most widely adopted strategy due to their practical and fast application. The most recent methods involve the best fitting of thinning data using various exponential seg- ments or a power-law curve on semilog plots of thickness (or mass/area) versus square root of isopach area. The exponential method is mainly sensitive to the number and the choice of straight segments, whereas the power-law method can better reproduce the natural thinning of tephra deposits but is strongly sensitive to the proximal or distal extreme of integration. We analyze a large data set of tephra deposits and propose a new empirical method for the deter- mination of tephra-deposit volumes that is based on the integration of the Weibull function. The new method shows a better agreement with observed data, reconciling the debate on the use of the exponential versus power-law method. In fact, the Weibull best fitting only depends on three free parameters, can well reproduce the gradual thinning of tephra deposits, and does not depend on the choice of arbitrary segments or of arbitrary extremes of integration.

Numerical simulation of tephra transport and deposition of the 1982 El Chichón eruption and implications for hazard assessment

El Chichón volcano, Chiapas, Mexico, erupted explosively on March 29th, 1982, after a repose period of about 550 years. Amongst ten eruptive episodes documented between March 29th and April 4th, only the three that occurred on March 29th and April 4th produced significant pyroclastic tephra deposits. Here we use analytical (HAZMAP) and numerical (FALL3D) tephra transport models to reconstruct the deposits and the atmospheric plume dispersal associated with the three main fallout units of the 1982 eruption. On the basis of such a reconstruction, we produce hazard maps of tephra fallout associated to a Plinian eruption and discuss the implications of such a severe eruption scenario.

Desert dust deposition on Mt. Elbrus, Caucasus Mountains, Russia in 2009-2012 as recorded in snow and shallow ice core: high-resolution "provenancing", transport patterns, physical properties and soluble ionic composition

A record of dust deposition events between 2009 and 2012 on Mt. Elbrus, Caucasus Mountains derived from a snow pit and a shallow ice core is presented for the first time for this region. A combination of isotopic analysis, SEVIRI red-green-blue composite imagery, MODIS atmospheric optical depth fields derived using the Deep Blue algorithm, air mass trajectories derived using the HYSPLIT model and analysis of meteorological data enabled identification of dust source regions with high temporal (hours) and spatial (cf. 20–100 km) resolution. Seventeen dust deposition events were detected; fourteen occurred in March–June, one in February and two in October. Four events originated in the Sahara, predominantly in north-eastern Libya and eastern Algeria. Thirteen events originated in the Middle East, in the Syrian Desert and northern Mesopotamia, from a mixture of natural and anthropogenic sources. Dust transportation from Sahara was associated with vigorous Saharan depressions, strong surface winds in the source region and mid-tropospheric south-westerly flow with daily winds speeds of 20–30 m s−1 at 700 hPa level and, although these events were less frequent, they resulted in higher dust concentrations in snow. Dust transportation from the Middle East was associated with weaker depressions forming over the source region, high pressure centered over or extending towards the Caspian Sea and a weaker southerly or south-easterly flow towards the Caucasus Mountains with daily wind speeds of 12–18 m s−1 at 700 hPa level. Higher concentrations of nitrates and ammonium characterise dust from the Middle East deposited on Mt. Elbrus in 2009 indicating contribution of anthropogenic sources. The modal values of particle size distributions ranged between 1.98 μm and 4.16 μm. Most samples were characterised by modal values of 2.0–2.8 μm with an average of 2.6 μm and there was no significant difference between dust from the Sahara and the Middle East.

Quantification of chemical and physical processes influencing ozone during long-range transport using a trajectory ensemble

During long-range transport, many distinct processes – including photochemistry, deposition, emissions and mixing – contribute to the transformation of air mass composition. Partitioning the effects of different processes can be useful when considering the sensitivity of chemical transformation to, for example, a changing environment or anthropogenic influence. However, transformation is not observed directly, since mixing ratios are measured, and models must be used to relate changes to processes. Here, four cases from the ITCT-Lagrangian 2004 experiment are studied. In each case, aircraft intercepted a distinct air mass several times during transport over the North Atlantic, providing a unique dataset and quantifying the net changes in composition from all processes. A new framework is presented to deconstruct the change in O3 mixing ratio (Δ O3) into its component processes, which were not measured directly, taking into account the uncertainty in measurements, initial air mass variability and its time evolution. The results show that the net chemical processing (Δ O3chem) over the whole simulation is greater than net physical processing (Δ O3phys) in all cases. This is in part explained by cancellation effects associated with mixing. In contrast, each case is in a regime of either net photochemical destruction (lower tropospheric transport) or production (an upper tropospheric biomass burning case). However, physical processes influence O3 indirectly through addition or removal of precursor gases, so that changes to physical parameters in a model can have a larger effect on Δ O3chem than Δ O3phys. Despite its smaller magnitude, the physical processing distinguishes the lower tropospheric export cases, since the net photochemical O3 change is −5 ppbv per day in all three cases. Processing is quantified using a Lagrangian photochemical model with a novel method for simulating mixing through an ensemble of trajectories and a background profile that evolves with them. The model is able to simulate the magnitude and variability of the observations (of O3, CO, NOy and some hydrocarbons) and is consistent with the time-average OH following air-masses inferred from hydrocarbon measurements alone (by Arnold et al., 2007). Therefore, it is a useful new method to simulate air mass evolution and variability, and its sensitivity to process parameters.

Tracing troposphere-to-stratosphere transport above a mid-latitude deep convective system

Within the project SPURT (trace gas measurements in the tropopause region) a variety of trace gases have been measured in situ in order to investigate the role of dynamical and chemical processes in the extra-tropical tropopause region. In this paper we report on a flight on 10 November 2001 leading from Hohn, Germany (52�N) to Faro, Portugal (37�N) through a strongly developed deep stratospheric intrusion. This streamer was associated with a large convective system over the western Mediterranean with potentially significant troposphere-to-stratosphere transport. Along major parts of the flight we measured unexpectedly high NOy mixing ratios. Also H2O mixing ratios were significantly higher than stratospheric background levels confirming the extraordinary chemical signature of the probed air masses in the interior of the streamer. Backward trajectories encompassing the streamer enable to analyze the origin and physical characteristics of the air masses and to trace troposphere-to-stratosphere transport. Near the western flank of the streamer features caused by long range transport, such as tropospheric filaments characterized by sudden drops in the O3 and NOy mixing ratios and enhanced CO and H2O can be reconstructed in great detail using the reverse domain filling technique. These filaments indicate a high potential for subsequent mixing with the stratospheric air. At the south-western edge of the streamer a strong gradient in the NOy and the O3 mixing ratios coincides very well with a sharp gradient in potential vorticity in the ECMWF fields. In contrast, in the interior of the streamer the observed highly elevated NOy and H2O mixing ratios up to a potential temperature level of 365K and potential vorticity values of maximum 10 PVU cannot be explained in terms of resolved troposphere-to-stratosphere transport along the backward trajectories. Also mesoscale simulations with a High Resolution Model reveal no direct evidence for convective H2O injection up to this level. Elevated H2O mixing ratios in the ECMWF and HRM are seen only up to about tropopause height at 340 hPa and 270 hPa, respectively, well below flight altitude of about 200 hPa. However, forward tracing of the convective influence as identified by satellite brightness temperature measurements and counts of lightning strokes shows that during this part of the flight the aircraft was closely following the border of an air mass which was heavily impacted by convective activity over Spain and Algeria. This is evidence that deep convection at mid-latitudes may have a large impact on the tracer distribution of the lowermost stratosphere reaching well above the thunderstorms anvils as claimed by recent studies using cloud-resolving models.

General circulation model estimates of aerosol transport and radiative forcing during the Indian Ocean Experiment

Aerosol sources, transport, and sinks are simulated, and aerosol direct radiative effects are assessed over the Indian Ocean for the Indian Ocean Experiment (INDOEX) Intensive Field Phase during January to March 1999 using the Laboratoire de Me´te´orologie Dynamique (LMDZT) general circulation model. The model reproduces the latitudinal gradient in aerosol mass concentration and optical depth (AOD). The model-predicted aerosol concentrations and AODs agree reasonably well with measurements but are systematically underestimated during high-pollution episodes, especially in the month of March. The largest aerosol loads are found over southwestern China, the Bay of Bengal, and the Indian subcontinent. Aerosol emissions from the Indian subcontinent are transported into the Indian Ocean through either the west coast or the east coast of India. Over the INDOEX region, carbonaceous aerosols are the largest contributor to the estimated AOD, followed by sulfate, dust, sea salt, and fly ash. During the northeast winter monsoon, natural and anthropogenic aerosols reduce the solar flux reaching the surface by 25 W m�2, leading to 10–15% less insolation at the surface. A doubling of black carbon (BC) emissions from Asia results in an aerosol single-scattering albedo that is much smaller than in situ measurements, reflecting the fact that BC emissions are not underestimated in proportion to other (mostly scattering) aerosol types. South Asia is the dominant contributor to sulfate aerosols over the INDOEX region and accounts for 60–70% of the AOD by sulfate. It is also an important but not the dominant contributor to carbonaceous aerosols over the INDOEX region with a contribution of less than 40% to the AOD by this aerosol species. The presence of elevated plumes brings significant quantities of aerosols to the Indian Ocean that are generated over Africa and Southeast and east Asia.

Sources, sinks, and transatlantic transport of North African dust aerosol: a multimodel analysis and comparison with remote sensing data

This study evaluates model-simulated dust aerosols over North Africa and the North Atlantic from five global models that participated in the Aerosol Comparison between Observations and Models phase II model experiments. The model results are compared with satellite aerosol optical depth (AOD) data from Moderate Resolution Imaging Spectroradiometer (MODIS), Multiangle Imaging Spectroradiometer (MISR), and Sea-viewing Wide Field-of-view Sensor, dust optical depth (DOD) derived from MODIS and MISR, AOD and coarse-mode AOD (as a proxy of DOD) from ground-based Aerosol Robotic Network Sun photometer measurements, and dust vertical distributions/centroid height from Cloud Aerosol Lidar with Orthogonal Polarization and Atmospheric Infrared Sounder satellite AOD retrievals. We examine the following quantities of AOD and DOD: (1) the magnitudes over land and over ocean in our study domain, (2) the longitudinal gradient from the dust source region over North Africa to the western North Atlantic, (3) seasonal variations at different locations, and (4) the dust vertical profile shape and the AOD centroid height (altitude above or below which half of the AOD is located). The different satellite data show consistent features in most of these aspects; however, the models display large diversity in all of them, with significant differences among the models and between models and observations. By examining dust emission, removal, and mass extinction efficiency in the five models, we also find remarkable differences among the models that all contribute to the discrepancies of model-simulated dust amount and distribution. This study highlights the challenges in simulating the dust physical and optical processes, even in the best known dust environment, and stresses the need for observable quantities to constrain the model processes.