Geology and landslide geomorphology of the Burpee Hills, Skagit County, Washington, USA

Landforms within the Skagit Valley record a complex history of land evolution from Late Pleistocene to the present. Late Pleistocene glacial deposits and subsequent incision by the Skagit River formed the Burpee Hills terrace. The Burpee Hills comprises an approximately 205-m-thick sequence of sediments, including glacio-lacustrine silts and clays, overlain by sandy advance outwash and capped by coarse till, creating a sediment-mantled landscape where mass wasting occurs in the form of debris flows and deep-seated landslides (Heller, 1980; Skagit County, 2014). Landslide probability and location are necessary metrics for informing citizens and policy makers of the frequency of natural hazards. Remote geomorphometric analysis of the site area using airborne LiDAR combined with field investigation provide the information to determine relative ages of landslide deposits, to classify geologic units involved, and to interpret the recent hillslope evolution. Thirty-two percent of the 28-km2 Burpee Hills landform has been mapped as landslide deposits. Eighty-five percent of the south-facing slope is mapped as landslide deposits. The mapped landslides occur predominantly within the advance outwash deposits (Qgav), this glacial unit has a slope angle ranging from 27 to 36 degrees. Quantifying surface roughness as a function of standard deviation of slope provides a relative age of landslide deposits, laying the groundwork for frequency analysis of landslides on the slopes of the Burpee Hills. The south-facing slopes are predominately affected by deep-seated landslides as a result of Skagit River erosion patterns within the floodplain. The slopes eroded at the toe by the Skagit River have the highest roughness coefficients, suggesting that areas with more frequent disturbance at the toe are more prone to sliding or remobilization. Future work including radiocarbon dating and hydrologic-cycle investigations will provide a more accurate timeline of the Burpee Hills hillslope evolution, and better information for emergency management and planners in the future.

Linking surficial geomorphology with vertical structure in high and low energy marine environments

Senior thesis written for Oceanography 445

Vasa Creek Channel Conditions and Slope Stability, Bellevue, Washington

The focus of this report is is the channel conditions at Vasa Creek, Bellevue, Washington, with regard to kokanee habitat and slope stability. This required a geomorphic and geologic assessment of the stream and riparian corridor along Vasa Creek. I focused my efforts in a 720m study-reach just south of I-90 in which City of Bellevue had no information. My assessment is divided into 3 categories: channel morphology, geology, and landslide hazards. I described the channel morphology by determining the gradient of the channel, longitudinal and cross-channel geometries, grain size distribution, embeddedness observations, type of channel reaches present, and the locations of significant in-channel woody-debris, landslides, scarps, landslide debris, and erosional features. This was done by conducting a longitudinal survey, 7 cross-channel surveys, pebble counts, and visual observations with the aid of a GPS device for mapping. I completed my geological assessment using both field observations and borehole data provided by GeoMapNW. Borehole data provided logs of the subsurface material at specific locations. In the field, I interpreted local geology using material in the channel as well as exposures in the adjacent slope. I completed the landslide hazard assessment using GIS methods supplemented by field observations. GIS methods included the use of aerial LiDAR to discern slope values and locations of features. Features of interest include the locations of scarps, landslides, landslide debris, and erosional features which were observed in the field. I classified 4 slope classes using ArcMap10 along with the locations of previously mapped landslides, scarps, and landslide debris. I describe the risk of slope failure according to the Washington Administration Code definition of critical areas (WAC 365-190-120 6a-i). My results are presented in the form of a map suite containing a channel morphology map, geology map, and landslide hazard map. The channel is a free-formed alluvial plane-bed reach with infrequent step-pools with riffles associated with landslide debris that chokes the channel. Overall I found that there is not the potential for kokanee habitat due flashy behavior (sudden high flow events), landslide inundation, and a lack of favorable conditions within the channel. The updated geologic map displays advance outwash deposits and alluvium present within the study-reach, as opposed to exposures of the Blakeley Formation along with other corrections from borehole data interpretations. The landslide hazard map shows that there are areas at high risk for slope failure along the channel that should be looked into further.

Locating the Seattle Fault in Bellevue, WA: Combining geomorphic indicators with borehole data

The Seattle Fault is an active east-west trending reverse fault zone that intersects both Seattle and Bellevue, two highly populated cities in Washington. Rupture along strands of the fault poses a serious threat to infrastructure and thousands of people in the region. Precise locations of fault strands are still poorly constrained in Bellevue due to blind thrusting, urban development, and/or erosion. Seismic reflection and aeromagnetic surveys have shed light on structural geometries of the fault zone in bedrock. However, the fault displaces both bedrock and unconsolidated Quaternary deposits, and seismic data are poor indicators of the locations of fault strands within the unconsolidated strata. Fortunately, evidence of past fault strand ruptures may also be recorded indirectly by fluvial processes and should also be observable in the subsurface. I analyzed hillslope and river geomorphology using LiDAR data and ArcGIS to locate surface fault traces and then compare/correlate these findings to subsurface offsets identified using borehole data. Geotechnical borings were used to locate one fault offset and provide input to a cross section of the fault constructed using Rockworks software. Knickpoints, which may correlate to fault rupture, were found upstream of this newly identified fault offset as well as upstream of a previously known fault segment.

