SEARCH · Geology Reports

Results for “Simulation & Gaming”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

29 records · Page 2Linked to original sources

USGS2017

Water-quality models to assess algal community dynamics, water quality, and fish habitat suitability for two agricultural land-use dominated lakes in Minnesota, 2014

Fish habitat can degrade in many lakes due to summer blue-green algal blooms. Predictive models are needed to better manage and mitigate loss of fish habitat due to these changes. The U.S. Geological Survey (USGS), in cooperation with the Minnesota Department of Natural Resources, developed predictive water-quality models for two agricultural land-use dominated lakes in Minnesota—Madison Lake and Pearl Lake, which are part of Minnesota’s sentinel lakes monitoring program—to assess algal community dynamics, water quality, and fish habitat suitability of these two lakes under recent (2014) meteorological conditions. The interaction of basin processes to these two lakes, through the delivery of nutrient loads, were simulated using CE-QUAL-W2, a carbon-based, laterally averaged, two-dimensional water-quality model that predicts distribution of temperature and oxygen from interactions between nutrient cycling, primary production, and trophic dynamics. The CE-QUAL-W2 models successfully predicted water temperature and dissolved oxygen on the basis of the two metrics of mean absolute error and root mean square error. For Madison Lake, the mean absolute error and root mean square error were 0.53 and 0.68 degree Celsius, respectively, for the vertical temperature profile comparisons; for Pearl Lake, the mean absolute error and root mean square error were 0.71 and 0.95 degree Celsius, respectively, for the vertical temperature profile comparisons. Temperature and dissolved oxygen were key metrics for calibration targets. These calibrated lake models also simulated algal community dynamics and water quality. The model simulations presented potential explanations for persistently large total phosphorus concentrations in Madison Lake, key differences in nutrient concentrations between these lakes, and summer blue-green algal bloom persistence. Fish habitat suitability simulations for cool-water and warm-water fish indicated that, in general, both lakes contained a large proportion of good-growth habitat and a sustained period of optimal growth habitat in the summer, without any periods of lethal oxythermal habitat. For Madison and Pearl Lakes, examples of important cool-water fish, particularly game fish, include northern pike ( Esox lucius ), walleye ( Sander vitreus ), and black crappie ( Pomoxis nigromaculatus ); examples of important warm-water fish include bluegill ( Lepomis macrochirus ), largemouth bass ( Micropterus salmoides ), and smallmouth bass ( Micropterus dolomieu ). Sensitivity analyses were completed to understand lake response effects through the use of controlled departures on certain calibrated model parameters and input nutrient loads. These sensitivity analyses also operated as land-use change scenarios because alterations in agricultural practices, for example, could potentially increase or decrease nutrient loads.

Minnesota↗

USGS2009

Hydrologic characterization for Spring Creek and hydrologic budget and model scenarios for Sheridan Lake, South Dakota, 1962-2007

