USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geologic Notes”

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.

701 records · Page 39Linked to original sources

Lateral and vertical channel movement and potential for bed-material movement on the Madison River downstream from Earthquake Lake, Montana

The 1959 Hebgen Lake earthquake caused a massive landslide (Madison Slide) that dammed the Madison River and formed Earthquake Lake. The U.S. Army Corps of Engineers excavated a spillway through the Madison Slide to permit outflow from Earthquake Lake. In June 1970, high streamflows on the Madison River severely eroded the spillway channel and damaged the roadway embankment along U.S. Highway 287 downstream from the Madison Slide. Investigations undertaken following the 1970 flood events concluded that substantial erosion through and downstream from the spillway could be expected for streamflows greater than 3,500 cubic feet per second (ft 3 /s). Accordingly, the owners of Hebgen Dam, upstream from Earthquake Lake, have tried to manage releases from Hebgen Lake to prevent streamflows from exceeding 3,500 ft 3 /s measured at the U.S. Geological Survey (USGS) gaging station 0638800 Madison River at Kirby Ranch, near Cameron, Montana. Management of flow releases from Hebgen Lake to avoid exceeding the threshold streamflow at USGS gaging station 06038800 is difficult, and has been questioned for two reasons. First, no road damage was reported downstream from the Earthquake Lake outlet in 1993, 1996, and 1997 when streamflows exceeded the 3,500-ft 3 /s threshold. Second, the 3,500-ft 3 /s threshold generally precludes releases of higher flows that could be beneficial to the blue-ribbon trout fishery downstream in the Madison River. In response to concerns about minimizing streamflow downstream from Earthquake Lake and the possible armoring of the spillway, the USGS, in cooperation with the Madison River Fisheries Technical Advisory Committee (MADTAC; Bureau of Land Management; Montana Department of Environmental Quality; Montana Fish, Wildlife and Parks; PPL-Montana; U.S. Department of Agriculture Forest Service - Gallatin National Forest; and U.S. Fish and Wildlife Service), conducted a study to determine movement of the Madison River channel downstream from Earthquake Lake and to investigate the potential for bed material movement along the same reach. The purpose of this report is to present information about the lateral and vertical movement of the Madison River from 1970 to 2006 for a 1-mile reach downstream from Earthquake Lake and for Raynolds Pass Bridge, and to provide an analysis of the potential for bed-material movement so that MADTAC can evaluate the applicability of the previously determined threshold streamflow for initiation of damaging erosion. As part of this study channel cross sections originally surveyed by the USGS in 1971 were resurveyed in 2006. Incremental channel-movement distances were determined by comparing the stream centerlines from 14 aerial photographs taken between 1970 and 2006. Depths of channel incision and aggregation were determined by comparing the 2006 and 1971 cross-section and water-surface data. Particle sizes of bed and bank materials were measured in 2006 and 2008 using the pebble-count method and sieve analyses. A one-dimensional hydraulic-flow model (HEC-RAS) was used to calculate mean boundary-shear stresses for various streamflows; these calculated boundary-shear stresses were compared to calculated critical-shear stresses for the bed materials to determine the potential for bed-material movement. A comparison of lateral channel movement distances with annual peak streamflows shows that streamflows higher than the 3,500-ft 3 /s threshold were followed by lateral channel movement except from 1991 to 1992 and possibly from 1996 to 1997. However, it was not possible to discern whether the channel moved gradually or suddenly, or in response to one peak flow, to several peak flows, or to sustained flows. The channel moved between 2002 and 2005 even when streamflows were less than the threshold streamflow of 3,500 ft 3 /s. Comparisons of cross sections and aerial photographs show that the channel has moved laterally and incised and aggraded to varying degrees. The channel has developed meander bends and has incised as much as 5–12 feet (ft) through the upstream part of the Madison Slide (cross sections 1400–800). Near cross section 800, the stream has eroded into the steep right bank between the stream and the road where fill was mechanically placed after 1970. Channel movement also was noted downstream from the Madison Slide. Near Raynolds Pass Bridge, about 3 miles (mi) downstream from Earthquake Lake, elevations across the channel have changed by -1.4 ft to +1.9 ft, but these changes were local in nature and could represent a few rocks or depressions in the bed. Overall, it does not appear that the materials eroded from the Madison Slide are causing aggradation in the subreach near the Raynolds Pass Bridge. Comparisons of critical shear stresses to mean boundary-shear stresses indicate that the D50 particle sizes (median size) along the right side of the bed between cross sections 400 and 500 and along the right side of the bed between cross sections 1300 and 1400 could move at the threshold streamflow. In contrast, most of the D84 particle sizes at those two locations probably will not move at the threshold streamflow. This lack of movement for the larger particles at the threshold streamflow could lead to further armoring of the bed as the D50 and smaller-sized particles are removed from the bed and transported downstream. The Shields parameter values from 0.04 to 0.08 that were used to calculate critical shear stresses could be conservative for a high-gradient stream such as the Madison. A higher, less conservative, Shields parameter would result in higher critical shear stresses, meaning that higher streamflows would be required to move material than those reported herein. In addition, because materials in the channel thalweg are exposed to higher boundary-shear stresses than the materials along the sides of the channel, larger, more erosion-resistant materials likely exist in the deeper parts of the channel where high-flow depths and velocities prevented sediment sampling. Movement of these materials might require higher critical shear stresses than estimated in this report. Characterization of sediment sizes in the center of the stream and observation of bed-material movement for a range of streamflows could provide information to help refine the Shields parameter and critical-shear stress estimates for bed materials in the Madison River downstream from Earthquake Lake. Furthermore, resurveying cross sections and water-surface elevations more frequently (either annually or after high streamflows) could better define the relation between streamflow and lateral and vertical channel movement.

Montana↗

Aquifers of Arkansas: protection, management, and hydrologic and geochemical characteristics of groundwater resources in Arkansas

Sixteen aquifers in Arkansas that currently serve or have served as sources of water supply are described with respect to existing groundwater protection and management programs, geology, hydrologic characteristics, water use, water levels, deductive analysis, projections of hydrologic conditions, and water quality. State and Federal protection and management programs are described according to regulatory oversight, management strategies, and ambient groundwater-monitoring programs that currently (2013) are in place for assessing and protecting groundwater resources throughout the State. Physical attributes, groundwater geochemistry, and groundwater quality are described for each of the 16 aquifers of the State. Information in regard to the hydrology and geochemistry of each of the aquifers is summarized from about 550 historical and recent publications. Additionally, more than 8,000 sites with groundwater-quality data were obtained from the U.S. Geological Survey National Water Information System and the Arkansas Department of Environmental Quality databases and entered into a spatial database to investigate distribution and trends in chemical constituents for each of the aquifers. The 16 aquifers of the State were divided into two major physiographic regions of the State: the Coastal Plain Province (referred to as Coastal Plain) of eastern and southern Arkansas, which includes 11 of the 16 aquifers, and the Interior Highlands Division (referred to as Interior Highlands) of western Arkansas, which includes the remaining 5 aquifers. The 11 aquifers in the Coastal Plain consist of various geologic units that are Cenozoic in age and consist primarily of Cretaceous, Tertiary, and Quaternary sands, gravels, silts, and clays. Groundwater in the Coastal Plain represents one of the most valuable natural resources in the State, driving the economic engines of agriculture, while also supplying abundant water for commercial, industrial, and public-supply use. In terms of age from youngest to oldest, the aquifers of the Coastal Plain include Quaternary alluvial aquifers, including the Mississippi River Valley alluvial aquifer (the most important aquifer in Arkansas in terms of volume of use and economic benefits), the Jackson Group (a regional confining unit that served for decades as an important source of domestic supply), and the Cockfield, Sparta, Cane River, Carrizo, Wilcox, Nacatoch, Ozan, Tokio, and Trinity aquifers. The Mississippi River Valley alluvial aquifer accounts for approximately 94 percent of all groundwater used in the State, and the aquifer is used primarily for irrigation purposes. The Sparta aquifer is the second most important aquifer in terms of use, and the aquifer was used in the past dominantly as a source of public and industrial supply, although increasing irrigation use is occurring because of critically declining water levels in the Mississippi River Valley alluvial aquifer. Other aquifers of the Coastal Plain generally are used as important local sources of domestic, industrial, and public supply, in addition to other minor uses. Water quality generally is good for all aquifers of the Coastal Plain, except for elevated iron concentrations and localized areas of high salinity. The high salinity results from intrusion from underlying formations, evapotranspiration processes in areas of low recharge, and inadequate flushing in downgradient areas of residual salinity from deposition in marine environments. Trends in the spatial distribution of individual chemical constituents are related to position along the flow path for most aquifers of the Coastal Plain. These trends include elevated iron and nitrate concentrations with lower pH values and dissolved solids in groundwater from the outcrop areas, transitioning to lower iron and nitrate (related to changes in redox) and higher pH and dissolved solids (dominantly from the dissolution of carbonate minerals) in groundwater downgradient from outcrop areas. Groundwater generally trended from a calcium- to a sodium-bicarbonate water type with increasing cation exchange along the flow path. The Interior Highlands of western Arkansas has less reported groundwater use than other areas of the State, reflecting a combination of factors. These factors include prevalent and increasing use of surface water, less intensive agricultural uses, lower population and industry densities, lesser potential yield of the resource, and lack of detailed reporting. The overall low yields of aquifers of the Interior Highlands result in domestic supply as the dominant use, with minor industrial, public, and commercial-supply use. Where greater volumes are required for growth of population and industry, surface water is the greatest supplier of water needs in the Interior Highlands. The various aquifers of the Interior Highlands generally occur in shallow, fractured, well-indurated, structurally modified bedrock of this mountainous region of the State, as compared to the relatively flat-lying, unconsolidated sediments of the Coastal Plain. In terms of age from youngest to oldest, the aquifers of the Interior Highlands include: the Arkansas River Valley alluvial aquifer, the Ouachita Mountains aquifer, the Western Interior Plains confining system, the Springfield Plateau aquifer, and the Ozark aquifer. Spatial trends in groundwater geochemistry in the Interior Highlands differ greatly from trends noted for aquifers of the Coastal Plain. In the Coastal Plain, the prevalence of long regional flow paths results in regionally predictable and mappable geochemical changes along the flow paths. In the Interior Highlands, short, topographically controlled flow paths (from hilltops to valleys) within small watersheds represent the predominant groundwater-flow system. As such, dense data coverage from numerous wells would be required to effectively characterize these groundwater basins and define small-scale geochemical changes along any given flow path for aquifers of the Interior Highlands. Changes in geochemistry generally were related to rock type and residence time along individual flow paths. Dominant changes in geochemistry for the Ouachita Mountains aquifer and the Western Interior Plains confining system are attributed to rock/water interaction and changes in redox zonation along the flow path. In these areas, groundwater evolves along flow paths from a calcium- to a sodium-bicarbonate water type with increasing reducing conditions resulting in denitrification, elevated iron and manganese concentrations, and production of methane in the more geochemically evolved and strongest reducing conditions. In the Ozark and Springfield Plateau aquifers, rapid influx of surface-derived contaminants, especially nitrogen, coupled with few to no attenuation processes was attributed to the karst landscape developed on Mississippian- and Ordovician-age carbonate rocks of the Ozark Plateaus. Increasing nitrate concentrations are related to increasing agricultural land use, and areas of mature karst development result in higher nitrate concentrations than areas with less karst features.

