Search USGS⌕ Search

SEARCH · Search USGS

Results for “Final Report”

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.

At least 1,549 records · Page 86Linked to original sources

The coral reef of South Moloka'i, Hawai'i— Portrait of a sediment-threatened fringing reef

Moloka‘i, with the most extensive coral reef in the main Hawaiian Islands, is especially sacred to Hina, the Goddess of the Moon. As Hinaalo, she is the Mother of the Hawaiian people; as Hinapuku‘a, she is the Goddess of Fishermen; and in the form Hina‘opuhalako‘a, she is the Goddess who gave birth to coral, coral reefs, and all spiny marine organisms. Interdependence between the reef’s living resources, the people, and their cosmology was the basis for management of Moloka‘i’s coastal waters for over a thousand years. The ancient residents of Moloka‘i built the greatest concentration of fishponds known anywhere, but their mastery of mariculture, something needed now more than ever, was lost after near genocide from exotic Western diseases. Subsequent destruction of the native vegetation for exotic cattle, goats, pigs, sugar cane, and pineapple caused soil erosion and sedimentation on the reef flat. This masterful volume clearly documents that soil washing into the sea is the major threat to the reef today. Abandoned fishponds, choked with sediment, now act as barriers and mud traps, making damage to corals less than it would otherwise would have been. The role of mud and freshwater from land in preventing coral reef growth, clearly articulated in Charles Darwin’s first book, The Structure and Distribution of Coral Reefs, is the major theme of this book. All around the tropics, coral reefs have died from huge increases in terrestrial sedimentation that resulted from destruction of hillside forests for cash-crop agriculture and pastures in the colonial era, especially in Latin America, Asia, and the islands of the Caribbean and Indo-Pacific. It is obvious that one cannot manage the coastal zone as a unit separate from the watersheds that drain into it. Yet there has been surprisingly little comprehensive scientific study of these impacts. In this landmark volume, U.S. Geological Survey researchers and their colleagues have developed and applied a remarkably integrated approach to the reefs of Moloka‘i, combining geology, oceanography, and biology to provide an in-depth understanding of the processes that have made these reefs grow and that now limit them. They have joined old fashioned natural history of marine animals and plants with study of the geological evolution of the island, hydrology, meteorology, and land-use history, to an arsenal of new methods of remote sensing, including aerial photography, laser ranging, infrared thermal mapping, seismic reflection, in-situ instrumentation to measure chemical parameters of water quality, and direct measurements of the physical driving forces affecting them—such as wave energy, currents, sedimentation, and sediment transport. They provide a level of documentation and insight that has never been available for any reef before. A remarkable feature of this book is that it is aimed at the people of Moloka‘i to inform them of what is happening to their reef and what they might do to preserve their vital resources. The scientific data and interpretations are expressed in unusually clear and comprehensible language, free of the professional jargon that makes most technical publications impenetrable to the public that most needs to know about them, yet without loss of scientific rigor. Here readers will see clearly explained the whole path of soil loss, from the impacts of wild pigs and goats at higher elevations, deforestation of the hills for cattle pasture at lower levels, and denudation of low lands for cash crops. The resulting biological impoverishment has bared the soils, which wash away in flash storms, smothering the inshore reefs, whose growth was already limited because they had grown right up to sea level. The data in this book show that the mud doesn’t get far if it is washed into the sea during a big storm with heavy waves. Afterwards this mud keeps getting stirred up by every succeeding storm, spreading and affecting corals over wider areas until it is finally washed out of the system—and that only happens if there is no more new mud washing onto the reef. I saw this myself a few years ago in Pila‘a Bay on Kaua‘i, where a bulldozed hillside of abandoned sugar cane fields had slumped right on top of a coral reef following exceptional rains. Years later, the algae species were zoned in a way that clearly mapped the distribution of nutrients washed into the bay, most likely from fertilizers bound to the eroded soils. That pattern closely mimics, on a small scale, that shown in Moloka‘i in this volume, where the inner reef is covered with algae, zoned by species in a way that points to land-based sources of nutrients, while the outermost reef slope is still coral dominated, and the deep algae seem to indicate deep-water nutrient upwelling. What of the future? The Hawaiian Islands have been exceptionally fortunate to be spared the worst coral heatstroke death from high temperatures, at least to date. So far, the worst global warming impacts have luckily been small in this region, and the small number of people on Moloka‘i has kept population densities, and sewage pollution, low compared to the more developed islands. Nutrients from years of sugar and pineapple fertilization, and the washing of this soil onto the reefs, show clear influences on the pattern of algae on the reef. Even at very low levels of nutrients, well below that which drives algae to smother and kill coral reefs, more algae is present. Soil erosion control is therefore the key to better management of both nutrients and turbidity on Moloka‘i reefs. To that end land management actions mentioned in this book, such as suppressing wild fires and eliminating wild goats and pigs, could be made even more effective if supplemented by active erosion control using plants whose roots bind the soil effectively in place. Through all of these efforts, Hina and the people of Moloka‘i could be happy again!

Hawaii↗

Behavior and movements of adult spring Chinook salmon ( Oncorhynchus tshawytscha ) in the Chehalis River Basin, southwestern Washington, 2015

Recent interest in flood control and restoration strategies in the Chehalis River Basin has increased the need to understand the current status and ecology of spring Chinook salmon. Based on the extended period between freshwater entry and spawn timing, spring Chinook salmon have the longest exposure of all adult Chinook salmon life histories to the low-flow and high water temperature conditions that typically occur during summer. About 100 adult spring Chinook salmon were found dead in the Chehalis River in July and August 2009. Adult Chinook salmon are known to hold in cool-water refugia during warm summer months, but the extent to which spring Chinook salmon might use thermal refugia in the Chehalis River is unknown. The movements and temperature exposures of adult spring Chinook salmon following their return to the Chehalis River were investigated using radiotelemetry and transmitters equipped with temperature sensors, combined with water temperature monitoring throughout the basin. A total of 23 spring Chinook salmon were radio-tagged between April and early July 2015; 11 were captured and released in the main-stem Chehalis River, and 12 were captured and released in the South Fork Newaukum River. Tagged fish were monitored with a combination of fixed-site monitoring locations and regular mobile tracking, from freshwater entry through the spawning period. Water temperature and flow conditions in the main-stem Chehalis River during 2015 were atypical compared to historical averages. Mean monthly water temperatures between March and July 2015 were higher than any decade since 1960 and mean daily flows were 30–70 percent of the flows in previous years. Overall, 96 percent of the tagged fish were detected, with a mean of 62 d in the detection history of tagged fish. Of the 11 fish released in the main-stem Chehalis River, six fish (55 percent) moved upstream, either shortly after release (2–7 d, 50 percent), or following a short delay (12–18 d, 50 percent). One fish released in the main-stem Chehalis River remained near the release location for 64 d before moving upstream. The final fates for the seven fish that moved upstream in the main-stem Chehalis River included two fish with unknown fates, two fish with a fate of pre-spawn mortality, and three fish that were assigned a fate of spawner. Four (36 percent) of the radio-tagged Chinook salmon released in the main-stem Chehalis River showed limited movement from their release sites, and were assigned fates of unknown (one fish), pre-spawn mortality (one fish), and spit/mortality (2 fish). The 12 spring Chinook salmon released in the South Fork Newaukum River remained in the South Fork Newaukum River throughout the study period. Five (42 percent) of these fish were actively moving through the spawning period and were assigned a fate of spawner. Seven (58 percent) of these fish were detected for a period following release, but their detection histories ended prior to the spawning period. The fates assigned to these seven fish included two fish with spit/mortality fates and five fish with fates of pre-spawn mortality. Tagged fish in both the Chehalis River and the South Fork Newaukum River showed limited movements during the peak water temperatures in July and August, and were not frequently detected at sites where water temperatures were greater than 21 °C. Pre-spawn mortality due to predation or harvest may be an important factor in the Chehalis River Basin as it was the assigned fate for 27 percent of the fish released in the main-stem Chehalis River and 42 percent of the fish released in the South Fork Newaukum River. This study represents a substantial contribution to the understanding of spring Chinook salmon in the Chehalis River Basin. The water temperatures and flow conditions during the 2015 study period were not typical of the historical conditions in the basin and the numbers of tagged fish monitored was relatively low, so results should be interpreted with those cautions in mind.

Washington↗

Hydrogeology and simulated groundwater availability in reaches 3 and 4 of the Washita River aquifer, southern Oklahoma, 1980–2017

