An Integrated Assessment of the Red Sea Urchin (Strongylocentrotus franciscanus) Fisheries in Southern California and Washington State: Addressing Sustainability and Vulnerability in the face of Climate Change

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

Cultivating Care: Understanding Intimate Partner Violence Experiences of Undocumented Latinas in Washington State

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

The generation and propagation of internal waves in Nootka Sound

senior thesis written for Oceanography 445

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.

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.

The Ledgewood-Bonair Landslide toe and where did it go? A case study in littoral sediment budgets from a deep-seated landslide in Island County, Washington

On the morning of March 27th, 2013, a small portion of a much larger landslide complex failed on the western shoreline of central Whidbey Island, Island County, Washington. This landslide, known as the Ledgewood-Bonair Landslide (LB Landslide), mobilized as much as 150,000 cubic meters of unconsolidated glacial sediment onto the coastline of the Puget Sound (Slaughter et al., 2013, Geotechnical Engineering Services, 2013). This study aims to determine how sediment from the Ledgewood-Bonair Landslide has acted on the adjacent beaches 400 meters to the north and south, and specifically to evaluate the volume of sediment contributed by the slide to adjacent beaches, how persistent bluff-derived accretion has been on adjacent beaches, and how intertidal grain sizes changed as a result of the bluff-derived sediment, LiDAR imagery from 2013 and 2014 were differenced and compared to beach profile data and grain size photography. Volume change results indicate that of the 41,850 cubic meters of sediment eroded at the toe of the landslide, 8.9 percent was redeposited on adjacent beaches within 1 year of the landslide. Of this 8.9 percent, 6.3 percent ended up on the north beach and 2.6 percent ended up on the south beach. Because the landslide deposit was primarily sands, silts, and clays, it is reasonable to assume that the remaining 91.1 percent of the sediment eroded from the landslide toe was carried out into the waters of the Puget Sound. Over the course of the two-year study, measurable accretion is apparent up to 150 meters north and 100 meters south of the landslide complex. Profile data also suggests that the most significant elevation changes occurred within the first two and half months since the landslides occurrence. The dominant surficial grain size of the beach soon after the landslide was coarse-sand; in the years following the landslide, 150 meters north of the toe the beach sediment became finer while 100 meters south of the toe the beach sediment became coarser. Overall, the LB Landslide has affected beach profile and grain size only locally, within 150 meters of the landslide toe.

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.

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.

Fetch-limited wave growth in Nootka Sound

Senior thesis written for Oceanography 445

Tidal Energy Resource Characterization in Alberni Inlet, B.C.

Senior thesis written for Oceanography 444

A Meteorological and Snow observational data set from Snoqualmie Pass (921 m), Washington Cascades, U.S.

We introduce a quality controlled observational atmospheric, snow, and soil data set from Snoqualmie Pass, Washington, U.S.A., to enable testing of hydrometeorological and snow process representations within a rain-snow transitional climate where existing observations are sparse and limited. Continuous meteorological forcing (including air temperature, total precipitation, wind speed, specific humidity, air pressure, short- and longwave irradiance) are provided at hourly intervals for a 24-year historical period (water years 1989-2012) and at half-hourly intervals for a more-recent period (water years 2013-2015), separated based on the availability of observations. Additional observations include 40-years of snow board new snow accumulation, multiple measurements of total snow depth, and manual snow pits, while more recent years include sub-daily surface temperature, snowpack drainage, soil moisture and temperature profiles, and eddy co-variance derived turbulent heat flux. This data set is ideal for testing hypotheses about energy balance, soil and snow processes in the rain-snow transition zone. Plots of live data can be found here: http://depts.washington.edu/mtnhydr/cgi/plot.cgi

Analyzing drifter designs and data gathering methods in complex estuarine environments

Senior thesis written for Oceanography 445

Growing biofuel outdoors in Seattle: effects of natural light level and temperature on marine microalgal culture density and lipid content in contained natural exposure systems

Senior thesis written for Oceanography 445

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.