USGS Science⌕ Search

Geology topics

Theresa D. Olsen

Publications and source records attributed to Theresa D. Olsen.

18 recordsLinked to original sources
USGS2018

Simulation of groundwater storage changes in the Quincy Basin, Washington

The Miocene Columbia River Basalt Group and younger sedimentary deposits of lacustrine, fluvial, eolian, and cataclysmic-flood origins compose the aquifer system of the Quincy Basin in eastern Washington. Irrigation return flow and canal leakage from the Columbia Basin Project have caused groundwater levels to rise substantially in some areas. Water resource managers are considering extraction of additional stored groundwater to supply increasing demand. To help address these concerns, the transient groundwater model of the Quincy Basin documented in this report was developed to quantify the changes in groundwater flow and storage. The model based on the U.S. Geological Survey modular three-dimensional finite-difference numerical code MODFLOW uses a 1-kilometer finite-difference grid and is constrained by logs from 698 wells in the study area. Five model layers represent two sedimentary hydrogeologic units and underlying basalt formations. Head-dependent flux boundaries represent the Columbia River and other streams, lakes and reservoirs, underflow to and (or) from adjacent areas, and discharge to agricultural drains and springs. Specified flux boundaries represent recharge from precipitation and anthropogenic sources, including irrigation return flow and leakage from water-distribution canals and discharge through groundwater withdrawal wells. Transient conditions were simulated from 1920 to 2013 using annual stress periods. The model was calibrated with the parameter-estimation code PEST to a total of 4,064 water levels measured in 710 wells. Increased recharge since predevelopment resulted in an 11.5 million acre-feet increase in storage in the Quincy Groundwater Management Subarea of the Quincy Basin. Four groundwater-management scenarios were formulated with input from project stakeholders and were simulated using the calibrated model to provide representative examples of how the model could be used to evaluate the effect on groundwater levels as a result of potential changes in recharge, groundwater withdrawals, or increased flow in Crab Creek. Decreased recharge and increased groundwater withdrawals both resulted in declines in groundwater levels over 2013 conditions, whereas increasing the flow in Crab Creek resulted in increased groundwater levels over 2013 conditions.

Washington↗
USGS2017

Groundwater data collection for the Quinault Indian Nation, Grays Harbor and Jefferson Counties, Washington

Groundwater data were collected on the Quinault Indian Reservation to provide the Quinualt Indian Nation (QIN) with basic knowledge of the existing wells and springs on the reservation, and to establish a water-level network to be monitored by QIN to begin building a long-term groundwater dataset. The 327 mi 2 Quinault Indian Reservation is located within the heavily forested Queets-Quinault watershed along the west-central coast of Washington and includes the coastal communities of Taholah and Queets, and the inland community of Amanda Park. Groundwater data were collected or compiled for 87 sites—82 wells and 5 springs. In October 2016, a field inventory was done to locate the sites and acquire site data. Groundwater levels were measured in 15 of the field-inventoried wells and 3 of those wells were observed as flowing (artesian). A monthly groundwater‑level monitoring network of 13 wells was established by the U.S. Geological Survey in March 2017, and the network was transferred to QIN in June 2017 for continued measurements. Several data needs were identified that would provide a more complete understanding of the groundwater system of the Quinault Indian Reservation. The collection of monthly water-level data for multiple years is an important first step in understanding seasonal and long term changes in water levels. Additionally, the collection of baseline groundwater chemistry and quality data across the reservation would help with future efforts to monitor existing and potentially changing groundwater quality conditions. Development of a water budget of the Queets-Quinault Watershed and the reservation within that area would provide water users with a better understanding of this important resource and provide needed information about the competing demands on local water sources.

Washington↗
USGS2016

Numerical simulation of the groundwater-flow system of the Kitsap Peninsula, west-central Washington

A groundwater-flow model was developed to improve understanding of water resources on the Kitsap Peninsula. The Kitsap Peninsula is in the Puget Sound lowland of west-central Washington, is bounded by Puget Sound on the east and by Hood Canal on the west, and covers an area of about 575 square miles. The peninsula encompasses all of Kitsap County, Mason County north of Hood Canal, and part of Pierce County west of Puget Sound. The peninsula is surrounded by saltwater, and the hydrologic setting is similar to that of an island. The study area is underlain by a thick sequence of unconsolidated glacial and interglacial deposits that overlie sedimentary and volcanic bedrock units that crop out in the central part of the study area. Twelve hydrogeologic units consisting of aquifers, confining units, and an underlying bedrock unit form the basis of the groundwater-flow model. Groundwater flow on the Kitsap Peninsula was simulated using the groundwater-flow model, MODFLOW‑NWT. The finite difference model grid comprises 536 rows, 362 columns, and 14 layers. Each model cell has a horizontal dimension of 500 by 500 feet, and the model contains a total of 1,227,772 active cells. Groundwater flow was simulated for transient conditions. Transient conditions were simulated for January 1985–December 2012 using annual stress periods for 1985–2004 and monthly stress periods for 2005–2012. During model calibration, variables were adjusted within probable ranges to minimize differences between measured and simulated groundwater levels and stream baseflows. As calibrated to transient conditions, the model has a standard deviation for heads and flows of 47.04 feet and 2.46 cubic feet per second, respectively. Simulated inflow to the model area for the 2005–2012 period from precipitation and secondary recharge was 585,323 acre-feet per year (acre-ft/yr) (93 percent of total simulated inflow ignoring changes in storage), and simulated inflow from stream and lake leakage was 43,905 acre-ft/yr (7 percent of total simulated inflow). Simulated outflow from the model primarily was through discharge to streams, lakes, springs, seeps, and Puget Sound (594,595 acre-ft/yr; 95 percent of total simulated outflow excluding changes in storage) and through withdrawals from wells (30,761 acre-ft/yr; 5 percent of total simulated outflow excluding changes in storage). Six scenarios were formulated with input from project stakeholders and were simulated using the calibrated model to provide representative examples of how the model could be used to evaluate the effects on water levels and stream baseflows of potential changes in groundwater withdrawals, in consumptive use, and in recharge. These included simulations of a steady-state system, no-pumping and return flows, 15-percent increase in current withdrawals in all wells, 80-percent decrease in outdoor water to simulate effects of conservation efforts, 15-percent decrease in recharge from precipitation to simulate a drought, and particle tracking to determine flow paths. Changes in water-level altitudes and baseflow amounts vary depending on the stress applied to the system in these various scenarios. Reducing recharge by 15 percent between 2005 and 2012 had the largest effect, with water-level altitudes declining throughout the model domain and baseflow amounts decreasing by as much as 18 percent compared to baseline conditions. Changes in pumping volumes had a smaller effect on the model. Removing all pumping and resulting return flows caused increased water-level altitudes in many areas and increased baseflow amounts of between 1 and 3 percent.

