USGS Science⌕ Search

SEARCH · USGS Science

Results for “Case Studies in the Environment”

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

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

447 records · Page 25Linked to original sources

Occam's shadow: levels of analysis in evolutionary ecology - where to next?

Evolutionary ecology is the study of evolutionary processes, and the ecological conditions that influence them. A fundamental paradigm underlying the study of evolution is natural selection. Although there are a variety of operational definitions for natural selection in the literature, perhaps the most general one is that which characterizes selection as the process whereby heritable variation in fitness associated with variation in one or more phenotypic traits leads to intergenerational change in the frequency distribution of those traits. The past 20 years have witnessed a marked increase in the precision and reliability of our ability to estimate one or more components of fitness and characterize natural selection in wild populations, owing particularly to significant advances in methods for analysis of data from marked individuals. In this paper, we focus on several issues that we believe are important considerations for the application and development of these methods in the context of addressing questions in evolutionary ecology. First, our traditional approach to estimation often rests upon analysis of aggregates of individuals, which in the wild may reflect increasingly non-random (selected) samples with respect to the trait(s) of interest. In some cases, analysis at the aggregate level, rather than the individual level, may obscure important patterns. While there are a growing number of analytical tools available to estimate parameters at the individual level, and which can cope (to varying degrees) with progressive selection of the sample, the advent of new methods does not reduce the need to consider carefully the appropriate level of analysis in the first place. Estimation should be motivated a priori by strong theoretical analysis. Doing so provides clear guidance, in terms of both (i) assisting in the identification of realistic and meaningful models to include in the candidate model set, and (ii) providing the appropriate context under which the results are interpreted. Second, while it is true that selection (as defined) operates at the level of the individual, the selection gradient is often (if not generally) conditional on the abundance of the population. As such, it may be important to consider estimating transition rates conditional on both the parameter values of the other individuals in the population (or at least their distribution), and population abundance. This will undoubtedly pose a considerable challenge, for both single- and multi-strata applications. It will also require renewed consideration of the estimation of abundance, especially for open populations. Thirdly, selection typically operates on dynamic, individually varying traits. Such estimation may require characterizing fitness in terms of individual plasticity in one or more state variables, constituting analysis of the norms of reaction of individuals to variable environments. This can be quite complex, especially for traits that are under facultative control. Recent work has indicated that the pattern of selection on such traits is conditional on the relative rates of movement among and frequency of spatially heterogeneous habitats, suggesting analyses of evolution of life histories in open populations can be misleading in some cases.

Journal of Applied Statistics↗

Dynamics of national forests assessed using the Landsat record: Case studies in eastern United States

The national forests (NFs) in the United States are protected areas managed for multiple purposes, and therefore are subject to both natural and anthropogenic disturbances. Monitoring forest changes arising from such disturbances and the post-disturbance recovery processes is essential for assessing the conditions of the NFs and the effectiveness of management approaches. In this study, we used time series stacks of Landsat images (LTSS) to evaluate the dynamics of seven NFs in eastern United States, including the De Soto NF, the Talladega NF, the Francis Marion NF, and the Uwharrie NF in southeastern U.S., and the Chequamegon NF, the Hiawatha NF, and the Superior NF in northern U.S. Each LTSS consisted of 12–14 Landsat images acquired for the same location, spanning from 1984 to 2006 with a nominal interval of one image every 2 years. Each LTSS was analyzed using a vegetation change tracker (VCT) algorithm to map forest disturbance. Accuracy assessments of the derived disturbance maps revealed that they had overall accuracy values of about 80%, with most of the disturbance classes having user's accuracies ranging from 70% to 95%. The producer's accuracies were generally lower, with the majority being in the range between 50% and 70%. While this may suggest that the disturbance maps could slightly underestimate disturbances, a more detailed assessment of the omission errors revealed that the majority of the disagreements were due to minor disturbances like thinning or storm damages that were identified by the image analysts but were not captured by the VCT algorithm. The derived disturbance year maps revealed that while each of the seven NFs consisted of 90% or more forest land, significant portions of the forests were disturbed since 1984. Mapped disturbances accounted for about 30%–45% of total land area in the four NFs in southeastern U.S. and about 10%–20% in the three NFs in northern U.S. The disturbance rates were generally higher in the buffer zones surrounding each NF, and varied considerably over time. The time series approach employed in this study represents a new approach for monitoring forest resources using the Landsat or similar satellite data records. The disturbance products derived using this approach were spatially explicit and contained much more temporal details than conventional bi-temporal change products, and likely will be found more useful by many users including ecologists and resources managers. The high disturbance rates found in the southeastern U.S. suggest that this region may have a more significant role in modulating the atmospheric carbon budget than currently recognized.

Remote Sensing of Environment↗

Natural offshore oil seepage and related tarball accumulation on the California coastline — Santa Barbara Channel and the southern Santa Maria Basin; source identification and inventory