The 1973 Oklahoma Groundwater Law (Oklahoma Statutes §82–1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s aquifers to determine the maximum annual yield for each groundwater basin. Because more than 20 years have elapsed since the final order was issued, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, conducted an updated hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater flow and availability in reaches 3 and 4 of the Washita River aquifer in southern Oklahoma for a study period spanning 1980–2017. A hydrogeologic framework and conceptual model were developed to guide the construction and calibration of a numerical model of the Washita River aquifer. The numerical model was calibrated to water-table-altitude observations at selected wells, base-flow observations at selected U.S. Geological Survey streamgages, and the conceptual-model recharge. Three types of groundwater-availability scenarios were run using the calibrated numerical model. These scenarios were used to (1) estimate equal-proportionate-share pumping rates, (2) quantify the potential effects of projected well withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical 10-year drought. With Washita River main-stem inflows, the 20-, 40-, and 50-year equal-proportionate-share pumping rates under normal recharge conditions were about 3.08 acre-feet per acre per year for reach 3 and about 3.80 acre-feet per acre per year for reach 4. Projected 50-year pumping scenarios were used to simulate the effects of modified well withdrawal rates. Because well withdrawals were less than 1 percent of the calibrated numerical-model water budget, changes to the well pumping rates had little effect on Washita River base flows and groundwater storage in the Washita River aquifer. A hypothetical 10-year drought scenario was used to simulate the potential effects of a prolonged period of reduced recharge on groundwater storage. Groundwater storage at the end of the drought period was 4.6 percent less than the groundwater storage of the calibrated numerical model at the end of the drought period.

Oklahoma↗

Geology and mining history of the Southeast Missouri Barite District and the Valles Mines, Washington, Jefferson, and St. Francois Counties, Missouri

The Southeast Missouri Barite District and the Valles Mines are located in Washington, Jefferson, and St. Francois Counties, Missouri, where barite and lead ore are present together in surficial and near-surface deposits. Lead mining in the area began in the early 1700’s and extended into the early 1900’s. Hand mining of lead in the residuum resulted in widespread pits (also called shafts or diggings), and there was some underground mining of lead in bedrock. By the 1860’s barite was recovered from the residuum by hand mining, also resulting in widespread diggings, but generally not underground mines in bedrock. Mechanized open-pit mining of the residuum for barite began in the 1920’s. Barite production slowed by the 1980’s, and there has not been any barite mining since 1998. Mechanized barite mining resulted in large mined areas and tailings ponds containing waste from barite mills. The U.S. Environmental Protection Agency (EPA) has determined that lead is present in surface soils in Washington and Jefferson Counties at concentrations exceeding health-based screening levels. Also, elevated concentrations of barium, arsenic, and cadmium have been identified in surface soils, and lead concentrations exceeding the Federal drinking-water standard of 15 micrograms per liter have been identified in private drinking-water wells. Potential sources of these contaminants are wastes associated with barite mining, wastes associated with lead mining, or unmined natural deposits of barium, lead, and other metals. As a first step in helping EPA determine the source of soil and groundwater contamination, the U.S. Geological Survey (USGS), in cooperation with the EPA, investigated the geology and mining history of the Southeast Missouri Barite District and the Valles Mines. Ore minerals are barite (barium sulfate), galena (lead sulfide), cerussite (lead carbonate), anglesite (lead sulfate), sphalerite (zinc sulfide), smithsonite (zinc carbonate), and chalcopyrite (copper-iron sulfide). The Cambrian Potosi Dolomite is the most important formation for the ore deposits, followed by the Eminence Dolomite. Because galena, sphalerite, and barite are less soluble than dolomite, chemical weathering of the ore-bearing dolomite bedrock resulted in the concentration of ore minerals in the residuum. Most of the barite and lead mining was in the residuum, which averages 10 to 15 feet thick. Lead mining by French explorers may have begun in 1719 along Old Mines Creek at Cabanage de Renaudiere, which was followed shortly by the discovery of lead and the development of lead mines at Mine Renault (also called Forche a Renault Mine), Old Mines, and at other places along the Big River, Mineral Fork, and Forche a Renault Creek. Lead mining began sometime between 1775 and 1780 at Mine a Breton, the name of which was later changed to Potosi. Other mining areas were developed in the early part of the 19th century, including Fourche a Courtois (Palmer Mines), the French Diggings, and the Richwoods Mines. Zinc became a valuable resource after the Civil War, and the Valles Mines was an important supplier of zinc as well as lead, with at least some production up until the 1920’s. Lead mining declined in the early part of the 20th century as mining in the Old Lead Belt, Mine La Motte, and the Tri-State District expanded. The earliest lead mines were diggings in the residuum and were round holes (shafts) about 4 feet in diameter dug with pick and shovel about 15–20 feet deep, with drifts dug a short distance laterally from the bottom of the shafts. This mining process was repeated a short distance away until a large area was covered with pits. Some mining in bedrock began by about 1800, with shafts as deep as 170 feet and as much as several hundred feet of lateral drifts. Smelting of the lead ore to elemental lead was first done using a log furnace, which was inefficient; estimates have been made that only about 50 percent of the lead was recovered, and the remainder was lost to the ashes (slags) and to volatilization. Starting in 1798, ash furnaces were used to smelt the ashes from the log furnaces. These two furnaces were worked in tandem for many years but were gradually replaced by other furnaces, including the Scotch hearth. Estimates of lead recovery as high as 80–90 percent have been made for the Scotch hearth. By the mid-1870’s the air furnace was being used, also with estimated lead recovery as high as 80–90 percent. Zinc furnaces were built when zinc became a valuable commodity, but much of the zinc ore was shipped out of the area, either to a smelter in St. Louis, Missouri, or to other smelters. The total lead and zinc production from the Southeast Missouri Barite District and the Valles Mines is estimated at 180,000 tons of lead and 60,000 tons of zinc. An estimated 97,000 tons of lead and an estimated 120,000 tons of zinc were lost during smelting. The estimated losses do not include losses at the mine site during mining and preparation for smelting, such as the loss of fine-grained galena during hand cleaning or the discarding of zinc ore before its value was known, for which no estimates are available. Hand mining for barite in the residuum was active by at least the 1860’s and peaked from 1905 to the 1930’s when several thousand people were engaged in barite mining. Hand mining (diggings) and cleaning of the ore was done in much the same way as earlier lead mining, with the additional use of a rattle box to further clean the barite. Mechanized open-pit mining of old barite diggings began in 1924 to recover barite left behind by hand mining, and washing plants were used to clean the clay from the barite. Hand mining, however, continued to thrive, and washer plants began to close temporarily in 1931; nearly all of the barite produced before 1937 was by hand mining. By the 1940’s, however, all barite mining was mechanized. Mechanized mining used shovels powered by steam, gasoline, or electricity (and by the 1950’s draglines and front-end loaders) to mine the residuum. The ore was loaded onto rail cars (and by the 1940’s, trucks) for shipment to washer plants. Clay was removed from the barite using a log washer, and a jig was used to concentrate the barite. Overflow from the log washers was waste and went to a mud (tailings) pond. The coarse jig tailings went to tailings piles or were used as railroad ballast and, later, to create roads within the mine pit. Some barite was ground, depending on its final use, and some ground barite was bleached using a hot solution of sulfuric acid to remove impurities such as iron minerals and lead sulfide (galena). An earlier bleaching process used lead-lined tanks. Large quantities of water were required for milling the barite; some was recirculated water and the remainder came from dammed streams or was pumped from wells. Tailings and wastewater were impounded behind dikes that were built across small valleys and were increased in height as necessary using washer waste and any overburden that had been stripped. In some cases, dikes were built across valleys that had already been mined for barite. The total production of barite from the Southeast Missouri Barite District and the Valles Mines is estimated to have been about 13.1 million tons. Most of the barite production was from Washington County. Hand mining and processing of barite was inefficient. Estimates of barite recovery range from less than one-fourth to about one-half because pillars between the shafts in the residuum needed to be left unmined for stability. With mechanized mining, large amounts of barite were lost during the milling process. It has been estimated that about 30 percent of the barite was lost and that about two-thirds of the lost barite was fine-grained and was discharged to the tailings ponds. Some galena was lost to the tailings ponds. A 1972 inventory of tailings ponds by the Missouri Geological Survey identified 67 ponds in the Southeast Missouri Barite District (there are more than this currently documented). Results from samples from four ponds that were drilled were used to estimate that the 67 ponds contained almost 39 million tons (or cubic yards) of tailings averaging about 5 percent barite, for a potential reserve of 1.935 million tons of barite. It is not known how much lead was removed during barite mining, either by hand or mechanized mining and processing, how much lead was recovered, or how much lead went as fines to the tailing ponds or as coarse material to mine roads or was otherwise lost.

Missouri↗

Assessing the relationship between Section 404 and wetland losses: a feasibility study

The primary objective of the Clean Water Act of 1977 (33 U.S.C. 1251) is to restore and maintain the physical, chemical, and biological integrity of the Nation's waters. Section 404 of the Act regulates the discharge of dredged or fill materials into wetlands and represents the primary Federal authority for regulation of wetland alterations. Since its inception, the Section 404 program has been controversial in regard to the extent to which it was intended to provide wetlands regulation. Section 404 requires those who wish to discharge dredged or fill material into waters of the United States, which include many wetlands, to first obtain a Federal permit. The Environmental Protection Agency (EPA) has overall responsibility for administration of the Section 404 program and promulgates guidelines that must be followed in issuing permits. In addition, EPA has the final authority to prohibit specific discharges if the environmental impacts are unacceptable. The U.S. Army Corps of Engineers (Corps) issues Section 404 permits, which can be of two types. Individual Permits are issued following case-by-case reviews of proposed discharges. General Permits, which can be either nationwide or regional in scope, are authorized by the Corps for categories of activities that are similar in nature and that have only minimal individual and cumulative adverse environmental impacts. EPA, the National Marine Fisheries Service (NMFS), the Fish and Wildlife Service (Service), and State natural resource agencies review and comment on permit applications and offer recommendations on appropriate mitigation measures. Although comments from the Service and other natural resource agencies are advisory in nature (EPA's veto authority excepted), they can serve as the basis for modifying, conditioning, or denying a Section 404 permit. In 1986, in a survey conducted by the National Ecology Research Center, Service personnel indicated interest in additional information concerning both wetland trends and the impacts of activities authorized by the Corps under Section 404 (Roelle 1986). Although there is some information concerning wetland losses for certain geographic areas and for the Nation as a whole (Frayer et al. 1983; Tiner 1984), there appears to be little information on how these losses relate to the Section 404 permitting process. The primary objective of this study was to determine the feasibility of estimating wetland losses in relationship to Individual and General Permits issued under Section 404. A secondary objective was to assemble data on acceptance and implementation of specific mitigation recommendations offered by the Service and other natural resource agencies in connection with development activities on wetlands examined. At present, 26 categories of Nationwide permits have been authorized by the Corps. Nationwide permit 26 was of particular interest in this study because it specifically authorizes discharges into wetlands under certain circumstances (see Figure 1 and related text for a more complete discussion of circumstances under which Nationwide Permit 26 is applicable). All subsequent references to Nationwide or General permits pertain to Nationwide Permit 26.

