USGS Science⌕ Search

SEARCH · USGS Science

Results for “Internal 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,639 records · Page 91Linked to original sources
USGS2006

Application of municipal biosolids to dry-land wheat fields - A monitoring program near Deer Trail, Colorado (USA). A presentation for an international conference: "The Future of Agriculture: Science, Stewardship, and Sustainability", August 7-9, 2006, Sacramento, CA

Since late 1993, Metro Wastewater Reclamation District of Denver (Metro District), a large wastewater treatment plant in Denver, Colorado, has applied Grade I, Class B biosolids to about 52,000 acres of non-irrigated farmland and rangeland near Deer Trail, Colorado. In cooperation with the Metro District in 1993, the U.S. Geological Survey (USGS) began monitoring ground water at part of this site. In 1999, the USGS began a more comprehensive study of the entire site to address stakeholder concerns about the chemical effects of biosolids applications. This more comprehensive monitoring program has recently been extended through 2010. Monitoring components of the more comprehensive study included biosolids collected at the wastewater treatment plant, soil, crops, dust, alluvial and bedrock ground water, and stream bed sediment. Streams at the site are dry most of the year, so samples of stream bed sediment deposited after rain were used to indicate surface-water effects. This presentation will only address biosolids, soil, and crops. More information about these and the other monitoring components are presented in the literature (e.g., Yager and others, 2004a, b, c, d) and at the USGS Web site for the Deer Trail area studies at http://co.water.usgs.gov/projects/CO406/CO406.html. Priority parameters identified by the stakeholders for all monitoring components, included the total concentrations of nine trace elements (arsenic, cadmium, copper, lead, mercury, molybdenum, nickel, selenium, and zinc), plutonium isotopes, and gross alpha and beta activity, regulated by Colorado for biosolids to be used as an agricultural soil amendment. Nitrogen and chromium also were priority parameters for ground water and sediment components. In general, the objective of each component of the study was to determine whether concentrations of priority parameters (1) were higher than regulatory limits, (2) were increasing with time, or (3) were significantly higher in biosolids-applied areas than in a similar farmed area where biosolids were not applied. Where sufficient samples could be collected, statistical methods were used to evaluate effects. Rigorous quality assurance was included in all aspects of the study. The roles of hydrology and geology also were considered in the design, data collection, and interpretation phases of the study. Study results indicate that the chemistry of the biosolids from the Denver plant was consistent during 1999-2005, and total concentrations of regulated trace elements were consistently lower than the regulatory limits. Plutonium isotopes were not detected in the biosolids. Leach tests using deionized water to simulate natural precipitation indicate arsenic, molybdenum, and nickel were the most soluble priority parameters in the biosolids. Study results show no significant difference in concentrations of priority parameters between biosolids-applied soils and unamended soils where no biosolids were applied. However, biosolids were applied only twice during 1999-2003. The next soil sampling is not scheduled until 2010. To date concentrations of most of the priority parameters were not much greater in the biosolids than in natural soil from the sites. Therefore, many more biosolids applications would need to occur before biosolids effects on the soil priority constituents can be quantified. Leach tests using deionized water to simulate precipitation indicate that molybdenum and selenium were the priority parameters that were most soluble in both biosolids-applied soil and natural or unamended soil. Study results do not indicate significant differences in concentrations of priority parameters between crops grown in biosolids-applied areas and crops grown where no biosolids were applied. However, crops were grown only twice during 1999-2003, so only two crop samples could be collected. The wheat-grain elemental data collected during 1999-2003 for both biosolids-applied areas and unamended areas are similar

Open-File Report↗
USGS2003

Preliminary geologic map of the Murrieta 7.5' quadrangle, Riverside County, California

The Murrieta quadrangle is located in the northern part of the Peninsular Ranges Province and includes parts of two structural blocks, or structural subdivisions of the province. The quadrangle is diagonally crossed by the active Elsinore fault zone, a major fault zone of the San Andreas fault system, and separates the Santa Ana Mountains block to the west from the Perris block to the east. Both blocks are relatively stable internally and within the quadrangle are characterized by the presence of widespread erosional surfaces of low relief. The Santa Ana Mountains block, in the Murrieta quadrangle, is underlain by undifferentiated, thick-layered, granular, impure quartzite and well-layered, fissile, phyllitic metamorphic rock of low metamorphic grade. Both quartzite and phyllitic rocks are Mesozoic. Unconformably overlying the metamorphic rocks are remnants of basalt flows having relatively unmodified flow surfaces. The age of the basalt is about 7-8Ma. Large shallow depressions on the surface of the larger basalt remnants form vernal ponds that contain an endemic flora. Beneath the basalt the upper part of the metamorphic rocks is deeply weathered. The weathering appears to be the same as the regional Paleocene saprolitic weathering in southern California. West of the quadrangle a variable thickness sedimentary rock, physically resembling Paleogene rocks, occurs between the basalt and metamorphic rock. Where not protected by the basalt, the weathered rock has been removed by erosion. The dominant feature on the Perris block in the Murrieta quadrangle is the south half of the Paloma Valley ring complex, part of the composite Peninsular Ranges batholith. The complex is elliptical in plan view and consists of an older ring-dike with two subsidiary short-arced dikes that were emplaced into gabbro by magmatic stoping. Small to large stoped blocks of gabbro are common within the ring-dikes. A younger ring-set of hundreds of thin pegmatite dikes occur largely within the central part of the complex. These pegmatite dikes were emplaced into a domal fracture system, apparently produced by cauldron subsidence, and include in the center of the complex, a number of flat-floored granophyre bodies. The granophyre is interpreted to be the result of pressure quenching of pegmatite magma. Along the eastern edge of the quadrangle is the western part of a large septum of medium metamorphic grade Mesozoic schist. A dissected basalt flow caps the Hogbacks northeast of Temecula, and represents remnants of a channel filling flow. Beneath the basalt is a thin deposit of stream gravel. Having an age of about 10Ma, this basalt is about 2-3Ma older than the basalt flows in the Santa Ana Mountains. The Elsinore fault zone forms a complex of pull-apart basins. The west edge of the fault zone, the Willard Fault, is marked by the high, steep eastern face of the Santa Ana Mountains. The east side of the zone, the Wildomar Fault, forms a less pronounced physiographic step. In the center of the quadrangle a major splay of the fault zone, the Murrieta Hot Springs Fault, strikes east. Branching of the fault zone causes the development of a broad alluvial valley between the Willard Fault and the Murrieta Hot Springs Fault. All but the axial part of the zone between the Willard and Wildomar Faults consist of dissected Pleistocene sedimentary units. The axial part of the zone is underlain by Holocene and latest Pleistocene sedimentary units.

California↗
USGS1994

Principal reference section for part of the Eocene Ghazij Formation, Moghal Mine area, Mach coal field, Balochistan, Pakistan

The information presented on this sheet was collected as part of a joint U.S. Geological Survey-Geological Survey of Pakistan program sponsored by the U.S. Agency for International Development. As a project within this program, the coal-bearing Ghazij Formation (Eocene) was investigated in the northeastern part of Balochistan cast and south of the provincial capital of Quetta. Strata exposed in this area range in age from Permian to Holocene and crop out as a belt of folded and thrusted rocks that form a southeast-facing orocline. In this region of Pakistan, the Ghazij can usually be divided into three parts. The lower part is the thickest (probably more than 1,000 m) and consists of gray-weathering calcareous mudrock (shale, mudstone, and impure claystone) and a few tabular bodies of fine-to medium-grained calcareous sandstone. The middle part (27-300 m) consists of gray-weathering calcareous mudrock and tabular to lenticular bodies of fine-to medium-grained calcareous sandstone; beds of carbonaceous shale and coal are common. The upper part (as thick as 533 m) contains reddish-weathering calcareous mudrock that contains scattered lenticular bodies of fine- to medium-grained calcareous sandstone. Fossil plant debris is common in mudrock of the lower and middle parts of the Ghazij, and bivalves and gastropods are common in the middle part; the upper part of the Ghazij is usually unfossiliferous. This three-fold division of the Ghazij is less distinct in the Johan area. Here, the upper part of the formation is clearly identifiable, but rocks below it are poorly exposed and assigning a stratigraphic level that separates the middle and lower parts of the formation is problematic. Below the upper part of the formation is a thick sequence of greenish-gray calcareous mudrock that contains locally abundant plant debris and isolated bodies of brown-weathering sandstone. Rare carbonaceous shale and even rarer coal are present in the upper part of this sequence, and this interval of the formation might correspond to the middle part of the Ghazji exposed in areas to the north. We propose that, in the Johan area, those rocks below the upper part of the formation be referred to as the main body of the Ghazij (for example, main-body Ghazij). Underlying the Ghazij are the carbonate rocks of the Paleocene Dungan Formation (or its equivalent), and overlying the Ghazij are the mostly carbonate rocks of the Eocene Kirthar Formation (or its equivalent). Both contacts can be conformable or unconformable. All of the pre-Neogene rocks in Balochistan are greatly deformed by the collision of India and Asia. The Ghazij is especially susceptible to regional compressional tectonics because it contains a large amount of shale and is sandwiched between two thick carbonate units. As a result, bedding-plane faults and isoclinal folds are common. As part of our study of the Ghazij Formation, five stratigraphic sections were measured: one near Pir Ismail Ziarat, one in the Sor Range, two in the vicinity of Mach, and one near Johan. Each area's section is published separately.