Channel Migration and Avulsion Hazards at Sunshine Point Campground, Mount Rainier National Park, Washington

In November 2006, the flood of record on the upper Nisqually River destroyed part of Sunshine Point Campground in Mount Rainier National Park, Washington. The Nisqually River migrated north and reoccupied five acres of its floodplain; Tahoma Creek partially avulsed into the west floodplain, topping banks of an undersized channel and flooding the campground. I assessed hazards to infrastructure at the old campground location, where the Park proposes to rebuild the remaining campground roads and sites. This assessment focuses on two major hazards: northward Nisqually River migration, which may reincorporate the floodplain into the river destroying infrastructure; and Tahoma Creek avulsions, which may flood the campgroud and deposit sediment burying campground infrastructure. I quantify northward migration by: estimating migration rates and changes to channel width; evaluating river occupation of the pre- and post-2006 campground; and estimating scour depths at revetments protecting the campground. I digitized the Nisqually River channels and channel centerlines from maps and images between 1955 and 2013 into a GIS, which I used to estimate migration rate and river width changes. Centerline migration rates average 9 ft/yr along the length of the Nisqually River study reach; at Sunshine Point lateral migration rates average 11 ft/yr. Maximum migration along the study reach was 19 ft/yr between 2006 and 2009. Greater than average migration rates and channel widths correspond to river confluences and include the Tahoma Creek confluence at Sunshine Point. To determine historical channel locations and the frequency that the river occupied different parts of its floodplain, I digitized the river from maps and images between 1903 and 2013. The Nisqually River flows through Sunshine Point Campground in eight out of 15 historical images. I assess scour at revetments protecting infrastructure from the Nisqually River during a 100-year recurrence interval flood using measured cross-sections. During a 100-year flood, the Nisqually River may scour up to 10 feet below the bed elevation. These scour depths can destabilize critical revetments leaving loose unconsolidated riverbanks exposed to Nisqually River flows. To determine the causes, locations, and frequency of flood hazards from Tahoma Creek avulsions, I field map avulsion channels and compare the results with imagery and channel width changes between 1955 and 2013. Mapped avulsion channels occur with swaths of dead vegetation or nascent vegetation; both dead and recent vegetation are visibly distinct from surrounding vegetation in aerial images. Times of changes to these vegetation anomalies correspond to increases in Tahoma Creek channel width. Avulsions have occurred at least three times in the study period: pre-1955, between 1979 and 1984, and in 2006. The 1984 and 2006 avulsions both occur after increases in Tahoma Creek reach averaged width. The NPS is considering two options to rebuild Sunshine Point Campground, both at the same location. The hazards posed by the Nisqually River and Tahoma Creek at Sunshine Point will affect both construction options equally. Migration hazards to the campground may be reduced by limiting the proposed campground infrastructure to an elevated ridge that has not been occupied by the Nisqually River since 1903. The hazards of damage from migration may be reduced by revetments, which were effective in preventing northward Nisqually River migration in 1959 and 1965. Tahoma Creek avulsions are related increased of Tahoma Creek reach averaged widths, which are near a 58- year maximum, and occurred during a 10-year flood in 1984. The campground may be as susceptible to flooding from avulsions during as little as a 10-year flood. A large avulsion may occur with the next significant Tahoma Creek width increase. Glacial retreat has been shown to increase debris flow activity and increase sediment delivery to Mount Rainier rivers. Increased sediment discharge has been correlated with aggradation, which will further encourage Tahoma Creek avulsions.

Shallow Landslide Risk Assessment of the Westerly Bluffs of the 7.5-minute Freeland Quadrangle, Whidbey Island, Washington

Landslides often occur on slopes rendered unstable by underlying geology, geomorphology, hydrology, weather-climate, slope modifications, or deforestation. Unfortunately, humans commonly exacerbate such unstable conditions through careless or imprudent development practices. Due to local geology, geography, and climatic conditions, Puget Sound of western Washington State is especially landslide-prone. Despite this known issue, detailed analyses of landslide risks for specific communities are few. This study aims to classify areas of high landslide risk on the westerly bluffs of the 7.5 minute Freeland quadrangle based on a combined approach: mapping using LiDAR imagery and the Landform Remote Identification Model (LRIM) to identify landslides, and implementation of the Shallow Slope Stability Model (SHALSTAB) to establish a landslide exceedance probability. The objective is to produce a risk assessment from two shallow landslide scenarios: (1) minimum bluff setback and runout and (2) maximum bluff setback and runout. A simple risk equation that takes into account the probability of hazard occurrence with physical and economic vulnerability (van Westen, 2004) was applied to both scenarios. Results indicate an possible total loss as much as $32.6b from shallow landslides, given a setback of 12 m and a runout of 235 m.