Washington↗
USGS2016

Simulation of groundwater storage changes in the eastern Pasco Basin, Washington

The Miocene Columbia River Basalt Group and younger sedimentary deposits of lacustrine, fluvial, eolian, and cataclysmic-flood origins compose the aquifer system of the Pasco Basin in eastern Washington. Irrigation return flow and canal leakage from the Columbia Basin Project have caused groundwater levels to rise substantially in some areas, contributing to landslides along the Columbia River. Water resource managers are considering extraction of additional stored groundwater to supply increasing demand and possibly mitigate problems caused by the increased water levels. To help address these concerns, the transient groundwater model of the Pasco Basin documented in this report was developed to quantify the changes in groundwater flow and storage. The MODFLOW model uses a 1-kilometer finite-difference grid and is constrained by logs and water levels from 846 wells in the study area. Eight model layers represent five sedimentary hydrogeologic units and underlying basalt formations. Head‑dependent flux boundaries represent the Columbia and Snake Rivers to the west and south, respectively, underflow to and (or) from adjacent areas to the northeast, and discharge to agricultural drains, springs, and groundwater withdrawal wells. Specified flux boundaries represent recharge from infiltrated precipitation and anthropogenic sources, including irrigation return flow and leakage from water-distribution canals. The model was calibrated with the parameter ‑ estimation code PEST++ to groundwater levels measured from 1907 through 2013 and measured discharge to springs and estimated discharge to agricultural drains. Increased recharge since pre-development resulted in a 6.8 million acre-feet increase in storage in the 508-14 administrative area of the Pasco Basin. Four groundwater-management scenarios simulate the 7-year drawdown resulting from withdrawals in different locations. Withdrawals of 2 million gallons per day (Mgal/d) from a hypothetical well field in the upper Ringold Formation along the Columbia River could generate 30–70 feet of drawdown, which may reduce landslide susceptibility along the White Bluffs. Drawdowns resulting from a 1 Mgal/d withdrawal from wells screened in either Pasco gravels, upper Ringold Formation, or both Ringold Formation and underlying basalt are simulated in the other three scenarios, and differ because of the contrasting hydraulic conductivities within the screened intervals.

Washington↗
USGS2014

Hydrogeologic framework, groundwater movement, and water budget of the Kitsap Peninsula, west-central Washington

This report presents information used to characterize the groundwater-flow system on the Kitsap Peninsula, and includes descriptions of the geology and hydrogeologic framework, groundwater recharge and discharge, groundwater levels and flow directions, seasonal groundwater-level fluctuations, interactions between aquifers and the surface‑water system, and a water budget. The Kitsap Peninsula is in the Puget Sound lowland of west-central Washington, is bounded by Puget Sound on the east and by Hood Canal on the west, and covers an area of about 575 square miles. The peninsula encompasses all of Kitsap County, the part of Mason County north of Hood Canal, and part of Pierce County west of Puget Sound. The peninsula is surrounded by saltwater and the hydrologic setting is similar to that of an island. The study area is underlain by a thick sequence of unconsolidated glacial and interglacial deposits that overlie sedimentary and volcanic bedrock units that crop out in the central part of the study area. Geologic units were grouped into 12 hydrogeologic units consisting of aquifers, confining units, and an underlying bedrock unit. A surficial hydrogeologic unit map was developed and used with well information from 2,116 drillers’ logs to construct 6 hydrogeologic sections and unit extent and thickness maps. Unconsolidated aquifers typically consist of moderately to well-sorted alluvial and glacial outwash deposits of sand, gravel, and cobbles, with minor lenses of silt and clay. These units often are discontinuous or isolated bodies and are of highly variable thickness. Unconfined conditions occur in areas where aquifer units are at land surface; however, much of the study area is mantled by glacial till, and confined aquifer conditions are common. Groundwater in the unconsolidated aquifers generally flows radially off the peninsula in the direction of Puget Sound and Hood Canal. These generalized flow patterns likely are complicated by the presence of low-permeability confining units that separate discontinuous bodies of aquifer material and act as local groundwater-flow barriers. Groundwater-level fluctuations observed during the monitoring period (2011–12) in wells completed in unconsolidated hydrogeologic units indicated seasonal variations ranging from 1 to about 20 feet. The largest fluctuation of 33 feet occurred in a well that was completed in the bedrock unit. Streamgage discharge measurements made during 2012 indicate that groundwater discharge to creeks in the area ranged from about 0.41 to 33.3 cubic feet per second. During 2012, which was an above-average year of precipitation, the groundwater system received an average of about 664,610 acre-feet of recharge from precipitation and 22,122 acre-feet of recharge from return flows. Most of this annual recharge (66 percent) discharged to streams, and only about 4 percent was withdrawn from wells. The remaining groundwater recharge (30 percent) left the groundwater system as discharge to Hood Canal and Puget Sound.