Balochistan↗
USGS1994

Principal reference section for part of the Eocene Ghazij Formation, Sarawan River area, Johan coal field, Balochistan, Pakistan

The information presented on this sheet was collected as part of a joint U.S. Geological Survey-Geological Survey of Pakistan program sponsored by the U.S. Agency for International Development. As a project within this program, the coal-bearing Ghazij Formation (Eocene) was investigated in the northeastern part of Balochistan cast and south of the provincial capital of Quetta. Strata exposed in this area range in age from Permian to Holocene and crop out as a belt of folded and thrusted rocks that form a southeast-facing orocline. In this region of Pakistan, the Ghazij can usually be divided into three parts. The lower part is the thickest (probably more than 1,000 m) and consists of gray-weathering calcareous mudrock (shale, mudstone, and impure claystone) and a few tabular bodies of fine-to medium-grained calcareous sandstone. The middle part (27-300 m) consists of gray-weathering calcareous mudrock and tabular to lenticular bodies of fine-to medium-grained calcareous sandstone; beds of carbonaceous shale and coal are common. The upper part (as thick as 533 m) contains reddish-weathering calcareous mudrock that contains scattered lenticular bodies of fine- to medium-grained calcareous sandstone. Fossil plant debris is common in mudrock of the lower and middle parts of the Ghazij, and bivalves and gastropods are common in the middle part; the upper part of the Ghazij is usually unfossiliferous. This three-fold division of the Ghazij is less distinct in the Johan area. Here, the upper part of the formation is clearly identifiable, but rocks below it are poorly exposed and assigning a stratigraphic level that separates the middle and lower parts of the formation is problematic. Below the upper part of the formation is a thick sequence of greenish-gray calcareous mudrock that contains locally abundant plant debris and isolated bodies of brown-weathering sandstone. Rare carbonaceous shale and even rarer coal are present in the upper part of this sequence, and this interval of the formation might correspond to the middle part of the Ghazji exposed in areas to the north. We propose that, in the Johan area, those rocks below the upper part of the formation be referred to as the main body of the Ghazij (for example, main-body Ghazij). Underlying the Ghazij are the carbonate rocks of the Paleocene Dungan Formation (or its equivalent), and overlying the Ghazij are the mostly carbonate rocks of the Eocene Kirthar Formation (or its equivalent). Both contacts can be conformable or unconformable. All of the pre-Neogene rocks in Balochistan are greatly deformed by the collision of India and Asia. The Ghazij is especially susceptible to regional compressional tectonics because it contains a large amount of shale and is sandwiched between two thick carbonate units. As a result, bedding-plane faults and isoclinal folds are common. As part of our study of the Ghazij Formation, five stratigraphic sections were measured: one near Pir Ismail Ziarat, one in the Sor Range, two in the vicinity of Mach, and one near Johan. Each area's section is published separately.

Balochistan↗
USGS2007

Ground-water recharge in the arid and semiarid southwestern United States

Ground-water recharge in the arid and semiarid southwestern United States results from the complex interplay of climate, geology, and vegetation across widely ranging spatial and temporal scales. Present-day recharge tends to be narrowly focused in time and space. Widespread water-table declines accompanied agricultural development during the twentieth century, demonstrating that sustainable ground-water supplies are not guaranteed when part of the extracted resource represents paleorecharge. Climatic controls on ground-water recharge range from seasonal cycles of summer monsoonal and winter frontal storms to multimillennial cycles of glacial and interglacial periods. Precipitation patterns reflect global-scale interactions among the oceans, atmosphere, and continents. Large-scale climatic influences associated with El Niño and Pacific Decadal Oscillations strongly, but irregularly, control weather in the study area, so that year-to-year variations in precipitation and ground-water recharge are large and difficult to predict. Proxy data indicate geologically recent periods of naturally occurring multidecadal droughts unlike any in the modern instrumental record. Any anthropogenically induced climate change will likely reduce ground-water recharge through diminished snowpack at higher elevations. Future changes in El Niño and monsoonal patterns, both crucial to precipitation in the study area, are highly uncertain in current models. Current land-use modifications influence ground-water recharge through vegetation, irrigation, and impermeable area. High mountain ranges bounding the study area—the San Bernadino Mountains and Sierra Nevada to the west, and the Wasatch and southern Colorado Rocky Mountains to the east—provide external geologic controls on ground-water recharge. Internal geologic controls stem from tectonic processes that led to numerous, variably connected alluvial-filled basins, exposure of extensive Paleozoic aquifers in mountainous recharge areas, and distinct modes of recharge in the Colorado Plateau and Basin and Range subregions. The chapters in this professional paper present (first) an overview of climatic and hydrogeologic framework (chapter A), followed by a regional analysis of ground-water recharge across the entire study area (chapter B). These are followed by an overview of site-specific case studies representing different subareas of the geographically diverse arid and semiarid southwestern United States (chapter C); the case studies themselves follow in chapters D–K. The regional analysis includes detailed hydrologic modeling within the framework of a high-resolution geographic-information system (GIS). Results from the regional analysis are used to explore both the distribution of ground-water recharge for mean climatic conditions as well as the influence of two climatic patterns—the El Niño-Southern Oscillation and Pacific Decadal Oscillation—that impart a high degree of variability to the hydrologic cycle. Individual case studies employ a variety of geophysical and geochemical techniques to investigate recharge processes and relate the processes to local geologic and climatic conditions. All of the case studies made use of naturally occurring tracers to quantify recharge. Thermal and geophysical techniques that were developed in the course of the studies are presented in appendices. The quantification of ground-water recharge in arid settings is inherently difficult due to the generally low amount of recharge, its spatially and temporally spotty nature, and the absence of techniques for directly measuring fluxes entering the saturated zone from the unsaturated zone. Deep water tables in arid alluvial basins correspond to thick unsaturated zones that produce up to millennial time lags between changes in hydrologic conditions at the land surface and subsequent changes in recharge to underlying ground water. Recent advances in physical, chemical, isotopic, and modeling techniques have fostered new types of recharge assessments. Chemical and isotopic techniques include an increasing variety of environmental tracers that are useful and robust. Physically based techniques include the use of heat as a tracer and computationally intensive geophysical imaging tools for characterizing hydrologic conditions in the unsaturated zone. Modeling-based techniques include spatially distributed water-budget computations using high-resolution remotely sensed and ground-based geographic data. Application of these techniques to arid and semiarid settings in the southwestern United States reveals distinct patterns of recharge corresponding to geologic setting, climatic and vegetative history, and land use. Analysis of recharge patterns shows that large expanses of alluvial basin floors are drying out under current climatic conditions, with little to no recharge to underlying ground water. Ground-water recharge occurs mainly beneath upland catchments in which thin soils overlie permeable bedrock, ephemeral channels in which flow may average only several hours per year, and active agricultural areas. The chapters in this professional paper represent a coordinated attempt to develop a better understanding of one of the Nation's most critical yet difficult-to-quantify renewable resources.

Professional Paper↗
USGS2014

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

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

Nevada↗
USGS2019

Assessing potential effects of highway and urban runoff on receiving streams in total maximum daily load watersheds in Oregon using the stochastic empirical loading and dilution model

