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.

Determining Spatial Distribution and Physical Properties of the Vashon Advance Outwash near Mountlake Terrace, Washington

This study is aimed at determining the spatial distribution, physical properties, and groundwater conditions of the Vashon advance outwash (Qva) in the Mountlake Terrace, WA area. The Qva is correlative with the Esperance Sand, as defined at its type section; however, local variations in the Qva are not well-characterized (Mullineaux, 1965). While the Qva is a dense glacial unit with low compressibility and high frictional shear strength (Gurtowski and Boirum, 1989), the strength of this unit can be reduced when it becomes saturated (Tubbs, 1974). This can lead to caving or flowing in excavations, and on a larger scale, can lead to slope failures and mass-wasting when intersected by steep slopes. By studying the Qva, we can better predict how it will behave under certain conditions, which will be beneficial to geologists, hydrogeologists, engineers, and environmental scientists during site assessments and early phases of project planning. In this study, I use data from 27 geotechnical borings from previous field investigations and C-Tech Corporation’s EnterVol software to create three-dimensional models of the subsurface geology in the study area. These models made it possible to visualize the spatial distribution of the Qva in relation to other geologic units. I also conducted a comparative study between data from the borings and generalized published data on the spatial distribution, relative density, soil classification, grain-size distribution, moisture content, groundwater conditions, and aquifer properties of the Qva. I found that the elevation of the top of the Qva ranges from 247 to 477 ft. I found that the Qva is thickest where the modern topography is high, and is thinnest where the topography is low. The thickness of the Qva ranges from absent to 242 ft. Along the northern, east-west trending transect, the Qva thins to the east as it rises above a ridge composed of Pre- Vashon glacial deposits. Along the southern, east-west trending transect, the Qva pinches out against a ridge composed of pre-Vashon interglacial deposits. Two plausible explanations for this ridge are paleotopography and active faulting associated with the Southern Whidbey Fault Zone. Further investigations should be done using geophysical methods and the modeling methods described in this study to determine the nature of this ridge. The relative density of the Qva in the study area ranges from loose to very dense, with the loose end of the spectrum probably relating to heave in saturated sands. I found subtle correlations between density and depth. Volumetric analysis of the soil groups listed in the boring logs indicate that the Qva in the study area is composed of approximately 9.5% gravel, 89.3% sand, and 1.2% silt and clay. The natural moisture content ranges from 3.0 to 35.4% in select samples from the Qva. The moisture content appears to increase with depth and fines content. The water table in the study area ranges in elevation from 231.9 to 458 ft, based on observations and measurements recorded in the boring logs. The results from rising-head and falling-head slug tests done at a single well in the study area indicate that the geometric mean of hydraulic conductivity is 15.93 ft/d (5.62 x 10-03 cm/s), the storativity is 3.28x10-03, and the estimated transmissivity is 738.58 ft2/d in the vicinity of this observation well. At this location, there was 1.73 ft of seasonal variation in groundwater elevation between August 2014 and March 2015.

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.

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.

Observations of Joints Present in the Lawton Clay Member of the Vashon Stade, Exposed at a Coastal Bluff in Discovery Park, Seattle, Washington

The southwest-facing coastal bluff present at Discovery Park, Seattle, Washington, displays distinctive joints throughout the exposed Lawton Clay Member. Exhibiting a characteristic local stratigraphy of permeable advance outwash over the impermeable proglacial lacustrine clay, this bluff is located in an area of Seattle at high risk from landslides. This project addressed the relationship between the joints observed at this coastal bluff and the coherency of the bluff as a whole, through remote sensing and field measurements. Aerial drone photography taken of the bluff was processed through a photogrammetry software to produce a 3-dimensional Structure from Motion model, allowing for a digital manipulation and broad examination of the bluff not possible by foot. Stereonet plots produced from these measurements provided insight into patterns of varying joint strike along a horizontal transect of the observed bluff face. Taken together, these two visualizations provided a better picture of the possible chicken-and-egg interaction of the joints and bluff topography; they demonstrated the likelihood that the joint formation at the bluff was most likely to be primarily influenced by the local topography of the bluff over other sources of possible tensional stress in the immediate area.

Analysis of the Hydrogeological Conditions and the Effects of Development on the Hydraulic Head at Redmond Ridge East

Redmond Ridge East (RRE) is a large-scale master plan community in East King County, WA. In this report, I evaluate the spatial variability of the Quaternary Advance Outwash (Qva) at RRE and the time-series data for 16 water wells with the intent to better understand groundwater below the RRE area. I investigate changes between pre- and post-development conditions through the determination of temporal changes in annual water level, annual water level fluctuations, hydraulic head response to precipitation, and ambient drainage of the aquifer. I also perform a basic analysis of the annual aquifer recharge and a determination for the storage through the implementation of the water table fluctuation (WTF) method. Associated Earth Sciences (AESI) was tasked with monitoring the geological and environmental impacts during the development of RRE and collected the data I use in this report. AESI involvement in monitoring began in 1998 and extends to the present. Sixteen wells were identified in the RRE area with adequate temporal data to conduct the analysis. A comparison of the well logs and aquifer testing data allowed local variations in the Qva to be mapped. The WTF was used to determine a range of reasonable specific yield values for locations where the Qva was unconfined. Yearly average of the seasonal water level high and lows, and the fluctuations were quantified. Temporal relationships were established through linear regression. The average water level was found to be increasing in some locations, and the corresponding fluctuations were found to decrease. However, no clear change between pre- and post-development was observed. The response of hydraulic head to precipitation was investigated through an analysis of hydrographs for ten wells. Periods of consistent response and the corresponding precipitation during each period were delineated. A linear relationship between precipitation and water level change was determined. The threshold precipitation under which there is a positive response in the hydraulic head was established. No observable changes were apparent between pre- and post-development conditions. The ambient drainage for the Qva was calculated using recessional periods on the hydrograph. The transmissivity of Qva varies with thickness of the overlying lodgment till and thickness of the Qva, itself. Water level fluctuations observed in the Qva are consistent with regional observations. Localized areas in the Qva display the large 10 foot fluctuations and these anomalies are likely due to a combination of the local variability in the storativity as well as the concentration and channeling of water due to geographical variations in the Qva and the overlying topography. All trends seen in the RRE area remained relatively constant through time. There was no evidence showing an effect of development on the hydraulic head at RRE. This implies that the style and distribution of infiltration has not changed as a result of development, and that any measures in place are properly mitigating the effects of development on the RRE region.