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.

A Comparative Analysis of Two Slug Test Methods in Puget Lowland Glacio-Fluvial Sediments near Coupeville, WA

Two different slug test field methods are conducted in wells completed in a Puget Lowland aquifer and are examined for systematic error resulting from water column displacement techniques. Slug tests using the standard slug rod and the pneumatic method were repeated on the same wells and hydraulic conductivity estimates were calculated according to Bouwer & Rice and Hvorslev before using a non-parametric statistical test for analysis. Practical considerations of performing the tests in real life settings are also considered in the method comparison. Statistical analysis indicates that the slug rod method results in up to 90% larger hydraulic conductivity values than the pneumatic method, with at least a 95% certainty that the error is method related. This confirms the existence of a slug-rod bias in a real world scenario which has previously been demonstrated by others in synthetic aquifers. In addition to more accurate values, the pneumatic method requires less field labor, less decontamination, and provides the ability to control the magnitudes of the initial displacement, making it the superior slug test procedure.

Physical Characteristics and Estimated Hydraulic Conductivity from Grain-Size Distribution for Vashon Till, South-Central Puget Lowland, Washington

Population growth, urban development, and increased commercial and industrial activity in the south-central Puget Lowlands of Washington State has led to an increased demand for groundwater. The Vashon till is a glacially consolidated, low-permeability unit comprising unstratified clay, silt, cobbles and boulders with ubiquitous coarse-grained lenses and is an extensive surficial unit throughout the south-central Puget Lowland. Thus, understanding the physical and hydrological characteristics – specifically, the hydraulic conductivity – of this unit is a necessary component of a groundwater model. This study provides (1) a record of the physical characteristics of Vashon till deposits within the study area; and (2) an estimate of the highest, lowest, and average value of saturated hydraulic conductivity based on the grain-size distribution of Vashon till samples collected from six field sites in the Puyallup River Watershed. Analysis shows that the average moisture content ranges between about 1 and 6%, average dry bulk density is about 2.20 g/cm3, and average porosity is about 17%. Grain-size distributions show that half of the samples analyzed are well graded, while the other half is poorly graded. Grain-size distributions also show an average d10 value of about 0.20 mm, and average ff values ≤ 16%, which are key values in estimating the saturated hydraulic conductivity of over-consolidated glacial deposits. Based on these observed values, the estimates of hydraulic conductivity range from a minimum of 0.02 m/d to a maximum of 1.38 m/d in within the general Vashon till.

Stratigraphic analysis of the southwestern quadrant of the Cedar Hills Regional Landfill: Incorporating recent borehole explorations and radiometric dating with previous investigations.

This report presents the results of stratigraphic analysis of the southwestern quadrant of the Cedar Hills Regional Landfill (CHRLF). My report was intended to incorporate the recent Area 8 borehole data into the pre-existing analyses. This analysis was conducted during the preparation of the Area 8 Hydrogeologic Report, but is my independent investigation and does not represent the opinion of UEC or their associates. The CHRLF, in Maple Valley, WA, south of Squak Mountain, is a municipal solid waste landfill that has been in operation since the 1960s. A network of borings, the product of previous investigations, exists for the study area. I utilized the compiled boring logs, previous investigations, and the recently acquired data to produce a series of interpretative cross-sections for the study area. I recognized 9 distinct stratigraphic units, including fill. My interpreted stratigraphic units are similar to those identified in previous investigations such as the Area 7 Hydrogeologic investigation (HDR Engineering and Associates, 2008). These units include pre-Olympia aged non-glacial alluvium, glacial alluvium, and glacial till. Additionally, younger, Vashon-aged deposits of glacial till, recessional outwash, recessional lacustrine, and ice-contact were observed. An isolated “till-like” deposit was observed below the Vashon till. This could possibly represent an older till as mapped by Sweet Edwards (1985) and Booth (1995). I cite the continuity of the lower contact of the Vashon till (Unit 5, Table 2) and the upper contact pre-Vashon non-glacial fluvial deposits (Unit 9, Table 2) as evidence that faults or other structural features do not offset the deposits in the study area. This conclusion supports the findings of the pre-existing body of work within the landfill property and the nearby Queen City Farms property.

A Study of Shallow Groundwater Migration through Anthropogenic Fill along a Seawall

