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.

Water chemistry controls on minor elemental concentrations and strontium isotopic compositions of freshwater biogenic carbonate in northeastern Washington streams

Tracking the movement of migratory freshwater fish is essential to those invested in rebuilding declining fish populations. Using strontium isotopic signatures to match calcified fish tissues to streams where fish spawn is a useful method of tracking migratory fish where physical tracking methods such as radio, acoustic, or external tags, have proven unsuccessful. In this study, we develop tools to practice this method of tracking fish in Lake Roosevelt and its upstream tributaries in Washington State by analyzing the elemental concentrations and 87Sr/86Sr ratios of water samples, and mussel shell samples. This study evaluates whether mussel shells act as an appropriate proxy for water chemistry by comparing the 87Sr/86Sr isotope ratios of water samples to the 87Sr/86Sr isotope ratios of mussel shells sampled from the same, or nearby, locations. We compare concentrations of Ba, Ca, Cd, Cu, Fe, Mg, Pb, Sr, and U in the water and mussel shell samples to determine the feasibility of using mussel shells as a proxy for water chemistry. If it is determined that the concentrations of these elements in mussel shells reflect that of the surrounding water composition, the elemental composition of mussel shells can be compared to that of calcified tissues in fish, such as otoliths, to infer the location of the natal stream. We report analyses of water and mussel shell samples collected from Lake Roosevelt, Sanpoil River, Spokane River, Colville River, Kettle River, Pend Oreille River, Kootenay River, and Columbia River in Washington State. Each of these rivers is a tributary to Lake Roosevelt, and each flows through different geologic units. We hypothesize that the differences in the rock units of each stream’s watershed are reflected in the elemental concentrations and strontium isotopic ratios of water in each stream and in the lake. We also hypothesize that the composition of the mussel shells will match the composition of the water samples, therefore allowing us to use the mussel shells as a proxy for local water chemistry. Additionally, we hypothesize that the composition of the mussel shells will vary by location, and that we will be able to then infer where a fish is from by matching the composition of the fish in question to the mussels we have analyzed. We found that 87Sr/86Sr values for water and mussel hinge samples collected from tributaries east of Lake Roosevelt are significantly higher than the 87Sr/86Sr values for samples collected from tributaries west of Lake Roosevelt with averages of 0.7235 and 0.7089, respectively. The average 87Sr/86Sr ratios for water and mussel hinge samples collected within Lake Roosevelt is 0.7158, which is between the averages for samples collected east and west of the lake. Generally, older rocks are exposed on the east side of the lake, and younger rocks on the west side of the lake, so our 87Sr/86Sr values support the hypothesis that geologic units are a primary control on water chemistry, and that tributary compositions mix to form an average weighed by flow in Lake Roosevelt. The 87Sr/86Sr values for water and mussel shell samples collected from the same locations have a strong, positive linear correlation, suggesting that mussel shell 87Sr/86Sr ratios reflect the 87Sr/86Sr ratios of the ambient water. With these data, we can distinguish between different streams and the lake, but cannot distinguish between samples from within the same stream or within Lake Roosevelt. The Sr:Ca and Fe:Ca ratios of water samples show positive correlations with mussel shell compositions, with R2 values of 0.82 and 0.52, respectively. Ratios of Mg, Ba, Cu, Cd, Pb, and U to Ca showed little or no positive correlation between water and mussel shell samples. The elemental concentration data collected for this study do not demonstrate whether a correlation between elemental ratios in water samples and elemental ratios in mussel shell samples collected from the same location exists. Positive Sr:Ca and Fe:Ca correlations for water versus mussel shell samples indicate that perhaps for some elements, the composition of mussel shells are representative of the composition of ambient water. Using elemental concentration ratios to complement 87Sr/86Sr isotopic data may enhance our ability to identify correlations between water and mussel shell samples, and ultimately between mussel shell and otolith samples. The hinge part of a mussel shell may be used as a proxy for local water composition because the mussel shell composition reflects that of the local ambient water. The hinge of the mussel has the same composition as the whole mussel shell. We measured variation of 87Sr/86Sr ratios in the water among different streams and Lake Roosevelt. The 87Sr/86Sr values for samples collected in tributaries east of Lake Roosevelt, which erode older rocks, are higher for mussel shell and water samples than the average 87Sr/86Sr values for mussel shell and water samples collected in tributaries west of Lake Roosevelt, which flow through younger rocks.

Surface Elevation Change Processes in the Snohomish River Estuary, Washington

