Catastrophic Flooding In North-Central Pennsylvania: Geomorphic Findings From The September 2011 Event And Their Significance To

Tropical Storm Lee produced 25-36 cm of rainfall in north-central Pennsylvania on September 4th through 8th of 2011. Loyalsock Creek, Muncy Creek, and Fishing Creek experienced catastrophic flooding resulting in new channel formation, bank erosion, scour of chutes, deposition/reworking of point bars and chute bars, and reactivation of the floodplain. This study was created to investigate aspects of both geomorphology and sedimentology by studying the well-exposed gravel deposits left by the flood, before these features are removed by humans or covered by vegetation. By recording the composition of gravel bars in the study area and creating lithofacies models, it is possible to understand the 2011 flooding. Surficial clasts on gravel bars are imbricated, but the lack of imbrication and high matrix content of sediments at depth suggests that surface imbrication of the largest clasts took place during hyperconcentrated flow (40-70% sediment concentration). The imbricated clasts on the surface are the largest observed within the bars. The lithofacies recorded are atypical for mixed-load stream lithofacies and more similar to glacial outburst flood lithofacies. This paper suggests that the accepted lithofacies model for mixed-load streams with gravel bedload may not always be useful for interpreting depositional systems. A flume study, which attempted to duplicate the stratigraphy recorded in the field, was run in order to better understand hyperconcentrated flows in the study area. Results from the study in the Bucknell Geology Flume Laboratory indicate that surficial imbrication is possible in hyperconcentrated conditions. After flooding the flume to entrain large amounts of sand and gravel, deposition of surficially imbricated gravel with massive or upward coarsening sedimentology occurred. Imbrication was not observed at depth. These experimental flume deposits support our interpretation of the lithofacies discovered in the field. The sizes of surficial gravel bar clasts show clear differences between chute and point bars. On point bars, gravels fine with increasing distance from the channel. Fining also occurs at the downstream end of point bars. In chute deposits, dramatic fining occurs down the axis of the chute, and lateral grain sizes are nearly uniform. Measuring the largest grain size of sandstone clasts at 8-11 kilometer intervals on each river reveals anomalies in the downstream fining trends. Gravel inputs from bedrock outcrops, tributaries, and erosion of Pleistocene outwash terraces may explain observed variations in grain size along streams either incised into the Appalachian Plateau or located near the Wisconsinan glacial boundary. Atomic Mass Spectrometry (AMS) radiocarbon dating of sediment from recently scoured features on Muncy Creek and Loyalsock Creek returned respective ages of 500 BP and 2490 BP. These dates suggest that the recurrence interval of the 2011 flooding may be several hundred to several thousand years. This geomorphic interval of recurrence is much longer then the 120 year interval calculated by the USGS using historical stream gauge records.

GIS based Watershed Classification in Lao PDR

Every inclined land surface has a potential for soil and water degradation, the seriousness depends on a multitude of parameters such as slope, soil type, geomorphology, rainfall, land use and natural vegetation cover. In Laos this intensified land use leads to reduced vegetation cover, to increased soil erosion, decreasing yield, and finally is likely to influence the hydrological regime. Against this background the Mekong River Commission (MRC) elaborated a spatial explicit Watershed Classification (WSC) for the Lower Mekong Basin. Based on topographic factors derived from a high-resolution Digital Terrain Model, five watershed classes are calculated, giving indication about the sensitivity to resource degradation by soil erosion. The WSC allows spatial priority setting for watershed management and generally supports informed decision making on reconnaissance level. In the conclusions focus is laid on general considerations when GIS techniques are used for spatial decision support in a development context.

Modeling and classifying variable width riparian zones utilizing digital elevation models, flood height data, digital soil data and national wetlands inventory : a new approach for riparian zone delineation

Riparian zones are dynamic, transitional ecosystems between aquatic and terrestrial ecosystems with well defined vegetation and soil characteristics. Development of an all-encompassing definition for riparian ecotones, because of their high variability, is challenging. However, there are two primary factors that all riparian ecotones are dependent on: the watercourse and its associated floodplain. Previous approaches to riparian boundary delineation have utilized fixed width buffers, but this methodology has proven to be inadequate as it only takes the watercourse into consideration and ignores critical geomorphology, associated vegetation and soil characteristics. Our approach offers advantages over other previously used methods by utilizing: the geospatial modeling capabilities of ArcMap GIS; a better sampling technique along the water course that can distinguish the 50-year flood plain, which is the optimal hydrologic descriptor of riparian ecotones; the Soil Survey Database (SSURGO) and National Wetland Inventory (NWI) databases to distinguish contiguous areas beyond the 50-year plain; and land use/cover characteristics associated with the delineated riparian zones. The model utilizes spatial data readily available from Federal and State agencies and geospatial clearinghouses. An accuracy assessment was performed to assess the impact of varying the 50-year flood height, changing the DEM spatial resolution (1, 3, 5 and 10m), and positional inaccuracies with the National Hydrography Dataset (NHD) streams layer on the boundary placement of the delineated variable width riparian ecotones area. The result of this study is a robust and automated GIS based model attached to ESRI ArcMap software to delineate and classify variable-width riparian ecotones.

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.

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.