Volcanic processes and human exposure as elements to build a risk model for Volcan de Fuego, Guatemala

The activity of Fuego volcano during the 1999 - 2013 eruptive episode is studied through field, remote sensing and observatory records. Mapping of the deposits allows quantifying the erupted volumes and areas affected by the largest eruptions during this period. A wide range of volcanic processes results in a diversity of products and associated deposits, including minor airfall tephra, rockfall avalanches, lava flows, and pyroclastic flows. The activity can be characterized by long term, low level background activity, and sporadic larger explosive eruptions. Although the background activity erupts lava and ash at a low rate (~ 0.1 m3/s), the persistence of such activity over time results in a significant contribution (~ 30%) to the eruption budget during the studied period. Larger eruptions produced the majority of the volume of products during the studied period, mainly during three large events (May 21, 1999, June 29, 2003, and September 13, 2012), mostly in the form of pyroclastic flows. A total volume of ~ 1.4 x 108 m3 was estimated from the mapped deposits and the estimated background eruption rate. Posterior remobilization of pyroclastic flow material by stream erosion in the highly confined Barranca channels leads to lahar generation, either by normal rainfall, or by extreme rainfall events. A reassessment of the types of products and volumes erupted during the decade of 1970's allows comparing the activity happening since 1999 with the older activity, and suggests that many of the eruptive phenomena at Fuego may have similar mechanisms, despite the differences in scale between. The deposits of large pyroclastic flows erupted during the 1970's are remarkably similar in appearance to the deposit of pyroclastic flows from the 1999 - 2013 period, despite their much larger volume; this is also the case for prehistoric eruptions. Radiocarbon dating of pyroclastic flow deposits suggests that Fuego has produced large eruptions many times during the last ~ 2 ka, including larger eruptions during the last 500 years, which has important hazard implications. A survey was conducted among the local residents living near to the volcano, about their expectations of possible future crises. The results show that people are aware of the risk they could face in case of a large eruption, and therefore they are willing to evacuate in such case. However, their decision to evacuate may also be influenced by the conditions in which the evacuation could take place. If the evacuation represents a potential loss of their livelihood or property they will be more hesitant to leave their villages during a large eruption. The prospect of facing hardship conditions during the evacuation and in the shelters may further cause reluctance to evacuate. A short discussion on some of the issues regarding risk assessment and management through an early warning system is presented in the last chapter.

Observations of volcanic emissions using satellite remote sensing

Volcanoes pose a threat to the human population at regional and global scales and so efficient monitoring is essential in order to effectively manage and mitigate the risks that they pose. Volcano monitoring from space has been possible for over thirty years and now, more than ever, a suite of instruments exists with the capability to observe emissions of gas and ash from a unique perspective. The goal of this research is to demonstrate the use of a range of satellite-based sensors in order to detect and quantify volcanic sulphur dioxide, and to assess the relative performances of each sensor against one another. Such comparisons are important in order to standardise retrievals and permit better estimations of the global contribution of sulphur dioxide to the atmosphere from volcanoes for climate modelling. In this work, retrievals of volcanic sulphur dioxide from a number of instruments are compared, and the individual performances at quantifying emissions from large, explosive volcanic eruptions are assessed. Retrievals vary widely from sensor to sensor, and often the use of a number of sensors in synergy can provide the most complete picture, rather than just one instrument alone. Volcanic emissions have the ability to result significant economic loses by grounding aircraft due to the high risk associated with ash encountering aircraft. As sulphur dioxide is often easier to measure than ash, it is often used as a proxy. This work examines whether this is a reasonable assumption, using the Icelandic eruption in early 2010 as a case study. Results indicate that although the two species are for the most part collocated, separation can occur under some conditions, meaning that it is essential to accurately measure both species in order to provide effective hazard mitigation. Finally, the usefulness of satellite remote sensing in quantifying the passive degassing from Turrialba, Costa Rica is demonstrated. The increase in activity from 2005 – 2010 can be observed in satellite data prior to the phreatic phase of early 2010, and can therefore potentially provide a useful indication of changing activity at some volcanoes.

La Yeguada volcanic complex, western Panama : an assessment of the geologic hazards using new 40AR/39Ar Ages

La Yeguada volcanic complex is one of three Quaternary volcanic centers in Panama, and is located on the southern slope of the Cordillera Central mountain range in western Panama, province of Veraguas. To assess potential geologic hazards, this study focused on the main dome complex near the village of La Laguna and also examined a cinder cone 10 km to the northwest next to the village of Media Luna. Based on newly obtained 40Ar/39Ar ages, the most recent eruption occurred approximately 32 000 years ago at the Media Luna cinder cone, while the youngest dated eruption at the main dome complex occurred 0.357 ± 0.019 Ma, producing the Castillo dome unit. Cerro Picacho is a separate dome located 1.5 km east of the main complex with a date of 4.47 ± 0.23 Ma, and the El Satro Pyroclastic Flow unit surrounds the northern portion of the volcanic complex and has an age of 11.26 ± 0.17 Ma. No Holocene (10 000 years ago to present) activity is recorded at the La Yeguada volcanic complex and therefore, it is unlikely to produce another eruption. The emergence of a new cinder cone is a possibility, but the associated hazards tend to be low and localized, and this does not pose a significant threat to the small communities scattered throughout the area. The main geologic hazard at the La Yeguada volcanic complex is from landslides coming off the many steep slopes.

