Geology topics

Tamara M. Wood

Publications and source records attributed to Tamara M. Wood.

At least 19 recordsLinked to original sources

USGS2023

Implications of water, sediment, and nutrient budgets for the restoration of a shallow, turbid lake in semiarid southeastern Oregon

Malheur Lake is the largest lake in the endorheic Harney Basin in southeastern Oregon. Since the 1990s, Malheur Lake—which averages depths of about 1 meter—has been in a degraded, turbid state lacking submergent and emergent vegetation. The goals of this study were to identify the major sources of sediment and nutrients to Malheur Lake to determine the importance of managing nutrients for lake restoration. Discrete water samples were analyzed for nutrient (total phosphorus, total nitrogen, orthophosphate, nitrate+nitrite, and ammonia) and suspended-sediment concentrations, and additional parameters including chlorophyll- a and phytoplankton biomass were measured in lake samples. Lake area fluctuated from a minimum of 3,300 hectares (ha) to a maximum of 11,300 ha in water years 2019 and 2020. In water year 2019, inflow from the tributaries created a 1,400-hectare area in the lake with low turbidity that persisted for multiple months. Land-use practices and water diversions along the tributaries affected the hydrographs and nutrient and suspended-sediment concentrations reaching the lake. As lake area increased, storage of sediment-associated constituents in the water column increased in excess of external loads because of resuspension. In 2019, 69 percent of the increase in suspended-sediment storage in the water column was attributed to internal resuspension and 31 percent was from external loading. Sediment was deposited as lake area decreased, and water-column storage decreased even as positive external loading continued. The internal resuspension, deposition, and external loading of suspended sediment likely is decreasing topographic heterogeneity in the lake. Concentrations of total phosphorus and orthophosphate are substantially higher than in the 1980s, and the lake is eutrophic. Phytoplankton in the lake was light limited in 2019–20, and restoration actions that prioritize vegetation establishment would reduce bioavailable nutrients for phytoplankton while increasing light in the water column.

Oregon↗

USGS2023

Future climate-induced changes in mixing and deep oxygen content of a caldera lake with hydrothermal heat and salt inputs

Vertical profiles of temperature, salinity and dissolved oxygen in Crater Lake, a caldera lake in the Oregon Cascade Range that receives hydrothermal inputs of heat and salt, were simulated with a 1-dimensional model. Twelve Global Circulation Models and two Representative Concentration Pathways (RCPs) were used to develop boundary conditions from 1950 to 2099. The model simulated the ventilation of deep water initiated by reverse stratification and subsequent thermobaric instability. All models predicted a reduction in the frequency of deep ventilation events, from an ensemble median frequency of 5.4 winters decade−1 during 1950–2005 to 4.3 (RCP4.5) or 2.5 (RCP8.5) winters decade−1 during 2045–2099. Favorable conditions for thermobaric instability-induced mixing currently occur infrequently and will become rare in the future. The salinity gradient resulting from hydrothermal inputs presents an additional barrier to thermobaric instability that will continue through 2099. A redistribution of salt to the deep lake may prevent ventilation all the way to the bottom in the future. Hypolimnetic dissolved oxygen percent saturation remained above 75% within the 21st century, consistent with oligotrophy and very small oxygen demands. The rate of change in all variables accelerated approaching 2099, coincident with elimination of winter reverse stratification. Historically, about half of the hydrothermal heat added to Crater Lake has been vented to the atmosphere. In the RCP8.5 scenario, the hydrothermal heat will cease to be vented to the atmosphere by the end of the 21st century, and then the temperature of the deep waters will increase rapidly.

Oregon↗

USGS2022

Light attenuation and erosion characteristics of fine sediments in a highly turbid, shallow, Great Basin Lake—Malheur Lake, Oregon, 2017–18