Oil spillage from natural sources is very common in the waters of southern California. Active oil extraction and shipping is occurring concurrently within the region and it is of great interest to resource managers to be able to distinguish between natural seepage and anthropogenic oil spillage. The major goal of this study was to establish the geologic setting, sources, and ultimate dispersal of natural oil seeps in the offshore southern Santa Maria Basin and Santa Barbara Basins. Our surveys focused on likely areas of hydrocarbon seepage that are known to occur between Point Arguello and Ventura, California. Our approach was to 1) document the locations and geochemically fingerprint natural seep oils or tar; 2) geochemically fingerprint coastal tar residues and potential tar sources in this region, both onshore and offshore; 3) establish chemical correlations between offshore active seeps and coastal residues thus linking seep sources to oil residues; 4) measure the rate of natural seepage of individual seeps and attempt to assess regional natural oil and gas seepage rates; and 5) interpret the petroleum system history for the natural seeps. To document the location of sub-sea oil seeps, we first looked into previous studies within and near our survey area. We measured the concentration of methane gas in the water column in areas of reported seepage and found numerous gas plumes and measured high concentrations of methane in the water column. The result of this work showed that the seeps were widely distributed between Point Conception east to the vicinity of Coal Oil Point, and that they by in large occur within the 3-mile limit of California State waters. Subsequent cruises used sidescan and high resolution seismic to map the seafloor, from just south of Point Arguello, east to near Gaviota, California. The results of the methane survey guided the exploration of the area west of Point Conception east to Gaviota using a combination of seismic instruments. The seafloor was mapped by sidescan sonar, and numerous lines of high -resolution seismic surveys were conducted over areas of interest. Biomarker and stable carbon isotope ratios were used to infer the age, lithology, organic matter input, and depositional environment of the source rocks for 388 samples of produced crude oil, seep oil, and tarballs mainly from coastal California. These samples were used to construct a chemometric fingerprint (multivariate statistics) decision tree to classify 288 additional samples, including tarballs of unknown origin collected from Monterey and San Mateo County beaches after a storm in early 2007. A subset of 9 of 23 active offshore platform oils and one inactive platform oil representing a few oil reservoirs from the western Santa Barbara Channel were used in this analysis, and thus this model is not comprehensive and the findings are not conclusive. The platform oils included in this study are from west to east: Irene, Hildago, Harvest, Hermosa, Heritage, Harmony, Hondo, Holly, Platform A, and Hilda (now removed). The results identify three “tribes” of 13 C-rich oil samples inferred to originate from thermally mature equivalents of the clayey-siliceous, carbonaceous marl, and lower calcareous-siliceous members of the Monterey Formation. Tribe 1 contains four oil families having geochemical traits of clay-rich marine shale source rock deposited under suboxic conditions with substantial higher-plant input. Tribe 2 contains four oil families with intermediate traits, except for abundant 28,30-bisnorhopane, indicating suboxic to anoxic marine marl source rock with hemipelagic input. Tribe 3 contains five oil families with traits of distal marine carbonate source rock deposited under anoxic conditions with pelagic but little or no higher-plant input. Tribes 1 and 2 occur mainly south of Point Conception in paleogeographic settings where deep burial of the Monterey Formation source rock favored generation from all three members or their equivalents. In this area, oil from the clayey-siliceous and carbonaceous marl members (Tribes 1 and 2) may overwhelm that from the lower calcareous-siliceous member (Tribe 3) because the latter is thinner and less oil-prone than the overlying members. Tribe 3 occurs mainly north of Point Conception, where shallow burial caused preferential generation from the underlying lower calcareous-siliceous member or another unit with similar characteristics. It is very desirable to be able to clearly distinguish the naturally occurring seep oils from the anthropogenically derived platform oils. Within the “training set” of oils and tars (388 samples), the biomarker parameters are sometimes sufficient to allow unique discrimination of individual platform oils. More often however, platform samples and seep samples with sources geographically close to each other are too similar to each other, with respect to the biomarker parameters, to definitively differentiate them on that basis alone. In some cases other parameters can be helpful. These other parameters are related to the degree of biogeochemical degradation or weathering that the oils or tars have experienced. These components include the typical oil distribution of n-alkane hydrocarbons and isoprenoids pristane and phytane. All of the platform oils in our sample set contain these components. On the other hand, the seep oils or tars have been exposed to significant biodegradation while in the near subsurface. The majority, but not all of seep oils or tars have been biodegraded up to or beyond the loss of n-alkanes and isoprenoids. Seep oils found in the vicinity of Coal Oil Point or Arroyo Burro are apparently the least weathered and are particularly likely to retain significant n-alkanes and isoprenoids. Therefore the combination of chemometric fingerprinting and the presence or absence of n-alkanes and isoprenoids help to differentiate anthropogenic production oils versus natural seeps oils and tars. The differentiation is not always definitive because of the close chemical similarity of some samples and the variability in the biodegradation progression. This is the case near Coal Oil Point, and near Platform A (Dos Cuadros Field) where seep oils and Platform Holly and Platform A oils are genetically very similar and cannot be definitively distinguished after a period of a few days of weathering. In contrast, oils from the Point Conception platforms can be distinguished on the basis of chemometric fingerprinting alone. In the middle of this spectrum are oils from Platforms Harmony, Heritage, and Hondo, where it is expected that oil weathering would take on the order of two weeks to a month to produce tarballs similar to those seen near Point Conception. In this case there is a much greater degree of weathering needed to proceed from produced oil to the biodegraded tar characteristic of tarball stranded on the beach. Tar deposition on beaches was monitored as part of cooperative with the County of Santa Barbara Energy Division and the U.S. Geological Survey during 2001-2003. We found tar deposition varies on a seasonal basis. In general, tarballs accumulate at a faster rate or remain longer on all beaches during the summer and fall months. The reasons for this are unclear based on our limited observations, however we speculate that factors such as prevailing winds and currents combined with more quiescent wave conditions favors the accumulation and preservation of tarballs on the beach during the summer and fall months. In contrast, winter storms, with much greater wave action remove beach sand and other materials, and stormy seas tend to break up oil that might weather into tarballs. Natural seepage is affected by the spring/neap tidal cycle; however, the link to tar deposition is unclear. Longer periods of monitoring are needed to address the variability in the data and provide a more robust statistical analysis.

California↗

Determination of total mercury in biological and geological samples

The analytical chemist is faced with several challenges when determining mercury in biological and geological materials. These challenges include widespread mercury contamination, both in the laboratory and the environment, possible losses of mercury during sample preparation and digestion, the wide range of mercury values commonly observed, ranging from the low nanogram per gram or per liter for background areas to hundreds of milligrams per kilogram in contaminated or ore-bearing areas, great matrix diversity, and sample heterogeneity1. These factors can be naturally occurring or anthropogenic, but must be addressed to provide a precise and accurate analysis. Although there are many instrumental methods available for the successful determination of mercury, no one technique will address all problems or all samples all of the time. The approach for the determination of mercury used at the U.S. Geological Survey, Crustal Imaging and Characterization Team, Denver Laboratories, utilizes a suite of complementary instrumental methods when approaching a study requiring mercury analyses. Typically, a study could require the analysis of waters, leachates or selective digestions of solids, vegetation, and biological materials such as tissue, bone, or shell, soils, rocks, sediments, coals, sludges, and(or) ashes. No one digestion or sample preparation method will be suitable for all of these matrices. The digestions typically employed at our laboratories include: (i) a closed-vessel microwave method using nitric acid and hydrogen peroxide, followed by digestion/dilution with a nitric acid/sodium dichromate solution, (ii) a robotic open test-tube digestion with nitric acid and sodium dichromate, (iii) a sealed Teflon? vessel with nitric acid and sodium dichromate, (iv) a sealed glass bottle with nitric acid and sodium dichromate, or (v) open test tube digestion with nitric and sulfuric acids and vanadium pentoxide. The common factor in all these digestions is that they are very oxidative to ensure the conversion of all mercury forms into Hg (II). Each method of digestion has its advantages and limitations. The method of detection used in our laboratories involves a combination of an in-house, custom, classic continuous-flow cold-vapor atomic absorption spectrometry (CVAAS), a commercially available, automated, flow-injection and a continuous flow cold-vapor atomic fluorescence spectrometry (CV-AFS) systems, and a relatively new, automated and integrated approach where solid or liquid samples are thermally decomposed under an oxygen atmosphere (a nitrogen atmosphere is used for coals) and the released mercury vapor trapped onto a gold gauze and then thermally released into an AAS system. Other less frequently used instrumental methods available for the determination of mercury include inductively coupled plasma ? optical emission spectrometry (ICP-OES), inductively couple plasma ? mass spectrometry (ICP-MS) (both solution nebulization and laser ablation), and instrumental neutron activation analysis (INAA). Results from two case studies involving the determination of mercury in the challenging matrices of biological materials will be presented. These will include fillet, liver and stomach-content samples from grayling for a baseline/background study in Alaska, and samples of meat tissue and shell material from Tanner crabs from Glacier Bay, Alaska. These studies show that the method of digestion is more important than a very sensitive detection limit for mercury.

Open-File Report↗

Multiple lines of evidence point to pesticides as stressors affecting invertebrate communities in small streams in five United States regions

Multistressor studies were performed in five regions of the United States to assess the role of pesticides as stressors affecting invertebrate communities in wadable streams. Pesticides and other chemical and physical stressors were measured in 75 to 99 streams per region for 4 weeks, after which invertebrate communities were surveyed (435 total sites). Pesticides were sampled weekly in filtered water, and once in bed sediment. The role of pesticides as a stressor to invertebrate communities was assessed by evaluating multiple lines of evidence: toxicity predictions based on measured pesticide concentrations, multivariate models and other statistical analyses, and previously published mesocosm experiments. Toxicity predictions using benchmarks and species sensitivity distributions and statistical correlations suggested that pesticides were present at high enough concentrations to adversely affect invertebrate communities at the regional scale. Two undirected techniques—boosted regression tree models and distance-based linear models—identified which pesticides were predictors of (respectively) invertebrate metrics and community composition. To put insecticides in context with known, influential covariates of invertebrate response, generalized additive models were used to identify which individual pesticide(s) were important predictors of invertebrate community condition in each region, after accounting for natural covariates. Four insecticides were identified as stressors to invertebrate communities at the regional scale: bifenthrin, chlordane, fipronil and its degradates, and imidacloprid. Fipronil was particularly important in the Southeast region, and imidacloprid, bifenthrin, and chlordane were important in multiple regions. For imidacloprid, bifenthrin, and fipronil, toxicity predictions were supported by mesocosm experiments that demonstrated adverse effects on naïve aquatic communities when dosed under controlled conditions. These multiple lines of evidence do not prove causality—which is challenging in the field under multistressor conditions—but they make a strong case for the role of insecticides as stressors adversely affecting invertebrate communities in streams within the five sampled regions.