Probabilistic-based hurricane risk assessment and mitigation considering the potential impacts of climate change

Studies are suggesting that hurricane hazard patterns (e.g. intensity and frequency) may change as a consequence of the changing global climate. As hurricane patterns change, it can be expected that hurricane damage risks and costs may change as a result. This indicates the necessity to develop hurricane risk assessment models that are capable of accounting for changing hurricane hazard patterns, and develop hurricane mitigation and climatic adaptation strategies. This thesis proposes a comprehensive hurricane risk assessment and mitigation strategies that account for a changing global climate and that has the ability of being adapted to various types of infrastructure including residential buildings and power distribution poles. The framework includes hurricane wind field models, hurricane surge height models and hurricane vulnerability models to estimate damage risks due to hurricane wind speed, hurricane frequency, and hurricane-induced storm surge and accounts for the timedependant properties of these parameters as a result of climate change. The research then implements median insured house values, discount rates, housing inventory, etc. to estimate hurricane damage costs to residential construction. The framework was also adapted to timber distribution poles to assess the impacts climate change may have on timber distribution pole failure. This research finds that climate change may have a significant impact on the hurricane damage risks and damage costs of residential construction and timber distribution poles. In an effort to reduce damage costs, this research develops mitigation/adaptation strategies for residential construction and timber distribution poles. The costeffectiveness of these adaptation/mitigation strategies are evaluated through the use of a Life-Cycle Cost (LCC) analysis. In addition, a scenario-based analysis of mitigation strategies for timber distribution poles is included. For both residential construction and timber distribution poles, adaptation/mitigation measures were found to reduce damage costs. Finally, the research develops the Coastal Community Social Vulnerability Index (CCSVI) to include the social vulnerability of a region to hurricane hazards within this hurricane risk assessment. This index quantifies the social vulnerability of a region, by combining various social characteristics of a region with time-dependant parameters of hurricanes (i.e. hurricane wind and hurricane-induced storm surge). Climate change was found to have an impact on the CCSVI (i.e. climate change may have an impact on the social vulnerability of hurricane-prone regions).

Risk analysis of light-frame wood construction due to multiple hazards

Light-frame wood buildings are widely built in the United States (U.S.). Natural hazards cause huge losses to light-frame wood construction. This study proposes methodologies and a framework to evaluate the performance and risk of light-frame wood construction. Performance-based engineering (PBE) aims to ensure that a building achieves the desired performance objectives when subjected to hazard loads. In this study, the collapse risk of a typical one-story light-frame wood building is determined using the Incremental Dynamic Analysis method. The collapse risks of buildings at four sites in the Eastern, Western, and Central regions of U.S. are evaluated. Various sources of uncertainties are considered in the collapse risk assessment so that the influence of uncertainties on the collapse risk of lightframe wood construction is evaluated. The collapse risks of the same building subjected to maximum considered earthquakes at different seismic zones are found to be non-uniform. In certain areas in the U.S., the snow accumulation is significant and causes huge economic losses and threatens life safety. Limited study has been performed to investigate the snow hazard when combined with a seismic hazard. A Filtered Poisson Process (FPP) model is developed in this study, overcoming the shortcomings of the typically used Bernoulli model. The FPP model is validated by comparing the simulation results to weather records obtained from the National Climatic Data Center. The FPP model is applied in the proposed framework to assess the risk of a light-frame wood building subjected to combined snow and earthquake loads. The snow accumulation has a significant influence on the seismic losses of the building. The Bernoulli snow model underestimates the seismic loss of buildings in areas with snow accumulation. An object-oriented framework is proposed in this study to performrisk assessment for lightframe wood construction. For home owners and stake holders, risks in terms of economic losses is much easier to understand than engineering parameters (e.g., inter story drift). The proposed framework is used in two applications. One is to assess the loss of the building subjected to mainshock-aftershock sequences. Aftershock and downtime costs are found to be important factors in the assessment of seismic losses. The framework is also applied to a wood building in the state of Washington to assess the loss of the building subjected to combined earthquake and snow loads. The proposed framework is proven to be an appropriate tool for risk assessment of buildings subjected to multiple hazards. Limitations and future works are also identified.

Understanding volcanic processes using UV camera measurements of sulfur dioxide and coincident infrasound and seismicity

