Geology topics

Susan A. Wherry

Publications and source records attributed to Susan A. Wherry.

9 recordsLinked to original sources

USGS2023

Modeling flow and water quality in reservoir and river reaches of the Mahoning River Basin, Ohio

Executive Summary The U.S. Army Corps of Engineers (USACE) is considering changes to the management of water surface elevation in four lakes in the Mahoning River Basin. These changes would affect the timing and amounts of water released to the Mahoning River and could affect the water quality of those releases. To provide information on possible water-quality effects from these operational changes, flow and water-quality models were constructed for Berlin Lake, Lake Milton, Michael J Kirwan Reservoir, Mosquito Creek Lake, Mosquito Creek, and the Mahoning River from the dams downstream to Lowellville, Ohio. The models were calibrated for two calendar years each, with model years selected depending on the availability of water-quality data. Models were developed with CE-QUAL-W2 version 4.2 (Wells, S.A., 2020, CE-QUAL-W2—A two-dimensional, laterally averaged, hydrodynamic and water quality model [version 4.2]: Portland State University, variously paged), a two-dimensional, laterally averaged hydrodynamic and water-quality model. Modeled constituents included flow, velocity, ice cover, water temperature, total dissolved solids (TDS), sulfate, chloride, inorganic suspended sediment, nitrate, ammonia, total Kjeldahl nitrogen, orthophosphate, total phosphorus, dissolved and particulate organic matter, algae, and dissolved oxygen. Iron was included for the lake models, but not the river. A whole-basin model, with the four lake models and river model, was used to run model scenarios to examine the effects of altered lake water surface elevations on flow and water quality in the lakes, the lake outflows, and the Mahoning River. The initial whole-basin model, with calendar year 2013 hydrology and measured or typical water quality, was designated as scenario 0. Mahoning River flows for calendar year 2013 were close to a 20-year median flow. Four additional scenarios were constructed based on reservoir operations model (RES-SIM) model water surface elevations for the four lakes as provided by USACE. Scenario 1 was the RES-SIM base case, scenario 2 kept Berlin Lake water surface elevations higher in summer, scenario 3 allowed 25 percent of summer flood storage to extend the guide curve, and scenario 4 allowed more flexibility in lake management by removing any downstream Mahoning River minimum flow requirements. The Mahoning River model was not changed in any scenarios but received altered flows from the lakes. Significant findings from this study include the following: In two of the four lakes (Berlin and Mosquito Creek Lakes), development of lake model grids using recent bathymetric surveys suggests that sedimentation in these lakes has occurred since they were constructed, altering volume-elevation curves. Tests of model parameter sensitivity showed that modeled water temperature, TDS, and dissolved oxygen were relatively insensitive to model parameter values. Modeled chlorophyll a , a measure of algal concentration, was most sensitive to parameter values; nitrate and total phosphorus concentrations were affected by a few of the parameters tested. As a group, the lake model results were more sensitive to model parameter values compared to the Mahoning River model. Data gaps were identified for inflows, both for water quantity and water quality, that could be filled through future sampling programs. Ample data were available from within the waterbodies for model calibration. The model simulated the general spatial and temporal patterns of water temperature, TDS, chloride, sulfate, nutrients, suspended sediment, organic matter, chlorophyll a , and dissolved oxygen in the lakes and Mahoning River. From late spring to autumn in the years modeled (2006, 2013, 2017–19 depending on the lake), all lakes developed thermal stratification and periods of anoxia in bottom waters. Stratification was most stable in Michael J Kirwan Reservoir and least stable in Mosquito Creek Lake. The stratification and anoxia in Berlin Lake, Lake Milton, and Mosquito Creek Lake could be interrupted by high-flow inputs moving through those lakes. The model predicted the release of ammonia and iron during anoxic periods in the lake hypolimnions. Concentrations of TDS, nitrate, orthophosphate, and total phosphorus increased in the Mahoning River down to Lowellville, the end of the river model, in the years modeled. These concentrations were greater than those in upstream lake releases. Chloride and sulfate concentrations were underpredicted in the Mahoning River, suggesting the presence of unreported chloride and sulfate inputs to the river, at least in the years modeled. Model scenario 4 kept water surface elevations the highest in all lakes in the April to mid-December period, compared to scenarios 1–3. Model scenario 2 kept water surface elevations in Berlin Lake higher in summer and late autumn, compared to scenarios 1 and 3, but to satisfy downstream minimum flow requirements, water surface elevations in the other lakes had periods of lower water surface elevation. As a group, scenarios 1–3 had largely similar effects on flow and water surface elevation in the Mahoning River because the lake releases in those scenarios still met downstream Mahoning River flow targets. Modeling the removal of downstream flow targets, scenario 4 had periods of lower flow in the Mahoning River from April to mid-September as water was held in the lakes, and periods of higher Mahoning River flow from mid-September through November as the lakes were drawn down to prepare for winter flood-risk management. In the four scenarios, all the lakes and lake outflows had generally similar seasonal cycles of water quality, though some differences were predicted. For instance, higher concentrations of iron and ammonia in the Lake Milton hypolimnion were modeled during a period of both low inflows from Berlin Lake and low outflows at Lake Milton dam. It is possible that those changes could be minimized by maintaining more flow or pulses of higher flow through the lake. Compared to the scenario 1 base case, changes to Mahoning River water quality were relatively minor for scenarios 2 and 3, which maintained downstream flows but shifted the flow source among the upstream lakes. The largest changes in Mahoning River water quality were predicted between Leavittsburg and Lowellville for scenario 4. The periods of lower lake outflows between April and mid-September led to correspondingly higher concentrations of TDS, orthophosphate, total phosphorus, and nitrate in the river, compared to the base case scenario 1. Conversely, the overall greater lake outflows from mid-September through November in scenario 4 led to periods of lower concentrations of TDS and nutrients in that portion of the river, at that time of year.