The U.S. Geological Survey cooperated with South Dakota Game, Fish and Parks to characterize hydrologic information relevant to management of water resources associated with Sheridan Lake, which is formed by a dam on Spring Creek. This effort consisted primarily of characterization of hydrologic data for a base period of 1962 through 2006, development of a hydrologic budget for Sheridan Lake for this timeframe, and development of an associated model for simulation of storage deficits and drawdown in Sheridan Lake for hypothetical release scenarios from the lake. Historically, the dam has been operated primarily as a 'pass-through' system, in which unregulated outflows pass over the spillway; however, the dam recently was retrofitted with an improved control valve system that would allow controlled releases of about 7 cubic feet per second (ft3/s) or less from a fixed depth of about 60 feet (ft). Development of a hydrologic budget for Sheridan Lake involved compilation, estimation, and characterization of data sets for streamflow, precipitation, and evaporation. The most critical data need was for extrapolation of available short-term streamflow records for Spring Creek to be used as the long-term inflow to Sheridan Lake. Available short-term records for water years (WY) 1991-2004 for a gaging station upstream from Sheridan Lake were extrapolated to WY 1962-2006 on the basis of correlations with streamflow records for a downstream station and for stations located along two adjacent streams. Comparisons of data for the two streamflow-gaging stations along Spring Creek indicated that tributary inflow is approximately proportional to the intervening drainage area, which was used as a means of estimating tributary inflow for the hydrologic budget. Analysis of evaporation data shows that sustained daily rates may exceed maximum monthly rates by a factor of about two. A long-term (1962-2006) hydrologic budget was developed for computation of reservoir outflow from Sheridan Lake for the historical pass-through operating system. Two inflow components (stream inflow and precipitation) and one outflow component (evaporation) were considered. The hydrologic budget uses monthly time steps within a computational year that includes two 6-month periods - May through October, for which evaporation is accounted for, and November through April, when evaporation is considered negligible. Results indicate that monthly evaporation rates can substantially exceed inflow during low-flow periods, and potential exists for outflows to begin approaching zero-flow conditions substantially prior to the onset of zero-inflow conditions, especially when daily inflow and evaporation are considered. Results also indicate that September may be the month for greatest potential benefit for enhancing fish habitat and other ecosystem values in downstream reaches of Spring Creek with managed releases of cool water. Computed monthly outflows from Sheridan Lake for September are less than 1.0 ft3/s for 8 of the 44 years (18 percent) and are less than 2.0 ft3/s for 14 of the 44 years (32 percent). Conversely, none of the computed outflows for May are less than 2.0 ft3/s. A short-term (July through September 2007) data set was used to calculate daily evaporation from Sheridan Lake and to evaluate the applicability of published pan coefficients. Computed values of pan coefficients of approximately 1.0 and 1.1 for two low-flow periods are larger than the mean annual pan coefficient of 0.74 for the area that is reported in the literature; however, the computed values are consistent with pan coefficients reported elsewhere for similar late summer and early fall periods. Thus, these results supported the use of variable monthly pan coefficients for the long-term hydrologic budget. A hydrologic model was developed using the primary components of the hydrologic budget and was used to simulate monthly storage deficits and drawdown for Sheridan Lake using hypothetical

South Dakota↗

USGS2011

White-tailed deer age ratios as herd management and predator impact measures in Pennsylvania

A review of the Pennsylvania Game Commission's (PGC) deer management program and public concern about predator impacts on deer ( Odocoileus virginianus ) populations compelled the PGC to investigate the role of age ratios in developing management recommendations. Age ratios, such as proportion of juveniles in the antlerless harvest, may provide an index to population productivity and predator impacts. We estimated proportion of juveniles in the antlerless harvest from hunter-killed deer, population trends using the Pennsylvania (USA) sex–age–kill model, and reproduction from road-killed females. Using these estimates and a simulation model, we concluded that no single age-ratio value would serve as a reliable measure of population status. Wildlife Management Unit-specific trends in proportion of juveniles in the antlerless harvest and population trends provided the most relevant management information. We also provide an example decision chart to guide management actions in response to declining age ratios in the harvest. Although predator management activities and juvenile survival studies are often desired by the public, our decision-chart example indicated a number of deer management options exist before investing resources in predator management activities and juvenile survival studies.

Pennsylvannia↗

USGS2019

Water-balance modeling of selected lakes for evaluating viability as long-term fisheries in Kidder, Logan, and Stutsman Counties, North Dakota

