SEARCH · USGS Science

Results for “Computer contribution”

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.

At least 469 records · Page 26Linked to original sources

USGS1979

Storm and flood of July 31-August 1, 1976, in the Big Thompson River and Cache la Poudre River basins, Larimer and Weld Counties, Colorado

PART A: Devastating flash floods swept through the canyon section of Larimer County in north-central Colorado during the night of July 31-August I, 1976, causing 139 deaths, 5 missing persons, and more than $35 million in total damages. The brunt of the storms occurred over the Big Thompson River basin between Drake and Estes Park with rainfall amounts as much as 12 inches being reported during the storm period. In the Cache la Poudre River basin to the north, a rainfall amount of 10 inches was reported for one locality while 6 inches fell over a widespread area near the central part of the basin. The storms developed when strong low-level easterly winds to the rear of a polar front pushed a moist, conditionally unstable airmass upslope into the Front Range of the Rocky Mountains. Orographic uplift released the convective instability, and light south-southeasterly winds at middle and upper levels allowed the storm complex to remain nearly stationary over the foothills for several hours. Minimal entrainment of relatively moist air at middle and upper levels, very low cloud bases, and a slightly tilted updraft structure contributed to a high precipitation efficiency. Intense rainfall began soon after 1900 MDT (Mountain Daylight Time) in the Big Thompson River and the North Fork Cache la Poudre River basins. A cumulative rainfall curve developed for Glen Comfort from radar data indicates that 7.5 inches of rain fell during the period 1930-2040 MDT on July 31. In the central part of the storm area west of Fort Collins, the heaviest rainfall began about 2200 MDT on July 31 and continued until 0100 MDT on August 1. Peak discharges were extremely large on many streams in the storm area-exceeding previously recorded maximum discharges at several locations. The peak discharge of the Big Thompson River at the gaging station at the canyon mouth, near Drake was 31,200 cubic feet per second or more than four times the previous maximum discharge of 7,600 cubic feet per second at the site during 88 years of flood history. At the gaging station on the North Fork Big Thompson River at Drake, the peak discharge on July 31 was 8,710 cubic feet per second as compared to the previous maximum discharge during 29 years of record of 1,290 cubic feet per second. Peak discharges for three small tributaries near the area of heaviest rainfall northeast of Estes Park exceeded previously recorded maximum discharges for basins of less than 4 square miles in Colorado. Stream velocities were rapid along the tributaries near the storm center and on the Big Thompson River in the canyon section, with average velocities of 20-25 feet per second being common. The flood crest on the Big Thompson River moved through the 7.7-mile reach between Drake and the canyon mouth in about 30 minutes for an average travel rate of 15 miles per hour, or about 23 feet per second. The peak discharge of the flood on the Big Thompson River at the canyon mouth exceeded the 100-year flood discharge for the site by a ratio of 1.8. Upstream in the Big Thompson River basin, the flood was even more rare being 3.8 times the estimated 100-year flood discharge at the site on the Big Thompson River just upstream from Drake. In the Cache la Poudre River basin, recurrence intervals were computed to be 100 years for the flood on Deadman Creek and 16 years for Rist Canyon and the Cache la Poudre River at the canyon mouth near Fort Collins. Although the rainfall and flood discharges were unusually large, they are not unprecedented for some areas along the eastern foothills and plains of Colorado. The May 1935 and June 1965 floods on some streams along the eastern plains greatly exceeded the 1976 flood peaks in the storm area. Prior floods on several other streams in the foothills have approximately equaled the 1976 peak discharges. PART B: Intense rainfall from the Big Thompson thunderstorm complex on the evening of July 31,1976, and the ensuing floods that evening and the fol

Professional Paper↗

USGS2025

Host responses and viral traits interact to shape the impacts of climate warming on highly pathogenic avian influenza in migratory waterfowl

Emerging infectious diseases pose threats to wildlife populations, as exemplified by recent outbreaks of avian influenza viruses in wild birds. Climate change can affect infection dynamics in wildlife through direct effects on pathogens (e.g., environmental decay rates) and changes to host ecology, including shifting migration patterns. Here, we adapt an existing mechanistic model that couples migration and infection to study how traits of highly pathogenic avian influenza (HPAI) viruses contribute to HPAI outcomes in migratory waterfowl, then apply this model to explore potential impacts of climate change on HPAI dynamics. We find that the simulated impacts of HPAI on the host population under baseline climate conditions varied from no impact to 100% mortality, depending on viral traits. In most cases, traits related to transmission (i.e., contact rates, shedding rates) were more important for HPAI establishment probability, infection prevalence, and mortality than were other viral traits (e.g., environmental temperature sensitivity, cross-protective immunity). We then simulated the effects of climate change (i.e., altered temperature regimes) on HPAI dynamics both via viral environmental decay and via changes in bird migration phenology. In these simulations, we found that a 9-day advancement in spring migration timing increased the duration of HPAI outbreaks by increasing time birds spent at their breeding grounds, leading to higher mortality and fewer infections. In contrast, increased viral decay in warmer years had a smaller, but opposite impact. These patterns depended on the primary transmission mode of HPAI (i.e., direct vs. environmental) and its sensitivity to environmental temperatures. Together, these results suggest that climate change is likely to increase the impacts of HPAI on waterfowl populations if HPAI relies strongly on direct transmission and birds advance their spring migration. Further integrating host-viral co-evolution and other climatic changes (e.g., salinity, humidity) could provide more precise predictions of how HPAI dynamics could change in the future.

Alaska, California, Oregon↗

USGS2010

Methods for estimating flow-duration and annual mean-flow statistics for ungaged streams in Oklahoma

Flow statistics can be used to provide decision makers with surface-water information needed for activities such as water-supply permitting, flow regulation, and other water rights issues. Flow statistics could be needed at any location along a stream. Most often, streamflow statistics are needed at ungaged sites, where no flow data are available to compute the statistics. Methods are presented in this report for estimating flow-duration and annual mean-flow statistics for ungaged streams in Oklahoma. Flow statistics included the (1) annual (period of record), (2) seasonal (summer-autumn and winter-spring), and (3) 12 monthly duration statistics, including the 20th, 50th, 80th, 90th, and 95th percentile flow exceedances, and the annual mean-flow (mean of daily flows for the period of record). Flow statistics were calculated from daily streamflow information collected from 235 streamflow-gaging stations throughout Oklahoma and areas in adjacent states. A drainage-area ratio method is the preferred method for estimating flow statistics at an ungaged location that is on a stream near a gage. The method generally is reliable only if the drainage-area ratio of the two sites is between 0.5 and 1.5. Regression equations that relate flow statistics to drainage-basin characteristics were developed for the purpose of estimating selected flow-duration and annual mean-flow statistics for ungaged streams that are not near gaging stations on the same stream. Regression equations were developed from flow statistics and drainage-basin characteristics for 113 unregulated gaging stations. Separate regression equations were developed by using U.S. Geological Survey streamflow-gaging stations in regions with similar drainage-basin characteristics. These equations can increase the accuracy of regression equations used for estimating flow-duration and annual mean-flow statistics at ungaged stream locations in Oklahoma. Streamflow-gaging stations were grouped by selected drainage-basin characteristics by using a k-means cluster analysis. Three regions were identified for Oklahoma on the basis of the clustering of gaging stations and a manual delineation of distinguishable hydrologic and geologic boundaries: Region 1 (western Oklahoma excluding the Oklahoma and Texas Panhandles), Region 2 (north- and south-central Oklahoma), and Region 3 (eastern and central Oklahoma). A total of 228 regression equations (225 flow-duration regressions and three annual mean-flow regressions) were developed using ordinary least-squares and left-censored (Tobit) multiple-regression techniques. These equations can be used to estimate 75 flow-duration statistics and annual mean-flow for ungaged streams in the three regions. Drainage-basin characteristics that were statistically significant independent variables in the regression analyses were (1) contributing drainage area; (2) station elevation; (3) mean drainage-basin elevation; (4) channel slope; (5) percentage of forested canopy; (6) mean drainage-basin hillslope; (7) soil permeability; and (8) mean annual, seasonal, and monthly precipitation. The accuracy of flow-duration regression equations generally decreased from high-flow exceedance (low-exceedance probability) to low-flow exceedance (high-exceedance probability) . This decrease may have happened because a greater uncertainty exists for low-flow estimates and low-flow is largely affected by localized geology that was not quantified by the drainage-basin characteristics selected. The standard errors of estimate of regression equations for Region 1 (western Oklahoma) were substantially larger than those standard errors for other regions, especially for low-flow exceedances. These errors may be a result of greater variability in low flow because of increased irrigation activities in this region. Regression equations may not be reliable for sites where the drainage-basin characteristics are outside the range of values of independent vari

Scientific Investigations Report↗

USGS2019

Comparison of groundwater-model construction methods, representations of glacial geology, model designs, and groundwater-model flow simulations within Elkhart County, Indiana