Report↗

A field guide to the central, creeping section of the San Andreas fault and the San Andreas Fault Observatory at Depth

This field trip is along the central section of the San Andreas fault and consists of eight stops that illustrate surface evidence of faulting, in general, and features associated with active fault creep, in particular. Fault creep is slippage along a fault that occurs either in association with small-magnitude earthquakes or without any associated large-magnitude earthquakes. Another aspect of the trip is to highlight where there are multiple fault traces along this section of the San Andreas fault zone in order to gain a better understanding of plate-boundary processes. The first stop is along the Calaveras fault, part of the San Andreas fault system, at a location where evidence of active fault creep is abundant and readily accessible. The stops that follow are along the San Andreas fault and at convenient locations to present and discuss rock types juxtaposed across the fault that have been transported tens to hundreds of kilometers by right-lateral motion along the San Andreas fault. Stops 6 and 7 are examples of recent studies of different aspects of the fault: drilling into the fault at the depth of repeating magnitude (M) 2 earthquakes with the San Andreas Fault Observatory at Depth (SAFOD) and the geological, geophysical, and seismological study of M 6 earthquakes near the town of Parkfield. Along with the eight official stops on this field trip are 12 “rolling stops”—sites of geologic interest that add to the understanding of features and processes in the creeping section of the fault. Many of the rolling stops are located where stopping is difficult to dangerous; some of these sites are not appropriate for large vehicles (buses) or groups; some sites are not appropriate for people at all. We include photographs of or from many of these sites to add to the reader's experience without adding too many stops or hazards to the trip. An extensive set of literature is available for those interested in the San Andreas fault or in the creeping section, in particular. For more scientifically oriented overviews of the fault, see Wallace ( 1990 ) and Irwin ( 1990 ); for a more generalized overview with abundant, colorful illustrations, see Collier ( 1999 ). Although the presence of small sections of the San Andreas fault was known before the great 1906 San Francisco earthquake, it was only after that event and subsequent geologic investigations reported in Lawson ( 1908 ) that showed the fault as a long structure, extending all the way from east of Los Angeles into northern California. Prentice ( 1999 ) described the importance of the 1908 “Lawson report” and how it pivotally influenced the understanding of the San Andreas. Hill ( 1981 ) presented a wonderful introduction to the evolution of thought on the San Andreas. Geologic maps and maps of the most recently active fault trace in the creeping section, or large parts of it, include those by Brown ( 1970 ), Dibblee ( 1971 , 1980 ), and Wagner et al. ( 2002 ); detailed geologic maps are discussed at various stops in this guide. Various aspects of the creeping section of the San Andreas fault have been the focus of many geologic field trips in the past few decades. Guidebooks for some of those trips include those by Gribi ( 1963a , 1963b ), Brabb et al. ( 1966 ), Rogers ( 1969 ), Bucknam and Haller ( 1989 ), Harden et al. ( 2001 ), and Stoffer ( 2005 ). The creeping section of the San Andreas fault zone lies between areas that experienced large-displacement surface breakage during great earthquakes in 1857 and 1906 (Fig. 1 inset). Burford and Harsh ( 1980 ) divided the creeping section into three segments: (1) a northwest section where the creep rate increases to the southeast in step-like increments, (2) a central section where the creep rate is relatively constant at a maximum value of ∼30 mm/yr (∼1.2 in/yr), and (3) a southeast section where the creep rate decreases to the southeast (Fig. 2 ). The rate of slip along the creeping section of the fault zone has been measured using creepmeters, alignment arrays, and laser distance-measuring devices. The aperture of measurements over which these measurements are made ranges from 10 m (∼33 ft) (creepmeters) to 100 m (∼330 ft) (alignment arrays) to kilometers and tens of kilometers (laser measuring devices). Creepmeter and alignment-array measurements are here termed “near-fault” measurements; laser measurements over distances of 1–2 km (∼0.6–1.2 mi) are termed “intermediate-scale” measurements; laser measurements over tens of kilometers (miles) are termed “broadscale” measurements. Comparisons among near-fault, intermediate-scale, and broadscale measurements and geologic maps show that the northwest part of the creeping section of the fault is composed of two narrow zones of active deformation, one along the San Andreas fault and one along the Calaveras-Paicines fault, whereas the central and southeast sections are both composed of a single relatively narrow zone of deformation. The southeast section is transitional to a locked zone southeast of Cholame; a locked fault is one that slips only in association with a moderate to large earthquake. Throughout the creeping section of the San Andreas fault zone, broadscale measurements generally indicate more deformation than near-fault and intermediate-scale measurements, which are in reasonably close agreement except at Monarch Peak (Mustang Ridge), near the center of the creeping section and our Stop 5 ((Figs. 1 ) and 2 ). Features that we see on this trip include offset street curbs, closed depressions (sag ponds), fault scarps (steep slopes formed by movement along a fault), a split and displaced tree, offset fence lines, fresh fractures, and offset road lines (Fig. 3 is a sketch showing some of the landforms that represent deformation by an active fault). We also see evidence of long-term maturity of the San Andreas fault, as indicated by fault features and displaced rock types (Fig. 4 ). Finally, we will visit sites of ongoing research into the processes associated with earthquakes and their effects. Discussions include drilling into the San Andreas fault at the SAFOD drill site and the 2004 Parkfield earthquake and its effects and implications.

California↗

Magnitude of flood flows for selected annual exceedance probabilities for streams in Massachusetts

