LARGE SCALE LANDSLIDE ON THE ONTONAGON RIVER, MICHIGAN

In 2003, a large landslide occurred along the Ontonagon River, located in the Upper Peninsula of Michigan, and adjacent to US-45 in Ontonagon County. The failure took place during the springtime, when the river reached a peak discharge that was the second highest on record. The volume of the slide has been estimated to be approximately 1,400,000 cubic yards. The colluvium blocked the river, forcing a new channel to be carved around the debris. The landslide consisted of a silt layer at its base, overlain by a coarsening upward sand sequence, and finally a varved glacio-lacustrine clay with sparse dropstone inclusions making up the upper section of hillside.

PROBABILISTIC MODELING OF RAINFALL INDUCED LANDSLIDE HAZARD ASSESSMENT IN SAN JUAN LA LAGUNA, SOLOLÁ, GUATEMALA

The municipality of San Juan La Laguna, Guatemala is home to approximately 5,200 people and located on the western side of the Lake Atitlán caldera. Steep slopes surround all but the eastern side of San Juan. The Lake Atitlán watershed is susceptible to many natural hazards, but most predictable are the landslides that can occur annually with each rainy season, especially during high-intensity events. Hurricane Stan hit Guatemala in October 2005; the resulting flooding and landslides devastated the Atitlán region. Locations of landslide and non-landslide points were obtained from field observations and orthophotos taken following Hurricane Stan. This study used data from multiple attributes, at every landslide and non-landslide point, and applied different multivariate analyses to optimize a model for landslides prediction during high-intensity precipitation events like Hurricane Stan. The attributes considered in this study are: geology, geomorphology, distance to faults and streams, land use, slope, aspect, curvature, plan curvature, profile curvature and topographic wetness index. The attributes were pre-evaluated for their ability to predict landslides using four different attribute evaluators, all available in the open source data mining software Weka: filtered subset, information gain, gain ratio and chi-squared. Three multivariate algorithms (decision tree J48, logistic regression and BayesNet) were optimized for landslide prediction using different attributes. The following statistical parameters were used to evaluate model accuracy: precision, recall, F measure and area under the receiver operating characteristic (ROC) curve. The algorithm BayesNet yielded the most accurate model and was used to build a probability map of landslide initiation points. The probability map developed in this study was also compared to the results of a bivariate landslide susceptibility analysis conducted for the watershed, encompassing Lake Atitlán and San Juan. Landslides from Tropical Storm Agatha 2010 were used to independently validate this study’s multivariate model and the bivariate model. The ultimate aim of this study is to share the methodology and results with municipal contacts from the author's time as a U.S. Peace Corps volunteer, to facilitate more effective future landslide hazard planning and mitigation.

A STUDY OF A LANDSLIDE IN BIRMINGHAM, ALABAMA

In 1988 a landslide occurred at a construction site in Birmingham, Alabama in which a portion of the construction site required excavating a rock slope with a group of apartments that were located at the top of the slope. During construction, two separate landslides occurred causing one and half of the apartment buildings to collapse downslope. The slope failure was investigated by two firms. One firm investigated the site conditions and the second firm investigated the design of the cut slope. The main concerns in the investigation were (1) the lack of consideration for the existing joint system, (2) using averaged the strength parameters, (3) the possibility of damaging the slope with blasting, and (4) the potential that there were underground mines at the site. The Rocscience program RocPlane was used to model the in situ conditions and the excavation. The model showed that the joint system’s pore water pressure was most likely the main factor in the failure.

Eruptive history of Maderas volcano using new 40Ar/39Ar ages and geochemical analyses

Maderas volcano is a small, andesitic stratovolcano located on the island of Ometepe, in Lake Nicaragua, Nicaragua with no record of historic activity. Twenty-one samples were collected from lava flows from Maderas in 2010. Selected samples were analyzed for whole-rock geochemical data using ICP-AES and/or were dated using the 40Ar/39Ar method. The results of these analyses were combined with previously collected data from Maderas as well as field observations to determine the eruptive history of the volcano and create a geologic map. The results of the geochemical analyses indicate that Maderas is a typical Central American andesitic volcano similar to other volcanoes in Nicaragua and Costa Rica and to its nearest neighbor, Concepción volcano. It is different from Concepción in one important way – higher incompatible elements. Determined age dates range from 176.8 ± 6.1 ka to 70.5 ± 6.1 ka. Based on these ages and the geomorphology of the volcano which is characterized by a bisecting graben, it is proposed that Maderas experienced two clear generations of development with three separate phases of volcanism: initial build-up of the older cone, pre-graben lava flows, and post-graben lava flows. The ages also indicate that Maderas is markedly older than Concepción which is historically active. Results were also analyzed regarding geologic hazards. The 40Ar/39Ar ages indicate that Maderas has likely been inactive for tens of thousands of years and the risk of future volcanic eruptions is low. However, earthquake, lahar and landslide hazards exist for the communities around the volcano. The steep slopes of the eroded older cone are the most likely source of landslide and lahar hazards.

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.

A COMPARISON OF PIXEL-BASED VERSUS OBJECT ORIENTED ANALYSIS OF LANDSLIDES USING HISTORICAL REMOTE SENSING DATA

With recent advances in remote sensing processing technology, it has become more feasible to begin analysis of the enormous historic archive of remotely sensed data. This historical data provides valuable information on a wide variety of topics which can influence the lives of millions of people if processed correctly and in a timely manner. One such field of benefit is that of landslide mapping and inventory. This data provides a historical reference to those who live near high risk areas so future disasters may be avoided. In order to properly map landslides remotely, an optimum method must first be determined. Historically, mapping has been attempted using pixel based methods such as unsupervised and supervised classification. These methods are limited by their ability to only characterize an image spectrally based on single pixel values. This creates a result prone to false positives and often without meaningful objects created. Recently, several reliable methods of Object Oriented Analysis (OOA) have been developed which utilize a full range of spectral, spatial, textural, and contextual parameters to delineate regions of interest. A comparison of these two methods on a historical dataset of the landslide affected city of San Juan La Laguna, Guatemala has proven the benefits of OOA methods over those of unsupervised classification. Overall accuracies of 96.5% and 94.3% and F-score of 84.3% and 77.9% were achieved for OOA and unsupervised classification methods respectively. The greater difference in F-score is a result of the low precision values of unsupervised classification caused by poor false positive removal, the greatest shortcoming of this method.