The Stochastic Empirical Loading and Dilution Model (SELDM) was developed by the U.S. Geological Survey (USGS) in cooperation with the Federal Highway Administration to simulate stormwater quality. To assess the effects of runoff, SELDM uses a stochastic mass-balance approach to estimate combinations of pre-storm streamflow, stormflow, highway runoff, event mean concentrations (EMCs) and stormwater constituent loads from a site of interest. In addition, SELDM can be used to assess the effects of stormwater Best Management Practices (BMPs), which are designed to mitigate the adverse effects of runoff into a waterbody. Adverse effects of stormwater on receiving waters are one of the greatest unsolved water-quality problems Nationwide. State DOTs, municipalities, Federal facilities, and private property owners who manage impervious surfaces need information about the potential magnitude of their contributions and the potential effectiveness of methods to mitigate the adverse effects of runoff. Because the efficacy of at-site controls are limited, information about the potential effectiveness of alternative strategies is needed. The USGS, in cooperation with the Oregon Department of Transportation (ODOT), conducted a study to research methods in which SELDM can be used to enhance the efficiency of ODOT’s stormwater program, support the development of a stormwater banking program, and meet environmental goals. Results can be used to develop a strategic, systems-level approach to stormwater management by considering entire watersheds instead of individual road crossings. Two watersheds, Bear Creek and Mill Creek, in western Oregon were selected for analysis. Within each watershed, seven road crossings were selected for demonstrating the utility of SELDM in nested basins. Precipitation statistics, pre-storm streamflow, runoff coefficients, and hydrograph recession factors were calculated for each location and used in SELDM to simulate flow, water-quality concentrations, and constituent loads in the upstream basin, from the highway (or developed area), and downstream from the road crossing. Three water-quality constituents were selected for modeling: suspended-sediment concentration (SSC), total phosphorus (TP), and total copper (TCu). Using water-quality transport curves, the relations between streamflow and SSC and between streamflow and TP were simulated. Concentrations of TCu were simulated by configuring a linear relation between SSC and TCu. A generic BMP was simulated using the median treatment statistics for flow reductions, hydrograph extensions, concentration reductions, and minimum irreducible concentrations from nine BMP categories with data from the 2012 International BMP database. Five simulation scenarios were modeled for demonstrative purposes. These simulations were used to evaluate potential effects of different watershed properties, water-quality inputs, and stormwater mitigation measures. Instream EMCs were compared to hypothetical water-quality criteria for suspended sediment, total phosphorus, and total copper to demonstrate the concept of water-quality risk analysis. For all five scenarios, it was assumed that highway runoff concentrations were independent of location or average annual daily traffic. These five scenarios are as follows: • Simulation Scenario 1—Natural Conditions (hereafter Simulation Scenario 1) represents conditions in an undeveloped watershed. This scenario demonstrates that the strategic placement of a hypothetical road crossing within a watershed could be used to avoid exceeding water-quality standards of TP and SSC, but that no location choice results in meeting TCu standards. Implementation of BMP had the most pronounced effects on downstream water-quality constituent EMCs at road crossings with the highest ratio of highway catchment area to upstream drainage area, but the largest effect of BMP treatment on mean annual load is based on highway catchment area alone. • Simulation Scenario 2—Current Conditions (hereafter Simulation Scenario 2) represents current watershed conditions, where all developed area upstream from the road crossing was modeled as a highway and combined with the undeveloped part of the upstream drainage area (scenario 2A) and where the output from scenario 2A is used for the upstream area (developed area and the undeveloped area), and where the road crossing is added as usual (scenario 2B). Scenario 2 results indicate that attaining water-quality standards is more difficult with upstream developed areas. Specific road-crossing sites can be selected to achieve the fewest water-quality exceedances per year, but water-quality targets are not met without BMP implementation, and in some instances are not achievable even with BMP implementation. Results from this scenario also serve to quantify the upper limit of constituent reduction if funding were available to implement BMPs to large areas of development, and to quantify how much area would need BMP implementation to achieve water-quality targets. • Simulation Scenario 3—Alternative Road Layouts (hereafter Simulation Scenario 3) was designed to assess the sensitivity of SELDM to various road layouts. In this scenario, different highway configurations were superimposed at one road crossing. Results indicate that downstream waterquality constituent EMCs did not exhibit much variation, but annual water-quality constituent loads varied considerably. • Simulation Scenario 4—Varying Road Width (hereafter Simulation Scenario 4) was designed to assess the sensitivity of SELDM to road width. Similar to scenario 3, the results indicate little variation in downstream water-quality constituent EMCs, but annual water-quality constituent loads increased in proportion to road width. • Simulation scenario 5—Changes to Impervious Area (hereafter Simulation Scenario 5) was designed to investigate the effects of changing amounts of imperviousness upstream from the road crossing. Results indicate that the downstream water-quality constituent EMCs are highly correlated with the percentage of impervious area upstream.

Oregon↗
USGS2012

Parking lot runoff quality and treatment efficiencies of a hydrodynamic-settling device in Madison, Wisconsin, 2005-6

A hydrodynamic-settling device was installed in 2004 to treat stormwater runoff from a roof and parking lot located at the Water Utility Administration Building in Madison, Wis. The U.S. Geological Survey, in cooperation with the Wisconsin Department of Natural Resources, the City of Madison, cities in the Waukesha Permit Group, Hydro International, Earth Tech, Inc., National Sanitation Foundation International, and the U.S. Environmental Protection Agency, monitored the device from November 2005 through September 2006 to evaluate it as part of the U.S. Environmental Protection Agency's Environmental Technology Verification Program. Twenty-three runoff events monitored for flow volume and water quality at the device's inlet and outlet were used to calculate the percentage of pollutant reduction for the device. The geometric mean concentrations of suspended sediment (SS), "adjusted" total suspended solids (TSS), total phosphorus (TP), dissolved phosphorus (DP), total recoverable zinc (TZn), and total recoverable copper (TCu) measured at the inlet were 107 mg/L (milligrams per liter), 92 mg/L, 0.17 mg/L, 0.05 mg/L, 38 μg/L (micrograms per liter), and 12 μg/L, respectively, and these concentrations are in the range of values observed in stormwater runoff from other parking lots in Wisconsin and Michigan. Efficiency of the settling device was calculated using the efficiency ratio and summation of loads (SOL) methods. Using the efficiency ratio method, the device reduced concentrations of SS, and DP, by 19, and 15, percent, respectively. Using the efficiency ratio method, the device increased "adjusted" TSS and TZn concentrations by 5 and 19, respectively. Bypass occurred for 3 of the 23 runoff events used in this assessment, and the bypass flow and water-quality concentrations were used to determine the efficiency of the bypass system. Concentrations of SS, "adjusted" TSS, and DP were reduced for the system by 18, 5, and 18, respectively; however, TZn increased by 5 percent. Some of the TSS concentrations were "adjusted" to add the particles that remained on the sieves during sample processing. The loads of SS, "adjusted" TSS, and DP were reduced using the SOL method for the settling device by 38, 9, and 19 percent, respectively, and TZn increased by 13 percent. For the bypass system, the loads of SS, "adjusted" TSS, and DP had percentage reductions of 39, 12, 22, respectively, however TZn increased by 4 percent. The SOL method produced percentage reductions for SS and 'adjusted" TSS that were twice those for the efficiency ratio method. Removing the two large runoff events on August 23 and 24, 2006, from the SOL calculation brought the reduction for SS down to 16 and increased "adjusted" TSS by 4 percent. The two large runoff events were anomalies in that the runoff volumes and dissolved solids concentrations were greatly increased by overflow from an adjacent recycling facility. The SOL method was used to determine the percentage of SS load reduction for six different particle sizes for both the settling device and bypass system. Essentially no load reduction was observed for particles less than 125 micrometers (μm) in diameter, and about a 90-percent reduction occurred for particle sizes greater than 250 μm in diameter. The large removal efficiencies for particle sizes greater than 250 μm in diameter were further supported by the fact that more than 80 percent of the particle sizes trapped in the sump were greater than 250 μm in diameter. These results support the claim by the manufacturer of achieving a large percentage load reduction for particle sizes greater than 250 μm in diameter.

Wisconsin↗
USGS2002

Topographic map of Mars M 25M RKN