Science of the Total Environment↗

A 16-year record of eolian dust in Southern Nevada and California, USA: Controls on dust generation and accumulation

An ongoing project monitors modern dust accumulation in the arid southwestern United States to (1) determine the rate and composition of dust inputs to soils and (2) relate dust accumulation to weather patterns to help predict the effects of climate change on dust production and accumulation. The 16-year records of 35 dust-trap sites in the eastern Mojave Desert and southern Great Basin reveal how generation and accumulation of dust, including the silt-clay, carbonate, and soluble-salt fractions, is affected by the amount and seasonal distribution of rainfall and the behavior of different source types (alluvium, dry playas, and wet playas). Accumulation rates (fluxes) of the silt-clay fraction of dust, including carbonates, range from about 2-20 g/m2/yr. Average rates are higher in the southern part of the study area (south of latitude 36.5??N) and annually fluctuate over a larger range than rates in the northern part of the area. Sites throughout the study area show peaks in dust flux in the 1984-1985 sampling period and again in 1997-1999; northern sites also show increased flux in 1987-1988 and southern sites in 1989-1991. These peaks of dust flux correspond with both La Nina (dry) conditions and with strong El Nino (wet) periods. The accumulation rates of different components of mineral dusts fluctuate differently. For example, soluble-salt flux increases in 1987-1988, coincident with a moderate El Nino event, and increases very strongly in 1997-1999, overlapping with a strong El Nino event. Both of these high-rainfall winters were preceded and accompanied by strong summer rains. In contrast, little or no change in soluble-salt flux occurred during other periods of high winter rainfall but little summer rain, e.g. 1992-1995. The differences between northern vs. southern sites and between sites with playa dust sources vs. alluvial dust sources indicate that regional differences in the response of precipitation and vegetation growth to ENSO influence and differences in the response of source types control dust production and accumulation. A major factor is the hydrologic condition of surface sediments. The silt-clay and soluble-salt fluxes increased during the El Nino events of 1987-1988 and 1997-1998 at sites close to "wet" playas with shallow depths to groundwater (<10 m), consistent with the concept that active evaporative concentration of salts disrupts surface crusts and increases the susceptibility of surface sediment to deflation. The silt-clay flux also increased during drought periods (1989-1991, 1995-1997) at sites downwind of alluvial sources and "dry" playas with deeper groundwater (<10 m). These increases are probably related to the die-off of drought-stressed vegetation on alluvial sediments, and in some cases to local runoff events that deliver fresh sediment to playa margins and distal portions of alluvial fans. ?? 2006 Elsevier Ltd. All rights reserved.

Journal of Arid Environments↗

Pesticide residues in birds and mammals

SUMMARY: Residues of organochlorine pesticides and their breakdown products are present in the tissues of essentially all wild birds throughout the world. These chemicals accumulate in fat from a relatively small environmental exposure. DDE and dieldrin are most prevalent. Others, such as heptachlor epoxide, chlordane, endrin, and benzene hexachloride also occur, the quantities and kinds generally reflecting local or regional use. Accumulation may be sufficient to kill animals following applications for pest control. This has occurred in several large-scale programmes in the United States. Mortality has also resulted from unintentional leakage of chemical from commercial establishments. Residues may persist in the environment for many years, exposing successive generations of animals. In general, birds that eat other birds, or fish, have higher residues than those that eat seeds and vegetation. The kinetic processes of absorption, metabolism, storage, and output differ according to both kind of chemical and species of animal. When exposure is low and continuous, a balance between intake and excretion may be achieved. Residues reach a balance at an approximate animal body equilibrium or plateau; the storage is generally proportional to dose. Experiments with chickens show that dieldrin and heptachlor epoxide have the greatest propensity for storage, endrin next, then DDT, then lindane. The storage of DDT was complicated by its metabolism to DDE and DDD, but other studies show that DDE has a much greater propensity for storage than either DDD or DDT. Methoxychlor has little cumulative capacity in birds. Residues in eggs reflect and parallel those in the parent bird during accumulation, equilibrium, and decline when dosage is discontinued. Residues with the greatest propensity for storage are also lost most slowly. Rate of loss of residues can be modified by dietary components and is speeded by weight loss of the animal. Under sublethal conditions of continuous exposure to an organochlorine pesticide, the concentrations of residues in the different tissues are ordinarily directly correlated with each other. When the dosage is at lethal levels, or when stored residues are mobilised to lethal levels, the balanced relationship is disrupted. The concentrations of residues in the brain provide the most rigorous criteria for diagnosis of death due to these chemicals, and levels are generally similar across a wide range of species of birds and mammals. Residues in liver are closely correlated with recent dose, either from direct intake or from mobilisation from storage, and so reflect hazardous exposure. Residues in the whole carcass show the storage reserve, and so indicate the potential for adverse effects from lethal mobilisation or from the continuous slow mobilisation that occurs during the normal processes of metabolism and excretion. A synchronous, rapid, and widespread decline in weight and thickness of shells of eggs laid by many species of wild birds occurred in the late 1940's and has persisted. Birds of prey were primarily affected; exceptions apparently are the result of lesser exposure because of different food habits. Many species of fish-eating birds are also affected. Others, however, appear to be more resistant and to accumulate much higher residues before shell-thinning occurs. Seed-eating birds do not appear to have been generally affected; their exposure is ordinarily lower, but physiological factors also seem to be involved. A relationship between shell-thinning and population decline has been established for many species. In exceptional cases, such as the herring gull, persistent re-nesting and other population reactions have overcome adverse effects at the population level. The discovery of shell-thinning among natural populations, and the hypothesis that this thinning was related to the occurrence of organochlorine pesticides, stimulated experimental studies to determine wheth

Book chapter↗

Potential minability and economic viability of the Antaramut-Kurtan-Dzoragukh coal field, north-central Armenia; a prefeasibility study

The U. S. Geological Survey (USGS) conducted a coal resource assessment of several areas in Armenia from 1997 to 1999. This report, which presents a prefeasibility study of the economic and mining potential of one coal deposit found and studied by the USGS team, was prepared using all data available at the time of the study and the results of the USGS exploratory work, including core drilling, trenching, coal quality analyses, and other ongoing field work. On the basis of information currently available, it is the authors? opinion that a small surface coal mine having about a 20-year life span could be developed in the Antaramut-Kurtan-Dzoragukh coal field, specifically at the Dzoragukh site. The mining organization selected or created to establish the mine will need to conduct necessary development drilling and other work to establish the final feasibility study for the mine. The company will need to be entrepreneurial, profit oriented, and sensitive to the coal consumer; have an analytical management staff; and focus on employee training, safety, and protection of the environment. It is anticipated that any interested parties will be required to submit detailed mining plans to the appropriate Armenian Government agencies. Further development work will be required to reach a final decision regarding the economic feasibility of the mine. However, available information indicates that a small, economic surface mine can be developed at this locality. The small mine suggested is a typical surface-outcropstripping, contour mining operation. In addition, auger mining is strongly suggested, because the recovery of these low-cost mining reserves will help to ensure that the operation will be a viable, economic enterprise. (Auger mining is a system in which large-diameter boreholes are placed horizontally into the coal seam at the final highwall set as the economic limit for the surface mining operation). A special horizontal boring machine, which can be imported from Russia, is required for auger mining. Although auger-mining coal reserves do exist, the necessary development work will further verify the extent of these reserves and all of the other indicated reserves. The following items are based on the detailed study reported in this publication. Initial investment.?Following an investment of US $85,000 over a 12-month period in mine development drilling and other activities, a decision must be taken regarding further investment in an ongoing mining operation. If the new data support the opening of the surface mine, __________________________ 1Consultant, 6024 Morning Dew Drive, Austin, TX 78749. 2 U.S. Geological Survey, 956 National Center, Reston, VA 20192 1 2 MINABILITY AND ECONOMIC VIABILITY, ANTARAMUT-KURTAN-DZORAGUKH COAL FIELD the $85,000 development cost is amortized over the first 10 years of mine production. If the new data do not support the opening of the mine, the $85,000 is considered a business development expense that may be written off against profits from other operations for income or other tax purposes or simply as a business loss. Total capital required.?The equipment costs will reach a total of $900,500 which will be amortized over a 7-year period to establish estimated coal mining costs. Estimated working capital costs are $300,000, which will be borrowed. Surface mining reserves.?Approximately 840,200 metric tonnes of surface minable coal reserves at 9.3 m3 of overburden per metric tonne of minable coal is indicated. Recovery of the minable coal at 85 percent will yield 714,000 recoverable metric tonnes of marketable as-mined coal. Auger mining reserves.?Auger-mining reserves of 576,000 metric tonnes are indicated. Recoverable auger-mining reserves of 202,000 metric tonnes (at 35-percent recovery) can be expected. Auger-mining production will vary according to the hole size being used, but, in either case, augering is a very profitable addition to the mining oper