Water levels in lakes and wetlands in the central North Dakota Missouri Coteau region that were either dry or only sporadically held water since before the 1930s have been rising since the early 1990s in response to an extended wet period. The lakes have remained full since the mid-1990s, which has provided benefits to migratory waterfowl, fisheries, and wildlife. A small shift in climate conditions, either to drier or wetter conditions, can have a large effect on the lake levels of these water bodies. The North Dakota Game and Fish Department identified five lakes as candidates for sustaining long-term fisheries. The lakes are in Kidder, Stutsman, and Logan Counties, and some lakes might receive inflow from mostly freshwater aquifers, such as the Central Dakota and Streeter aquifers, and were mostly dry during the early 1990s. After about 1995, the lakes had filled up and were deep enough to sustain populations of game fish such as walleye, perch, and northern pike. Before investing in development of permanent fisheries and associated infrastructure, such as campgrounds and boat ramps, fisheries biologists needed to know if the lake levels are likely to remain high in coming decades. The U.S. Geological Survey, in cooperation with the North Dakota Game and Fish Department, developed a water-balance model to determine the effects of precipitation, evapotranspiration, and groundwater interaction on lake volumes. The model was developed using climate input data and lake volumes for the calibration period 1992 through 2016, during which historical lake volumes could be estimated using land surface elevation data and Landsat images. Long-term (1940–2018) climate input data were used with the water-balance model to reconstruct historical lake volumes prior to the calibration period, and block-bootstrapping was used to simulate potential future climate input data and lake volumes for 2017 through 2067. The simulated future lake volumes were used to estimate the likelihood of annual lake volumes remaining consistent, increasing, or decreasing through the year 2067. Of the five lakes, Sibley Lake was the most likely to sustain a long-term fishery for a period longer than 50 years. The simulated lake volumes for Alkaline Lake, Big Mallard Marsh, and Remmick Lake indicated the lakes have a 50-percent chance to fall below 75 percent of their 2016 volume by about 2030, 2067, and 2025, respectively. Simulation results for Marvin Miller Lake were substantially different compared to the other four lakes and indicated the lake has a 50-percent chance to fall below 75 percent of its 2016 volume prior to 2025.

North Dakota↗

USGS2010

Preliminary assessment of factors influencing riverine fish communities in Massachusetts

The U.S. Geological Survey, in cooperation with the Massachusetts Department of Conservation and Recreation (MDCR), Massachusetts Department of Environmental Protection (MDEP), and the Massachusetts Department of Fish and Game (MDFG), conducted a preliminary investigation of fish communities in small- to medium-sized Massachusetts streams. The objective of this investigation was to determine relations between fish-community characteristics and anthropogenic alteration, including flow alteration and impervious cover, relative to the effect of physical basin and land-cover (environmental) characteristics. Fish data were obtained for 756 fish-sampling sites from the Massachusetts Division of Fisheries and Wildlife fish-community database. A review of the literature was used to select a set of fish metrics responsive to flow alteration. Fish metrics tested include two fish-community metrics (fluvial-fish relative abundance and fluvial-fish species richness), and five indicator species metrics (relative abundance of brook trout, blacknose dace, fallfish, white sucker, and redfin pickerel). Streamflows were simulated for each fish-sampling site using the Sustainable Yield Estimator application (SYE). Daily streamflows and the SYE water-use database were used to determine a set of indicators of flow alteration, including percent alteration of August median flow, water-use intensity, and withdrawal and return-flow fraction. The contributing areas to the fish-sampling sites were delineated and used with a Geographic Information System (GIS) to determine a set of environmental characteristics, including elevation, basin slope, percent sand and gravel, percent wetland, and percent open water, and a set of anthropogenic-alteration variables, including impervious cover and dam density. Two analytical techniques, quantile regression and generalized linear modeling, were applied to determine the association between fish-response variables and the selected environmental and anthropogenic explanatory variables. Quantile regression indicated that flow alteration and impervious cover were negatively associated with both fluvial-fish relative abundance and fluvial-fish species richness. Three generalized linear models (GLMs) were developed to quantify the response of fish communities to multiple environmental and anthropogenic variables. Flow-alteration variables are statistically significant for the fluvial-fish relative-abundance model. Impervious cover is statistically significant for the fluvial-fish relative-abundance, fluvial-fish species richness, and brook trout relative-abundance models. The variables in the equations were demonstrated to be significant, and the variability explained by the models, as measured by the correlation between observed and predicted values, ranges from 39 to 65 percent. The GLM models indicated that, keeping all other variables the same, a one-unit (1 percent) increase in the percent depletion or percent surcharging of August median flow would result in a 0.4-percent decrease in the relative abundance (in counts per hour) of fluvial fish and that the relative abundance of fluvial fish was expected to be about 55 percent lower in net-depleted streams than in net-surcharged streams. The GLM models also indicated that a unit increase in impervious cover resulted in a 5.5-percent decrease in the relative abundance of fluvial fish and a 2.5-percent decrease in fluvial-fish species richness.