NOTES ON BASE This map is based on data from the Mars Orbiter Laser Altimeter (MOLA) (Smith and others, 2001), an instrument on NASA's Mars Global Surveyor (MGS) spacecraft (Albee and others, 2001). The image used for the base of this map represents more than 600 million measurements gathered between 1999 and 2001, adjusted for consistency (Neumann and others, 2001, 2002) and converted to planetary radii. These have been converted to elevations above the areoid as determined from a martian gravity field solution GMM2 (Lemoine and others, 2001), truncated to degree and order 50, and oriented according to current standards (see below). The average accuracy of each point is originally ~100 meters in horizontal position, and ~1 meter in radius (Neumann, 2001). However, the total elevation uncertainty is at least +3 m due to the global error in the areoid (+1.8 meters according to Lemoine and others, 2001) and regional uncertainties in its shape (communication from Neumann, 2002). The measurements were converted into a digital elevation model (DEM) (communication from Neumann, 2002; Neumann and others, 2001; Smith and others, 2001) using Generic Mapping Tools software (Wessel and Smith, 1998), with a resolution of 0.015625 degrees per pixel or 64 pixels per degree. In projection, the pixels are 926.17 meters in size at the equator. PROJECTION The Mercator projection is used between latitudes +57°, with a central meridian at 0° and latitude equal to the nominal scale at 0°. The Polar Stereographic projection is used for the polar regions north of the +55° parallel and south of the –55° parallel with a central meridian set for both at 0°. The adopted equatorial radius is 3,396.19 km (Duxbury and others 2002; Seidelmann and others 2002). COORDINATE SYSTEM The MOLA data were initially referenced to an internally consistent inertial coordinate system, derived from tracking of the MGS spacecraft. By adopting appropriate values for the orientation of Mars as defined by the International Astronomical Union (IAU) and the International Association of Geodesy (IAG) (Seidelmann and others, 2002), these inertial coordinates were converted into the planet-fixed coordinates (longitude and latitude) used on this map. These values include the orientation of the north pole of Mars (including the effects of precession), the rotation rate of Mars, and a value for W0 of 176.630°, where W0 is the angle along the equator to the east, between the 0° meridian and the equator's intersection with the celestial equator at the standard epoch J2000.0 (Seidelmann and others, 2002). This value of W0 was chosen (Duxbury and others, 2002) in order to place the 0° meridian through the center of the small (~500 m) crater Airy-0, located in the crater Airy (de Vaucouleurs and others, 1973; Seidelmann and others, 2002). Longitude increases to the east, and latitude is planetocentric as allowed by IAU/IAG standards (Seidelmann and others, 2002) and in accordance with current NASA and USGS standards (Duxbury and others, 2002). A secondary grid (printed in red) has been added to the map as a reference to the west longitude/planetographic latitude system that is also allowed by IAU/IAG standards (Seidelmann and others, 2002) and has also been used for Mars. The figure adopted to compute this secondary grid is an oblate spheroid with an equatorial radius of 3,396.19 km and a polar radius of 3,376.2 km (Duxbury and others, 2002; Seidelmann and others, 2002). MAPPING TECHNIQUES To create the topographic base image, the original DEM produced by the MOLA team in Simple Cylindrical projection with a resolution of 64 pixels per degree was projected into the Mercator and Polar Stereographic pieces. A shaded relief was generated from each DEM with a sun angle of 30° from horizontal and a sun azimuth of 270°, as measured clockwise from north, and a vertical exaggeration of 100%. Illumination is from the west, which follows a long-standing USGS tradition for planetary maps. This allows for continuity in the shading between maps and quadrangles, and most closely resembles lighting conditions found on imagery. The DEM values were then mapped to a smooth global color look-up table. Note that the chosen color scheme simply represents elevation changes and is not intended to imply anything about surface characteristics (e.g., past or current presence of water or ice). These two files were then merged and scaled to 1:25 million for the Mercator portion and 1:15,196,708 for the two Polar Stereographic portions, with a resolution of 300 dots per inch. The projections have a common scale of 1:13,923,113 at +56° latitude. NOMENCLATURE Names on this sheet are approved by the IAU and have been applied for features clearly visible at the scale of this map. For a complete list of the IAU-approved nomenclature for Mars, see the Gazetteer of Planetary Nomenclature at http://planetarynames.wr.usgs.gov/. Font color was chosen for readability. Names followed by an asterisk are provisionally approved. *M 25M RKN: Abbreviation for Mars, 1:25,000,000 series, shaded relief (R), with color (K) and nomenclature (N) (Greeley and Batson, 1990)

Open-File Report↗
USGS2007

Map and Database of Probable and Possible Quaternary Faults in Afghanistan

The U.S. Geological Survey (USGS) with support from the U.S. Agency for International Development (USAID) mission in Afghanistan, has prepared a digital map showing the distribution of probable and suspected Quaternary faults in Afghanistan. This map is a key component of a broader effort to assess and map the country's seismic hazards. Our analyses of remote-sensing imagery reveal a complex array of tectonic features that we interpret to be probable and possible active faults within the country and in the surrounding border region. In our compilation, we have mapped previously recognized active faults in greater detail, and have categorized individual features based on their geomorphic expression. We assigned mapped features to eight newly defined domains, each of which contains features that appear to have similar styles of deformation. The styles of deformation associated with each domain provide insight into the kinematics of the modern tectonism, and define a tectonic framework that helps constrain deformational models of the Alpine-Himalayan orogenic belt. The modern fault movements, deformation, and earthquakes in Afghanistan are driven by the collision between the northward-moving Indian subcontinent and Eurasia. The patterns of probable and possible Quaternary faults generally show that much of the modern tectonic activity is related to transfer of plate-boundary deformation across the country. The left-lateral, strike-slip Chaman fault in southeastern Afghanistan probably has the highest slip rate of any fault in the country; to the north, this slip is distributed onto several fault systems. At the southern margin of the Kabul block, the style of faulting changes from mainly strike-slip motion associated with the boundary between the Indian and Eurasian plates, to transpressional and transtensional faulting. North and northeast of the Kabul block, we recognized a complex pattern of potentially active strike-slip, thrust, and normal faults that form a conjugate shear system in a transpressional region of the Trans-Himalayan orogenic belt. The general patterns and orientations of faults and the styles of deformation that we interpret from the imagery are consistent with the styles of faulting determined from focal mechanisms of historical earthquakes. Northwest-trending strike-slip fault zones are cut and displaced by younger, southeast-verging thrust faults; these relations define the interaction between northwest-southeast-oriented contraction and northwest-directed extrusion in the western Himalaya, Pamir, and Hindu Kush regions. Transpression extends into north-central Afghanistan where north-verging contraction along the east-west-trending Alburz-Marmul fault system interacts with northwest-trending strike-slip faults. Pressure ridges related to thrust faulting and extensional basins bounded by normal faults are located at major stepovers in these northwest-trending strike-slip systems. In contrast, young faulting in central and western Afghanistan indicates that the deformation is dominated by extension where strike-slip fault zones transition into regions of normal faults. In addition to these initial observations, our digital map and database provide a foundation that can be expanded, complemented, and modified as future investigations provide more detailed information about the location, characteristics, and history of movement on Quaternary faults in Afghanistan.

Open-File Report↗
USGS2016

Flood-inundation maps for Lake Champlain in Vermont and in northern Clinton County, New York

