Numerical and granulometric approaches to geophysical granular flows

Volcanic ash clouds can be fed by an upward-directed eruption column (Plinian column) or by elutriation from extensive pyroclastic-flows (coignimbrite cloud). For large-scale eruptions, there is considerable uncertainty about which mechanism is dominant. Here we analyze in a novel way a comprehensive grainsize database for pyroclastic deposits. We demonstrate that the Mount Pinatubo climactic eruption deposits were substantially derived from coignimbrite clouds, and not only by a Plinian cloud as generally thought. Coignimbrite ash-fall deposits are much richer in breathable <10 m ash (5–25 wt%) than pure Plinian ash at most distances from the source volcano. We also show that coignimbrite ash clouds, as at Pinatubo, are expected to be more water rich than Plinian clouds, leading to removal of more HCl prior to stratospheric injection, thereby reducing their atmospheric impact.

GROUND-BASED AND SATELLITE REMOTE SENSING OF PAROXYSMAL ERUPTIONS AT ETNA VOLCANO, 2011-2012.

Mt Etna's activity has increased during the last decade with a tendency towards more explosive eruptions that produce paroxysmal lava fountains. From January 2011 to April 2012, 25 lava fountaining episodes took place at Etna's New South-East Crater (NSEC). Improved understanding of the mechanism driving these explosive basaltic eruptions is needed to reduce volcanic hazards. This type of activity produces high sulfur dioxide (SO2) emissions, associated with lava flows and ash fall-out, but to date the SO2 emissions associated with Etna's lava fountains have been poorly constrained. The Ultraviolet (UV) Ozone Monitoring Instrument (OMI) on NASA's Aura satellite and the Atmospheric Infrared Sounder (AIRS) on Aqua were used to measure the SO2 loadings. Ground-based data from the Observatoire de Physique du Globe de Clermont-Ferrand (OPGC) L-band Doppler radar, VOLDORAD 2B, used in collaboration with the Italian National Institute of Geophysics and Volcanology in Catania (INGV-CT), also detected the associated ash plumes, giving precise timing and duration for the lava fountains. This study resulted in the first detailed analysis of the OMI and AIRS SO2 data for Etna's lava fountains during the 2011-2012 eruptive cycle. The HYSPLIT trajectory model is used to constrain the altitude of the observed SO2 clouds, and results show that the SO2 emission usually coincided with the lava fountain peak intensity as detected by VOLDORAD. The UV OMI and IR AIRS SO2 retrievals permit quantification of the SO2 loss rate in the volcanic SO2 clouds, many of which were tracked for several days after emission. A first attempt to quantitatively validate AIRS SO2 retrievals with OMI data revealed a good correlation for high altitude SO2 clouds. Using estimates of the emitted SO2 at the time each paroxysm, we observe a correlation with the inter-paroxysm repose time. We therefore suggest that our data set supports the collapsing foam (CF) model [1] as driving mechanism for the paroxysmal events at the NSEC. Using VOLDORAD-based estimates of the erupted magma mass, we observe a large excess of SO2 in the eruption clouds. Satellite measurements indicate that SO2 emissions from Etnean lava fountains can reach the lower stratosphere and hence could pose a hazard to aviation. [1] Parfitt E.A (2004). A discussion of the mechanisms of explosive basaltic eruptions. J. Volcanol. Geotherm. Res. 134, 77-107.

CHARACTERIZING THE FIRST HISTORIC ERUPTION OF NABRO, ERITREA: INSIGHTS FROM THERMAL AND UV REMOTE SENSING

June 2011 saw the first historic eruption of Nabro volcano, one of an ongoing sequence of eruptions in the Afar-Red Sea region since 2005. It halted air travel in northern Africa, contaminated food and water sources, and displaced thousands from their homes. Due to its remote location, little was known about this event in terms of the quantity of erupted products and the timing and mechanisms of their emplacement. Geographic isolation, previous quiescence and regional civil unrest meant that this volcano was effectively unmonitored at the time of eruption, and opportunities for field study are limited. Using free, publicly available satellite data, I examined rates of lava effusion and SO2 emission in order to quantify the amount of erupted products and understand the temporal evolution of the eruption, as well as explore what information can be gleaned about eruption mechanisms using remote sensing data. These data revealed a bimodal eruption, beginning with explosive activity marked by high SO2 emission totalling 1824 - 2299 KT, and extensive ash fall of 270 - 440 km2. This gave way to a period of rapid effusion, producing a ~17 km long lava flow, and a volume of ~22.1 x 106 m3. Mass balance between the SO2 and lava flows reveals no sulfur 'excess', suggesting that nearly all of the degassed magma was extruded. The 2011 eruption of Nabro lasted nearly 6 weeks, and may be considered the second largest historic eruption in Africa. Work such as this highlights the importance of satellite remote sensing for studying and monitoring volcanoes, particularly those in remote regions that may be otherwise inaccessible.