Bulletin↗

Physical, chemical, and biological aspects of subsurface organic waste injection near Wilmington, North Carolina

From May 1968 to December 1972, an industrial organic waste was injected at rates of 100 to 200 gallons per minute (6.3 to 12.6 litres per second) into a sand, gravel, and limestone aquifer of Late Cretaceous age by Hercules Inc. located near Wilmington, North Carolina. This report presents both field and laboratory data pertaining to the physical, chemical, and biological effects of waste injection into the subsurface at this particular site, a case history of the operation, predictions of the reactions between certain organic wastes and the aquifer components, and descriptions of the effects of these reactions on the subsurface movement of the wastes. The case history documents a situation in which subsurface waste injection could not be considered a successful means of waste disposal. The first injection well was used only for 1 year due to excessive wellhead pressure build-up above the specified pressure limit of 150 pounds per square inch (10.3 bars). A second injection well drilled as a replacement operated for only 5 months before it too began to have problems with plugging. Upward leakage of waste into shallower aquifers was also detected at several wells in the injection-observation well system. The multiple problems of plugging, high pressures, and waste leakage suggested that the reactive nature of the waste with the aquifer into which it was injected was the primary reason for the difficulties experienced with waste injection. A site study was initiated in June 1971 to investigate waste-aquifer interactions. The first stage of the study determined the hydrogeologic conditions at the site, and characterized the industrial waste and the native ground water found in the injection zone and other aquifers. The injection zone consisted of multiple permeable zones ranging in depth from about 850 to 1,000 feet (259 to 305 metres) below land surface. In addition to the injection zone, aquifers were found near depths of 60, 300, 500, and 700 feet (18, 91, 152, and 213 metres) below land surface. The aquifers from 300 feet (91 metres) down to the injection zone were flowing artesian with the natural pressure of the injection zone being 65 feet (20 metres) above land surface at the site. The dissolved solids concentration in the native ground water increased with depth to an average value of 20,800 mg/l (milligram per litre) (two-thirds that of seawater) in the water from the injection zone. Sodium chloride was the major dissolved solid, and all of the ground water below 300-feet (91-metres) depth was slightly alkaline. Dissolved organic carbon of the industrial waste averaged 7,100 mg/l and 95 percent of the organic carbon was identified and quantified. The major organic waste constituents in order of decreasing abundance were acetic acid, formic acid, p-toluic acid, formaldehyde, methanol, terephthalic acid, phthalic acid, and benzoic acid. Prior to injection, the waste was neutralized with lime to pH 4 so that the major inorganic waste constituent was calcium at a concentration of 1,300 mg/l. The second stage of the site study involved the observation of waste-aquifer interactions at various wells as the waste arrived and passed by the wells. Water samples obtained from three observation wells located 1,500 to 2,000 feet (457 to 607 metres) from the original injection well gave evidence for biochemical waste transformations at low waste concentrations. Gas that effervesced from these water samples contained up to 54 percent methane by volume. Ferrous iron concentrations as high as 35 mg/l, hydrogen sulfide gas, and sulfide precipitates were additional indicators of biochemical reductive processes in the subsurface environment. Approximately 3,000 organisms per millilitre were found in uncontaminated ground water from the injection zone whereas in waste-contaminated wells, the number increased to levels as high as 1,000,000 organisms per millilitre. High concentrations of waste were found to be toxic to microo

Professional Paper↗

Selenium in ecosystems within the mountaintop coal mining and valley-fill region of southern West Virginia-assessment and ecosystem-scale modeling

Coal and associated waste rock are among environmental selenium (Se) sources that have the potential to affect reproduction in fish and aquatic birds. Ecosystems of southern West Virginia that are affected by drainage from mountaintop coal mines and valleys filled with waste rock in the Coal, Gauley, and Lower Guyandotte watersheds were assessed during 2010 and 2011. Sampling data from earlier studies in these watersheds (for example, Upper Mud River Reservoir) and other mining-affected watersheds also are included to assess additional hydrologic settings and food webs for comparison. Basin schematics give a comprehensive view of sampled species and Se concentration data specific to location and date. Food-web diagrams document the progression of Se trophic transfer across suspended particulate material, invertebrates, and fish for each site to serve as the basis for developing an ecosystem-scale model to predict Se exposure within the hydrologic conditions and food webs of southern West Virginia. This approach integrates a site-specific predator’s dietary exposure pathway into modeling to ensure an adequate link to Se toxicity and, thus, to species vulnerability. Site-specific fish abundance and richness data in streams documented various species of chub, shiner, dace, darters, bass, minnow, sunfish, sucker, catfish, and central stoneroller (Campostoma anomalum), mottled sculpin (Cottus bairdii), and least brook lamprey (Lampetra aepyptera). However, Se assessment species for streams, and hence, model species for streams, were limited to creek chub (Semotilus atromaculatus) and central stoneroller. Both of these species of fish are generally considered to have a high tolerance for environmental stress based on traditional comparative fish community assessment, with creek chub being present at all sites. Aquatic insects (mayfly, caddisfly, stonefly, dobsonfly, chironomid) were the main invertebrates sampled in streams. Collection of suspended particulate material acted as an integrator of organic-rich, fine-grained biomass present in streams. The base-case food web modeled for streams was suspended particulate material to aquatic insect to creek chub, with comparative modeling of a direct particulate-to-stoneroller food web. Model species for a reservoir setting were based on an earlier study of bluegill sunfish (Lepomis macrochirus), green sunfish (Lepomis cyanellus), and largemouth bass (Micropterus salmoides). Several reservoir food webs were considered based on a variety of invertebrates (insect, snail, clam). For stream and reservoir settings, predicted Se concentrations in exposure scenarios showed a high degree of correlation (r 2 = 0.91 for invertebrates and 0.75 for fish) with field observations of Se concentrations when modeling was initiated from suspended-particulate-material Se concentrations and model transfer parameters defined previously in the literature were used. These strong correlations validate the derived site-specific model and establish sufficient confidence that the predictions from the developed model can be quantitatively applied to the ecosystems in southern West Virginia. An application of modeling used a metric describing the partitioning of Se between particulate material and dissolved phases (K d ) to allow determination of a dissolved Se concentration that would be necessary to attain a site-specific Se fish body burden. The operationally defined K d quantifies the complex process of transformation at the base of a food web on a site-specific basis. The magnitude of this metric is known to vary with such factors as Se speciation, particulate-material type, and hydrology. This application (1) ties dissolved Se concentrations to fish tissue concentrations; (2) allows consideration of different choices for intervening site-specific exposure steps that set Se bioaccumulation, partitioning, and bioavailability; and (3) generates implications for management decisions that define protection through different regulatory pathways and guidelines. The range of model outcomes accounts for critical sources of variability and establishes whether site and food-web characterization were adequate to represent the dynamics of the system with certainty. This is especially true in terms of particulate-material phases at the base of the food web and utilization of K d in different hydrologic settings. For streams, a range of field-derived K d ds were applied to food-web exposure scenarios within a framework of locational and hydrologic variables (area of stream basin; stream gradient and discharge) that may affect the magnitude of K d . Overlaying even a coarse temporal scale that acknowledges variability in stream dissolved Se and Se speciation, such as through seasonal derivation of K d , can substantially narrow model uncertainty. Modeling that constrains the place and time of greatest ecosystem Se sensitivity within a specified food web gives insight into Se risk and identifies controlling management alternatives within a watershed or stream basin. If there is a range of hydrologic settings, specificity is needed to establish a hierarchy of in-stream and off-stream habitats for a watershed approach that takes into account Se-enriched water moving through different K d and food web environments. If there is a range of predator vulnerabilities (measured as a combination of food-web Se biodynamics and response in Se toxicity tests) within the site-specific community of fish species to be protected, then choice of fish species is critical to protection because it determines the food web and, hence, the magnitude of biotransfer through which Se is modeled. Whether creek chub is representative of the vulnerability to Se of all fish species encountered within the study-site ecosystems will require additional species-specific data and analysis. A range of site-specific scenarios illustrated here set model outcomes, but the final quantitative evaluation of alternatives and their implications will be those generated through choices and guidance formulated by state and other agencies in their decisionmaking processes. Proposed additions and refinements to the ecosystem-scale site-specific approach developed here include consideration of: measurement of temporally matched pairs of dissolved and suspended-particulate-material Se concentrations across a broader range of stream sites to expand the stream K d database and to test the representativeness of a suspended-particulate-material sample within a stream; characterization of different phases of particulate material across seasons to better define the base of the food web and connect to invertebrate feeding; refinement of model assumptions concerning dietary preferences and composition for fish to develop additional trophic transfer factors (TTFs) (for example, calculation of TTFinvertebrate composite for mixed diets); expansion of modeling of fish species and their food webs to include Se-vulnerable species; temporal characterization of a predator’s life cycle and habitat use as additional model layers to integrate with Se biodynamics in streams; investigation of the effect of stream gradient on K d based on a finer scale than presented here in terms of such variables as residence time, watershed dilution, and physical habitat attributes (for example, amount of ponding versus run or riffle within a stream); and linkage to discharge through use of stream gaging to record variability and enable model organization within water-year types and discharge seasons. Investigating the presence and variability of prey and predator species in demographically open systems such as streams also is key to model outcomes given the overall environmental stressors (for example, general landscape change, food-web disruption, recolonization potential) imposed on the composition of biological communities in coal mining and valley-fill affected watersheds