Massachusetts↗

USGS2019

Multi-species duck harvesting using dynamic programming and multi-criteria decision analysis

1.Multiple species are often exposed to a common hunting season, but harvest and population objectives may not be fully achieved if harvest potential varies among species and/or species abundances are not correlated through time. Our goal was to develop an approach for setting a common hunting season that would recognize heterogeneity in species productivity and would select annual hunting seasons conditioned on the status of individual species. 2.We first used stochastic dynamic programming to generate optimal, state‐dependent harvest strategies for 18 candidate regulatory scenarios. We simulated the performance of these strategies, and then used multi‐criteria decision making to identify preferred regulatory scenarios for duck hunting seasons in the Atlantic Flyway of the U.S. 3.Generally, estimates of annual population size were not correlated among species. Mallards had the highest estimated intrinsic rate of growth, green‐winged teal, wood ducks, and ring‐necked ducks had intermediate values, and goldeneyes were the least productive. Estimated carrying capacity was highest for mallards and lowest for green‐winged teal. 4.Managers had greatest interest in maximizing season length (33%) and aggregate duck abundance (28%), and less interest in maximizing aggregate harvest (19%) and the number of years between a change in hunting season regulations (19%). Several regulatory scenarios provided acceptable tradeoffs among these objectives. 5.Synthesis and applications. Separate hunting seasons for various species of game may be untenable, either due to the added cost and regulatory complexity, or because selective harvesting of stocks may be difficult due to problems in species identification. Rather than averaging species‐specific productivities, or basing hunting seasons on the least (or most) productive species, we describe an approach in which productivity and annual population status are considered explicitly for each species. By combining stochastic dynamic programming with multi‐criteria decision analysis, we can identify a regulatory strategy that can address a diverse set of objectives and explicitly recognize the tradeoffs among them. To meet the Atlantic Flyway's objectives as identified by waterfowl managers, our results suggest a regulatory strategy in which the harvest is targeted at 98% of aggregate maximum sustainable yield, most emphasis is placed on accumulating harvests of mallards and wood ducks, and by using a set of regulatory options that are more conservative than those currently in use.

Journal of Applied Ecology↗

USGS2014

A test of the compensatory mortality hypothesis in mountain lions: a management experiment in West-Central Montana

Mountain lions ( Puma concolor ) are widely hunted for recreation, population control, and to reduce conflict with humans, but much is still unknown regarding the effects of harvest on mountain lion population dynamics. Whether human hunting mortality on mountain lions is additive or compensatory is debated. Our primary objective was to investigate population effects of harvest on mountain lions. We addressed this objective with a management experiment of 3 years of intensive harvest followed by a 6-year recovery period. In December 2000, after 3 years of hunting, approximately 66% of a single game management unit within the Blackfoot River watershed in Montana was closed to lion hunting, effectively creating a refuge representing approximately 12% (915 km 2 ) of the total study area (7,908 km 2 ). Hunting continued in the remainder of the study area, but harvest levels declined from approximately 9/1,000 km 2 in 2001 to 2/1,000 km 2 in 2006 as a result of the protected area and reduced quotas outside. We radiocollared 117 mountain lions from 1998 to 2006. We recorded known fates for 63 animals, and right-censored the remainder. Although hunting directly reduced survival, parameters such as litter size, birth interval, maternity, age at dispersal, and age of first reproduction were not significantly affected. Sensitivity analysis showed that female survival and maternity were most influential on population growth. Life-stage simulation analysis (LSA) demonstrated the effect of hunting on the population dynamics of mountain lions. In our non-hunted population, reproduction (kitten survival and maternity) accounted for approximately 62% of the variation in growth rate, whereas adult female survival accounted for 30%. Hunting reversed this, increasing the reliance of population growth on adult female survival (45% of the variation in population growth), and away from reproduction (12%). Our research showed that harvest at the levels implemented in this study did not affect population productivity (i.e., maternity), but had an additive effect on mountain lion mortality, and therefore population growth. Through harvest, wildlife managers have the ability to control mountain lion populations.

