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.

The Mineralogy and Strength Characteristics of Selected Glaciolacustrine Clays in the Puget Sound Region

The combination of rainy climate, glaciolacustrine clays, and steep topography of the Puget Lowland creates slope stability issues for the regional population. Several glaciolacustrine deposits of laminated silt and clay of different ages contribute to the likelihood of slope failure. The glaciolacustrine deposits are generally wet, range in thickness from absent to >30m, and consist of laminated silt and clay with sand interbeds at the tops and bottoms, sandy laminae throughout the deposits, occasional dropstones and shear zones. The glaciolacustrine deposits destabilize slopes by 1) impeding groundwater flow percolating through overlying glacial outwash sediments, 2) having sandy laminae that lower strength by increasing pore pressure during wet seasons, and 3) increasing the potential for block-style failure because of secondary groundwater pathways such as laminae and vertical fractures. Eight clay samples from six known landslide deposits were analyzed in this study for their mineralogy, clay fraction and strength characteristics. The mineralogy was determined using X-ray Diffractometry (XRD) which revealed an identical mineralogic suite among all eight samples consisting of chlorite, illite and smectite. Nonclay minerals appearing in the X-ray diffractogram include amphibole and plagioclase after removal of abundant quartz grains. Hydrometer tests yielded clay-size fraction percentages of the samples ranging from 10% to 90%, and ring shear tests showed that the angle of residual shear resistance (phi_r) ranged from 11° to 31°. Atterberg limits of the samples were found to have liquid limits ranging from 33 to 83, with plastic limits ranging from 25 to 35 and plasticity indices ranging from 6 to 48. The results of the hydrometer and residual shear strength tests suggest that phi_r varies inversely with the clay-size fraction, but that this relationship was not consistent among all eight samples. The nature of the XRD analysis only revealed the identity of the clay minerals present in the samples, and provided no quantitative information. Thus, the extent to which the mineralogy influenced the strength variability among the samples cannot be determined given that the mineral assemblages are identical. Additional samples from different locations within each deposit along with quantitative compositional analyses would be necessary to properly account for the observed strength variability.

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.

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.