Geology topics

Kenneth H. Johnson

Publications and source records attributed to Kenneth H. Johnson.

11 recordsLinked to original sources

USGS2016

Groundwater contaminant plume maps and volumes, 100-K and 100-N Areas, Hanford Site, Washington

This study provides an independent estimate of the areal and volumetric extent of groundwater contaminant plumes which are affected by waste disposal in the 100-K and 100-N Areas (study area) along the Columbia River Corridor of the Hanford Site. The Hanford Natural Resource Trustee Council requested that the U.S. Geological Survey perform this interpolation to assess the accuracy of delineations previously conducted by the U.S. Department of Energy and its contractors, in order to assure that the Natural Resource Damage Assessment could rely on these analyses. This study is based on previously existing chemical (or radionuclide) sampling and analysis data downloaded from publicly available Hanford Site Internet sources, geostatistically selected and interpreted as representative of current (from 2009 through part of 2012) but average conditions for groundwater contamination in the study area. The study is limited in scope to five contaminants—hexavalent chromium, tritium, nitrate, strontium-90, and carbon-14, all detected at concentrations greater than regulatory limits in the past. All recent analytical concentrations (or activities) for each contaminant, adjusted for radioactive decay, non-detections, and co-located wells, were converted to log-normal distributions and these transformed values were averaged for each well location. The log-normally linearized well averages were spatially interpolated on a 50 × 50-meter (m) grid extending across the combined 100-N and 100-K Areas study area but limited to avoid unrepresentative extrapolation, using the minimum curvature geostatistical interpolation method provided by SURFER ® data analysis software. Plume extents were interpreted by interpolating the log-normally transformed data, again using SURFER ® , along lines of equal contaminant concentration at an appropriate established regulatory concentration . Total areas for each plume were calculated as an indicator of relative environmental damage. These plume extents are shown graphically and in tabular form for comparison to previous estimates. Plume data also were interpolated to a finer grid (10 × 10 m) for some processing, particularly to estimate volumes of contaminated groundwater. However, hydrogeologic transport modeling was not considered for the interpolation. The compilation of plume extents for each contaminant also allowed estimates of overlap of the plumes or areas with more than one contaminant above regulatory standards. A mapping of saturated aquifer thickness also was derived across the 100-K and 100–N study area, based on the vertical difference between the groundwater level (water table) at the top and the altitude of the top of the Ringold Upper Mud geologic unit, considered the bottom of the uppermost unconfined aquifer. Saturated thickness was calculated for each cell in the finer (10 × 10 m) grid. The summation of the cells’ saturated thickness values within each polygon of plume regulatory exceedance provided an estimate of the total volume of contaminated aquifer, and the results also were checked using a SURFER ® volumetric integration procedure. The total volume of contaminated groundwater in each plume was derived by multiplying the aquifer saturated thickness volume by a locally representative value of porosity (0.3). Estimates of the uncertainty of the plume delineation also are presented. “Upper limit” plume delineations were calculated for each contaminant using the same procedure as the “average” plume extent except with values at each well that are set at a 95-percent upper confidence limit around the log-normally transformed mean concentrations, based on the standard error for the distribution of the mean value in that well; “lower limit” plumes are calculated at a 5-percent confidence limit around the geometric mean. These upper- and lower-limit estimates are considered unrealistic because the statistics were increased or decreased at each well simultaneously and were not adjusted for correlation among the well distributions (i.e., it is not realistic that all wells would be high simultaneously). Sources of the variability in the distributions used in the upper- and lower-extent maps include time varying concentrations and analytical errors. The plume delineations developed in this study are similar to the previous plume descriptions developed by U.S. Department of Energy and its contractors. The differences are primarily due to data selection and interpolation methodology. The differences in delineated plumes are not sufficient to result in the Hanford Natural Resource Trustee Council adjusting its understandings of contaminant impact or remediation.

Washington↗

USGS2016

Numerical simulation of groundwater flow at Puget Sound Naval Shipyard, Naval Base Kitsap, Bremerton, Washington