Montana↗

USGS2013

Interactions between striped bass and other gamefish in reservoirs

Competitive interactions among reservoir fishes may be pronounced because fish assemblages in these artificial environments have had little time to develop niche-partitioning strategies that alleviate negative interspecific interactions. Such interactions may at times have been intensified by introductions of predators such as striped bass Morone saxatilis , introduced to create additional fisheries and control pelagic clupeids. Possible interactions between existing fish assemblages and striped bass include predation and competition. While there is a perception among angler groups that predation by striped bass on co-existing game fish is significant, most studies have reported little or no predation on game fish my striped bass and have considered predation rare and inconsequential. Moreover, predation that occurs will likely be compensatory and fail to reduce overall game fish survival. Any indirect effect of striped bass predation by restricting prey-sized game fish to limited refuge sites remains unknown. Exploitative competition may be more common. Although infrequently, introduced striped bass have depleted prey resources shared with other piscivores, particularly when stocking rates have been high, when there is a high rate of natural reproduction, or when prey supply has plunged in response to environmental fluxes. Fluctuation in prey supply, associated with ordinary environmental variability, and associated time lages in prey supply and predator demand, preclude adjusting predator densities to exactly balance demand with supply. The frequency of low supply-demand rations varies across systems and exhibits seasonal trends. Nevertheless, chronic supply-demand imbalances are manageable where the predator assemblage is at least partially controlled through stocking, harvest regulations, or both. Because of the poor state of knowledge concerning the parameters defining balance and because uncontrollable annual fluctuations preclude exact management of alternating prey levels, we suggest adjusting stocking to manage demand to that it equals the median historical prey supply. Simulating the removal of striped bass and predicting the aftermath may be the most cost-efficient way to provide decision support for stakeholders involved in determining if a striped bass stocking program is beneficial to most users.

American Fisheries Society Symposium↗

USGS2011

Factors influencing riverine fish assemblages in Massachusetts

The U.S. Geological Survey, in cooperation with the Massachusetts Department of Conservation and Recreation, Massachusetts Department of Environmental Protection, and the Massachusetts Department of Fish and Game, conducted an investigation of fish assemblages in small- to medium-sized Massachusetts streams. The objective of this study was to determine relations between fish-assemblage characteristics and anthropogenic factors, including impervious cover and estimated flow alteration, relative to the effects of environmental factors, including physical-basin characteristics and land use. The results of this investigation supersede those of a preliminary analysis published in 2010. Fish data were obtained for 669 fish-sampling sites from the Massachusetts Division of Fisheries and Wildlife fish-community database. A review of the literature was used to select fish metrics - species richness, abundance of individual species, and abundances of species grouped on life history traits - responsive to flow alteration. The contributing areas to the fish-sampling sites were delineated and used with a geographic information system to determine a set of environmental and anthropogenic factors that were tested for use as explanatory variables in regression models. Reported and estimated withdrawals and return flows were used together with simulated unaltered streamflows to estimate altered streamflows and indicators of flow alteration for each fish-sampling site. Altered streamflows and indicators of flow alteration were calculated on the basis of methods developed in a previous U.S. Geological Survey study in which unaltered daily streamflows were simulated for a 44-year period (water years 1961-2004), and streamflow alterations were estimated by use of water-withdrawal and wastewater-return data previously reported to the State for the 2000-04 period and estimated domestic-well withdrawals and septic-system discharges. A variable selection process, conducted using principal components analysis and Spearman rank correlation, was used to select a set of 15 non-redundant environmental and anthropogenic factors to test for use as explanatory variables in the regression analyses. Twenty-one fish species were used in a multivariate analysis of fish-assemblage patterns. Results of nonmetric multidimensional scaling and hierarchical cluster analysis were used to group fish species into fluvial and macrohabitat generalist habitat-use classes. Two analytical techniques, quantile regression and generalized linear modeling, were applied to characterize the association between fish-response variables and environmental and anthropogenic explanatory variables. Quantile regression demonstrated that as percent impervious cover and an indicator of percent alteration of August median flow from groundwater withdrawals increase, the relative abundance and species richness of fluvial fish decrease. The quantile regression plots indicate that (1) as many as seven fluvial fish species are expected in streams with little flow alteration or impervious cover, (2) no more than four fluvial fish species are expected in streams where flow alterations from groundwater withdrawals exceed 50 percent of the August median flow or the percent area of impervious cover exceeds 15 percent, and (3) few fluvial fish remain at high rates of withdrawal (approaching 100 percent) or high rates of impervious cover (between 25 and 30 percent). Three generalized linear models (GLMs) were developed to quantify the response of fluvial fish to multiple environmental and anthropogenic variables. All variables in the GLM equations were demonstrated to be significant (p less than 0.05, with most less than 0.01). Variables in the fluvial-fish relative-abundance model were channel slope, estimated percent alteration of August median flow from groundwater withdrawals, percent wetland in a 240-meter buffer strip, and percent impervious cover. Variables in the fluvial-fish species-richness model were drainage area, channel slope, total undammed reach length, percent wetland in a 240-meter buffer strip, and percent impervious cover. Variables in the brook trout relativeabundance model were drainage area, percent open water, and percent impervious cover. The variability explained by the GLM models, as measured by the pseudo R2, ranged from 18.2 to 34.6, and correlations between observed and predicted values ranged from 0.50 to 0.60. Results of GLM models indicated that, keeping all other variables the same, a one-unit (1 percent) increase in the percent depletion of August median flow would result in a 0.9-percent decrease in the relative abundance (in counts per hour) of fluvial fish. The results of GLM models also indicated that a unit increase in impervious cover (1 percent) resulted in a 3.7-percent decrease in the relative abundance of fluvial fish, a 5.4-percent decrease in fluvial-fish species richness, and an 8.7-percent decrease in brook trout relative abundance.