Arkasas↗

Discharge and dissolved-solids characteristics and trends of Snake River above Jackson Lake at Flagg Ranch, Wyoming, 1986–2018

The headwaters of the Snake River are in the mountains of northwestern Wyoming. Maintaining the recognized high quality of water in Grand Teton National Park is a National Park Service (NPS) priority. To characterize and understand the water resources of Grand Teton National Park, the NPS established a monitoring program to monitor the quality of area surface waters. Beginning in 2006, water was sampled by the NPS and analyzed for a range of chemical species at the Snake River above Jackson Lake at Flagg Ranch streamgage 13010065 (hereafter referred to as “Snake River at Flagg Ranch”), a site where the U.S. Geological Survey (USGS) previously sampled and analyzed water from 1986 through 2004. The USGS, in cooperation with the NPS, evaluated water-quality data collected by both entities to determine if discharge and total dissolved solids (referred to as dissolved solids) have changed in the Snake River at the Flagg Ranch. To understand potential changes with time in dissolved solids, discharge was analyzed between January 1986 and December 2018, which corresponds with the time period when water-quality data were collected. Mean annual discharge varied during this time, with high, low, mean, and median flows generally increasing from 1986 through 1998, decreasing through 2005, and then generally increasing through 2018. Combining water-quality data collected by the USGS and NPS provides a longer, more complete dataset for analyses. During the period of time when NPS was the sampling agency, specific conductance data were collected, but dissolved-solids data were not. The specific conductance data from both agencies were evaluated to determine if combining the data was justified. The interquartile ranges of data collected by both agencies are similar, and rapid, large changes in values during the period of transition between USGS and NPS sampling do not occur. The USGS and NPS datasets are not statistically different in the spring, summer, or fall, but are statistically different in the winter. The winter differences could be a function of the lack of wintertime NPS sampling, which excludes higher-concentration, lower-discharge data or a function of changes in the actual concentration in the stream. Although there is some difference in the winter datasets, the similarity in sampling methods and general overall data characteristics justifies combining the data for trend analyses. Because the dissolved-solids parameter is useful for managers, it is often calculated from specific conductance using a linear regression model when dissolved-solids data are absent. For this study, creating a modeled dataset of dissolved solids for the NPS data collection period of time provided a longer, more complete dataset of dissolved-solids concentrations. The concentrations of dissolved solids over time are identified by season and indicate that samples collected in the fall and winter have higher concentrations than samples collected in spring and summer. Specifically, the mean dissolved-solids concentrations in fall and winter are around 188 milligrams per liter (mg/L), whereas the mean concentrations are around 130 mg/L in spring and summer. This difference is generally attributed to the dilution of spring and summer samples by snowmelt generated runoff during the high-flow period of the year. Trend analyses of dissolved-solids concentrations and loads indicate that an upward trend in concentration from 1986 to 2018 is likely, and a downward trend in load is highly likely. Comparing 1986 to 2018, dissolved-solids concentration is estimated to have increased by 2.25 mg/L (1.4 percent). During that same period, the dissolved-solids load is estimated to have decreased 11.8 million kilograms per year (12-percent decrease). This decrease is consistent with the estimated decrease in annual mean of daily mean discharge. Because 10 percent of the total change in dissolved-solids load is related to a change in the concentration-discharge relationship and 2 percent is related to changes in discharge, the decreased load is related less to changes in discharge and more to landscape scale processes that are affecting the concentration-discharge relationship. As noted above, the data collected by the USGS and NPS are generally comparable with regards to sampling and analytical methods, and data collected by both agencies were used as one dataset for trend analyses. The current NPS sampling schedule, however, is creating a dataset biased towards lower concentration dissolved-solids data, which occurs during higher summer flows, by only sampling during April through November. From 1986 to 2018, the percentage of NPS samples is small enough that the effect on trends is expected to be minimal. Because of the importance of low flow (winter season) data, it is likely that an April through November sampling regime may affect the ability to detect trends or determine seasonality in the future. Collection of winter data in particular is important based on the findings that the changes in the modeled concentration-discharge relationship over time have been most pronounced during the winter (represented by February) months.

Wyoming↗

Surface faults in the gulf coastal plain between Victoria and Beaumont, Texas

Displacement of the land surface by faulting is widespread in the Houston-Galveston region, an area which has undergone moderate to severe land subsidence associated with fluid withdrawal (principally water, and to a lesser extent, oil and gas). A causative link between subsidence and fluid extraction has been convincingly reported in the published literature. However, the degree to which fluid withdrawal affects fault movement in the Texas Gulf Coast, and the mechanism(s) by which this occurs are as yet unclear. Faults that offset the ground surface are not confined to the large (>6000-km 2 ) subsidence “bowl” centered on Houston, but rather are common and characteristic features of Gulf Coast geology. Current observations and conclusions concerning surface faults mapped in a 35,000-km 2 area between Victoria and Beaumont, Texas (which area includes the Houston subsidence bowl) may be summarized as follows: (1) Hundreds of faults cutting the Pleistocene and Holocene sediments exposed in the coastal plain have been mapped. Many faults lie well outside the Houston-Galveston region; of these, more than 10% are active, as shown by such features as displaced, fractured, and patched road surfaces, structural failure of buildings astride faults, and deformed railroad tracks. (2) Complex patterns of surface faults are common above salt domes. Both radial patterns (for example, in High Island, Blue Ridge, Clam Lake, and Clinton domes) and crestal grabens (for example, in the South Houston and Friendswood-Webster domes) have been recognized. Elongate grabens connecting several known and suspected salt domes, such as the fault zone connecting Mykawa, Friendswood-Webster, and Clear Lake domes, suggest fault development above rising salt ridges. (3) Surface faults associated with salt domes tend to be short (<5 km in length), numerous, curved in map view, and of diverse trend. Intersecting faults are common. In contrast, surface faults in areas unaffected by salt diapirism are frequently mappable for appreciable distances (>10 km), occur singly or in simple grabens, have gently sinuous traces, and tend to lie roughly parallel to the ENE-NE “coastwise” trend common to regional growth faults identified in subsurface Tertiary sediments. (4) Evidence to support the thesis that surface scarps are the shallow expression of faults extending downward into the Tertiary section is mostly indirect, but nonetheless reasonably convincing. Certainly the patterns of crestal grabens and radiating faults mapped on the surface above salt domes are more than happenstance; analogous fault patterns have been documented around these structures at depth. Similarly, some of the long surface faults not associated with salt domes seem to have subsurface counterparts among known regional growth faults documented through well logs and seismic data. Correlations between surface scarps and faults offsetting subsurface data are not conclusive because of the large vertical distances (1900- 3800 m) involved in making the most of the inferred connections. Nevertheless, the large number of successful correlations - in trend, movement sense, and position - suggests that many surface scarps represent merely the most recent displacements on faults formed during the Tertiary. (5) Upstream-facing fault scarps in this region of low relief can be significant impediments to streams. Locally, both abandoned, mud-filled Pleistocene distributary channels and, more commonly, Holocene drainage lines still occupied by perennial streams reflect the influence of faulting on their development. Some bend sharply near faults and have tended to flow along or pond against the base of scarps; others meander within topographically expressed grabens. Such evidence for Quaternary displacement of the ground surface is widespread in the Texas Gulf coast. In the general, however, streams in areas now offset by faulting show no disruption of their courses where they cross fault scarps. Such scarps are probably very young, and where they can be demonstrated to partly or wholly predate fluid withdrawal, very recent natural fault activity is indicated. (6) Early aerial photographs (1930) of the entire region and topographic maps (1915-16 surveys) of Harris County (Houston and vicinity) show that many faults had already displaced the land surface at a time when appreciable pressure declines in subjacent strata were localized to relatively few areas of large-scale pumping. Prehistoric faulting of the land surface, as noted above, appears to have affected much of the Texas Gulf Coast. (7) A relation between groundwater extraction and current motion on active faults is suspected because of the increased incidence of ground failure in the Houston-Galveston subsidence bowl. This argument is weakened somewhat by recognition of numerous surface faults, some of them active today, far beyond the periphery of the strongly subsiding area. Moreover, tilt beam records from two monitored faults in northwest Houston and accounts of fault damage from local residents demonstrate a complex, episodic nature of fault creep which can only partially be correlated with groundwater production. Nevertheless, although specific mechanisms are in doubt, the extraction of groundwater from shallow (<800-m) sands is probably a major factor in contributing to current displacement of the ground surface in the Houston-Galveston region. Within this large area, the number of faults recognizable from aerial photographs has increased at least tenfold between 1930 and 1970. Elsewhere in the Texas Gulf Coast only a moderate increase has been noted, some of which is possibly attributable to oil and gas production. Surface fault density in the Houston-Galveston region is far greater than in any other area of the Texas Gulf Coast investigated to date. A plausible explanation for these differences is that large overdrafts of groundwater over an extended period of time in the Houston-Galveston region have stimulated fault activity there. Throughout the Texas Gulf Coast, however, a natural contribution to fault motion remains a distinct possibility.