Geomorphic Impacts of the 2013 Colorado Front Range Flood on Black Canyon Creek and North Fork Big Thompson River

In September 2013, the Colorado Front Range experienced a five-day storm that brought record-breaking precipitation to the region. As a consequence, many Front Range streams experienced flooding, leading to erosion, debris flows, bank failures and channel incision. I compare the effects that debris flows and flooding have on the channel bar frequency, frequency and location of wood accumulation, and on the shape and size of the channel along two flood impacted reaches located near Estes Park and Glen Haven, Colorado within Rocky Mountain National Park and Arapaho-Roosevelt National Forest: Black Canyon Creek (BCC) and North Fork Big Thompson River (NFBT). The primary difference between the two study areas is that BCC was inundated by multiple debris flows, whereas NFBT only experienced flooding. Fieldwork consisted of recording location and size of large wood and channel bars and surveying reaches to produce cross-sections. Additional observations were made on bank failures in NFBT and the presence of boulders in channel bars in BCC to determine sediment source. The debris flow acted to scour and incise BCC causing long-term alteration. The post-flood channel cross-sectional area is as much as 7 to 23 times larger than the pre-flood channel, caused by the erosion of the channel bed to bedrock and the elimination of riparian vegetation. Large wood was forced out of the stream channel and deposited outside of the bankfull channel. Flooding in NFBT caused bank erosion and widening that contributed sediment to channel bars, but accomplished little stream-bed scour. As a result, there was relatively little damage to mid-channel and riparian vegetation, and most large wood remained within the wetted channel.

River and Hillslope Response to Localized Uplift Along a Left Bend of the San Andreas Fault, Santa Cruz Mountains, CA

Studying landscape evolution of the Earthís surface is difficult because both tectonic forces and surface processes control its response to perturbation, and ultimately, its shape and form. Researchers often use numerical models to study erosional response to deformation because there are rarely natural settings in which we can evaluate both tectonic activity and topographic response over appropriate time scales (103-105 years). In certain locations, however, geologic conditions afford the unique opportunity to study the relationship between tectonics and topography. One such location is along the Dragonís Back Pressure Ridge in California, where the landscape moves over a structural discontinuity along the San Andreas Fault and landscape response to both the initiation and cessation of uplift can be observed. In their landmark study, Hilley and Arrowsmith (2008) found that geomorphic metrics such as channel steepness tracked uplift and that hillslope response lagged behind that of rivers. Ideal conditions such as uniform vegetation density and similar lithology allowed them to view each basin as a developmental stage of response to uplift only. Although this work represents a significant step forward in understanding landscape response to deformation, it remains unclear how these results translate to more geologically complex settings. In this study, I apply similar methodology to a left bend along the San Andreas Fault in the Santa Cruz Mountains, California. At this location, the landscape is translated through a zone of localized uplift caused by the bend, but vegetation, lithology, and structure vary. I examine the geomorphic response to uplift along the San Andreas Fault bend in order to determine whether predicted landscape patterns can be observed in a larger, more geologically complex setting than the Dragonís Back Pressure Ridge. I find that even with a larger-scale and a more complex setting, geomorphic metrics such as channel steepness index remain useful tools for evaluating landscape evolution through time. Steepness indices in selected streams of study record localized uplift caused by the restraining bend, while hillslope adjustment in the form of landsliding occurs over longer time scales. This project illustrates that it is possible to apply concepts of landscape evolution models to complex settings and is an important contribution to the body of geomorphological study.

The role of large woody debris and riparian forest in channel avulsion in the Carbon River, Mount Rainier National Park, Washington

A specific type of natural log jam in the upper alluvial reach of the Carbon River was found to influence secondary channel avulsion, causing flooding hazards to the adjacent Carbon River Road in the northwest quadrant of Mount Rainier National Park, Washington. The fence-like natural log jam was characterized by large woody debris buttressed horizontally against standing riparian trees (i.e. ìfence railsî and ìfence postî). The objectives of this report are two-fold. First, physical characteristics and spatial distribution were documented to determine the geomorphic controls on the fence-like log jams. Second, the function and timing of the natural log jam in relation to channel avulsion was determined to provide insight into flooding hazards along the Carbon River Road. The fence-like log jams are most abundant in the upper reaches of the Carbon River between 3.0 and 5.5 kilometers from the Carbon Glacier terminus, where longitudinal gradient significantly decreases from about 0.06 to 0.03. Sediment impoundment can occur directly upstream of the fence-like log jam, creating vertical bed elevation difference as high as 1.32 meters, and can form during low magnitude, high frequency flood event (3.5-year recurrence interval). In some locations, headcuts and widening of secondary channel were observed directly to the side of the log jams, suggesting its role in facilitating secondary channel avulsions. Areas along the Carbon River Road more prone to damages from avulsion hazards were identified by coupling locations of the log jams and Relative Water Surface Elevation map created using the 1-meter 2012 Light Detection and Ranging Digital Elevation Map. Ultimately, the results of this report may provide insight to flooding hazards along the Carbon River Road from log jam-facilitated channel avulsion.