Massachusetts↗

USGS2024

Wolf harvest management strategy evaluation: Annual Report, 2024

Wolf harvest season setting is complicated and controversial. State law requires Montana Fish, Wildlife and Parks (MFWP) to both reduce the wolf population and avoid federal relisting under the Endangered Species Act (Montana Fish, Wildlife and Parks, 2002). Disparate stakeholder groups each have different objectives for wolf management. For instance, big game advocates want to see improved big game populations and hunting opportunities in northwest Montana, while wolf advocates want to see regulations that minimize wolf mortality. Decision making about season setting tries to balance these objectives. Wolf hunting and trapping season decisions are made by the Montana Fish and Wildlife Commission and are informed by annual wolf abundance estimates from an integrated patch occupancy model (iPOM, Sells et al., 2022c) as well as the predictions of wolf abundance into the future under potential constant harvest levels. Parametric uncertainty (uncertainty surrounding the value of a parameter) from the iPOM estimates is propagated through to future projections, providing the Commission with plausible and worst-case outcomes of different levels of public harvest over the short term, i.e., five years into the future, on the wolf population in Montana (Parks et al., 2024). An alternative approach to inform wolf management and harvest decisions is through adaptive management. Adaptative management is appropriate for decisions that are made iteratively and when monitoring data are collected to learn about the outcomes from decisions, where monitoring data help to reduce critical uncertainties regarding ecosystem function or management outcomes (Walters, 1986; Williams, 2011). Management strategy evaluation (MSE) is one way to develop an adaptive management framework. MSE was developed by fisheries managers and scientists to more accurately and fully incorporate various forms of uncertainty, consider long-term time horizons, and add more transparency in a fisheries context (Punt et al., 2016). It has been used routinely and has become a standard approach for complicated and contentious marine fisheries management situations, yet it has been underutilized in wildlife management (but see Bunnefeld et al., 2013, 2011). MSE is a forward simulation approach for testing prospective management options or strategies over a wide range of possible states (Punt et al., 2016). A MSE framework captures the ‘truth’ or what is happening in the system (termed the ‘the operating model’) and the information available to the decision makers (termed ‘the estimation model’ or ‘management strategy’). More precisely, there are four main processes modeled. First, models are constructed based on current understanding and data to represent ‘truth’. Second, the collection of monitoring data is simulated from the ‘truth’ model. Third, the simulated monitoring data are fit to an estimation model and the next time step’s population metrics are predicted from the estimated parameters. Fourth, based on the estimation model results and the predictions, the decision-making process is simulated following a management strategy, whereby a decision is made and the implementation of this decision feeds back into the ‘truth’ model (Figure 1). This process continues through time. Additionally, each simulation through time is repeated to capture the full range of stochasticity and uncertainty.