The U.S. Geological Survey, in cooperation with the Massachusetts Department of Transportation, determined the magnitude of flood flows at selected annual exceedance prob­abilities (AEPs) at streamgages in Massachusetts and from these data developed equations for estimating flood flows at ungaged locations in the State. Flood magnitudes were deter­mined for the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent AEPs at 220 streamgages, 125 of which are in Massachusetts and 95 are in the adjacent States of Connecticut, New Hamp­shire, New York, Rhode Island, and Vermont. AEP flood flows were computed for streamgages using the expected moments algorithm weighted with a recently computed regional skew­ness coefficient for New England. Regional regression equations were developed to estimate the magnitude of floods for selected AEP flows at ungaged sites from 199 selected streamgages and for 60 potential explanatory basin characteristics. AEP flows for 21 of the 125 streamgages in Massachusetts were not used in the final regional regression analysis, primarily because of regulation or redundancy. The final regression equations used general­ized least squares methods to account for streamgage record length and correlation. Drainage area, mean basin elevation, and basin storage explained 86 to 93 percent of the variance in flood magnitude from the 50- to 0.2-percent AEPs, respec­tively. The estimates of AEP flows at streamgages can be improved by using a weighted estimate that is based on the magnitude of the flood and associated uncertainty from the at-site analysis and the regional regression equations. Weighting procedures for estimating AEP flows at an ungaged site on a gaged stream also are provided that improve estimates of flood flows at the ungaged site when hydrologic characteristics do not abruptly change. Urbanization expressed as the percentage of imperviousness provided some explanatory power in the regional regression; however, it was not statistically significant at the 95-percent confidence level for any of the AEPs examined. The effect of urbanization on flood flows indicates a complex interaction with other basin characteristics. Another complicating factor is the assumption of stationarity, that is, the assumption that annual peak flows exhibit no significant trend over time. The results of the analysis show that stationarity does not prevail at all of the streamgages. About 27 percent of streamgages in Massachusetts and about 42 percent of streamgages in adjacent States with 20 or more years of systematic record used in the study show a significant positive trend at the 95-percent confidence level. The remaining streamgages had both positive and negative trends, but the trends were not statistically significant. Trends were shown to vary over time. In particular, during the past decade (2004–2013), peak flows were persistently above normal, which may give the impression of positive trends. Only continued monitoring will provide the information needed to determine whether recent increases in annual peak flows are a normal oscillation or a true trend. The analysis used 37 years of additional data obtained since the last comprehensive study of flood flows in Massa­chusetts. In addition, new methods for computing flood flows at streamgages and regionalization improved estimates of flood magnitudes at gaged and ungaged locations and better defined the uncertainty of the estimates of AEP floods.

Massachusetts↗

Methods for estimating selected spring and fall low-flow frequency statistics for ungaged stream sites in Iowa, based on data through June 2014

A statewide study was led to develop regression equations for estimating three selected spring and three selected fall low-flow frequency statistics for ungaged stream sites in Iowa. The estimation equations developed for the six low-flow frequency statistics include spring (April through June) 1-, 7-, and 30-day mean low flows for a recurrence interval of 10 years and fall (October through December) 1-, 7-, and 30-day mean low flows for a recurrence interval of 10 years. Estimates of the three selected spring statistics are provided for 241 U.S. Geological Survey continuous-record streamgages, and estimates of the three selected fall statistics are provided for 238 of these streamgages, using data through June 2014. Because only 9 years of fall streamflow record were available, three streamgages included in the development of the spring regression equations were not included in the development of the fall regression equations. Because of regulation, diversion, or urbanization, 30 of the 241 streamgages were not included in the development of the regression equations. The study area includes Iowa and adjacent areas within 50 miles of the Iowa border. Because trend analyses indicated statistically significant positive trends when considering the period of record for most of the streamgages, the longest, most recent period of record without a significant trend was determined for each streamgage for use in the study. Geographic information system software was used to measure 63 selected basin characteristics for each of the 211streamgages used to develop the regional regression equations. The study area was divided into three low-flow regions that were defined in a previous study for the development of regional regression equations. Because several streamgages included in the development of regional regression equations have estimates of zero flow calculated from observed streamflow for selected spring and fall low-flow frequency statistics, the final equations for the three low-flow regions were developed using two types of regression analyses—left-censored and generalized-least-squares regression analyses. A total of 211 streamgages were included in the development of nine spring regression equations—three equations for each of the three low-flow regions. A total of 208 streamgages were included in the development of nine fall regression equations—three equations for each of the three low-flow regions. A censoring threshold was used to develop 15 left-censored regression equations to estimate the three fall low-flow frequency statistics for each of the three low-flow regions and to estimate the three spring low-flow frequency statistics for the southern and northwest regions. For the northeast region, generalized-least-squares regression was used to develop three equations to estimate the three spring low-flow frequency statistics. For the northeast region, average standard errors of prediction range from 32.4 to 48.4 percent for the spring equations and average standard errors of estimate range from 56.4 to 73.8 percent for the fall equations. For the northwest region, average standard errors of estimate range from 58.9 to 62.1 percent for the spring equations and from 83.2 to 109.4 percent for the fall equations. For the southern region, average standard errors of estimate range from 43.2 to 64.0 percent for the spring equations and from 78.1 to 78.7 percent for the fall equations. The regression equations are applicable only to stream sites in Iowa with low flows not substantially affected by regulation, diversion, or urbanization and with basin characteristics within the range of those used to develop the equations. The regression equations will be implemented within the U.S. Geological Survey StreamStats Web-based geographic information system application. StreamStats allows users to click on any ungaged stream site and compute estimates of the six selected spring and fall low-flow statistics; in addition, 90-percent prediction intervals and the measured basin characteristics for the ungaged site are provided. StreamStats also allows users to click on any Iowa streamgage to obtain computed estimates for the six selected spring and fall low-flow statistics.

Iowa↗

Genomics of Arctic cod

The Arctic cod (Boreogadus saida) is an abundant marine fish that plays a vital role in the marine food web. To better understand the population genetic structure and the role of natural selection acting on the maternally-inherited mitochondrial genome (mitogenome), a molecule often associated with adaptations to temperature, we analyzed genetic data collected from 11 biparentally-inherited nuclear microsatellite DNA loci and nucleotide sequence data from from the mitochondrial DNA (mtDNA) cytochrome b (cytb) gene and, for a subset of individuals, the entire mitogenome. In addition, due to potential of species misidentification with morphologically similar Polar cod (Arctogadus glacialis), we used ddRAD-Seq data to determine the level of divergence between species and identify species-specific markers. Based on the findings presented here, Arctic cod across the Pacific Arctic (Bering, Chukchi, and Beaufort Seas) comprise a single panmictic population with high genetic diversity compared to other gadids. High genetic diversity was indicated across all 13 protein-coding genes in the mitogenome. In addition, we found moderate levels of genetic diversity in the nuclear microsatellite loci, with highest diversity found in the Chukchi Sea. Our analyses of markers from both marker classes (nuclear microsatellite fragment data and mtDNA cytb sequence data) failed to uncover a signal of microgeographic genetic structure within Arctic cod across the three regions, within the Alaskan Beaufort Sea, or between near-shore or offshore habitats. Further, data from a subset of mitogenomes revealed no genetic differentiation between Bering, Chukchi, and Beaufort seas populations for Arctic cod, Saffron cod (Eleginus gracilis), or Walleye pollock (Gadus chalcogrammus). However, we uncovered significant differences in the distribution of microsatellite alleles between the southern Chukchi and central and eastern Beaufort Sea samples of Arctic cod. Finally, using ddRAD-Seq data, we identified species-specific markers and in conjunction with mitogenome data, identified an Arctic cod x Polar cod hybrid in western Canadian Beaufort Sea. Overall, the lack of genetic structure among Arctic cod within the Bering, Chukchi and Beaufort seas of Alaska is concordant with the absence of geographic barriers to dispersal and typical among marine fishes. Arctic cod may exhibit a genetic pattern of isolation-by-distance, whereby populations in closer geographic proximity are more genetically similar than more distant populations. As this signal is only found between our two fartherest localities, data from populations elsewhere in the species’ global range are needed to determine if this is a general characteristic. Further, tests for selection suggested a limited role for natural selection acting on the mitochondrial genome of Arctic cod, but do not exclude the possibility of selection on genes involved in nuclear-mitogenome interactions. Unlike previous genetic assessment of Arctic cod sampled from the Chukchi Sea, the high levels of genetic diversity found in Arctic cod assayed in this study, across regions, suggests that the species in the Beaufort and Chukchi seas does not suffer from low levels of genetic variation, at least at neutral genetic markers. The large census size of Arctic cod may allow this species to retain high levels of genetic diversity. In addition, we discovered the presence of hybridization between Arctic and Polar cod (although low in frequency). Hybridization is expected to occur when environmental changes modify species distributions that result in contact between species that were previously separated. In such cases, hybridization may be an evolutionary mechanism that promotes an increase in genetic diversity that may provide species occupying changing environments with locally-adapted genotypes and, therefore, phenotypes. Natural selection can only act on the standing genetic variation present within a population. Therefore, given its higher levels of genetic diversity in combination with a large population size, Arctic cod may be resilient to current and future environmental change, as high genetic diversity is expected to increase opportunities for positive selection to act on genetic variants beneficial in different environments, regardless of the source of that genetic variation.