Automated data-processing methods allow hydrologists to efficiently incorporate digital well-record datasets into the construction of hydrostratigraphic frameworks for groundwater-flow models. The method selected to construct the hydrostratigraphic framework can affect the extent of geologic heterogeneity that can be included in the model. The detail generated from a hydrostratigraphic framework can affect groundwater simulation results. The effects of detail on model accuracy, groundwater-flow simulations, and particle-tracking simulations are described in this study. This report compares differences in hydrostratigraphic frameworks and results of groundwater models using (1) a method that incorporates more hydrologic judgment at the expense of using limited lithologic data and (2) a method that is more automated and uses all available lithologic data. The study additionally evaluates the effect of model discretization and inclusion of more (or less) geologic detail on simulation results. Two methods were used to create hydrostratigraphic frameworks of glacial deposits in the St. Joseph River Basin. One method, referred to as the subjective method, manually identifies stratigraphic boundaries using a sample of well logs from State databases and uses two-dimensional kriging to create three model layers of the study area. Indicator kriging is used to define aquifer extent in each layer. The second method, referred to as the objective method, uses three-dimensional kriging to automatically create a detailed heterogeneous model of the study area using all wells logs from the State database. The objective method increases detail in the vertical by greatly increasing the number of computer groundwater model layers from 3 to 30. In Elkhart County, Indiana, a previously published model represents the product of the subjective method, and a newly calibrated model of the same area represents the product of the objective method. An automated calibration procedure was used with the objective model (derived from the objective method) for Elkhart County. The two most-sensitive parameters for the Elkhart County objective model are horizontal hydraulic conductivity of the sand and the combined sand and gravel/gravel deposits. Vertical hydraulic conductivity of the fine-grained and intermediate-sized deposits could not be estimated, possibly indicating major flow paths are along a continuously connected series of sand and gravel deposits and not through a confining layer. The statistics measuring model calibration accuracy for the objective model were slightly better than statistics for the subjective model (model derived from the subjective method) of Elkhart County, but the hydraulic conductivities and flow rates for the two models were different. The mean absolute errors between simulated and measured groundwater levels are 2.04 and 2.16 feet for the objective and subjective models, respectively. Simulated seepage losses from and groundwater discharges to measured stream reaches in the objective model were evenly balanced in terms of over and under simulations of measured values; the subjective model tended to overpredict measured groundwater discharge to streams. The overprediction may be related to the 58 percent greater total inflow and outflow through the subjective model. The greater flow rate through the subjective model results from higher horizontal hydraulic conductivities in the subjective model than in the objective model. Horizontal hydraulic conductivity ranged from 23.9 to 111 feet per day in the objective model and generally ranged from 170 to 370 feet per day in the subjective model. The improvement in calibration statistics for the objective model relative to the subjective model may be from increased detail in how the objective model represents the distribution of fine- and coarse-grained deposits. The improvement also could be associated with the difference in methods used to represent the continuity of the confining unit. The effect of differences in horizontal hydraulic conductivity distributions between the two models for Elkhart County is evident in the groundwater-flow paths simulated by the objective and subjective models. At a withdrawal well location, the flow lines produced by the objective model indicate a wider contributing area than that for the subjective model. The discontinuous confining unit represented in the objective model provided the opportunity for groundwater flow to split into an upper and lower path. The split in flow simulated by the objective model at one location was independently supported by bromide concentrations in groundwater; the subjective model did not duplicate the split in flow.

Indiana↗

USGS1995

Simulated ground-water flow and sources of water in the Killbuck Creek Valley near Wooster, Wayne County, Ohio

The stratified-drift aquifer in the 3,000-ft (feet)-wide and 100-ft-deep buried valley of Killbuck Creek near Wooster in northeastern Ohio was studied. The stratified drift with adjacent sandstone and shale bedrock produce a system of ground-water flow representative of the western part of the glaciated north-eastern United States. The stratified-drift aquifer is an excellent source of water for municipal and industrial wells. The aquifer is recharged locally by water from precipitation on the valley floor and uplands, by infiltration from streams, and by lateral flow to the valley from the uplands. As a result, the aquifer is vulnerable to surface or subsurface spills of contaminants in the valley or the adjacent uplands. Quality of water in the stratified drift is affected by influx of water from bedrock lateral to or beneath the valley. This influx is controlled, in part, by the pumping stress placed on the stratified-drift aquifer. Hydrogeologic and aqueous-geochemical data were analyzed to establish the framework necessary for stead-state and transient simulations of ground-water flow in stratified drift and bedrock with a three-layer ground-water-flow model. A new model routine, the Variable-Recharge procedure, was developed to simulate areal recharge and the contribution of the uplands to the drift system. This procedure allows for water applied to land surface to infiltrate or to be rejected. Rejected recharge and ground water discharged when the water table is at land surface form surface runoff-this excess upland water can be redirected as runoff to other parts of the model. Infiltration of streamwater, areal recharge to uplands and valley, and lateral subsurface flow from the uplands to the valley are sources of water to the stratufued0druft aquifer. Water is removed from the stratified-drift aquifer at Wooster primarily by production wells pumping at a rate of approximately 8.5 ft3/s (cubic feet per second). The ground-water budget resulting from two types of simulations of ground-water flow in this study indicates the primary sources of water to the wells are recharge at or near land surface and lateral subsurface flow from the shale and sandstone bedrock. Components of recharge at land surface include induced infiltration from streams, precipitation on the valley floor, and infiltration of unchanneled upland runoff that reaches the valley floor. The steady-state simulation was designed to represent conditions during the fall of 1984. The transient simulation was designed to represent an 11-day snowmelt event, 23 February to 5 March 1985, that caused water levels to rise significantly throughout the valley. Areal recharge to the valley and flow from the uplands to the valley were determined through the Variable-Recharge procedure. The total steady-state recharge to the valley was 12.5 ft3/s. Upland sources, areal valley recharge, and induced infiltration from Killnuck Creek accounted for 63, 23, and 8 percent, respectively, of the valley recharge. An analysis of the simulated vertical flow to the buried stratified drift through surficial slit, clay, and fine sand indicates that about 75 percent of the total recharge to the buried deposits is the sum of areally extensive, relatively small flows less than about 0.01 ft? /s per model node), whereas about 25 percent of the recharge results from a really restricted, relatively large flows (greater than about 0.01 ft? /s per model node). The large-magnitude flows are located primarily beneath Clear and Little Killbuck Creeks where seepage provides abundant recharge and the surficial sediments grade into coarser alluvial-fan deposits. Chemical and isotopic studies of ground water and streamwater combined with measurements of stream infiltration provide independent support for the conclusions derived from computer simulation of ground-water flow. In addition, the chemical and isotopic studies helped quantity the rate and pathways of infiltrating water from

Water-Resources Investigations Report↗

USGS2001

Simulation of flow in the upper North Coast Limestone Aquifer, Manati-Vega Baja area, Puerto Rico

A two-dimensional computer ground-water model was constructed of the Manati-Vega Baja area to improve the understanding of the unconfined upper aquifer within the North Coast Province of Puerto Rico. The modeled area covers approximately 79 square miles within the municipios of Manati and Vega Baja and small portions of Vega Alta and Barceloneta. Steady-state two-dimensional ground-water simulations were correlated to conditions prior to construction of the Laguna Tortuguero outlet channel in 1940 and calibrated to the observed potentiometric surface in March 1995. At the regional scale, the unconfined Upper North Coast Limestone aquifer is a diffuse ground-water flow system through the Aguada and Aymamon limestone units. The calibrated model input parameters for aquifer recharge varied from 2 inches per year in coastal areas to 18 inches per year in the upland areas south of Manati and Vega Baja. The calibrated transmissivity values ranged from less than 500 feet squared per day in the upland areas near the southern boundary to 70,000 feet squared per day in the areas west of Vega Baja. Increased ground-water withdrawals from 1.0 cubic foot per second for 1940 conditions to 26.3 cubic feet per second in 1995, has reduced the natural ground-water discharge to springs and wetland areas, and induced additional recharge from the rivers. The most important regional drainage feature is Laguna Tortuguero, which is the major ground-water discharge body for the upper aquifer, and has a drainage area of approximately 17 square miles. The discharge to the sea from Laguna Tortuguero through the outlet channel has been measured on a bi-monthly basis since 1974. The outflow represents a combination of ground- and surface-water discharge over the drainage area. Hydrologic conditions, prior to construction of the Laguna Tortuguero outlet channel in 1943, can be considered natural conditions with minimal ground-water pumpage (1.0 cubic foot per second), and heads in the lagoon were 2.4 feet higher. The model was calibrated to March 1995 conditions during a dry period of minimal aquifer recharge and relatively constant water levels in the upper aquifer. For the steady-state 1995 model simulation, however, ground-water pumpage had been increased to 26.3 cubic foot per second, due to increased demand for public water supply, the heads at 0.9 feet, and the outflow to the sea at Laguna Tortuguero had been lowered considerably. Simulated ground-water inflow for 1940 hydrologic conditions included 35.9 cubic feet per second from areal recharge, contributions from streamflow along the southern boundary of 1.6 cubic feet per second, and streamflow infiltration to the upper aquifer of 4.2 cubic feet per second. Simulated ground-water outflow for 1940 hydrologic conditions are discharge to springs of 17.4 cubic feet per second, total ground-water withdrawals of 1.0 cubic feet per second, and aquifer contribution to streamflow or wetland areas of 23.4 cubic feet per second. Simulated ground-water inflow for hydrologic conditions of March 1995 include d contributions from streamflow along the southern boundary of 1.6 cubic feet per second, areal recharge of 35.9 cubic feet per second, and streamflow infiltration to the upper aquifer of 11 cubic feet per second. Simulated ground-water outflow for hydrologic conditions of March 1995 are ground-water withdrawals of 26.3 cubic feet per second, discharge from springs of 7.3 cubic feet per second, and aquifer contribution to streamflow or wetland areas of 14 .9 cubic feet per second. The overall ground-water budget increased from 41.8 cubic feet per second for 1940 conditions to 48.6 cubic feet per second for the hydrologic conditions of March 1995. The increase in ground-water budget is a direct result of increased ground-water withdrawals, which induced greater streamflow infiltration. Simulated ground-water flux to Laguna Tortuguero for 1940 conditions was 11 cubic feet per second, which drop