Digital flood-inundation maps for an approximately100-mile length of Lake Champlain in Addison, Chittenden, Franklin, and Grand Isle Counties in Vermont and northern Clinton County in New York were created by the U.S. Geological Survey (USGS) in cooperation with the International Joint Commission (IJC). The flood-inundationmaps, which can be accessed through the International Joint Commission (IJC) Web site at http://www.ijc.org/en_/ , depict estimates of the areal extent flooding correspondingto selected water levels (stages) at the USGS lake gage on the Richelieu River (Lake Champlain) at Rouses Point, N.Y. (station number 04295000). In this study, wind and seiche effects (standing oscillating wave with a long wavelength) were not taken into account and the flood-inundation mapsreflect 11 stages (elevations) for Lake Champlain that are static for the study length of the lake. Near-real-time stages at this lake gage, and others on Lake Champlain, may be obtained on the Internet from the USGS National Water Information System at http://waterdata.usgs.gov/ or the National Weather Service Advanced Hydrologic Prediction Service at http:/water.weather.gov/ahps/ , which also forecasts flood hydrographs at the Richelieu River (Lake Champlain) at Rouses Point. Static flood boundary extents were determined for LakeChamplain in Addison, Chittenden, Franklin, and Grand Isle Counties in Vermont and northern Clinton County in New York using recently acquired (2013–2014) lidar (light detection and ranging) and may be referenced to any of the five USGS lake gages on Lake Champlain. Of these five lakgages, USGS lake gage 04295000, Richelieu River (Lake Champlain) at Rouses Point, N.Y., is the only USGS lake gage that is also a National Weather Service prediction location. Flood boundary extents for the Lake Champlain static flood-inundation map corresponding to the May 201 flood(103.2 feet [ft], National Geodetic Vertical Datum [NGVD] 29) were evaluated by comparing these boundary extents against the inundation area extents determined for the May 2011 flood (which incorporated documented high-water marksfrom the flood of May 201) (Bjerklie and others, 2014). A digital elevation model (DEM) was created by USGS, within a geographic information system (GIS), from the recently flown and processed light detection and ranging(lidar) data (2013–2014) in Vermont and the lake shore area of northern Clinton County in New York. The lidar data have a vertical accuracy of 0.3 to 0.6-ft (9.6 to 18.0-centimeters [cm]) and a horizontal resolution of 2.3 to 4.6 ft (0.7 to 1.4 meters). This DEM was used in determining the floodboundary for 11 flood stages at 0.5-ft intervals from 100.0 to104.0 ft (NGVD 29) and 1-ft intervals from 104.0 to 106.0 ft (NGVD 29) as referenced to the USGS lake gage 04295000, Richelieu River (Lake Champlain) at Rouses Point, N.Y. In addition, the May 2011 flood-inundation area for elevation103.20 ft (NGVD 29) (102.77 ft, North American Vertical Datum [NAVD] 88) was determined from this DEM. The May 2011 flood is the highest recorded lake water level (stage)at the Rouses Point, N.Y., lake gage. Flood stages greater than 101.5 ft (NGVD 29) exceed the “major flood stage”as defined by the NationalWeather Service for USGS lake gage 04295000. The availability of these maps, along with Internet information regarding current stage from the USGS lake gage and forecasted high-flow stages from the NationalWeather Service, will provide emergency management personnel and residents with information that is critical for flood responseactivities such as evacuations and road closures, as well as for post-flood recovery eforts.

New York, Vermont↗
USGS2013

Sediment distribution and hydrologic conditions of the Potomac aquifer in Virginia and parts of Maryland and North Carolina

Sediments of the heavily used Potomac aquifer broadly contrast across major structural features of the Atlantic Coastal Plain Physiographic Province in eastern Virginia and adjacent parts of Maryland and North Carolina. Thicknesses and relative dominance of the highly interbedded fluvial sediments vary regionally. Vertical intervals in boreholes of coarse-grained sediment commonly targeted for completion of water-supply wells are thickest and most widespread across the central and southern parts of the Virginia Coastal Plain. Designated as the Norfolk arch depositional subarea, the entire sediment thickness here functions hydraulically as a single interconnected aquifer. By contrast, coarse-grained sediment intervals are thinner and less widespread across the northern part of the Virginia Coastal Plain and into southern Maryland, designated as the Salisbury embayment depositional subarea. Fine-grained intervals that are generally avoided for completion of water-supply wells are increasingly thick and widespread northward. Fine-grained intervals collectively as thick as several hundred feet comprise two continuous confining units that hydraulically separate three vertically spaced subaquifers. The subaquifers are continuous northward but merge southward into the single undivided Potomac aquifer. Lastly, far southeastern Virginia and northeastern North Carolina are designated as the Albemarle embayment depositional subarea, where both coarse- and fine-grained intervals are of only moderate thickness. The entire sediment thickness functions hydraulically as a single interconnected aquifer. A substantial hydrologic separation from overlying aquifers is imposed by the upper Cenomanian confining unit. Potomac aquifer sediments were deposited by a fluvial depositional complex spanning the Virginia Coastal Plain approximately 100 to 145 million years ago. Westward, persistently uplifted granite and gneiss source rocks sustained a supply of coarse-grained sand and gravel. Immature, high-gradient braided streams deposited longitudinal bars and channel fills across the Norfolk arch subarea. By contrast, across the Salisbury and Albemarle embayment subareas, mature, medium- to low-gradient meandering streams deposited medium- to coarse-grained channel fills and point bars segregated from fine-grained overbank deposits. The Virginia depositional complex merged northward across the Salisbury embayment subarea with another complex in Maryland. Here, additional sediments were received from schist source rocks that underwent three cycles of initial uplift and rapid erosion followed by crustal stability and erosional leveling. Because of the predominance of coarse-grained sediments, transmissivity, hydraulic conductivity, and regional velocities of lateral flow through the Potomac aquifer are greatest across the Norfolk arch depositional subarea, but decrease progressively northward with increasingly fine-grained sediments. Confining units hydraulically separate the Potomac aquifer from overlying aquifers, as indicated by large vertical hydraulic gradients. By contrast, most of the Potomac aquifer internally functions hydraulically as a single interconnected aquifer, as indicated by uniformly small vertical gradients. Most fine-grained sediments within the aquifer do not hydraulically separate overlying and underlying coarse-grained sediments. Across the Salisbury embayment depositional subarea, however, hydraulic separation among the vertically spaced subaquifers is imposed by the intervening confining units. The Potomac aquifer is the largest and most heavily used source of groundwater in the Virginia Coastal Plain. Water-level declines as great as 200 feet create the potential for saltwater intrusion. Conventional stratigraphic correlation has been generally ineffective at accurately characterizing complexly distributed fluvial sediments that compose the Potomac aquifer. Consequently, the aquifer’s internal hydraulic connectivity and overall hydrologic function have not been well understood. Water-supply planning and development efforts have been hampered, and interpretations of regulatory criteria for allowable water-level declines have been ambiguous. An investigation undertaken during 2010–11 by the U.S. Geological Survey, in cooperation with the Virginia Department of Environmental Quality, provides a comprehensive regional description of the spatial distribution of Potomac aquifer sediments and their relation to hydrologic conditions. Altitudes and thicknesses of 2,725 vertical sediment intervals represent the spatial distribution of Potomac aquifer sediments in the Virginia Coastal Plain and adjacent parts of Maryland and North Carolina. Sediment intervals are designated as either dominantly coarse or fine grained and were determined by interpretation of geophysical logs and ancillary information from 456 boreholes. Sediment-interval and borehole summary statistical data indicate regional trends in sediment lithology and stratigraphic continuity, upon which three structurally based and hydrologically distinct sediment depositional subareas are designated. Broad patterns of sediment deposition over time are inferred from published sediment pollen-age data. Discrepancies in previously drawn hydrostratigraphic relations between southeastern Virginia and northeastern North Carolina are partly resolved based on borehole geophysical logs and a recently documented geologic map and corehole. A conceptual model theorizes the depositional history of the sediments and geologically accounts for their distribution. Documented pumping tests of the Potomac aquifer at 197 locations produced 336 values of transmissivity and 127 values of storativity. Based on effective aquifer thicknesses, 296 values of sediment hydraulic conductivity and 113 values of sediment specific storage are calculated. Vertical hydraulic gradients are calculated from 9,479 pairs of water levels measured between November 17, 1953, and October 4, 2011, in 129 closely spaced pairs of wells. Borehole sediment-interval and related data provide a means to achieve high yielding production wells in the Potomac aquifer by site-specific targeting of drilling operations toward water-bearing coarse-grained sand and gravel. Advance knowledge of the potential of different parts of the aquifer also aids in planning optimal groundwater-development areas. Depositional subareas further provide a possible context for resource management. Current (2013) regulatory limits on water-level declines are relative to top surfaces of subdivided upper, middle, and lower Potomac aquifers across the entire Virginia Coastal Plain, but have the potential to exceed the same limit relative to a single undivided Potomac aquifer. By contrast, designation of the sediments as a single aquifer in the Norfolk arch and Albemarle embayment subareas—and as a series of vertically spaced subaquifers and intervening confining units in the Salisbury embayment subarea—best reflects understanding of the Potomac aquifer and can avoid the potential for excessive water-level declines. Simulation modeling to evaluate effects of groundwater withdrawals could be designed similarly, including vertical discretization and (or) zonation of the Potomac aquifer based on depositional subareas and a geostatistical distribution of aquifer properties derived from borehole sediment-interval data. Further resource-management information needs extend beyond the developed part of the Potomac aquifer, particularly across the Northern Neck and Middle Peninsula where only the shallowest part of the aquifer is known, and include structural aspects such as faults, basement bedrock, and the Chesapeake Bay impact crater.

Maryland, North Carolina, Virginia↗
USGS2007

Geophysical analysis of the Salmon Peak Formation near Amistad Reservoir Dam, Val Verde County, Texas, and Coahuila, Mexico, March 2006, to aid in piezometer placement