Texas↗

Petrology of Newberry Volcano, central Oregon

Note: This paper is dedicated to Aaron and Elizabeth Waters on the occasion of Dr. Waters' retirement. The eastern flank of the central and southern Cascade Mountains is bordered by a belt of shield volcanoes that appears to be a subprovince of the Oregon high-alumina plateau basalt petrologic province. Most of the volcanoes in this belt are low shields in which differentiation from the parent high-alumina basalt magma has been relatively slight, but several are large complex shield centers where differentiation has been extreme. The location of these large centers, and of some of the smaller volcanoes as well, was largely determined by intersecting concentrations of faults and fault-fissures of three regional fault systems. One of the largest of the complex volcanic centers is Newberry Volcano in central Oregon, a shield volcano with a big caldera at its summit. The stratigraphy of the caldera walls and of features on the caldera floor at Newberry allows detailed interpretation of the history of the younger parts of the volcano and caldera. The formation of Newberry Caldera was apparently a slow process controlled largely by faulting along the three regional fault systems. The magma conduits were probably a gridlike plexus of intersecting dikes and fissures, with larger “magma pockets” at the grid intersections. The magma was trapped in shallow chambers and periodically released by faulting. The entrapment of the magma allowed differentiation in the shallow chambers. The stratigraphy and petrology of the wall sequence also allows determination of the relative time at which the caldera had grown large enough to hold a caldera lake. On differentiation plots, chemical analyses of the Newberry rocks show two trends: rocks erupted before the presence of a lake in the caldera trend toward slight iron enrichment, whereas rocks erupted after water was present in the caldera generally trend toward alkali enrichment. These different trends are attributed to differences in the oxygen fugacity of the magma which, in turn, are related to the presence or absence of large volumes of water in the caldera lake. The interpretation is supported by field, petrographic, petrologic, chemical, trace-element, and isotopic data. Plots of existing data for the Medicine Lake Highland Volcano, another large complex shield center in the belt, show the same type of two-trend relation as those of Newberry Volcano.

Oregon↗

Loess studies in central United States: Evolution of concepts

Few words in the realm of earth science have caused more debate than "loess". It is a common term that was first used as a name of a silt deposit before it was defined in a scientific sense. Because this "loose" deposit is easily distinguished from other more coherent deposits, it was recognized as a matter of practical concern and later became the object of much scientific scrutiny. Loess was first recognized along the Rhine Valley in Germany in the 1830s and was first noted in the United States in 1846 along the lower Mississippi River where it later became the center of attention. The use of the name eventually spread around the world, but its use has not been consistently applied. Over the years some interpretations and stratigraphic correlations have been validated, but others have been hotly contested on conceptual grounds and semantic issues. The concept of loess evolved into a complex issue as loess and loess-like deposits were discovered in different parts of the US. The evolution of concepts in the central US developed in four indefinite stages: the eras of (1) discovery and development of hypotheses, (2) conditional acceptance of the eolian origin of loess, (3) "bandwagon" popularity of loess research, and (4) analytical inquiry on the nature of loess. Toward the end of the first era around 1900, the popular opinion on the meaning of the term loess shifted from a lithological sense of loose silt to a lithogenetic sense of eolian silt. However, the dual use of the term fostered a lingering skepticism during the second era that ended in 1944 with an explosion of interest that lasted for more than a decade. In 1944, R.J. Russell proposed and H.N. Fisk defended a new non-eolian, property-based, concept of loess. The eolian advocates reacted with surprise and enthusiasm. Each side used constrained arguments to show their view of the problem, but did not examine the fundamental problem, which was not in the proofs of their hypothesis, but in the definition of the term. Between 1944 and about 1950, the debates about loess reached a maximum level of complexity. The main semantic problem was submersed in peripheral arguments about physical properties and genetic interpretations. The scholarly treatment of the subject by Fisk and Russell stimulated quality responses from a diversity of earth scientists interested in academic and applied studies, particularly geo-history, pedology, soil mechanics and stratigraphy. The long-lasting popularity of loess studies during the bandwagon era lasted to about 1970. By that time, the analytical and technical interests had attracted the mainstream into the fourth era with a focus beyond the old arguments. Although Fisk and Russell found themselves defending an unpopular theory, they stimulated a scientific interest in the late Quaternary history of the Mississippi Valley that may never be exceeded.

Engineering Geology↗

Miocene tholeiitic basalts of coastal Oregon and Washington and their relations to coeval basalts of the Columbia Plateau

Note: This paper is dedicated to Aaron and Elizabeth Waters on the occasion of Dr. Waters' retirement. Tholeiitic basalt flows and breccias of Miocene age in western Oregon and Washington form three distinct stratigraphic units. Each unit was erupted from coastal vents marked by dikes and sills of the same composition as associated extrusive rocks. The three coastal basalt units are interbedded with predominantly marine sedimentary rocks of middle to late Miocene age. These units are here named, from older to younger, the Depoe Bay Basalt, Cape Foulweather Basalt, and basalt of Pack Sack Lookout. The three units can be distinguished by their petrographic characteristics. The Depoe Bay Basalt is nonporphyritic; Cape Foulweather Basalt has sparse large labradorite phenocrysts; and Pack Sack basalt has labradorite phenocrysts with numerous pyroxene and glass inclusions as well as augite and olivine phenocrysts. Chemical analyses of basalts from these three units show that each has a distinct and uniform composition. The Depoe Bay Basalt is characterized by high SiO 2 content; the Cape Foulweather Basalt has high content of total iron, TiO 2 , and P 2 O 5 ; and Pack Sack basalt is marked by relatively high MgO and CaO content. The Depoe Bay Basalt, Cape Foulweather Basalt, and basalt of Pack Sack Lookout on the coast occur in the same stratigraphic order and are essentially the same ages as three basalt units that erupted on the Columbia Plateau. The plateau-derived units are the Yakima and late-Yakima petrographic types of Waters (1961) and the Pomona flow of Schmincke (1967). The virtual identity in chemical composition of the Depoe Bay Basalt and Yakima-type basalt, the Cape Foulweather Basalt and the late-Yakima–type basalt, and the Pack Sack basalt and the Pomona basalt flow indicate that each pair is consanguineous. Fissure vents for the plateau basalt are located in eastern Oregon and Washington and western Idaho more than 500 km east of the coastal vent areas. Thus, a regional mechanism of magma generation or emplacement is required. Three models of magma genesis considered in this report are: (1) partial melting of the subducted Juan de Fuca plate; (2) partial melting along a nearly horizontal shear zone at the base of the American plate; and (3) partial melting within the asthenosphere and fractionation during ascent of the magma.

Oregon, Washington↗

Geographic information system (GIS)-based maps of Appalachian basin oil and gas fields