The exsolution of volatiles from magma maintains an important control on volcanic eruption styles. The nucleation, growth, and connectivity of bubbles during magma ascent provide the driving force behind eruptions, and the rate, volume, and ease of gas exsolution can affect eruptive activity. Volcanic plumes are the observable consequence of this magmatic degassing, and remote sensing techniques allow us to quantify changes in gas exsolution. However, until recently the methods used to measure volcanic plumes did not have the capability of detecting rapid changes in degassing on the scale of standard geophysical observations. The advent of the UV camera now makes high sample rate gas measurements possible. This type of dataset can then be compared to other volcanic observations to provide an in depth picture of degassing mechanisms in the shallow conduit. The goals of this research are to develop a robust methodology for UV camera field measurements of volcanic plumes, and utilize this data in conjunction with seismoacoustic records to illuminate degassing processes. Field and laboratory experiments were conducted to determine the effects of imaging conditions, vignetting, exposure time, calibration technique, and filter usage on the UV camera sulfur dioxide measurements. Using the best practices determined from these studies, a field campaign was undertaken at Volcán de Pacaya, Guatemala. Coincident plume sulfur dioxide measurements, acoustic recordings, and seismic observations were collected and analyzed jointly. The results provide insight into the small explosive features, variations in degassing rate, and plumbing system of this complex volcanic system. This research provides useful information for determining volcanic hazard at Pacaya, and demonstrates the potential of the UV camera in multiparameter studies.

SLOPE STABILITY ANALYSIS OF MOUNT MEAGER, SOUTH-WESTERN BRITISH COLUMBIA, CANADA

The Mount Meager Volcanic Complex (MMVC) in south-western British Columbia is a potentially active, hydrothermally altered massif comprising a series of steep, glaciated peaks. Climatic conditions and glacial retreat has led to the further weathering, exposure and de-buttressing of steep slopes composed of weak, unconsolidated material. This has resulted in an increased frequency of landslide events over the past few decades, many of which have dammed the rivers bordering the Complex. The breach of these debris dams presents a risk of flooding to the downstream communities. Preliminary mapping showed there are numerous sites around the Complex where future failure could occur. Some of these areas are currently undergoing progressive slope movement and display features to support this such as anti-scarps and tension cracks. The effect of water infiltration on stability was modelled using the Rocscience program Slide 6.0. The main site of focus was Mount Meager in the south- east of the Complex where the most recent landslide took place. Two profiles through Mount Meager were analysed along with one other location in the northern section of the MMVC, where instability had been detected. The lowest Factor of Safety (FOS) for each profile was displayed and an estimate of the volume which could be generated was deduced. A hazard map showing the inundation zones for various volumes of debris flows was created from simulations using LAHARZ. Results showed the massif is unstable, even before infiltration. Varying the amount of infiltration appears to have no significant impact on the FOS annually implying that small changes of any kind could also trigger failure. Further modelling could be done to assess the impact of infiltration over shorter time scales. The Slide models show the volume of material that could be delivered to the Lillooet River Valley to be of the order of 109 m3 which, based on the LAHARZ simulations, would completely inundate the valley and communities downstream. A major hazard of this is that the removal of such a large amount of material has the potential to trigger an explosive eruption of the geothermal system and renew volcanic activity. Although events of this size are infrequent, there is a significant risk to the communities downstream of the complex.

STABILITY ANALYSIS AND HAZARD ASSESSMENT OF THE NORTHERN SLOPES OF SAN VICENTE VOLCANO IN CENTRAL EL SALVADOR

Geologic hazards affect the lives of millions of people worldwide every year. El Salvador is a country that is regularly affected by natural disasters, including earthquakes, volcanic eruptions and tropical storms. Additionally, rainfall-induced landslides and debris flows are a major threat to the livelihood of thousands. The San Vicente Volcano in central El Salvador has a recurring and destructive pattern of landslides and debris flows occurring on the northern slopes of the volcano. In recent memory there have been at least seven major destructive debris flows on San Vicente volcano. Despite this problem, there has been no known attempt to study the inherent stability of these volcanic slopes and to determine the thresholds of rainfall that might lead to slope instability. This thesis explores this issue and outlines a suggested method for predicting the likelihood of slope instability during intense rainfall events. The material properties obtained from a field campaign and laboratory testing were used for a 2-D slope stability analysis on a recent landslide on San Vicente volcano. This analysis confirmed that the surface materials of the volcano are highly permeable and have very low shear strength and provided insight into the groundwater table behavior during a rainstorm. The biggest factors on the stability of the slopes were found to be slope geometry, rainfall totals and initial groundwater table location. Using the results from this analysis a stability chart was created that took into account these main factors and provided an estimate of the stability of a slope in various rainfall scenarios. This chart could be used by local authorities in the event of a known extreme rainfall event to help make decisions regarding possible evacuation. Recommendations are given to improve the methodology for future application in other areas as well as in central El Salvador.