West Virginia↗

Spectral mixture analysis for surveillance of harmful algal blooms (SMASH): A field-, laboratory-, and satellite-based approach to identifying cyanobacteria genera from remotely sensed data

Algal blooms around the world are increasing in frequency and severity, often with the possibility of adverse effects on human and ecosystem health. The health and economic impacts associated with harmful algal blooms, or HABs, provide compelling rationale for developing new methods for monitoring these events via remote sensing. Although concentrations of chlorophyll- a and key pigments like phycocyanin are routinely estimated from satellite images and used to infer algal or cyanobacterial cell counts, current methods are unable to provide information on the taxonomic composition of a bloom. This study introduced a new approach capable of differentiating among genera based on their reflectance characteristics: Spectral Mixture Analysis for Surveillance of HABs, or SMASH. The foundation of SMASH is a multiple endmember spectral mixture analysis (MESMA) algorithm that takes a library of cyanobacteria endmembers and a hyperspectral image as input and estimates the fractional abundance of each genus, plus water, on a per-pixel basis. Importantly, we assume that the water column consists of only pure water and cyanobacteria, implying that our linear spectral unmixing models do not account for other optically active constituents such as suspended sediment and colored dissolved organic matter (CDOM). We used reflectance spectra for 12 genera measured under a microscope to populate an algal spectral library and applied the SMASH workflow to satellite images from four waterbodies across the United States. Normalized spectral separability scores indicated that the 12 genera were distinct from one another and the MESMA algorithm reproduced known input fractions for simulated mixtures that included all pairwise combinations of genera and water. We used Upper Klamath Lake as an example to illustrate data products generated via SMASH: maps of the normalized difference chlorophyll index and cyanobacterial index, a MESMA-based classification of algal genera, fraction images for each endmember, and a root mean square error (RMSE) image that summarizes uncertainty. For Upper Klamath Lake, these outputs highlighted a complex algal bloom featuring several genera, primarily Aphanizomenon , and intricate spatial patterns associated with gyres. The maximum RMSE constraint imposed on the MESMA algorithm provided a means of avoiding false positive detection of genera not present in a waterbody but must not be set so low as to leave much of an image unclassified in cases where genera included in the library are present. Comparison of endmember fractions with relative biovolumes calculated from field samples indicated that taxonomic information from SMASH was consistent with field observations. For example, the algorithm successfully identified Microcystis in Owasco Lake but avoided misclassifying Asterionella , a genus not yet included in our library, in Detroit Lake. This proof-of-concept investigation demonstrates the potential of SMASH to enhance our understanding of algal blooms, particularly with respect to their spatial and temporal dynamics.

New York, Oregon, Texas↗

New U-Pb geochronology and geochemistry of Paleozoic metaigneous rocks from western Yukon and eastern Alaska, cross-border synthesis, and implications for tectonic models