One of the more recent maps of Appalachian basin oil and gas fields (and the adjoining Black Warrior basin) is the U.S. Geological Survey (USGS) compilation by Mast and others (1998) (see Trippi and others, this volume, chap. I.1). This map is part of a larger oil and gas field map for the conterminous United States that was derived by Mast and others (1998) from the Well History Control System (WHCS) database of Petroleum Information, Inc. (now IHS Energy Group). Rather than constructing the map from the approximately 500,000 proprietary wells in the Appalachian and Black Warrior part of the WHCS database, Mast and others (1998) subdivided the region into a grid of 1-mi 2 (square mile) cells and allocated an appropriate type of hydrocarbon production (oil production, gas production, oil and gas production, or explored but no production) to each cell. Each 1-mi 2 cell contains from 0 to 5 or more exploratory and (or) development wells. For example, if the wells in the 1-mi 2 cell consisted of three oil wells, one gas well, and one dry well, then the cell would be characterized on the map as an area of oil and gas production. The map by Mast and others (1998) accurately shows the distribution and types of hydrocarbon accumulation in the Appalachian and Black Warrior basins, but it does not show the names of individual fields. To determine the locality and name of individual oil and gas fields, one must refer to State oil and gas maps (for example, Harper and others, 1982), which are generally published at scales of 1:250,000 or 1:500,000 (see References Cited), and (or) published journal articles. Other recent USGS Appalachian basin oil and gas field maps show the distribution of oil and gas production with a cell size as small as 0.25 mi 2 , such as the maps converted by Trippi and others (this volume, chap. I.1) from proprietary well-location maps used in the USGS 2002 assessment of oil and gas resources of the Appalachian basin (Milici and others, 2003). Another set of Appalachian basin oil and gas cell maps (based on a cell size of 0.25 mi 2 ) was created for the USGS 1995 National Assessment of United States Oil and Gas Resources (Gautier and others, 1995; Beeman and others, 1996). Between 1991 and 1994, R.T. Ryder (with R.E. Mattick, J.B. Roen, and J.R. San Filipo, USGS, Reston, Va.) compiled oil and gas fields on stable-base mylar greenline base maps (scale 1:500,000) for selected plays in the Appalachian basin. These map compilations included field names and field numbers where assigned by State agencies. The purpose of the maps was to provide supporting data for the USGS 1995 National Assessment of United States Oil and Gas Resources (Gautier and others, 1995). In particular, the greenline oil and gas field maps were linked, where possible, with production data from State records and (or) published literature in order to determine ultimate sizes for conventional fields and estimated ultimate recovery (EUR) values for wells in continuous accumulations (for definitions of the conventional and continuous terminology, see USGS National Oil and Gas Assessment Team, 1995; Schmoker, 1997; Schenk and Pollastro, 2002). This approach was used in the 1995 national oil and gas assessment because ultimate field size and EUR data were unavailable in the Appalachian region from Petroleum Information, Inc., and other commercial sources. In 2006 and 2007, the greenline Appalachian basin field maps were digitized under the supervision of Scott Kinney and converted to geographic information system (GIS) files for chapter I.1 (this volume). By converting these oil and gas field maps to a digital format and maintaining the field names where noted, they are now available for a variety of oil and gas and possibly carbon-dioxide sequestration projects. Having historical names assigned to known digitized conventional fields provides a convenient classification scheme into which cumulative production and ultimate field-size databases can be organized. Moreover, as exploratory and development drilling expands across the basin, many previously named fields that were originally treated as conventional fields have evolved into large, commonly unnamed continuous-type accumulations. These new digital maps will facilitate a comparison between EUR values from recently drilled, unnamed parts of continuous accumulations and EUR values from named fields discovered early during the exploration cycle of continuous accumulations.

Appalachian basin↗

Simulation of the interaction of karstic lakes Magnolia and Brooklyn with the upper Floridan Aquifer, southwestern Clay County, Florida

The stage of Lake Brooklyn, in southwestern Clay County, Florida, has varied over a range of 27 feet since measurements by the U.S. Geological Survey began in July 1957. The large stage changes have been attributed to the relation between highly transient surface-water inflow to the lake and subsurface conduits of karstic origin that permit a high rate of leakage from the lake to the Upper Floridan aquifer. After the most recent and severe stage decline (1990-1994), the U.S. Geological Survey began a study that entailed the use of numerical ground-water flow models to simulate the interaction of the lake with the Upper Floridan aquifer and the large fluctuations of stage that were a part of that process. A package (set of computer programs) designed to represent lake/aquifer interaction in the U.S. Geological Survey Modular Finite-Difference Ground-Water Flow Model (MODFLOW-96) and the Three-Dimensional Method-of-Characteristics Solute-Transport Model (MOC3D) simulators was prepared as part of this study, and a demonstration of its capability was a primary objective of the study. (Although the official names are Brooklyn Lake and Magnolia Lake (Florida Geographic Names), in this report the local names, Lake Brooklyn and Lake Magnolia, are used.) In the simulator of lake/aquifer interaction used in this investigation, the stage of each lake in a simulation is updated in successive time steps by a budget process that takes into account ground-water seepage, precipitation upon and evaporation from the lake surface, stream inflows and outflows, overland runoff inflows, and augmentation or depletion by artificial means. The simulator was given the capability to simulate both the division of a lake into separate pools as lake stage falls and the coalescence of several pools into a single lake as the stage rises. This representational capability was required to simulate Lake Brooklyn, which can divide into as many as 10 separate pools at sufficiently low stage. In the first of two calibrated models, recharge to the water table, specified as a monthly rate, was set equal to 40 percent of the monthly rainfall rate. The specified rate of inflow to the uppermost stream segment was set equal to outflows from Lake Lowry estimated from lake stage and the 1994-97 rating table. Leakage to the intermediate and Upper Floridan aquifers was assumed to occur from the surficial aquifer system through the confining layers directly beneath deeper parts of the lake bottom. A leakance coefficient value of 0.001 feet per day per foot of thickness was used beneath Lake Magnolia, and a value of 0.005 feet per day per foot of thickness was used beneath most of Lake Brooklyn. With these values, the conductance through the confining layers beneath Lake Brooklyn was about 19 times that beneath Lake Magnolia. The simulated stages of Lake Brooklyn matched the measured stages reasonably well in the early (1957-72) and later (1990-98) parts of the simulation time period, but the match was unsatisfactory in an intermediate time period (1973-89). To resolve this discrepancy, the hypothesis was proposed that undocumented losses of water from Alligator Creek upstream from Lake Brooklyn or from the lake itself occurred between 1973 and 1989 when there was sufficient streamflow. The resulting simulation of lake stages matched the measured lake stages accurately during the entire simulation time period. The model was then revised to incorporate the assumption that only 20 percent of precipitation recharged the water table (the second calibrated model). Recalibration of the model required that leakance values for the confining units under deeper parts of the lakes also be reduced by nearly 50 percent. The stages simulated with the new parameter assumptions, but retaining the assumption of surface-water losses, were an excellent match of the measured values. The stage of Lake Magnolia was also simulated accurately. The results of sensitivity analyses show that simulated streamflow between Lakes Magnolia and Brooklyn tends to be water-budget controlled, and is not appreciably affected by the specified outflow altitude or channel characteristics of the receiving stream. To match heads measured in observation wells of the surficial aquifer network, the assigned hydraulic conductivity values were zoned, and ranged from a minimum of 4 feet per day to a maximum of 400 feet per day in the first calibrated model. These values were reduced by about 50 percent in the second calibrated model. Differences between observation wells were noted in the abruptness of changes of measured head values, and in the relation of the timing of peak measured heads and simulated peak heads. These differences seemed to be correlated with the depth of the water table below land surface. Spatially uniform values of transmissivity were specified for the intermediate (10,000 feet squared per day) and Upper Floridan (100,000 feet squared per day) aquifers. Simulated heads in the Upper Floridan aquifer layer follow the trend of the heads measured in a long-term observation well with data beginning in 1960. This result suggests that the observed head decline could be explained entirely in terms of the stage decline in Lake Brooklyn and may not indicate a regional trend.

Florida↗

Zoning of ore deposits in and adjoining the Leadville District, Colorado