This report provides the findings and opinions of a historical document review, hydraulic balance calculation, and proposed additional study for a property that was historically used as a bulk petroleum storage and distribution facility. The property lies along the base, west, of a heavily vegetated bluff with a tidally influenced body of water west-adjacent to the property. The western portion of the property is bounded by a seawall spanning approximately 3,200 linear feet trending north-south. The seawall’s construction details are not known, save for a 225-foot section of driven sheet pile wall located within the northern portion of the property’s seawall. Due to the presence of petroleum hydrocarbons in soil and groundwater at the property, a cleanup action for the property will likely be overseen by the state regulatory agency. The property is currently undergoing remedial investigation in an effort to identify the lateral and vertical extent in which contaminants at the property have come to be located, also known as the “site” as defined by the Model Toxics Control Act (MTCA). The majority of the property bounded within the seawall area has been characterized; however, the shoreline sediments located immediately west-adjacent of the seawall have not been properly delineated. Identifying the bounds of the site to the west within sediment is pivotal for the purposes of the remedial investigation. Since the west-adjacent shoreline is so extensive, conducting a complete sediment sampling event along the entire shoreline would be cost-prohibitive due to analytical costs and logistical issues at the property. Because of the extensive nature of the shoreline, it would greatly benefit the client and project to focus sampling efforts at areas of greater risk for contaminants along the shoreline by identifying potential preferential pathways for contaminants to migrate off of the property and into adjacent shoreline sediments. The review of historical studies of the property yielded some useful information; however much of the findings included within the historical studies were lacking original raw data, therefore limiting the information obtained. The calculated hydraulic balance for the property yielded a relatively large surplus of recharge to groundwater after precipitation events, reinforcing the concept that contaminant have potentially historically, and currently, been migrating into the adjacent shoreline through preferential pathways along the seawall. Due to the limitation within the historical studies for the property as well as the groundwater recharge identified in the hydraulic balance, an additional study was proposed in an effort to provide additional aquifer characteristics along the seawall, and the ability to observe flow propagation at and proximate to the seawall in two-dimensions through time without the need to piece separate studies together. This proposed study includes a single simultaneous tidal study which comprises select monitoring points along the seawall. This report has identified the need for additional data that can be collected through available avenues for the property based upon the client’s desires and project needs. Ultimately, the proposed additional study is suggested based upon its relatively low capital investment, and ability meet the requirements relevant to the specific project needs and scope. Assuming preferential pathways are identified through the additional study proposed within this report, a representative and cost-effective sediment sampling plan can then be put in place in an effort to define the site.

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.

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.

Development of a Conceptual Site Model for Evaluating Seawater Intrusion under Sea Level Rise Scenarios Using Analytical Methods

Island County is located in the Puget Sound of Washington State and includes several islands, the largest of which is Whidbey Island. Central Whidbey Island was chosen as the project site, as residents use groundwater for their water supply and seawater intrusion near the coast is known to contaminate this resource. In 1989, Island County adopted a Saltwater Intrusion Policy and used chloride concentrations in existing wells in order to define and map “risk zones.” In 2005, this method of defining vulnerability was updated with the use of water level elevations in conjunction with chloride concentrations. The result of this work was a revised map of seawater intrusion vulnerability that is currently in use by Island County. This groundwater management strategy is defined as trigger-level management and is largely a reactive tool. In order to evaluate trends in the hydrogeologic processes at the site, including seawater intrusion under sea level rise scenarios, this report presents a workflow where groundwater flow and discharge to the sea are quantified using a revised conceptual site model. The revised conceptual site model used several simplifying assumptions that allow for first-order quantitative predictions of seawater intrusion using analytical methods. Data from water well reports included lithologic and well construction information, static water levels, and aquifer tests for specific capacity. Results from specific capacity tests define the relationship between discharge and drawdown and were input for a modified Theis equation to solve for transmissivity (Arihood, 2009). Components of the conceptual site model were created in ArcGIS and included interpolation of water level elevation, creation of groundwater basins, and the calculation of net recharge and groundwater discharge for each basin. The revised conceptual site model was then used to hypothesize regarding hydrogeologic processes based on observed trends in groundwater flow. Hypotheses used to explain a reduction in aquifer thickness and hydraulic gradient were: (1) A large increase in transmissivity occurring near the coast. (2) The reduced aquifer thickness and hydraulic gradient were the result of seawater intrusion. (3) Data used to create the conceptual site model were insufficient to resolve trends in groundwater flow. For Hypothesis 2, analytical solutions for groundwater flow under Dupuit assumptions were applied in order to evaluate seawater intrusion under projected sea level rise scenarios. Results indicated that a rise in sea level has little impact on the position of a saltwater wedge; however, a reduction in recharge has significant consequences. Future work should evaluate groundwater flow using an expanded monitoring well network and aquifer recharge should be promoted by reducing surface water runoff.

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.