Washington↗
USGS2013

Hydrogeology of the Little Spokane River basin, Spokane, Stevens, and Pend Oreille Counties, Washington

A study of the hydrogeologic framework of the Little Spokane River Basin was conducted to identify and describe the principal hydrogeologic units in the study area, their hydraulic characteristics, and general directions of groundwater movement. The Little Spokane River Basin includes an area of 679 square miles in northeastern Washington State covering parts of Spokane, Stevens, and Pend Oreille Counties. The groundwater system consists of unconsolidated sedimentary deposits and isolated, remnant basalt layers overlying crystalline bedrock. In 1976, a water resources program for the Little Spokane River was adopted into rule by the State of Washington, setting instream flows for the river and closing its tributaries to further uses. Spokane County representatives are concerned about the effects that additional groundwater development within the basin might have on the Little Spokane River and on existing groundwater resources. Information provided by this study will be used in future investigations to evaluate the effects of potential increases in groundwater withdrawals on groundwater and surface-water resources in the basin. The hydrogeologic framework consists of eight hydrogeologic units: the Upper aquifer, Upper confining unit, Lower aquifers, Lower confining unit, Wanapum basalt unit, Latah unit, Grande Ronde basalt unit, and Bedrock. The Upper aquifer is composed mostly of sand and gravel and varies in thickness from 4 to 360 ft, with an average thickness of 70 ft. The aquifer is generally finer grained in areas farther from main outwash channels. The estimated horizontal hydraulic conductivity ranges from 4.4 to 410,000 feet per day (ft/d), with a median hydraulic conductivity of 900 ft/d. The Upper confining unit is a low-permeability unit consisting mostly of silt and clay, and varies in thickness from 5 to 400 ft, with an average thickness of 100 ft. The estimated horizontal hydraulic conductivity ranges from 0.5 to 5,600 ft/d, with a median hydraulic conductivity of 8.2 ft/d. The Lower aquifers unit consists of localized confined aquifers or lenses consisting mostly of sand that occur at depth in various places in the basin; thickness of the unit ranges from 8 to 150 ft, with an average thickness of 50 ft. The Lower confining unit is a low-permeability unit consisting mostly of silt and clay; thickness of the unit ranges from 35 to 310 ft, with an average thickness of 130 ft. The Wanapum basalt unit includes the Wanapum Basalt of the Columbia River Basalt Group, thin sedimentary interbeds, and, in some places, overlying loess. The unit occurs as isolated remnants on the basalt bluffs in the study area and ranges in thickness from 7 to 140 ft, with an average thickness of 60 ft. The Latah unit is a mostly low-permeability unit consisting of silt, clay, and sand that underlies and is interbedded with the basalt units. The Latah unit ranges in thickness from 10 to 700 ft, with an average thickness of 250 ft. The estimated horizontal hydraulic conductivity ranges from 0.19 to 15 ft/d, with a median hydraulic conductivity of 0.56 ft/d. The Grande Ronde unit includes the Grande Ronde Basalt of the Columbia River Basalt Group and sedimentary interbeds. Unit thickness ranges from 30 to 260 ft, with an average thickness of 140 ft. The estimated horizontal hydraulic conductivity ranges from 0.03 to 13 ft/d, with a median hydraulic conductivity of 2.9 ft/d. The Bedrock unit is the only available source of groundwater where overlying sediments are absent or insufficiently saturated. The estimated horizontal hydraulic conductivity ranges from 0.01 to 5,000 ft/d, with a median hydraulic conductivity of 1.4 ft/d. The altitude of the buried bedrock surface ranges from about 2,200 ft to about 1,200 ft. Groundwater movement in the Little Spokane River Basin mimics the surface-water drainage pattern of the basin, moving from the topographically high tributary-basin areas toward the topographically lower valley floors. Water-level altitudes range from more than 2,700 ft to about 1,500 ft near the basin’s outlet.

Washington↗
USGS2012

Changes in sediment volume in Alder Lake, Nisqually River Basin, Washington, 1945-2011

The Nisqually River drains the southwest slopes of Mount Rainier, a glaciated stratovolcano in the Cascade Range of western Washington. The Nisqually River was impounded behind Alder Dam when the dam was completed in 1945 and formed Alder Lake. This report quantifies the volume of sediment deposited by the Nisqually and Little Nisqually Rivers in their respective deltas in Alder Lake since 1945. Four digital elevation surfaces were generated from historical contour maps from 1945, 1956, and 1985, and a bathymetric survey from 2011. These surfaces were used to compute changes in sediment volume since 1945. Estimates of the volume of sediment deposited in Alder Lake between 1945 and 2011 were focused in three areas: (1) the Nisqually River delta, (2) the main body of Alder Lake, along a 40-meter wide corridor of the pre-dam Nisqually River, and (3) the Little Nisqually River delta. In each of these areas the net deposition over the 66-year period was 42,000,000 ± 4,000,000 cubic meters (m 3 ), 2,000,000 ± 600,000 m 3 , and 310,000 ± 110,000 m 3 , respectively. These volumes correspond to annual rates of accumulation of 630,000 ± 60,000 m 3 /yr, 33,000 ± 9,000 m 3 /yr, and 4,700 ± 1,600 m 3 /yr, respectively. The annual sediment yield of the Nisqually (1,100 ± 100 cubic meters per year per square kilometer [(m 3 /yr)/km 2 ]) and Little Nisqually River basins [70 ± 24 (m 3 /yr)/km 2 ] provides insight into the yield of two basins with different land cover and geomorphic processes. These estimates suggest that a basin draining a glaciated stratovolcano yields approximately 15 times more sediment than a basin draining forested uplands in the Cascade Range. Given the cumulative net change in sediment volume in the Nisqually River delta in Alder Lake, the total capacity of Alder Lake since 1945 decreased about 3 percent by 1956, 8 percent by 1985, and 15 percent by 2011.

