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.

An investigation of contaminant distribution in subsurface marine wood waste at Mill-A South Terminal, Port of Everett, Washington

This project investigates the correlation between contaminants and the wood waste present in marine sediments off the shore of the Port of Everett in the former Weyerhaeuser Mill-A pulp mill site. The investigation includes the results of two field studies, which tested contaminant levels in 22 boreholes as well as several surface samples. The contaminants include heavy metals and wood waste byproducts. These results, along with 14 other bore logs, provide the framework for a three-dimensional site model, interpolating the full extent of the depositional units and organic and inorganic chemicals found at Mill-A. The sediments of interest are divided into five depositional units defined by the percent wood content and type of wood: native material (<5% wood), intermediate (<30% wood), sawdust (<30% wood), woodchips (<30% wood), and poorly sorted sands with silt (SM-SP) (0% wood). The contaminants include arsenic, 2,4-dimethylphenol, and total organic carbon. Three-dimensional modeling software, RockWorks, interpolated the discrete borehole data of sediment and contaminants assuming horizontal continuity between sampling locations. The sediment distribution was calculated within concentration ranges for each contaminant of concern. The lowest detection limits, the screening levels, and the cleanup levels defined these ranges. Total organic carbon served as a proxy to estimate the quantity of wood waste in the sediment. As a known byproduct of wood decomposition, 2,4-dimethylphenol was expected to be more prevalent in the depositional units with more wood waste. Finally, arsenic was a proxy for other contaminants to determine if contaminants at Mill-A are dominant in sediments with high percentages of wood waste. The volumetric distribution established that high levels of total organic carbon are present in the sediment with higher percentages of wood waste. This correlation was stronger in the decomposing sawdust-rich sediment than the woodchip-rich sediment. The 2,4-dimethylphenol concentrations above cleanup standards were dominant in the sawdust-rich, intermediate and native sediments. Concentrations of 2,4-dimethylphenol below cleanup levels characterized the native sediment. The distribution of arsenic showed no statistically significant correlation to wood content in sediment. These results do not support the hypothesis of contaminant-rich wood waste, as many of the high concentrations of contaminants were not in the wood-rich sediments. This suggests that the contaminants are more distributed among all depositional units at Mill-A rather than focused within sediments with a high percent of wood waste. Understanding the distribution of potentially toxic compounds with wood waste is important for restoring the Puget Sound waterways to a more habitable environment. Future studies should include new data to validate these results and to limit the uncertainty of the extent of contaminants. Future studies may also find motive in looking for a correlation between contaminants and grain size based on previous studies linking these characteristics. These investigations will benefit the current cleanup effort as well as future cleanup efforts at similarly contaminated waterways.

Environmental controls on zooplankton composition and distribution in Muchalat Inlet, Nootka Sound, B.C., Canada

Senior thesis written for Oceanography 445

Modeling collisional processes in plasmas using discontinuous numerical methods

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

Modeling of Distributed Generation under Next Generation Interconnection Requirements

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

Stochastic Modeling of Reversible Biochemical Reaction-Diffusion Systems and High-Resolution Shock-Capturing Methods for Fluid Interfaces

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

The spectral albedo of sea ice and salt crusts on the tropical ocean of Snowball Earth: II. Optical modeling

During the Snowball Earth events of the Neoproterozoic, tropical regions of the ocean could have developed a precipitated salt lag deposit left behind by sublimating sea ice. The major salt would have been hydrohalite, NaCl•2H2O. The crystals in such a deposit can be small and highly scattering, resulting in an allwave albedo similar to that of snow. The snow-free sea ice from which such a crust could develop has a lower albedo, around 0.5, so the development of a crust would substantially increase the albedo of tropical regions on Snowball Earth. Hydrohalite crystals are much less absorptive than ice in the near- infrared part of the solar spectrum, so their presence at the surface would increase the overall albedo as well as altering its spectral distribution. In this paper, we use laboratory measurements of the spectral albedo of a hydrohalite lag deposit, in combination with a radiative transfer model, to infer the inherent optical properties of hydrohalite as functions of wavelength. Using this result, we model mixtures of hydrohalite and ice representing both artificially created surfaces in the laboratory and surfaces relevant to Snowball Earth. The model is tested against sequences of laboratory measurements taken during the formation and the dissolution of a lag deposit of hydrohalite. We present a parameterization for the broadband albedo of cold, sublimating sea ice as it forms and evolves a hydrohalite crust, for use in climate models of Snowball Earth.