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.

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.

Federal and state policy influence on woody biomass utilization

A considerable portion of public lands in the United States is at risk of uncharacteristically severe wildfires due to a history of fire suppression. Wildfires already have detrimental impacts on the landscape and on communities in the wildland-urban interface (WUI) due to unnatural and overstocked forests. Strategies to mitigate wildfire risk include mechanical thinning and prescribed burning in areas with high wildfire risk. The material removed is often of little or no economic value. Woody biomass utilization (WBU) could offset the costs of hazardous fuel treatments if removed material could be used for wood products, heat, or electricity production. However, barriers due to transportation costs, removal costs, and physical constraints (such as steep slopes) hinder woody biomass utilization. Various federal and state policies attempt to overcome these barriers. WBU has the potential to aid in wildfire mitigation and meet growing state mandates for renewable energy. This research utilizes interview data from individuals involved with on-the-ground woody biomass removal and utilization to determine how federal and state policies influence woody biomass utilization. Results suggest that there is not one over-arching policy that hinders or promotes woody biomass utilization, but rather woody biomass utilization is hindered by organizational constraints related to time, cost, and quality of land management agencies’ actions. However, the use of stewardship contracting (a hybrid timber sale and service contract) shows promise for increased WBU, especially in states with favorable tax policies and renewable energy mandates. Policy recommendations to promote WBU include renewal of stewardship contracting legislations and a re-evaluation of land cover types suited for WBU. Potential future policies to consider include the indirect role of carbon dioxide emission reduction activities to promote wood energy and future impacts of air quality regulations.

Wolf habitat selection at the territory level : seasonal and interannual variation and influence on reproductive success

Habitat selection has been one of the main research topics in ecology for decades. Nevertheless, many aspects of habitat selection still need to be explored. In particular, previous studies have overlooked the importance of temporal variation in habitat selection and the value of including data on reproductive success in order to describe the best quality habitat for a species. We used data collected from radiocollared wolves in Yellowstone National Park (USA), between 1996 and 2008, to describe wolf habitat selection. In particular, we aimed to identify i) seasonal differences in wolf habitat selection, ii) factors influencing interannual variation in habitat selection, and iii) the effect of habitat selection on wolf reproductive success. We used probability density functions to describe wolf habitat use and habitat coverages to represent the habitat available to wolves. We used regression analysis to connect habitat use with habitat characteristics and habitat selection with reproductive success. Our most relevant result was discovering strong interannual variability in wolf habitat selection. This variability was in part explained by pack identity and differences in litter size and leadership of a pack between two years (summer) and in pack size and precipitation (winter). We also detected some seasonal differences. Wolves selected open habitats, intermediate elevations, intermediate distances from roads, and avoided steep slopes in late winter. They selected areas close to roads and avoided steep slopes in summer. In early winter, wolves selected wetlands, herbaceous and shrub vegetation types, and areas at intermediate elevation and distance from roads. Surprisingly, the habitat characteristics selected by wolves were not useful in predicting reproductive success. We hypothesize that interannual variability in wolf habitat selection may be too strong to detect effects on reproductive success. Moreover, prey availability and competitor pressure may also have an influence on wolf reproductive success, which we did not assess. This project demonstrated how important temporal variation is in shaping patterns of habitat selection. We still believe in the value of running long-term studies, but the effect of temporal variation should always be taken into account.

Topographic and meteorological influences on spatial scaling of heavy convective rainfall in mountainous regions

Characterizing the spatial scaling and dynamics of convective precipitation in mountainous terrain and the development of downscaling methods to transfer precipitation fields from one scale to another is the overall motivation for this research. Substantial progress has been made on characterizing the space-time organization of Midwestern convective systems and tropical rainfall, which has led to the development of statistical/dynamical downscaling models. Space-time analysis and downscaling of orographic precipitation has received less attention due to the complexities of topographic influences. This study uses multiscale statistical analysis to investigate the spatial scaling of organized thunderstorms that produce heavy rainfall and flooding in mountainous regions. Focus is placed on the eastern and western slopes of the Appalachian region and the Front Range of the Rocky Mountains. Parameter estimates are analyzed over time and attention is given to linking changes in the multiscale parameters with meteorological forcings and orographic influences on the rainfall. Influences of geographic regions and predominant orographic controls on trends in multiscale properties of precipitation are investigated. Spatial resolutions from 1 km to 50 km are considered. This range of spatial scales is needed to bridge typical scale gaps between distributed hydrologic models and numerical weather prediction (NWP) forecasts and attempts to address the open research problem of scaling organized thunderstorms and convection in mountainous terrain down to 1-4 km scales.

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.

Cut slope design for the access to an underground copper mine

Slope stability analysis is a major area of research in geotechnical engineering. That being said, very little is written in the geotechnical engineering literature on the design of box-cuts. The goal of this thesis will be to investigate the proper design of a boxcuts, and to design a box-cut for access to an underground copper mine. Issues that need to be considered in the box-cut design include, long term dewatering design, slope stability analysis, and erosion control. The soils at the project site were extremely low permeability, as a result a system of ejectors was designed both to improve the stability of the slopes and prevent flooding. Based on the results of limit equilibrium analysis and finite element analysis, a slope design of two horizontal on one vertical was selection, with a rock fill buttress providing reinforcement. Finally, Michigan DOT standards for seeding were used to provide erosion control

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.