The ore deposits in the heart of the Leadville district , Colorado , are divisible, as shown by Emmons, Irving, and Loughlin,2 into two main classes on the basis of mineral composition: "contact-metamorphic" silicates accompanied by iron oxides, which adjoin an obscure igneous stock, and mixed sulphides accompanied by manganosiderite and quartz gangue, which cut the "contact-metamorphic" deposits and form a broad zone extending for 2 1/2 miles to the west and for shorter distances in other directions. A detailed consideration of the mineralogy of the mixed sulphide ore deposits clearly shows that they in turn are divisible into groups, referred to as hotter and intermediate mesothermal deposits . These groups, however, because of complications in structural control, are not ideally arranged from the standpoint of zoning .Emmons, Irving, and Loughlin noted that barite was a conspicuous gangue mineral in certain outlying deposits on the west slope of the Mosquito Range. The junior author during 1928-1931 had an opportunity to study the outlying deposits over an area of twenty square miles, which includes mines near the heads of Iowa Gulch and Empire Gulch, and the Weston Pass district (Fig. 1), and found that these baritic deposits form a somewhat distinct outer zone in which the primary ore consists of galena and light-colored zinc blende, with small amounts of silver and practically no gold, in a gangue of barite, jasperoid, and crystalline quartz. There is also recognized a still more remote though poorly defined zone in which the same sulphides are present but jasperoid is less conspicuous and barite is practically absent, whereas dolomite is the characteristic gangue mineral. The deposits in these two outer zones are fissure-fillings and comparatively small replacement bodies. They contrast strongly with the sulphide deposits in the heart of the Leadville district and are assigned to the outer mesothermal and epithermal zones. As they have not been described before, general features of these deposits are given in some detail, followed by their relations to the deposits within the Leadville district . Analogies are also drawn with similar deposits in other districts. In the preparation of this paper the authors have been impressed with the confusion that may result from too simple a conception of zoning . For a proper conception of Leadville , the geologist must have sufficient regard for structural conditions, must make due allowance for successive stages of mineral deposition within a small block of ground or even within a single ore shoot, and must also allow for the reactions between minerals formed at an early stage and solutions introduced during a late stage. The effects of such complications in the Leadville district might give the impression that certain minerals found in one of the inner zones were formed at relatively high temperatures, but close study of paragenesis indicates late deposition at relatively low temperatures. Again, as already stated, there is sufficient inconsistency in the distribution of certain varieties of deposits to give a confused impression of zonal arrangement unless the influence of structure on the circulation of ore -forming solutions is appreciated. The field facts give a far more complex picture than might be anticipated from a simple statement of the theory of zoning .

Colorado↗

Geochronology of Precambrian rocks of the Teton Range, Wyoming

Note: This paper is dedicated to Aaron and Elizabeth Waters on the occasion of Dr. Waters' retirement. The oldest rocks in the Teton Range are complexly deformed interlayered biotite gneiss, plagioclase gneiss, amphibole gneiss, and amphibolite. Also, within these rocks, there are concordant bodies of strongly lineated quartz monzonite gneiss, here named the Webb Canyon Gneiss, which may be of volcanic origin. Coarse metagabbro, here named the Rendezvous Metagabbro, is intrusive into the layered gneiss sequence and was metamorphosed and deformed along with the enclosing rocks. These older rocks are cut by discordant plutons and swarms of undeformed dikes of quartz monzonite and associated pegmatite. The quartz monzonite, which makes up much of the central part of the Teton Range, is here named the Mount Owen Quartz Monzonite. The youngest Precambrian rocks are undeformed dikes of slightly metamorphosed tholeiitic diabase. A Rb-Sr whole-rock isochron on the Webb Canyon Gneiss and the Rendezvous Metagabbro indicates that these rocks were metamorphosed 2,875 ± 150 m.y. ago. The initial Sr. ratio of 0.700 suggests that the original rocks are probably not significantly older than the metamorphism. The Mount Owen Quartz Monzonite has a whole-rock isochron age of 2,495 ± 75 m.y. and an unusually high initial ratio of 0.732. Plagioclase-microcline isochrons from two samples of the quartz monzonite indicate partial re-equilibration of the Rb-Sr system during a thermal event 1,800 m.y. ago. The age of the diabase dikes has not been definitely determined, but biotite in the wall rocks of one major dike has a K-Ar age of 1,450 m.y. This suggests that the dike was emplaced during or prior to a thermal event 1,300 to 1,500 m.y. ago that was responsible for resetting many of the previously reported K-Ar mineral ages throughout the range. The geochronologic record in the Teton Range is very similar to that elsewhere in the Wyoming Precambrian province. Major metamorphic events with ages between 2,700 and 2,900 m.y. have been identified in the Bighorn, Beartooth, Little Belt, and Granite Mountains. Post-tectonic granitic rocks with ages of 2,500 to 2,700 m.y. have been found in the Wind River Range and the Granite Mountains. Later thermal events have affected Rb-Sr systematics of rocks in the Beartooth Mountains, Wind River Range, and Granite Mountains, as well as in the Teton Range at about the same time as major episodes of regional metamorphism in terranes flanking the Wyoming province in southwestern Montana and in the Front Range in Colorado.

Wyoming↗

Notes on the origin of copromacrinite based on nitrogen functionalities and δ13C and δ15N determined on samples from the Peach Orchard coal bed, southern Magoffin County, Kentucky

This paper represents the first attempt to show, by means other than just petrographic ones, that one type of macrinite, herein designated copromacrinite, may result from macrofauna feces. For that purpose a combination of coal petrography, X-ray photoelectron spectroscopy, and elemental-analysis continuous-flow isotope ratio mass spectrometry methods were used to determine nitrogen functionalities and &delta; 13 C and &delta; 15 N compositions in 1) vitrinite-rich, 2) fusinite + semifusinite-rich, and 3) macrinite-rich (with a possible coprolitic origin) samples of the high volatile A bituminous Peach Orchard coal (Bolsovian; Middle Pennsylvanian) from Magoffin County, Kentucky. There were no significant differences between pyridinic-N and quaternary-N abundance in the three samples, however, pyrrolic-N was higher (~ 54%) in the macrinite-rich sample than in the other two samples (~ 38%). The data suggest that pyridinic-N and quaternary-N are independent of maceral group composition and that pyrrolic-N is dependent on maceral composition (fusinite + semifusinite versus macrinite). &delta; 13 C values obtained for bulk and demineralized coal of the vitrinite- and fusinite + semifusinite-rich samples are similar with &delta; 13 C values of &minus; 24.80 &plusmn; 0.01&permil; VPDB and &minus; 24.61 &plusmn; 0.09&permil; VPDB for bulk samples and &minus; 24.81 &plusmn; 0.07&permil; VPDB and &minus; 24.52 &plusmn; 0.04&permil; VPDB for demineralized samples. These values are within the expected range for vitrinite-rich samples and the slightly higher &delta; 13 C value of the fusinite + semifusinite-rich sample is expected as &delta; 13 C values for inertinite are higher than for vitrinite. However, there was a significant shift to a lower &delta; 13 C value (&minus; 26.80 &plusmn; 0.01&permil; VPDB for the bulk sample value) for the macrinite-rich sample. Because the samples are basically isorank, and &delta; 13 C (and &delta; 15 N) shifts do not occur during maturation until anthracite rank, the difference may be related to the presence or composition of the macrinite within the sample which lacks heat-effect indicators, such as devolatilization vacuoles and distorted pores. &delta; 15 N values are also similar for bulk and demineralized coal of the vitrinite- and fusinite + semifusinite-rich samples, and the bulk values were heavier in this samples (3.07 &plusmn; 0.03&permil; Air and 2.92 &plusmn; 0.10&permil; Air, respectively), and much lighter (&minus; 2.83 &plusmn; 0.09&permil; Air) for the macrinite-rich sample. The study of Peach Orchard coal samples using reflected-light microscopy, isotopic composition, and nitrogen-forms analyses revealed that the macrinite-rich sample contains macrinite with coprolitic features (e.g. oxidation rind, mix of undigested palynomorphs, frequent and randomly located funginite, agglutination pulp of semifusinite reflectance, internal lack of bedding fabric, and suggestion of structures resulting from intestines and stomach walls), more pyrrolic-N (~ 16%), and lower &delta; 13 C (~ 2&permil; VPDB) and &delta; 15 N (~ 4&permil; Air) values than the vitrinite and semifusinite + fusinite rich samples. These findings suggest that the maceral macrinite has multiple origins based on petrography and measurable chemical differences between the macrinite, vitrinite, and semifusinite + fusinite fractions within the coal. Assuming that copromacrinite observed is an excretion then the anomalies observed may result from the symbiotic relations between the macrofauna (e.g. cockroaches) and microbiota during the digestive processes, and the nitrogen balance mechanisms inside macrofauna body.

Kentucky↗

Loess stratigraphy of the Lower Mississippi Valley