Information about groundwater-flow paths and locations where groundwater discharges at and near Puget Sound Naval Shipyard is necessary for understanding the potential migration of subsurface contaminants by groundwater at the shipyard. The design of some remediation alternatives would be aided by knowledge of whether groundwater flowing at specific locations beneath the shipyard will eventually discharge directly to Sinclair Inlet of Puget Sound, or if it will discharge to the drainage system of one of the six dry docks located in the shipyard. A 1997 numerical (finite difference) groundwater-flow model of the shipyard and surrounding area was constructed to help evaluate the potential for groundwater discharge to Puget Sound. That steady-state, multilayer numerical model with homogeneous hydraulic characteristics indicated that groundwater flowing beneath nearly all of the shipyard discharges to the dry-dock drainage systems, and only shallow groundwater flowing beneath the western end of the shipyard discharges directly to Sinclair Inlet. Updated information from a 2016 regional groundwater-flow model constructed for the greater Kitsap Peninsula was used to update the 1997 groundwater model of the Puget Sound Naval Shipyard. That information included a new interpretation of the hydrogeologic units underlying the area, as well as improved recharge estimates. Other updates to the 1997 model included finer discretization of the finite-difference model grid into more layers, rows, and columns, all with reduced dimensions. This updated Puget Sound Naval Shipyard model was calibrated to 2001–2005 measured water levels, and hydraulic characteristics of the model layers representing different hydrogeologic units were estimated with the aid of state-of-the-art parameter optimization techniques. The flow directions and discharge locations predicted by this updated model generally match the 1997 model despite refinements and other changes. In the updated model, most groundwater discharge recharged within the boundaries of the shipyard is to the dry docks; only at the western end of the shipyard does groundwater discharge directly to Puget Sound. Particle tracking for the existing long-term monitoring well network suggests that only a few wells intercept groundwater that originates as recharge within the shipyard boundary.

Washington↗

USGS2015

Hydrogeologic framework, groundwater movement, and water budget in the Puyallup River Watershed and vicinity, Pierce and King Counties, Washington

This report presents information used to characterize the groundwater-flow system in the Puyallup River Watershed and vicinity, 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 study area covers about 1,220 square miles in northern Pierce and southern King Counties, Washington; extends north to the Green River and Auburn Valley and southwest to the Puyallup River and adjacent uplands; and is bounded on the south and east by foothills of the Cascade Range and on the west by Puget Sound. The area is underlain by a northwest-thickening sequence of unconsolidated glacial and interglacial deposits, which overlie sedimentary and volcanic bedrock units that crop out in the foothills along the southern and eastern margin of the study area. Geologic units were grouped into 13 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 1,012 drillers’ logs to construct 8 hydrogeologic sections, and unit extent and thickness maps. Groundwater in unconsolidated glacial and interglacial aquifers generally flows to the northwest towards Puget Sound, and to the north and northeast towards the Puyallup River, White River, and Green River valleys. These generalized flow patterns are complicated by the presence of low permeability confining units and bedrock that separate discontinuous bodies of aquifer material and act as local groundwater-flow barriers. Water levels in wells completed in the unconsolidated hydrogeologic units show seasonal variations ranging from less than 1 to about 32 feet during the monitoring period (March 2011–March 2013). Synoptic streamflow measurements made in October 2011 and October 2012 indicated a total groundwater discharge to streams in the water-budget area (520 square miles located within the larger study area) of at least 349,000 and 280,000 acre-feet per year, respectively. Annual groundwater discharge to streams likely exceeds these values because streamflow measurements were made during the dry, late-summer and early-autumn period when groundwater levels typically are at annual lows. Most stream reaches in the study area either gain flow from groundwater discharge or exhibit near-neutral conditions with no substantial gain or loss of flow. Groundwater discharge occurs at numerous springs in the area; the total reported discharge of springs in the area is approximately 80,300 acre-feet per year. The water-budget area received about 1,428,000 acre-feet or about 52 inches of precipitation per year (January 1, 2011, to December 31, 2012). About 41 percent of precipitation enters the groundwater system as recharge. Seven percent of this recharge is withdrawn from wells and the remainder leaves the groundwater system as discharge to rivers, discharge to springs, or submarine discharge to Puget Sound, or exits the study area through subsurface flow in the Green River valley.

Washington↗

USGS2013

Transient calibration of a groundwater-flow model of Chimacum Creek Basin and vicinity, Jefferson County, Washington: a supplement to Scientific Investigations Report 2013-5160