Malheur Lake is a large, shallow, turbid lake in southeastern Oregon that fluctuates widely in surface area in response to yearly precipitation and climatic cycles. High suspended-sediment concentrations (SSCs) likely are negatively affecting the survival of aquatic plants by reducing the intensity of solar radiation reaching the plants, thus inhibiting photosynthesis. This study was designed to determine the types of suspended material, the erodibility of the lakebed, the attenuation of photosynthetically active radiation (PAR) through the water column, and the effects of wind and precipitation on SSC. Two sites in the lake were monitored for approximately 5 months during the summer growing season each year (2017–19). At these sites, turbidity, chlorophyll a fluorescence (a surrogate for concentration), and underwater PAR measurements were collected continuously, and discrete samples were collected every 2 weeks and analyzed for SSC, loss on ignition, and chlorophyll a concentration. Underwater PAR profile measurements were collected during site visits, and a nearby meteorological station recorded terrestrial PAR and wind speeds. About 18 percent of suspended material in the water was organic and mostly detrital. Nearly 100 percent of all suspended material was fine material (less than 63 micrometers), and more than 90 percent of the surficial lakebed material was fine material. The high concentrations of fine material in the water column can be expected to strongly attenuate light. SSC was significantly higher at both sites in 2018 compared to 2017 and 2019; the interannual differences were mostly due to the lower amount of precipitation in 2018, which resulted in shallower lake depths. Three years of SSC values multiplied by water depth showed a seasonal pattern: concentrations were often highest in early spring, lowest in summer, and intermediate in autumn. Episodic wind events with speeds of 5–10 meters per second caused rapid increases in turbidity above background that lasted for a few days. However, a baseline SSC value multiplied by water depth (estimated to be 0.11 kilograms per square meter) was present between wind events and even under ice, suggesting a persistent suspension of very fine, highly erodible material. Terrestrial and underwater PAR measurements were used to develop a relation between PAR attenuation and turbidity that can be used in modeling restoration scenarios. Calculated bottom shear stress caused by wind-generated waves ranged from 0 to 0.4 pascals (Pa). Erosion experiments indicated variability in the bottom sediments from the two lake sites, but much of the lakebed is highly erodible at a threshold of 0.05 to 0.1 Pa. Restoration actions may target the persistent turbidity (for example, the use of flocculation) or transient turbidity (for example, construction of wave-reduction barriers), with a goal of attaining approximately 36 micromoles photons per square meter per second of PAR at the lakebed to promote emergence of sago pondweed and other desirable plants. Currently, that threshold often is reached from 4 to 34 centimeters (cm) below the water surface in 1 meter water depth, depending on wind conditions, but halving persistent turbidity would increase the upper end of the range to 55 cm. Additional studies regarding the effects of (1) sediment drying on resuspension and (2) nutrient inputs and internal cycling on phytoplankton populations would help determine the most appropriate restoration strategies.

Oregon↗

USGS2020

Use of boosted regression trees to quantify cumulative instream flow resulting from curtailment of irrigation in the Sprague River basin, Oregon

A boosted regression trees (BRT) approach was used to estimate the amount by which streamflow is increased when irrigation is regulated (curtailed) upstream of a streamgage on the Sprague River in southern-central Oregon. The BRT approach differs from most other approaches that require baseline conditions for comparison, where those baseline conditions are determined from past observations by searching for hydrologically similar years when irrigation was not regulated. Such baseline conditions are always imperfect estimates of the true baseline conditions. The BRT approach instead estimates unique baseline conditions for any year in which irrigation is regulated by calculating the baseline condition based on measurements of precipitation and weather observations that determine evapotranspiration, and other measurements that are proxies for the effects of climate and regional groundwater pumping on water-table elevation, using a model that has been trained in years of no regulation. The amount by which streamflow is increased by regulation is then calculated by subtracting the estimated baseline conditions from the measured streamflow. The approach is challenged by the fact that the streamflow increase may be a small fraction of the total streamflow; nonetheless, during 2 years in which regulation was started early and was implemented consistently through the season, the increased flow made up about one third of the flow past the streamgage during the regulation period. An advantage of this approach is that with rigorous model testing with holdout data, the threshold for detecting streamflow increase and intervals around the estimates of increase at a desired level of confidence can be quantified. The model relies on datasets that are readily available and updated continuously and therefore can be used operationally to inform resource management.

Oregon↗

USGS2019

Annual variations in microcystin occurrence in Upper Klamath Lake, Oregon, based on high-throughput DNA sequencing, qPCR, and environmental parameters

Cyanobacteria-dominated blooms in Upper Klamath Lake, Oregon, create poor water quality and produce microcystins that may be detrimental to local wildlife and human health. Genetic tools, including high-throughput DNA sequencing and quantitative polymerase chain reaction (qPCR), have been shown to improve the identification and quantification of key groups associated with these blooms over more traditional techniques. We examined the seasonal and interannual variations in nutrient (nitrogen and phosphorus) concentrations between 2013 and 2014 to describe the relations between these factors and the growth dynamics of Aphanizomenon and toxigenic Microcystis as described by DNA sequencing and qPCR. Although total nutrients and chlorophyll a concentrations were similar between years, qPCR results showed the cyanobacterial populations to be 40 times larger in 2014 and indicated a large shift from an Aphanizomenon -dominant, low microcystin-level regime in 2013 to one dominated later in the season by microcystin-producing Microcystis in 2014. In both years, the transition from Aphanizomenon to Microcystis was coincident with a late-season increase in nitrite-plus-nitrate concentrations and in dissolved inorganic nitrogen to dissolved inorganic phosphorus (DIN:DIP) ratios. However, these increases did not explain the large interannual differences in total cyanobacteria abundance. Rather, we hypothesized that year-to-year differences in bioavailable phosphorus, which also manifested as lower total nitrogen to total phosphorus (TN:TP) ratios, were responsible.

Oregon↗