Water-Resources Investigations Report↗

USGS2020

Ecological status of aquatic communities in selected streams in the Milwaukee Metropolitan Sewerage District planning area of Wisconsin, 2004–13

A total of 14 wadable streams in urban or urbanizing watersheds near Milwaukee, Wisconsin, were sampled in 2004, 2007, 2010, and 2013 to assess the ecological status of aquatic communities (biota), including benthic algae and invertebrates, and fish. To assess temporal variation, additional community sampling was also done at a subset of three sites in 2011 and 2012. Relative abundances of each type of organism were used to calculate biological metrics, such as richness and diversity, percentages of intolerant and tolerant organisms, and indexes of biotic integrity for invertebrates and fish. Selected environmental (physical and chemical) data in the streams were collected to evaluate potential relations to the biota and the ecological health of the stream. Physical and chemical data included land use/land cover, stream discharge from U.S. Geological Survey (USGS) streamgages (except at 2 creeks that were not gaged), stream habitat, microhabitat at invertebrate collection locations, water quality (except at 2 creeks that were not gaged), field measurements of several water-quality constituents, measures of benthic algal biomass, and toxicity and chemical tests on extracts from passive samplers deployed at a subset of 6 sites. Relative abundances of organisms and biological metrics were compared among sampling years and with environmental metrics to evaluate the ecological status of these streams and determine primary stressors on the aquatic communities, with the aim of helping resource managers understand and work toward improving the ecological health of these and other urban and urbanizing rivers in the study area. Biological metrics for most sites indicated some level of diminished ecological status when compared across all sampled sites and when compared with rating scales for selected metrics. The least degraded sites among all those sampled—indicated by aggregate bioassessments for algae, invertebrates, and fish metrics and in order starting with the best overall condition—were the Milwaukee River near Cedarburg, Menomonee River at Menomonee Falls, Jewel Creek, and Milwaukee River at Milwaukee. The most degraded sites were Menomonee River at Wauwatosa, Root River at Greenfield, Lincoln Creek, and the Kinnickinnic River. Differences in aggregate bioassessments indicate that aquatic communities at the Menomonee River at Wauwatosa site and the Root River at Greenfield site were worse in 2013 than in 2004; however, Oak Creek and Honey Creek sites were better. In 2013, several sites had less than 30-percent pollution-sensitive diatoms indicating degraded algal assemblages. Invertebrate metrics for most of the 14 sites in 2013 were lower than in 2004 and indicate that invertebrate assemblages at most sampled sites were more degraded in 2013. Tolerant fish taxa made up more than 40 percent of assemblages at most sites and nearly 100 percent of assemblages at four sites. At times, in some smaller streams, too few fish were captured to compute an Index of Biotic Integrity with confidence, and invertebrates provided a better means for assessing the ecological status and water quality. With these few exceptions, the use of all three groups of biota provided the most robust assessments at the 14 sites in 2004–13. Physical and chemical stressors were correlated to adverse effects on aquatic biota at the sampled streams. Passive samplers were deployed at a subset of six sites in 2013. Microtox results indicated there was little or no toxicity at the Milwaukee River near Cedarburg site and at the Oak Creek site, slight toxicity at the Lincoln Creek and Honey Creek sites, and moderate toxicity at the Milwaukee River at Milwaukee site and the Little Menomonee River site; however, based on cytochrome-P450 reporter gene system toxicity tests, potential toxicity from hydrophobic organic contaminants was measured at all six sites. For all 14 sites, physical and chemical stressors related to urbanization correlated with biological metrics for algae, invertebrates, and fish. Most stressors for aquatic biota reflected an urban signature. Stressors related to ecological condition in our study were chemical and physical, such as developed land, impervious surface in the watershed, urban land in a buffer area around the stream (a 100-foot [30-meter]-wide area on each side of the stream, and maximum instantaneous discharge normalized by drainage area (a measure of flood and scour effects). Chemical stressors included low waterborne concentrations of dissolved oxygen and high concentrations of chloride, zinc and other metals, nutrients (nitrite and phosphorus), and fecal coliform bacteria. Although algae, invertebrates, and fish did not always demonstrate a significant response to the same stressors, higher abundances of high total phosphorus-indicator diatoms, lower ratings for invertebrate biotic integrity indexes and percentages of mayflies-stoneflies-caddisflies, and lower values for fish biotic integrity indexes underscored possible adverse effects of even low levels of developed land. Developed land is typically associated with more rapid runoff, which washes chemicals from impervious surfaces into area waterways and degrades stream habitat for aquatic communities. However, with respect to at least chloride from road salt, diatoms tolerant to dissolved salts were significantly lower with as little as 1-percent mixed forest in the watershed. Lower percentages of urban land in the stream buffer correlated with healthier aquatic assemblages of algae, invertebrates, and fish. The assessment of algal, invertebrate, and fish assemblages coupled with physical and chemical data were highly useful in evaluating the ecological status of aquatic communities at the 14 sites and for determining environmental stressors that may be contributing to reduced stream condition. Some of the stressors could potentially be removed or lessened with stream rehabilitation or changes in watershed management.

Wisconsin↗

USGS2020

Groundwater quality and geochemistry of West Virginia’s southern coal fields

Coal mining has been the dominant industry and land use in West Virginia’s southern coal fields since the mid-1800s. Mortality rates for a variety of serious chronic conditions, such as diabetes, heart disease, and some forms of cancer in Appalachian coal mining regions, are higher than in areas lacking substantial coal mining activity within the Appalachian Region or elsewhere in the United States. Causes of the increased mortality and morbidity are not clear, but poor diet, high rates of smoking, socioeconomic factors, and the quality of groundwater used by area residents are all possible contributing factors. This study was conducted by the U.S. Geological Survey in cooperation with the West Virginia Department of Health and Human Resources and the West Virginia Department of Environmental Protection, with grant support from the Centers for Disease Control and Prevention (CDC) to assess the quality of groundwater in southern West Virginia. The data from this assessment of groundwater quality may be used by the CDC and other agencies to potentially investigate the role or lack thereof of groundwater quality with respect to mortality and morbidity rates in the region. The study was conducted in a region where a high density of current or past coal mining combined with a lack of advanced sewage treatment could affect concentrations of commonly occurring constituents plus contaminants, including nitrate, trace metals, major ions, indicator bacteria, radon, hydrogen sulfide, and dissolved hydrocarbons. Because rural residential wells and mine outfalls are considered private sources of water in the region, and are therefore unregulated and unmonitored, water-quality data are sparse. To fill the data gap and assess the groundwater quality in the region, water-quality samples were collected from 60 sites in a 10-county area. The 60 sites sampled included 46 rural residential homeowner wells and 14 mine outfall discharges used for residential supply. For this study, all samples were collected prior to any filtration or other treatments, typically at the pressure tank, and are indicative of total and dissolved constituents in the untreated water. Generally, data for the 60 sites indicate that most waters sampled do not exceed thresholds for most U.S. Environmental Protection Agency (EPA) drinking-water standards and U.S. Geological Survey (USGS) drinking-water screening criteria. However, there were several notable exceptions. Turbidity exceeded the 5-Nephelometric Turbidity Unit (NTU) EPA treatment technique (TT) drinking-water standard in 14 of 60 (23 percent) sites sampled and exceeded the 1-NTU TT standard in 51 of 60 (85 percent) sites sampled. Turbidity is common in many wells in southern West Virginia and may be attributed to iron oxyhydroxide precipitates, sediment carried into the aquifers from the shallow soil zone due to improperly constructed or cased wells or transported to the aquifer in shallow stress-relief fracture zones or through permeable bedding-plane partings. For the sites sampled, 31 of 60 (52 percent) had pH values at, above, or below the upper and lower range of the EPA Secondary Maximum Contaminant Level (SMCL, 6.5–8.5 standard units). Of those 31 sites, 28 (90 percent) were indicative of acidic corrosive water and 3 (10 percent) were indicative of alkaline water. The Langelier Saturation Index (LSI), which is a measure of the corrosivity of the water, was computed for all sites sampled for the study. Eighty-two percent of the sites sampled had waters that were classified as corrosive, based on a LSI less than −0.5. Corrosive water has the potential to leach lead, copper, and other metals from lead, copper, galvanized, or lead-tin soldered connections in water lines. The chloride to sulfate mass ratio also was assessed with the alkalinity to indicate the potential to promote galvanic corrosion (PPGC) of water lines and plumbing fixtures. Only one of the sites (1.7 percent) classified as a corrosive water site, had a PPGC considered high; the remaining sites were classified as having either a moderate (53.3 percent) or low (45 percent) PPGC. Therefore, the type of plumbing systems sampled for this study may be affected by corrosive water, but the potential for leaching trace metals and other constituents from residential plumbing systems containing older galvanized pipes or lead-tin soldered copper pipes is moderate to low. The indicator bacteria total coliform and Escherichia coli ( E. coli ) also were detected in groundwater samples to varying degrees. Total coliforms, which are a broad class of indicator bacteria, are common in groundwater in southern West Virginia and were detected in 39 of the 60 sites (65 percent) sampled. The presence of total coliform bacteria is a potential indicator of surface contamination, due to improperly constructed or cased wells, or infiltration of soil or other surface contaminants into the aquifer or well bore. E. coli bacteria, however, are much more indicative of fecal contamination of groundwater from either human or animal sources, and 14 of the 60 (23 percent) sites sampled had detections of E. coli . Although only a few strains of E. coli are known pathogens, their presence in groundwater may be an indicator of other related pathogens such as viruses and should be regarded as a serious potential issue. Water treatment such as chlorination, ozonation, or ultraviolet light may be appropriate to kill potential pathogenic bacteria or viruses in the source water. Manganese and iron were prevalent contaminants in the groundwater samples collected for this study, with 30 of 60 (50 percent) sites analyzed for manganese and 25 of 60 (42 percent) sites analyzed for iron exceeding the proposed 50- and 300-micrograms per liter (µg/L) SMCL drinking-water standards, respectively, for aesthetic criteria such as taste, odor, or staining of plumbing fixtures. Fourteen of the 60 sites sampled (23 percent) had concentrations of manganese that exceeded the 300-µg/L USGS health-based screening level, and 1 site exceeded the 1,600-µg/L EPA drinking-water equivalent level, which is based on a lifetime exposure level. Sodium is another common constituent in groundwater within the study area. Sodium has an EPA health-based value (HBV) of 20 milligrams per liter (mg/L) for individuals who are on a sodium-restricted diet for blood pressure or other health reasons. Sodium concentrations exceeded the 20-mg/L EPA HBV in 27 of 60 (45 percent) samples. Radon, a naturally occurring carcinogenic radioactive gas known to cause lung cancer, was detected at concentrations at or exceeding the proposed 300-picocuries per liter (pCi/L) EPA Maximum Contaminant Level (MCL) in 12 of the 60 (20 percent) sites sampled. Sites with radon gas concentrations exceeding the 300-pCi/L proposed MCL have the potential for airborne concentrations of radon to exceed the 4-pCi/L indoor air standard. Inhalation of radon can cause lung cancer, and the 4-pCi/L indoor air standard is based on an inhalation standard. Therefore, homeowners whose wells have radon gas concentrations exceeding 300 pCi/L may be advised to have their indoor air tested to determine if indoor air concentrations exceed the 4-pCi/L indoor air standard established by the EPA. Various factors were analyzed statistically and graphically to determine whether they have an influence on groundwater quality within the study area, including topographic setting, well depth, type of mining (surface or underground), type of site (well or mine outfall), and geologic formation. Only geologic formation and the type of site sampled had strong statistical correlations with one or more of the constituents of concern for this study. The overall chemistry of outfalls (mine outfalls) and wells was significantly different, with a much higher dissolved oxygen content in outfalls than in wells. The dissolved oxygen content is the primary component driving the oxidation and reduction of minerals, and the precipitation of minerals that are saturated or super saturated with respect to various cations and anions. Median dissolved oxygen concentrations for the outfalls sampled was 8.75 mg/L, and only 0.4 mg/L for the wells sampled. Median concentrations of sulfate and selenium were much higher in waters from the outfalls sampled, with median concentrations of 73.75 mg/L and 2.35 µg/L, respectively, compared to the wells sampled, which had median concentrations of 18.3 mg/L and less than (<) the 0.05-µg/L method detection limit, respectively. The maximum selenium concentration was for a well, with a concentration of 16.6 µg/L. The geochemical processes that control sulfate and selenium concentrations in groundwater are similar and are the result of the oxidation of sulfide minerals such as pyrite and ferroselite. Iron and manganese concentrations were elevated in most of the wells sampled, with median concentrations of 269.5 and 124.5 µg/L, respectively, but were rarely detected in the outfalls sampled, with median concentrations of < 4.0 and < 0.4 µg/L, respectively. The difference in iron and manganese between wells and outfalls is indicative of the role of dissolved oxygen on processes controlling groundwater chemistry in the region. Three principal geologic formations were assessed for the study, and the overall chemistry for the Pocahontas, New River, and Kanawha Formations varied substantially with respect to several constituents. Concentrations of calcium, magnesium, and total dissolved solids were highest for sites sampled in the Pocahontas Formation, with median concentrations of 41.9, 18.6, and 312 mg/L, respectively. For constituents that are commonly associated with mining activity, the highest concentrations were for sites sampled in the New River Formation, with median concentrations of iron and manganese of 2,450 µg/L and 482 µg/L, respectively, and a median pH of 6.35 standard units. Concentrations of barium also were elevated in samples collected from sites in the New River Formation, with a median barium concentration of 184 µg/L. The source of the barium is not fully known but may be associated with commingling of shallow groundwater with deeper brines or dissolution of the mineral barite. The highest median sulfate concentrations were from sites sampled in the Pocahontas Formation, with a median concentration of 64.0 mg/L. Of the 12 sites at or exceeding the 300-pCi/L proposed drinking-water standard for radon, 8 (67 percent of MCL exceedances) were for sites deriving water from the Kanawha Formation, 3 (25 percent of MCL exceedances) were for sites deriving water from the New River Formation, and only 1 site was for water from the Pocahontas Formation (8 percent of proposed MCL exceedances). Dissolved hydrocarbons, including methane, ethane, propane, propene, n - and i -butane, 1-butene, n - and i -pentane, pentane, 2- and 3-ethyl pentane, hexane, and benzene were analyzed in samples collected from 59 of the 60 sites to assess the potential occurrence and sources of these trace gases in groundwater within the study area. Results of the analysis indicate that most of the gas is of shallow biogenic origin, possibly associated with coal-bed methane, but a subset of samples has a gas signature and a chloride to bromide ratio indicative of potential mixing with deeper thermogenic gases. Only 2 of the 59 (3.3 percent) sites sampled had concentrations of methane gas, which is a highly combustible and explosive gas, exceeding the 10 milligrams per kilogram level of concern established by the U.S. Office of Surface Mining Reclamation and Enforcement. Principal components analysis was used to assess the primary geochemical processes occurring in the aquifers sampled. The first principal component had significant positive loadings for bromide, chloride, silica, ammonia, barium, iron, manganese, and arsenic, and significant negative loadings for dissolved oxygen, potassium, nitrate, and uranium, and reflects reduction and oxidation (redox) processes occurring in deeper anoxic groundwater or shallow oxic groundwater. The strong positive loadings for iron, manganese, barium, and arsenic are correlated with reducing conditions often found deeper in the aquifer. More oxic water is correlated with oxidation of nitrogen species to nitrate and environmental mobilization of uranium and sulfate in shallow wells and mine outfalls.

West Virginia↗

USGS2024

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri River between Kansas City and St. Louis, Missouri, May 19–26, 2021

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, near nine bridges at eight highway crossings of the Missouri River between Kansas City and St. Louis, Missouri, from May 19 to 26, 2021. A multibeam echosounder mapping system was used to obtain channel-bed elevations for river reaches about 1,640 to 1,840 feet (ft) longitudinally and generally extending laterally across the active channel from bank to bank during low to moderate flood-flow conditions. These surveys provided channel geometry and hydraulic conditions at the time of the surveys and provided characteristics of scour holes that may be useful in developing or verifying predictive guidelines or equations for computing potential scour depth. These data also may be useful to the Missouri Department of Transportation as a low to moderate flood-flow assessment of the bridges for stability and integrity issues with respect to bridge scour during floods. Bathymetric data were collected around every in-channel pier. Scour holes were present at most piers for which bathymetry could be obtained, except those on banks or surrounded by riprap. Occasionally, scour holes were minor and difficult to discern from nearby dunes and ripples. All the bridge sites in this study were previously surveyed and documented in previous studies. Comparisons between bathymetric surfaces from the previous surveys and those of the current (2021) study do not indicate any consistent correlation between channel-bed elevations and streamflow conditions. The average difference between the bathymetric surfaces varied from 1.59 ft higher to 0.95 ft lower in 2021 than 2017, which corresponds to a gain of 100,200 cubic yards and a loss of 55,800 cubic yards, respectively. The average difference between the bathymetric surfaces varied from 2.74 ft higher to 3.05 ft lower in 2021 than 2013, which corresponds to a gain of 111,500 cubic yards and a loss of 169,200 cubic yards, respectively. The average difference between the bathymetric surfaces varied from 4.52 ft higher to 1.38 ft lower in 2021 than 2011, which corresponds to a gain of 221,100 cubic yards and a loss of 90,300 cubic yards, respectively. The most substantial overall net gain was 221,100 cubic yards between 2011 and 2021 at structures L0550 and A4497 at Jefferson City (site 20). The large net gain likely results from a combination of the mitigation of the scour holes near pier 4 of both bridges and the substantially lower flow in 2021 than in 2011. Alternatively, the most substantial overall net loss was 169,200 cubic yards between 2013 and 2021 at structure A6288 at Hermann (site 21), despite comparable streamflows. Pier size, nose shape, and skew to approach flow had a substantial effect on the size of the scour hole observed at a given pier. Larger and deeper scour holes were present at piers with wide or blunt noses caused by exposed footings or caissons. When a pier was skewed to primary approach flow, the scour hole was generally deeper and larger than at a similar pier without skew; furthermore, the shape of the scour hole near skewed piers in this study generally was longer and deeper on the side with impinging flow. At structure A6288 at Hermann (site 21), the scour hole near pier 5 was difficult to discern from nearby dunes and ripples, whereas the upstream edge of the footing was visible at pier 4, which likely contributes to the larger scour hole near that pier; the top of the footing may blunt the horseshoe vortex at pier 5, but the exposed front of the footing may exacerbate the vortex at pier 4.

Missouri↗

USGS2026

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of Virginia’s Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis conducted on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore shows that groundwater conditions vary substantially throughout the Coastal Plain aquifer system due to local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have fundamentally altered groundwater flow from pre-development conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers supplying groundwater is a relatively small portion of the total groundwater water budget for the aquifer system. Analyses of groundwater budgets from the Virginia Coastal Plain model show that groundwater flow is generally outward from the surficial aquifer to rivers and coastal water bodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer currently is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from areas outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2025, rates of withdrawal substantially lower than the historical maximum have resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or in adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. Downward groundwater flow continues under current conditions, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, those rates are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Estimated groundwater inflow from coastal regions could be contributing to saltwater intrusion, though that was not measured directly in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, show that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up over 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. Currently, estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The current withdrawal rate from the confined Yorktown-Eastover system is near the highest reported rate for the Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and small rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

EarthArXiv↗

USGS2026

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of the Virginia Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis based on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore indicates that groundwater conditions vary substantially throughout the Coastal Plain aquifer system because of local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have altered groundwater flow from predevelopment conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers is a relatively small part of the total groundwater budget for the aquifer system compared to the rate of recharge at the land surface. Analyses of groundwater budgets from the Virginia Coastal Plain model indicate that groundwater flow is generally outward from the surficial aquifer to rivers and coastal waterbodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2023, rates of withdrawal substantially lower than the historical maximum resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. For the most recent simulated conditions (2023) downward groundwater flow continues, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, rates of downward flow are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Groundwater inflow from coastal regions could be contributing to saltwater intrusion, even though that was not measured in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, indicate that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up more than 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. For the most recent simulated period (2023), estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The reported withdrawal rate for this period from the confined Yorktown-Eastover system is near the highest reported rate for the Virginia Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and low rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

Virginia↗

USGS2014

Simulation of the Lower Walker River Basin hydrologic system, west-central Nevada, using PRMS and MODFLOW models

Walker Lake is a terminal lake in west-central Nevada with almost all outflow occurring through evaporation. Diversions from Walker River since the early 1900s have contributed to a substantial reduction in flow entering Walker Lake. As a result, the lake is receding, and salt concentrations have increased to a level in which Oncorhynchus clarkii henshawi (Lahontan Cutthroat trout) are no longer present, and the lake ecosystem is threatened. Consequently, there is a concerted effort to restore the Walker Lake ecosystem and fishery to a level that is more sustainable. However, Walker Lake is interlinked with the lower Walker River and adjacent groundwater system which makes it difficult to understand the full effect of upstream water-management actions on the overall hydrologic system including the lake level, volume, and dissolved-solids concentrations of Walker Lake. To understand the effects of water-management actions on the lower Walker River Basin hydrologic system, a watershed model and groundwater flow model have been developed by the U.S. Geological Survey in cooperation with the Bureau of Reclamation and the National Fish and Wildlife Foundation. The watershed model was developed using the precipitation runoff modeling system (PRMS) and the groundwater flow model was constructed using the MODular groundwater FLOW model (MODFLOW) and both were calibrated for the lower Walker River Basin. These models can be incorporated in an integrated Groundwater and Surface-water FLOW (GSFLOW) model of the lower Walker River Basin. Additionally, the MODFLOW model developed for this study is useful for efficiently simulating long-term and large-scale effects of water-management actions on groundwater hydrology, streamflow, and Walker Lake level, volume, and dissolved-solids concentrations. The lower Walker River Basin PRMS model (LWR_PRMS) was constructed using a subbasin approach to aid in development and calibration, and simulates a 30-year period from 1978 to 2007 using daily time steps. The LWR_PRMS was used to estimate the distribution of groundwater recharge specified in the MODFLOW model. The highest rates of groundwater recharge occur in the Wassuk Range beneath perennial and ephemeral stream channels, whereas lower rates of recharge occur beneath alluvial fans along mountain fronts. The total groundwater recharge estimated using PRMS was about 25,000 acre-feet per year. The lower Walker River Basin MODFLOW (LWR_MF) model simulates an 89-year period using monthly time steps. The LWR_MF was constructed with an initial steady-state simulation to represent dynamic equilibrium conditions from 1908 to 1918 and then a transient simulation representing the period 1919–2007. The model was calibrated using a combination of manual and automated methods of adjusting model parameters to minimize errors between model simulated results and weighted observations of groundwater levels, streamflows, and lake level. Hydrologic conditions simulated with the LWR_MF include the movement and change in storage of groundwater, and the water budgets for Walker River, Walker Lake, and the groundwater system. The LWR_MF computed dissolved-solids concentrations for Walker Lake using simulated lake volume and an assumed constant internal salt mass of 37.2 million tons. Effects of potential changes in water management on future conditions (scenarios) of the lower Walker River Basin hydrologic system and Walker Lake from 2011 to 2070 were evaluated. Several water-management scenarios were considered, including a baseline scenario that represents no changes in system management, improved irrigation efficiencies for the Walker River Indian Irrigation Project (WRIIP), a range of increased streamflows entering the lower Walker River Basin, and, the fallowing of fields on the WRIIP. For the baseline scenario, it was assumed that streamflow conditions from 1981 to 2010 will be repeated in the future. Results indicate that Walker Lake level and volume continue to decline but at a slower rate as the surface area of the lake becomes smaller and lake evaporation decreases. Dissolved-solids concentrations in Walker Lake continue to increase and increase much more rapidly during periods when minimal flows reach the lake due to a diminished lake volume. Alternatively, in years with high runoff, lake level increases are greater and dissolved-solids decreases are greater, compared with equivalent runoffs experienced during 1981–2010. The simulated effects of improving WRIIP efficiencies on Walker River streamflows, Walker Lake inflow, level, and dissolved-solids concentrations, and crop consumptive use, are compared with the baseline reference scenario for a range of irrigation efficiency improvements from 0 to 25 percent over 60 years. Results indicate that water is conserved through a reduction in irrigation-induced groundwater recharge and subsequent groundwater discharge through evapotranspiration. The conserved water mostly goes to increased streamflow to Walker Lake, followed by increased crop consumptive use, then increased evaporation from Weber Reservoir. The simulated effects of increased streamflows at Walker River at Wabuska streamgage (10301500) on Walker Lake inflow, level, and dissolved-solids concentrations, and crop consumptive use, are compared with the baseline scenario after 60 years under two different management methods for Weber Reservoir. Results indicate Walker Lake level and dissolved-solids concentrations stabilized with increased irrigation-season streamflow of about 40,000 acre-feet per year at the Walker River at Wabuska streamgage. Walker Lake level increased, and dissolved-solids concentration decreased, with increased flows of 50,000 acre-feet per year or more. After 60 years with additional irrigation-season streamflows of 50,000 acre-feet per year, Walker Lake level increased by about 48 feet, and lake dissolved-solids concentrations decreased by about 3,000 milligrams per liter (mg/L). With 75,000 acre-feet per year of additional streamflow, Walker Lake level increased by 70 feet, and dissolved-solids concentration decreased by 7,600 milligrams per liter. The effects of fallowing of Walker River Indian Irrigation Project fields from 2007 to 2010 on Walker Lake inflow, level, and dissolved solids were evaluated. Fallowing resulted in a near doubling of Walker River inflow to Walker Lake during this period, an increase in Walker Lake level of about 1.4 feet, and a decrease in dissolved-solids concentration of about 540 mg/L.

Nevada↗

USGS2024

U.S. Geological Survey Earthquake Hazards Program decadal science strategy, 2024–33

Executive Summary Earthquakes represent one of our Nation’s most significant and costly natural hazards, with estimated annual loses from earthquakes close to $15 billion in 2023. Over the past two centuries, 37 U.S. States have experienced an earthquake exceeding a magnitude of 5, and 50 percent of States have a significant potential for future damaging shaking; these statistics speak to the need for nationwide interest and investment in earthquake hazard characterization and risk reduction. Authorized under the Earthquake Hazards Reduction Authorization Act, the U.S. Geological Survey (USGS) Earthquake Hazards Program (EHP) provides the scientific information, situational awareness, and knowledge necessary to reduce deaths, injuries, and economic losses from earthquakes and earthquake-induced tsunamis, landslides, and soil liquefaction. The EHP supports activities in three focused topical areas: (1) earthquake monitoring, (2) hazard assessment, and (3) applied research, using the results of each—and the coordination among them—to further support risk translation and communication in regions at risk nationwide. For earthquake monitoring, the Advanced National Seismic System (ANSS), a cooperative effort of USGS networks, university partner regional seismic networks, and real-time geodetic networks, collects and analyzes data on earthquakes; issues timely, reliable notifications of their occurrence and impacts; and provides data for earthquake research, hazard, and risk assessment as a foundation for building an earthquake-resilient Nation. The USGS-operated ShakeAlert Earthquake Early Warning system is a recent addition to EHP’s ANSS infrastructure. In the realm of earthquake hazard assessment, the EHP contributes to earthquake risk mitigation strategies by developing the National Seismic Hazard Model and maps, and other related products, that describe the likelihood and potential effects of earthquakes nationwide, especially in the urban areas of highest risk. The EHP also conducts research on the causes, characteristics, and effects of earthquakes and prioritizes work that directly increases the accuracy and precision of earthquake hazards assessments, earthquake forecasts, and earthquake monitoring and situational-awareness products and that supports the Nation’s earthquake mitigation practices. Bridging the EHP’s efforts across research, hazard assessments, and earthquake monitoring is a broad and comprehensive collection of earthquake information products, including the National Seismic Hazard Model, ShakeAlert, and other products describing impact, such as ShakeMap and PAGER (Prompt Assessment of Global Earthquakes for Response), which have been developed and integrated into EHP’s real-time monitoring systems. EHP funds external partners to carry out many important collaborative activities through an active external grants program—one of the largest in the USGS—and through cooperative agreements with other partners such as the university-operated regional seismic networks, funded as part of the ANSS. To continue its support of earthquake hazard characterization and risk reduction, the EHP aims to strengthen its foundational products and practices while positioning itself to respond to the evolving needs of the Nation and follow best practices of the scientific community. This document describes a strategy for the program to ensure it can meet these demands. The foundational priorities outlined in this strategy represent those activities that remain critical to the core functionality of the program and those that can be supported under current fiscal year 2024-level appropriations. Priorities described as aspirational are important for future growth, and to maintain the program’s position as a leading global resource in earthquake science, but would require increases in appropriated funding to be fully realized. Across the program’s portfolio of activities, several major themes have been identified as the most critical activities to advance EHP science over the coming decade. Together, these activities provide the framework necessary to integrate critical hazard characterization and risk reduction activities across the program. They provide the structure for research to advance the understanding of where, when, and why earthquakes occur and how we can use improved knowledge to drive short-term and actionable forecasts of seismic activity. They expand the usefulness of critical earthquake products and advance the sophistication of those products to keep pace with the rapidly evolving needs of an ever-expanding user base while maintaining the position of the USGS as a global leader in earthquake science. Focus on system-level science. Establish an automated earthquake-processing pipeline. Enhance the accuracy and reliability of the ShakeAlert Earthquake Early Warning system and plan for extension to other regions. Implement time-dependent earthquake forecasting. Develop physically realistic models. Expand computational capacity. This science strategy is organized into three primary sections. The first section provides an overview of the EHP and its budget, governance, and program council. Readers familiar with the program may wish to focus on the second section, which describes the core of the science strategy, including priorities across each of the EHP’s major program activities in monitoring, hazard assessment, and targeted research. The third section outlines science priorities that cut across program activities, including those involving collaborations external to the EHP.

American Samoa, Guam, Northern Marianas Islands, U↗

USGS1997

Influence of evaporation, ground water, and uncertainty in the hydrologic budget of Lake Lucerne, a seepage lake in Polk County, Florida

Evaporation losses and the interaction of ground water with Lake Lucerne were studied to determine the influence of these two processes on the hydrologic budget of a seepage lake. Lake Lucerne is representative of the numerous seepage lakes of sinkhole origin in the karst terrain of central Florida. Because of permeable surficial deposits, ground-water inflow is the only significant contribution from the surrounding watershed. The lake recharges the underlying Upper Floridan aquifer and, as a result, is susceptible to increased leakage induced by pumping from this aquifer. Ground-water fluxes determined in the study were analyzed to define the proportion of the total lake leakage induced by pumping from the Upper Floridan aquifer. A hydrologic budget is analyzed for the 1-year period from October 1985 to September 1986. Ground-water inflow and leakage are significant components of the hydrologic budget. Changes in the quantity of either of these fluxes can substantially alter lake stage. Ground-water inflow contributed from 20 to 37 percent of the total annual inflow to the lake. Leakage from the lake accounted for 18 to 23 percent of the total annual outflow. Water withdrawals from the Upper Floridan aquifer increased annual lake leakage by 22 percent over nonpumping conditions. Most of the increase (92 percent) in leakage occurred during April, May, and June 1986, when local citrus irrigation was highest. For the study year, ground-water inflow and leakage volumes were calculated by flow-net analysis to be equal to 10.5 and 12.6 inches, respectively, of water depth above (or below) the lake surface. These estimates were revised upward on the basis of an analysis of the error in the hydrologic-budget equation. Revised ground water inflow exceeded annual leakage from the lake. Ground-water inflow rates were increased by 120 percent to 23.6 inches, and leakage was increased by 40 percent to 17.5 inches. Differences between the two estimates probably reflect the uncertainty in the hydraulic conductivity estimates of the porous media around the lake and the unaccounted effect of transient ground-water inflow. The geometry of the sinkhole complex beneath Lake Lucerne and pumping in the Upper Floridan aquifer are primary controls on ground water interactions with the lake and, in particular, lake leakage. A numerical ground-water model was used to test the effects of these two factors on ground-water interactions with the lake. Results indicate that the intermediate confining unit below Lake Lucerne has been breached and replaced by materials about two orders of magnitude more conductive. Anisotropy in the surficial aquifer is approximately 100 and controls the depth of the ground-water flow intercepted by the lake. Lake sediments having low permeability may control the distribution of leakage through the lakebed but did not appreciably reduce total leakage rates in these simulations. Evaporation loss was the major outflow component of the hydrologic budget. Annual lake evaporation determined by the energy-budget method was 57.9 inches, about 8 inches greater than long-term estimates for the region. The greater rate was attributed to drier than normal conditions: rainfall totaled 40.9 inches during the study year, about 10 inches less than the long-term average. Similar annual evaporation rates were determined by the energy-budget method, the simpler mass-transfer method, and by corrected pan evaporation from an onsite pan; however, the probable errors associated with these other two methods were greater than for the energy budget. Weekly energy-budget evaporation rates ranged from 0.04 inch per day in early January 1986 to 0.26 inch per day in early May 1986. The largest monthly energy-budget evaporation rates occurred in April and May 1986, 7.16 and 7.12 inches per month, respectively. Monthly evaporation estimated from corrected pan evaporation generally was within 10 percent of the energy-budget estimate but differed by as much as 35 percent. Daily energy-budget evaporation also was computed for 321 days during the year.

Florida↗

USGS1947

Age and growth of the kiyi, Leucichthys kiyi Koelz, in Lake Michigan

Ages were determined and individual growth histories were calculated from the examination and measurement of the scales of 1,649 kiyis captured at seven localities in Lake Michigan in 1931 and 1932. The numbers of individuals employed for the investigation of other phases of the life history (such as length-frequency distributions, length-weight relationship, and sex ratio) varied according to the amount of data available or required. Age-group IV was dominant in the 1931 collections from Racine, Port Washington, and Kewaunee, Wisconsin, and age-group V dominated the 1932 samples from the Fox Islands and from three localities southward of Manistique, Michigan. A trend was noticeable toward an increase in average age from south to north. Among the explanations suggested for the observed differences in age composition were: Variation with latitude in the natural span of life; differences in fishing intensity; fluctuations in the strength of year classes (to account possibly for the shift in the dominant age group from 1931 to 1932). The oldest male kiyi belonged to the VII group and the oldest female was a member of the X group. The possible distorting effects of such factors as gear selection traceable to differences in the mesh sizes of nets fished in 1931 and 1932, selection by nets on the basis of the condition (K) of the fish, and local variations in fishing intensity and hence in the selective destruction of rapidly growing individuals in the fishery were held to be sufficiently great to render doubtful the significance of most of the observed local differences in growth rate. Kiyis from all samples were combined to determine the general growth in length. The growth in weight of the Fox Islands fish, however, was considered separately as these fish were consistently lighter than kiyis of corresponding length from other localities. The Lake Michigan kiyi grows slowly, with the females growing slightly more rapidly than the males. The grand average calculated lengths indicated, for example, that the females did not attain a total length of 10 inches until the fifth year of life or the males until the sixth. Similarly, the calculated weight of 4 ounces was not reached until the fifth or sixth year (with the actual time varying with sex and locality). The season's growth of the kiyi probably begins sometime in May and most or all of the growth is completed by the end of August. The calculated lengths of the age groups exhibited large discrepancies that differed from “Lee's phenomenon” as ordinarily observed in that the data for the later rather than the earlier years of life were affected most severely. Chief among the factors held responsible for these discrepancies were gear selection and the selective destruction of the more rapidly growing individuals in the fishery. Errors inherent in the (direct-proportion) method of computing growth from scale measurements were considered to have been unimportant. The Lake Michigan kiyi exhibits growth compensation–the tendency for the smaller of the young fish to have the more rapid growth in the later years of life. Comparisons with the average lengths and weights of the age groups of the Lake Ontario kiyi given by Pritchard (1931) indicated the Lake Michigan fish to be the larger at the earlier ages (age-groups II and III) and the smaller at the later ages (age-groups IV to VI). The length-frequency distributions of the age groups exhibited extensive overlap. As many as eight age groups were represented in a single centimeter interval of length. The length frequencies and average lengths of all fish collected, arranged according to the mesh sizes of the gill nets by which they were captured, revealed that the selective action of these nets in the taking of kiyis was much more obvious in the numbers of fish in the catch than in their average size. As an illustration, in 1930–1931, the 2 3/4-inch mesh nets took fish that were only 0.1 inch longer than those in 2 1/2-inch meshes but captured less than one fourth as many. Gill nets fished in northern Lake Michigan in 1932 captured kiyis that averaged 0.2 to 0.4 inch longer than those taken in the same meshes in southern Lake Michigan in 1930–1931. Because of the more slender form of kiyis from the northeastern island region of Lake Michigan, data on the general length weight relationship were compiled separately for fish of that area and for those of the great central basins of the lake. In both regions the weight increased to a power slightly greater than the cube of the length. Available information on condition indicated that the coefficient (K) was higher in August and early September than in May, June, and July. Condition declined from early September to October and early November–the latter period the time of most active spawning. Spawning itself was accompanied by an additional loss of about 12 per cent of the body weight of females and of somewhat less than 2 per cent of the weight of males. Analysis of the variations of K within a group that was homogeneous with respect to age, sex, maturity, and time of collection revealed that a net of a particular mesh size tends to take the heavier of the shorter fish and the lighter of the longer fish within its range of effectiveness. Among fish of the same length the values of K tended to increase with increase in the mesh size of the nets employed for their capture. Practically all fish in the samples were mature (only 11 immature in more than 6,000). These “immature” fish were probably “non-functional” since all of them approached or exceeded the average length of the mature kiyis. Females were strongly predominant in the collections at all seasons but were relatively more plentiful during the summer (90 per cent of the total) than during the spawning period (75 per cent). Possible factors contributing to this predominance of females and to the change in the sex ratio at the spawning season were discussed. A decrease in the relative abundance of males with increase in age appears to be characteristic of the kiyi. This decrease indicates a differential mortality of th sexes (greater relative destruction of males in the spawning period when they are unusually abundant or a greater natural mortality rate for the males). Current fishery regulations on mesh size and closed seasons afford the kiyi good protection but offer no guarantee against depletion from too intensive fishing.

Transactions of the American Fisheries Society↗

USGS2004

Geohydrology of the French Creek Basin and simulated effects of drought and ground-water withdrawals, Chester County, Pennsylvania

This report describes the results of a study by the U.S. Geological Survey, in cooperation with the Delaware River Basin Commission, to develop a regional ground-water-flow model of the French Creek Basin in Chester County, Pa. The model was used to assist water-resource managers by illustrating the interconnection between ground-water and surface-water systems. The 70.7-square mile French Creek Basin is in the Piedmont Physiographic Province and is underlain by crystalline and sedimentary fractured-rock aquifers. Annual water budgets were calculated for 1969-2001 for the French Creek Basin upstream of streamflow-measurement station French Creek near Phoenixville (01472157). Average annual precipitation was 46.28 in. (inches), average annual streamflow was 20.29 in., average annual base flow determined by hydrograph separation was 12.42 in., and estimated average annual ET (evapotranspiration) was 26.10 in. Estimated average annual recharge was 14.32 in. and is equal to 31 percent of the average annual precipitation. Base flow made up an average of 61 percent of streamflow. Ground-water flow in the French Creek Basin was simulated using the finite-difference MODFLOW-96 computer program. The model structure is based on a simplified two-dimensional conceptualization of the ground-water-flow system. The modeled area was extended outside the French Creek Basin to natural hydrologic boundaries; the modeled area includes 40 square miles of adjacent areas outside the basin. The hydraulic conductivity for each geologic unit was calculated from reported specific-capacity data determined from aquifer tests and was adjusted during model calibration. The model was calibrated for above-average conditions by simulating base-flow and water-level measurements made on May 1, 2001, using a recharge rate of 20 in/yr (inches per year). The model was calibrated for below-average conditions by simulating base-flow and water-level measurements made on September 11 and 17, 2001, using a recharge rate of 6.2 in/yr. Average conditions were simulated by adjusting the recharge rate until simulated streamflow at streamflow-measurement station 01472157 matched the long-term (1968-2001) average base flow of 54.1 cubic feet per second. The recharge rate used for average conditions was 15.7 in/yr. The effect of drought in the French Creek Basin was simulated using a drought year recharge rate of 8 in/yr for 3 months. After 3 months of drought, the simulated streamflow of French Creek at streamflow-measurement station 01472157 decreased 34 percent. The simulations show that after 6 months of average recharge (15.7 in/yr) following drought, streamflow and water levels recovered almost to pre-drought conditions. The effect of increased ground-water withdrawals on stream base flow in the South Branch French Creek Subbasin was simulated under average and drought conditions with pumping rates equal to 50, 75, and 100 percent of the Delaware River Basin Commission Ground Water Protected Area (GWPA) withdrawal limit (1,393 million gallons per year) with all pumped water removed from the basin. For average recharge conditions, the simulated streamflow of South Branch French Creek at the mouth decreased 18, 28, and 37 percent at a withdrawal rate equal to 50, 75, and 100 percent of the GWPA limit, respectively. After 3 months of drought recharge conditions, the simulated streamflow of South Branch French Creek at the mouth decreased 27, 40, and 52 percent at a withdrawal rate equal to 50, 75, and 100 percent of the GWPA limit, respectively. The effect of well location on base flow, water levels, and the sources of water to the well was simulated by locating a hypothetical well pumping 200 gallons per minute in different places in the Beaver Run Subbasin with all pumped water removed from the basin. The smallest reduction in the base flow of Beaver Run was from a well on the drainage divide between the French Creek Basin and the Marsh Creek Basin to the south; the simulated base flow of Beaver Run at the mouth was reduced 1 percent. The greatest reduction in the base flow of Beaver Creek was from a well close to Beaver Run; the simulated base flow of Beaver Run at the mouth was reduced 8 percent. The simulations showed that (1) if the contributing area of a well is in a basin, pumping will affect stream base flow and water levels in that basin whether the well is inside or outside that basin; (2) wells in different areas of a basin away from a divide produce a similar reduction in base flow; (3) a well within a basin will derive more water from diverted base flow and less water from storage than a well on or near a basin divide; and (4) the reduction in base flow at the mouth of the stream is the same for a well in the headwaters and a well downstream near the confluence. Model simulations illustrate some of the typical analyses and results that can be produced. The model was calibrated using annual values for recharge and ground-water ET and then was run using the annual values in a seasonally independent transient mode to show changes with time. The timing and relative magnitude of some of the changes simulated with the model when viewed in terms of a normal climatic year may be subject to considerable uncertainty because of the variability in seasonal recharge and ground-water ET rates. Transient model simulations for short-term periods are indicative of possible hydrologic system response and are considered an approximation.

Pennsyvlania↗

USGS2024

Status and understanding of groundwater quality in the Mojave Basin Domestic-Supply Aquifer study unit, 2018—California GAMA Priority Basin Project

Groundwater quality in the western part of the Mojave Desert in San Bernardino County, California, was investigated in 2018 as part of the California State Water Resources Control Board Groundwater Ambient Monitoring and Assessment Program Priority Basin Project. The Mojave Basin Domestic-Supply Aquifer study unit (MOBS) region was divided into two study areas—floodplain and regional—to assess differences between the two major aquifers used for drinking water supply in the area. This assessment characterized the quality of ambient groundwater and not the quality of treated drinking water. The study included three components: (1) a status assessment, which characterized the quality of groundwater resources used for domestic drinking-water supply in the floodplain and regional study areas; (2) a brief understanding assessment, which evaluated factors that could potentially affect the quality of groundwater used by domestic wells in the region; and (3) a comparative assessment between the groundwater resources used by domestic wells and public-supply wells in the two study areas. The domestic-well assessment was based on data collected by the U.S. Geological Survey from 48 domestic wells in January–May 2018. The public-supply assessment was based on data for samples from 322 public-supply wells in 2008–18, either collected by the U.S. Geological Survey or compiled from the California State Water Resources Control Boards Division of Drinking Water publicly available database. Concentrations of water-quality constituents in ambient groundwater were compared to regulatory and non-regulatory benchmarks typically used by the State of California and Federal agencies as health-based or aesthetic standards for public drinking water. Relative concentrations, defined as the measured concentration divided by the benchmark concentration, were classified as high (greater than 1.0), moderate (greater than 0.5 for inorganic constituents or 0.1 for organic and special-interest constituents, and not high), or low (concentrations lower than moderate). The floodplain and regional study areas were divided into 15 and 35 grid cells, respectively, and grid-based methods were used to compute the areal proportions of the two study areas with high, moderate, or low relative concentrations of individual constituents and classes of constituents. For the domestic-supply assessment, one or more inorganic constituents with health-based benchmarks were detected at high relative concentrations in 58 percent of the regional study area and 13 percent of the floodplain study area. The inorganic constituents with health-based benchmarks detected at high relative concentrations in the regional study area were arsenic, chromium and hexavalent chromium, fluoride, adjusted gross alpha particle activity, uranium, molybdenum, strontium, and nitrate; only arsenic was detected at high relative concentrations in the floodplain study area. One or more inorganic constituents with secondary maximum contaminant level benchmarks were detected at high concentrations in 15 and 6.7 percent of the regional and floodplain study areas, respectively. The constituents detected at high relative concentrations in the regional study area were total dissolved solids, chloride, sulfate, and iron; only total dissolved solids and sulfate were detected at high relative concentrations in the floodplain study area. Organic constituents were not detected at moderate or high relative concentrations in either the regional or floodplain study areas. Volatile organic compounds were detected at low relative concentrations in 21 and 27 percent of the regional and floodplain study areas, respectively, and pesticides were detected at low relative concentrations in 9.1 and 20 percent of the regional and floodplain study areas, respectively. The only individual organic constituent detected in more than 10 percent of either study area was the trihalomethane trichloromethane. Total coliform bacteria were detected in 15 and 27 percent of the grid wells in the regional and floodplain study areas, respectively. The greater prevalence of high relative concentrations of many inorganic constituents in the regional study area compared to the floodplain area likely indicates the greater diversity of geologic material at depth in aquifer material and generally finer-grained alluvium compared to the floodplain study area combined with generally older groundwater that has had more contact time with aquifer materials. In general, trace element concentrations (1) increased with increasing groundwater age, (2) increased with distance from recharge sources in the mountains, and (3) increased with closer proximity to some types of geological units. In general, groundwater from domestic wells in the floodplain study area is young, with most samples containing a component of modern groundwater based on tritium and unadjusted carbon-14 activities, whereas groundwater from domestic wells in the regional study area generally is old, with most samples having unadjusted carbon-14 ages of 5,000–40,000 years. Public-supply wells in MOBS generally were deeper than domestic wells and presumably are in contact with older, more weathered alluvium that may have more mobile trace elements, such as arsenic or uranium. However, only 26 percent of the public-supply regional study area had high relative concentrations of inorganic constituents, compared to 58 percent for the domestic regional study area. The percentages of the public-supply and domestic floodplain study areas with high relative concentrations of inorganic constituents were 11 and 13 percent, respectively. The ages of groundwater used by public-supply and domestic wells in each study area were similar, which was not expected given the greater depth of the public-supply wells. Three potential factors may contribute to these results: (1) greater spatial footprint of domestic well network, which may result in domestic wells pumping groundwater from fractured bedrock or mineralized areas not used by public-supply wells; (2) greater pumping rates in public-supply wells, resulting in more water being withdrawn from coarse-grained, heterogeneous alluvium than finer-grained layers, which may have higher concentrations of (or more mobile) inorganic constituents; and (3) a greater degree of well management with public-supply wells, which may include pausing use of or decommissioning wells if treating or blending water is not feasible to lower constituent concentrations.

California↗

USGS2007

Mercury in precipitation in Indiana, January 2001–December 2003

Mercury in precipitation was monitored during 2001 through 2003 at four locations in Indiana as part of the National Atmospheric Deposition Program-Mercury Deposition Network (NADP-MDN). Monitoring stations were operated at Roush Lake near Huntington, Clifty Falls State Park near Madison, Monroe County Regional Airport near Bloomington, and Indiana Dunes National Lakeshore near Porter. At these monitoring stations, precipitation amounts were measured and weekly samples were collected for analysis of total mercury and methylmercury by low-level methods. Wet deposition was computed with the total mercury and methylmercury concentrations and the precipitation amounts. In 3 years of weekly samples collected at the four monitoring stations, the volume-weighted total mercury concentration was 11.5 ng/L (nanograms per liter). As a reference for comparison, the total mercury concentration in 47 percent of the samples analyzed was greater than the Indiana water-quality standard for mercury (12 ng/L, protecting aquatic life) and nearly all of the concentrations exceeded the Indiana water-quality standards for mercury in the Great Lakes system (1.8 ng/L, protecting human health, and 1.3 ng/L, protecting wild mammals and birds). The precipitation-weighted concentrations at three of the monitoring stations in Indiana in 2003 were in the top 40 percent of all monitoring stations in the NADP-MDN and the concentration at Indiana Dunes was the eighth highest in the NADP-MDN for 2003. At the four monitoring stations during the study period, the mean weekly total mercury deposition was 243 ng/m 2 (nanograms per square meter) and mean annual total mercury deposition was 12,623 ng/m 2 . The annual mercury deposition at the four monitoring stations in Indiana in 2003 was in the top 40 percent of all monitoring stations in the NADP-MDN and the annual mercury deposition at the Clifty Falls station was the tenth highest in the NADP-MDN for 2003. For the 3-year period, the median methylmercury concentration in weekly samples was 0.058 ng/L with a maximum of 5.77 ng/L. Normalized methylmercury deposition was 2.09 ng/m 2 per inch of precipitation and methylmercury deposition was 0.7 percent of the total mercury deposition. The annual and mean weekly methylmercury deposition was highest at the Roush Lake station. Among the monitoring stations in the NADP-MDN with methylmercury data, methylmercury deposition at the monitoring stations in Indiana appeared to be higher than at eight stations in Wisconsin and Minnesota for that same time period, although methylmercury concentrations in Indiana were similar to or lower than those in Wisconsin and Minnesota. Geographically, the weekly total mercury concentrations at Indiana Dunes and Clifty Falls were statistically higher than concentrations at Bloomington, although a statistical difference in weekly total mercury deposition was not found among the four monitoring stations. Annual mercury emissions from sources in the vicinity of Indiana Dunes and Clifty Falls in 2001 were more than 10 times those at Bloomington, although other factors may help explain the differences in total mercury concentrations, such as the types of mercury emissions, mercury transport from sources outside Indiana, and meteorological conditions. Mercury concentrations and deposition varied at the four monitoring stations during the 3-year period. Total mercury concentrations in weekly samples ranged from 1.54 to 77 ng/L and weekly mercury deposition ranged from 0.8 to 2,456 ng/m 2 . Data from weekly samples exhibited seasonal patterns. Total mercury concentrations and deposition were highest in spring and summer and lowest in winter. Methylmercury concentrations were highest in winter and methylmercury deposition was highest in spring. Annual precipitation at the four monitoring stations was highest in 2003, exceeding the precipitation normals in spring and summer 2003. Annual mercury deposition in 2003 at the Roush Lake, Clifty Falls, and Indiana Dunes was as much as 41 to 67 percent higher in 2003 than in 2001 or in 2002 at those stations. Total mercury deposition that was more than 10 percent of the mean annual deposition (1,262 ng/m 2 ) was recorded in 11 of 551 weekly samples from the study period. These samples contained approximately 3 inches or more of rain and most were collected in spring and summer 2003. The highest deposition (2,456 ng/m 2 in a sample from Roush Lake) was 15.7 percent of the annual deposition at that station and approximately 10 times the mean weekly deposition for Indiana. High deposition recorded in three weekly samples at Clifty Falls contributed 31 percent of the annual deposition at that station in 2003. Weekly samples with high mercury deposition may help to explain the differences in annual mercury deposition among the four monitoring stations in Indiana.

Indiana↗