Washington↗
USGS2012

Geomorphic analysis of the river response to sedimentation downstream of Mount Rainier, Washington

A study of the geomorphology of rivers draining Mount Rainier, Washington, was completed to identify sources of sediment to the river network; to identify important processes in the sediment delivery system; to assess current sediment loads in rivers draining Mount Rainier; to evaluate if there were trends in streamflow or sediment load since the early 20th century; and to assess how rates of sedimentation might continue into the future using published climate-change scenarios. Rivers draining Mount Rainier carry heavy sediment loads sourced primarily from the volcano that cause acute aggradation in deposition reaches as far away as the Puget Lowland. Calculated yields ranged from 2,000 tonnes per square kilometer per year [(tonnes/km 2 )/yr] on the upper Nisqually River to 350 (tonnes/km 2 )/yr on the lower Puyallup River, notably larger than sediment yields of 50–200 (tonnes/km2)/yr typical for other Cascade Range rivers. These rivers can be assumed to be in a general state of sediment surplus. As a result, future aggradation rates will be largely influenced by the underlying hydrology carrying sediment downstream. The active-channel width of rivers directly draining Mount Rainier in 2009, used as a proxy for sediment released from Mount Rainier, changed little between 1965 and 1994 reflecting a climatic period that was relatively quiet hydrogeomorphically. From 1994 to 2009, a marked increase in geomorphic disturbance caused the active channels in many river reaches to widen. Comparing active-channel widths of glacier-draining rivers in 2009 to the distance of glacier retreat between 1913 and 1994 showed no correlation, suggesting that geomorphic disturbance in river reaches directly downstream of glaciers is not strongly governed by the degree of glacial retreat. In contrast, there was a correlation between active-channel width and the percentage of superglacier debris mantling the glacier, as measured in 1971. A conceptual model of sediment delivery processes from the mountain indicates that rockfalls, glaciers, debris flows, and main-stem flooding act sequentially to deliver sediment from Mount Rainier to river reaches in the Puget Lowland over decadal time scales. Greater-than-normal runoff was associated with cool phases of the Pacific Decadal Oscillation. Streamflow-gaging station data from four unregulated rivers directly draining Mount Rainier indicated no statistically significant trends of increasing peak flows over the course of the 20th century. The total sediment load of the upper Nisqually River from 1945 to 2011 was determined to be 1,200,000±180,000 tonnes/yr. The suspended-sediment load in the lower Puyallup River at Puyallup, Washington, was 860,000±300,000 tonnes/yr between 1978 and 1994, but the long-term load for the Puyallup River likely is about 1,000,000±400,000 tonnes/yr. Using a coarse-resolution bedload transport relation, the long-term average bedload was estimated to be about 30,000 tonnes/yr in the lower White River near Auburn, Washington, which was four times greater than bedload in the Puyallup River and an order of magnitude greater than bedload in the Carbon River. Analyses indicate a general increase in the sediment loads in Mount Rainier rivers in the 1990s and 2000s relative to the time period from the 1960s to 1980s. Data are insufficient, however, to determine definitively if post-1990 increases in sediment production and transport from Mount Rainier represent a statistically significant increase relative to sediment-load values typical from Mount Rainier during the entire 20th century. One-dimensional river-hydraulic and sediment-transport models simulated the entrainment, transport, attrition, and deposition of bed material. Simulations showed that bed-material loads were largest for the Nisqually River and smallest for the Carbon River. The models were used to simulate how increases in sediment supply to rivers transport through the river systems and affect lowland reaches. For each simulation, the input sediment pulse evolved through a combination of translation, dispersion, and attrition as it moved downstream. The characteristic transport times for the median sediment-size pulse to arrive downstream for the Nisqually, Carbon, Puyallup, and White Rivers were approximately 70, 300, 80, and 60 years, respectively.

Washington↗
USGS2011

Conceptual model and numerical simulation of the groundwater-flow system of Bainbridge Island, Washington