Loesses of the Lower Mississippi Valley (LMV) are world-famous. Sir Charles Lyell (1847), Hilgard (1860), Stafford (1869), Call (1891) and Mabry (1898), thought the LMV loess was a single water deposit although "double submergence" was noted by Call (1891) and Salisbury (1891). Shimek (1902) and Emerson (1918) recognized LMV loess as a wind deposit which came from the valley. Although wind-deposited loess gained wide acceptance, Russell (1944a) published his controversial theory of "loessification" which entailed weathering of backswamp deposits, downslope movement and recharge by carbonates to form loess. Wascher et al. (1947) identified three LMV loesses, mapped distributions and strongly supported eolian deposition. Leighton and Willman (1950), identified four loesses and supported eolian deposition as did Krinitzsky and Turnbull (1967) and Snowden and Priddy (1968), but Krinitzsky and Turnbull questioned the deepest loess. Daniels and Young (1968) and Touchet and Daniels (1970) studied the distribution of loesses in south-central Louisiana. West et al. (1980) and Rutledge et al. (1985) studied the source areas and wind directions which deposited the loesses on and adjoining Crowley's Ridge. B.J. Miller and co-workers (Miller et al., 1985, 1986, Miller and Alford, 1985) proposed that the Loveland Silt was Early Wisconsin rather than Illinoian age and advanced the name Sicily Island loess. They proposed the underlying loess was Illinoian and advanced the name Crowley's Ridge. We termed the loesses, from the surface downward, Peoria Loess, Roxana Silt, Loveland/Sicily Island loess, Crowley's Ridge Loess and Marianna loess. Researchers agree that the surfical Peoria Loess is Late Wisconsin and the Roxana Silt is Late to Middle Wisconsin, but little agreement exists on the age of the older loesses. Pye and Johnson (1988) proposed Early Wisconsin for the Loveland/Sicily Island. McKay and Follmer (1985) suggested this loess correlated with a loess under Illinoian till. Clark et al. (1989) agreed on Crowley's Ridge, but suggested the Loveland/Sicily Island loess on Sicily Island was older. Mirecki and Miller (1994) and Millard and Maat (1994) suggested an Illinoian age for the Loveland/Sicily Island loess. Miller and co-workers suggested, as did Pye and Johnson (1988), an Illinoian age for the Crowley's Ridge loess. McKay and Follmer (1985) suggested it correlated with a loess under "Kansan" till. Stratigraphy indicates the Marianna is the older of the five loesses. Researchers identified loess on both the east and west side of the LMV as well as on higher terraces within the valley. Many researchers assumed unaltered loesses were commonly yellowish brown, and silts or silt loams (West et al., 1980; Miller et al., 1986). The nonclay fraction of unweathered LMV loesses was dominated by quartz followed by carbonates, mainly dolomites, followed by feldspars, and micas. Clays were dominated by montmorillonite followed by micaceous minerals, kaolinite and vermiculite (Miller et al., 1986). Soils in the Crowley's Ridge loess are most developed, followed by the soils in the Loveland/Sicily Island which are more developed than the modern soils in the Peoria Loess. Soils in the Roxana and Marianna loesses are least developed and the Farmdale Soil of the Roxana is the weaker of the two (Miller et al., 1986). There is certainly overlapping range in the degree of soil development in the various loesses.

Engineering Geology↗

Factors affecting water quality in selected carbonate aquifers in the United States, 1993-2005

Carbonate aquifers are an important source of water in the United States; however, these aquifers can be particularly susceptible to contamination from the land surface. The U.S. Geological Survey National Water-Quality Assessment (NAWQA) Program collected samples from wells and springs in 12 carbonate aquifers across the country during 1993–2005; water-quality results for 1,042 samples were available to assess the factors affecting ground-water quality. These aquifers represent a wide range of climate, land-use types, degrees of confinement, and other characteristics that were compared and evaluated to assess the effect of those factors on water quality. Differences and similarities among the aquifers were also identified. Samples were analyzed for major ions, radon, nutrients, 47 pesticides, and 54 volatile organic compounds (VOCs). Geochemical analysis helped to identify dominant processes that may contribute to the differences in aquifer susceptibility to anthropogenic contamination. Differences in concentrations of dissolved oxygen and dissolved organic carbon and in ground-water age were directly related to the occurrence of anthropogenic contaminants. Other geochemical indicators, such as mineral saturation indexes and calcium-magnesium molar ratio, were used to infer residence time, an indirect indicator of potential for anthropogenic contamination. Radon exceeded the U.S. Environmental Protection Agency proposed Maximum Contaminant Level (MCL) of 300 picocuries per liter in 423 of 735 wells sampled, of which 309 were drinking-water wells. In general, land use, oxidation-reduction (redox) status, and degree of aquifer confinement were the most important factors affecting the occurrence of anthropogenic contaminants. Although none of these factors individually accounts for all the variation in water quality among the aquifers, a combination of these characteristics accounts for the majority of the variation. Unconfined carbonate aquifers that had high percentages of urban or agricultural land, or a combination of both, had higher concentrations and higher frequency of detections for most of the anthropogenic contaminants than areas with other combinations of land use and degree of aquifer confinement. Redox status is an indicator of more recently recharged water and affects the fate of some contaminants. Median concentrations of nitrate were highest in the Valley and Ridge and Piedmont aquifers and lowest in the Biscayne and Silurian-Devonian/Upper carbonate aquifers. Nitrate concentrations were significantly higher in unconfined aquifers than in confined aquifers and semiconfined/mixed confined aquifers (wells in aquifers with breached confining layers or wells open to both a confined and an unconfined aquifer). Water recharged after 1953 had significantly higher concentrations of nitrate than water recharged prior to 1953. Redox status was also a key factor affecting nitrate concentrations; in recently recharged waters, samples in oxic waters had significantly higher concentrations of nitrate than anoxic waters, regardless of land use in the area around the well. Samples from 54 wells (5 percent) exceeded the U.S. Environmental Protection Agency MCL of 10 mg/L for nitrate in drinking water. Most of the samples exceeding the drinking-water standard (52 samples, or 5 percent) were in domestic supply wells in agricultural areas. The Piedmont and Valley and Ridge aquifers had the largest number of samples (45) exceeding the MCL; in the remaining aquifers only 9 samples had concentrations of nitrate that exceeded the MCL (about 1 percent). None of the water recharged prior to 1953 and only a single sample from a confined aquifer had nitrate concentrations that exceeded 10 mg/L as N. Wells were sampled for a minimum of 47 pesticides. Detection frequencies and comparisons varied depending on the assessment level used. At least 1 of the 47 pesticides was detected at 510 (50 percent) of the 1,027 sites where pesticide data were available using the ‘all detections’ assessment level—that is, including any quantified detection as well as any estimated values where the compound was definitively detected. Multiple pesticides were frequently detected in a sample of water from a site; 34 percent of the samples had two to five pesticides detected in the same sample, and 4 percent of the samples had six or more pesticides detected. Dieldrin was detected at 20 sites, 9 of which were from either domestic or public supply wells, at a concentration above the Health-Based Screening Level (HBSL) of 0.002 µg/L. Diazinon was detected at a concentration greater than the HBSL of 1 µg/L at a single site, which was also a domestic supply well. These are the only samples where a pesticide exceeded a human-health benchmark. The most frequently occurring pesticide compounds were four herbicides—atrazine, simazine, metolachlor, and prometon—and deethylatrazine, a degradate of atrazine. These pesticides typically were detected at concentrations that were less than 10 percent of a human-health benchmark. Of the four frequently occurring pesticides, only samples for atrazine (3 percent) and simazine (0.1 percent) had concentrations that exceeded 10 percent of the human-health benchmark; most of these cases were in agricultural areas. It is important to note, however, that the most frequently occurring pesticide degradate compound—deethylatrazine—has no human-health benchmark. Using a common assessment level of 0.01 µg/L, four of the aquifers—Biscayne, Mississippian, Piedmont, and Valley and Ridge—had at least one of these five compounds detected in more than 30 percent of the wells sampled. These four aquifers, along with the Ordovician, Ozark Plateaus, and Prairie du Chien aquifers were the aquifers or aquifer systems that had concentrations of pesticides that exceeded 10 percent of a human-health benchmark. Water recharged after 1953 had a significantly higher percentage of detections of pesticides than water recharged before 1953, and water from unconfined aquifers had a significantly higher percentage of detections of pesticides than water from confined or semiconfined/mixed confined aquifers. Water from sites in unconfined aquifers, where land use was agricultural or urban, accounted for the vast majority of detections of pesticides. Dissolved oxygen concentration was positively related to pesticide occurrence, which likely reflects the positive association between dissolved oxygen concentration and recently recharged water. Water samples were collected for analysis of VOCs at 793 sites—154 samples were analyzed for 54 VOCs from 1993 through 1995 and 639 samples were analyzed for 86 VOCs from 1996 through 2005. Twenty percent of samples contained one or more VOCs at concentrations greater than or equal to 0.2 µg/L (159 of 793 samples). The aquifers with the highest percentage of samples containing one or more VOCs were the Castle Hayne (about 41 percent of samples) and Biscayne aquifers (34 percent). The most frequently detected VOCs were chloroform, tetrahydrofuran, tetrachloroethene (PCE), toluene, acetone, ethylmethylketone, methyl tert-butyl ether (MTBE), and trichloroethene (TCE). Low-level concentrations of VOCs occurred in a much larger percentage of a subset of the data (the 639 samples analyzed using a low-level analytical method). In these samples, 69 percent of the 639 samples contained 1 or more VOCs, indicating the vulnerability of the carbonate aquifers to low-level VOC contamination. Four VOCs were detected at concentrations exceeding their respective MCLs in five samples, all of which were from drinking-water wells. Vinyl chloride concentrations exceeded the MCL of 2 µg/L in two samples from urban areas in the unconfined Biscayne aquifer. PCE, TCE, and 1,2-dichloropropane each had one sample with a concentration greater than their MCLs of 5 µg/L; these samples were from agricultural and urban areas in the unconfined Mississippian aquifer. Water quality in the 12 carbonate aquifers was highly variable. Most of the samples met drinking-water standards. The occurrence of anthropogenic contaminants was related to contaminant sources but also was affected by degree of aquifer confinement, ground-water age, and redox status. Areas with higher amounts of agricultural or urban land in unconfined aquifers were the most likely to have elevated concentrations of anthropogenic contaminants.

