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.

Interpreting submarine structures for glacial retreat periodicity within Muchalat Inlet in Nootka Sound

Senior thesis written for Oceanography 445

Evaluation of Potential Post-Glacial Faults Using Bare-Earth LiDAR Topographic Data Near Ridley Island, British Columbia

I present results of my evaluation to identify topographic lineaments that are potentially related to post-glacial faulting using bare-earth LiDAR topographic data near Ridley Island, British Columbia. The purpose of this evaluation has been to review bare-earth LiDAR data for evidence of post-glacial faulting in the area surrounding Ridley Island and provide a map of the potential faults to review and possibly field check. My work consisted of an extensive literature review to understand the tectonic, geologic, glacial and sea level history of the area and analysis of bare-earth LiDAR data for Ridley Island and the surrounding region. Ridley Island and the surrounding north coast of British Columbia have a long and complex tectonic and geologic history. The north coast of British Columbia consists of a series of accreted terranes and some post-accretionary deposits. The accreted terranes were attached to the North American continent during subduction of the Pacific Plate between approximately 200 Ma and 10 Ma. The terrane and post-accretionary deposits are metamorphosed sedimentary, volcanic and intrusive rocks. The rocks have experienced significant deformation and been intruded by plutonic bodies. Approximately 10 Ma subduction of the Pacific Plate beneath the North America Plate ceased along the central and north coast of British Columbia and the Queen Charlotte Fault Zone was formed. The Queen Charlotte Fault Zone is a transform-type fault that separates the Pacific Plate from the North America Plate. Within the past 1 million years, the area has experienced multiple glacial/interglacial cycles. The most recent glacial cycle occurred approximately 23,000 to 13,500 years ago. Few Quaternary deposits have been mapped in the area. The vast majority of seismicity around the northwest coast of British Columbia occurs along the Queen Charlotte Fault Zone. Numerous faults have been mapped in the area, but there is currently no evidence to suggest these faults are active (i.e. have evidence for post-glacial surface displacement or deformation). No earthquakes have been recorded within 50 km of Ridley Island. Several small earthquakes (less than magnitude 6) have been recorded within 100 km of the island. These earthquakes have not been correlated to active faults. GPS data suggests there is ongoing strain in the vicinity of Ridley Island. The strain has the potential to be released along faults, but the calculated strain may be a result of erroneous data or accommodated aseismically. Currently, the greatest known seismic hazard to Ridley Island is the Queen Charlotte Fault Zone. LiDAR data for Ridley Island, Digby Island, Lelu Island and portions of Kaien Island, Smith Island and the British Columbia mainland were reviewed and analyzed for evidence of postglacial faulting. The data showed a strong fabric across the landscape with a northwest-southeast trend that appears to mirror the observed foliation in the area. A total of 80 potential post-glacial faults were identified. Three lineaments are categorized as high, forty-one lineaments are categorized as medium and thirty-six lineaments are categorized as low. The identified features should be examined in the field to further assess potential activity. My analysis did not include areas outside of the LiDAR coverage; however faulting may be present there. LiDAR data analysis is only useful for detecting faults with surficial expressions. Faulting without obvious surficial expressions may be present in the study area.

Evaluation of the Relationship Between Saturated Hydraulic Conductivity and Grain-Size Distribution of Fluvio-Glacial Deposits, Puget Lowland, Washington

The relationship between saturated hydraulic conductivity (Ks) and grain-size distribution was evaluated for 49 sites underlain by either glacially over consolidated or normally consolidated fluvio-glacial deposits in the Puget Lowland. A linear regression comprising pairs of grain-size analyses and pilot infiltration tests predicts Ks with a 1 sigma uncertainty of a factor of about 3.5 with 70% of the population variance accounted for. The correlation coefficient R^2 of about 0.90 shows that there is a strong correlation between the grain-size distribution and Ks. In contrast, a widely applied analysis proposed by Massmann (2003) explains only 20% of the population variance for normally consolidated materials with an R^2 of only 0.15. That analysis entirely fails to explain the population variance for over consolidated materials. The method developed in this study is recommended for determination of Ks for fluvio-glacial deposits of the Puget Lowland.

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.

Computing Hydraulic Conductivity from Borehole Infiltration Tests and Grain-Size Distribution in Advance Glacial Outwash Deposits, Puget Lowland, Washington

This study evaluates two methods for estimating a soilís hydraulic conductivity: in-situ infiltration tests and grain-size analyses. There are numerous formulas in the literature that relate hydraulic conductivity to various parameters of the infiltrating medium, but studies have shown that these formulas do not perform well when applied to depositional environments that differ from those used to derive the formulas. Thus, there exists a need to specialize infiltration tests and related grain-size analyses for the Vashon advance outwash in the Puget Lowland. I evaluated 134 infiltration tests and 119 soil samples to find a correlation between grain-size parameters and hydraulic conductivity. This work shows that a constant-head borehole infiltration test that accounts for capillarity with alpha approximately 5m^-1 is an effective method for calculating hydraulic conductivity from our flow tests. Then, by conducting grain-size analysis and applying a multiple linear regression, I show that the hydraulic conductivity can also be estimated by log(K) = 1.906 + 0.102D_10 + 0.039D_60 - 0.034D_90 - 7.952F_fines. This result predicts the infiltration rate with a 95% confidence interval of 20 ft/day. The results of study are for application in the Puget Lowland.

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

Senior thesis written for Oceanography 445

Geological structure and acoustic properties of sediments in Effingham Inlet

Thesis written for Oceanography 445

Identifying a historic shoreline using fine-scale slope variation

Senior thesis written for Oceanography 445