A steady-state groundwater-flow model described in Scientific Investigations Report 2013-5160, ”Numerical Simulation of the Groundwater-Flow System in Chimacum Creek Basin and Vicinity, Jefferson County, Washington” was developed to evaluate potential future impacts of growth and of water-management strategies on water resources in the Chimacum Creek Basin. This supplement to that report describes the unsuccessful attempt to perform a calibration to transient conditions on the model. The modeled area is about 64 square miles on the Olympic Peninsula in northeastern Jefferson County, Washington. The geologic setting for the model area is that of unconsolidated deposits of glacial and interglacial origin typical of the Puget Sound Lowlands. The hydrogeologic units representing aquifers are Upper Aquifer (UA, roughly corresponding to recessional outwash) and Lower Aquifer (LA, roughly corresponding to advance outwash). Recharge from precipitation is the dominant source of water to the aquifer system; discharge is primarily to marine waters below sea level and to Chimacum Creek and its tributaries. The model is comprised of a grid of 245 columns and 313 rows; cells are a uniform 200 feet per side. There are six model layers, each representing one hydrogeologic unit: (1) Upper Confining unit (UC); (2) Upper Aquifer unit (UA); (3) Middle Confining unit (MC); (4) Lower Aquifer unit (LA); (5) Lower Confining unit (LC); and (6) Bedrock unit (OE). The transient simulation period (October 1994–September 2009) was divided into 180 monthly stress periods to represent temporal variations in recharge, discharge, and storage. An attempt to calibrate the model to transient conditions was unsuccessful due to instabilities stemming from oscillations in groundwater discharge to and recharge from streamflow in Chimacum Creek. The model as calibrated to transient conditions has mean residuals and standard errors of 0.06 ft ±0.45 feet for groundwater levels and 0.48 ± 0.06 cubic feet per second for flows. Although the expected seasonal trends were observed in model results, the typical observed annual variation of groundwater levels of about 2 feet was not. Streamflow at the most downstream observation point was about three times larger than simulated streamflow. Because the transient version of the model proved inherently unstable, it was not used to simulate forecast conditions for alternate hydrologic or anthropogenic changes. Adaptation of alternate stream simulation packages, such as RIV, or newer versions of MODFLOW, such as MODFLOW-NWT, could possibly assist with achieving calibration to transient conditions.

Washington↗

USGS2013

Effect of simulated tree canopy removal on a municipal wellfield in the Puget Sound aquifer system, Thurston County, Washington

Effects of tree canopy removal on a wellfield were simulated using a groundwater flow model characteristic of hydrogeologic settings in the Puget Sound aquifer system. Effects were estimated according to simulated changes in flow patterns that may result from tree canopy removal associated with varying degrees of residential development. The flow model used was a modified version of a model of the hydrogeologic setting in Thurston County, Washington; the wellfield was one planned for Olympia, Washington, and the canopy modifications spanned a range of possible land use change scenarios. The relative effects of tree canopy removal were estimated in terms of potential changes in capture zones for the wellfield and groundwater levels. Because of the depth of the wellfield and the dispersal of the effects from changes in recharge at ground surface, potential changes in wellfield capture zones and groundwater levels were discernible but small compared to other possible influences.

Washington↗

USGS2013

Numerical simulation of the groundwater-flow system in Chimacum Creek Basin and vicinity, Jefferson County, Washington

A groundwater-flow model was developed to evaluate potential future effects of growth and of water-management strategies on water resources in the Chimacum Creek Basin. The model covers an area of about 64 square miles (mi 2 ) on the Olympic Peninsula in northeastern Jefferson County, Washington. The Chimacum Creek Basin drains an area of about 53 mi 2 and consists of Chimacum Creek and its tributary East Fork Chimacum Creek, which converge near the town of Chimacum and discharge to Port Townsend Bay near the town of Irondale. The topography of the model area consists of north-south oriented, narrow, regularly spaced parallel ridges and valleys that are characteristic of fluted glaciated surfaces. Thick accumulations of peat occur along the axis of East Fork Chimacum Creek and provide rich soils for agricultural use. The study area is underlain by a north-thickening sequence of unconsolidated glacial (till and outwash) and interglacial (fluvial and lacustrine) deposits, and sedimentary and igneous bedrock units that crop out along the margins and the western interior of the model area. Six hydrogeologic units in the model area form the basis of the groundwater-flow model. They are represented by model layers UC (upper confining), UA (upper aquifer), MC (middle confining), LA (lower aquifer), LC (lower confining), and OE (bedrock). Groundwater flow in the Chimacum Creek Basin and vicinity was simulated using the groundwater-flow model, MODFLOW-2005. The finite-difference model grid comprises 245 columns, 313 rows, and 6 layers. Each model cell has a horizontal dimension of 200 × 200 feet (ft). The thickness of model layers varies throughout the model area and ranges from 5 ft in the non-bedrock units to more than 2,400 ft in the bedrock. Groundwater flow was simulated for steady-state conditions, which were simulated for calibration of the model using average recharge, discharge, and water levels for the 180-month period October 1994–September 2009. The model as calibrated has a mean residual of 4.5 ft and a standard error on the mean of 2.1 ft for heads, and 0.64±0.42 cubic feet per second for streamflows. After the model was calibrated, a Current Conditions simulation was developed to reflect current (October 2008–September 2009) hydrologic conditions, with representative pumping, return flows, and “normal” recharge (based on National Weather Service average precipitation for 1981 to 2010). The Current Conditions simulation was used to estimate current flow quantities, and as a basis to compare other simulations.Simulated steady-state inflow to the model area from precipitation and secondary recharge, or “return flow,” was 16,347 acre-feet per year (acre-ft/yr); groundwater inflow from other basins to the north of the model boundary was 1,518 acre-ft/yr (net, 3,114 acre-ft/yr in and 1,596 acre-ft/yr out) and simulated inflow from lake leakage was 613 acre-ft/yr (net, 684 acre-ft/yr in and 71 acre-ft/yr out). Simulated outflow from the model primarily was through discharge to Puget Sound (10,022 acre-ft/yr), streams (5,424 acre-ft/yr ), springs and seeps (1,521 acre-ft/yr), and through withdrawals from wells (1,506 acre-ft/yr). Four simulations were formulated using the calibrated model—one to represent current conditions (2009, the end of the period used for calibration) and three to provide representative examples of how the model can be used to evaluate the relative effects of potential changes in groundwater withdrawals and consumptive use on groundwater levels and stream base flows: Probable Future Use, based on population projections; Full Beneficial Use, based on Jefferson County Public Utility District #1 water rights; Sanitary Sewer, based on eliminating septic return flows in the Urban Growth Area. Particle tracking was used to assess flowpaths from sources and to sinks, and the effects of the presence of irrigation wells and their depths was assessed.

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↗