Scientific Investigations Report↗

Groundwater-flow model and effects of projected groundwater use in the Ozark Plateaus Aquifer System in the vicinity of Greene County, Missouri — 1907-2030

Recent and historical periods of rapid growth have increased the stress on the groundwater resources in the Ozark aquifer in the Greene County, Missouri area. Historical pumpage from the Ozark aquifer has caused a cone of depression beneath Springfield, Missouri. In an effort to ease its dependence on groundwater for supply, the city of Springfield built a pipeline in 1996 to bring water from Stockton Lake to the city. Rapid population growth in the area coupled with the expanding cone of depression raised concern about the sustainability of groundwater as a resource for future use. A groundwater-flow model was developed by the U.S. Geological Survey in cooperation with Greene County, Missouri, the U. S. Army Corps of Engineers, and the Missouri Department of Natural Resources to assess the effect that increased groundwater demand is having on the long-term availability of groundwater in and around Greene County, Missouri. Three hydrogeologic units were represented in the groundwater-flow model: the Springfield Plateau aquifer, the Ozark confining unit, and the Ozark aquifer. The Springfield Plateau aquifer is less than 350 feet thick in the model area and generally is a low yield aquifer suitable only for domestic use. The Ozark aquifer is composed of a more than 900-foot thick sequence of dolomite and sandstone in the model area and is the primary aquifer throughout most of southern Missouri. Wells open to the entire thickness of the Ozark aquifer typically yield 1,000 gallons per minute or more. Between the two aquifers is the Ozark confining unit composed of as much as 98 feet of shale and limestone. Karst features such as sinkholes, springs, caves, and losing streams are present in both aquifers, but the majority of these features occur in the Springfield Plateau aquifer. The solution-enlarged fracture and bedding plane conduits in the karst system, particularly in the Springfield Plateau aquifer, are capable of moving large quantities of groundwater through the aquifer in relatively short periods of time. Pumpage rates in the model area increased from 1,093,268 cubic feet per day in 1962 to 2,693,423 cubic feet per day in 1987 to 4,330,177 cubic feet per day in 2006. Annual precipitation ranged from 25.21 inches in 1953 to 62.45 inches in 1927 from 1915 to 2006 in the model area. Recharge to the model was calculated as 2.53 percent of the annual precipitation and was varied annually. Recharge was distributed over the model area based on land slope and was adjusted in the city limits of Springfield to account for the impervious surface. A groundwater model with annual stress periods from 1907 to 2030 was developed using a transient calibration period from 1987 to 2006 and a prediction period from 2007 to 2030 to simulate flow in the Springfield Plateau aquifer and the Ozark aquifer. For the model area of approximately 2,870 square miles, the model hydrogeologic units and hydraulic properties were discretized into 253 rows, 316 columns, and 3 layers with the layer boundaries crossing hydrogeologic unit boundaries in some areas. The horizontal cell spacing was 1,000 feet by 1,000 feet. The model was calibrated by minimizing the difference between simulated head and observed water levels and simulated and observed flows in rivers and springs. Population and the associated groundwater use were estimated for 12 communities and the unincorporated area of Greene County based on past growth. Each was analyzed individually, and a low and high annual rate of growth relative to the 2006 population was computed for each community or group. Low growth rates ranged from 0.215 percent per year in Springfield to 6.997 percent per year in Rogersville. Total growth from 2006 to 2030 at the low growth rate ranged from 5.2 percent in Springfield to 167.9 percent in Rogersville. High growth rates ranged from 0.236 percent per year in Springfield to 7.345 percent per year in Rogersville. Total growth from 2006 to 2030 at the high growth rate ranged from 5.7 percent in Springfield to 176.3 percent in Rogersville. Response of the flow system to selected hypothetical pumping stresses and recharge conditions was simulated using the calibrated model. Seven hypothetical scenarios were simulated from 2007 to 2030 to test the effects of various stresses on the head in the Ozark aquifer. Hypothetical scenario 1 continued the 2006 pumping rates without change to the end of 2030. Scenario 2 assumed a low population growth rate with a 4-year drought at the beginning of the prediction period. Scenario 3 assumed a low population growth rate with a 4-year drought at the end of the prediction period. Scenario 4 assumed a high population growth rate with a 4-year drought at the beginning of the prediction period. Scenario 5 assumed a high population growth rate with a 4-year drought at the end of the prediction period. Scenario 6 and 7 had one new industrial well installed within the city limits of Springfield and one new industrial well installed about 3.5 miles east of Rogersville. Scenario 6 assumed a low population growth rate and scenario 7 assumed a high population growth rate. Results were compared by examining differences in head at the end of the simulation period. All scenarios examined resulted in potentiometric-surface declines from 2006 levels. Results from scenario 1 indicated that even with no increase in pumping, the potentiometric surface in the Springfield area continued to decline. The maximum decline of approximately 62 feet from the 2006 potentiometric surface occurred in Springfield. The maximum decline from the 2006 potentiometric surface in scenarios 2 and 3 was approximately 203 feet and in scenarios 4 and 5 was approximately 207 feet. The drought occurring at the end of the simulation period tended to broaden the drawdown area relative to the drought at the beginning. Drought timing did not substantially affect the potentiometric surface in the Ozark aquifer except for where the Ozark aquifer was exposed. Although not a substantial difference, the high population growth rate scenarios tended to have larger declines than the low population growth rate scenarios. As in the previous scenarios, little difference was noted between the low and high growth rate in scenario 6 and 7. Scenarios 6 and 7 showed declines of more than 640 feet from the 2006 potentiometric surface at the new well located in Springfield. The drawdown at the new wells decreased relatively quickly with increased distance from the well. Simulated head in the nearby cities of Nixa, Ozark, and Republic was nearly the same for scenarios 2 through 7 and was lower than the head predicted for scenario 1. Results from scenarios 2 through 7 indicate that the potentiometric surface in 2030 near these cities could decline 100 feet or more from the 2006 levels. Because model layers 2 and 3, representing the Ozark confining unit and most of the thickness of the Ozark aquifer, were simulated as confined, drawdown in the wells in the area of the Ozark aquifer that is unconfined or becomes unconfined during the simulation period will likely be under predicted.

Missouri↗

U.S. Geological Survey science for the Wyoming Landscape Conservation Initiative—2014 annual report

