Geology ReportsSearch

SEARCH · Geology Reports

Results for “Project Completion 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 613 records · Page 34Linked to original sources

Preliminary Earthquake Hazard Map of Afghanistan

Introduction Earthquakes represent a serious threat to the people and institutions of Afghanistan. As part of a United States Agency for International Development (USAID) effort to assess the resource potential and seismic hazards of Afghanistan, the Seismic Hazard Mapping group of the United States Geological Survey (USGS) has prepared a series of probabilistic seismic hazard maps that help quantify the expected frequency and strength of ground shaking nationwide. To construct the maps, we do a complete hazard analysis for each of ~35,000 sites in the study area. We use a probabilistic methodology that accounts for all potential seismic sources and their rates of earthquake activity, and we incorporate modeling uncertainty by using logic trees for source and ground-motion parameters. See the Appendix for an explanation of probabilistic seismic hazard analysis and discussion of seismic risk. Afghanistan occupies a southward-projecting, relatively stable promontory of the Eurasian tectonic plate (Ambraseys and Bilham, 2003; Wheeler and others, 2005). Active plate boundaries, however, surround Afghanistan on the west, south, and east. To the west, the Arabian plate moves northward relative to Eurasia at about 3 cm/yr. The active plate boundary trends northwestward through the Zagros region of southwestern Iran. Deformation is accommodated throughout the territory of Iran; major structures include several north-south-trending, right-lateral strike-slip fault systems in the east and, farther to the north, a series of east-west-trending reverse- and strike-slip faults. This deformation apparently does not cross the border into relatively stable western Afghanistan. In the east, the Indian plate moves northward relative to Eurasia at a rate of about 4 cm/yr. A broad, transpressional plate-boundary zone extends into eastern Afghanistan, trending southwestward from the Hindu Kush in northeast Afghanistan, through Kabul, and along the Afghanistan-Pakistan border. Deformation here is expressed as a belt of major, north-northeast-trending, left-lateral strike-slip faults and abundant seismicity. The seismicity intensifies farther to the northeast and includes a prominent zone of deep earthquakes associated with northward subduction of the Indian plate beneath Eurasia that extends beneath the Hindu Kush and Pamirs Mountains. Production of the seismic hazard maps is challenging because the geological and seismological data required to produce a seismic hazard model are limited. The data that are available for this project include historical seismicity and poorly constrained slip rates on only a few of the many active faults in the country. Much of the hazard is derived from a new catalog of historical earthquakes: from 1964 to the present, with magnitude equal to or greater than about 4.5, and with depth between 0 and 250 kilometers. We also include four specific faults in the model: the Chaman fault with an assigned slip rate of 10 mm/yr, the Central Badakhshan fault with an assigned slip rate of 12 mm/yr, the Darvaz fault with an assigned slip rate of 7 mm/yr, and the Hari Rud fault with an assigned slip rate of 2 mm/yr. For these faults and for shallow seismicity less than 50 km deep, we incorporate published ground-motion estimates from tectonically active regions of western North America, Europe, and the Middle East. Ground-motion estimates for deeper seismicity are derived from data in subduction environments. We apply estimates derived for tectonic regions where subduction is the main tectonic process for intermediate-depth seismicity between 50- and 250-km depth. Within the framework of these limitations, we have developed a preliminary probabilistic seismic-hazard assessment of Afghanistan, the type of analysis that underpins the seismic components of modern building codes in the United States. The assessment includes maps of estimated peak ground-acceleration (PGA), 0.2-second spectral acceleration (SA), and 1.0-secon

Open-File Report

Status and understanding of groundwater quality in the northern Sierra Nevada foothills domestic-supply aquifer study units, 2015–17—California GAMA Priority Basin Project

Groundwater quality in the northern Sierra Nevada foothills region of California was investigated as part of California State Water Resources Control Board (SWRCB) Groundwater Ambient Monitoring Assessment Priority Basin Project (GAMA-PBP). The region was divided into two study units: the Yuba-Bear watersheds (YBW) study unit and the American-Cosumnes-Mokelumne watersheds (ACMW) study unit. The GAMA-PBP made a spatially unbiased assessment of aquifer systems used for domestic drinking-water supply in the study region, which are predominantly composed of fractured, hard-rock aquifers of varying lithology. These assessments characterized the quality of raw groundwater to evaluate ambient conditions in the domestic-supply aquifer and not the quality of treated drinking water. The study included three components: (1) a status assessment , which characterized the quality of groundwater resources used for domestic drinking-water supply in the YBW and ACMW study units; (2) an understanding assessment , which evaluated natural and anthropogenic explanatory factors that could potentially affect groundwater quality in the study region; and (3) a comparative assessment between the groundwater resources used for domestic and public drinking-water supply in the study region. The status assessment was based on data collected by the GAMA-PBP from 74 sites in the YBW study unit during 2015–16 and 67 sites in the ACMW study unit from 2016 to 2017. To contextualize water-quality results, concentrations of water-quality constituents in ambient groundwater were compared to regulatory and non-regulatory benchmarks typically used by the State of California and Federal agencies as health-based or aesthetic standards for public drinking water. The status assessment used a grid-based method to estimate proportions of groundwater resources with concentrations approaching or exceeding benchmark thresholds. This method provides spatially unbiased results and allows inter-comparability with similar groundwater-quality assessments. Inorganic constituents with health-based benchmarks were present at high relative concentration (RC), meaning they exceeded the benchmark threshold, in 5.4 and 10 percent of domestic-supply aquifer systems in the YBW and ACMW study units, respectively. Inorganic constituents with aesthetic-based benchmarks were detected at high-RCs in 20 and 28 percent of the YBW and ACMW study units, respectively. The inorganic constituents present at high RC were arsenic, barium, boron, molybdenum, strontium, nitrate, adjusted gross-alpha particle activity, chloride, total dissolved solids, specific conductance, iron, manganese, and hardness. Groundwater samples were tested for presence or absence of three microbial indicators (total coliform, Escherichia coli , and Enterococci ). At least one microbial indicator was present in 26 and 28 percent of the YBW and ACMW study units, respectively. At least one organic constituent was detected in 30 and 42 percent of the YBW and ACMW study units, respectively. Organic constituents were not present at high RC, but tetrachloroethene (PCE), trichloroethene (TCE), and toluene were detected in the YBW study unit at moderate RC (between the benchmark concentration and one-tenth of the benchmark concentration). Methyl tert -butyl ether (MTBE) and chloroform were present at low RC (less than one-tenth of the benchmark concentration) in the YBW and ACMW study units with detection frequencies greater than 10 percent. Perchlorate, a constituent of special interest, was detected in 31 and 41 percent of the YBW and ACMW study units, respectively, at either low or moderate RCs. Relations among select water-quality constituents and potential explanatory factors were evaluated using statistical and graphical approaches. Nitrate, microbial indicators, and perchlorate were all correlated to elevation-dependent variables relating to climate, land use, and recharge condition. Isotopic and dissolved noble-gas tracers indicated these water-quality constituents are associated with recharge conditions associated with irrigation during the summer dry-season, which is common in areas of rural-residential or agricultural land uses. Higher concentrations of iron and manganese were primarily associated with anoxic groundwater in aquifers of metasedimentary lithology. Increased hardness was primarily associated with anoxic groundwater in aquifers of mafic-ultramafic or metavolcanics lithologies at lower elevations in the study region in the Melones fault zone. Chloroform and MTBE were associated with shallow groundwater (wells depths less than 130 m) under oxic and anoxic redox conditions, respectively. The comparative assessment evaluated differences between the aquifer systems used for domestic- and public-supply in study region based on (1) well-construction characteristics, and (2) water quality. Analysis of over 60,000 well-completion reports in the study region showed that although domestic-supply wells span the deepest depth zones in regional aquifers, median depths for public-supply wells were significantly greater than those of domestic-supply wells in both study units. Water-quality data from more than 300 public-supply wells in the study region were assessed using a spatially weighted method for calculation aquifer-scale proportions and compared with the domestic-supply assessment results. Detections of inorganic constituents at high RC and detection frequencies for organic constituents were generally similar between the domestic- and public-supply aquifer systems in both study units, with a few notable exceptions in the ACMW study unit: nitrate was greater for the public- compared to domestic-supply aquifer system and both manganese, hardness, and MTBE were greater in the domestic- compared to public-supply aquifer system. These differences are likely related to contrasting land uses, aquifer lithologies, landscape positions, and depths characterizing domestic- and public-supply wells in the ACMW study unit. Overall, fewer samples from domestic-supply wells in the northern Sierra Nevada foothills exceeded health-based benchmarks compared to aesthetic-based benchmarks for groundwater quality. Exceedences of health-based benchmarks were primarily caused by nitrate and coliform bacteria, which were associated with recharge from diverted surface water used primarily for irrigation. Exceedences of aesthetic-based benchmarks were primarily caused by iron, managanese, and hardness, which were associated with geologic factors. Regional irrigation practices and aquifer lithology can affect groundwater quality in fractured-rock aquifers in the northern Sierra Nevada foothills used for domestic drinking-water supply.