Since 1992, numerous sinkholes have developed northwest of the Amistad Reservoir dam on the Rio Grande. Increases in the discharge of springs south of the dam, on the western side of the Rio Grande, in Coahuila, Mexico, have been documented. In 1995 the Mexico Section of the International Boundary and Water Commission (IBWC) completed a study of the western embankment (Coahuila, Mexico) of the dam that included surface geophysics, borehole geophysics, and installation of piezometers to learn more about subsurface conditions. As part of a 5-year safety inspection in 2005, technical advisors recommended that one line of similarly constructed piezometers be installed on the eastern embankment (Val Verde County, Texas) of the dam for comparison of water levels (potentiometric head) on both the western and eastern embankments of Amistad Reservoir dam. To provide technical assistance for the horizontal and vertical placement of piezometers on the eastern embankment of Amistad Reservoir dam, the U.S. Geological Survey, in cooperation with the U.S. Section of the IBWC, conducted a study along both the western and eastern embankments of Amistad Reservoir dam. The study involved an integrated approach using surface and borehole geophysical methods. In the western embankment investigation, geological and geophysical characteristics that indicate relatively large water-yielding properties of the Salmon Peak Formation were identified. The direct-current (DC) resistivity method was selected as the surface geophysical reconnaissance technique to correlate relatively large water-yielding properties of the Salmon Peak Formation, identified from analysis of borehole geophysical logs, with variations in subsurface resistivity. The dipole-dipole array and the reciprocal Schlumberger array were selected as the most applicable DC resistivity arrays. Two resistivity units were identified in both the dipole-dipole array data and the reciprocal Schlumberger array data along DC resistivity profiles on both embankments. Resistivity unit 1 generally is of relatively low resistivity, ranging from 45 to 150 ohm-meters compared with resistivity unit 2, which ranges from 120 to 345 ohm-meters (depending on the DC array type). The presence of mapped sinkholes in the reservoir north of the western embankment study area and the zone of increased water content (as indicated by zones of low neutron log count rates in nearby piezometers) leads to the conclusion that resistivity unit 1 is a preferential flow path where surface water from Amistad Reservoir is forced into the ground-water system (because of increased head from the reservoir). In the eastern embankment investigation, trends in the spatial distribution of sinkholes and the occurrence of weathered zones were identified from geologic descriptions of cores. The correlation of surface geophysical DC resistivity, historical lithologic data, and general trend of documented sinkholes along the eastern end of the eastern embankment profile were used to justify further exploration (drilling of piezometers) in the eastern expression of resistivity unit 1. The spatial location of the piezometers and the screened intervals were selected to best match the locations of the screened intervals of the western embankment piezometers. Six piezometers were installed on the eastern embankment and logged using borehole geophysical techniques. Surface DC resistivity sections superimposed on the resistivity logs for two piezometers indicate three discernible resistivity units that correlate with resistivity units 2, 1, and 2, respectively, identified in the western embankment study area. Resistivity units 1 and 2 in the DC resistivity profiles generally correspond with low and high resistivity zones, respectively, on the normal and lateral resistivity logs collected in the nearby piezometers at the time of installation.

Scientific Investigations Report↗
USGS2017

Fluorine

Fluorine compounds are essential in numerous chemical and manufacturing processes. Fluorspar is the commercial name for fluorite (isometric CaF 2 ), which is the only fluorine mineral that is mined on a large scale. Fluorspar is used directly as a fluxing material and as an additive in different manufacturing processes. It is the source of fluorine in the production of hydrogen fluoride or hydrofluoric acid, which is used as the feedstock for numerous organic and inorganic chemical compounds. The United States was the world’s leading producer of fluorspar until the mid-1950s. In the mid-1970s, the U.S. fluorspar mining industry began to decline because of foreign competition. By 1982, there was essentially only a single U.S. producer left, and that company ceased mining in 1996. Consumption of fluorspar in the United States peaked in the early 1970s, which was also the peak period of U.S. steel production. Since then, U.S. fluorspar consumption has decreased substantially; the United States has nonetheless increased its imports of downstream fluorine compounds, such as, in order of tonnage imported, hydrofluoric acid, aluminum fluoride, and cryolite. This combination of no U.S. production (until recently) and high levels of consumption has made the United States the world’s leading fluorspar-importing country, in all its various forms. The number of fluorspar-exporting countries has decreased substantially in recent decades, and, as a result, the United States has become dependent on just a few countries to supply its needs. In 2013, the United States imported the majority of its fluorspar from three countries, which were, in descending order of the amount imported, Mexico, China, and South Africa. Geologically, in igneous systems, fluorine is one of a number of elements that are “incompatible.” These incompatible elements become concentrated in the residual magma while the common silicates crystallize upon magma ascent and cooling, leading to relatively high fluorine concentrations in the more evolved or differentiated igneous rocks and in hydrothermal deposits associated with those evolved igneous rocks. In sedimentary rocks, fluorine’s highest concentrations are found in phosphorites because fluorine substitutes for hydroxyl ions in apatite, which leads to fluorine concentrations of, typically, from 2 to 4 weight percent in phosphorites. Because of the presence of fluorine, phosphate fertilizer manufacturers can produce a fluorosilicic acid byproduct. Most deposits mined for fluorine are hydrothermal, however, and consist of fluorine minerals that precipitated from hot water. Magmatic brines and brines from deep within sedimentary basins that have high concentrations of dissolved fluoride are the mineralizing fluids for various types of hydrothermal fluorspar deposits. Relatively dilute hydrothermal fluids that formed in some volcanic rocks can also transport sufficient fluoride to form a high-grade fluorspar deposit. Fluorite has low solubility in a common range of hydrothermal temperatures, particularly from about 160 degrees Celsius (°C) down to 60 °C. The increasing fluorite solubility below 60 °C partly explains why some water with exceptionally high levels of dissolved fluorine are found even at ambient temperatures in evaporitic lake basins in some East African Rift valleys in Kenya and Tanzania. The geologic conditions that led to the high concentrations there are known to exist in a number of other places in the world as well, including, perhaps, places in the Basin and Range province of the United States. Eight minerals or mineral groups have sufficient fluorine in their structures to be considered as possible ores of the element; they are bastnaesite (also spelled bastnäsite; and other fluorocarbonates), cryolite, sellaite, villiaumite, fluorite, fluorapatite (in phosphorites), various phyllosilicates, and topaz. Fluorite is currently the only mineral that is mined for fluorine, and nomineral except fluorite is likely to become a source of commercially produced fluorine as a primary product as long as supplies from relatively thick and high-grade fluorite deposits continue to be available. At least seven classes (which include one subclass) of hydrothermal fluorite deposits are recognized; they are classified according to their tectonic and (or) magmatic settings, as follows: (1) carbonatite-related fluorspar deposits; (2) alkaline-intrusion-related fluorspar deposits; (3) alkaline-volcanic-related epithermal fluorspar deposits; (4) Mississippi Valley-type fluorspar deposits (and a subclass of salt-related carbonate-hosted fluorspar deposits); (5) fluorspar deposits related to strongly differentiated granites; (6) subalkaline-volcanic-related epithermal fluospar deposits; and (7) fluorspar deposits that appear to be conformable within tuffaceous limy lacustrine sediments. An eighth class (not hydrothermal) is that of fluorspar deposits concentrated in soils and weathered zones; that is, residual fluorspar deposits. Generally, fluorspar deposits related to strongly differentiated granites have larger tonnages and lower grades than carbonatite-related fluorspar deposits, which, in turn, have larger tonnages and lower grades than fluorspar vein deposits from various other classes. The United States has a few identified resources of fluorspar, most notably the Klondike II property in the Illinois- Kentucky fluorspar district located about 8 kilometers southwest of Salem, Kentucky, which has a large vein that contains at least 1.6 million metric tons at a grade of 60 percent CaF 2 (Feytis, 2009). Additional fluorspar resources of lower grade but larger tonnage have been identified at Hicks Dome in the Illinois-Kentucky fluorspar district and at Lost River near the western tip of the Seward Peninsula in Alaska, along with a couple of dozen smaller, higher grade resources. Internationally, new mines that either opened before the beginning of 2013 or were scheduled to open soon after that time include the Nui Phao tungsten-fluorspar-bismuth-copper-gold deposit in northern Vietnam; the St. Lawrence project in Newfoundland, Canada, which is located in a well-known fluorspar district; the Bamianshan deposit, which is related to a strongly differentiated granite in northwestern Zhejiang Province, China, near some of that Province’s large, subalkaline-volcanic-related epithermal veins; and the Nokeng project in South Africa, which is also related to a strongly differentiated granite. Other deposits in northwestern Australia, Nevada (United States), Norway, South Africa, and Sweden have been identified and could be put into production within just a few years. Among undiscovered resources, an interesting possibility might be to produce a fluorine product from evaporitic, high-fluorine, high-pH sodium-carbonate brines like Lake Magadi (Kenya) and Lake Natron (Tanzania) in Africa’s Eastern Rift Valley. In addition, apparently conformable fluorspar deposits in tuffaceous limy lacustrine sediments, such as those in Italy, are likely to occur in similar young alkalic volcanic settings elsewhere in the world. Modern geophysical and geochemical exploration techniques have typically not been brought to bear in exploration for new fluorspar deposits, although such techniques are likely to be used in future exploration. The tendency for fluorine to dissolve in significant concentrations in water at low temperature allows both surface water and groundwater to be used as sampling media in geochemical exploration. Evolved granite-related fluorspar deposits may be particularly susceptible to geophysical exploration methods because crystalline rocks that form a basement to sedimentary sections can be approximately defined with gravity and magnetic methods, and magnetite-bearing skarns can be directly detected with magnetic surveys. Environmental considerations of fluorine mining focus especially on drinking water, where high fluorine concentrations can lead to tooth decay; dental and skeletal fluorosis; and bone and cartilage conditions, including genu valgum, which is the crippling bone deformity more commonly known as knock knee. Trace amounts of other elements in fluorspar ores are a concern at some deposits; for example, high beryllium concentrations in alkaline-volcanic-related epithermal deposits or high cadmium concentrations associated with Mississippi Valley-type and salt-related carbonate-hosted fluorspar deposits. Future research might include testing whether fluorine can be extracted economically from high-pH, sodium-carbonate brines and exploring for new occurrences of apparently conformable fluorspar deposits in tuffaceous limy lacustrine sediments outside of the Latium Region of Italy. Other promising new areas of research could be studies of fluorspar deposit fluid inclusion compositions by quadrupole mass spectrometry, by noble gas mass spectrometry on irradiated fluid inclusions, or by chlorine isotopes, while also measuring the chemistry of the same fluid inclusions either by bulk crush-and-leach methods or by laser ablation-inductively coupled plasma mass spectrometry. Advanced studies of fluid inclusion chemistry could be applied beneficially to some of the enigmatic large epithermal fluorspar veins at various places in the world, where they might determine those deposits’ possible relationships to igneous intrusions, or to dissolved salt, or to heated meteoric water in volcanic sections, or perhaps to all three. This knowledge could help focus new exploration.