Invasion patterns of emerald ash borer and European earthworms in forested ecosystems

Invasive and exotic species present a serious threat to the health and sustainability of natural ecosystems. These species often benefit from anthropogenic activities that aid their introduction and dispersal. This dissertation focuses on invasion dynamics of the emerald ash borer, native to Asia, and European earthworms. These species have shown detrimental impacts in invaded forest ecosystems across the Great Lakes region, and continue to spread via human-assisted long distance dispersal and by natural modes of dispersal into interior forests from areas of introduction. Successful forest management requires that the impact and effect of invasive species be considered and incorporated into management plans. Understanding patterns and constraints of introduction, establishment, and spread will aid in this effort. To assist in efforts to locate introduction points of emerald ash borer, a multicriteria risk model was developed to predict the highest risk areas. Important parameters in the model were road proximity, land cover type, and campground proximity. The model correctly predicted 85% of known emerald ash borer invasion sites to be at high risk. The model’s predictions across northern Michigan can be used to focus and guide future monitoring efforts. Similar modeling efforts were applied to the prediction of European earthworm invasion in northern Michigan forests. Field sampling provided a means to improve upon modeling efforts for earthworms to create current and future predictions of earthworm invasion. Those sites with high soil pH and high basal area of earthworm preferred overstory species (such as basswood and maples) had the highest likelihood of European earthworm invasion. Expanding beyond Michigan into the Upper Great Lakes region, earthworm populations were sampled across six National Wildlife Refuges to identify potential correlates and deduce specific drivers and constraints of earthworm invasion. Earthworm communities across all refuges were influenced by patterns of anthropogenic activity both within refuges and in surrounding ecoregions of study. Forest composition, soil pH, soil organic matter, anthropogenic cover, and agriculture proximity also proved to be important drivers of earthworm abundance and community composition. While there are few management options to remove either emerald ash borer or European earthworms from forests after they have become well established, prevention and early detection are important and can be beneficial. An improved understanding the factors controlling the distribution and invasion patterns of exotic species across the landscape will aid efforts to determine their consequences and generate appropriate forest management solutions to sustain ecosystem health in the presence of these invaders.

Detection and landing behavior of emerald ash borer, Agrilus planipennis, at low population density

The exotic emerald ash borer, Agrilus planipennis Fairmaire (Coleoptera: Buprestidae), was first discovered in North America in southeastern Michigan, USA, and Windsor, Ontario, Canada in 2002. Significant ash (Fraxinus spp.) mortality has been caused in areas where this insect has become well established, and new infestations continue to be discovered in several states in the United States and in Canada. This beetle is difficult to detect when it invades new areas or occurs at low density. Girdled trap tree and ground surveys have been important tools for detecting emerald ash borer populations, and more recently, purple baited prism traps have been used in detection efforts. Girdled trap trees were found to be more effective than purple prism traps at detecting emerald ash borer as they acted as sinks for larvae in an area of known low density emerald ash borer infestation. The canopy condition of the trap trees was not predictive of whether they were infested or not, indicating that ground surveys may not be effective for detection in an area of low density emerald ash borer population. When landing rates of low density emerald ash borer populations were monitored on non-girdled ash trees, landing rates were higher on larger, open grown trees with canopies that contain a few dead branches. As a result of these studies, we suggest that the threshold for emerald ash borer detection using baited purple prism traps hung at the canopy base of trees is higher than for girdled trap trees. In addition, detection of developing populations of EAB may be possible by selectively placing sticky trapping surfaces on non-girdled trap trees that are the larger and more open grown trees at a site.

Finite element study on fracture patterns in human scaphoid wrist bone due to free fall