OCS Study↗

Death Valley regional groundwater flow system, Nevada and California: Hydrogeologic framework and transient groundwater flow model

A numerical three-dimensional (3D) transient groundwater flow model of the Death Valley region was developed by the U.S. Geological Survey for the U.S. Department of Energy programs at the Nevada Test Site and at Yucca Mountain, Nevada. Decades of study of aspects of the groundwater flow system and previous less extensive groundwater flow models were incorporated and reevaluated together with new data to provide greater detail for the complex, digital model. A 3D digital hydrogeologic framework model (HFM) was developed from digital elevation models, geologic maps, borehole information, geologic and hydrogeologic cross sections, and other 3D models to represent the geometry of the hydrogeologic units (HGUs). Structural features, such as faults and fractures, that affect groundwater flow also were added. The HFM represents Precambrian and Paleozoic crystalline and sedimentary rocks, Mesozoic sedimentary rocks, Mesozoic to Cenozoic intrusive rocks, Cenozoic volcanic tuffs and lavas, and late Cenozoic sedimentary deposits of the Death Valley regional groundwater flow system (DVRFS) region in 27 HGUs. Information from a series of investigations was compiled to conceptualize and quantify hydrologic components of the groundwater flow system within the DVRFS model domain and to provide hydraulic-property and head-observation data used in the calibration of the transient-flow model. These studies reevaluated natural groundwater discharge occurring through evapotranspiration (ET) and spring flow; the history of groundwater pumping from 1913 through 1998; groundwater recharge simulated as net infiltration; model boundary inflows and outflows based on regional hydraulic gradients and water budgets of surrounding areas; hydraulic conductivity and its relation to depth; and water levels appropriate for regional simulation of prepumped and pumped conditions within the DVRFS model domain. Simulation results appropriate for the regional extent and scale of the model were provided by acquiring additional data, by reevaluating existing data using current technology and concepts, and by refining earlier interpretations to reflect the current understanding of the regional groundwater flow system. Groundwater flow in the Death Valley region is composed of several interconnected, complex groundwater flow systems. Groundwater flow occurs in three subregions in relatively shallow and localized flow paths that are superimposed on deeper, regional flow paths. Regional groundwater flow is predominantly through a thick Paleozoic carbonate rock sequence affected by complex geologic structures from regional faulting and fracturing that can enhance or impede flow. Spring flow and ET are the dominant natural groundwater discharge processes. Groundwater also is withdrawn for agricultural, commercial, and domestic uses. Groundwater flow in the DVRFS was simulated using MODFLOW-2000, the U.S. Geological Survey 3D finitedifference modular groundwater flow modeling code that incorporates a nonlinear least-squares regression technique to estimate aquifer parameters. The DVRFS model has 16 layers of defined thickness, a finite-difference grid consisting of 194 rows and 160 columns, and uniform cells 1,500 meters (m) on each side. Prepumping conditions (before 1913) were used as the initial conditions for the transient-state calibration. The model uses annual stress periods with discrete recharge and discharge components. Recharge occurs mostly from infiltration of precipitation and runoff on high mountain ranges and from a small amount of underflow from adjacent basins. Discharge occurs primarily through ET and spring discharge (both simulated as drains) and water withdrawal by pumping and, to a lesser amount, by underflow to adjacent basins simulated by constant-head boundaries. All parameter values estimated by the regression are reasonable and within the range of expected values. The simulated hydraulic heads of the final calibrated transient model generally fit observed heads reasonably well (residuals with absolute values less than 10 meters) with two exceptions: in most areas of nearly flat hydraulic gradient the fit is considered moderate (residuals with absolute values of 10 to 20 meters), and in areas of steep hydraulic gradient along the Eleana Range and western part of Yucca Flat, southern part of the Owlshead Mountains, southern part of the Bullfrog Hills, and the north-northwestern part of the model domain (residuals with absolute values greater than 20 meters). Groundwater discharge residuals are fairly random, with as many areas where simulated flows are less than observed flows as areas where simulated flows are greater. The highest unweighted groundwater discharge residuals occur at Death Valley, Sarcobatus Flat (northeastern area), Tecopa, and early observations at Manse Spring in Pahrump Valley. High weighted-discharge residuals were computed in Indian Springs Valley and parts of Death Valley. Most of these inaccuracies in head and discharge can be attributed to insufficient representation of the hydrogeology in the HFM and(or) discharge estimates, misrepresentation of water levels, and(or) model error associated with grid-cell size. The model represents the large and complex groundwater flow system of the Death Valley region at a greater degree of refinement and accuracy than has been possible previously. The representation of detail provided by the 3D digital hydrogeologic framework model and the numerical groundwater flow model enabled greater spatial accuracy in every model parameter. The lithostratigraphy and structural effects of the hydrogeologic framework; recharge estimates from simulated net infiltration; discharge estimates from ET, spring flow, and pumping; and boundary inflow and outflow estimates all were reevaluated, some additional data were collected, and accuracy was improved. Uncertainty in the results of the flow model simulations can be reduced by improving on the quality, interpretation, and representation of the water-level and discharge observations used to calibrate the model and improving on the representation of the HGU geometries, the spatial variability of HGU material properties, the flow model physical framework, and the hydrologic conditions.

California, Nevada↗

Late Holocene volcanism at Medicine Lake Volcano, northern California Cascades

Late Holocene volcanism at Medicine Lake volcano in the southern Cascades arc exhibited widespread and compositionally diverse magmatism ranging from basalt to rhyolite. Nine well-characterized eruptions have taken place at this very large rear-arc volcano since 5,200 years ago, an eruptive frequency greater than nearly all other Cascade volcanoes. The lavas are widely distributed, scattered over an area of ~300 km 2 across the >2,000-km 2 volcano. The eruptions are radiocarbon dated and the ages are also constrained by paleomagnetic data that provide strong evidence that the volcanic activity occurred in three distinct episodes at ~1 ka, ~3 ka, and ~5 ka. The ~1-ka final episode produced a variety of compositions including west- and north-flank mafic flows interspersed in time with fissure rhyolites erupted tangential to the volcano’s central caldera, including the youngest and most spectacular lava flow at the volcano, the ~950-yr-old compositionally zoned Glass Mountain flow. At ~3 ka, a north-flank basalt eruption was followed by an andesite eruption 27 km farther south that contains quenched basalt inclusions. The ~5-ka episode produced two caldera-focused dacitic eruptions. Quenched magmatic inclusions record evidence of intrusions that did not independently reach the surface. The inclusions are present in five andesitic, dacitic, and rhyolitic host lavas, and were erupted in each of the three episodes. Compositional and mineralogic evidence from mafic lavas and inclusions indicate that both tholeiitic (dry) and calcalkaline (wet) parental magmas were present. Petrologic evidence records the operation of complex, multi-stage processes including fractional crystallization, crustal assimilation, and magma mixing. Experimental evidence suggests that magmas were stored at 3 to 6 km depth prior to eruption, and that both wet and dry parental magmas were involved in generating the more silicic magmas. The broad distribution of eruptive events and the relative accessibility and good exposure of lavas, combined with physical and petrologic evidence for multiple and varied mafic inputs, has created an unusual opportunity to understand the workings of this large magmatic system. A combined total of more than 25 intrusive and extrusive events are indicated for late Holocene time. Plutonic inclusions, some with ages as young as Holocene, were also brought to the surface in five of the eruptions. All eruptions took place along northwest- to northeast-trending alignments of vents, reflecting the overall east-west extensional tectonic environment. The interaction of tectonism and volcanism is a dominant influence at this subduction-related volcano, located where the west edge of the extensional Basin and Range Province impinges on the Cascades arc. Ongoing subsidence focused at the central caldera has been documented along with geophysical evidence for a small magma body. This evidence, combined with the frequency of eruptive and intrusive activity in late Holocene time, an active geothermal system, and intermittent long-period seismic events indicate that the volcano is likely to erupt again.