USGS2019

Assessment of Columbia and Willamette River flood stage on the Columbia Corridor Levee System at Portland, Oregon, in a future climate

To support Levee Ready Columbia’s (LRC’s) effort to re-certify levees along the Columbia and Willamette Rivers and remain accredited, two 2-dimensional hydraulic models, Adaptive Hydraulics and Delft3D-Flexible Mesh, were used to simulate the effects of plausible extreme high water during the 2030 to 2059 period. The Columbia River was simulated from Bonneville Dam, situated at river mile (RM) 145, to the mouth of Columbia River, and the Willamette River was simulated from Willamette Falls, RM 26.2, to the Columbia River confluence. Inputs to the models included light detection and ranging (lidar) and bathymetric mapping data to determine bed level, and boundary conditions in the form of daily inflow hydrographs and water levels in the ocean offshore of the mouth of the Columbia River. Future conditions were based on climate science data developed by the U.S. Army Corps of Engineers and others. These conditions included future streamflow and coastal ocean water levels. The hypothetical, extreme but plausible, upstream boundary was based on scaling up the hydrographs from the 1996 flood. Scaling factors were determined by comparing the peak flow rankings determined from flood frequency analyses of historical unregulated periods and 2040s simulated unregulated winter streamflow. The comparison resulted in scaling up the Columbia River hydrograph by 40-percent and scaling up the Willamette River and Lower Columbia River tributaries hydrographs by 20-percent. The downstream ocean boundary was based on a combination of sea-level change, high tide, and storm surge. The models were calibrated for two historical periods: (1) from January 15 to February 28, 1996, and (2) from April 12 to July 12, 1997. The two models compared well to the measured water-surface elevation over the historical periods and had good performance statistics, with root-mean square error ranging from 0.085 to 0.32 meters, Nash-Sutcliffe values greater than 0.96, and bias ranging from -0.03 to 0.28 meters. The simulated peak stage in the Columbia River at Vancouver, Washington, for 1996 was 9.60 and 9.98 meters (31.5 and 32.7 feet) compared to the measured peak of 9.89 meters (32.5 feet). Future peak stage then was simulated with boundary conditions representing extreme but plausible future conditions at the inflow sites and the ocean boundary. The two calibrated models compared well in their simulations of extreme but plausible future conditions. For the 0-meter sea-level change scenario, the simulated peak stage in the Columbia River at Vancouver was 11.15 and 11.39 meters (36.6 and 37.4 feet); and for the 1-meter sea-level change scenario, the simulated peak stage in the Columbia River was 11.25 and 11.54 meters (36.9 and 37.9 feet). The total increase in stage as compared to the 1996 measured peak stage ranged from 1.26 to 1.65 meters (4.13 to 5.40 feet).

Oregon↗

USGS2018

A metabolism-based whole lake eutrophication model to estimate the magnitude and time scales of the effects of restoration in Upper Klamath Lake, south-central Oregon

A whole lake eutrophication (WLE) model approach for phosphorus and cyanobacterial biomass in Upper Klamath Lake, south-central Oregon, is presented here. The model is a successor to a previous model developed to inform a Total Maximum Daily Load (TMDL) for phosphorus in the lake, but is based on net primary production (NPP), which can be calculated from dissolved oxygen, rather than scaling up a small-scale description of cyanobacterial growth and respiration rates. This phase 3 WLE model is a refinement of the proof-of-concept developed in phase 2, which was the first attempt to use NPP to simulate cyanobacteria in the TMDL model. The calibration of the calculated NPP WLE model was successful, with performance metrics indicating a good fit to calibration data, and the calculated NPP WLE model was able to simulate mid-season bloom decreases, a feature that previous models could not reproduce. In order to use the model to simulate future scenarios based on phosphorus load reduction, a multivariate regression model was created to simulate NPP as a function of the model state variables (phosphorus and chlorophyll a ) and measured meteorological and temperature model inputs. The NPP time series was split into a low- and high-frequency component using wavelet analysis, and regression models were fit to the components separately, with moderate success. The regression models for NPP were incorporated in the WLE model, referred to as the “scenario” WLE (SWLE), and the fit statistics for phosphorus during the calibration period were mostly unchanged. The fit statistics for chlorophyll a , however, were degraded. These statistics are still an improvement over prior models, and indicate that the SWLE is appropriate for long-term predictions even though it misses some of the seasonal variations in chlorophyll a . The complete whole lake SWLE model, with multivariate regression to predict NPP, was used to make long-term simulations of the response to 10-, 20-, and 40-percent reductions in tributary nutrient loads. The long-term mean water column concentration of total phosphorus was reduced by 9, 18, and 36 percent, respectively, in response to these load reductions. The long-term water column chlorophyll a concentration was reduced by 4, 13, and 44 percent, respectively. The adjustment to a new equilibrium between the water column and sediments occurred over about 30 years.