The tectonic evolution of and relation between the Yukon-Tanana terrane and the Lake George assemblage, as well as other associated tectonic assemblages in western Yukon and eastern Alaska, have been debated for decades. The Yukon-Tanana terrane is widely considered to be an allochthonous rifted fragment derived from the Laurentian continental margin, whereas the Lake George assemblage and associated assemblages are currently interpreted to be part of the parautochthonous continental margin of western North America (Laurentia). To address these topics, we present 40 new U-Pb zircon ages and 20 new whole-rock geochemical analyses. We incorporate these data into a new compilation of available geological mapping for a large area that straddles the Alaska-Yukon border, together with 34 previously published U-Pb age determinations and an extensive geochemical database of metaigneous rocks from Late Devonian to Early Mississippian and middle to late Permian assemblages in this area. Magmatism in the Lake George assemblage and related assemblages occurred in two pulses from about 371 to 360 and from about 358 to 347 million years ago (Ma); geochemical discrimination diagrams indicate a large crustal component, possibly indicative of arc magmatism, for felsic metaigneous rocks and a range of tectonic environments for mafic rocks. Magmatism in the Fortymile River and related assemblages, and parts of the Nasina assemblage—all parts of the Yukon-Tanana terrane—are mainly Early Mississippian and span a crystallization age range from about 361 to 343 Ma; geochemical discrimination diagrams for these rocks indicate primarily arc geochemical signatures for both mafic and felsic rocks. Middle to late Permian crystallization ages (about 261–253 Ma) are indicated for felsic metaigneous rocks in the Klondike assemblage and some of the felsic metaigneous rocks in the Nasina assemblage. Based on our mapping, we propose the existence of a possible unconformity between the Mississippian and Permian felsic metavolcanic rocks within the Nasina assemblage that is marked by sporadic occurrences of stretched-pebble conglomerate. Our combined database supports the well-established model of a magmatic arc comprising the Fortymile River and Finlayson assemblages of the rifted Yukon-Tanana terrane continental fragment on which a middle to late Permian arc (Klondike assemblage) was later built. The assemblages of the Yukon-Tanana terrane were subsequently intruded by Late Triassic to Early Jurassic granitoids, presumably during reaccretion of the Yukon-Tanana terrane to the continental margin. Permian and Late Triassic to Early Jurassic intrusions have not been mapped in the now structurally lower plate Lake George assemblage; their absence is one of the lines of evidence that have been used to support the parautochthonous, rather than allochthonous, origin of the Lake George assemblage and related assemblages. Our new data, together with previously published ranges of igneous crystallization ages and geochemical tectonic signatures of the Late Devonian to Early Mississippian magmatic rocks in the Lake George assemblage and associated assemblages and in the Fortymile River, Nasina, and correlated assemblages of the Yukon-Tanana terrane, indicate that the currently accepted interpretation of the Lake George assemblage and associated rocks being part of parauthochthonous North America is not the only possible interpretation of this tectonic entity. Approximately half of the dated intrusive rocks in the Lake George assemblage are contemporaneous with the metaigneous rocks of the Yukon-Tanana terrane arc (<361 Ma). We speculate that our approximately 361 Ma U-Pb age for quartz syenite in part of the North American continental margin in south-central Yukon defines the beginning of rifting of the Laurentian margin. Although the currently favored model of prolonged middle Paleozoic subduction and extension in both the Yukon-Tanana terrane and parautochthonous North America allows for simultaneous middle Paleozoic magmatism on both sides of the Slide Mountain Ocean, we now propose an alternative hypothesis in which the Lake George assemblage represents a deeper part of the rifted Yukon-Tanana terrane arc. If this is the case, the absence of Permian and Late Triassic to Early Jurassic arc rocks in the Lake George assemblage could be explained either by the arcs of these ages not being wide enough to have affected the Lake George assemblage or by tectonic displacement of these arc rocks away from the Lake George assemblage. Our approximately 259 Ma U-Pb zircon age and geochemical analyses of metarhyolite in the Seventymile terrane in Alaska, which comprises remnants of the back-arc basin that separated the Yukon-Tanana terrane from the Laurentian continental margin, confirm the presence of a late middle Permian volcanic arc component to the terrane. Our approximately 319 Ma U-Pb zircon age from the Chicken assemblage (as redefined in this study) in eastern Alaska, combined with previously reported fossil ages and a U-Pb zircon age from this assemblage, indicate that it is a Late Mississippian to Early Pennsylvanian arc assemblage. We propose several other relatively young, locally developed arc assemblages outboard of the ancient continental margin of Laurentia that may correlate with the Chicken assemblage, but we consider its origin to remain an enigma.

Alaska↗

SAR studies in the Yuma Desert, Arizona: Sand penetration, geology, and the detection of military ordnance debris

Synthetic Aperture Radar (SAR) images acquired over part of the Yuma Desert in southwestern Arizona demonstrate the ability of C-band (5.7-cm wavelength), L-band (24.5 cm), and P-band (68 cm) AIRSAR signals to backscatter from increasingly greater depths reaching several meters in blow sand and sandy alluvium. AIRSAR images obtained within the Barry M. Goldwater Bombing and Gunnery Range near Yuma, Arizona, show a total reversal of C- and P-band backscatter contrast (image tone) for three distinct geologic units. This phenomenon results from an increasingly greater depth of radar imaging with increasing radar wavelength. In the case of sandy- and small pebble-alluvium surfaces mantled by up to several meters of blow sand, backscatter increases directly with SAR wavelength as a result of volume scattering from a calcic soil horizon at shallow depth and by volume scattering from the root mounds of healthy desert vegetation that locally stabilize blow sand. AIRSAR images obtained within the military range are also shown to be useful for detecting metallic military ordnance debris that is located either at the surface or covered by tens of centimeters to several meters of blow sand. The degree of detectability of this ordnance increases with SAR wavelength and is clearly maximized on P-band images that are processed in the cross-polarized mode (HV). This effect is attributed to maximum signal penetration at P-band and the enhanced PHV image contrast between the radar-bright ordnance debris and the radar-dark sandy desert. This article focuses on the interpretation of high resolution AIRSAR images but also Compares these airborne SAR images with those acquired from spacecraft sensors such as ERS-SAR and Space Radar Laboratory (SIR-C/X-SAR).Synthetic Aperture Radar (SAR) images acquired over part of the Yuma Desert in southwestern Arizona demonstrate the ability of C-band (5.7-cm wavelength), L-band (24.5 cm), and P-band (68 cm) AIRSAR signals to backscatter from increasingly greater depths reaching several meters in blow sand and sandy alluvium. AIRSAR images obtained within the Barry M. Goldwater Bombing and Gunnery Range near Yuma, Arizona, show a total reversal of C- and P-band backscatter contrast (image tone) for three distinct geologic units. This phenomenon results from an increasingly greater depth of radar imaging with increasing radar wavelength. In the case of sandy- and small pebble-alluvium surfaces mantled by up to several meters of blow sand, backscatter increases directly with SAR wavelength as a result of volume scattering from a calcic soil horizon at shallow depth and by volume scattering from the root mounds of healthy desert vegetation that locally stabilize blow sand. AIRSAR images obtained within the military range are also shown to be useful for detecting metallic military ordnance debris that is located either at the surface or covered by tens of centimeters to several meters of blow sand. The degree of detectability of this ordnance increases with SAR wavelength and is clearly maximized on P-band images that are processed in the cross-polarized mode (HV). This effect is attributed to maximum signal penetration at P-band and the enhanced PHV image contrast between the radar-bright ordnance debris and the radar-dark sandy desert. This article focuses on the interpretation of high resolution AIRSAR images but also compares these airborne SAR images with those acquired from spacecraft sensors such as ERS-SAR and Space Radar Laboratory (SIR-C/X-SAR).

Remote Sensing of Environment↗

Occurrence of fecal-indicator bacteria and protocols for identification of fecal-contamination sources in selected reaches of the West Branch Brandywine Creek, Chester County, Pennsylvania