Groundwater is the sole source of drinking water for the population of Bainbridge Island. Increased use of groundwater supplies on Bainbridge Island as the population has grown over time has created concern about the quantity of water available and whether saltwater intrusion will occur as groundwater usage increases. A groundwater-flow model was developed to aid in the understanding of the groundwater system and the effects of groundwater development alternatives on the water resources of Bainbridge Island. Bainbridge Island is underlain by unconsolidated deposits of glacial and nonglacial origin. The surficial geologic units and the deposits at depth were differentiated into aquifers and confining units on the basis of areal extent and general water-bearing characteristics. Eleven principal hydrogeologic units are recognized in the study area and form the basis of the groundwater-flow model. A transient variable-density groundwater-flow model of Bainbridge Island and the surrounding area was developed to simulate current (2008) groundwater conditions. The model was calibrated to water levels measured during 2007 and 2008 using parameter estimation (PEST) to minimize the weighted differences or residuals between simulated and measured hydraulic head. The calibrated model was used to make some general observations of the groundwater system in 2008. Total flow through the groundwater system was about 31,000 acre-ft/ yr. The recharge to the groundwater system was from precipitation and septic-system returns. Groundwater flow to Bainbridge Island accounted for about 1,000 acre-ft/ yr or slightly more than 5 percent of the recharge amounts. Groundwater discharge was predominately to streams, lakes, springs, and seepage faces (16,000 acre-ft/yr) and directly to marine waters (10,000 acre-ft/yr). Total groundwater withdrawals in 2008 were slightly more than 6 percent (2,000 acre-ft/yr) of the total flow. The calibrated model was used to simulate predevelopment conditions, during which no groundwater pumping or secondary recharge occurred and currently developed land was covered by conifer forests. Simulated water levels in the uppermost aquifer generally were slightly higher at the end of 2008 than under predevelopment conditions, likely due to increased recharge from septic returns and reduced evapotranspiration losses due to conversion of land cover from forests to current conditions. Simulated changes in water levels for the extensively used sea-level aquifer were variable, although areas with declines between zero and 10 feet were common and generally can be traced to withdrawals from public-supply drinking wells. Simulated water-level declines in the deep (Fletcher Bay) aquifer between predevelopment and 2008 conditions ranged from about 10 feet in the northeast to about 25 feet on the western edge of the Island. These declines are related to groundwater withdrawals for public-supply purposes. The calibrated model also was used to simulate the possible effects of increased groundwater pumping and changes to recharge due to changes in land use and climactic conditions between 2008 and 2035 under minimal, expected, and maximum impact conditions. Drawdowns generally were small for most of the Island (less than 10 ft) for the minimal and expected impact scenarios, and were larger for the maximum impact scenario. No saltwater intrusion was evident in any scenario by the year 2035. The direction of flow in the deep Fletcher Bay aquifer was simulated to reverse direction from its predevelopment west to east direction to an east to west direction under the maximum impact scenario.

Scientific Investigations Report↗
USGS2011

Hydrogeologic framework, groundwater movement, and water budget in the Chimacum Creek basin and vicinity, Jefferson County, Washington

This report presents information used to characterize the groundwater flow system in the Chimacum Creek basin. It includes descriptions of the geology and hydrogeologic framework; groundwater recharge and discharge; groundwater levels and flow directions; seasonal fluctuations in groundwater level; interactions between aquifers and the surface-water system; and a groundwater budget. The study area covers 124 square miles in northeastern Jefferson County, Washington, and includes the Chimacum Creek basin, which drains an area of about 37 square miles. The area is underlain by a north-thickening sequence of unconsolidated glacial and interglacial deposits that overlie sedimentary and igneous bedrock units that crop out along the margins and western interior of the study area. Six hydrogeologic units consisting of unconsolidated aquifers and confining units, along with an underlying bedrock unit, were identified. A surficial hydrogeologic map was developed and used with well information from 187 drillers' logs to construct 4 hydrogeologic sections, and maps showing the extent and thickness of the units. Natural recharge was estimated using precipitation-recharge relation regression equations developed for western Washington, and estimates were calculated for return flow from data on domestic indoor and outdoor use and irrigated agriculture. Results from synoptic streamflow measurements and water table elevations determined from monthly measurements at monitoring wells are presented and compared with those from a study conducted during 2002-03. A water budget was calculated comprising long-term average recharge, domestic public-supply withdrawals and return flow, self-supplied domestic withdrawals and return flow, and irrigated agricultural withdrawals and return flow.

Washington↗
USGS2010

Hydrogeologic framework, groundwater and surface-water systems, land use, pumpage, and water budget of the Chamokane Creek basin, Stevens County, Washington

A study of the water resources of the unconsolidated groundwater system of the Chamokane Creek basin was conducted to determine the hydrogeologic framework, interactions of shallow and deep parts of the groundwater system with each other and the surface-water system, changes in land use and land cover, and water-use estimates. Chamokane Creek basin is a 179 mi2 area that borders and partially overlaps the Spokane Indian Reservation in southern Stevens County in northeastern Washington State. Aquifers within the Chamokane Creek basin are part of a sequence of glaciofluvial and glaciolacustrine sediment that may reach total thicknesses of about 600 ft. In 1979, most of the water rights in the Chamokane Creek basin were adjudicated by the United States District Court requiring regulation in favor of the Spokane Tribe of Indians' senior water right. The Spokane Tribe, the State of Washington, and the United States are concerned about the effects of additional groundwater development within the basin on Chamokane Creek. Information provided by this study will be used to evaluate the effects of potential increases in groundwater withdrawals on groundwater and surface-water resources within the basin. The hydrogeologic framework consists of six hydrogeologic units: The Upper outwash aquifer, the Landslide Unit, the Valley Confining Unit, the Lower Aquifer, the Basalt Unit, and the Bedrock Unit. The Upper outwash aquifer occurs along the valley floors of the study area and consists of sand, gravel, cobbles, boulders, with minor silt and (or) clay interbeds in places. The Lower aquifer is a confined aquifer consisting of sand and gravel that occurs at depth below the Valley confining unit. Median horizontal hydraulic conductivity values for the Upper outwash aquifer, Valley confining unit, Lower aquifer, and Basalt unit were estimated to be 540, 10, 19, and 3.7 ft/d, respectively. Many low-flow stream discharge measurements at sites on Chamokane Creek and its tributaries were at or near zero flow. The most notable exception is where Chamokane Creek is supported by discharge of large springs from the Upper outwash aquifer in the southern part of the basin. Most high-flow measurements indicated gains in streamflow (groundwater discharging to the stream). Large streamflow losses, however, were recorded near the north end of Walkers Prairie where streamflow directly recharges the Upper outwash aquifer. The similarity in seasonal water-level fluctuations in the Upper outwash aquifer and the Lower aquifer indicate that these systems may be fairly well connected. Land use and land cover change analysis indicates that Chamokane Creek basin has been dominated by forests with some pasture and agricultural lands with sparse residential development from the 1980s to present. Loss in forest cover represents the largest change in land cover in the basin between 1987 and 2009. This appears to be mostly due to forestry activities, especially in the northern part of the basin. Since 1987, more than 18,000 acres of evergreen forest have been logged and are at various stages of regrowth. Estimated average annual total groundwater pumpage in the basin increased from 224 million gallons per year (Mgal/yr) in 1980 to 1,330 Mgal/yr in 2007. The largest withdrawals during 2007 were to supply two fish hatcheries, with a combined total annual pumpage of about 1,150 Mgal. Annual groundwater pumpage values from 1980 through 2007 for the study area ranged from 21.1 to 28.9 Mgal/yr for domestic wells and 0.38 to 23.7 Mgal/yr for public supply. An approximate water budget for a typical year in the Chamokane Creek basin indicates that 19.6 in. of precipitation are balanced by 4.7 in. of streamflow discharge from the basin, and 14.9 in. of evapotranspiration.