USGS2010

Hydrogeologic framework, groundwater movement, and water budget in the Chambers-Clover Creek watershed and vicinity, Pierce County, Washington

This report presents information used to characterize the groundwater-flow system in the Chambers-Clover Creek Watershed and vicinity, 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 study area covers about 706 square miles in western Pierce County, Washington, and extends north to the Puyallup River, southwest to the Nisqually River, and is bounded on the south and east by foothills of the Cascade Range and on the west by Puget Sound. The area is underlain by a northwest-thickening sequence of unconsolidated glacial and interglacial deposits which overlie sedimentary and volcanic bedrock units that crop out in the foothills along the southern and southeastern margin of the study area. Geologic units were grouped into 11 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 450 drillers’ logs to construct 6 hydrogeologic sections, and unit extent and thickness maps. Groundwater in unconsolidated glacial and interglacial aquifers generally flows to the northwest towards Puget Sound, and to the north and northeast towards the Puyallup River. 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. Water levels in wells completed in the unconsolidated hydrogeologic units show seasonal variations ranging from less than 1 to about 50 feet. The largest groundwater-level fluctuation (78 feet) observed during the monitoring period (March 2007–September 2008) was in a well completed in the bedrock unit. Synoptic streamflow measurements made in September 2007 and July 2008 indicated a total groundwater discharge to streams in the study area of 87,310 and 92,160 acre-feet per year, respectively. The synoptic streamflow measurements show a complex pattern of gains and losses to streamflows that varies throughout the study area, and appears to be influenced in places by local topography. Groundwater discharge occurs at numerous springs in the area and the total previously reported discharge of springs in the area is approximately 80,000 acre-feet per year. There are, in addition, many unmeasured springs and the total spring discharge in the area is unknown. The water-budget area (432 mi 2 located within the larger study area) received an annual average (September1, 2006, to August 31, 2008) of about 1,025,000 acre-ft or about 45 inches of precipitation a year. About 44 percent of precipitation enters the groundwater system as recharge. Almost one-half of this recharge (49 percent) discharges to the Puyallup and Nisqually Rivers and leaves the groundwater system as submarine groundwater discharge to Puget Sound. The remaining groundwater recharge discharges to streams (20 percent) and springs (18 percent) or is withdrawn from wells (13 percent)

Washington↗

USGS2010

Numerical simulation of the groundwater-flow system in tributary subbasins and vicinity, lower Skagit River basin, Skagit and Snohomish Counties, Washington