This is the seventh report produced by the U.S. Geological Survey (USGS) for the Wyoming Landscape Conservation Initiative (WLCI) to detail annual activities conducted by the USGS for addressing specific management needs identified by WLCI partners. In FY2014, there were 26 projects, including a new one that was completed, two others that were also completed, and several that entered new phases or directions. The 26 projects fall into several categories: (1) synthesizing and analyzing existing data to identify current conditions on the landscape and using the data to develop models for projecting past and future landscape conditions; (2) monitoring indicators of ecosystem conditions and the effectiveness of on-the-ground habitat projects; (3) conducting research to elucidate the mechanisms underlying wildlife and habitat responses to changing land uses; (4) managing and making accessible the large number of databases, maps, and other products being developed; and (5) coordinating efforts among WLCI partners, helping them use USGS-developed decision-support tools, and integrating WLCI outcomes with future habitat enhancement and research projects. The new (completed) project was the development and publication of a public outreach piece for visitors of Fossil Butte National Monument. The final product was a USGS Fact Sheet that capitalized on previously collected elk-monitoring data to interpret the ecology of the Monument&rsquo;s elk population and the importance of the Monument&rsquo;s habitats to this highly visible wildlife species. One of the completed projects entailed developing and evaluating a synthetic approach to high-resolution satellite imagery for use in effectiveness monitoring, which culminated in a journal article. The other completed project was a coalescing of two similar tasks under data and information management that pertain to Web application development and the development of outreach and graphic products into a single integrated project that focuses on developing and maintaining/upgrading Web applications and other tools for visualizing, mapping, and using geospatial data. Major accomplishments for FY2014 included several publications, including Part B of an energy resources map that (with Part A) depicts coal, wind, oil, gas, oil shale, uranium, and solar energy production in the WLCI region. Two published works associated with sage-grouse included a Wildlife Monograph on prioritizing species&rsquo; habitats across large landscapes, multiple seasons, and novel areas (using sage-grouse in Wyoming as an example), and a USGS Data Series report that includes both the data used in the habitat-prioritization models and the habitat prioritization models developed for sage-grouse. Our Science Team also published a framework for conducting large, collaborative projects that rely on geospatial data, and a paper that describes the efficacy of fusing satellite data collected at various resolutions for measuring and monitoring vegetation changes. These products are all invaluable tools for maximizing the efficiency and effectiveness of managing species of concern, conducting future landscape-scale assessments, and monitoring status and trends of landscape conditions. Other highlights of FY2014 included a renewed effort to gather and analyze wildlife and habitat status and trend data for the WLCI Interagency Monitoring Database (IAMD) to assess long-term trends and cumulative effects associated with land-use and climate changes. Water-monitoring efforts included drilling four new groundwater-monitoring wells in the Green and New Fork River basins near the proposed Normally Pressured Lance Formation energy development, and continued data collection at established water-monitoring sites. Three additional wells were sampled as part of the Wyoming Groundwater Monitoring Network, bringing the total to 19 Network wells sampled in the WLCI region since 2010. Combined, these water-monitoring efforts can help to identify potential changes in water quality or levels that may result from land-use changes. Major terrestrial monitoring accomplishments included processing satellite imagery from 1985&minus;2010 to develop a historical perspective of long-term vegetation changes, which can serve as a basis for monitoring current and future trends in sagebrush steppe. Such data are crucial tools for agencies tasked with sage-grouse management and conservation. The USGS WLCI Science Team also continued monitoring and testing methods for evaluating WLCI habitat treatments designed to promote aspen regeneration and enhance sage-grouse habitat, and to assess how those treatments influence invasive species distributions and ungulate herbivory. Highlights included analyzing field data collected to elucidate the relationships between sage-grouse habitat use and the proximity of energy infrastructure, and using new instruments to measure productivity responses of aspen woodlands to various factors. Numerous FY2014 accomplishments specifically addressed agency needs to manage and conserve Wyoming&rsquo;s wildlife species of concern. A pygmy rabbit habitat model and Wyoming distribution map were completed to identify factors associated with rabbit habitat occupancy. Previous work on sage-grouse population dynamics was expanded to better understand the factors that drive long-term viability of sage-grouse populations and to develop a tool that helps to identify key factors limiting sage-grouse persistence in Wyoming. Field work and data analyses continued for elucidating the relationships between sagebrush songbird abundance and productivity, the intensity of energy development, and community dynamics of nest predators. For the mule deer study, mixed mountain shrublands important to migrating and wintering mule deer were mapped and delivered to WLCI partners. Additionally, the relationships between energy development and crucial winter habitat for mule deer were evaluated, and a new phase of work was implemented to better understand relationships between plant phenology and mule deer migration movements. Finally, initial analyses of data collected to evaluate fish-community composition in relation to habitat quality indicate that water quality, as measured by concentrations of hydrocarbons, water temperature, and others parameters, has been diminished in subwatersheds with higher levels of energy development. Overall, the outcomes and products of these wildlife studies contribute significantly to the information and tools needed for addressing effects of land-use changes on Wyoming&rsquo;s species of concern. Finally, capabilities of the WLCI Web site and the USGS ScienceBase infrastructure were maintained and upgraded to help ensure access to and efficient use of all the WLCI data, products, assessment tools, and outreach materials that have been developed. Of particular note is the completion of three Web applications developed for mapping (1) the 1900&minus;2008 progression of oil and gas development;(2) the predicted distributions of Wyoming&rsquo;s Species of Greatest Conservation Need; and (3) the locations of coal and wind energy production, sage-grouse distribution and core management areas, and alternative routes for transmission lines within the WLCI region. Collectively, these applications tools provide WLCI planners and managers with powerful tools for better understanding the distributions of wildlife species and potential alternatives for energy development.

Wyoming↗

Spatial and temporal variability of harmful algal blooms in Milford Lake, Kansas, May through November 2016

The U.S. Geological Survey, in cooperation with the Kansas Department of Health and Environment (KDHE), completed a study to quantify the spatial and temporal variability of cyanobacterial blooms in Milford Lake, Kansas, over a range of environmental conditions at various time scales (hours to months). A better understanding of the spatial and temporal variability of cyanobacteria and microcystin will inform sampling and management strategies for Milford Lake and for other lakes with cyanobacterial harmful algal bloom (CyanoHAB) issues throughout the Nation. Spatial and temporal variability were assessed in the upstream one-third of Milford Lake (designated as “Zone C” by KDHE) during May through November 2016 using a combination of time-lapse photography, continuous water-quality monitors, discrete phytoplankton, chlorophyll, and microcystin samples, and spatially dense near-surface data. Combined, these data were used to characterize variability of cyanobacterial abundance, algal biomass, and microcystin concentrations in Zone C of Milford Lake before, during, and after cyanobacterial blooms in 2016. Temporal patterns were evaluated during May through November 2016 using time-lapse photography at six locations in Zone C and at a single point location (the Wakefield site) using a combination of discrete and continuously measured water-quality data (including the cyanobacterial pigment phycocyanin). Based on time-lapse photography, CyanoHABs developed in Zone C of Milford Lake in early July and persisted through the end of November. Bloom accumulations at individual sites were dependent on wind direction. After a change in wind direction, it would take about 1 day for accumulations to become visible at different locations. During periods with low wind, accumulations were widespread and visible at all sites. Cyanobacteria were absent from the algal community at the Wakefield site in late May and were a minor component of the community in June; however, by mid-July the cyanobacteria were dominant and remained dominant until early November. Chlorophyll and microcystin concentrations at the Wakefield site were estimated using sensor-measured phycocyanin based on regression models developed for Zone C. Regression-estimated concentrations likely are more indicative of seasonal patterns in algal biomass (as indicated by chlorophyll concentrations) and microcystin than discretely collected samples because regression-estimated data have a much higher temporal resolution. Based on regression estimates, algal biomass and microcystin concentrations at the Wakefield site steadily increased from May through August. After August, concentrations decreased but remained relatively high compared to May and June. Daily chlorophyll maxima were as much as 400 times higher than daily minima, and daily microcystin maxima were as many as several orders of magnitude higher than daily minima. The extreme variability in algal biomass and microcystin concentrations at the Wakefield site reflects the development and dissipation of blooms, as indicated by the time-lapse cameras. Based on regression-estimated microcystin concentrations, the KDHE watch and warning thresholds for microcystin were exceeded during mid-June through late November. Exceedance of KDHE advisory thresholds often changed from no advisory to watch or warning over the course of the day because of the variability in algal biomass and microcystin concentrations caused by bloom development and dissipation. Continuous water-quality monitors may be useful in informing public-health decisions in lakes with variable CyanoHAB conditions; however, site-specific models need to be developed, and best practices for using continuous water-quality monitors to inform CyanoHAB management strategies need to be established. Spatial data were collected on May 26, July 21, and September 15, 2016, using a combination of a boat-mounted array and discrete water-quality samples analyzed for phytoplankton community composition and chlorophyll and microcystin concentrations. Spatial patterns were described using regression-estimated chlorophyll and microcystin concentrations. During the May 26, 2016, spatial surveys, cyanobacterial abundances were relatively low throughout Zone C and did not exceed KDHE guidance values compared to spatial surveys on July 21 and September 15. Regression-estimated chlorophyll concentrations were indicative of higher algal biomass uplake in Zone C, and decreases in the downlake direction towards Zone B. Regression-estimated chlorophyll concentrations also were more variable uplake than downlake. Based on regression estimates, microcystin concentrations did not exceed KDHE guidance values anywhere in Zone C on May 26. Spatial patterns in microcystin throughout Zone C did not match patterns in regression-estimated chlorophyll concentrations, likely because the algal community was not dominated by cyanobacteria at most locations in May. During the July 21, 2016, spatial surveys, cyanobacterial abundances in Zone C exceeded KDHE guidance values in 50 percent of samples. The algal community in Zone C was dominated by cyanobacteria at all locations except two, where cyanobacteria codominated with diatoms. Both locations where cyanobacteria and diatoms codominated were north of the causeway. Regression-estimated chlorophyll concentrations were indicative of higher algal biomass north of the causeway and on the eastern shore of Zone C. On July 21, algal biomass did not always decrease in the downlake direction. There was a decrease just south of the causeway but an increase shortly after with higher concentrations into Zone B. Spatial maps indicated changes in algal distribution at a 0.5-meter depth, with algae moving to the central part of the lake north of the causeway and along the eastern shore south of the causeway. Most regression-estimated microcystin concentrations on July 21 exceeded KDHE guidance values, reflecting the pervasive bloom conditions in Zone C during this period. Spatial patterns in regression-estimated microcystin concentrations throughout Zone C were similar to patterns seen in discrete samples and regression-estimated chlorophyll concentrations, with higher concentrations north of the causeway and on the east shore of Zone C. During the September 15, 2016, spatial surveys, cyanobacterial abundances did not exceed KDHE guidance values. The algal community north of the causeway was dominated by diatoms. The algal community throughout the rest of Zone C was dominated by cyanobacteria. Of regression-estimated microcystin concentrations on September 15, 80 percent did not exceed KDHE guidance values. Spatial patterns indicated northward movement of the cyanobacterial bloom consistent with a wind shift noted the previous day. On September 14, winds were generally from the north to northwest, shifting to the south by September 15. There was a northward progression of chlorophyll and microcystin during the spatial surveys. These data, along with the camera data and spatial and wind data from May and July, indicate that wind can be a major driver of the spatial and temporal variability of cyanobacterial blooms in Milford Lake and likely plays a role in the extent and duration of near-shore accumulations.

Kansas↗