Washington↗
USGS2009

Geologic Setting and Hydrogeologic Units of the Columbia Plateau Regional Aquifer System, Washington, Oregon, and Idaho

The Columbia Plateau Regional Aquifer System (CPRAS) covers approximately 44,000 square miles of northeastern Oregon, southeastern Washington, and western Idaho. The area supports a $6 billion per year agricultural industry, leading the Nation in production of apples and nine other commodities (State of Washington Office of Financial Management, 2007; U.S. Department of Agriculture, 2007). Groundwater availability in the aquifers of the area is a critical water-resource management issue because the water demand for agriculture, economic development, and ecological needs is high. The primary aquifers of the CPRAS are basalts of the Columbia River Basalt Group (CRBG) and overlying basin-fill sediments. Water-resources issues that have implications for future groundwater availability in the region include (1) widespread water-level declines associated with development of groundwater resources for irrigation and other uses, (2) reduction in base flow to rivers and associated effects on temperature and water quality, and (3) current and anticipated effects of global climate change on recharge, base flow, and ultimately, groundwater availability. As part of a National Groundwater Resources Program, the U.S. Geological Survey began a study of the CPRAS in 2007 with the broad goals of (1) characterizing the hydrologic status of the system, (2) identifying trends in groundwater storage and use, and (3) quantifying groundwater availability. The study approach includes documenting changes in the status of the system, quantifying the hydrologic budget for the system, updating the regional hydrogeologic framework, and developing a groundwater-flow simulation model for the system. The simulation model will be used to evaluate and test the conceptual model of the system and later to evaluate groundwater availability under alternative development and climate scenarios. The objectives of this study were to update the hydrogeologic framework for the CPRAS using the available geologic mapping and well information and to develop a digital, three-dimensional hydrogeologic model that could be used as the basis of a groundwater-flow model. This report describes the principal geologic and hydrogeologic units of the CPRAS and geologic map and well data that were compiled as part of the study. The report also describes simplified regional hydrogeologic sections and unit extent maps that were used to conceptualize the framework prior to development of the digital 3-dimensional framework model.

Idaho, Oregon, Washington↗
USGS2009

Navigability potential of Washington rivers and streams determined with hydraulic geometry and a geographic information system

Using discharge and channel geometry measurements from U.S. Geological Survey streamflow-gaging stations and data from a geographic information system, regression relations were derived to predict river depth, top width, and bottom width as a function of mean annual discharge for rivers in the State of Washington. A new technique also was proposed to determine bottom width in channels, a parameter that has received relatively little attention in the geomorphology literature. These regression equations, when combined with estimates of mean annual discharge available in the National Hydrography Dataset, enabled the prediction of hydraulic geometry for any stream or river in the State of Washington. Predictions of hydraulic geometry can then be compared to thresholds established by the Washington State Department of Natural Resources to determine navigability potential of rivers. Rivers with a mean annual discharge of 1,660 cubic feet per second or greater are 'probably navigable' and rivers with a mean annual discharge of 360 cubic feet per second or less are 'probably not navigable'. Variance in the dataset, however, leads to a relatively wide range of prediction intervals. For example, although the predicted hydraulic depth at a mean annual discharge of 1,660 cubic feet per second is 3.5 feet, 90-percent prediction intervals indicate that the actual hydraulic depth may range from 1.8 to 7.0 feet. This methodology does not determine navigability - a legal concept determined by federal common law - instead, this methodology is a tool for predicting channel depth, top width, and bottom width for rivers and streams in Washington.

Washington↗
USGS2009

Estimating Low-Flow Frequency Statistics and Hydrologic Analysis of Selected Streamflow-Gaging Stations, Nooksack River Basin, Northwestern Washington and Canada