Oregon↗

USGS2016

Simulation of deep ventilation in Crater Lake, Oregon, 1951–2099

The frequency of deep ventilation events in Crater Lake, a caldera lake in the Oregon Cascade Mountains, was simulated in six future climate scenarios, using a 1-dimensional deep ventilation model (1DDV) that was developed to simulate the ventilation of deep water initiated by reverse stratification and subsequent thermobaric instability. The model was calibrated and validated with lake temperature data collected from 1994 to 2011. Wind and air temperature data from three general circulation models and two representative concentration pathways were used to simulate the change in lake temperature and the frequency of deep ventilation events in possible future climates. The lumped model air2water was used to project lake surface temperature, a required boundary condition for the lake model, based on air temperature in the future climates. The 1DDV model was used to simulate daily water temperature profiles through 2099. All future climate scenarios projected increased water temperature throughout the water column and a substantive reduction in the frequency of deep ventilation events. The least extreme scenario projected the frequency of deep ventilation events to decrease from about 1 in 2 years in current conditions to about 1 in 3 years by 2100. The most extreme scenario considered projected the frequency of deep ventilation events to be about 1 in 7.7 years by 2100. All scenarios predicted that the temperature of the entire water column will be greater than 4 °C for increasing lengths of time in the future and that the conditions required for thermobaric instability induced mixing will become rare or non-existent. The disruption of deep ventilation by itself does not provide a complete picture of the potential ecological and water quality consequences of warming climate to Crater Lake. Estimating the effect of warming climate on deep water oxygen depletion and water clarity will require careful modeling studies to combine the physical mixing processes affected by the atmosphere with the multitude of factors affecting the growth of algae and corresponding water clarity.

Oregon↗

USGS2015

Revision and proposed modification for a total maximum daily load model for Upper Klamath Lake, Oregon

This report presents Phase 2 of the review and development of the mass balance water-quality model, originally developed in 2001, that guided establishment of the phosphorus (P) total maximum daily load (TMDL) for Upper Klamath and Agency Lakes, Oregon. The purpose of Phase 2 was to incorporate a longer (19-year) set of external phosphorus loading data into the lake TMDL model than had originally been available, and to develop a proof-of-concept method for modeling algal mortality and the consequent decrease in chlorophyll a that had not been possible with the 2001 TMDL model formulation. Using the extended 1991–2010 external phosphorus loading dataset, the lake TMDL model was recalibrated following the same procedures outlined in the Phase 1 review. The version of the model selected for further development incorporated an updated sediment initial condition, a numerical solution method for the chlorophyll a model, changes to light and phosphorus factors limiting algal growth, and a new pH-model regression, which removed Julian day dependence in order to avoid discontinuities in pH at year boundaries. This updated lake TMDL model was recalibrated using the extended dataset in order to compare calibration parameters to those obtained from a calibration with the original 7.5-year dataset. The resulting algal settling velocity calibrated from the extended dataset was more than twice the value calibrated with the original dataset, and, because the calibrated values of algal settling velocity and recycle rate are related (more rapid settling required more rapid recycling), the recycling rate also was larger than that determined with the original dataset. These changes in calibration parameters highlight the uncertainty in critical rates in the Upper Klamath Lake TMDL model and argue for their direct measurement in future data collection to increase confidence in the model predictions.

Oregon↗

USGS2014

Statistical analysis of the water-quality monitoring program, Upper Klamath Lake, Oregon, and optimization of the program for 2013 and beyond

Upper Klamath Lake in south-central Oregon has become increasingly eutrophic over the past century and now experiences seasonal cyanobacteria-dominated and potentially toxic phytoplankton blooms. Growth and decline of these blooms create poor water-quality conditions that can be detrimental to fish, including two resident endangered sucker species. Upper Klamath Lake is the primary water supply to agricultural areas within the upper Klamath Basin. Water from the lake is also used to generate power and to enhance and sustain downstream flows in the Klamath River. Water quality in Upper Klamath Lake has been monitored by the Klamath Tribes since the early 1990s and by the U.S. Geological Survey (USGS) since 2002. Management agencies and other stakeholders have determined that a re-evaluation of the goals for water-quality monitoring is warranted to assess whether current data-collection activities will continue to adequately provide data for researchers to address questions of interest and to facilitate future natural resource management decisions. The purpose of this study was to (1) compile an updated list of the goals and objectives for long-term water-quality monitoring in Upper Klamath Lake with input from upper Klamath Basin stakeholders, (2) assess the current water-quality monitoring programs in Upper Klamath Lake to determine whether existing data-collection strategies can fulfill the updated goals and objectives for monitoring, and (3) identify potential modifications to future monitoring plans in accordance with the updated monitoring objectives and improve stakeholder cooperation and data-collection efficiency. Data collected by the Klamath Tribes and the USGS were evaluated to determine whether consistent long-term trends in water-quality variables can be described by the dataset and whether the number and distribution of currently monitored sites captures the full range of environmental conditions and the multi-scale variability of water-quality parameters in the lake. Also, current monitoring strategies were scrutinized for unnecessary redundancy within the overall network.