The presence of fecal-indicator bacteria indicates the potential presence of pathogens originating from the fecal matter of warm-blooded animals. These pathogens are responsible for numerous human diseases ranging from common diarrhea to meningitis and polio. The detection of fecal-indicator bacteria and interpretation of the resultant data are, therefore, of great importance to water-resource managers. Current (2005) techniques used to assess fecal contamination within the fluvial environment primarily assess samples collected from the water column, either as grab samples or as depth- and (or) width-integrated samples. However, current research indicates approximately 99 percent of all bacteria within nature exist as attached, or sessile, bacteria. Because of this condition, most current techniques for the detection of fecal contamination, which utilize bacteria, assess only about 1 percent of the total bacteria within the fluvial system and are, therefore, problematic. Evaluation of the environmental factors affecting the occurrence and distribution of bacteria within the fluvial system, as well as the evaluation and modification of alternative approaches that effectively quantify the larger population of sessile bacteria within fluvial sediments, will present water-resource managers with more effective tools to assess, prevent, and (or) eliminate sources of fecal contamination within pristine and impaired watersheds. Two stream reaches on the West Branch Brandywine Creek in the Coatesville, Pa., region were studied between September 2002 and August 2003. The effects of sediment particle size, climatic conditions, aquatic growth, environmental chemistry, impervious surfaces, sediment and soil filtration, and dams on observed bacteria concentrations were evaluated. Alternative approaches were assessed to better detect geographic sources of fecal contamination including the use of turbidity as a surrogate for bacteria, the modification and implementation of sandbag bacteria samplers, and the use of optical brighteners. For the purposes of this report, sources of bacteria were defined as geographic locations where elevated concentrations of bacteria are observed within, or expected to enter, the main branch of the West Branch Brandywine Creek. Biologic sources (for example, waterfowl) were noted where applicable; however, no specific study of biologic sources (such as bacterial source tracking) was conducted. Data indicated that specific bacterial populations within fluvial sediments could be related to specific particle-size ranges. This relation is likely the result of the reduced porosity and permeability associated with finer sediments and the ability of specific bacteria to tolerate particular environments. Escherichia coli (E. coli) showed a higher median concentration (2,160 colonies per gram of saturated sediment) in the 0.125 to 0.5-millimeter size range of natural sediments than in other ranges, and enterococcus bacteria showed a higher median concentration (61,830 colonies per gram of saturated sediment) in the 0.062 to 0.25-millimeter size range of natural sediments than in other ranges. There were insufficient data to assess the particle-size relation to fecal coliform bacteria and (or) fecal streptococcus bacteria. Climatic conditions were shown to affect bacteria concentrations in both the water column and fluvial sediments. Drought conditions in 2002 resulted in lower overall bacteria concentrations than the more typically wet year of 2003. E. coli concentrations in fluvial sediment along the Coatesville study reach in 2002 had a median concentration of 92 colonies per gram of saturated sediment; in 2003, the median concentration had risen to 4,752 colonies per gram of saturated sediment. Symbiotic relations between bacteria and aquatic growth were likely responsible for increased bacteria concentrations observed within an impoundment area on the Coatesville study reach. This reach showed evidence of elevated aquatic growth and sharp increases in E. coli concentrations from upstream to downstream through the impoundment area in both 2002 and 2003. In 2003, E. coli concentrations within the waters column increased from 940 colonies per 100 milliliters upstream to 6,000 colonies per 100 milliliters at the dam crest. Given that these bacteria likely resulted from natural bacterial regrowth, the use of E. coli as an indicator of fecal contamination was severely impaired. Variable environmental conditions along the West Branch Brandywine Creek made the common field-chemical parameters of specific conductance, temperature, pH, and dissolved oxygen ineffective and (or) impossible to use for the determination of inputs of fecal contamination. Extreme variations in chemical gradients commonly were related to the urban/industrial signature of the watershed. For example, during base-flow sampling in 2002, specific-conductance values exceeding 1,000 microsiemens per centimeter observed in effluent from a local steel mill. This effluent raised the specific conductance within the West Branch Brandywine from just above 200 microsiemens per centimeter upstream from the outfall to just below 500 microsiemens per centimeter downstream from the outfall. These chemical gradients also, likely, had an effect on the initial colonization of bacteria, the formation of biofilms, and the persistence of certain types of bacteria along the study reach. Data collected in 2003 indicated that nutrients increased during both base-flow and stormflow conditions along the Coatesville study reach. For example, during base-flow sampling in 2003, 20 pounds of phosphorus was shown to enter the West Branch Brandywine Creek along the Coatesville study reach. The largest contributors to this base-flow nutrient load were likely two wastewater-treatment facilities adjacent to the study reach. During stormflow sampling in 2003, 480 pounds of phosphorus was shown to enter the West Branch Brandywine Creek along the Coatesville study reach. Data, along with other research, indicated the largest contributor to this stormflow nutrient load was likely remobilized sediment originating from a large dam impoundment. These elevated nutrient concentrations were considered sufficient to promote accelerated aquatic growth along the reach. Data collected in 2003 showed that wastewater constituents entered the West Branch Brandywine Creek largely from urban storm-sewer systems. Samples from the primary storm sewer for the city of Coatesville had detections for 20 of 69 wastewater constituents. These constituents included both strong and weak indicators of fecal contamination and generally indicated the storm-sewer system along the Coatesville study reach was a likely source of fecal-indicator bacteria and fecal contamination under base-flow conditions. By comparison, 5 constituents were detected in samples from the upstream end of the reach, and 10 constituents were detected in samples from the downstream end of the reach. During stormflow, numbers of detections were similar along the entire length of the study reach-five in samples from the upstream end, eight in samples from the center of the reach, and seven in samples from the downstream end of the reach. These data indicate that point sources (such as culverts and pipes, septic systems, and wastewater-treatment facilities) are not likely the origin of bacteria contamination during stormflow. The bacteria concentrations observed during stormflow events probably result from remobilized sessile bacteria stored within fluvial sediments. In this case, these bacteria should not be considered indicators of current fecal contamination. Impervious surfaces were found to increase bacteria concentrations along the West Branch Brandywine Creek because contaminated runoff from impervious areas generally flows into, and is concentrated within, the confines of the local storm-sewer system. During 2002, storm-sewer outfalls draining impervious areas were associated with all major locations of elevated bacterial concentrations (greater than 1,200 colonies per gram of saturated sediment) in fluvial sediments. During 2003, wetter conditions and overall bacteria concentrations higher than in 2002 resulted in point sources of bacterial contamination becoming less pronounced; however, the storm-sewer system, draining adjacent impervious areas, was still observed to be the primary source of bacteria along the reach. Where stormwater and (or) other runoff from these areas was allowed to infiltrate and (or) flow through wetland and riparian buffers, bacteria concentrations were not observed to be elevated above background levels commonly observed throughout similar areas of the same reach. Two run-of-the-river dams along the Coatesville study reach were evaluated for their effects on observed bacterial concentrations. These dams were shown to have greater or lesser effects on bacterial concentrations depending on the size of the structure and the capacity of the structure to impede flows. The smaller upstream dam had an approximate height of 3 feet and showed little observed effect on measured turbidity values; these data indicated that the dam did not effectively impede the flow of water or sediment within the West Branch Brandywine Creek. Consequently, this small dam did not show any observed effect on bacterial concentrations either upstream or downstream of the structure. The larger dam, near the middle of the reach, had an approximate height of 20 feet and showed greater effects on both turbidity and bacteria concentrations. The capacity of the larger dam to impede flows, combined with nutrients entering the reach, resulted in increased biologic activity throughout the impoundment area. Within this larger impoundment, enterococcus bacteria populations were observed to decrease sharply and E. coli bacteria populations were observed to increase sharply as flow approached the dam crest. All bacteria levels were then observed to drop to background levels, in both the water column and fluvial sediment, immediately downstream from the dam crest. Additional study is required to determine the cause for this rapid die off. Turbidity was assessed as a potential surrogate for E. coli bacteria. Regression analysis indicated higher turbidity levels usually can indicate higher concentrations of bacteria (R2 = 0.67), but the relation was too sporadic on the West Branch Brandywine Creek to use turbidity as a surrogate for estimated bacteria concentrations. Evaluation of data from individual base-flow and stormflow events resulted in variable and generally poor statistical relations between E. coli bacteria and turbidity (R2 values ranged from 0.02 to 0.94). Sandbag samplers were used in 2003 to determine their suitability for the assessment of fecal contamination. Sandbag samplers rely on the ability of bacteria to attach to surfaces and use the larger sessile bacteria populations instead of the more commonly used planktonic bacteria populations. E. coli bacteria concentrations observed in the sandbag samplers, after 1 week in place, were similar to those found within natural sediments collected concurrently. Enterococcus bacteria concentrations within the same sandbag samplers were not similar, and were generally lower, than those observed within the natural sediments. This discrepancy was likely because sand within the samplers was sieved to a size that was likely too coarse for enterococcus bacteria to persist. Optical-brightener samplers were installed along with each sandbag sampler. Optical brighteners are additives used in common household detergents; therefore, detection of optical brighteners, along with elevated fecal-indicator bacteria concentrations, strongly indicates a link to humans. Positive results for optical brighteners were detected only at the outfalls of two sewage-treatment facilities; because of treatment of the effluent from these facilities, these samples did not have elevated bacteria concentrations. The lack of additional positive results was largely because this method is not sensitive to low concentrations of optical brighteners.

Pennsylvania↗

Flood-frequency analyses from paleoflood investigations for Spring, Rapid, Boxelder, and Elk Creeks, Black Hills, western South Dakota