Professional Paper↗
USGS1993

SPECtrum Processing Routines User's Manual Version 3 (program SPECPR)

Introduction The SPECtrum Processing Routines (SPECPR) is a large-scale interactive program for general one-dimensional array processing and optimized for reflectance spectroscopy data and analysis. The program processes one-dimensional arrays up to 4852 data points and the operations include addition, subtraction, multiplication, division, trigonometric functions, logarithmic and exponential functions, and many more specialized routines. The Specpr user interface is a menu driven, character command system with all user input entered as ascii characters to the program. User commands are thoroughly checked for the context in which they will be used. This provides for essentially all user input mistakes to be caught and appropriate error messages to be issued. The user interface also allows for command aliasing, variable substitutions, command history and batch command processing to occur at any point in the program. The user can also control his or her own variables, increment and decrement them, and monitor certain internal variables. This flexibility allows for sophisticated programming, e.g. from specific computations, to management of a database. All arrays are treated as one dimensional lists of numbers, each with an appropriate header for identification and history. Any array can be plotted versus another array, and each axis labeled appropriately. For example, you may plot a reflectance spectrum versus wavelength or wavelength versus reflectance. The array type is simply a label, so data like temperature versus time can be treated just as easily as a reflectance spectrum. Some routines are specific to certain purposes, however. For example, the Planck black body generator is specific to intensity versus wavelength, whereas a smoothing routine could operate equally on reflectance versus wavelength or temperature versus time. Each array element can have an error bar associated with it and errors are propagated through all appropriate routines. Data points can be marked deleted and deleted points are tracked appropriately. These features allow the program to be effectively used on real world data. This philosophy has resulted in a flexible system for which a user can manipulate data arrays efficiently. The program was designed for analysis of laboratory, field, telescopic, and spacecraft spectroscopic data, and although general in nature, there are biases built in. For example, horizontal axis labels default to wavelength in micrometers, but the defaults can be changed. Because of the orientation in Specpr to process planetary and terrestrial data, and the need to access such data, Specpr has been extended to have access to non-Specpr default file types. In particular, Specpr can "skewer" an 3-dimensional (3D) data file along any of the three axes. Further, block skewers can be done with the standard deviation of the mean computed for each channel in the block. Currently, Specpr can access any of the standard file types common in the terrestrial and planetary remote sensing communities. Details on 3D file I/O are given in Chapter 6. Specpr has multiple record types available within a single Specpr data file. Currently defined record types are data and text. In a data record, a standard one-dimensional array is held, along with its header information. The header information includes a title, history, dates and time of data acquisition and when the data were last processed, the user who processed the data, information typical of a spectrum like temperature and viewing geometry, as well as pointers to wavelengths, resolution, and text. If a spectrum is more than 256 channels in length, then the data gets put in succeeding records in the file, and the following records are continuation records. The second record type is the text, where a title and a block of text (up to 19 kbytes) can be stored. This record type is typically used for a description of samples, experiments, instruments, or data processing for a particular data set. It could also store the actual commands used to create a spectrum. The commands for a complicated plot for a publication might be stored in a text record. That way, you can easily regenerate the plot, or even write the commands to a non-Specpr file, modify them with any editor and re-execute them. The file types, text and command processing, combined with the math and special functions, provide for very powerful and general analysis tool. These facilities also provide the framework for database management.

Open-File Report↗
USGS2011

A promising tool for subsurface permafrost mapping-An application of airborne geophysics from the Yukon River Basin, Alaska

Permafrost is a predominant physical feature of the Earth's Arctic and Subarctic clines and a major consideration encompassing ecosystem structure to infrastructure engineering and placement. Perennially frozen ground is estimated to cover about 85 percent of the state of Alaska where northern reaches are underlain with continuous permafrost and parts of interior Alaska are underlain by areas of discontinuous and (or) sporadic permafrost (fig. 1). The region of Interior Alaska, where permafrost is scattered among unfrozen ground, is a complex mosaic of terrains and habitats. Such diversity creates arrays of lakes and surface-water and groundwater patterns that continental populations of migratory waterfowl and internationally significant fisheries have adapted to over time. A road or pipeline might pass over frozen and unfrozen ground, affecting the types of materials and engineering approaches needed to sustain the infrastructure. Effective mapping of discontinuous permafrost at scales meaningful ecologically and (or) from an engineering perspective has been a long-standing challenge. Using techniques such as borehole logging for site-specific assessments or botanical techniques that can suggest underlying permafrost distributions can be labor intensive and difficult to accomplish at the scale and remoteness of much of Alaska. The climate is changing in the Arctic and Subarctic regions. The warming observed throughout much of Alaska could create widespread changes in permafrost. How the warming of the permafrost affects near-surface processes, ecosystems, and community infrastructure and ecosystems is not clear. A better understanding of the dynamic distribution and physical properties of permafrost, from continuous to discontinuous, will provide knowledge of how the permafrost environment may change in the future and help inform engineering and natural resource response strategies. Here we discuss an application of an airborne remote sensing methodology for mapping and shade imaging permafrost characteristics at various scales. This work provides the first look into three-dimensional distribution of permafrost in the areas around Fort Yukon and is a demonstration of the application of airborne electromagnetic (AEM) to permafrost mapping. Such an approach is attractive, particularly in Arctic and Subarctic studies, where ground access is difficult and ecosystems are fragile. In June 2010, the U.S. Geological Survey (USGS) conducted an AEM survey near Fort Yukon, Alaska. The primary focus of this survey was to map the distribution of permafrost in selected areas in order to supply information for the development of groundwater models of the Yukon River Basin. However, the methodologies have more far-reaching ecological and engineering applications. Approximately 1,800 line kilometers were acquired in a combination of typical block style surveying in the immediate area of Fort Yukon and in long reconnaissance lines over a broader area. The widely spaced lines were flown to cross the modern Yukon River in "X" like patterns with intersections at features that have been previously studied (fig. 2). AEM is used to gather data on the electrical resistivity of materials in the subsurface below the flight path of the helicopter, which are then analyzed to interpret the subsurface lithology and the location and extent of permafrost. For this survey, the electrical resistivity was imaged to depths on the order of 50-100 meters. Images from the survey can be qualitatively compared with known permafrost features and suggest new permafrost features. Electrical properties of earth materials are affected by lithology as well as temperature and the presence of ice; frozen materials become substantially more resistive. This allows for the identification of permafrost from the resistivity image (Abraham and others, 2011). In the area of Fort Yukon, the AEM survey shows elevated resistivities extending to depth, likely indicative of thick permafrost. This depth corresponds well to observations from a borehole drilled in the area in the late 1990s, which detected permafrost to a depth of about 100 meters (Clark and others, 2009). In contrast to the area of Fort Yukon, the Yukon River and its floodplain are not associated with deep resistive sediments, suggesting a lack of deep permafrost, at least within the depth range of the AEM mapping (fig. 3).