Estuaries provide crucial ecosystem functions and contain significant socio-economic value. Within Washington State, estuaries supply rearing habitat for juvenile salmon during their transition period from freshwater to open sea. In order to properly manage wetland resources and restore salmon habitat, the mechanisms through which estuaries evolve and adapt to pressures from climate change, most notably eustatic sea level rise, must be understood. Estuaries maintain elevation relative to sea level rise through vertical accretion of sediment. This report investigates the processes that contribute to local surface elevation change in the Snohomish Estuary, conveys preliminary surface elevation change results from RTK GPS monitoring, and describes how surface elevation change will be monitored with a network of RSET-MH’s. Part of the tidal wetlands within the Snohomish River Estuary were converted for agricultural and industrial purposes in the 1800’s, which resulted in subsidence of organic soils and loss of habitat. The Tulalip Tribes, the National Oceanic and Atmospheric Administration (NOAA), Northwest Indian Fisheries Commission (NWIFC), and the Environmental Protection Agency (EPA) are conducting a large-scale restoration project to improve ecosystem health and restore juvenile salmon habitat. A study by Crooks et al. (2014) used 210Pb and carbon densities within sediment cores to estimate wetland re-building capacities, sediment accretion rates, and carbon sequestration potential within the Snohomish Estuary. This report uses the aforementioned study in combination with research on crustal movement, tidal patterns, sediment supply, and sea level rise predictions in the Puget Sound to project how surface elevation will change in the Snohomish Estuary with respect to sea level rise. Anthropogenic modification of the floodplain has reduced the quantity of vegetation and functional connectivity within the Snohomish Estuary. There have been losses up to 99% in vegetation coverage from historic extents within the estuary in both freshwater and mesohaline environments. Hydrographic monitoring conducted by NOAA and the Tulalip Tribe shows that 85% of the historic wetland area is not connected to the main stem of the Snohomish (Jason Hall 2014, unpublished data, NOAA). As vegetation colonization and functional connectivity of the floodplains of the Snohomish estuary is re-established through passive and active restoration, sediment transport and accretion is expected to increase. Under the Intergovernmental Panel on Climate Change (IPCC) “medium- probability” scenario sea level is projected to rise at a rate of 4.28 mm/year in the Puget Sound. Sea level rise in the Snohomish Estuary will be exacerbated from crustal deformation from subsidence and post-glacial rebound, which are measured to be -1.4 mm/year and -0.02 mm/year, respectively. Sediment accretion rates calculated by Crooks et al. (2014) and RTK GPS monitoring of surface elevation change of the Marysville Mitigation site from 2011-2014 measured vertical accretion rates that range from -48-19 mm/year and have high spatial variability. Sediment supply is estimated at 490 thousand tons/year, which may be an under-estimate because of the exclusion of tidal transport in this value. The higher rates of sediment accretion measured in the Snohomish Estuary suggest that the Snohomish will likely match or exceed the pace of sea level rise under “medium-probability” projections. The network of RSET-MH instruments will track surface elevation change within the estuary, and provide a more robust dataset on rates of surface elevation change to quantify how vertical accretion and subsidence are contributing to surface elevation change on a landscape scale.

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.

Exploring the spatial variation of tides for seafloor mapping data corrections

Senior thesis written for Oceanography 445

Assessing Historical Thinking at a State History Museum