Montana↗

USGS2020

Hydrology of Haskell Lake and investigation of a groundwater contamination plume, Lac du Flambeau Reservation, Wisconsin

Haskell Lake is a shallow, 89-acre drainage lake in the headwaters of the Squirrel River, on the Lac du Flambeau Reservation in northern Wisconsin. The lake has long been valued by the Lac du Flambeau Band of Lake Superior Chippewa Indians (LDF Tribe) for abundant wild rice and game fish. In recent decades, however, wild rice has mostly disappeared from the lake and the fishery has declined. A petroleum contamination plume discovered in the 1990s in the shallow aquifer upgradient from the northern end of the lake poses a threat to the ecological health of the lake and the aquifer, which is the sole drinking water source for nearby residents and businesses. Understanding of the lake’s hydrology is important to the LDF Tribe as they seek to restore wild rice and maintain the ecological health of the Haskell Lake/Tower Creek watershed. An improved understanding of lithology in the area of the contamination plume, documentation of a contamination pathway from groundwater in the plume source area to Haskell Lake, and an understanding of the plume extent beneath the lake are needed to advance remediation efforts. Evaluation of the fraction of groundwater discharge that is contaminated relative to the overall lake water budget is desired as a first step towards determining the extent of ecological effects from the plume. A cooperative study between the U.S. Geological Survey and the LDF Tribe was initiated to quantify the lake water budget and the sources of water to the lake, to provide a rough estimate of the maximum quantity of groundwater discharge to the lake that may be contaminated, and to improve the conceptual understanding of the plume extent and subsurface materials in the area of contamination. The results of this study can help inform natural resource management of the Haskell Lake/Tower Creek watershed, including planned wild rice restoration and cleanup of the contaminant plume. During 2016–17, field data on lake and groundwater levels, gradients, fluxes, and subsurface lithology were collected using a variety of techniques that ranged from basic measurement of water levels and streamflows to distributed temperature sensing, vertical temperature profiling, and several shallow geophysical methods. The data were used to inform a MODFLOW–NWT model that simulated the contributing groundwatershed, including the water budget for Haskell Lake and Tower Creek using the Lake, Streamflow-Routing, and Unsaturated Zone-Flow Packages. Particle tracking with the MODFLOW solution (using MODPATH 6) was used to improve understanding of the downgradient extent of the contamination plume, estimate groundwater flux through the plume area, and delineate the groundwater contributing area (groundwatershed) for the lake/creek system. Linear uncertainty estimates for model results were computed during model parameter estimation using the software package PEST++. Results indicate groundwater discharge along the perimeter of Haskell Lake, with groundwater accounting for about 22 (± 11.5) percent of the lake water budget. Field data and particle tracking results indicate discharge of the entire contamination plume to Haskell Lake. Although the exact locations where contaminated groundwater enters the lake are unknown, the downgradient extent of the plume beneath Haskell Lake is likely limited to within about 700 feet from the shore. Groundwater flux through the plume accounts for at most about 1.4 percent of total groundwater discharge to Haskell Lake, or about 0.3 percent of the lake water budget. Most groundwater discharging to Haskell Lake and Tower Creek originates as terrestrial recharge. A lesser amount originates in or passes through neighboring lakes, including Buckskin, Crawling Stone, Broken Bow, Tippecanoe, and Jerms Lakes, as well as several unnamed kettles. The average age of simulated groundwater discharge to the lake is about 20 years.

Wisconsin↗