Scaphoid is one of the 8 carpal bones found adjacent to the thumb supported proximally by Radius bone. During the free fall, on outstretched hand, the impact load gets transferred to the scaphoid at its free anterior end. Unique arrangement of other carpal bones in the palm is also one of the reasons for the load to get transferred to scaphoid. About half of the total load acting upon carpal bone gets transferred to scaphoid at its distal pole. There are about 10 to 12 clinically observed fracture pattern in the scaphoid due to free fall. The aim of the study is to determine the orientation of the load, magnitude of the load and the corresponding fracture pattern. This study includes both static and dynamic finite element models validated by experiments. The scaphoid model has been prepared from CT scans of a 27 year old person. The 2D slices of the CT scans have been converted to 3D model by using MIMICS software. There are four cases of loading studied which are considered to occur clinically more frequently. In case (i) the load is applied at the posterior end at distal pole whereas in case (ii), (iii) and (iv), the load is applied at anterior end at different directions. The model is given a fixed boundary condition at the region which is supported by Radius bone during the impact. Same loading and boundary conditions have been used in both static and dynamic explicit finite element analysis. The site of fracture initiation and path of fracture propagation have been identified by using max principal stress / gradient and max principal strain / gradient criterion respectively in static and dynamic explicit finite element analysis. Static and dynamic impact experiments were performed on the polyurethane foam specimens to validate the finite element results. Experimental results such as load at fracture, site of fracture initiation and path of fracture propagation have been compared with the results of finite element analysis. Four different types of fracture patterns observed in clinical studies have been identified in this study.

Characterization of a Coal Fly Ash-Cement Slurry by the Absolute Foam Index

This dissertation established a standard foam index: the absolute foam index test. This test characterized a wide range of coal fly ash by the absolute volume of air-entraining admixture (AEA) necessary to produce a 15-second metastable foam in a coal fly ash-cement slurry in a specified time. The absolute foam index test was used to characterize fly ash samples having loss on ignition (LOI) values that ranged from 0.17 to 23.3 %wt. The absolute foam index characterized the fly ash samples by absolute volume of AEA, defined as the amount of undiluted AEA solution added to obtain a 15-minute endpoint signified by 15-second metastable foam. Results were compared from several foam index test time trials that used different initial test concentrations to reach termination at selected times. Based on the coefficient of variation (CV), a 15-minute endpoint, with limits of 12 to 18 minutes was chosen. Various initial test concentrations were used to accomplish consistent contact times and concentration gradients for the 15-minute test endpoint for the fly ash samples. A set of four standard concentrations for the absolute foam index test were defined by regression analyses and a procedure simplifying the test process. The set of standard concentrations for the absolute foam index test was determined by analyzing experimental results of 80 tests on coal fly ashes with loss on ignition (LOI) values ranging from 0.39 to 23.3 wt.%. A regression analysis informed selection of four concentrations (2, 6, 10, and 15 vol.% AEA) that are expected to accommodate fly ashes with 0.39 to 23.3 wt.% LOI, depending on the AEA type. Higher concentrations should be used for high-LOI fly ash when necessary. A procedure developed using these standard concentrations is expected to require only 1-3 trials to meet specified endpoint criteria for most fly ashes. The AEA solution concentration that achieved the metastable foam in the foam index test was compared to the AEA equilibrium concentration obtained from the direct adsorption isotherm test with the same fly ash. The results showed that the AEA concentration that satisfied the absolute foam index test was much less than the equilibrium concentration. This indicated that the absolute foam index test was not at or near equilibrium. Rather, it was a dynamic test where the time of the test played an important role in the results. Even though the absolute foam index was not an equilibrium condition, a correlation was made between the absolute foam index and adsorption isotherms. Equilibrium isotherm equations obtained from direct isotherm tests were used to calculate the equilibrium concentrations and capacities of fly ash from 0.17 to 10.5% LOI. The results showed that the calculated fly ash capacity was much less than capacities obtained from isotherm tests that were conducted with higher initial concentrations. This indicated that the absolute foam index was not equilibrium. Rather, the test is dynamic where the time of the test played an important role in the results. Even though the absolute foam index was not an equilibrium condition, a correlation was made between the absolute foam index and adsorption isotherms for fly ash of 0.17 to 10.5% LOI. Several batches of mortars were mixed for the same fly ash type increasing only the AEA concentration (dosage) in each subsequent batch. Mortar air test results for each batch showed for each increase in AEA concentration, air contents increased until a point where the next increase in AEA concentration resulted in no increase in air content. This was maximum air content that could be achieved by the particular mortar system; the system reached its air capacity at the saturation limit. This concentration of AEA was compared to the critical micelle concentration (CMC) for the AEA and the absolute foam index.