California↗

Large carbon isotope fractionation associated with oxidation of methyl halides by methylotrophic bacteria

The largest biological fractionations of stable carbon isotopes observed in nature occur during production of methane by methanogenic archaea. These fractionations result in substantial (as much as ≈70‰) shifts in δ 13 C relative to the initial substrate. We now report that a stable carbon isotopic fractionation of comparable magnitude (up to 70‰) occurs during oxidation of methyl halides by methylotrophic bacteria. We have demonstrated biological fractionation with whole cells of three methylotrophs (strain IMB-1, strain CC495, and strain MB2) and, to a lesser extent, with the purified cobalamin-dependent methyltransferase enzyme obtained from strain CC495. Thus, the genetic similarities recently reported between methylotrophs, and methanogens with respect to their pathways for C 1 -unit metabolism are also reflected in the carbon isotopic fractionations achieved by these organisms. We found that only part of the observed fractionation of carbon isotopes could be accounted for by the activity of the corrinoid methyltransferase enzyme, suggesting fractionation by enzymes further along the degradation pathway. These observations are of potential biogeochemical significance in the application of stable carbon isotope ratios to constrain the tropospheric budgets for the ozone-depleting halocarbons, methyl bromide and methyl chloride. Methyl bromide (MeBr) and methyl chloride (MeCl) are, respectively, the most abundant volatile bromo- and chlorocarbons in the troposphere and are major contributors to stratospheric ozone destruction (1). Both compounds have natural and human-influenced sources and a predominant sink by reaction with OH in the troposphere (2–4). MeBr also has a bacterial soil sink (5) that represents about 20% of the estimated total removal from the troposphere, and it is likely that a soil sink of similar magnitude exists for MeCl (6). Hence, if an isotopic fractionation is associated with the soil sink, it will influence the isotopic compositions of MeBr and MeCl in the lower atmosphere (7). The δ 13 C value of industrially produced MeBr ranges between −43.5‰ and −66.4‰ (7), but δ 13 C values of tropospheric MeBr and natural sources are not yet known. The δ 13 C of atmospheric MeCl has been measured from −22‰ to −45‰ (8, 9). If carbon isotope ratios are to be used to constrain the budgets of these methyl halides, it is essential to determine the extent of carbon isotope fractionation that occurs during biological degradation of these compounds. Methylotrophic bacteria use C 1 compounds, which are simple organic molecules that contain no carbon–carbon bonds. Strains IMB-1, CC495, and MB2 are as-yet-unnamed facultative methylotrophs isolated from agricultural soil, woodland leaf litter, and coastal seawater, respectively (10–13), environments where methyl halides are produced. They are members of the α subgroup of the Proteobacteria . On the basis of 16S rRNA gene sequences, strains IMB-1 and CC495 show some phylogenetic alignment with the genus Rhizobium (10, 11) and are very closely related to the new genus Pseudoaminobacte r (I. McDonald, personal communication). Strain MB2 aligns within the Ruegeria clade [J. K. Schaefer, K. D. Goodwin, I. R. McDonald, J. C. Murrell and R.S.O., unpublished work]. All of these aerobic bacteria are methylotrophs in that they can grow by using MeBr or MeCl as their sole carbon source, but they do not metabolize methane. They oxidize MeBr, MeCl, and methyl iodide (MeI) to CO 2 . Soil bacteria are known to consume MeBr at the ambient tropospheric mixing ratio of around 10 parts per trillion by volume (5). Preliminary experiments with strain IMB-1 indicate that it can oxidize MeBr at these mixing ratios** and is therefore likely to be characteristic of bacteria associated with MeBr uptake by soils. We examined δ 13 C of MeCl, MeBr, and MeI during oxidation by whole-cell suspensions of IMB-1 and CC495 and also the change in δ 13 C values of the three methyl halides during oxidation by the marine strain MB2. In addition, we measured the fractionation of carbon isotopes during formation of methane thiol (MeSH) from MeCl by the purified cobalamin-dependent enzyme, halomethane:bisulphide/halide ion methyltransferase (11) from CC495, to determine whether this initial step in MeCl degradation could account for the observed fractionation by whole cells. Finally, we determined the fractionation associated with the degradation of MeBr during field studies with agricultural soil by monitoring MeBr concentration and δ 13 C of MeBr in the headspace of flux chambers under fumigation conditions.

Proceedings of the National Academy of Sciences of↗

Handbook of acute toxicity of chemicals to fish and aquatic invertebrates : summaries of toxicity tests conducted at Columbia National Fisheries Research Laboratory, 1965-78

Acute toxicity is a major subject of research at Columbia National Fisheries Research Laboratory for evaluating the impact of toxic chemicals on fishery resources. The Laboratory has played a leading role in developing research technology for toxicity testing and data interpretation. In 1965-78, more than 400 chemicals were tested against a variety of invertebrates and fish species representative of both cold- and warm-water climates. The use of acute toxicity tests for assessing the potential hazard of chemical contaminants to aquatic organisms is well documented (Boyd 1957; Henderson et al. 1960; Sanders and Cope 1966; Macek and McAllister 1970). Static acute toxicity tests provide rapid and (within limits) reproducible concentration-response curves for estimating toxic effects of chemicals on aquatic organisms. These tests provide a database for determining relative toxicity of a large number of chemicals to a variety of species and for estimating acute effects of chemical spills on natural aquatic systems; they also assist in determining priority and design of additional toxicity studies. Acute toxicity tests usually provide estimates of the exposure concentration causing 50% mortality (LC50) to test organisms during a specified period of time. For certain invertebrates, the effective concentration is based on immobilization, or some other identifiable endpoint, rather than on lethality. The application of the LC50 has gained acceptance among toxicologists and is generally the most highly rated test for assessing potential adverse effects of chemical contaminants to aquatic life (Brungs and Mount 1978; American Institute for Biological Sciences 1978 a ). The literature contains numerous papers dealing with the acute toxicity of chemicals to freshwater organisms. However, there is a tremendous need for a concise compendium of toxicity data covering a large variety of chemicals and test species. This Handbook is a compilation of a large volume of acute toxicity data from the Columbia Laboratory and its field laboratories. It presents definitive acute toxicity data on 271 chemicals tested against a variety of freshwater invertebrates and fishes. The chemicals represent all major groups of pesticides, as well as numerous industrial chemicals. This compilation should serve as a useful database for the many agencies and organizations dealing with research and management programs concerned with the impact of chemicals on aquatic resources. The Columbia Laboratory has played a major role in developing currently used standard methodology for static acute toxicity testing. The use of standardized methodology greatly reduces variation in results. The data presented here have been carefully scrutinized to eliminate tests that failed to follow acceptable procedures. Handling of test organisms and procedures for static toxicity tests followed those described by Lennon and Walker (1964) and Macek and McAllister (1970), and conform well with those recommended by Brauhn and Schoettger (1975) and the Committee on Methods for Toxicity Tests with Aquatic Organisms (1975). The species of fish and invertebrates that were tested are listed in phylogenetic order in Tables 1 and 2. Fish were obtained from Federal and State hatcheries as either eggs or fry. Original stocks of invertebrates were collected and cultured from wild populations with no known source of contamination; these populations were replenished regularly. The invertebrates were cultured in the Laboratory by methods similar to those described by Sanders and Cope (1966). Test chemicals usually consisted of technical or analytical grade samples of known purity. Formulations of the chemicals were also tested when available. When purity of test chemicals was known, all calculated concentrations were based on percent active ingredients. Stock solutions were prepared immediately before each test, with commercial grade acetone as the carrier solvent. Occasionally, ethanol or dimethyl-formamide was substituted. Solvent concentrations did not exceed 0.5 mL/L in final dilution water. Test water (dilution water) was reconstituted from deionized water of at least 10 6 ohms resistivity by the addition of appropriate reagent grade chemicals (Marking 1969). Water was buffered to maintain a pH of 7.2 to 7.5, an alkalinity of 30 to 35 mg/L, and a hardness of 40 to 50 mg/L as CaCO 3 . Test water was mixed thoroughly and aerated before transfer into test chambers. Fish were acclimated to dilution water by gradually changing the water in acclimated tanks from 100% well water to 100% reconstituted water over a 1- to 3-day period at the desired testing temperature. Invertebrates were acclimated from well water to dilution water over a 4- to 6-h period. Toxicity tests were conducted under static conditions without aeration, and the organisms were not fed during acclimation or testing. Temperature of test solutions was maintained within ± 1°C of that required for a given test. Toxicity tests with fish were conducted in 18.9-liter (5-gal) wide-mouthed jars containing 15 liters of test solution. Fingerling fish weighing 0.2 to 1.5 g were tested at each concentration. Caution was taken not to exceed 0.8 g of test organisms per liter of solution. Duplicate test chambers were used to accommodate larger fish. Test chambers varied in size for invertebrates, depending on the species used; volume of test solution ranged from 0.25 to 4 liters. At least 10 organisms were exposed to each concentration for all definitive tests. At least six concentrations were used per toxicity test. The tests began upon initial exposure to the toxicant and continued for 96 h. Immobilization tests with invertebrates were conducted for only 48 h. The number of dead or affected organisms in each test chamber were recorded and the dead organisms were removed every 24 h; general observations on the condition of test organisms were also recorded at these times. Toxicity data were analyzed by a statistical method described by Litchfield and Wilcoxon (1949) to determine LC50 (theoretical estimate of the concentration lethal to 50% of the test animals) and 95% confidence intervals. This method is recommended by the American Public Health Association (1971) and by Sprague (1969) for determining median lethal concentrations. The procedure is easily modified for computing a single LC50 when replicate tests are performed.