Flood-frequency analyses for the Black Hills area are important because of severe flooding of June 9-10, 1972, that was caused by a large mesoscale convective system and caused at least 238 deaths. Many 1972 peak flows are high outliers (by factors of 10 or more) in observed records that date to the early 1900s. An efficient means of reducing uncertainties for flood recurrence is to augment gaged records by using paleohydrologic techniques to determine ages and magnitudes of prior large floods (paleofloods). This report summarizes results of paleoflood investigations for Spring Creek, Rapid Creek (two reaches), Boxelder Creek (two subreaches), and Elk Creek. Stratigraphic records and resulting long-term flood chronologies, locally extending more than 2,000 years, were combined with observed and adjusted peak-flow values (gaged records) and historical flood information to derive flood-frequency estimates for the six study reaches. Results indicate that (1) floods as large as and even substantially larger than 1972 have affected most of the study reaches, and (2) incorporation of the paleohydrologic information substantially reduced uncertainties in estimating flood recurrence. Canyons within outcrops of Paleozoic rocks along the eastern flanks of the Black Hills provided excellent environments for (1) deposition and preservation of stratigraphic sequences of late-Holocene flood deposits, primarily in protected slack-water settings flanking the streams; and (2) hydraulic analyses for determination of associated flow magnitudes. The bedrock canyons ensure long-term stability of channel and valley geometry, thereby increasing confidence in hydraulic computations of ancient floods from modern channel geometry. Stratigraphic records of flood sequences, in combination with deposit dating by radiocarbon, optically stimulated luminescence, and cesium-137, provided paleoflood chronologies for 29 individual study sites. Flow magnitudes were estimated from elevations of flood deposits in conjunction with hydraulic calculations based on modern channel and valley geometry. Reach-scale paleoflood chronologies were interpreted for each study reach, which generally entailed correlation of flood evidence among multiple sites, chiefly based on relative position within stratigraphic sequences, unique textural characteristics, or results of age dating and flow estimation. The FLDFRQ3 and PeakfqSA analytical models (assuming log-Pearson Type III frequency distributions) were used for flood-frequency analyses for as many as four scenarios: (1) analysis of gaged records only; (2) gaged records with historical information; (3) all available data including gaged records, historical flows, paleofloods, and perception thresholds; and (4) the same as the third scenario, but ?top fitting? the distribution using only the largest 50 percent of gaged peak flows. The PeakfqSA model is most consistent with procedures adopted by most Federal agencies for flood-frequency analysis and thus was (1) used for comparisons among results for study reaches, and (2) considered by the authors as most appropriate for general applications of estimating low-probability flood recurrence. The detailed paleoflood investigations indicated that in the last 2,000 years all study reaches have had multiple large floods substantially larger than in gaged records. For Spring Creek, stratigraphic records preserved a chronology of at least five paleofloods in approximately (~) 1,000 years approaching or exceeding the 1972 flow of 21,800 cubic feet per second (ft 3 /s). The largest was ~700 years ago with a flow range of 29,300-58,600 ft 3 /s, which reflects the uncertainty regarding flood-magnitude estimates that was incorporated in the flood-frequency analyses. In the lower reach of Rapid Creek (downstream from Pactola Dam), two paleofloods in ~1,000 years exceeded the 1972 flow of 31,200 ft 3 /s. Those occurred ~440 and 1,000 years ago, with flows of 128,000-256,000 and 64,000-128,000 ft 3 /s, respectively. Five smaller paleofloods of 9,500-19,000 ft 3 /s occurred between ~200 and 400 years ago. In the upper reach of Rapid Creek (above Pactola Reservoir), the largest recorded floods are substantially smaller than for lower Rapid Creek and all other study reaches. Paleofloods of ~12,900 and 12,000 ft 3 /s occurred ~1,000 and 1,500 years ago. One additional paleoflood (~800 years ago) was similar in magnitude to the largest gaged flow of 2,460 ft 3 /s Boxelder Creek was treated as having two subreaches because of two tributaries that affect peak flows. During the last ~1,000 years, paleofloods of ~39,000-78,000 ft 3 /s and 40,000-80,000 ft 3 /s in the upstream subreach have exceeded the 1972 peak flow of 30,800 ft 3 /s. One other paleoflood was similar to the second largest gaged flow (16,400 ft 3 /s in 1907). For the downstream subreach, paleofloods of 61,300-123,000 ft 3 /s and 52,500-105,000 ft 3 /s in the last ~1,000 years have substantially exceeded the 1972 flood (50,500 ft 3 /s). Four additional paleofloods had flows between 14,200 and 33,800 ft 3 /s. The 1972 flow on Elk Creek (10,400 ft 3 /s) has been substantially exceeded at least five times in the last 1,900 years. The largest paleoflood (41,500-124,000 ft 3 /s) was ~900 years ago. Three other paleofloods between 37,500 and 120,000 ft 3 /s occurred between 1,100 and 1,800 years ago. A fifth paleoflood of 25,500-76,500 ft 3 /s was ~750 years ago. Considering analyses for all available data (PeakfqSA model) for all six study reaches, the 95-percent confidence intervals about the low-probability quantile estimates (100-, 200-, and 500-year recurrence intervals) were reduced by at least 78 percent relative to those for the gaged records only. In some cases, 95-percent uncertainty intervals were reduced by 99 percent or more. For all study reaches except the two Boxelder Creek subreaches, quantile estimates for these long-term analyses were larger than for the short-term analyses. The 1972 flow for the Spring Creek study reach (21,800 ft 3 /s) corresponds with a recurrence interval of ~400 years. Recurrence intervals are ~500 years for the 1972 flood magnitudes along the lower Rapid Creek reach and the upstream subreach of Boxelder Creek. For the downstream subreach of Boxelder Creek, the large 1972 flood magnitude (50,500 ft 3 /s) exceeds the 500-year quantile estimate by about 35 percent. The recurrence interval of ~100 years for 1972 flooding along the Elk Creek study reach is small relative to other study reaches along the eastern margin of the Black Hills. All of the paleofloods plot within the bounds of a national envelope curve, indicating that the national curve represents exceedingly rare floods for the Black Hills area. Elk Creek, lower Rapid Creek, and the downstream subreach of Boxelder Creek all have paleofloods that plot above a regional envelope curve; in the case of Elk Creek, by a factor of nearly two. The Black Hills paleofloods represent some of the largest known floods, relative to drainage area, for the United States. Many of the other largest known United States floods are in areas with physiographic and climatologic conditions broadly similar to the Black Hills-semiarid and rugged landscapes that intercept and focus heavy precipitation from convective storm systems. The 1972 precipitation and runoff patterns, previous analyses of peak-flow records, and the paleoflood investigations of this study support a hypothesis of distinct differences in flood generation within the central Black Hills study area. The eastern Black Hills are susceptible to intense orographic lifting associated with convective storm systems and also have high relief, thin soils, and narrow and steep canyons-factors favoring generation of exceptionally heavy rain-producing thunderstorms and promoting runoff and rapid concentration of flow into stream channels. In contrast, storm potential is smaller in and near the Limestone Plateau area, and storm runoff is further reduced by substantial infiltration into the limestone, gentle topography, and extensive floodplain storage. Results of the paleoflood investigations are directly applicable only to the specific study reaches and in the case of Rapid Creek, only to pre-regulation conditions. Thus, approaches for broader applications were developed from inferences of overall flood-generation processes, and appropriate domains for application of results were described. Example applications were provided by estimating flood quantiles for selected streamgages, which also allowed direct comparison with results of at-site flood-frequency analyses from a previous study. Several broad issues and uncertainties were examined, including potential biases associated with stratigraphic records that inherently are not always complete, uncertainties regarding statistical approaches, and the unknown applicability of paleoflood records to future watershed conditions. The results of the paleoflood investigations, however, provide much better physically based information on low-probability floods than has been available previously, substantially improving estimates of the magnitude and frequency of large floods in these basins and reducing associated uncertainty.

South Dakota↗