Thesis (Master's)--University of Washington, 2016-06

A Study of the Association of Maternal Occupations and Deliveries by Cesarean Section for Infants Born in Washington, 2011-2013

Thesis (Master's)--University of Washington, 2016-06

Economic Evaluations of State Policies

Thesis (Ph.D.)--University of Washington, 2016-08

I-82 Prosser Landslide Investigation Benton City, WA

The scope of this technical report is to establish the mechanisms by which the eastbound lanes of Interstate 82 at mile post (MP) 91.9 near Benton City continue to deform. Within the Washington State Department of Transportation (WSDOT), the area is known as the Prosser Landslide and has been an ongoing concern since the 1980s. Results from previous technical investigations have been conflicted or inconclusive as to whether landslide movement persists beneath or through the shear key-buttress or that pavement distress is related to swelling of a clay-rich unit that underlies the slope and interstate. For this report, the following steps were taken. First, I conducted a desk review of archived reports, memos, data, and drill logs from the original construction of I-82 and previous geotechnical investigations commissioned by WSDOT. Findings of this desk review are reported in Part III. Second, WSDOT drillers drilled two new boreholes at the Prosser Landslide site above the buttress and instrumentation was installed within the boreholes. Borehole logs produced from the 2013 drilling can be found in Appendix A of this report. Material retrieved from the suspected failure zone during drilling was tested at the WSDOT Materials Lab by WSDOT personnel for its mechanical properties including Atterberg limits, grain-size analysis, and residual shear strength (Appendix B). Samples were also analyzed for mineral content using X -ray powder diffraction (XRD). These data and observations are reported in Part III and Appendix C. Finally, using drill logs produced by WSDOT from the latest drilling and from historic drilling campaigns, I constructed a 2-D geologic model of the landslide site. This model is the basis for slope stability analysis reported in Part IV and Appendix D. This study concludes that the deformation observed in the eastbound lanes of I-82 could be the result of continued landslide movement, despite previous remediation efforts.

A comparison of landslide inventories along the shoreline of Steamboat Island Peninsula, Thurston County, WA: office based protocol vs. field based protocol

The mountain ranges and coastlines of Washington State have steep slopes, and they are susceptible to landslides triggered by intense rainstorms, rapid snow melts, earthquakes, and rivers and waves removing slope stability. Over a 30-year timespan (1984-2014 and includes State Route (SR) 530), a total of 28 deep-seated landslides caused 300 million dollars of damage and 45 deaths (DGER, 2015). During that same timeframe, ten storm events triggered shallow landslides and debris flows across the state, resulting in nine deaths (DGER, 2015). The loss of 43 people, due to the SR 530 complex reactivating and moving at a rate and distance unexpected to residents, highlighted the need for an inventory of the stateís landslides. With only 13% of the state mapped (Lombardo et al., 2015), the intention of this statewide inventory is to communicate hazards to citizens and decision makers. In order to compile an accurate and consistent landslide inventory, Washington needs to adopt a graphic information system (GIS) based mapping protocol. A mapping protocol provides consistency for measuring and recording information about landslides, including such information as the type of landslide, the material involved, and the size of the movement. The state of Oregon shares similar landslide problems as Washington, and it created a GIS-based mapping protocol designed to inform its residents, while also saving money and reducing costly hours in the field (Burns and Madin, 2009). In order to determine if the Oregon Department of Geology and Mineral Industries (DOGAMI) protocol, developed by Burns and Madin (2009), could serve as the basis for establishing Washingtonís protocol, I used the office-based DOGAMI protocol to map landslides along a 40-50 km (25-30 mile) shoreline in Thurston County, Washington. I then compared my results to the field-based landslide inventory created in 2009 by the Washington Division of Geology and Earth Resources (DGER) along this same shoreline. If the landslide area I mapped reasonably equaled the area of the DGER (2009) inventory, I would consider the DOGAMI protocol useful for Washington, too. Utilizing 1m resolution lidar flown for Thurston County in 2011 and a GIS platform, I mapped 36 landslide deposits and scarp flanks, covering a total area of 879,530 m2 (9,467,160 ft2). I also found 48 recent events within these deposits. With an exception of two slides, all of the movements occurred within the last fifty years. Along this same coastline, the DGER (2009) recorded 159 individual landslides and complexes, for a total area of 3,256,570 m2 (35,053,400 ft2). At a first glance it appears the DGER (2009) effort found a larger total number and total area of landslides. However, in addition to their field inventory, they digitized landslides previously mapped by other researchers, and they did not field confirm these landslides, which cover a total area of 2,093,860 m2 (22,538,150 ft2) (DGER, 2009). With this questionable landslide area removed and the toes and underwater landslides accounted for because I did not have a bathymetry dataset, my results are within 6,580 m2 (70,840 ft2) of the DGERís results. This similarity shows that the DOGAMI protocol provides a consistent and accurate approach to creating a landslide inventory. With a few additional modifications, I recommend that Washington State adopts the DOGAMI protocol. Acquiring additional 1m lidar and adopting a modified DOGAMI protocol poises the DGER to map the remaining 87% of the state, with an ultimate goal of informing citizens and decision makers of the locations and frequencies of landslide hazards on a user-friendly GIS platform.

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.

Macro and micro plastics in an urbanized and non-urbanized fjord estuary in the Northeast Pacific Ocean

Senior thesis written for Oceanography 445

Hydroaccoustic pinging as a diver tool for underwater navigation in surveying Puget Sound Eelgrass (Zostera marina) densities

Senior thesis written for Oceanography 445

Mapping Whidbey Island and Possession Sound landslide susceptibility using an additive linear model

Senior thesis written for Oceanography 445

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

Senior thesis written for Oceanography 445