Ohio↗

USGS2022

Summary of extreme water-quality conditions in Upper Klamath Lake, Oregon, 2005–19

This study used the complete set of continuous water-quality (WQ) data and discrete measurements of total ammonia collected by the U.S. Geological Survey from 2005 to 2019 at the four core sites in Upper Klamath Lake, Oregon, to examine relations between variables and extreme conditions that may be harmful for endemic Lost River suckers ( Deltistes luxatus ) and shortnose suckers ( Chasmistes brevirostris ). Several graphical and tabular approaches were used to compare variables, sites, and years to better understand the factors contributing to and timing of extreme WQ in the lake. Extreme WQ thresholds were defined as the 1st or 99th percentiles of the daily average dataset of water temperature, pH, and dissolved oxygen (DO) concentration, and the weekly estimated un-ionized ammonia (NH 3 ) from 2005 to 2019. Extreme WQ days were defined as those when at least 12 hours of measurements exceeded the extreme WQ threshold. The core site at Mid-Trench, which was also the deepest measurement site with a full-pool depth of 15 meters and at which water-quality sondes were deployed at the top and bottom of the water column, had the most extreme conditions of high water temperature, low DO, and high NH 3 . The upper sonde at Mid-Trench represented 40 percent of all days of extremely high water temperature (days with at least 12 hours exceeding 24.38 degrees Celsius) in the lake and 71 percent of all weekly estimates of extremely high NH 3 (greater than 264 micrograms per liter) in the lake. The lower sonde at Mid-Trench represented 85 percent of all days of extremely low DO (days with at least 12 hours of DO concentrations less than 1.76 milligrams per liter) in the lake. In each of the study years, poor water quality at Mid-Trench, as represented by several metrics, lasted for multiple days. The shallowest site at the Williamson River outlet represented 54 percent of all days of extremely high pH (days with at least 12 hours of pH measurements exceeding 10.04) in the lake. The seasonality of extreme WQ during the summer sampling period (limited to June through September) was evaluated and most days of extremely high water temperature (83 percent) and extremely high pH (54 percent) occurred in July, whereas most days of extremely low DO (57 percent) and extremely high NH 3 (57 percent) occurred in August. The years with the most days of extreme WQ accumulated for all variables (high water temperature, low DO, high pH, and high NH 3 ) were 2012–15 and 2017, which all occurred in the latter half of the study period. The years with the fewest accumulated days of extreme WQ were 2010 and 2011.

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↗

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↗