Resource Publication↗

Two-dimensional streamflow simulations of the Jordan River, Midvale and West Jordan, Utah

The Jordan River in Midvale and West Jordan, Utah, flows adjacent to two U.S. Environmental Protection Agency Superfund sites: Midvale Slag and Sharon Steel. At both sites, geotechnical caps extend to the east bank of the river. The final remediation tasks for these sites included the replacement of a historic sheet-pile dam and the stabilization of the river banks adjacent to the Superfund sites. To assist with these tasks, two hydraulic modeling codes contained in the U.S. Geological Survey (USGS) Multi-Dimensional Surface-Water Modeling System (MD_SWMS), System for Transport and River Modeling (SToRM) and Flow and Sediment Transport and Morphological Evolution of Channels (FaSTMECH), were used to provide predicted water-surface elevations, velocities, and boundary shear-stress values throughout the study reach of the Jordan River. A SToRM model of a 0.7 mile subreach containing the sheet-pile dam was used to compare water-surface elevations and velocities associated with the sheet-pile dam and a proposed replacement structure. Maps showing water-surface elevation and velocity differences computed from simulations of the historic sheet-pile dam and the proposed replacement structure topographies for streamflows of 500 and 1,000 cubic feet per second (ft 3 /s) were created. These difference maps indicated that the velocities associated with the proposed replacement structure topographies were less than or equal to those associated with the historic sheet-pile dam. Similarly, water-surface elevations associated with the proposed replacement structure topographies were all either greater than or equal to water-surface elevations associated with the sheet-pile dam. A FaSTMECH model was developed for the 2.5-mile study reach to aid engineers in bank stabilization designs. Predicted water-surface elevations, velocities and shear-stress values were mapped on an aerial photograph of the study reach to place these parameters in a spatial context. Profile plots of predicted cross-stream average water-surface elevations and cross-stream maximum and average velocities showed how these parameters change along the study reach for two simulated discharges of 1,040 ft 3 /s and 2,790 ft 3 /s. The profile plots for the simulated streamflow of 1,040 ft 3 /s show that the highest velocities are associated with the constructed sheet-pile replacement structure. Results for the simulated streamflow of 2,790 ft 3 /s indicate that the geometry of the 7800 South Bridge causes more backwater and higher velocities than the constructed sheet-pile replacement structure.

Utah↗

Spatial distribution, temporal variability, and chemistry of the salt wedge in the lower Charles River, Massachusetts, June 1998 to July 1999

The Charles River is of great recreational and ecological value to the Boston metropolitan region and the Commonwealth of Massachusetts. It is also the focus of the U.S. Environmental Protection Agency (USEPA) Region I, Clean Charles 2005 Task Force. The main goal of the Task Force is to make the Charles River 'fishable and swimmable' by the year 2005. Achieving 'fishable and swimmable' conditions will require continued progress in addressing a range of environmental conditions now degrading water quality, including the infiltration of saltwater from Boston Harbor into the freshwater Charles River.To better understand the pattern of saltwater intrusion, the U.S. Geological Survey (USGS), in cooperation with the U.S. Environmental Protection Agency (USEPA), Massachusetts Department of Environmental Management (MADEM), and New England Interstate Water Pollution Control Commission (NEIWPCC), collected data on the spatial distribution, temporal variability, and chemistry of the saltwater that entered the lower Charles River from June 1998 to July 1999. The purpose of this investigation is to extend and complement a regional-scale study of Charles River water quality conducted in 1996 (T. Faber, U.S. Environmental Protection Agency, written commun., 1997), and the ongoing water monitoring activities of the Massachusetts Water Resources Authority (MWRA) and the Charles River Watershed Association (CRWA). The data collected by this investigation supports the Clean Charles 2005 Task Force by providing detailed information concerning a major factor limiting 'fishable and swimmable' conditions in the lower Charles River. Finally, the study will be used to assist current planning efforts of the Metropolitan District Commission (MDC) to restore the historic parklands of the lower Charles River.The 'Basin' is the local term for the reach of the Charles River that begins at the Watertown Dam in Watertown, Mass., and extends about 8 mi through suburban and urban areas to Boston Harbor. Discharge to the harbor is controlled by the 'new' Charles River Dam in Boston (fig. 1). The Basin was created by construction of the 'old' Charles River Dam in 1908 to solve Boston's sanitary problems. Prior to the building of the old Charles River Dam, the lower Charles River was a tidal estuary in which the water levels rose and fell twice daily with the tidal cycle. Low tide would expose untreated sewage that was discharged directly into the river. Exposed sewage created noxious odors and served as a breeding ground for mosquitoes that caused sporadic epidemics of malaria and yellow fever (Jobin, 1998). Damming of the river interrupted the normal tidal cycle and flooded the estuary by creating a freshwater pool (the Basin) that had a constant water elevation of about 0.8 meters (m) above mean sea level. Flooding of the estuary initially improved sanitary conditions and the Basin became a source of enjoyment for the local population and the focus of a large waterfront park in Boston and Cambridge (Jobin, 1998).Although the infiltration of saltwater from the harbor into the Basin was anticipated when the old Charles River Dam was built, neither the magnitude nor the consequences of the infiltration was considered. By 1975, the Metropolitan District Commission (MDC) determined that harbor water covered about 80 percent of the river bottom within the Basin and composed about 50 percent of its depth. The MDC also concluded that fish kills and odors in the spring of 1975 were likely the result of the sulfide-rich saltwater mixing with the overlying freshwater (Metropolitan District Commission, 1975).Saltwater from Boston Harbor that enters the Basin is known as the 'salt wedge' because of the shape it assumes as it moves upstream. Freshwater discharge from upstream pushes against the intruding harbor water until the density differences cause stratification to occur; the freshwater then overrides the denser harbor water (Fischer and others, 1979). The depth from

Water-Resources Investigations Report↗

Reconnaissance map of the Cenozoic geology in the Carlin basin area, Elko and Eureka counties, Nevada