California

The Pajarito Plateau: A bibliography

This bibliography is the result of two initially independent projects. As the consulting archaeologist at Los Alamos National Laboratory (LANL), Charlie R. Steen collected entries at the suggestion of the staff of the Environmental Surveillance Group of the Health, Safety, and Environmental Division, HSE-8. The primary purpose was to aid the staff in evaluating cultural resources on LANL lands. In addition to works that related to the archaeology and history of the area, Steen included notations of a few books and articles in other fields such as geology and natural history. It was hoped that they also would be of value to other organizations and to students of past human activities on the Pajarito Plateau. At the same time, the National Park Service (NPS) was planning a major survey of Bandelier National Monument (BNM). As part of this plan, the author was asked to prepare a background document that described research previously carried out in the area, including an annotated bibliography. Although the survey would be limited to the park boundaries, the larger Pajarito Plateau is a more logical study area from physiographic, environmental, and cultural perspectives; hence the focus was on this larger region. Mathien (1986) also included some references to natural resources studies, particularly those initiated by NPS within Bandelier National Monument. Both bibliographies were made available to Colleen Olinger and Beverly Larson of the Health and Environmental Services Group at Los Alamos. They realized that while neither was complete, each included entries missing from the other. Larson suggested the two bibliographies be combined. (At this time, Craig Allen was studying the landscape of the Jemez Mountains [Allen 1984c, 1989]. His investigations included much detailed information on natural resource studies and were added in 1991 and 1992.) To limit the scope of their work, Steen and Mathien had chosen their parameter: the Pajarito Plateau. Geographically, the Pajarito Plateau is described as the high tableland that lies between the Jemez Mountains on the west and the Rio Grande on the east. From north to south, it extends from the Chama Valley to La Cañada de Cochiti (Hewett 1906:14)(Figure 1). Because human activity rarely stops at such definite boundaries, major ethnographic studies of Tewa (San Ildefonso and Santa Clara) and Keres (Cochiti) linguistic groups are included. (Even though most of the historic pueblos occupied by the Tewa and Keres are not located on the Pajarito Plateau, oral traditions and archaeological data suggest that these groups once occupied sites on the plateau.) Towa studies are not included because Steen believed Towa ancestors were not involved in major cultural developments of the Pajarito Plateau. In addition, a bibliography of the Jemez area (home of Towa people) has been prepared by Michael Elliott (1982) and included with his nomination of large Pueblo sites near Jemez Springs to the National Register of Historic Places that is on file at the Museum of New Mexico, Laboratory of Anthropology, in Santa Fe. Both Steen and Mathien included references to geographically and historically related material that does not focus on the Pajarito Plateau but, nonetheless, is important to understanding the area's archaeology and physical environment, for example, lithic resources available from Cerro Pedernal or in the Sangre de Cristo Mountains.

New Mexico

Water and sediment quality of the Lake Andes and Choteau Creek basins, South Dakota, 1983-2000

The Bureau of Reclamation has proposed construction of the Lake Andes/Wagner Irrigation Demonstration Project to investigate environmental effects of irrigation of glacial till soils substantially derived from marine shales. During 1983-2000, the U.S. Geological Survey collected hydrologic, water-quality, and sediment data in the Lake Andes and Choteau Creek Basins, and on the Missouri River upstream and downstream from Choteau Creek, to provide baseline information in support of the proposed demonstration project. Lake Andes has a drainage area of about 230 mi2 (square miles). Tributaries to Lake Andes are ephemeral. Water-level fluctuations in Lake Andes can be large, and the lake has been completely dry on several occasions. The outlet aqueduct from Lake Andes feeds into Garden Creek, which enters Lake Francis Case just upstream from Fort Randall Dam on the Missouri River. For Lake Andes tributary stations, calcium, magnesium, and sodium are approximately codominant among the cations, and sulfate is the dominant anion. Dissolved-solids concentrations typically range from about 1,000 mg/L (milligrams per liter) to about 1,700 mg/L. Major-ion concentrations for Lake Andes tend to be higher than the tributaries and generally increase downstream in Lake Andes. Proportions of major ions are similar among the different lake units (with the exception of Owens Bay), with calcium, magnesium, and sodium being approximately codominant among cations, and sulfate being the dominant anion. Owens Bay is characterized by a calcium sulfate water type. Dissolved-solids concentrations for Lake Andes typically range from about 1,400 to 2,000 mg/L. Whole-water nitrogen and phosphorus concentrations are similar among the Lake Andes tributaries, with median whole-water nitrogen concentrations ranging from about 1.6 to 2.4 mg/L, and median whole-water phosphorus concentrations ranging from about 0.5 to 0.7 mg/L. Whole-water nitrogen concentrations in Lake Andes are similar among the different units, with medians that range from about 2.4 to 4.0 mg/L. Median whole-water phosphorus concentrations for the different Lake Andes units range from 0.2 to 0.5 mg/L, and decrease downstream through Lake Andes. Median selenium concentrations are substantially lower for Andes Creek (3 ?g/L (micrograms per liter)) than for the other tributary stations (34, 18, and 7 ?g/L). Median selenium concentrations for the lake stations (ranging from less than 1 to 2 ?g/L) are substantially lower than tributary stations. The pesticides 2,4-D and atrazine were the most commonly detected pesticides in Lake Andes. Median concentrations for 2,4-D for Lake Andes range from 0.07 to 0.11 ?g/L; the median concentration for Owens Bay is 0.04 ?g/L. Median concentrations for atrazine for Lake Andes range from 0.2 to 0.4 ?g/L; the median concentration for Owens Bay is less than 0.1 ?g/L. Concentrations of both 2,4-D and atrazine are largest for the most upstream part of Lake Andes that is most influenced by tributary inflow. Median suspended-sediment concentrations for Lake Andes tributaries range from 22 to 56 mg/L. Most of the suspended sediment transported in the Lake Andes tributaries consists of particles less than 63 ?m (micrometers) in diameter. Concentrations of most constituents in bottom sediments generally had similar ranges and medians for the Lake Andes tributaries. However, Andes Creek generally had lower concentrations of several metals. For Lake Andes, medians and ranges for most constituents generally were similar among the different units. However, selenium concentrations tended to be higher in the upstream part of the lake, and generally decreased downstream. Results of vertical sediment cores collected from a single site in the South Unit of Lake Andes in October 2000 indicate that selenium loading to Lake Andes increased during the period 1952 through 2000. Choteau Creek has a drainage area of 619 mi2. In the upstream part of the basin, Chotea

Water-Resources Investigations Report

Using environmental DNA (eDNA) to assess the presence of cavefish and cave crayfish populations in caves of the Ozark Highlands

Many cavefishes and cave crayfishes are considered of conservation concern; however, sampling these species is inherently difficult given their occupied environments. The goal of our project was to verify the presence of select karst organisms while developing the foundation for sampling approaches that might be useful to conservation and management agencies. Our project objectives were to develop assays to amplify deoxyribonucleic acid (DNA) from several species of Ozark cavefishes and cave crayfishes and complete an initial surveillance of locations across the Ozark Highlands using environmental DNA (eDNA). Using DNA either provided by agency cooperators or that we extracted from tissue samples, we PCR amplified and then sequenced the Cytochrome Oxidase 1 (CO1) gene for cave crayfishes and the NADH Dehydrogenase Subunit 2 (ND2) gene for cavefishes. We developed species-specific primers and probes for five cave crayfishes and two cavefishes. From February 2017 to May 2017, we sampled 1–5 sampling units from 42 caves, wells, and springs (i.e., sites) using eDNA and traditional visual surveys. We measured physicochemical parameters at each sampling unit to estimate detection probability associated with both techniques. We also calculated two occupancy covariates for each site using geospatial data. We successfully amplified Troglichthys rosae DNA from the environment and detected DNA representing this species at 24 of 40 sites. At 16 of the sites where we detected T. rosae DNA, we did not visually observe the species. Although our assay for Typlichthys eigenmanni successfully amplified the target DNA from the environment, it also resulted in false absences where the species was visually confirmed. Using eDNA to detect cave crayfishes was much more difficult. The assay for Cambarus subterraneus did not work for eDNA samples and we were unable to pick up DNA from the environment, even at locations where it was visually confirmed. Alternatively, the eDNA surveys worked well for C. tartarus and we were able to amplify DNA at every site where it was visually observed. Our assay for C. aculabrum was based on a single sample obtained from GenBank, and did not amplify eDNA from field samples. Lastly, our eDNA results from samples in the known range of Orconectes stygocaneyi suggested the species may be found at an additional cave. Detection using eDNA based on our O. stygocaneyi assay was likely low because it was designed from a pseudogene; however, positive eDNA samples were sequenced to confirm species-specific DNA. Detection probability of both cavefishes and cave crayfishes varied by survey technique and was influenced by water volume, water clarity, water velocity, and substrate. Detection of cavefishes and cave crayfishes via visual surveys decreased when water volume increased, whereas detection using eDNA increased with greater water volume. Detection between taxa using either sample method was highest in habitats classified by fine substrates, except for eDNA detection of crayfishes which was greatest in coarse substrates. Detection of cavefishes increased with water clarity, but detection of cave crayfishes increased with turbidity. Detection probability of both cavefishes and crayfishes using eDNA increased slightly with water velocity, but decreased with visual surveys as water velocity increased. Occupancy by both taxa was positively related to particular geologic series. Crayfish occupancy was negatively related to fine-scale anthropogenic disturbance (i.e., 500-m buffer around the site), whereas crayfish showed no relationship with disturbance. Our results suggest possible range extensions, provide insights to factors driving detection using both sample techniques, and suggest areas where recharge zones may be shared among caves. Future efforts focused on a comprehensive evaluation of genetic diversity among cave crayfishes to improve assay design could improve detection and the applicability of eDNA as a supplemental and non-invasive sampling approach.