A groundwater-flow model was developed to evaluate the effects of potential groundwater withdrawals and consumptive use on streamflows in tributary subbasins of the lower portion of the Skagit River basin. The study area covers about 155 square miles along the Skagit River and its tributary subbasins (East Fork Nookachamps Creek, Nookachamps Creek, Carpenter Creek, Fisher Creek) in southwestern Skagit County and northwestern Snohomish County, Washington. The Skagit River occupies a large, relatively flat alluvial valley that extends across the northern and western margins of the study area, and is bounded to the south and east by upland and mountainous terrain. The alluvial valley and upland are underlain by unconsolidated deposits of glacial and inter- glacial origin. Bedrock underlies the alluvial valley and upland areas, and crops out throughout the mountainous terrain. Nine hydrogeologic units are recognized in the study area and form the basis of the groundwater-flow model. Groundwater flow in tributary subbasins of the lower Skagit River and vicinity was simulated using the groundwater-flow model, MODFLOW-2000. The finite-difference model grid consists of 174 rows, 156 columns, and 15 layers. Each model cell has a horizontal dimension of 500 by 500 feet. The thickness of model layers varies throughout the model area. Groundwater flow was simulated for both steady-state and transient conditions. The steady-state condition simulated average recharge, discharge, and water levels for the period, August 2006-September 2008. The transient simulation period, September 2006-September 2008, was divided into 24 monthly stress periods. Initial conditions for the transient model were developed from a 6-year ?lead-in? period that used recorded precipitation and Skagit River levels, and extrapolations of other boundary conditions. During model calibration, variables were adjusted within probable ranges to minimize differences between measured and simulated groundwater levels and stream baseflows. The final calibrated steady-state and transient models have weighted mean residual of -10.1 and -2.2 feet, respectively (negative residuals indicate that measured value is less than simulated value). Simulated inflow to the model area was about 144,000 acre-feet per year (acre-ft/yr) (81 percent of simulated inflow) from precipitation and secondary recharge, and about 32,700 acre-ft/yr (19 percent of simulated inflow) from stream and lake leakage. Simulated outflow from the model primarily was through discharge to streams and lakes (about 166,500 acre-ft/yr; 94 percent of simulated outflow), and withdrawals from wells (about 9,800 acre-ft/yr; 6 percent of simulated outflow). Model simulations were conducted to demonstrate model performance and to provide representative examples of how the model may be used to evaluate the effects of potential changes in groundwater withdrawals, consumptive use, and recharge on groundwater levels and tributary stream baseflows.

Scientific Investigations Report↗

USGS2009

Shallow Groundwater Movement in the Skagit River Delta Area, Skagit County, Washington

Shallow groundwater movement in an area between the lower Skagit River and Puget Sound was characterized by the U.S. Geological Survey to assist Skagit County and the Washington State Department of Ecology with the identification of areas where water withdrawals from existing and new wells could adversely affect streamflow in the Skagit River. The shallow groundwater system consists of alluvial, lahar runout, and recessional outwash deposits composed of sand, gravel, and cobbles, with minor lenses of silt and clay. Upland areas are underlain by glacial till and outwash deposits that show evidence of terrestrial and shallow marine depositional environments. Bedrock exposures are limited to a few upland outcrops in the southwestern part of the study area, and consist of metamorphic, sedimentary, and igneous rocks. Water levels were measured in 47 wells on a quarterly basis (August 2007, November 2007, February 2008, and May 2008). Measurements from 34 wells completed in the shallow groundwater system were used to construct groundwater-level and flow-direction maps and perform a linear-regression analysis to estimate the overall, time averaged shallow groundwater-flow direction and gradient. Groundwater flow in the shallow groundwater system generally moves in a southwestward direction away from the Skagit River and toward the Swinomish Channel and Skagit Bay. Local groundwater flow towards the river was inferred during February 2008 in areas west and southwest of Mount Vernon. Water-level altitudes varied seasonally, however, and generally ranged from less than 3 feet (August 2007) in the west to about 15 feet (May 2008) in the east. The time-averaged, shallow groundwater-flow direction derived from regression analysis, 8.5 deg south of west, was similar to flow directions depicted on the quarterly water-level maps. Seasonal changes in groundwater levels in most wells in the Skagit River Delta follow a typical pattern for shallow wells in western Washington. Water levels rise from October through March, when precipitation is high, and decline from April through September, when precipitation is lower. Groundwater levels in wells along the eastern margin of the study area also are likely influenced by stage on the Skagit River. Water levels in these wells remained elevated through April, and did not seem to begin to decline until the end of May in response to declining river stage. Groundwater levels in a well equipped with a continuous water-level recorder exhibited periodic fluctuations that are characteristic of ocean tides. This well is less than 1 mile east of the tidally influenced Swinomish Channel, and exhibited water-level fluctuations that correspond closely to predicted tidal extremes obtained from a tide gage near La Conner, Washington.

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↗