Oregon↗

USGS2014

Particle-tracking investigation of the retention of sucker larvae emerging from spawning grounds in Upper Klamath Lake, Oregon

This study had two objectives: (1) to use the results of an individual-based particle-tracking model of larval sucker dispersal through the Williamson River delta and Upper Klamath Lake, Oregon, to interpret field data collected throughout Upper Klamath and Agency Lakes, and (2) to use the model to investigate the retention of sucker larvae in the system as a function of Williamson River flow, wind, and lake elevation. This is a follow-up study to work reported in Wood and others (2014) in which the hydrodynamic model of Upper Klamath Lake was combined with an individual-based, particle-tracking model of larval fish entering the lake from spawning areas in the Williamson River. In the previous study, the performance of the model was evaluated through comparison with field data comprising larval sucker distribution collected in 2009 by The Nature Conservancy, Oregon State University (OSU), and the U.S. Geological Survey, primarily from the (at that time) recently reconnected Williamson River Delta and along the eastern shoreline of Upper Klamath Lake, surrounding the old river mouth. The previous study demonstrated that the validation of the model with field data was moderately successful and that the model was useful for describing the broad patterns of larval dispersal from the river, at least in the areas surrounding the river channel immediately downstream of the spawning areas and along the shoreline where larvae enter the lake. In this study, field data collected by OSU throughout the main body of Upper Klamath Lake, and not just around the Williamson River Delta, were compared to model simulation results. Because the field data were collected throughout the lake, it was necessary to include in the simulations larvae spawned at eastern shoreline springs that were not included in the earlier studies. A complicating factor was that the OSU collected data throughout the main body of the lake in 2011 and 2012, after the end of several years of larval drift collection in the Williamson River by the U.S. Geological Survey. Those larval drift data provided necessary boundary-condition information for the earlier studies, but there were no measured boundary conditions for larval input into model simulations during the years of this study (2011−12). Therefore, we developed a method to estimate a time series of larval drift in the Williamson River, and of the emergence of larvae from the gravel at the eastern shoreline springs, that captured the approximate timing of the larval pulse of the Lost River sucker (Deltistes luxatus) and shortnose sucker (Chasmistes brevirostris) and the relative magnitude of the pulses by species and spawning location. The method is not able to predict larval drift on any given day, but it can reasonably predict the approximate temporal progression of the larval drift through the season, based on counts of adult suckers returning to spawn. The accuracy in the timing of the larval pulses is not better than about plus or minus 5 days. Model results and field data were consistent in the basic progression of both catch per unit effort (CPUE) and larval length through time. The model simulation results also duplicated some of the characteristics of the spatial patterns of density in the field data, notably the tendency for high larval densities closer to the eastern and western shorelines. However, the model simulations could not explain high densities in the northern part of the lake or far into Ball Bay, locations that are far from the source of larvae in the Williamson River or eastern shoreline springs (as measured along the predominant transport pathways simulated in the model). This suggests the possibility of unaccounted-for spawning areas in the northern part of the lake and also that the period during which larvae are transported passively by the currents is shorter than the 46 days simulated in the model. Similarly, the progression of larval lengths in the field data is not a simple progression from smaller to larger fish away from sources in the river and springs, as simulated by the particle-tracking model; the smallest fish were caught at different times near the Williamson River, in the northwestern part of the lake, and in the southernmost part of the lake. This again suggests that fish may be spawning at places other than the river and eastern springs, that our understanding of larval transport is incomplete, or both. The model was used to run 96 numerical “experiments” in which lake elevation, river discharge, and wind forcing were varied systematically in order to investigate the sensitivity of particle retention to each variable, and with particular emphasis on the idea of managing lake elevation to control emigration. The estimates of particle retention cannot be equated directly to retention of fish larvae, primarily because there was no mortality included in the simulations, but the relative comparison of retention and emigration around the matrix of experimental conditions provided several “big picture” results: - Variables that cannot be controlled—winds and discharge—had the largest effect on retention. For example, at the lowest river discharge (20 cubic meters per second), simulated retention was high regardless of wind or lake elevation, whereas at the highest river discharge (100 cubic meters per second), retention was low regardless of wind or lake elevation. - When river discharge and wind were held constant, a higher elevation delayed the onset of the most rapid exit of particles by 1 (from the springs) to 4 (from the river) days, but did not determine overall retention. Only under the combination of conditions consisting of low discharge (50 cubic meters per second or less) and strong wind reversals for several days was there a consistent effect of lake elevation on overall retention several weeks into the simulation, and, under those conditions, retention was at the high end of the possible range regardless of lake elevation. - Under most combinations of conditions tested, after particles had been in the system for several days, the complex interaction between wind, elevation, and river discharge resulted in particle pathways, and therefore retention, being highly variable and unpredictable, at which point controlling lake elevation could not produce a predictable result. Therefore, on the basis of the model predictions, managing lake elevation probably is not a way to reliably provide any particular level of retention.