Low-flow frequency statistics were computed at 17 continuous-record streamflow-gaging stations and 8 miscellaneous measurement sites in and near the Nooksack River basin in northwestern Washington and Canada, including the 1, 3, 7, 15, 30, and 60 consecutive-day low flows with recurrence intervals of 2 and 10 years. Using these low-flow statistics, 12 regional regression equations were developed for estimating the same low-flow statistics at ungaged sites in the Nooksack River basin using a weighted-least-squares method. Adjusted R2 (coefficient of determination) values for the equations ranged from 0.79 to 0.93 and the root-mean-squared error (RMSE) expressed as a percentage ranged from 77 to 560 percent. Streamflow records from six gaging stations located in mountain-stream or lowland-stream subbasins of the Nooksack River basin were analyzed to determine if any of the gaging stations could be removed from the network without significant loss of information. Using methods of hydrograph comparison, daily-value correlation, variable space, and flow-duration ratios, and other factors relating to individual subbasins, the six gaging stations were prioritized from most to least important as follows: Skookum Creek (12209490), Anderson Creek (12210900), Warm Creek (12207750), Fishtrap Creek (12212050), Racehorse Creek (12206900), and Clearwater Creek (12207850). The optimum streamflow-gaging station network would contain all gaging stations except Clearwater Creek, and the minimum network would include Skookum Creek and Anderson Creek.

Scientific Investigations Report↗
USGS2009

Hydrogeologic Framework, Groundwater Movement, and Water Budget in Tributary Subbasins and Vicinity, Lower Skagit River Basin, Skagit and Snohomish Counties, Washington

A study to characterize the groundwater-flow system in four tributary subbasins and vicinity of the lower Skagit River basin was conducted by the U.S. Geological Survey to assist Skagit County and the Washington State Department of Ecology in evaluating the effects of potential groundwater withdrawals and consumptive use on tributary streamflows. This report presents information used to characterize the groundwater and surface-water flow system in the subbasins, and includes descriptions of the geology and hydrogeologic framework of the subbasins; groundwater recharge and discharge; groundwater levels and flow directions; seasonal groundwater-level fluctuations; interactions between aquifers and the surface-water system; and a water budget for the subbasins. The study area covers about 247 mi2 along the Skagit River and its tributary subbasins (East Fork Nookachamps Creek, Nookachamps Creek, Carpenter Creek, and Fisher Creek) in southwestern Skagit County and northwestern Snohomish County, Washington. The geology of the area records a complex history of accretion along the continental margin, mountain building, deposition of terrestrial and marine sediments, igneous intrusion, and the repeated advance and retreat of continental glaciers. A simplified surficial geologic map was developed from previous mapping in the area, and geologic units were grouped into nine hydrogeologic units consisting of aquifers and confining units. A surficial hydrogeologic unit map was constructed and, with lithologic information from 296 drillers'logs, was used to produce unit extent and thickness maps and four hydrogeologic sections. Groundwater in unconsolidated aquifers generally flows towards the northwest and west in the direction of the Skagit River and Puget Sound. This generalized flow pattern is likely complicated by the presence of low-permeability confining units that separate discontinuous bodies of aquifer material and act as local groundwater-flow barriers. Groundwater-flow directions in the sedimentary aquifer likely reflect local topographic relief (radial flow from bedrock highs) and more regional westward flow from the mountains to the Puget Sound. The largest groundwater-level fluctuations observed during the monitoring period (October 2006 through September 2008) occurred in wells completed in the sedimentary aquifer, and ranged from about 3 to 27 feet. Water levels in wells completed in unconsolidated hydrogeologic units exhibited seasonal variations ranging from less than 1 to about 10 feet. Synoptic streamflow measurements made in August 2007 and June 2008 indicate a total groundwater discharge to creeks in the tributary subbasin area of about 13.15 and 129.6 cubic feet per second (9,520 and 93,830 acre-feet per year), respectively. Streamflow measurements illustrate a general pattern in which the upper reaches of creeks in the study area tended to gain flow from the groundwater system, and lower creek reaches tended to lose water. Large inflows from tributaries to major creeks in the study area suggest the presence of groundwater discharge from upland areas underlain by bedrock. The groundwater system within the subbasins received an average (September 1, 2006 to August 31, 2008) of about 92,400 acre-feet or about 18 inches of recharge from precipitation a year. Most of this recharge (65 percent) discharges to creeks, and only about 3 percent is withdrawn from wells. The remaining groundwater recharge (32 percent) leaves the subbasin groundwater system as discharge to the Skagit River and Puget Sound.

Scientific Investigations Report↗
USGS2006

Freshwater and saline loads of dissolved inorganic nitrogen to Hood Canal and Lynch Cove, western Washington