Cooperator Science Series

Accuracy assessment/validation methodology and results of 2010–11 land-cover/land-use data for Pools 13, 26, La Grange, and Open River South, Upper Mississippi River System

Introduction/Background The U.S. Geological Survey (USGS)-Upper Midwest Environmental Sciences Center (UMESC) was responsible for development of several land cover/land use (LCU) systemic datasets of the Upper Mississippi River System (UMRS). These efforts (1989 and 2000) were funded by the U.S. Army Corps of Engineers’ Upper Mississippi River Restoration Program (UMRR) Long Term Resource Monitoring (LTRM) element. Development of systemic datasets includes the acquisition, processing, and serving of high-resolution aerial photography and land cover/land use spatial datasets ( http://www.umesc.usgs.gov/data_library/land_cover_use/land_cover_use_data.html ). In 2008, the UMRR reached a collaborative agreement with the U.S. Fish and Wildlife Service-Region 3 to collect high-resolution digital imagery of the entire UMRS floodplain during 2010–11 for LTRM. The UMESC helped acquire, process, and serve this imagery, as well as produce and serve the 2010–11 LCU systemic dataset of the UMRS floodplain. Digital imagery for Pools 13, 26, La Grange, and Open River South was collected using an Applanix DSS 439 digital sensor system with a 40 millimeter lens and Color Infrared (CIR) filter. The imagery was collected at a resolution of 20 centimeters/pixel (8 inches/pixel) for Pool 13 and 40 centimeters/pixel (16 inches/pixel) for Pools 26, Open River South, and La Grange. All imagery was projected to Universal Transverse Mercator (UTM) Zone 15, North American Datum of 1983 (NAD 83). The General Wetland Vegetation Classification (GWVC) system used for mapping is hierarchical, and its 31 classes can be collapsed into broader categories using either a 15- or 7-class level. While the 1989 and 2000 LCU systemic datasets have not gone through a traditional thematic accuracy assessment (AA) in the past, nor have they undergone a validation analysis, the end products are of high quality. For each systemic dataset produced (1989, 2000, 2010–11), extensive field reconnaissance is performed before photointerpretation. The intent of this field reconnaissance is to learn, test, and verify image signatures as they relate to vegetation types. Questionable areas on the imagery are visited, and the plants or land features observed in the area are recorded for reference. This procedure verifies vegetation signatures on the imagery with those on the ground. In addition, once the photointerpretation is complete, the final LCU dataset undergoes extensive quality assurance/quality control to ensure the imagery is mapped correctly. Since the 2000 LCU systemic dataset was developed, there has been a growing interest in completing thematic AAs for the LTRM LCU spatial datasets. The objective of an AA is to measure the probability that a particular location has been assigned its correct vegetation class. An AA estimates thematic (map class) errors in the data, giving users information needed to determine data suitability for a particular application. At the same time, data producers are able to learn more about the nature of errors in the data. Thus, the two attributes of an AA are “producers’ accuracy,” which is the probability that an AA point has been mapped correctly (also referred to as an error of omission); and “users’ accuracy,” which is the probability that the map actually represents what was found on the ground (also referred to as error of commission). Producers’ and users’ accuracies can be obtained from the same set of data by using different analyses. Accuracy assessment is an extensive effort that requires seasonal field personnel and equipment, data entry, analyses, and post processing—tasks that are costly and time consuming. The geospatial team at the UMESC has suggested a validation process for understanding the accuracy of the spatial datasets, which will be tested on at least some areas of the UMRS. Validation is not a true verification of map-class type in the field; however, it can provide the user of the map with useful information that is similar to a field AA. Similar to an AA, validation involves generating random points based on the total area for each map class. However, instead of collecting field data, two or three individuals not involved with the photo-interpretative mapping separately review each of the points onscreen and record a best-fit vegetation type(s) for each site. Once the individual analyses are complete, results are joined together and a comparative analysis is performed. The objective of this initial analysis is to identify areas where the validation results were in agreement (matches) and areas where validation results were in disagreement (mismatches). The two or three individuals then perform an analysis, looking at each mismatched site, and agree upon a final validation class. (If two vegetation types at a specific site appear to be equally prevalent, the validation team is permitted to assign the site two best-fit vegetation types.) Following the validation team’s comparative analysis of vegetation assignments, the data are entered into a database and compared to the mappers’ vegetation assignments. Agreements and disagreements between the map and validation classes are identified, and a contingency table is produced. This document presents the AA processes/results for Pools 13 and La Grange, as well as the validation process/results for Pools 13 and 26 and Open River South.

Iowa, Illinois, Minnesota, Missouri, Wisconsin

Assessment of undiscovered hydrocarbon resources of sub-Saharan Africa

Introduction The main objective of the U.S. Geological Survey’s (USGS) National and Global Petroleum Assessment Project is to assess the potential for undiscovered, technically recoverable oil and natural gas resources of the United States and the world (U.S. Geological Survey World Conventional Resources Assessment Team, 2012). The USGS updated assessments that were completed during the USGS World Petroleum Assessment 2000 (U.S. Geological Survey World Energy Assessment Team, 2000) and conducted new assessments in areas around the world that were not previously examined (U.S. Geological Survey World Conventional Resources Assessment Team, 2012). These assessments used the latest geology-based assessment methodology for conventional oil and gas resources. The new assessments are available at the USGS website, ( http://energy.usgs.gov/OilGas/AssessmentsData/WorldPetroleumAssessment.aspx ) . As part of this project, the USGS assessed 13 geologic provinces located in sub-Saharan Africa (U.S. Geological Survey World Conventional Resources Assessment Team, 2012). Coastal provinces were extended offshore to water depths ranging from 2,000 to 4,000 meters (m). Within these 13 geologic provinces 18 assessment units (figs. 1, 2) were identified. The west Africa provinces are (1) the Senegal, containing the passive-margin Senegal Basin of Middle Jurassic to Holocene age; (2) the West African Coastal, characterized by rift, passive-margin, and transform tectonics; (3) the Gulf of Guinea, characterized by transform tectonics; (4) the Niger Delta, containing more than 9,100 m of sedimentary rock and recent sediments; (5) the West-Central Coastal, which contains the Aptian salt basin, is dominated by both rift and sag tectonics, and includes the Congo Basin; and (6) the Orange River Coastal, containing more than 7,000 m of syn-rift and post-rift sedimentary rock. The West African Coastal Province was assessed for the first time, whereas the other five west Africa provinces were reassessed for the 2012 World Oil and Wandrey Gas Resource Assessment (fig. 1 of U.S. Geological Survey World Conventional Resources Assessment Team, 2012). More than 275 new oil and gas fields have been discovered in the six west Africa provinces (IHS Energy, 2008, 2009) since the USGS World Petroleum Assessment in 2000 (U.S. Geological Survey World Energy Assessment Team, 2000). These provinces were assessed because of increased energy exploration activity and new oil and gas discoveries within the provinces. Seven provinces not assessed as part of the World Petroleum Assessment 2000 (U.S. Geological Survey World Energy Assessment Team, 2000) were assessed by the USGS as part of the World Assessment 2012 (U.S. Geological Survey World Conventional Resources Assessment Team, 2012). These provinces are (1) the Chad Province, containing Cretaceous and Cenozoic-age lacustrine, continental, and minor marine rocks; (2) the Sud Province, containing Cretaceous and Paleogene age lacustrine, continental, and minor marine rocks; (3) the South Africa Coastal Province, which contains rift, transform, and passive-margin rocks; (4) the Mozambique Coastal Province, containing rift, drift, and passive-margin rocks; (5) the Morondava Province, which contains failed rift, drift, and passive-margin rocks; (6) the Tanzania Coastal Province, containing rift, drift, and passive-margin rocks; and (7) the Seychelles Province, which contains rift and drift rocks. At the time of this assessment 157 oil and gas fields had been discovered in the seven provinces (IHS Energy, 2009). These provinces were assessed because of increased interest and new oil and gas discoveries within the provinces. The assessment was geology-based and used the total petroleum system (TPS) concept. The geologic elements of a TPS are hydrocarbon source rocks (source rock maturation and hydrocarbon generation and migration), reservoir rocks (quality and distribution), and traps where hydrocarbon accumulates. Using these geologic criteria, 16 conventional total petroleum systems and 18 assessment units in the 13 provinces were defined. The undiscovered, technically recoverable oil and gas resources were assessed for all assessment units.

Data Series

Geology, ground-water hydrology, geochemistry, and ground-water simulation of the Beaumont and Banning Storage Units, San Gorgonio Pass area, Riverside County, California

Ground water has been the only source of potable water supply for residential, industrial, and agricultural users in the Beaumont and Banning storage units of the San Gorgonio Pass area, Riverside County, California. Ground-water levels in the Beaumont area have declined as much as 100 feet between the early 1920s and early 2000s, and numerous natural springs have stopped flowing. In 1961, the San Gorgonio Pass Water Agency (SGPWA) entered into a contract with the California State Department of Water Resources to receive 17,300 acre-feet per year of water to be delivered by the California State Water Project (SWP) to supplement natural recharge. Currently (2005), a pipeline is delivering SWP water into the area, and the SGPWA is artificially recharging the ground-water system using recharge ponds located along Little San Gorgonio Creek in Cherry Valley with the SWP water. In addition to artificial recharge, SGPWA is considering the direct delivery of SWP water for the irrigation of local golf courses and for agricultural supply in lieu of ground-water pumpage. To better understand the potential hydrologic effects of different water-management alternatives on ground-water levels and movement in the Beaumont and Banning storage units, existing geohydrologic and geochemical data were compiled, new data from a basin-wide ground-water level and water-quality monitoring network were collected, monitoring wells were installed near the Little San Gorgonio Creek recharge ponds, geohydrologic and geochemical analyses were completed, and a ground-water flow simulation model was developed. The San Gorgonio Pass area was divided into several storage units on the basis of mapped or inferred faults. This study addresses primarily the Beaumont and Banning storage units. The geologic units in the study area were generalized into crystalline basement rocks and sedimentary deposits. The younger sedimentary deposits and the surficial deposits are the main water-bearing deposits in the San Gorgonio Pass area. The water-bearing deposits were divided into three aquifers: (1) the perched aquifer, (2) the upper aquifer, and (3) the lower aquifer based on lithologic and downhole geophysical logs. Natural recharge in the San Gorgonio Pass area was estimated using INFILv3, a deterministic distributed- parameter precipitation-runoff model. The INFILv3 model simulated that the potential recharge of precipitation and runoff in the Beaumont and Banning storage units was about 3,710 acre-feet per year and that the potential recharge in 28 sub-drainage basins upstream of the storage units was about 6,180 acre-feet per year. The water supply for the Beaumont and Banning storage units is supplied by pumping ground water from wells in the Canyon (Edgar and Banning Canyons), Banning Bench, Beaumont, and Banning storage units. Total annual pumpage from the Beaumont and Banning storage units ranged from about 1,630 acre-feet in 1936 to about 20,000 acre-feet in 2003. Ground-water levels declined by as much as 100 feet in the Beaumont storage unit from 1926-2003 in response to ground-water pumping of about 450,160 acre-feet during this period. Since ground-water development began in the San Gorgonio Pass area, there have been several sources of artificial recharge to the basin including return flow from applied water on crops, golf courses, and landscape; septic-tank seepage; and infiltration of storm runoff diversions and imported water into recharge ponds. Return flow from applied water and septic-tank seepage was estimated to reach a maximum of about 8,100 acre-feet per year in 2003. Owing to the great depth of water in much of study area (in excess of 150 feet), the return flow and septic-tank seepage takes years to decades to reach the water table. Stable-isotope data indicate that the source of ground-water recharge was precipitation from storms passing through the San Gorgonio Pass as opposed to runoff from the higher altitudes of the San Bernar

Scientific Investigations Report

Simulation of groundwater flow in the Edwards-Trinity and related aquifers in the Pecos County region, Texas

The Edwards-Trinity aquifer is a vital groundwater resource for agricultural, industrial, and public supply uses in the Pecos County region of western Texas. The U.S. Geological Survey completed a comprehensive, integrated analysis of available hydrogeologic data to develop a numerical groundwater-flow model of the Edwards-Trinity and related aquifers in the study area in parts of Brewster, Jeff Davis, Pecos, and Reeves Counties. The active model area covers about 3,400 square miles of the Pecos County region of Texas west of the Pecos River, and its boundaries were defined to include the saturated areas of the Edwards-Trinity aquifer. The model is a five-layer representation of the Pecos Valley, Edwards-Trinity, Dockum, and Rustler aquifers. The Pecos Valley aquifer is referred to as the alluvial layer, and the Edwards-Trinity aquifer is divided into layers representing the Edwards part of the Edwards-Trinity aquifer and the Trinity part of the Edwards-Trinity aquifer, respectively. The calibration period of the simulation extends from 1940 to 2010. Simulated hydraulic heads generally were in good agreement with observed values; 1,684 out of 2,860 (59 percent) of the simulated values were within 25 feet of the observed value. The average root mean square error value of hydraulic head for the Edwards-Trinity aquifer was 34.2 feet, which was approximately 4 percent of the average total observed change in groundwater-level altitude (groundwater level). Simulated spring flow representing Comanche Springs exhibits a pattern similar to observed spring flow. Independent geochemical modeling corroborates results of simulated groundwater flow that indicates groundwater in the Edwards-Trinity aquifer in the Leon-Belding and Fort Stockton areas is a mixture of recharge from the Barilla and Davis Mountains and groundwater that has upwelled from the Rustler aquifer. The model was used to simulate groundwater-level altitudes resulting from prolonged pumping to evaluate sustainability of current and projected water-use demands. Each of three scenarios utilized a continuation of the calibrated model. Scenario 1 extended recent (2008) irrigation and nonirrigation pumping values for a 30-year period from 2010 to 2040. Projected groundwater-level changes in and around the Fort Stockton area under scenario 1 change little from current conditions, indicating that the groundwater system is near equilibrium with respect to recent (2008) pumping stress. Projected groundwater-level declines in the eastern part of the model area ranging from 5.0 to 15.0 feet are likely the result of nonequilibrium conditions associated with recent increases in pumping after a prolonged water-level recovery period of little or no pumping. Projected groundwater-level declines (from 15.0 to 31.0 feet) occurred in localized areas by the end of scenario 1 in the Leon-Belding area. Scenario 2 evaluated the effects of extended recent (2008) pumping rates as assigned in scenario 1 with year-round maximum permitted pumping rates in the Belding area. Results of scenario 2 are similar in water-level decline and extent as those of scenario 1. The extent of the projected groundwater-level decline in the range from 5.0 to 15.0 feet in the Leon-Belding irrigation area expanded slightly (about a 2-percent increase) from that of scenario 1. Maximum projected groundwater-level declines in the Leon-Belding irrigation area were approximately 31.3 feet in small isolated areas. Scenario 3 evaluated the effects of periodic increases in pumping rates over the 30-year extended period. Results of scenario 3 are similar to those of scenario 2 in terms of the areas of groundwater-level decline; however, the maximum projected groundwater-level decline increased to approximately 34.5 feet in the Leon-Belding area, and the extent of the decline was larger in area (about a 17-percent increase) than that of scenario 2. Additionally, the area of projected groundwater-level declines in the eastern part of the model area increased from that of scenario 2—two individual areas of decline coalesced into one larger area. The localized nature of the projected groundwater-level declines is a reflection of the high degree of fractured control on storage and hydraulic conductivity in the Edwards-Trinity aquifer. Additionally, the finding that simulated spring flow is highly dependent on the transient nature of hydraulic heads in the underlying aquifer indicates the importance of adequately understanding and characterizing the entire groundwater system.

Texas

Hydrogeology of northwestern Nassau and northeastern Queens Counties, Long Island, New York

A detailed investigation of the geology and ground-water occurrence in northwestern Nassau and northeastern Queens Counties, N.Y., has been completed by the U.S. Geological Survey, in cooperation with the Nassau County Department of Public Works and the New York State Water Resources Commission. The area, about 63 square miles, includes the peninsulas of Great Neck and Manhasset Neck and the western part of Nassau County as far south as Garden City. Unconsolidated sediments of Pleistocene and Cretaceous aged 200 to 800 feet thick, constitute the ground-water reservoir of the area. The reservoir includes three discrete aquifers, not all found in the same locality, which transect geologic boundaries. The shallow unconfined aquifer consists of those permeable Pleistocene and Cretaceous deposits that lie within the upper part of the zone of saturation, which extends from slightly below sea level to more than 110 feet above sea level. It is an abundant source of water in much of the area, particularly in the southern part, where the aquifer consists of coarse sand and gravel in a glacial-outwash plain. The intermediate or principal aquifer consists largely of that part of the Magothy(?) formation that occurs from about 50 feet below sea level downward to the top of the clay member of the Raritan formation. Locally, however, the aquifer includes Pleistocene deposits which blanket the Magothy(?) formation or lie in channels cut into it. Although the principal aquifer has limited extent to the north, it is the chief source of water in most of the area, except the peninsulas. The lower or deep confined aquifer is found beneath the entire project area, and consists of the Lloyd sand member of the Raritan formation and the hydraulically contiguous Jameco gravel. The lower limit of the aquifer is the bedrock surface; its upper limit is defined by the overlying clay member of the Raritan formation and the Gardiners clay, the latter abutting against the clay member on the north or, in some valleys and embankments, lying directly upon the clay member. Thus, hydraulic continuity exists between the Lloyd sand member of the Raritan formation (Cretaceous) and the Jameco gravel (Pleistocene) in the northern part of Manhasset and Great Necks. The two contiguous clay bodies overlying the aquifer, the clay member of the Raritan formation and the Gardiners clay, form an effective confining bed which probably extends beyond the shore lines of the project area. The deep confined aquifer is particularly important on Manhasset and Great Necks where locally it is the only source of water available for large public-supply or industrial needs. In general, ground-water supplies in sufficient quantity and of excellent quality can be obtained from the three aquifers underlying northwestern Nassau and northeastern Queens Counties, N.Y. Ground-water withdrawals for public supply have increased with population growth and expanded use from an average of about 10 millions of gallons per day in 1940 to 30 mgd in 1957. In addition, about 10 mgd are pumped for various industrial, institutional and private uses. Much of the water pumped by industry is returned to the ground by diffusion wells and recharge basins. However, an increasing amount of water is lost from the ground-water reservoir due to the expanding network of server systems discharging directly to the sea. In the Manhasset Neck and Great Neck areas, ground-water resources are approaching full development in terms of the total available supply, if such development has not been reached already (1959). The chief limiting factor to further expansion lies in the geometry of the aquifers and the reduced thickness of permeable saturated sediments. In the northern part of both peninsulas, the ground-water reservoir is only about 200 to 400 feet thick, and, locally, as much as 75 percent of the strata in the reservoir consists of rather impermeable silt and clay of the Raritan formation or the Gardiners clay. I

Water Supply Paper

An initial abstraction and constant loss model, and methods for estimating unit hydrographs, peak streamflows, and flood volumes for urban basins in Missouri

Streamflow data, basin characteristics, and rainfall data from 39 streamflow-gaging stations for urban areas in and adjacent to Missouri were used by the U.S. Geological Survey in cooperation with the Metropolitan Sewer District of St. Louis to develop an initial abstraction and constant loss model (a time-distributed basin-loss model) and a gamma unit hydrograph (GUH) for urban areas in Missouri. Study-specific methods to determine peak streamflow and flood volume for a given rainfall event also were developed. Distinct basin characteristics were evaluated and selected for use on the basis of their theoretical relation to flow, results from previous studies, and the ability to reliably measure the basin characteristic using digital datasets and geographic information system (GIS) technology. The key basin characteristics determined or computed for each of the 39 basins upstream from the streamflow-gaging stations were drainage area, percent impervious area, main-channel slope based on the 10- and 85-percent length method, percentage of the basin area in storage (lakes, ponds, reservoirs, wetlands), the composite Natural Resources Conservation Service curve number estimated from a combination of the soil type data and land-use characteristics, and the streamflow variability index developed for the recently completed study of low-flow regression in Missouri. Characteristics of spatial and temporal rainfall distribution came from the next generation weather radar (NEXRAD) network. Procedures were developed for this study to convert the variable radar sweep rate into a 5-minute total rainfall hyetograph using data from the radar bin at the centroid of a given basin. Additional characteristics determined for each storm on the basin included the 5-day and 14-day antecedent rainfall, estimated from the mean of daily rainfall values from various rain gages in the area. The database of observed rainfall and runoff events for the 39 basins upstream from the streamflow-gaging stations was analyzed to compute the optimal storm-specific initial abstraction and constant loss values, as well as the time to peak, peak streamflow, and shape factor values of the GUH. The optimal storm-specific values were used to develop a regional regression equation for initial abstraction; conversely, the constant loss was estimated not by regression but from either a generalized or specific regional mean value. The optimal storm-specific values of GUH time to peak, GUH peak streamflow, and GUH shape factor were used to develop regression equations for the GUH. The regression equations for the GUH initially were tested alone, and then were combined with the appropriate regional regression equation for initial abstraction and both the generalized regional and specific regional mean constant loss values. For the GUH regression equations, the interquartile range was substantially smaller than the range spanned by the minimum and maximum values, which indicates most of the errors have much smaller variation, and the minimum and maximum values may be extreme outliers. The central tendency of the regressed errors for peak streamflow and runoff hydrograph volume were both approximately zero, which implies a generally unbiased estimation of these values. The mean and median of the regressed errors for time to peak streamflow were both small but greater than zero, which implies the GUH regression equations create a hydrograph that has a peak that is later in time than observed. Specifically, the regressed times indicate an offset of about 10 minutes, on average, from observed. The mean and median of the regressed errors for widths of the runoff hydrograph at 50 and 75 percent were less than zero, which implies the GUH tends to slightly underestimate these widths compared to the observed. The appropriate regional initial abstraction regression equation was combined with both the generalized and the specific regional mean constant loss values and the GUH regression equations. Both the generalized regional mean constant loss and specific regional mean constant loss forms of the basin-loss model worked equally well to model the observed runoff hydrograph based on the error analysis, and neither model seems to make a consistently better approximation. Both initial abstraction and constant loss models combined with the GUH regression equations were further validated using several storms available after the start of the project in early 2011 with similar but consistently higher error results. If these methods are used in an urban area in Missouri other than those examined in this study, advice to the user is given to consider using the generalized regional mean values. If these methods are used in an urban area that is a subbasin of one of the basins in this study, advice to the user is given to consider using the specific regional mean values. The rainfall-runoff pairs from the storm-specific GUH analysis were further analyzed against various basin and rainfall characteristics to develop equations to estimate the peak streamflow and flood volume based on a quantity of rainfall on the basin.

Missouri

Preliminary geologic map of the San Guillermo Mountain Quadrangle, Ventura County, California

New 1:24,000-scale geologic mapping in the Cuyama 30' x 60' quadrangle, in support of the USGS Southern California Areal Mapping Project (SCAMP), is contributing to a more complete understanding of the stratigraphy, structure, and tectonic evolution of the complex junction area between the NW-striking Coast Ranges and EW-striking western Transverse Ranges. The 1:24,000-scale geologic map of the San Guillermo Mountain quadrangle is one of six contiguous 7 1/2' quadrangle geologic maps in the eastern part of the Cuyama map area being compiled for a more detailed portrayal and reevaluation of geologic structures and rock units shown on previous geologic maps of the area (e.g., Dibblee, 1979). The following observations and interpretations are based on the new San Guillermo Mountain geologic compilation: (1) The new geologic mapping in the northern part of the San Guillermo Mountain quadrangle allows for reinterpretation of fault architecture that bears on potential seismic hazards of the region. Previous mapping had depicted the eastern Big Pine fault (BPF) as a northeast-striking, sinistral strike-slip fault that extends for 30 km northeast of the Cuyama River to its intersection with the San Andreas fault (SAF). In contrast the new mapping indicates that the eastern BPF is a thrust fault that curves from a northeast strike to an east strike, where it is continuous with the San Guillermo thrust fault, and dies out further east about 15 km south of the SAF. This redefined segment of the BPF is a south-dipping, north-directed thrust, with dominantly dip slip components (rakes > 60 deg.), that places Middle Eocene marine rocks (Juncal and Matilija Formations) over Miocene through Pliocene(?) nonmarine rocks (Caliente, Quatal, and Morales Formations). Although a broad northeast-striking fault zone, exhibiting predominantly sinistral components of slip (rakes < 45 deg.), extends to the SAF as previously mapped, the fault zone does not connect to the southwest with the BPF but instead curves into a southwest-directed thrust fault system a short distance north of the BPF. Oligocene to Pliocene(?) nonmarine sedimentary and volcanic rocks of the Plush Ranch, Caliente, and Morales(?) Formations are folded on both sides of this fault zone (informally named the Lockwood Valley fault zone [LVFZ] on the map). South-southeast of the LVFZ overturned folds have southward vergence. Several moderate-displacement (< 50 m), mainly northwest-dipping thrust and reverse faults, exhibiting mostly sinistral-oblique slip, flank and strike parallel to the overturned folds. The fold vergence and thrust direction associated with the LVFZ is opposite to that of the redefined BPF, providing further evidence that the two faults are distinct structures. These revised fault interpretations bring into question earlier estimates of net sinistral strike-slip displacement of as much as 13 km along the originally defined eastern BPF, which assumed structural connection with the LVFZ. Also, despite sparse evidence for repeated Quaternary movement on the LVFZ (e.g., Dibblee, 1982), the potential for a large earthquake involving coseismic slip on both the LVFZ and the central BPF to the southwest may not be as great as once believed. (2) Several generations of Pleistocene and younger dissected alluvial terrace and fan deposits sit at various levels above modern stream channels throughout the quadrangle. These deposits give testimony to the recent uplift and related fault deformation that has occurred in the area. (3) A vast terrane of Eocene marine sedimentary rocks (Juncal and Matilija Formations and Cozy Dell Shale) exposed south of the Big Pine fault forms the southern two-thirds of the San Guillermo Mountain quadrangle. Benthic foraminifers collected from various shale intervals within the Juncal Formation indicate a Middle Eocene age (Ulatisian) for the entire formation (K. McDougall, unpub. data, 1998) and deposition at paleodepths as great as 2,000 m (i.e., lowe

California

Assessment of hydrology, water quality, and trace elements in selected placer-mined creeks in the birch creek watershed near central, Alaska, 2001-05

Executive Summary The U.S. Geological Survey, in cooperation with the Bureau of Land Management, completed an assessment of hydrology, water quality, and trace-element concentrations in streambed sediment of the upper Birch Creek watershed near Central, Alaska. The assessment covered one site on upper Birch Creek and paired sites, upstream and downstream from mined areas, on Frying Pan Creek and Harrison Creek. Stream-discharge and suspended-sediment concentration data collected at other selected mined and unmined sites helped characterize conditions in the upper Birch Creek watershed. The purpose of the project was to provide the Bureau of Land Management with baseline information to evaluate watershed water quality and plan reclamation efforts. Data collection began in September 2001 and ended in September 2005. There were substantial geomorphic disturbances in the stream channel and flood plain along several miles of Harrison Creek. Placer mining has physically altered the natural stream channel morphology and removed streamside vegetation. There has been little or no effort to re-contour waste rock piles. During high-flow events, the abandoned placer-mine areas on Harrison Creek will likely contribute large quantities of sediment downstream unless the mined areas are reclaimed. During 2004 and 2005, no substantial changes in nutrient or major-ion concentrations were detected in water samples collected upstream from mined areas compared with water samples collected downstream from mined areas on Frying Pan Creek and Harrison Creek that could not be attributed to natural variation. This also was true for dissolved oxygen, pH, and specific conductance-a measure of total dissolved solids. Sample sites downstream from mined areas on Harrison Creek and Frying Pan Creek had higher median suspended-sediment concentrations, by a few milligrams per liter, than respective upstream sites. However, it is difficult to attach much importance to the small downstream increase, less than 10 milligrams per liter, in median suspended-sediment concentration for either basin. During low-flow conditions in 2004 and 2005, previously mined areas investigated on Harrison Creek and on Frying Pan Creek did not contribute substantial suspended sediments to sample sites downstream from the mined areas. No substantial mining-related water- or sediment-quality problems were detected at any of the sites investigated in the upper Birch Creek watershed during low-flow conditions. Average annual streamflow and precipitation were near normal in 2002 and 2003. Drought conditions, extreme forest fire impact, and low annual streamflow set apart the 2004 and 2005 summer seasons. Daily mean streamflow for upper Birch Creek varied throughout the period of record-from maximums of about 1,000 cubic feet per second to minimums of about 20 cubic feet per second. Streamflow increased and decreased rapidly in response to rainfall and rapid snowmelt events because the steep slopes, thin soil cover, and permafrost areas in the watershed have little capacity to retain runoff. Median suspended-sediment concentrations for the 115 paired samples from Frying Pan Creek and 101 paired samples from Harrison Creek were less than the 20 milligrams per liter total maximum daily load. The total maximum daily load was set by the U.S. Environmental Protection Agency for the upper Birch Creek basin in 1996. Suspended-sediment paired-sample data were collected using automated samplers in 2004 and 2005, primarily during low-flow conditions. Suspended-sediment concentrations in grab samples from miscellaneous sites ranged from less than 1 milligram per liter during low-flow conditions to 1,386 milligrams per liter during a high-flow event on upper Birch Creek. Streambed-sediment samples were collected at six sites on Harrison Creek, two sites on Frying Pan Creek, and one site on upper Birch Creek. Trace-element concentrations of mercury, lead, and zinc in streambed sedimen

Scientific Investigations Report

Comparison of simulations of land-use specific water demand and irrigation water supply by MF-FMP and IWFM

Two hydrologic models, MODFLOW with the Farm Process (MF-FMP) and the Integrated Water Flow Model (IWFM), are compared with respect to each model’s capabilities of simulating land-use hydrologic processes, surface-water routing, and groundwater flow. Of major concern among the land-use processes was the consumption of water through evaporation and transpiration by plants. The comparison of MF-FMP and IWFM was conducted and completed using a realistic hypothetical case study. Both models simulate the water demand for water-accounting units resulting from evapotranspiration and inefficiency losses and, for irrigated units, the supply from surface-water deliveries and groundwater pumpage. The MF-FMP simulates reductions in evapotranspiration owing to anoxia and wilting, and separately considers land-use-related evaporation and transpiration; IWFM simulates reductions in evapotranspiration related to the depletion of soil moisture. The models simulate inefficiency losses from precipitation and irrigation water applications to runoff and deep percolation differently. MF-FMP calculates the crop irrigation requirement and total farm delivery requirement, and then subtracts inefficiency losses from runoff and deep percolation. In IWFM, inefficiency losses to surface runoff from irrigation and precipitation are computed and subtracted from the total irrigation and precipitation before the crop irrigation requirement is estimated. Inefficiency losses in terms of deep percolation are computed simultaneously with the crop irrigation requirement. The seepage from streamflow routing also is computed differently and can affect certain hydrologic settings and magnitudes ofstreamflow infiltration. MF-FMP assumes steady-state conditions in the root zone; therefore, changes in soil moisture within the root zone are not calculated. IWFM simulates changes in the root zone in both irrigated and non-irrigated natural vegetation. Changes in soil moisture are more significant for non-irrigated natural vegetation areas than in the irrigated areas. Therefore, to facilitate the comparison of models, the changes in soil moisture are only simulated by IWFM for the natural vegetation areas, and soil-moisture parameters in irrigated regions in IWFM were specified at constant values . The IWFM total simulated changes in soil moisture that are related to natural vegetation areas vary from stress period to stress period but are small over the entire two-year period of simulation. In the hypothetical case study, IWFM simulates more evapotranspiration and return flows and less streamflow infiltration than MF-FMP. This causes more simulated surface-water diversions upstream and less simulated water available to downstream farms in IWFM compared to MF-FMP. The evapotranspiration simulated by the two models is well correlated even though the quantity is different. The different approaches used to simulate soil moisture, evapotranspiration, and inefficient losses yield different results for deep percolation and pumpage. In IWFM, deep percolation is a function of soil moisture; therefore, the constant soil-moisture requirement for irrigated regions, assumed for this comparison, results in a constant deep percolation rate. This led to poor correlation with the variable deep percolation rates simulated in MF-FMP, where the deep percolation rate, a fraction of inefficiency losses from precipitation and irrigation, is a function of quasi-steady state infiltration for each soil type and a function of groundwater head. Similarly, the larger simulated evapotranspiration in IWFM is mainly responsible for larger simulated groundwater pumpage demands and related lower groundwater levels in IWFM compared to MF-FMP. Because of the differences in features between MF-FMP and IWFM, the user may find that for certain hydrologic settings one model is better suited than the other. The performance of MF-FMP and IWFM in this particular hypothetical test case, with a fixed framework composed of common initial and boundary conditions and input parameter values, does not necessarily predict the performance of MF-FMP and IWFM in a real-world situation with variable framework and parameter values. These differences may affect the evaluation of policies, projects, or water-balance analysis for some hydrologic settings. Generally, both models are powerful tools that simulate a connected system of aquifer, stream networks, land surface, root zone, and runoff processes. MF-FMP simulated the hypothetical test case in about 4 minutes compared to about 58 minutes for IWFM.

Technical Information Record

Assessment of adult pallid sturgeon fish condition, Lower Missouri River—Application of new information to the Missouri River Recovery Program

During spring 2015, Nebraska Game and Parks Commission (NGPC) biologists noted that pallid sturgeon (Scaphirhynchus albus) were in poor condition during sampling associated with the Pallid Sturgeon Population Assessment Project and NGPC’s annual pallid sturgeon broodstock collection effort. These observations prompted concerns that reduced fish condition could compromise reproductive health and population growth of pallid sturgeon. There was a further concern that compromised condition could possibly be linked to U.S. Army Corps of Engineers management actions and increase jeopardy to the species. An evaluation request was made to the Missouri River Recovery Program and the Effects Analysis Team was chartered to evaluate the issue. Data on all Missouri River pallid sturgeon captures were requested and received from the National Pallid Sturgeon Database. All data were examined for completeness and accuracy; 12,053 records of captures between 200 millimeters fork length (mm FL) and 1,200 mm FL were accepted. We analyzed condition using (1) the condition formula (Kn) from Shuman and others (2011); (2) a second Kn formulation derived from the 12,053 records (hereafter referred to as “Alternative Kn”); and (3) an analysis of covariance (ANCOVA) approach that did not rely on a Kn formulation. The Kn data were analyzed using group (average annual Kn) and individual (percentage in low, normal, and robust conditions) approaches. Using the Shuman Kn formulation, annual mean Kn was fairly static from 2005 to 2011 (although always higher in the upper basin), declined from 2012 to 2015, then remained either static (lower basin) or increasing (upper basin) in 2016. Under the Alternative Kn formulation, the upper basin showed no decline in Kn, whereas the lower basin displayed the same trend as the Shuman Kn formulation. Using both formulations, the individual approach revealed a more complex situation; at the same times and locations that there are fish in poor condition, there are nearby fish in normal or robust condition. The ANCOVA approach revealed that fish condition at size changed between 400 and 600 mm and that some of the apparent trend in low condition was caused by differences in sample size across the size range of the population (that is, greater catch of intermediate-sized fish compared to large fish). We examined basin, year, origin (hatchery compared to wild), segment, and size class for effects on condition and concluded that, since 2012, there has been an increase in the percentage of pallid sturgeon in low condition. There are basin, year, and segment effects; origin and size class do not seem to have an effect. The lower basin, in particular segment 9 (Platte River to Kansas River), had a high percentage of low-condition fish. Within the segment, there were bend-level effects, but the bend effect was not spatially contiguous. We concluded that existing data confirm concerns about declining fish condition, especially in the segments between Sioux City, Iowa, and Kansas City, Missouri. Although the evidence is strong that fish condition has been in decline from 2011 to 2015, additional analysis of individual fish histories may provide more confidence in this conclusion; such analysis was beyond the scope of this effort but is part of our recommendations. The most recent data in 2016 indicate that decline of condition may have leveled off; however, the length of record is insufficient to determine whether recent declines are within the background range of variation. We recommend that monitoring of fish condition should be increased and enhanced with additional health metrics. We also recommend that, should condition continue to decline, processes are deployed to bring low-condition adult fish into the hatchery to improve nutrition and condition. We could not determine the cause of declining fish condition with available data, but we compiled information on several dominant hypotheses in two main categories: inter- or intraspecific competition for resources and habitat conditions. Data are insufficient to indicate a specific causation or solution, and it is possible that multiple causes apply. We make recommendations for additional research that can be pursued to address uncertainties in trends in fish health as well as potential causes.

Missouri River

Hydrological information products for the Off-Project Water Program of the Klamath Basin Restoration Agreement

The Klamath Basin Restoration Agreement (KBRA) was developed by a diverse group of stakeholders, Federal and State resource management agencies, Tribal representatives, and interest groups to provide a comprehensive solution to ecological and water-supply issues in the Klamath Basin. The Off-Project Water Program (OPWP), one component of the KBRA, has as one of its purposes to permanently provide an additional 30,000 acre-feet of water per year on an average annual basis to Upper Klamath Lake through "voluntary retirement of water rights or water uses or other means as agreed to by the Klamath Tribes, to improve fisheries habitat and also provide for stability of irrigation water deliveries." The geographic area where the water rights could be retired encompasses approximately 1,900 square miles. The OPWP area is defined as including the Sprague River drainage, the Sycan River drainage downstream of Sycan Marsh, the Wood River drainage, and the Williamson River drainage from Kirk Reef at the southern end of Klamath Marsh downstream to the confluence with the Sprague River. Extensive, broad, flat, poorly drained uplands, valleys, and wetlands characterize much of the study area. Irrigation is almost entirely used for pasture. To assist parties involved with decisionmaking and implementation of the OPWP, the U.S. Geological Survey (USGS), in cooperation with the Klamath Tribes and other stakeholders, created five hydrological information products. These products include GIS digital maps and datasets containing spatial information on evapotranspiration, subirrigation indicators, water rights, subbasin streamflow statistics, and return-flow indicators. The evapotranspiration (ET) datasets were created under contract for this study by Evapotranspiration, Plus, LLC, of Twin Falls, Idaho. A high-resolution remote sensing technique known as Mapping Evapotranspiration at High Resolution and Internalized Calibration (METRIC) was used to create estimates of the spatial distribution of ET. The METRIC technique uses thermal infrared Landsat imagery to quantify actual evapotranspiration at a 30-meter resolution that can be related to individual irrigated fields. Because evaporation uses heat energy, ground surfaces with large ET rates are left cooler as a result of ET than ground surfaces that have less ET. As a consequence, irrigated fields appear in the Landsat images as cooler than nonirrigated fields. Products produced from this study include total seasonal and total monthly (April-October) actual evapotranspiration maps for 2004 (a dry year) and 2006 (a wet year). Maps showing indicators of natural subirrigation were also provided by this study. "Subirrigation" as used here is the evapotranspiration of shallow groundwater by plants with roots that penetrate to or near the water table. Subirrigation often occurs at locations where the water table is at or above the plant rooting depth. Natural consumptive use by plants diminishes the benefit of retiring water rights in subirrigated areas. Some agricultural production may be possible, however, on subirrigated lands for which water rights are retired. Because of the difficulty in precisely mapping and quantifying subirrigation, this study presents several sources of spatially mapped data that can be used as indicators of higher subirrigation probability. These include the floodplain boundaries defined by stream geomorphology, water-table depth defined in Natural Resources Conservation Service (NRCS) soil surveys, and soil rooting depth defined in NRCS soil surveys. The two water-rights mapping products created in the study were "points of diversion" (POD) and "place of use" (POU) for surface-water irrigation rights. To create these maps, all surface-water rights data, decrees, certificates, permits, and unadjudicated claims within the entire 1,900 square mile study area were aggregated into a common GIS geodatabase. Surface-water irrigation rights within a 5-mile buffer of the study area were then selected and identified. The POU area was then totaled by water right for primary and supplemental water rights. The maximum annual volume (acre-feet) allowed under each water right also was calculated using the POU area and duty (allowable annual irrigation application in feet). In cases where a water right has more than one designated POD, the total volume for the water right was equally distributed to each POD listed for the water right. Because of this, mapped distribution of diversion rates for some rights may differ from actual practice. Water-right information in the map products was from digital datasets obtained from the Oregon Water Resources Department and was, at the time acquired, the best available compilation of water-right information available. Because the completeness and accuracy of the water-right data could not be verified, users are encouraged to check directly with the Oregon Water Resources Department where specific information on individual rights or locations is essential. A dataset containing streamflow statistics for 72 subbasins in the study area was created for the study area. The statistics include annual flow durations (5-, 10-, 25-, 50-, and 95-percent exceedances) and 7-day, 10-year (7Q10) and 7-day, 2-year (7Q2) low flows, and were computed using regional regression equations based on measured streamflow records in the region. Daily streamflow records used were adjusted as needed for crop consumptive use; therefore the statistics represent streamflow under more natural conditions as though irrigation diversions did not exist. Statistics are provided for flow rates resulting from streamflow originating from within the entire drainage area upstream of the subbasin pour point (referring to the outlet of the contributing drainage basin). The statistics were computed for the purpose of providing decision makers with the ability to estimate streamflow that would be expected after water conservation techniques have been implemented or a water right has been retired. A final product from the study are datasets of indicators of the potential for subsurface return flow of irrigation water from agricultural areas to nearby streams. The datasets contain information on factors such as proximity to surface-water features, geomorphic floodplain characteristics, and depth to water. The digital data, metadata, and example illustrations for the datasets described in this report are available on-line from the USGS Water Resources National Spatial Data Infrastructure (NSDI) Node Website http://water.usgs.gov/lookup/getgislist or from the U.S. Government website DATA.gov at http://www.data.gov with links provided in a Microsoft&reg; Excel&reg; workbook in appendix A.

California;Oregon

The potential effects of sodium bicarbonate, a major constituent from coalbed natural gas production, on aquatic life

The production water from coalbed natural gas (CBNG) extraction contains many constituents. The U.S. Environmental Protection Agency has established aquatic life criteria for some of these constituents, and it is therefore possible to evaluate their risk to aquatic life. However, of the major ions associated with produced waters, chloride is the only one with an established aquatic life criterion (U.S. Environmental Protection Agency, 1988). The focus of this research was NaHCO 3 , a compound that is a major constituent of coalbed natural gas produced waters in the Tongue and Powder River Basins. This project included laboratory experiments, field in situ experiments, a field mixing zone study, and a fishery presence/absence assessment. Though this investigation focuses on the Tongue and Powder River Basins, the information is applicable to other watersheds where sodium bicarbonate is a principle component of product water either from CBNG or from traditional or unconventional oil and gas development. These data can also be used to separate effects of saline discharges from those potentially posed by other constituents. Finally, this research effort and the additional collaboration with USGS Water Resources and Mapping, Bureau of Land Management, US Environmental Protection Agency, State of Montana, State of Wyoming, Montana State University, University of Wyoming, and others as part of a Powder River Aquatic Task Group, can be used as a model for successful approaches to studying landscapes with energy development. The laboratory acute toxicity experiments were completed with a suite of organisms, including 7 species of fish, 5 species of invertebrates, and 1 amphibian species. Experiments performed on these multiple species resulted in LC50s that ranged from 1,120 to greater than (>) 8,000 milligrams sodium bicarbonate per liter (mg NaHCO 3 /L) (also defined as 769 to >8,000 milligrams bicarbonate per liter (mg HCO 3 -/L) or total alkalinity expressed as 608 to >4,181 milligrams calcium carbonate per liter (mg CaCO 3 /L)) that varied across species and lifestage within a species. The age at which fish were exposed to NaHCO 3 significantly affected the severity of toxic responses for some organisms. The chronic toxicity of NaHCO 3 was defined in experiments that lasted from 7&mdash;60 days post-hatch. For these experiments, sublethal effects such as growth and reproduction, in addition to significant reductions in survival were included in the final determination of effects. Chronic toxicity was observed at concentrations that ranged from 450 to 800mg NaHCO 3 /L (also defined as 430 to 657 mg HCO 3 -/L or total alkalinity expressed as 354 to 539 mg CaCO 3 /L) and the specific concentration depended on the sensitivity of the four species of invertebrates and fish exposed. Sublethal investigations during chronic studies revealed percent decrease in the activity of sodium-potassium adenosine triphosphatase (Na/K ATPase, an enzyme involved in ionoregulation) and the age of the fish at the onset of the decrease may affect the ability of fathead minnow to survive exposures to NaHCO 3 . A database of toxicity evaluations of NaHCO 3 on aquatic life has been constructed. Using these data, sample acute and chronic criteria of 459 and 381 mg NaHCO 3 /L, respectively, can be calculated for the protection of aquatic life. The final derivation and implementation of such criteria is, of course, left to the discretion of the concerned management agencies. A combination of in situ experiments, static-renewal experiments performed simultaneously with in situ experiments, and static renewal experiments performed with site water in the laboratory, demonstrated that untreated coalbed natural gas (CBNG) product water from the Tongue and Powder River Basins reduces survival of fathead minnow and pallid sturgeon. More precisely, the survival of early-lifestage fathead minnow, especially those less than 6-days post hatch (dph), likely is reduced significantly in the field when concentrations of NaHCO 3 rise above 1,500 mg/L. However, age was not a factor for pallid sturgeon and they were sensitive to product water regardless of age. Treatment with the Higgins Loop&trade; technology and dilution of untreated water increased survival in the laboratory. Both of these situations reduced ammonia in addition to the concentrations of NaHCO 3 . These experiments addressed the acute toxicity of effluent waters being added to the main stem rivers, but did not address issues related to the volumes of water that may be added to the watershed. Mixing zones of the three outfalls studied ranged from approximately 800&mdash;1,200 m below the confluence and the areas within these mixing zones with acutely lethal concentrations of NaHCO 3 (as defined by the presence of concentrated dye) are limited. The areas with concentrations of NaHCO 3 more than the concentrations likely to cause significant mortality, and more than the calculated sample water-quality criteria in the Tongue and Powder River Basins appear to be limited to tributaries and parts of mixing zones with considerable additions of untreated effluent.

Montana;Wyoming

Results of 2018–19 water-quality and hydraulic characterization of aquifer intervals using packer tests and preliminary geophysical-log correlations for selected boreholes at and near the former Naval Air Warfare Center Warminster, Bucks County, Pennsylvania

The U.S. Geological Survey (USGS) collected data on the vertical distribution of hydraulic head, specific capacity, and water quality using aquifer-interval-isolation tests and other vertical profiling methods in 15 boreholes completed in fractured sedimentary bedrock in Northampton, Warminster, and Warwick Townships, Bucks County, Pennsylvania during 2018–19. This work was done, in cooperation with the U.S. Navy, to support detailed investigations at and near the former Naval Air Warfare Center (NAWC) Warminster, where groundwater contamination with per- and polyfluoroalkyl substances (PFAS) had become a concern since 2014. Two PFAS compounds, perfluorooctanesulfonic acid (PFOS) and perfluorooctanoic acid (PFOA), have been measured in groundwater samples from supply and monitoring wells at or near NAWC Warminster in concentrations above U.S. Environmental Protection Agency health advisory levels for drinking water. The area is underlain by the Triassic Stockton Formation, which predominantly consists of sandstone interbedded with shale and siltstone beds and forms a layered fractured-rock aquifer used for private, industrial, and public drinking water supply. The vertical distribution of aquifer properties and water quality was assessed through hydraulic tests and sampling of aquifer intervals using a straddle-packer system (13 boreholes) or depth-discrete point sampling under known borehole-flow conditions (2 boreholes). Geophysical and video logs collected by USGS during 2017–19 were used to identify potential water-bearing fractures in 15 boreholes, which ranged in depth from 210 to 604 feet (ft) and included 6 boreholes drilled in 2018 and 9 existing wells on or near the former NAWC Warminster. Measured borehole flow was predominantly downward in most of the deepest boreholes (greater than 400 ft), which were commonly located at the highest land-surface elevations, with inflow from fractures at relatively shallow depths and outflow through fractures near or below depths of 500 ft below land surface. Hydraulic head differences measured during packer tests were up to about 60 ft between shallow and deep intervals. Borehole flow was predominantly upward in most boreholes less than 400 ft in depth and farther from, and at lower land-surface elevations than, the former NAWC Warminster. Total borehole specific capacity ranged from about 0.07 to 41 gallons per minute per foot [(gal/min)/ft]. Specific-capacity values for individual intervals ranged from 0.02 to 40.0 (gal/min)/ft, with a median of 1.14 (gal/min)/ft and a large range in values at most depths. Differences in water quality of samples as indicated by field properties (pH, dissolved oxygen, and specific conductance) and concentrations of dissolved major ions, PFOA, and PFOS were apparent among isolated intervals in the boreholes. Summed concentrations of PFOA and PFOS ranged from about 11 to 10,780 nanograms per liter (ng/L) and were greater than the 2016 U.S. Environmental Protection Agency health advisory of 70 ng/L for summed PFOA and PFOS concentrations in 62 of 104 intervals and discrete depths tested. The mass ratio of PFOS to PFOA was generally higher than 1.0 in samples with summed PFOA and PFOS concentrations greater than 70 ng/L, with ratio values as high as 8.7. In many boreholes, summed concentrations of PFOA and PFOS were positively related to chloride concentrations, which were elevated above natural-background values [less than 10 milligrams per liter] in most samples and as high as 717 milligrams per liter. Sources of the elevated chloride other than, or in addition to, common rock salt (sodium chloride) were indicated by chloride to sodium molar ratios greater than 1.0. Water-quality data indicated that sampled water from some intervals with lower hydraulic heads may be affected by water from intervals with higher hydraulic heads because of vertical flow in open boreholes; samples from these intervals with lower hydraulic heads may not be fully representative due to some component of cross contamination and should be interpreted with caution. Through a preliminary correlation of natural gamma and resistivity logs of boreholes drilled at and near the former NAWC Warminster, 11 lithologic units were identified and interpreted to strike northeast and dip to the northwest. Hydraulic heads were generally highest in isolated intervals that intercepted beds which, when projected up dip, crop out at the highest land-surface elevation on the former NAWC Warminster, indicating that the dipping-bed structure and topography are factors affecting the distribution of hydraulic head in the aquifer. The hydrogeologic framework in conjunction with the vertical distribution of hydraulic heads and water quality may assist in evaluating the locations of various PFAS sources and potential migration pathways of PFAS in groundwater at and near NAWC Warminster.

Pennsylvania