Oregon↗

USGS2014

Simulation and validation of larval sucker dispersal and retention through the restored Williamson River Delta and Upper Klamath Lake system, Oregon

A hydrodynamic model with particle tracking was used to create individual-based simulations to describe larval fish dispersal through the restored Williamson River Delta and into Upper Klamath Lake, Oregon. The model was verified by converting particle ages to larval lengths and comparing these lengths to lengths of larvae in net catches. Correlations of simulated lengths with field data were moderate and suggested a species-specific difference in model performance. Particle trajectories through the delta were affected by wind speed and direction, lake elevation, and shoreline configuration. Once particles entered the lake, transport was a function of current speed and whether behavior enhanced transport (swimming aligned with currents) or countered transport through greater dispersal (faster random swimming). We tested sensitivity to swim speed (higher speeds led to greater dispersal and more retention), shoreline configuration (restoration increased retention relative to pre-restoration conditions), and lake elevation (retention was maximized at an intermediate elevation). The simulations also highlight additional biological questions, such as the extent to which spatially heterogeneous mortality or fish behavior and environmental cues could interact with wind-driven currents and contribute to patterns of dispersal.

Oregon↗

USGS2013

U.S. Geological Survey water science strategy—Observing, understanding, predicting, and delivering water science to the Nation

Executive Summary This report expands the Water Science Strategy that began with the USGS Science Strategy, “Facing Tomorrow’s Challenges—U.S. Geological Survey Science in the Decade 2007–2017” (U.S. Geological Survey, 2007). This report looks at the relevant issues facing society and develops a strategy built around observing, understanding, predicting, and delivering water science for the next 5 to 10 years by building new capabilities, tools, and delivery systems to meet the Nation’s water-resource needs. This report begins by presenting the vision of water science for the USGS and the societal issues that are influenced by, and in turn influence, the water resources of our Nation. The essence of the Water Science Strategy is built on the concept of “water availability,” defined as spatial and temporal distribution of water quantity and quality, as related to human and ecosystem needs, as affected by human and natural influences. The report also describes the core capabilities of the USGS in water science—the strengths, partnerships, and science integrity that the USGS has built over its 134-year history.

Circular↗

USGS2013

Technical evaluation of a total maximum daily load model for Upper Klamath and Agency Lakes, Oregon

We reviewed a mass balance model developed in 2001 that guided establishment of the phosphorus total maximum daily load (TMDL) for Upper Klamath and Agency Lakes, Oregon. The purpose of the review was to evaluate the strengths and weaknesses of the model and to determine whether improvements could be made using information derived from studies since the model was first developed. The new data have contributed to the understanding of processes in the lakes, particularly internal loading of phosphorus from sediment, and include measurements of diffusive fluxes of phosphorus from the bottom sediments, groundwater advection, desorption from iron oxides at high pH in a laboratory setting, and estimates of fluxes of phosphorus bound to iron and aluminum oxides. None of these processes in isolation, however, is large enough to account for the episodically high values of whole-lake internal loading calculated from a mass balance, which can range from 10 to 20 milligrams per square meter per day for short periods. The possible role of benthic invertebrates in lake sediments in the internal loading of phosphorus in the lake has become apparent since the development of the TMDL model. Benthic invertebrates can increase diffusive fluxes several-fold through bioturbation and biodiffusion, and, if the invertebrates are bottom feeders, they can recycle phosphorus to the water column through metabolic excretion. These organisms have high densities (1,822–62,178 individuals per square meter) in Upper Klamath Lake. Conversion of the mean density of tubificid worms (Oligochaeta) and chironomid midges (Diptera), two of the dominant taxa, to an areal flux rate based on laboratory measurements of metabolic excretion of two abundant species suggested that excretion by benthic invertebrates is at least as important as any of the other identified processes for internal loading to the water column. Data from sediment cores collected around Upper Klamath Lake since the development of the TMDL model also contributed to this review. Cores were sequentially extracted to determine the distribution of phosphorus associated with several matrices in the sediment (freely exchangeable, metal-oxides, acid-soluble minerals, and residual). The concentrations of phosphorus in these fractions varied around the lake in patterns that reflect transport processes in the lake and the ultimate deposition of organic and inorganic forms of phosphorus from the water column. Both organic and inorganic phosphorus had higher concentrations in the northern part of the lake, in and just west of Goose Bay. At the time that these cores were collected, prior to restoration of the Williamson River Delta, this area was close to the shoreline of the lake and east of the Williamson River mouth. This contrasts with erosional inputs, which, in addition to being high to the east of the pre-restoration Williamson River mouth, were higher in the middle of the lake than at the northern end. Organic forms of phosphorus had particularly high concentrations in the northern bays. When these cores were used to calculate a new estimate of the whole-lake-averaged concentration of total phosphorus in the top 10 centimeters of the lake sediments, the estimate was about one-third of the best estimate available when the TMDL model was developed.