Yukon River Basin↗
USGS2006

The Virginia Coastal Plain hydrogeologic framework

A refined descriptive hydrogeologic framework of the Coastal Plain of eastern Virginia provides a new perspective on the regional ground-water system by incorporating recent understanding gained by discovery of the Chesapeake Bay impact crater and determination of other geological relations. The seaward-thickening wedge of extensive, eastward-dipping strata of largely unconsolidated sediments is classified into a series of 19 hydrogeologic units, based on interpretations of geophysical logs and allied descriptions and analyses from a regional network of 403 boreholes. Potomac aquifer sediments of Early Cretaceous age form the primary ground-water supply resource. The Potomac aquifer is designated as a single aquifer because the fine-grained interbeds, which are spatially highly variable and inherently discontinuous, are not sufficiently dense across a continuous expanse to act as regional barriers to ground-water flow. Part of the Potomac aquifer in the outer part of the Chesapeake Bay impact crater consists of megablock beds, which are relatively undeformed internally but are bounded by widely separated faults. The Potomac aquifer is entirely truncated across the inner part of the crater. The Potomac confining zone approximates a transition from the Potomac aquifer to overlying hydrogeologic units. New or revised designations of sediments of Late Cretaceous age that are present only south of the James River include the upper Cenomanian confining unit, the Virginia Beach aquifer and confining zone, and the Peedee aquifer and confining zone. The Virginia Beach aquifer is a locally important ground-water supply resource. Sediments of late Paleocene to early Eocene age that compose the Aquia aquifer and overlying Nanjemoy-Marlboro confining unit are truncated along the margin of the Chesapeake Bay impact crater. Sediments of late Eocene age compose three newly designated confining units within the crater, which are from bottom to top, the impact-generated Exmore clast and Exmore matrix confining units, and the Chickahominy confining unit. Piney Point aquifer sediments of early Eocene to middle Miocene age overlie most of the Chesapeake Bay impact crater and beyond, but are a locally significant ground-water supply resource only outside of the crater across the middle reaches of the Northern Neck, Middle, and York-James Peninsulas. Sediments of middle Miocene to late Miocene age that compose the Calvert confining unit and overlying Saint Marys confining unit effectively separate the underlying Piney Point aquifer and deeper aquifers from overlying shallow aquifers. Saint Marys aquifer sediments of late Miocene age separate the Calvert and Saint Marys confining units across two limited areas only. Sediments of the Yorktown-Eastover aquifer of late Miocene to late Pliocene age form the second most heavily used ground-water supply resource. The Yorktown confining zone approximates a transition to the overlying late Pliocene to Holocene sediments of the surficial aquifer, which extends across the entire land surface in the Virginia Coastal Plain and is a moderately used supply. The Yorktown-Eastover aquifer and the eastern part of the surficial aquifer are closely associated across complex and extensive hydraulic connections and jointly compose a shallow, generally semiconfined ground-water system that is hydraulically separated from the deeper system. Vertical faults extend from the basement upward through most of the hydrogeologic units but may be more widespread and ubiquitous than recognized herein, because areas of sparse boreholes do not provide adequate spatial control. Hydraulic conductivity probably is decreased locally by disruption of depositional intergranular structure by fault movement in the generally incompetent sediments. Localized fluid flow in open fractures may be unique in the Chickahominy confining unit. Some hydrogeologic units are partly to wholly truncated where displacements are large relative to unit thickness, resulting in lateral flow barriers or flow conduits. The tops of the Saint Marys confining unit, YorktownEastover aquifer, and Yorktown confining zone are widely sculpted by erosion that reflects both the present-day topography and buried paleochannels. Fault displacements across the top surfaces of these hydrogeologic units probably have been beveled by erosion. Additionally, erosion has modified the margins of many hydrogeologic units by truncation along the valleys of major rivers and their tributaries, beneath which underlying hydrogeologic units are incised. As a result, the surficial aquifer is in contact with a “patchwork” of underlying hydrogeologic units that create a complex array of hydraulic connections between the confined and unconfined ground-water systems.

Virginia↗
USGS1984

High-temperature, large-volume, lavalike ash-flow tuffs without calderas in southwestern Idaho

Rhyolitic rocks were erupted from vents in and adjacent to the Owyhee Mountains and Owyhee Plateau of southwestern Idaho from 16 m.y. ago to about 10 m.y. ago. They were deposited on a highly irregular surface developed on a variety of basement rocks that include granitic rocks of Cretaceous age, quartz latite and rhyodacite tuffs and lava flows of Eocene age, andesitic and basaltic lava flows of Oligocene age, and latitic and basaltic lava flows of early Miocene age. The rhyolitic rocks are principally welded tuffs that, regardless of their source, have one feature in common-namely internal characteristics indicating en-masse, viscous lavalike flowage. The flowage features commonly include considerable thicknesses of flow breccia at the bases of various cooling units. On the basis of the tabular nature of the rhyolitic deposits, their broad areal extents, and the local preservation of pyroclastic textures at the bases, tops, and distal ends of some of the deposits, we have concluded that the rocks were emplaced as ash flows at extremely high temperatures and that they coalesced to liquids before final emplacement and cooling. Temperatures of l090?C and higher are indicated by iron-titanium oxide compositions. Rhyolites that are about 16 m.y. old are preserved mostly in the downdropped eastern and western flanks of the Silver City Range and they are inferred to have been erupted from the Silver City Range. They rarely contain more than about 2 percent phenocrysts that consist of quartz and subequal amounts of plagioclase and alkali feldspar; commonly, they contain biotite, and they are the only rhyolitic rocks in the area to do so. The several rhyolitic units that are 14 m.y. to about 10 m.y. old contain only pyroxene-principally ferriferous and intermediate pigeonites-as mafic constituents. The rhyolites of the Silver City Range comprise many cooling units, none of which can be traced for great distances. Rocks erupted from the Owyhee Plateau include two sequences that were traced over areas having diameters of about 100 km. These two sheets are the herein-named Swisher Mountain Tuff, which is about 13.8 m.y. old, and the Little Jacks Tuff, which is about 10 m.y. old. The Swisher Mountain Tuff was erupted from the Juniper Mountain volcanic center, a gentle dome that is not bounded by arcuate faults indicative of cauldron subsidence. The tuff is 200 m thick over a considerable area in and adjacent to its source. It apparently thins gradually toward its distal edges, and it is inferred to be uniformly distributed around its source at Juniper Mountain. The unit contains vitrophyres at various intervals from base to top, and, although the vitrophyres are, in general, flow layered and commonly flow brecciated, they occasionally contain well-defined pumice clasts. The vitrophyres indicate compound cooling, and, near the distal edges of the sheet, some of them probably represent complete cooling breaks. The Little Jacks Tuff onlaps the Swisher Mountain Tuff in expo sures east of Juniper Mountain, and it is inferred to have been erupted from a source on the part of the Owyhee Plateau that lies just east of the area studied. This inferred source area, like that at Juniper Mountain, is also expressed today as a gentle dome without structural features indicative of cauldron subsidence. The Little Jacks Tuff, in most exposures in the deep canyons of the Plateau, consists of at least four cooling units, and, in places in the eastern part of the studied area near the source area, it possibly comprises as many as six. Although there is no obvious evidence of erosion between the various cooling units, magnetic polarity measurements indicate that there were at least two magnetic reversals during the eruption interval of the Little Jacks Tuff. Like the Swisher Mountain Tuff, the Little Jacks has flattened pumice clasts in a few outcrops-principally at the bases of the various cooling units. The two tuff sequences are calc-a

Professional Paper↗