The middle Miocene Carlin sedimentary basin encompasses a large area between the Adobe Range to the east, the Piñon Range to the south, the southern Independence Mountains and Marys Mountain to the west, and Swales Mountain to the north. The town of Carlin is in the southern part of the basin. The geologic map includes detailed to more reconnaissance mapping of Cenozoic units in the main part of the basin, including different facies of the middle Miocene Humboldt Formation. The mapping was part of a broader study of the Miocene and younger paleogeographic evolution of the region. Earlier work obtained numerous 40 Ar/ 39 Ar and tephra correlation dates on sedimentary and volcanic units in the basin. The basin connected to the east into the Elko sedimentary basin and to the northwest and northeast into similar smaller basins between present-day mountain ranges. Early sediments, largely fluvial, began to accumulate in the lowlands between the surrounding ranges at about 16.5 Ma. The sediments were derived from Paleozoic sedimentary and middle Tertiary volcanic units in the nearby highlands, and flow patterns indicate a general southward flow towards present-day Pine Valley. The eruption of the Palisade Canyon–Marys Mountain rhyolite flows at the southwest end of the basin at 15.3 Ma blocked the southward flow, and a lake began to form in the basin. As the lake grew in extent, sedimentary units around the fringes of the lake included a mixture of inflowing fluvial sediments mixed with the pyroclastic-fall, ash-rich sediments deposited in the lake. The lake margin expanded, and stratigraphic sections record the progressive transition from fluvial to mixed fluvial and lacustrine, and finally to entirely lacustrine. The volcanic rock dam was breached at about 14.7 Ma, the lake drained, and fluvial sediments blanketed the entire basin for an unknown period of time after that. Sedimentation progressively buried existing highlands and bridged gaps between adjacent basins. For example, the Carlin and Elko basins connected across the southern Adobe Range. Normal faulting produced numerous, mostly north- to north-northeast-striking faults that cut the sedimentary units and surrounding highlands largely after sedimentation ceased. The largest fault formed in the eastern third of the basin and tilted all of the sedimentary units in the western two-thirds of the basin, as well as the eastern part of Marys Mountain, to the east. Some offset took place during sedimentation. Many other normal faults of smaller extent and offset cut the sedimentary units. The integration of streams draining the Elko and Carlin basins began after about 9.8 Ma. The streams, which together comprised the early stages of the Humboldt River, flowed regionally southwestward beyond the Carlin basin. As many as thirteen, downward-stepping strath terraces in the Carlin basin record the progressive downcutting into and removal of the middle Miocene sediments. Gravel deposits form a thin veneer on some of the higher terraces. Clasts in those gravel deposits, as well as the overall terrace pattern, indicate southward drainage towards the Humboldt River. The erosion gradually re-exposed the flanks of the surrounding highlands. A brief pause in downcutting allowed the formation of a small lake in the Hemphillian (late Miocene), represented by lacustrine units northwest of Carlin. The sedimentary rocks of the Carlin basin conceal a large segment of the world-class, late Eocene Carlin gold trend, which extends from the southern Independence Mountains south into the Piñon Range. Sedimentation largely buried the Gold Quarry and Mike gold deposits in the northwestern part of the basin. Later faulting and erosion re-exposed the Gold Quarry deposit, but the Mike deposit remains buried. The basin’s sedimentary units conceal potential Paleozoic host rocks, and the sedimentary facies and post-sedimentation faults shown on the map may help guide interpretations of geophysical and other exploration data in the Carlin basin. The current map publication was supported by the USGS National Cooperative Geologic Mapping Program under STATEMAP award number G19AC00383.

Nevada↗

Seismic refraction profile, Kingdom of Saudi Arabia: field operations, instrumentation, and initial results

In February 1978 a seismic deep-refraction profile was recorded by the USGS along a 1000-km line across the Arabian Shield in western Saudi Arabia. The line begins in Paleozoic and Mesozoic cover rocks near Riyadh on the Arabian Platform, leads southwesterly across three major Precambrian tectonic provinces, traverses Cenozoic rocks of the coastal plain near Jizan (Tihamat Asir), and terminates at the outer edge of the Farasan Bank in the southern Red Sea. More than 500 surveyed recording sites were occupied, including 19 in the Farasan Islands. Six shot points were used--five on land, with charges placed mostly below water table in drill holes, and one at sea, with charges placed on the sea floor and fired from a ship. The total charge consumed was slightly in excess of 61 metric tons in 21 discrete firings. Seismic energy was recorded by means of a set of 100 newly developed portable seismic stations. Each station consists of a standard 2-Hz vertical geophone coupled to a self-contained analog recording instrument equipped with a magnetic-tape cassette. The stations were deployed in groups of 20 by five observer teams, each generally consisting of two scientist-technicians and a surveyor-guide. On the day prior to deployment, the instruments were calibrated and programmed for automatic operation by means of a specially designed device called a hand-held tester. At each of ten pre-selected recording time windows on a designated firing day, the instruments were programmed to turn on, stabilize, record internal calibration signals, record the seismic signals at three levels of amplification, and then deactivate. After the final window in the firing sequence, all instruments were retrieved and their data tapes removed for processing. A specially designed, field tape- dubbing system was utilized at shot point camps to organize and edit data recorded on the cassette tapes. The main functions of this system are to concatenate all data from each shot on any given day onto a single shot tape, and to provide hard copy for monitoring recorder performance so that any problems can be corrected prior to the next deployment. Composite digital record sections were produced from the dubbed tapes for each shot point by a portable processing and plotting system. The heart of this system is a DEC PDP 11VO3 computer, which controls a cassette playback unit identical to those used in the recorders and dubbers, a set of discriminators, a time-code translator, a digitizer, and a digital plotter. The system was used to maintain various informational data sets and to produce tabulations and listings of various sorts during the field operations, in addition to its main task of producing digital record sections. Two master clocks, both set to time signals broadcast by the British Broadcasting Corporation, provided absolute time for the recording operations. One was located on the ship and the other was stationed at a base camp on the mainland. The land-based master clock was used to set three additional master clocks located at the other active shot points a few days in advance of each firing, and these clocks were then used to set the internal clocks in the portable seismic stations via the hand-held tester. A master clock signal was also linked to the firing system at each shot point for determination of the absolute shot instant. It is possible to construct a generalized crustal model from examination of the six shot point composite record sections obtained in the field. Such a model rests upon a number of simplifying assumptions and will almost certainly be modified at a later stage of interpretation. The main assumptions are that the crust consists of two homogeneous isotropic layers having no velocity inversion,, that the Mohorovicic discontinuity is sharp, and that effects of surface inhomogeneities and elevation changes can be ignored. The main characteristics of the tentative model are the following: (1) The thickness of th

Open-File Report↗

A spatially explicit suspended-sediment load model for western Oregon

We calibrated the watershed model SPARROW (Spatially Referenced Regressions on Watershed attributes) to give estimates of suspended-sediment loads for western Oregon and parts of northwestern California. Estimates of suspended-sediment loads were derived from a nonlinear least squares regression that related explanatory variables representing landscape and transport conditions to measured suspended-sediment loads at 68 measurement stations. The model gives estimates of model coefficients and their uncertainty within a spatial framework defined by the National Hydrography Dataset Plus hydrologic network. The resulting model explained 64 percent of the variability in suspended-sediment yield and had a root mean squared error value of 0.737. The predictor variables selected for the final model were (1) generalized lithologic province, (2) mean annual precipitation, and (3) burned area (by recent wildfire). Other landscape characteristics also were considered, but they were not significant predictors of sediment transport, were strongly correlated with another predictor variable, or were not as significant as the predictors selected for the final model. The northern Oregon coastal drainages had the highest predicted suspended sediment yields (median yield 475 kilograms per hectare per year) and the Klamath River Basin had the lowest (median yield 53 kilograms per hectare per year). Quaternary deposits were, on average, the largest contributor to incremental suspended-sediment yield even though this lithologic province only makes up 17 percent of the modeling domain. Coast Range sedimentary rocks and Coast Range volcanic rocks had high suspended-sediment yields whereas, in addition to the Klamath terrane, the Western Cascade and High Cascade lithologic provinces had low suspended-sediment yields. Precipitation and the area affected by recent wildfire both positively correlated with suspended-sediment load. Suspended-sediment transport rates predicted by this SPARROW model are less than historical (1956–73) and long‑term (thousands of years) geological rates. This difference likely results, in part, from biases in the data underlying the SPARROW model, probably resulting in predicted suspended-sediment estimates that underestimate actual transport rates. However, the differences also likely owe to natural and human-caused variation in suspended-sediment yields as they respond to changes in climate, vegetation, fire frequency, and land use. In particular, decreases in mean annual suspended-sediment yields within the Umpqua River Basin since 1956–73 may owe to less intense forest harvest, passage of the Oregon Forest Practices Act of 1971, and increased emphasis in habitat protection in recent decades. Such sensitivity may have implications for the spatial and temporal distributions of aquatic and riparian habitats. Knowledge of the regionally important patterns and factors in suspended-sediment sources and transport could support broad-scale, water-quality management objectives and priorities. Because of biases and limitations of this model, however, these results are most applicable for general comparisons and for broad areas such as large watersheds. For example, despite having similar area, precipitation, and land-use, the Umpqua River Basin generates 68 percent more suspended sediment than the Rogue River Basin, chiefly because of the large area of Coast Range sedimentary province in the Umpqua River Basin. By contrast, the Rogue River Basin contains a much larger area of Klamath terrane rocks, which produce significantly less suspended load, although recent fire disturbance (in 2002) has apparently elevated suspended sediment yields in the tributary Illinois River watershed. Fine-scaled analysis, however, will require more intensive, locally focused measurements.

Scientific Investigations Report↗