Oregon↗

USGS2012

Time scales of change in chemical and biological parameters after engineered levee breaches adjacent to Upper Klamath and Agency Lakes, Oregon

Eight sampling trips were coordinated after engineered levee breaches hydrologically reconnected both Upper Klamath Lake and Agency Lake, Oregon, to adjacent wetlands. The reconnection, by a series of explosive blasts, was coordinated by The Nature Conservancy to reclaim wetlands that had for approximately seven decades been leveed for crop production. Sets of nonmetallic porewater profilers (U.S. Patent 8,051,727 B1; November 8, 2011; http://www.uspto.gov/web/patents/patog/ week45/OG/html/1372-2/US08051727-20111108.html.) were deployed during these trips in November 2007, June 2008, May 2009, July 2009, May 2010, August 2010, June 2011, and July 2011 (table 1). Deployments temporally spanned the annual cyanophyte bloom of Aphanizomenon flos-aquae and spatially involved three lake and four wetland sites. Spatial and temporal variation in solute benthic flux was determined by the field team, using the profilers, over an approximately 4-year period beginning 3 days after the levee breaches. The highest flux to the water column of dissolved organic carbon (DOC) was detected in the newly flooded wetland, contrasting negative or insignificant DOC fluxes at adjacent lake sites. Over the multiyear study, DOC benthic fluxes dissipated in the reconnected wetlands, converging to values similar to those for established wetlands and to the adjacent lake (table 2). In contrast to DOC, benthic sources of soluble reactive phosphorus, ammonium, dissolved iron and manganese from within the reconnected wetlands were consistently elevated (that is, significant in magnitude relative to riverine and established-wetland sources) indicating a multi-year time scale for certain chemical changes after the levee breaches (table 2). Colonization of the reconnected wetlands by aquatic benthic invertebrates during the study trended toward the assemblages in established wetlands, providing further evidence of a multiyear transition of this area to permanent aquatic habitat (table 3). Both the lake and wetland benthic environments substantively contribute to macro- and micronutrients in the water column. Wetland areas undergoing restoration, and those being used for water storage, function very differently relatively to the established wetland within the Upper Klamath Lake National Wildlife Refuge, adjacent Upper Klamath Lake. Developing long-term management strategies for water quality in the Upper Klamath Basin requires recognition of the multi-year time scales associated with restoring wetlands that provide natural, seasonal ecosystem function and services.

Oregon↗

USGS2012

Strategic directions for U.S. Geological Survey water science, 2012-2022 - Observing, understanding, predicting, and delivering water science to the Nation

Executive Summary This report expands the Water Science Strategy that was begun in the USGS Science Strategy, “Facing Tomorrow’s Challenges—U.S. Geological Survey Science in the Decade 2007–2017” (U.S. Geological Survey, 2007). The report looks at the relevant issues facing society and develops a strategy built around observing, understanding, predicting, and delivering water science for the next 5 to 10 years by building new capabilities, tools, and delivery systems to meet the Nation’s water-resource needs. This report begins by presenting the vision of water science for the USGS and the societal issues that are influenced by, and in turn influence, the water resources of our Nation. The essence of the Water Strategic Science Plan is built on the concept of “water availability,” defined as spatial and temporal distribution of water quantity and quality, as related to human and ecosystem needs, as affected by human and natural influences . The report also describes the core capabilities of the USGS in water science—the strengths, partnerships, and science integrity that the USGS has built over its 130-year history. Nine priority actions are presented in the report, which combine and elevate the numerous specific strategic actions listed throughout the report. Priority actions were developed as a means of providing the audience of this report with a list for focused attention, even if resources and time limit the ability of managers to address all of the strategic actions in the report. Priority actions focus on the following: Improve integrated science planning for water. Expand and enhance water-resource monitoring networks. Characterize the water cycle through development of state-of-the-art 3-D/4-D hydrogeologic framework models at multiple scales. Clarify the linkage between human water use (engineered hydrology) and the water cycle (natural hydrology). Advance ecological flow science. Provide flood-inundation science and information. Develop rapid deployment teams for water-related emergencies. Conduct integrated watershed assessment, research, and modeling. Deliver water data and analyses to the Nation. The body of the report is presented as a hierarchal set of 5 goals, 14 objectives, and 27 strategic actions that the USGS should undertake to advance water science through year 2022. The goals deal with: Providing society the information it needs regarding the amount and quality of water in all components of the water cycle at high temporal and spatial resolution, nationwide; Advancing our understanding of processes that determine water availability; Predicting changes in the quantity and quality of water resources in response to changing climate, population, land use, and management scenarios; Anticipating and responding to water-related emergencies and conflicts; and Delivering timely hydrologic data, analyses, and decision-support tools seamlessly across the Nation to support water-resource decisions. Scientific information produced on water resources would be without value if it were not communicated to society in a fashion that can inform decisions and actions. Therefore, the chapter following the goals describes how the USGS should inform, involve, and educate society about the science it produces. This includes discussions on local outreach and the use of social media for effective communication. This report concludes with a chapter devoted to the crosscutting science issues of the Water Mission Area with the other USGS Mission Areas: Climate and Land Use Change, Core Science Systems, Ecosystems, Energy and Minerals, Environmental Health Science, and Natural Hazards. Not one of these Mission Areas stands alone—all must work together and integrate their actions to fulfill the USGS science mission for the future. This final chapter identifies the important linkages that must be realized and maintained for this integration to occur.

Open-File Report↗

USGS2012

Water-quality data from Upper Klamath and Agency Lakes, Oregon, 2009-10

The U.S. Geological Survey Upper Klamath Lake water-quality monitoring program collected data from multiparameter continuous water-quality monitors, weekly water-quality samples, and meteorological stations during 2009 and 2010 from May through November each year. The results of these measurements and sample analyses, as well as quality-control data for the water-quality samples, are presented in this report for 14 sites on Upper Klamath Lake and 2 sites on Agency Lake. These 2 years of data demonstrate a contrast in the seasonal bloom of the dominant cyanobacterium, Aphanizomenon flos-aquae , that can be related to differences in the measured water quality and meteorological variables. Some of the significant findings from 2009 and 2010 are listed below. * Both 2009 and 2010 were characterized by two cyanobacteria blooms, but the blooms differed in timing and intensity. The first bloom in 2009 peaked in late June and at higher chlorophyll a concentrations at most sites than the first bloom in 2010, which peaked in mid-July. A major decline in the first 2009 bloom occurred in late July and was followed by a second bloom that peaked at most sites in mid-August and persisted through September. The decline of the weaker first bloom in 2010 occurred in early August and was followed by a more substantial second bloom that peaked between late August and early September at most sites. * Dissolved oxygen minima associated with bloom declines occurred approximately 2 weeks earlier in 2009 (mid-July) than in 2010 (early August). pH maxima associated with rapid bloom growth occurred in late June and again in mid-August in 2009 and in mid-July and late August in 2010. * In both years, the maxima for total phosphorus and total nitrogen concentrations coincided with the chlorophyll a maximum. The maxima for dissolved nutrient concentrations (orthophosphate, ammonia, and nitrite plus nitrate) coincided with the declines of the first blooms. * Total particulate carbon, total particulate nitrogen, and total particulate phosphorus concentrations were measured in 2009 only. The ratios of carbon to phosphorus and nitrogen to phosphorus in particulates were the highest of the entire season during the rapid growth phase of the first bloom and were the lowest of the season during the decline of the first bloom. These ratios increased with the onset of the second bloom in that year, but to a lesser degree. * Meteorological data show that 2009 was warmer (particularly in June and July), less windy, and more humid early in the season than 2010. The difference in water temperatures reflected the difference in air temperatures in that the lakes were warmer in 2009 than in 2010 starting in early May, when the sensors were deployed, through most of June. Water temperature peaked at a higher value in 2009, and there were more clear days in June 2009 than in June 2010.

Oregon↗

USGS2012

Dependence of flow and transport through the Williamson River Delta, Upper Klamath Lake, Oregon, on wind, river inflow, and lake elevation

The hydrodynamic model of Upper Klamath and Agency Lakes, Oregon, was used to run 384 realizations of a numerical tracer experiment in order to understand the relative effects of wind, lake elevation, and Williamson River inflow on flow and transport (the movement of water and passively transported constituents) through the Williamson River Delta. Significant findings from this study include: * The replacement rate of water increased in Tulana and Goose Bay with increasing lake elevation, Williamson River inflow, and wind speed. * The fraction of Williamson River inflow passing through either side of the Delta increased with lake elevation and Williamson River inflow. * The partial replacement rate of water in Goose Bay with water from the Williamson River increased with wind speed. * The partial replacement rate of water in Tulana with water from the Williamson River decreased with wind speed. * Strong wind forcing at the water surface caused more of the Williamson River inflow to pass through Goose Bay than through Tulana. * Westerly to northwesterly winds result in more of the Williamson River inflow passing through the Goose Bay side of the Delta than through the Tulana side. * Regression models developed from the tracer experiments can be used to quantify the dependencies between transport and the independent variables to obtain rough estimates of useful quantities such as residence time and steady-state solute concentrations.

Oregon↗