Hood Canal is a long (110 kilometers), deep (175 meters) and narrow (2 to 4 kilometers wide) fjord of Puget Sound in western Washington. The stratification of a less dense, fresh upper layer of the water column causes the cold, saltier lower layer of the water column to be isolated from the atmosphere in the late summer and autumn, which limits reaeration of the lower layer. In the upper layer of Hood Canal, the production of organic matter that settles and consumes dissolved oxygen in the lower layer appears to be limited by the load of dissolved inorganic nitrogen (DIN): nitrate, nitrite, and ammonia. Freshwater and saline loads of DIN to Hood Canal were estimated from available historical data. The freshwater load of DIN to the upper layer of Hood Canal, which could be taken up by phytoplankton, came mostly from surface and ground water from subbasins, which accounts for 92 percent of total load of DIN to the upper layer of Hood Canal. Although DIN in rain falling on land surfaces amounts to about one-half of the DIN entering Hood Canal from subbasins, rain falling directly on the surface of marine waters contributed only 4 percent of the load to the upper layer. Point-source discharges and subsurface flow from shallow shoreline septic systems contributed less than 4 percent of the DIN load to the upper layer. DIN in saline water flowing over the sill into Hood Canal from Admiralty Inlet was at least 17 times the total load to the upper layer of Hood Canal. In September and October 2004, field data were collected to estimate DIN loads to Lynch Cove - the most inland marine waters of Hood Canal that routinely contain low dissolved-oxygen waters. Based on measured streamflow and DIN concentrations, surface discharge was estimated to have contributed about one-fourth of DIN loads to the upper layer of Lynch Cove. Ground-water flow from subbasins was estimated to have contributed about one-half of total DIN loads to the upper layer. In autumn 2004, the relative contribution of DIN from shallow shoreline septic systems to the upper layer was higher in Lynch Cove (23 percent) than in the entire Hood Canal. Net transport of DIN into Lynch Cove by marine currents was measured during August and October 2004-a time of high biological productivity. The net transport of lower-layer water into Lynch Cove was significantly diminished relative to the flow entering Hood Canal at its entrance. Even though the net transport of saline water into the lower layer of Lynch Cove was only 119 cubic meters per second, estuarine currents between 33 and 47 m were estimated to have carried more than 35 times the total freshwater load of DIN to the upper layer from surface and ground water, shallow shoreline septic systems, and direct atmospheric rainfall. The subsurface maximums in measured turbidity, chlorophyll a, particulate organic carbon, and particulate organic nitrogen strongly suggest that the upward mixing of nitrate-rich deeper water is a limiting factor in supplying DIN to the upper layer that enhances marine productivity in Lynch Cove. The presence of phosphate in the upper layer in the absence of dissolved inorganic nitrogen also suggests that the biological productivity that leads to low dissolved-oxygen concentrations in the lower layer of Lynch Cove is limited by the supply of nitrogen rather than by phosphate loads. Although the near-shore zones of the shallow parts of Lynch Cove were sampled, a biogeochemical signal from terrestrial nitrogen was not found. Reversals in the normal estuarine circulation suggest that if the relative importance of the DIN load of freshwater terrestrial and atmospheric sources and the DIN load from transport of saline water by the estuarine circulation in controlling dissolved-oxygen concentrations in Lynch Cove is to be better understood, then the physical forces driving Hood Canal circulation must be better defined.

Scientific Investigations Report↗
USGS1999

Hydrogeologic assessment of the Sequim-Dungeness area, Clallam County, Washington

The Sequim-Dungeness area covers 116 square miles (mi 2 ) on the northern part of the Olympic Peninsula in northwestern Washington. The central part of this area (74 mi 2 ) was designated as a primary study area. During the past two decades, the population has rapidly increased, land use has changed from mostly agricultural to residential, and salmon populations in the Dungeness River have appreciably declined. The increasing competition for water combined with a close relation between ground water, the Dungeness River, and an extensive irrigation system has created a need for a better understanding of ground water and the relation between ground water and surface water in the study area. The Sequim-Dungeness area is underlain by as much as 2,000 feet of unconsolidated Quaternary deposits that are mostly of glacial origin. Interpretation of 10 hydrogeologic cross sections and lithologic logs of about 600 wells led to the delineation of three aquifers, two confining beds, and a lower unit of undifferentiated deposits. A bedrock unit at the bottom is considered the base of the ground-water system. Ground water in the study area is recharged from infiltration and percolation of precipitaton, percolation of unconsumed irrigation water, leakage from irrigation ditches, subsurface inflow through the southern study-area boundary, and leakage from streams. Average annual recharge for the study period (December 1995 to September 1997) was estimated to be 17.7 inches (in.) ( 151 cubic feet per second (ft 3 /s)). The distribution of recharge was 8.6 in. (74 ft3fs) from precipitation, 2.7 in. (23 ft 3 /s) from subsurface inflow, 3.1 in. (26 ft 3 /s) from irrigation, and 3.3 in. (28 ft 3 /s) from leakage from the Dungeness River. The 8.6 in. of recharge from precipitation is much higher than would be expected in an average year because average annual precipitation during the study period was about 28 in., which is 1.35 times higher than long-term average annual precipitation. The long-term average annual recharge from precipitation was estimated to be 5.4 in. (48 ft 3 /s). Ground water discharges as subsurface flow to saltwater bodies, flow to streams, flow to springs, and as withdrawals from wells. Subsurface flow to saltwater bodies and flow to springs were not estimated in this study. Estimated average annual discharge was 3.2 inches (in.) (27 ft 3 /s)) to the Dungeness River and 4.6 in. (39 ft 3 /s) to other streams in the study area. Gross withdrawals from wells in 1996 were estimated to be 1.0 in. (8.4 ft 3 /s). There was a small but statistically significant increase in nitrate concentrations in ground water from 1980 to 1996. Median concentrations in the primary study area were 0.37 milligrams per liter (mg/L) in 1980 and 0.46 mg/L in 1996. The areal pattern of elevated nitrate concentrations has not changed appreciably during the past 15 years. Elevated concentrations were found in a large area east of the Dungeness River and at scattered locations west of the Dungeness River. About 543,200 pounds of nitrogen are estimated to enter the ground-water system in the primary study area each year. Four sources account for about 85 percent of the nitrogen; residential fertilizers, septic systems, mineralization of soil organic matter, and agricultural fertilizers .. It appears that the four major sources are approximately equivalent in amounts of nitrogen. Concentrations of nitrate in the shallow aquifer were significantly higher under residential areas than under natural grasslands or forests. Median nitrate concentrations were 1.3 mg/L under residential areas, 0.55 mg/L under agricultural areas, and 0.12 mg/L under natural grasslands or forests.

Washington↗