USGS Science⌕ Search

SEARCH · USGS Science

Results for “Surveying and Land Information Systems”

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 721 records · Page 40Linked to original sources

Groundwater status and trends for the Columbia Plateau Regional Aquifer System, Washington, Oregon, and Idaho

Well information and groundwater-level measurements for the Columbia Plateau Regional Aquifer System in Washington, Oregon, and Idaho, were compiled from data provided by the U.S. Geological Survey and seven other organizations. From the full set of about 60,000 wells and 450,000 water-level measurements a subset of 761 wells within the aquifers of the Columbia River Basalt Group (CRBG) then was used to develop a simple linear groundwater-level trend map for 1968–2009. The mean of the trends was a decline of 1.9 feet per year (ft/yr), with 72 percent of the water levels in wells declining. Rates of declines greater than 1.0 ft/yr were measured in 50 percent of wells, declines greater than 2.0 ft/yr in 38 percent of wells, declines greater than 4.0 ft/yr in 29 percent of wells, and declines greater than 8.0 ft/yr in 4 percent of wells. Water-level data were used to identify groups of wells with similar hydraulic heads and temporal trends to delineate areas of overall similar groundwater conditions. Discontinuities in hydraulic head between well groups were used to help infer the presence of barriers to groundwater flow such as changes in lithology or the occurrence of folds and faults. In areas without flow barriers, dissimilarities in response of well groups over time resulted from the formation of groundwater mounds caused by recharge from irrigation or regions of decline caused by pumping. The areas of focus for this analysis included the Umatilla area, Oregon, and the Palouse Slope/eastern Yakima Fold Belt in the Columbia Basin Ground Water Management Area (GWMA) consisting of Adams, Franklin, Grant, and Lincoln Counties, Washington. In the Umatilla area, water levels from 286 wells were used to identify multiple areas of high hydraulic gradient that indicate vertical and horizontal barriers to groundwater flow. These barriers divide the groundwater-flow system into several compartments with varying degrees of interconnection. Horizontal flow barriers commonly correspond to mapped geologic structure and result in horizontal hydraulic gradients that progressively become steeper from north to south corresponding to an increase in structural complexity that may be impeding recharge from the uplands into the heavily developed areas. Most CRBG aquifers in the Umatilla area are declining and since 1970, cumulative declines range from about 100 to 300 feet. Significant vertical hydraulic gradients are documented for relatively small areas near Umatilla, and since the 1970s, downward vertical gradients in these areas have been increasing as hydraulic heads in the deeper units have declined. The absence of vertical gradients over much of the area may be a consequence of flow through commingling wells that results in the equilibration of the heads between aquifers. On the Palouse Slope in the central GWMA, large groundwater declines occurred during 1968–2009 along a north-south swath in the middle of the region. An analysis of 1,195 wells along major flow paths and through the area of persistent groundwater-level declines indicates that barriers to flow are not as evident in this area as in Umatilla. This is consistent with the geologic interpretation of the Palouse Slope as being a gently folded structure created by voluminous sheet flows of CRBG lavas. Groundwater discharge into the sediment-filled coulees, where the upper aquifers are intersected at land surface by incised canyons, is proposed as an alternative to explain local steepening of the hydraulic gradient along the Palouse Slope previously attributed to the presence of a groundwater dam. Comparison of generalized potentiometric surface maps developed for pre-development conditions and post-2000 conditions indicate that pre-development groundwater flow was from the uplands toward the Columbia and Snake River and that post-2000 flow patterns in the area are controlled by irrigation practices that have resulted in broad regions of elevated or depressed hydraulic head. In some cases, irrigation-related changes in head have reversed groundwater flow directions. Evidence of significant vertical hydraulic gradients exists, although much of the aquifer thickness is affected by commingling of wells. The effect of commingling and its relative contribution to problems related to groundwater-level declines remains unclear.

Washington, Oregon, Idaho↗

Hydrologic Data for Deep Creek Lake and Selected Tributaries, Garrett County, Maryland, 2007-08

Introduction Recent and ongoing efforts to develop the land in the area around Deep Creek Lake, Garrett County, Maryland, are expected to change the volume of sediment moving toward and into the lake, as well as impact the water quality of the lake and its many tributaries. With increased development, there is an associated increased demand for groundwater and surface-water withdrawals, as well as boat access. Proposed dredging of the lake bottom to improve boat access has raised concerns about the adverse environmental effects such activities would have on the lake. The Maryland Department of Natural Resources (MDDNR) and the U.S. Geological Survey (USGS) entered into a cooperative study during 2007 and 2008 to address these issues. This study was designed to address several objectives to support MDDNR?s management strategy for Deep Creek Lake. The objectives of this study were to: Determine the current physical shape of the lake through bathymetric surveys; Initiate flow and sediment monitoring of selected tributaries to characterize the stream discharge and sediment load of lake inflows; Determine sedimentation rates using isotope analysis of sediment cores; Characterize the degree of hydraulic connection between the lake and adjacent aquifer systems; and Develop an estimate of water use around Deep Creek Lake. Summary of Activities Data were collected in Deep Creek Lake and in selected tributaries from September 2007 through September 2008. The methods of investigation are presented here and all data have been archived according to USGS policy for future use. The material presented in this report is intended to provide resource managers and policy makers with a broad understanding of the bathymetry, surface water, sedimentation rates, groundwater, and water use in the study area. The report is structured so that the reader can access each topic separately using any hypertext markup (HTML) language reader. In order to establish a base-line water-depth map of Deep Creek Lake, a bathymetric survey of the lake bottom was conducted in 2007. The data collected were used to generate a bathymetric map depicting depth to the lake bottom from a full pool elevation of 2,462 feet (National Geodetic Vertical Datum of 1929). Data were collected along about 90 linear miles across the lake using a fathometer and a differentially corrected global positioning system. As part of a long-term monitoring plan for all surface-water inputs to the lake, streamflow data were collected continuously at two stations constructed on Poland Run and Cherry Creek. The sites were selected to represent areas of the watershed under active development and areas that are relatively stable with respect to development. Twelve months of discharge data are provided for both streams. In addition, five water-quality parameters were collected continuously at the Poland Run station including pH, specific conductance, temperature, dissolved oxygen, and turbidity. Water samples collected at Poland Run were analyzed for sediment concentration, and the results of this analysis were used to estimate the annual sediment load into Deep Creek Lake from Poland Run. To determine sedimentation rates, cores of lake-bottom sediments were collected at 23 locations. Five of the cores were analyzed using a radiometric-dating method, allowing average rates of sedimentation to be estimated for the time periods 1925 to 2008, 1925 to 1963, and 1963 to 2008. Particle-size data from seven cores collected at locations throughout the study area were analyzed to provide information on the amount of fine material in lake-bed sediments. Groundwater levels were monitored continuously in four wells and weekly in nine additional wells during October, November, and December of 2008. Water levels were compared to recorded lake levels and precipitation during the same period to determine the effect of lake-level drawdown and recovery on the adjacent aquifer systems. Water use in the Deep Creek Lake wa

Open-File Report↗

Metallogeny of the Great Basin: Crustal evolution, fluid flow, and ore deposits

The Great Basin physiographic province in the Western United States contains a diverse assortment of world-class ore deposits. It currently (2006) is the world’s second leading producer of gold, contains large silver and base metal (Cu, Zn, Pb, Mo, W) deposits, a variety of other important metallic (Fe, Ni, Be, REE’s, Hg, PGE) and industrial mineral (diatomite, barite, perlite, kaolinite, gallium) resources, as well as petroleum and geothermal energy resources. Ore deposits are most numerous and largest in size in linear mineral belts with complex geology. U.S. Geological Survey (USGS) scientists are in the final year of a research project initiated in the fall of 2001 to increase understanding of relations between crustal evolution, fluid flow, and ore deposits in the Great Basin. Because of its substantial past and current mineral production, this region has been the focus of numerous investigations over the past century and is the site of ongoing research by industry, academia, and state agencies. A variety of geoinformatic tools was used to organize, reinterpret, and display, in space and time, the large amounts of geologic, geophysical, geochemical, and hydrologic information deemed pertinent to this problem. This information, in combination with concentrated research on (1) critical aspects of the geologic history, (2) an area in northern Nevada that encompasses the major mineral belts, and (3) important mining districts and deposits, is producing new insights about the interplay between key tectonic events, hydrothermal fluid flow, and ore genesis in mineral belts. The results suggest that the Archean to Holocene history of the Great Basin was punctuated by several tectonic events that caused fluids of different origins (sea water, basinal brine, meteoric water, metamorphic water, magmatic water) to move through the crust. Basement faults reactivated during these events localized deformation, sedimentation, magmatism, and hydrothermal fluid flow in overlying rocks to form mineral belts that contain ore deposits of different types and ages that are locally superimposed (demonstrating inheritance). Fluid flow in these systems also was influenced by the distribution of permeable lithologies and paleotopographic highs and lows. Hydrothermal fluids evolved from their initial chemistries towards compositions that reflect the ƒ O 2 and ƒ S 2 buffering capacity of, and the ligands and metals present in, the rocks (±older mineralization) through which they moved. In northern Nevada, where gold deposits are relatively common, carbonaceous, pyritic strata buffered fluids of diverse origins to H 2 S-rich compositions so they could transport gold repeatedly over Paleozoic-Cenozoic time (convergent evolution). Ore formed where metal-laden fluids encountered effective physicochemical traps. Maps of Neogene basin fill and erosion surfaces identify areas where preexisting ore deposits have been progressively exposed or concealed. Comparisons with analogous terrains and deposit types in other parts of the world provide global context. The initial findings and some of the databases, geologic maps, sections, reconstructions, hydrogeologic models, topical syntheses, regional overviews, short courses, field guides, and deposit comparisons produced by project staff and associated managers, contractors, and collaborators have been presented in numerous abstracts, symposia, USGS publications, and professional journals over the last 5 years (see the extensive bibliography). Notable among these was the 2005 Geological Society of Nevada symposium in Reno, Nevada, and the 2005 Geological Society of America annual meeting in Salt Lake City, Utah, where project results were presented to audiences from around the nation and world. The final results of the project will be submitted for publication in 2007 to appropriate USGS and professional journals. A special issue of GEOSPHERE, scheduled for publication in 2007, will be devoted to the results of this project and related work. This special issue will reach an international audience and be available worldwide on the internet. Much of the research for this project has concentrated on areas that will receive the focused attention of the mining industry in the future. As such, the data and interpretations generated by this project have direct use for land-use managers in Federal, State, and local agencies. Improved hydrogeologic models developed by this project will considerably enhance ongoing and future water resource investigations in the region. The increased understanding of when, where, and how hydrothermal systems produce significant economic deposits has direct uses for mineral exploration and for future USGS mineral resource assessments in the Great Basin and other parts of the world.

Great Basin↗

Insights and strategic opportunities from the USGS 2024 Per- and Polyfluoroalkyl Substances (PFAS) Interagency Workshop

Introduction In 2021, the U.S. Geological Survey (USGS) published Circular 1490 titled, “Integrated Science for the Study of Perfluoroalkyl and Polyfluoroalkyl Substances (PFAS) in the Environment: A Strategic Science Vision for the U.S. Geological Survey” (Tokranov and others, 2021). Circular 1490 was created to be a resource for USGS scientists prioritizing and planning research related to per- and polyfluoroalkyl substances (PFAS) and to be a guide for developing partnerships with other scientists, State and Federal agencies, and stakeholders engaged in PFAS research and management and mitigation of the environmental and human-health effects of PFAS. This USGS PFAS Strategic Science Vision document was intended to be the foundation for a “living strategic vision,” periodically providing updates on the state of USGS PFAS research, emerging PFAS data gaps and needs, and progress on interagency and stakeholder PFAS partnerships and priorities. To meet this objective, the USGS planned to host an Interagency and Stakeholder PFAS Workshop every 2–3 years. During September 10–12, 2024, the USGS hosted the first Interagency and Stakeholder PFAS Workshop in Reston, Virginia. The Workshop brought together experts from other Federal agencies (U.S. Environmental Protection Agency, National Institute of Environmental Health Sciences, Food and Drug Administration, Department of Defense [Air Force, Army]), State agencies (Washington Fish and Wildlife, Virginia Department of Transportation), and academia (Harvard University, University of Maryland) to address key challenges relating to the measurement and modeling of PFAS and the implications for environmental health. Participants engaged in in-depth discussions centered around six pivotal topics related to PFAS: (1) sampling protocols, methods and interpretation; (2) environmental sources, source apportionment, and occurrence; (3) environmental fate and transport; (4) human and wildlife exposure routes and risk; (5) bioconcentration, bioaccumulation, and biomagnification; and (6) ecotoxicology and effects. Each topic had three breakout sessions. A recurrent theme of workshop discussions was how data on a nationwide scale for PFAS occurrence in various environmental matrices, including air, water, food crops, biota, soil, and streambed sediment could help to advance scientific understanding. Participants noted significant geospatial data gaps, particularly in the midwestern and southern United States and the Pacific Northwest. PFAS data collection tends to be more robust along the eastern seaboard and in California. Participants stressed how enhancing the integration of large and small datasets across various agencies could help to support national scale understanding of PFAS. To address these gaps, attendees suggested leveraging datasets from Federal entities like the USGS and the U.S. Department of Defense, State agencies, and municipal utility services to develop predictive contaminant detection and transport models. Improved coordination between water quality programs and USGS research could help to facilitate access to valuable data, leading to comprehensive databases that inform PFAS point (wastewater treatment plants and landfills) and nonpoint (runoff from land, atmospheric deposition, food packaging) sources, environmental transport mechanisms, environmental detection and concentrations, potential exposure routes, and health effects on different biota, including humans. A specific request was made to develop a map demarking the depth of modern (1953 or later) groundwater, which is susceptible to surface-derived anthropogenic (that is, human-made) contamination, based on tritium-age dating. Emphasis was placed on incorporation of hydrology, groundwater flow paths, groundwater–surface water interactions, and landscape factors in predictive statistical models as a step to improve contaminant source identification and tracking. Molecular fingerprinting approaches garnered attention as techniques to link specific PFAS mixtures detected in a sample to environmental sources and levels in biota (Dávila-Santiago and others, 2022). Integrating data from abiotic (that is, water, soil, and air) and biotic (that is, living organisms) systems identified as a research opportunity. For example, understanding the composition of soils and sediments, which include a mixture of mineral, plant, and animal components, could advance understanding of exposure pathways. The discussions highlighted opportunities to explore and understand the potential redistribution and biotic exposures of PFAS from biosolid and wastewater treatment plant effluent land application practices, in addition to atmospheric releases and discharges from landfill and wastewater treatment plants. Participants identified research gaps surrounding how these sources may contribute to contamination and may affect surrounding ecosystems, including a better definition of anthropogenic background concentrations. Moving forward, the collection of co-occurrence data was noted as a means to improve understanding of complex mixtures and to leverage companion modeling efforts focused on areas with high and low contamination levels to identify areas of concern and unaffected resources. Participants emphasized how centralized USGS databases and the establishment of sample-metadata archives can help to ensure that samples are preserved and accessible for future research. In conclusion, the workshop participants identified opportunities to bridge data gaps and improve measurement techniques, modeling frameworks, databases, and communication, to enhance the understanding of PFAS and their effects on environmental and human health. Upon completion of the workshop, participants indicated an interest in developing strategic data collection, modeling, and analytical approaches to address these challenges.

Open-File Report↗

Geochemistry and microbiology of groundwater and solids from extraction and monitoring wells and their relation to well efficiency at a Federally operated confined disposal facility, East Chicago, Indiana

In cooperation with the U.S. Army Corps of Engineers, Chicago District, the U.S. Geological Survey investigated the processes affecting water quality, geochemistry, and microbiology in representative extraction and monitoring wells at a confined disposal facility (CDF) in East Chicago, Indiana. The CDF is a 140-acre Federally-managed facility that was the former location of an oil refinery and is now used for the long-term disposal and storage of dredge material from the Indiana Harbor and Indiana Harbor Canal. Residual petroleum hydrocarbons and leachate from the CDF are contained within the facility by use of a groundwater cutoff wall. The wall consists of a soil-bentonite slurry and a gradient control system made up of an automated network of 96 extraction wells, 42 monitoring wells, and 2 ultrasonic sensors that maintain an inward hydraulic gradient at the site. The pumps in the extraction wells require vigilant maintenance and must be replaced when unable to withdraw water at a rate sufficient to maintain the required inward gradient. The wells are screened in the Calumet aquifer, a coarse-grained sand and gravel unit that extends approximately 35 feet below the land surface and is not utilized for drinking-water supply at the CDF or in the surrounding area. This study was initiated to identify the cause of decreased pump discharges and to identify potential mitigation strategies. For this study, the U.S. Geological Survey collected groundwater and solids from monitoring and extraction wells. Groundwater samples were collected during June 2014 for precautionary health screening and on four occasions during September 2014 through November 2014. Groundwater samples collected from two extraction wells during June 2014 were analyzed for concentrations of anthropogenic organic constituents. During September through November 2014, groundwater samples were collected from one additional extraction well, and samples from three monitoring wells were analyzed for concentrations of inorganic and organic constituents, dissolved gases, and bacterial abundance and diversity. Solid samples were collected during April 2014, during September 2014 through November 2014, and during November 2016. Solid samples were collected from the exterior of extraction-well pumps and as flocculent from water samples. Solid samples were collected from 10 wells, including 1 extraction well and 3 monitoring wells sampled for water quality. Solid samples were analyzed for mineralogy, solid-phase habit, geochemistry, and organic composition. The following is a list of observations that were made during this study: (1) the water quality is substantially variable among the six well locations sampled as part of this study—lower (more negative) redox values and higher concentrations of many constituents (including calcium, magnesium, sodium, and sulfate) and properties (including dissolved solids, hardness, and turbidity) were detected in sampled wells located near the extraction wells with the highest frequency of failure; (2) water-level drawdown is variable between extraction wells—wells with the greatest drawdown may pull deeper groundwater into the borehole; (3) dissolved gas results indicate reducing oxidation-reduction processes in the aquifer material that can feasibly contribute iron, carbon dioxide, and other byproducts from hydrocarbon degradation to precipitates and solids that accumulate on and impair pump operation; (4) crystalline and amorphous solid-phase minerals are precipitating in the borehole; (5) several types of bacteria are present in water pumped from extraction wells and are likely responsible for bonding mineral and microbiologic matter to the pump (and other well components); and (6) bacteria may create microenvironments that facilitate precipitation of solids or inhibit dissolution of unstable minerals once the bacteria adhere to biofilm attached to the pump. Results of the study indicate that bacteria may be accumulating and entrapping solid material on the exterior of pumps. This accumulation reduces heat transfer and water discharge from the pump and may lead to decreased efficiency or mechanical failure. Observations could not be made on the well screen, gravel pack, or surrounding geologic formation; therefore, mitigating measures in the borehole may not solve well-productivity issues. Remedies for the pump fouling problems were derived from the review and interpretation of data collected during this study and from information documented in other sources about groundwater well fouling. Potential remedies to problems associated with pump fouling at the CDF may include the following: (1) reducing attractiveness of the extraction wells for microbiological growth by modifying the chemical or physical environment of the well, (2) modifying the pump exterior to decrease microbiological adherence, (3) changing the pumping regime to control the chemistry of water entering the well from the surrounding aquifer material, (4) modifying the pumps to be less physically and thermally attractive, and (5) removing hydrocarbons from groundwater and the aquifer material surrounding the wells or adding surfactants to make them more mobile. Pilot scale testing may be necessary to identify the most effective treatment or combination of treatments.

Indiana↗

Groundwater/surface-water interactions in the Partridge River Basin and evaluation of hypothetical future mine pits, Minnesota

The Partridge River Basin (PRB) covers 156 square miles in northeastern Minnesota with headwaters in the Mesabi Iron Range. The basin is characterized by extensive wetlands, lakes, and streams in poorly drained and often thin glacial material overlying Proterozoic bedrock. To better understand the interaction between these extensive surface water features and the groundwater system, a three-dimensional, steady-state, groundwater-flow model of the PRB was developed by the U.S. Geological Survey in cooperation with the Great Lakes Indian Fish & Wildlife Commission using the finite-difference computer code MODFLOW-NWT. The model simulates steady-state base flow in streams and groundwater interactions using the streamflow routing (SFR2) package. Existing mining features including tailings basins, stockpiles, pumped mine pits, and flooded mine pits were simulated using either high hydraulic conductivity zones or the drain (DRN) package. The unsaturated zone flow (UZF) package was used to better represent the groundwater system in areas with a high water table and for wetlands often associated with such areas. UZF typically is used to represent unsaturated zone processes but also can simulate the rejection of recharge and groundwater discharge to the land surface when the water table is near land surface. The steady-state model used data from the 2011 to 2013 period when 2011 high-resolution land surface (light detecting and ranging [lidar]) data were available that reflected land-surface and water elevations from mining activity in the basin. The parameter-estimation software suite PEST_HP was used to obtain a best fit of the modeled to measured groundwater levels, streamflow, pit inflow rates, and mapped peat deposits. The PEST calibration used the target residuals from two models with the same model parameters and targets from two separate periods: (1) a 1995–2015 calibration model, which provided a larger number of calibration targets, and (2) a 2011–2013 mining conditions model, which included calibration targets that reflected conditions consistent with the modeled mine-workings topography. Calibration of the PRB model resulted in ranges of glacial horizontal hydraulic conductivity parameters that generally agreed with literature values and other models of the region. Horizontal hydraulic conductivity of the bedrock was higher in the upper bedrock layers where numerous and continuous fractures have been observed and lower in the deeper bedrock layers. Average basin-wide calibrated infiltration was 5.3 inches per year. An average of 4.6 inches per year of infiltration crosses the water table and becomes recharge and 0.7 inch per year is rejected by UZF due to saturated conditions at the land surface. Simulated groundwater runoff (the sum of rejected recharge and groundwater seepage to the land surface) can either be routed to streams or removed from the model as evapotranspiration. The calibrated model indicates relatively shallow groundwater-flow paths dominating and approximately 50 percent of the stream base flow coming from groundwater runoff. The 2011–2013 mining conditions model was then used to develop five model scenarios simulating the response of the groundwater and surface-water system to potential hydrologic stress. The purpose of these mine pit scenarios is to present a possible workflow to quantify a model’s uncertainty for a given model forecast and serve as a possible guide for initial data collection that may improve a future model’s ability to make such a forecast. The scenarios included one scenario with the currently existing Peter Mitchell pit at final buildout and flooded to an elevation of 1,500 feet, and four scenarios with a hypothetical, new mine pit plus the flooded Peter Mitchell at final buildout. The five model scenarios were used to forecast streamflow at six locations in the PRB, pit inflow rates for the new mine pits and the flooded Peter Mitchell pit, and the average depth to water in 12 wetlands. A linear uncertainty analysis was performed using information from the PEST calibration and tools in the PyEMU python package to assess model uncertainty propagation to the model forecasts. Streamflows generally were reduced with future mining and the greatest streamflow reductions occurred from the flooded Peter Mitchell Pit, probably due to its large size. Average depth to groundwater in wetlands was most affected the closer the wetland was to a new mine pit. Linear uncertainty methods were also used to evaluate data worth, which is the ability for potential new groundwater elevation observations to reduce the uncertainty in scenario forecasts. Data worth was performed for a grid of new hydraulic head observations. Overall, areas with nonnegligible data worth generally corresponded to wetland areas with no groundwater seepage to land surface from UZF. These model behaviors indicated that the land-surface boundary condition simulated by the UZF package was pinning the groundwater elevations to the land surface in areas with groundwater seepage (33 percent of the 2011–2013 base conditions model) such that the sensitivity to new observations in these areas was minimal. Therefore, representing wetlands as boundary conditions minimized the usefulness of data worth calculations because wetland areas were present over a large part of the model domain. Probabilistic capture zones were estimated for each of the mines in the model scenarios. A capture zone represents the area contributing recharge to a model feature, like a well or a mine pit, and can be calculated by forward tracking particles from the water table. By using Monte Carlo techniques, it is possible to generate estimated capture zones that include the probability of recharge capture given the uncertainty present in the model. Monte Carlo techniques use randomly generated model parameter sets sampled from a plausible parameter range to create many possible realizations. The resulting capture zone arrays were calculated by tallying the total number of realizations in which a particle from a model cell was captured by the feature. Probabilities from the Monte Carlo runs ranged from 1 (captured in 100 percent of the runs) near the pits to 0 (captured in 0 percent of the runs) at the edges of the capture zone. Capture zones were not always spatially continuous; for example, the capture zone for the proposed mine pits south of the flooded Peter Mitchell pit was discontinuous with capture surrounding the proposed mine pit and north of the flooded Peter Mitchell pit. This northern section represents deeper groundwater flow paths that originate in the topographic high, move under the flooded pit, and discharge into the proposed pit. This pattern of capture indicates the possibility of some deeper flow through the upper fractured bedrock when the shallow groundwater flow system is modified. These results underscore that future site-specific applications of the base condition model require the input of site-specific data and recalibration to focus on the site of interest.

Minnesota↗

Initial characterization of the groundwater system near the Lower Colorado Water Supply Project, Imperial Valley, California

In 2009, the U.S. Geological Survey, in cooperation with the city of Needles, began a study of the hydrogeology along the All-American Canal, which conveys water from the Colorado River to the Imperial Valley. The focus of this study was to gain a better understanding of the effect of lining the All-American Canal, and other management actions, on future total dissolved solids concentrations in groundwater pumped by Lower Colorado Water Supply Project wells that is delivered to the All-American Canal. The study included the compilation and evaluation of previously published hydrogeologic and geochemical information, establishment of a groundwater-elevation and groundwater-quality monitoring network, results of monitoring groundwater elevations and groundwater quality from 2009 to 2011, site-specific hydrologic investigations of the Lower Colorado Water Supply Project area, examination of groundwater salinity by depth by using time-domain electromagnetic surveys, and monitoring of groundwater-storage change by using microgravity methods. Prior to the completion of the All-American Canal in 1940, groundwater in the study area flowed from east to west, and groundwater was recharged primarily by underflow from the Colorado River Valley. After construction of the All-American Canal, groundwater elevations were altered in the study area as seepage of Colorado River water from the All-American Canal and other canals became the dominant recharge source. By 2005, groundwater elevations had increased by as much as 50–70 feet along the All-American Canal. Superimposed on the east-to-west groundwater gradient was groundwater movement away from the All-American Canal to the north and, most likely, to the south into Mexico. After lining the All-American Canal, from 2007 to 2010, groundwater elevations declined as seepage from the All-American Canal decreased. Between 2005 (the last complete groundwater-elevation survey prior to lining the All-American Canal) and 2011, groundwater elevations declined 20–40 feet along the All-American Canal and as much as 40–45 feet in the vicinity of Lower Colorado Water Supply Project pumping wells. Water-quality and isotope data were used to differentiate historically recharged groundwater from groundwater more recently recharged by seepage of Colorado River surface water from the All-American Canal. Prior to the completion of the All-American Canal in 1940, groundwater in the southern part of the study area was primarily sodium-chloride/sulfate type water that had relatively low total dissolved solids concentrations (500–820 milligrams per liter). During 2007–11, groundwater in the southern part of the study area, near the All-American Canal, ranged from sodium-chloride type water to mixed-cation-sulfate type water that had total dissolved solids concentrations generally less than 879 milligrams per liter. The stable-isotopic signature of groundwater near the All-American Canal sampled in 2009–11 indicated inputs of Colorado River water that had been affected by evaporation, and radioactive isotopes indicated that a substantial fraction of water had been recharged recently, within the past 60 years. This contrasted with historically recharged groundwater near the All-American Canal, which had higher sodium and chloride concentrations, and lower calcium and sulfate concentrations, than recent recharge from the All-American Canal. Groundwater at a distance from the All-American Canal, in the East Mesa, Algodones Dunes, Pilot Knob Mesa, and Cargo Muchacho Mountains piedmont, was found to have higher total dissolved solids concentrations (generally greater than 1,000 milligrams per liter) than recently recharged groundwater near the All-American Canal. Time-domain electromagnetic data indicated that low-salinity groundwater was present down to about 377 feet below land surface near the All-American Canal; groundwater salinity at depth increased with distance north from the All-American Canal. Groundwater several miles or more from the canal also did not contain tritium and had a residence time on the order of thousands to tens of thousands of years. The groundwater in the piedmont of the Cargo Muchacho Mountains had a distinctly light stable-isotopic signature indicative of recharge by runoff from local precipitation, whereas the stable isotopic signature of groundwater in the East Mesa and the Algodones Dunes indicated a mixture of local precipitation and historic Colorado River recharge sources. During and after lining the All-American Canal (2007–11), groundwater elevations in the Lower Colorado Water Supply Project area declined, while total dissolved solids concentrations remained relatively constant. The total dissolved solids concentrations in well LCWSP-2 ranged from 650 to 800 milligrams per liter during this study. Depth-specific water-quality and isotope sampling at well LCWSP-2 indicated the groundwater pumped from the deeper part of the screened interval (240–280 feet below land surface) contained a greater proportion of historical groundwater than the groundwater pumped from the shallower part of the screened interval (350–385 feet below land surface). Age-tracer data at well LCWSP-2 indicated that all depths of the screened interval had received recent recharge from seepage of Colorado River water from the All-American Canal.

California↗

A brief history and summary of the effects of river engineering and dams on the Mississippi River system and delta

The U.S. Geological Survey Forecast Mekong project is providing technical assistance and information to aid management decisions and build science capacity of institutions in the Mekong River Basin. A component of this effort is to produce a synthesis of the effects of dams and other engineering structures on large-river hydrology, sediment transport, geomorphology, ecology, water quality, and deltaic systems. The Mississippi River Basin (MRB) of the United States was used as the backdrop and context for this synthesis because it is a continental scale river system with a total annual water discharge proportional to the Mekong River, has been highly engineered over the past two centuries, and the effects of engineering have been widely studied and documented by scientists and engineers. The MRB is controlled and regulated by dams and river-engineering structures. These modifications have resulted in multiple benefits including navigation, flood control, hydropower, bank stabilization, and recreation. Dams and other river-engineering structures in the MRB have afforded the United States substantial socioeconomic benefits; however, these benefits also have transformed the hydrologic, sediment transport, geomorphic, water-quality, and ecologic characteristics of the river and its delta. Large dams on the middle Missouri River have substantially reduced the magnitude of peak floods, increased base discharges, and reduced the overall variability of intraannual discharges. The extensive system of levees and wing dikes throughout the MRB, although providing protection from intermediate magnitude floods, have reduced overall channel capacity and increased flood stage by up to 4 meters for higher magnitude floods. Prior to major river engineering, the estimated average annual sediment yield of the Mississippi River Basin was approximately 400 million metric tons. The construction of large main-channel reservoirs on the Missouri and Arkansas Rivers, sedimentation in dike fields, and protection of channel banks by revetments throughout the basin, have reduced the overall sediment yield of the MRB by more than 60 percent. The primary alterations to channel morphology by dams and other engineering projects have been (1) channel simplification and reduced dynamism; (2) lowering of channel-bed elevation; and (3) disconnection of the river channel from the flood plain, except during extreme flood events. Freshwater discharge from the Mississippi River and its associated sediment and nutrient loads strongly influence the physical and biological components in the northern Gulf of Mexico. Ninety percent of the nitrogen load reaching the Gulf of Mexico is from nonpoint sources with about 60 percent coming from fertilizer and mineralized soil nitrogen. Much of the phosphorus is from animal manure from pasture and rangelands followed by fertilizer applied to corn and soybeans. Increased nutrient enrichment in the northern Gulf of Mexico has resulted in the degradation of water quality as more phytoplankton grow, which increases turbidity and depletes oxygen in the lower depths creating what is known as the "dead zone." In 2002, the dead zone was 22,000 square kilometers (km2), an area similar to the size of the State of Massachusetts. Changes in the flow regime from engineered structures have had direct and indirect effects on the fish communities. The navigation pools in the upper Mississippi River have aged, and these overwintering habitats, which were created when the pools filled, have declined as sedimentation reduces water depth. Reproduction of paddlefish may have been adversely affected by dams, which impede access to suitable spawning habitats. Fishes that inhabit swift-current habitats in the unimpounded lower Mississippi River have not declined as much as in the upper Mississippi River. The decline of the pallid sturgeon may be attributable to channelization of the Missouri River above St. Louis, Missouri. The Missouri River supports a rich fish community and remains relatively intact. Nevertheless, the widespread and long history of human intervention in river discharge has contributed to the declines of about 25 percent of the species. The Mississippi River Delta Plain is built from six delta complexes composed of a massive area of coastal wetlands that support the largest commercial fishery in the conterminous United States. Since the early 20th century, approximately 4,900 km2 of coastal lands have been lost in Louisiana. One of the primary mechanisms of wetland loss on the Plaquemines-Balize complex is believed to be the disconnection of the river distributary network from the delta plain by the massive system of levees on the delta top, which prevent overbank flooding and replenishment of the delta top by sediment and nutrient deliveries. Efforts by Federal and State agencies to conserve and restore the Mississippi River Delta Plain began over three decades ago and have accelerated over the past decade. Regardless of these efforts, however, land losses are expected to continue because the reduced upstream sediment supplies are not sufficient to keep up with the projected depositional space being created by the combined forces of delta plain subsidence and global sea-level rise.

Alabama;Alberta;Arkansas;Colorado;Georgia;Illinois↗

2023-2024 Coastal sage scrub and chaparral community monitoring for western San Diego County

Western San Diego County is dominated by shrublands supporting biologically diverse native plant and animal communities. Widespread urbanization has led to regional habitat loss and fragmentation, and many species in these shrubland communities are rare, threatened, or endangered. Large-scale, multiple-species conservation planning has resulted in a regional preserve system that focuses on these shrubland communities. Several large-scale threats are leading to type conversion from shrub-dominated to non-native invasive annual grass-dominated vegetation. To understand the changes that are occurring to native shrublands, we have developed a vegetation monitoring program with several components at multiple spatial scales, focused on quantifying coastal sage scrub (CSS) and chaparral vegetation community characteristics. Several drivers of change associated with type conversion of native shrubland to non-native annual grassland have been identified by previous research including increasing fire frequency, nitrogen deposition from air pollution, and prolonged and intense drought associated with changing climate. Loss of ecological integrity indices have been proposed as useful measures of the threat of degradation and type conversion of shrublands in San Diego County. For this study, ecological integrity is defined as a system’s ability to maintain species’ relationships and functions comparable to natural habitat in the region. Previous studies have identified the percent cover of invasive non-native annual grasses as a proxy for overall ecological degradation (loss of integrity) that is consistent across native taxonomic groups. Increased cover of non-native grass is associated with lower integrity of the shrubland vegetation community as shrub-associated plant and animal species are replaced by species preferring grassy and disturbed habitats. The objectives of this CSS and chaparral vegetation community monitoring plan are to: 1) Determine the distribution, composition, structure, and integrity of CSS and chaparral vegetation communities on conserved lands in western San Diego County, 2) Identify whether these attributes of the vegetation communities are changing over time, and 3) Evaluate relationships of known drivers of change (threats) and environmental factors in association with changes in vegetation community attributes. The CSS and chaparral vegetation community monitoring program is divided into four components: 1) vegetation mapping, 2) GIS/remote sensing office analysis of landscape-scale data, 3) permanent field plots using Unmanned Aerial Systems (UAS) and field data collections, and 4) animal and target species surveys and rapid assessment protocols. The first component, which is not detailed in this vegetation monitoring plan, aims to map vegetation communities every 10-15 years based on a classification developed for western San Diego County. High resolution aerial imagery will be used to update the 2012 vegetation map and expand it across the entire study area. The second component uses remote sensing models to annually track ecological integrity of shrublands across the study area and will include a map of areas of change and areas of stability over several decades. These landscape-scale integrity classifications will be used to analyze ecological processes, threats, and abiotic factors relative to changes to shrubland ecological integrity over time and space. The third component includes field surveys of 100 permanent plots across areas historically mapped as shrublands. Surveys in the 1930s mapped vegetation types using plot data. By using this historical classification map, we included areas that have already type-converted from shrubland to non-native annual grassland. The plots were split between CSS (55 plots) and chaparral (45 plots) and stratified into four geographical eco-subregions to guarantee coverage over small patches of habitat along the heavily developed coast. Surveys will include the collection of UAS imagery at a very high resolution. Species-level identifications will be made from the imagery based on a plant list compiled of all species detected in the plot during a thorough field survey by botanists, combined with geo-referenced samples of plant species locations. In addition, herbaceous cover will be estimated in the field using nine 1-m diameter circles (one per subplot) to obtain ocular estimates of cover for each plant species within the circle frame. Soil samples will be collected and analyzed for important element compositions. These data will be analyzed to evaluate plot-level ecological integrity based on species composition and cover. Repeated monitoring will allow evaluation of changes in vegetation attributes over time with known drivers or threats and other environmental factors. In addition, analyses will focus on indicator species and various measures of biological diversity for the vegetation communities. Finally, animal species and rare plants will be assessed using either taxa-based rapid assessment protocols or specialized species-specific protocols for rarer species. The purpose of these assessments is to document the status, habitat, and threat covariates of specific species and confirm the species composition and diversity of animal taxonomic groups (e.g., pollinators). Diversity and abundance of animal species at vegetation plots will be used to refine measures and thresholds of ecological integrity. Rapid assessment protocols for animal taxonomic groups can include multiple detection methods such as, camera traps, cover boards, and bird point counts. Pollinators will be monitored at plots using a protocol currently being developed in conjunction with but separately from this plan. Target rare plant species will be monitored using the regional Inspect and Manage (IMG) protocol that measures the status of rare plant occurrences and habitat and threat covariates over time. Species-specific animal survey methods will be refined as these species are prioritized for future monitoring. The goal of this monitoring program is to classify CSS and chaparral vegetation community integrity, identify areas of degradation across western San Diego County, and characterize drivers, and environmental factors associated with loss of ecological integrity. A combination of vegetation mapping, landscape-scale remote sensing, and field plots will be used to address all the aspects of our research questions. Data compiled and collected will be available to conservation partners to help inform future management decisions.

California↗

Optical satellite data volcano monitoring: a multi-sensor rapid response system

In this chapter, the use of satellite remote sensing to monitor active geological processes is described. Specifically, threats posed by volcanic eruptions are briefly outlined, and essential monitoring requirements are discussed. As an application example, a collaborative, multi-agency operational volcano monitoring system in the north Pacific is highlighted with a focus on the 2007 eruption of Kliuchevskoi volcano, Russia. The data from this system have been used since 2004 to detect the onset of volcanic activity, support the emergency response to large eruptions, and assess the volcanic products produced following the eruption. The overall utility of such integrative assessments is also summarized. The work described in this chapter was originally funded through two National Aeronautics and Space Administration (NASA) Earth System Science research grants that focused on the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) instrument. A skilled team of volcanologists, geologists, satellite tasking experts, satellite ground system experts, system engineers and software developers collaborated to accomplish the objectives. The first project, Automation of the ASTER Emergency Data Acquisition Protocol for Scientific Analysis, Disaster Monitoring, and Preparedness , established the original collaborative research and monitoring program between the University of Pittsburgh (UP), the Alaska Volcano Observatory (AVO), the NASA Land Processes Distributed Active Archive Center (LP DAAC) at the U.S. Geological Survey (USGS) Earth Resources Observation and Science (EROS) Center, and affiliates on the ASTER Science Team at the Jet Propulsion Laboratory (JPL) as well as associates at the Earth Remote Sensing Data Analysis Center (ERSDAC) in Japan. This grant, completed in 2008, also allowed for detailed volcanic analyses and data validation during three separate summer field campaigns to Kamchatka Russia. The second project, Expansion and synergistic use of the ASTER Urgent Request Protocol (URP) for natural disaster monitoring and scientific analysis , has expanded the project to other volcanoes around the world and is in progress through 2011. The focus on ASTER data is due to the suitability of the sensor for natural disaster monitoring and the availability of data. The instrument has several unique facets that make it especially attractive for volcanic observations (Ramsey and Dehn, 2004). Specifically, ASTER routinely collects data at night, it has the ability to generate digital elevation models using stereo imaging, it can collect data in various gain states to minimize data saturation, it has a cross-track pointing capability for faster targeting, and it collects data up to ±85° latitude for better global coverage. As with any optical imaging-based remote sensing, the viewing conditions can negatively impact the data quality. This impact varies across the optical and thermal infrared wavelengths as well as being a function of the specific atmospheric window within a given wavelength region. Water vapor and cloud formation can obscure surface data in the visible and near infrared (VNIR)/shortwave infrared (SWIR) region due mainly to non-selective scattering of the incident photons. In the longer wavelengths of the thermal infrared (TIR), scattering is less of an issue, but heavy cloud cover can still obscure the ground due to atmospheric absorption. Thin clouds can be optically-transparent in the VNIR and TIR regions, but can cause errors in the extracted surface reflectance or derived surface temperatures. In regions prone to heavy cloud cover, optical remote sensing can be improved through increased temporal resolution. As more images are acquired in a given time period the chances of a clear image improve dramatically. The Advanced Very High Resolution Radiometer (AVHRR) routine monitoring, which commonly collects 4-6 images per day of any north Pacific volcano, takes advantage of this fact. The rapid response program described in this chapter also improves the temporal resolution of the ASTER instrument. ASTER has been acquiring images of volcanic eruptions since soon after its launch in December 1999. An early example included the observations of the large pyroclastic flow deposit emplaced at Bezymianny volcano in Kamchatka, Russia. The first images in March 2000, just weeks after the eruption, revealed the extent, composition, and cooling history of this large deposit and of the active lava dome (Ramsey and Dehn, 2004). The initial results from these early datasets spurred interest in using ASTER data for expanded volcano monitoring in the north Pacific. It also gave rise to the multi-year NASA-funded programs of rapid response scheduling and imaging throughout the Aleutian, Kamchatka and Kurile arcs. Since the formal establishment of the programs, the data have provided detailed descriptions of the eruptions of Augustine, Bezymianny, Kliuchevskoi and Sheveluch volcanoes over the past nine years (Wessels et al., in press; Carter et al., 2007, 2008; Ramsey et al., 2008; Rose and Ramsey, 2009). The initial research focus of this rapid response program was specifically on automating the ASTER sensor’s ability for targeted observational scheduling using the expedited data system. This urgent request protocol is one of the unique characteristics of ASTER. It provides a limited number of emergency observations, typically at a much-improved temporal resolution and quicker turnaround with data processing in the United States rather than in Japan. This can speed the reception of the processed data by several days to a week. The ongoing multi-agency research and operational collaboration has been highly successful. AVO serves as the primary source for status information on volcanic activity, working closely with the National Weather Service (NWS), Federal Aviation Administration (FAA), military and other state and federal emergency services. Collaboration with the Russian Institute of Volcanology and Seismology (IVS)/Kamchatka Volcanic Eruption Response Team (KVERT) is also maintained. Once a volcano is identified as having increased thermal output, ASTER is automatically tasked and the volcano is targeted at the next available opportunity. After the data are acquired, scientists at all the agencies have access to the images, with the primary science analysis carried out at the University of Pittsburgh and AVO. Results are disseminated to the responsible monitoring agencies and the global community through e-mail mailing lists.

Book chapter↗

Assessment of water levels, nitrate, and arsenic in the Carson Valley Alluvial Aquifer and the development of a data visualization tool for the Carson River Basin, Nevada

Residents of Carson Valley, Douglas County, Nevada, rely on the basin-fill alluvial aquifer underlying the valley for drinking water. Since the 1980s, groundwater levels and water-quality data have been collected to monitor the status of the aquifer system and to assist in planning efforts to address current (2024) and future demand. The U.S. Geological Survey (USGS), in cooperation with Douglas County, Nevada, evaluated trends in water levels, nitrate, and arsenic concentrations from a network of monitoring and domestic wells in Carson Valley. This work also assessed the monitoring well network to determine the suitability of wells for characterizing the occurrence of arsenic in the groundwater. Monitoring of constituents, such as nitrate and arsenic concentrations, is needed to assess changes in contaminant distribution and to evaluate the effect that changing land use and groundwater pumping has on their occurrence and transport. Results of the trend analysis indicate water levels are declining ( p <0.05) in 17 of 26 selected monitoring wells (65 percent). Areas with the largest change in water levels, with more than 20 feet of declines, were within the community areas of Johnson Lane, Ruhenstroth, South Agricultural, East Valley, and Fish Springs. Variations in water levels measured in wells from the Central Agricultural, Minden, Foothill, Alpine County (one well), and Gardnerville Ranchos areas show periods of increase and decrease over time, but they also maintain long-term declining trends ( p <0.05). Increases in nitrate concentrations in groundwater samples collected from 9 out of 14 selected monitoring wells (64 percent) are statistically significant ( p <0.05) within the Ruhenstroth, Gardnerville Ranchos, East Valley, Genoa, and Johnson Lane community areas. Samples collected from a well in Indian Hills/Jacks Valley indicated a decreasing trend in nitrate concentration over time. Nitrate concentrations in samples collected from wells in East Valley, Genoa, Johnson Lane, and Indian Hills/Jack Valley were consistently low (less than 3 milligrams per liter [mg/L]) and stable. Nitrate concentrations from selected wells in Johnson Lane and Garnerville Ranchos exceeded the U.S. Environmental Protection Agency (EPA) maximum contaminant level (MCL) of 10 mg/L (as nitrogen) and have trends that are increasing over time. In 2022, a sample collected from Johnson Lane had a concentration (7.3 mg/L) below the MCL with an increasing trend over time. Temporal trend analyses for groundwater arsenic concentrations in Carson Valley could not be done because of a lack of temporal data. However, using available historical data, arsenic concentrations seem to be greater in groundwater from wells located on the eastern and northern areas of the valley than in wells located on the western or southern areas. Groundwater arsenic concentrations exceed 5 micrograms per liter (μg/L) in most samples collected from wells in Johnson Lane, Airport, Central Agricultural, and East Valley areas and in many cases exceed the U.S. Environmental Protection Agency (EPA) MCL of 10 μg/L. Data indicate that groundwater from domestic wells screened at deeper intervals are likely more vulnerable to elevated arsenic concentrations than shallower wells. A groundwater network evaluation for Carson Valley identified potential modifications in the sampling locations and frequency to better understand the effect of groundwater pumping in communities where municipal and domestic demand are greatest, potentially enhancing understanding of contaminant transport in these areas. Potential modifications to the active well network include reducing the frequency of sample collection from existing network wells (6 out of 11) that have consistently shown low and stable nitrate concentrations, adding wells in areas where data are sparse, and increasing the number of wells in areas with elevated groundwater nitrate concentrations. Including the analysis of arsenic in samples from the active groundwater monitoring well network will provide more detail on the temporal and spatial variability of arsenic concentrations. A visualization tool for the Carson River Basin was developed to provide access to discrete and near real-time hydrologic and water-quality data. The Carson River Basin Hydro Mapper (CBH; U.S. Geological Survey, 2023b) shows active and historical discrete water levels measured by the USGS and the State of Nevada Division of Water Resources, discrete groundwater nitrate and arsenic concentration data collected by the USGS, near real-time streamflow, and surface water levels for select waterbodies. The hydrologic data in the CBH provides resource managers, the public, and the scientific community with an easily accessible tool to present and communicate the most up-to-date information available about local and basin-wide water resources.

Nevada↗

Water-quality conditions, and constituent loads and yields in the Cambridge drinking-water source area, Massachusetts, water years 2005–07

The source water area for the drinking-water supply of the city of Cambridge, Massachusetts, encompasses major transportation corridors, as well as large areas of light industrial, commercial, and residential land use. Because of ongoing development in the drinking-water source area, the Cambridge water supply has the potential to be affected by a wide variety of contaminants. The U.S. Geological Survey (USGS) has monitored surface-water quality in the Hobbs Brook and Stony Brook Basins, which compose the drinking-water source area, since 1997 (water year 1997) through continuous monitoring and discrete sample collection and, since 2004, through systematic collection of streamwater samples during base-flow and stormflow conditions at five primary sampling stations in the drinking-water source area. Four primary sampling stations are on small tributaries in the Hobbs Brook and Stony Brook Basins; the fifth primary sampling station is on the main stem of Stony Brook and drains about 93 percent of the Cambridge drinking-water source area. Water samples also were collected at six secondary sampling stations, including Fresh Pond Reservoir, the final storage reservoir for the raw water supply. Storm runoff and base-flow concentrations of calcium (Ca), chloride (Cl), sodium (Na), and sulfate (SO 4 ) were estimated from continuous records of streamflow and specific conductance for six monitoring stations, which include the five primary sampling stations. These data were used to characterize current water-quality conditions, estimate loads and yields, and describe trends in Cl and Na in the tributaries and main-stem streams in the Hobbs Brook and Stony Brook Basins. These data also were used to describe how streamwater quality is affected by various watershed characteristics and provide information to guide future watershed management. Water samples were analyzed for physical properties and concentrations of Ca, Cl, Na, and SO 4 , total nitrogen (TN), total phosphorus (TP), caffeine, and a suite of 59 polar pesticides. Values of physical properties and constituent concentrations varied widely, particularly in samples from tributaries. Median concentrations of Ca, Cl, Na, and SO4 in samples collected in the Hobbs Brook Basin (39.8, 392, 207, and 21.7 milligrams per liter (mg/L), respectively) were higher than those for the Stony Brook Basin (17.8, 87.7, 49.7, and 14.7 mg/L, respectively). These differences in major ion concentrations are likely related to the low percentages of developed land and impervious area in the Stony Brook Basin. Concentrations of dissolved Cl and Na in samples, and those estimated from continuous records of specific conductance (particularly during base flow), often were greater than the U.S. Environmental Protection Agency (USEPA) secondary drinking-water guideline for Cl (250 mg/L), the chronic aquatic-life guideline for Cl (230 mg/L), and the Commonwealth of Massachusetts, Executive Office of Energy and Environmental Affairs drinking-water guideline for Na (20 mg/L). Mean annual flow-weighted concentrations of Ca, Cl, and Na were generally positively correlated with the area of roadway land use in the subbasins. Correlations between mean annual concentrations of Ca and SO 4 in base flow and total roadway, total impervious, and commercial-industrial land uses were statistically significant. Concentrations of TN (range of 0.42 to 5.13 mg/L in all subbasins) and TP (range of 0.006 to 0.80 mg/L in all subbasins) in tributary samples did not differ substantially between the Hobbs Brook and Stony Brook Basins. Concentrations of TN and TP in samples collected during water years 2004–07 exceeded proposed reference concentrations of 0.57 and 0.024 mg/L, in 94 and 56 percent of the samples, respectively. Correlations between annual flow-weighted concentrations of TN and percentages of recreational land use and water-body area were statistically significant; however, no significant relation was found between TP and available land-use information. The volume of streamflow affected water-quality conditions at the primary sampling stations. Turbidity and concentrations of TP were positively correlated with streamflow. In contrast, concentrations of major ions were negatively correlated with streamflow, indicating that these constituents were diluted during stormflows. Concentrations of TN were not correlated with streamflow. Twenty-five pesticides and caffeine were detected in water samples collected in the drinking-water source area and in raw water collected from the Cambridge water-treatment facility intake at the Fresh Pond Reservoir. Imidacloprid, norflurazon, and siduron were the most frequently detected pesticides with the frequency of detections ranging from about 24 to 41 percent. Caffeine was detected in about 37 percent of water samples at concentrations ranging from 0.003 to 1.82 micrograms per liter (μg/L). Although some of the detected pesticides degrade rapidly, norflurazon and siduron are relatively stable and are able to immigrate though the serial reservoir system. Concentrations of 2,4-D, carbaryl, imazaquin, MCPA (2-methyl-4-chlorophenoxyacetic acid), metsulfuron-methyl, norflurazon, siduron, and caffeine were detected more frequently in stormflow samples than in base-flow samples. Concentrations of pesticides did not exceed USEPA drinking-water guidelines or other health standards and were several orders of magnitude less than the lethal exposure level established for several fish species common to the drinking-water source area. Imidacloprid, an insecticide, was the only pesticide with a concentration exceeding available long-term aquatic-life guidelines. Several pesticides correlated significantly with the amount of recreational, residential, and commercial area in the tributary subbasins. Mean annual base-flow concentrations of caffeine correlated significantly with parking-lot land use. For most tributaries, about 70 percent of the annual loads of Ca, Cl, Na, and SO 4 were associated with base flow. Upward temporal trends in annual loads of Cl and Na were identified on the basis of data for water years 1998 to 2008 for the outlet of the Cambridge Reservoir in the Hobbs Brook Basin; however, similar trends were not identified for the main stem of Stony Brook downstream from the reservoir. The proportions of the TN load attributed to base flow and stormflow were similar in each tributary. In contrast, more than 83 percent of the TP loads in the tributaries and about 73 percent of the TP load in main stem of Stony Brook were associated with stormflow. Mean annual yields of Ca, Cl, Na, and SO 4 in the Stony Brook Reservoir watershed, which represents most of the drinking-water source area, were 14, 85, 46, and 9 metric tons per square kilometer, respectively. Mean annual yields among the individual tributary subbasins varied extensively. Mean annual yields for the respective constituents increased with an increase in roadway and parking-lot area in the tributary subbasins. Mean annual yields of TN in the tributary subbasins ranged from about 740 to more than 1,200 kilograms per square kilometer and exceeded the yield for the main stem of Stony Brook at USGS station 01104460 upstream from the Stony Brook Reservoir. Mean annual yields estimated for the herbicides 2,4-D and imidacloprid ranged from 34 to 310 grams per square kilometer (g/km 2 ) and 3 to 170 g/km 2 , respectively. Annual loads for 2,4-D were entirely associated with stormflow. The largest annual load for imidacloprid was estimated for the main stem of Stony Brook; however, the highest annual yield for this pesticide, as well as for benomyl, carbaryl, metalaxyl, and propiconazole, was estimated for a tributary to the Stony Brook Reservoir that drains largely residential and recreational areas. Mean annual yields for the herbicide siduron ranged from 6.9 to 35 g/km 2 with most of the loads associated with stormflow. Mean annual yields for the insecticide diuron ranged from 2.1 to 4.4 g/km 2 . Annual yields of caffeine ranged from 11 to 410 g/km 2 .

Massachusetts↗

Characterization and occurrence of confined and unconfined aquifers in Quaternary sediments in the glaciated conterminous United States

The glacial aquifer system, which is a collection of aquifers within Quaternary sediments in the glaciated conterminous United States, is a principal aquifer that supplies groundwater that serves about 42 million people and accounts for about 5 percent of the Nation’s drinking water. This aquifer system (the area of maximum glacial advance) underlies parts of 25 States and covers 1.87×10 6 square kilometers. A hydrogeologic framework is presented that divides the glaciated United States into 17 distinct hydrogeologic terranes using a geologic approach based on previous mapping. Each hydrogeologic terrane contains Quaternary sediment that is derived from a common depositional history and can be characterized by similar texture and thickness. Characteristics of Quaternary sediments are described using attributes computed from a lithologic database of well logs compiled from 24 States (excluding Kentucky). The hydrogeologic framework presents a nationwide picture of the glacial aquifer system and provides generalizations concerning the nature of aquifers within it (for example, whether the aquifers are shallow or deep, and unconfined or confined). In this way insights can be gained from understanding the similarities and differences in distinct parts of the glacial aquifer system and how they relate to water use and quality and to aquifer vulnerability. Delineation of hydrogeologic terranes was based on an interpretation of existing geologic mapping of Quaternary sediments and the thickness of unconsolidated material. Overall thickness of Quaternary sediment was used to qualitatively rank the generalized complexity of the hydrogeologic framework in each terrane: “lower” complexity (assigned a terrane code 1), “moderate” complexity (terrane code 2), and “higher” complexity (terrane code 3). Letter designations appended to the terrane codes (for example, 1A, 1B, or 1C) differentiate terranes of similar complexity. Two unique areas, where thick, stratified, coarse-grained sediment dominates, were assigned terrane code 4. Elements of this hydrogeologic framework include a glacial environments and surficial sediments geodatabase, which includes lithologic, geomorphic, and stratigraphic characterization of Quaternary sediments based on previous mapping; a gridded database of sediment and aquifer characteristics computed from lithologic logs obtained from water-well driller records; a water-use database with information on public-water supply systems and sources of groundwater; and estimated recharge computed from a geologically based soil-water balance model. A generalized map of the bedrock geology based on previous State-level mapping is included as well. Quaternary sediment in the glaciated United States includes glacial, postglacial (Holocene) and nonglacial sediments. At land surface, 60 percent of the glacial sediment is till. Large areas of outwash and ice contact sediments are extensive throughout the Midwest but generally are confined to valleys in the Northeast and the Northwest. Lacustrine sediments were deposited in proglacial lakes adjacent to the present Great Lakes and in glacial Lake Agassiz in the eastern Dakotas and northwestern Minnesota. The median thickness of Quaternary sediment ranges from 6 to 45 meters across the 17 hydrogeologic terranes, but the maximum thickness is more than 500 meters in some areas. Quaternary sediments generally contain less than 10 percent coarse material; the median range is near zero percent under till to about 50 percent under ice contact and outwash sediments. About 80 percent of the coarse material lies within 25 to 40 meters of land surface. In most of the glaciated United States, there is a small likelihood of penetrating an aquifer-material interval containing coarse material at least 3 meters thick. A single aquifer-material interval was recorded in about 44 percent of lithologic logs, whereas about 11 percent of the logs penetrated multiple intervals. About 44 percent of water wells in the lithologic database are completed in Quaternary sediment, and many of these Quaternary water wells (42 percent) are confined by at least 7.5 meters of fine materials. About 33 percent of these Quaternary water wells are unconfined—the remainder are where only thin layers (less than 3 meters) of coarse material are present. The median depths of Quaternary water wells range from 13 to 40 meters among the 17 hydrogeologic terranes. Recharge ranges from more than 400 millimeters per year in the Northeast to 11 millimeters or less per year in the Dakotas and Montana (median value of 136 millimeters per year). Annual groundwater withdrawals compiled by county range on an areal basis from less than 1 to 370 millimeters per year, and the mean is 7.4 millimeters per year. About 36 percent of the withdrawals are for public-water supply, of which 70 percent are derived from Quaternary sediments. Groundwater withdrawals are less than 10 percent of recharge throughout most of the glaciated conterminous United States but are a larger proportion of recharge near urban areas in the Northeast and the Midwest, and in counties throughout drier parts of the Midwest. The salient characteristics of the 17 hydrogeologic terranes are presented through maps and a set of descriptive plots to facilitate visual comparisons between selected sediment and aquifer characteristics. The thickness of Quaternary sediment generally increases from the lower complexity terranes through the higher complexity terranes, consistent with their delineation. Median proportions of coarse material in Quaternary sediment and depths to aquifer-material intervals are highly variable (less than 10 to 50 percent, and 0 to 30 meters, respectively). Median thicknesses of aquifer-material intervals generally fall within a narrow range (10 to 20 meters), except in two terranes that contain thick coarse-grained sediment (30 to 35 meters). The source of water in wells varies from mostly bedrock wells in the lower complexity terranes to mostly Quaternary wells in the higher complexity terranes where the sediment is thickest. A tree diagram compiled from a hierarchical cluster analysis of a matrix composed of metrics based on sediment and aquifer characteristics, and the distribution of water wells in each terrane, indicates some groups of terranes that can be treated as comparable when analyzing groundwater flow and quality. Aquifer-material intervals indicated on maps prepared from the lithologic logs, including unconfined and confined conditions, correlate well with aquifer systems delineated on state maps for Illinois, Indiana, and North Dakota. The large scale of the study limits the resolution at which the maps can be interpreted, however, and alluvial units are not mapped correctly for some valleys in the Northeast and the Northwest. Lithologic logs used in the study are biased toward shallow depths because not all logs penetrate the entire thickness of Quaternary sediment, but this bias should not limit the utility of the sediment and aquifer descriptions because shallow depths are commonly exploited for water supply. The hydrogeologic framework will support ongoing studies of groundwater flow and quality in the U.S. Geological Survey National Water Quality Assessment program for the glaciated United States.

Scientific Investigations Report↗

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Katawas mineral district in Afghanistan: Chapter N in Local-area-enhanced, high-resolution natural-color and color-infrared satellite-image mosaics of mineral districts in Afghanistan

The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Defense Task Force for Business and Stability Operations, prepared databases for mineral-resource target areas in Afghanistan. The purpose of the databases is to (1) provide useful data to ground-survey crews for use in performing detailed assessments of the areas and (2) provide useful information to private investors who are considering investment in a particular area for development of its natural resources. The set of satellite-image mosaics provided in this Data Series (DS) is one such database. Although airborne digital color-infrared imagery was acquired for parts of Afghanistan in 2006, the image data have radiometric variations that preclude their use in creating a consistent image mosaic for geologic analysis. Consequently, image mosaics were created using ALOS (Advanced Land Observation Satellite; renamed Daichi) satellite images, whose radiometry has been well determined (Saunier, 2007a,b). This part of the DS consists of the locally enhanced ALOS image mosaics for the Katawas mineral district, which has gold deposits. ALOS was launched on January 24, 2006, and provides multispectral images from the AVNIR (Advanced Visible and Near-Infrared Radiometer) sensor in blue (420-500 nanometer, nm), green (520-600 nm), red (610-690 nm), and near-infrared (760-890 nm) wavelength bands with an 8-bit dynamic range and a 10-meter (m) ground resolution. The satellite also provides a panchromatic band image from the PRISM (Panchromatic Remote-sensing Instrument for Stereo Mapping) sensor (520-770 nm) with the same dynamic range but a 2.5-m ground resolution. The image products in this DS incorporate copyrighted data provided by the Japan Aerospace Exploration Agency (©AXA, 2008), but the image processing has altered the original pixel structure and all image values of the JAXA ALOS data, such that original image values cannot be recreated from this DS. As such, the DS products match JAXA criteria for value added products, which are not copyrighted, according to the ALOS end-user license agreement. The selection criteria for the satellite imagery used in our mosaics were images having (1) the highest solar-elevation angles (near summer solstice) and (2) the least cloud, cloud-shadow, and snow cover. The multispectral and panchromatic data were orthorectified with ALOS satellite ephemeris data, a process which is not as accurate as orthorectification using digital elevation models (DEMs); however, the ALOS processing center did not have a precise DEM. As a result, the multispectral and panchromatic image pairs were generally not well registered to the surface and not coregistered well enough to perform resolution enhancement on the multispectral data. Therefore, it was necessary to (1) register the 10-m AVNIR multispectral imagery to a well-controlled Landsat image base, (2) mosaic the individual multispectral images into a single image of the entire area of interest, (3) register each panchromatic image to the registered multispectral image base, and (4) mosaic the individual panchromatic images into a single image of the entire area of interest. The two image-registration steps were facilitated using an automated control-point algorithm developed by the USGS that allows image coregistration to within one picture element. Before rectification, the multispectral and panchromatic images were converted to radiance values and then to relative-reflectance values using the methods described in Davis (2006). Mosaicking the multispectral or panchromatic images started with the image with the highest sun-elevation angle and the least atmospheric scattering, which was treated as the standard image. The band-reflectance values of all other multispectral or panchromatic images within the area were sequentially adjusted to that of the standard image by determining band-reflectance correspondence between overlapping images using linear least-squares analysis. The resolution of the multispectral image mosaic was then increased to that of the panchromatic image mosaic using the SPARKLE logic, which is described in Davis (2006). Each of the four-band images within the resolution-enhanced image mosaic was individually subjected to a local-area histogram stretch algorithm (described in Davis, 2007), which stretches each band's picture element based on the digital values of all picture elements within a 315-m radius. The final databases, which are provided in this DS, are three-band, color-composite images of the local-area-enhanced, natural-color data (the blue, green, and red wavelength bands) and color-infrared data (the green, red, and near-infrared wavelength bands). All image data were initially projected and maintained in Universal Transverse Mercator (UTM) map projection using the target area's local zone (42 for Katawas) and the WGS84 datum. The final image mosaics are provided as embedded geotiff images, which can be read and used by most geographic information system (GIS) and image-processing software. The tiff world files (tfw) are provided, even though they are generally not needed for most software to read an embedded geotiff image. Within the Katawas study area, one subarea was designated for detailed field investigation (that is, the Gold subarea); this subarea was extracted from the area's image mosaic and is provided as a separate embedded geotiff image.

Katawas Mineral District↗

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Takhar mineral district in Afghanistan: Chapter Q in Local-area-enhanced, high-resolution natural-color and color-infrared satellite-image mosaics of mineral districts in Afghanistan

The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Defense Task Force for Business and Stability Operations, prepared databases for mineral-resource target areas in Afghanistan. The purpose of the databases is to (1) provide useful data to ground-survey crews for use in performing detailed assessments of the areas and (2) provide useful information to private investors who are considering investment in a particular area for development of its natural resources. The set of satellite-image mosaics provided in this Data Series (DS) is one such database. Although airborne digital color-infrared imagery was acquired for parts of Afghanistan in 2006, the image data have radiometric variations that preclude their use in creating a consistent image mosaic for geologic analysis. Consequently, image mosaics were created using ALOS (Advanced Land Observation Satellite; renamed Daichi) satellite images, whose radiometry has been well determined (Saunier, 2007a,b). This part of the DS consists of the locally enhanced ALOS image mosaics for the Takhar mineral district, which has industrial evaporite deposits. ALOS was launched on January 24, 2006, and provides multispectral images from the AVNIR (Advanced Visible and Near-Infrared Radiometer) sensor in blue (420–500 nanometer, nm), green (520–600 nm), red (610–690 nm), and near-infrared (760–890 nm) wavelength bands with an 8-bit dynamic range and a 10-meter (m) ground resolution. The satellite also provides a panchromatic band image from the PRISM (Panchromatic Remote-sensing Instrument for Stereo Mapping) sensor (520–770 nm) with the same dynamic range but a 2.5-m ground resolution. The image products in this DS incorporate copyrighted data provided by the Japan Aerospace Exploration Agency (©JAXA, 2008), but the image processing has altered the original pixel structure and all image values of the JAXA ALOS data, such that original image values cannot be recreated from this DS. As such, the DS products match JAXA criteria for value added products, which are not copyrighted, according to the ALOS end-user license agreement. The selection criteria for the satellite imagery used in our mosaics were images having (1) the highest solar-elevation angles (near summer solstice) and (2) the least cloud, cloud-shadow, and snow cover. The multispectral and panchromatic data were orthorectified with ALOS satellite ephemeris data, a process which is not as accurate as orthorectification using digital elevation models (DEMs); however, the ALOS processing center did not have a precise DEM. As a result, the multispectral and panchromatic image pairs were generally not well registered to the surface and not coregistered well enough to perform resolution enhancement on the multispectral data. Therefore, it was necessary to (1) register the 10-m AVNIR multispectral imagery to a well-controlled Landsat image base, (2) mosaic the individual multispectral images into a single image of the entire area of interest, (3) register each panchromatic image to the registered multispectral image base, and (4) mosaic the individual panchromatic images into a single image of the entire area of interest. The two image-registration steps were facilitated using an automated control-point algorithm developed by the USGS that allows image coregistration to within one picture element. Before rectification, the multispectral and panchromatic images were converted to radiance values and then to relative-reflectance values using the methods described in Davis (2006). Mosaicking the multispectral or panchromatic images started with the image with the highest sun-elevation angle and the least atmospheric scattering, which was treated as the standard image. The band-reflectance values of all other multispectral or panchromatic images within the area were sequentially adjusted to that of the standard image by determining band-reflectance correspondence between overlapping images using linear least-squares analysis. The resolution of the multispectral image mosaic was then increased to that of the panchromatic image mosaic using the SPARKLE logic, which is described in Davis (2006). Each of the four-band images within the resolution-enhanced image mosaic was individually subjected to a local-area histogram stretch algorithm (described in Davis, 2007), which stretches each band’s picture element based on the digital values of all picture elements within a 315-m radius. The final databases, which are provided in this DS, are three-band, color-composite images of the local-area-enhanced, natural-color data (the blue, green, and red wavelength bands) and color-infrared data (the green, red, and near-infrared wavelength bands). All image data were initially projected and maintained in Universal Transverse Mercator (UTM) map projection using the target area’s local zone (42 for Takhar) and the WGS84 datum. The final image mosaics for the Takhar area are provided as embedded geotiff images, which can be read and used by most geographic information system (GIS) and image-processing software. The tiff world files (tfw) are provided, even though they are generally not needed for most software to read an embedded geotiff image.

Takhar Mineral District↗

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Ghazni2 mineral district in Afghanistan

The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Defense Task Force for Business and Stability Operations, prepared databases for mineral-resource target areas in Afghanistan. The purpose of the databases is to (1) provide useful data to ground-survey crews for use in performing detailed assessments of the areas and (2) provide useful information to private investors who are considering investment in a particular area for development of its natural resources. The set of satellite-image mosaics provided in this Data Series (DS) is one such database. Although airborne digital color-infrared imagery was acquired for parts of Afghanistan in 2006, the image data have radiometric variations that preclude their use in creating a consistent image mosaic for geologic analysis. Consequently, image mosaics were created using ALOS (Advanced Land Observation Satellite; renamed Daichi) satellite images, whose radiometry has been well determined (Saunier, 2007a,b). This part of the DS consists of the locally enhanced ALOS image mosaics for the Ghazni2 mineral district, which has spectral reflectance anomalies indicative of gold, mercury, and sulfur deposits. ALOS was launched on January 24, 2006, and provides multispectral images from the AVNIR (Advanced Visible and Near-Infrared Radiometer) sensor in blue (420-500 nanometer, nm), green (520-600 nm), red (610-690 nm), and near-infrared (760-890 nm) wavelength bands with an 8-bit dynamic range and a 10-meter (m) ground resolution. The satellite also provides a panchromatic band image from the PRISM (Panchromatic Remote-sensing Instrument for Stereo Mapping) sensor (520-770 nm) with the same dynamic range but a 2.5-m ground resolution. The image products in this DS incorporate copyrighted data provided by the Japan Aerospace Exploration Agency ((c)JAXA, 2008, 2009), but the image processing has altered the original pixel structure and all image values of the JAXA ALOS data, such that original image values cannot be recreated from this DS. As such, the DS products match JAXA criteria for value added products, which are not copyrighted, according to the ALOS end-user license agreement. The selection criteria for the satellite imagery used in our mosaics were images having (1) the highest solar-elevation angles (near summer solstice) and (2) the least cloud, cloud-shadow, and snow cover. The multispectral and panchromatic data were orthorectified with ALOS satellite ephemeris data, a process which is not as accurate as orthorectification using digital elevation models (DEMs); however, the ALOS processing center did not have a precise DEM. As a result, the multispectral and panchromatic image pairs were generally not well registered to the surface and not coregistered well enough to perform resolution enhancement on the multispectral data. Therefore, it was necessary to (1) register the 10-m AVNIR multispectral imagery to a well-controlled Landsat image base, (2) mosaic the individual multispectral images into a single image of the entire area of interest, (3) register each panchromatic image to the registered multispectral image base, and (4) mosaic the individual panchromatic images into a single image of the entire area of interest. The two image-registration steps were facilitated using an automated control-point algorithm developed by the USGS that allows image coregistration to within one picture element. Before rectification, the multispectral and panchromatic images were converted to radiance values and then to relative-reflectance values using the methods described in Davis (2006). Mosaicking the multispectral or panchromatic images started with the image with the highest sun-elevation angle and the least atmospheric scattering, which was treated as the standard image. The band-reflectance values of all other multispectral or panchromatic images within the area were sequentially adjusted to that of the standard image by determining band-reflectance correspondence between overlapping images using linear least-squares analysis. The resolution of the multispectral image mosaic was then increased to that of the panchromatic image mosaic using the SPARKLE logic, which is described in Davis (2006). Each of the four-band images within the resolution-enhanced image mosaic was individually subjected to a local-area histogram stretch algorithm (described in Davis, 2007), which stretches each band's picture element based on the digital values of all picture elements within a 500-m radius. The final databases, which are provided in this DS, are three-band, color-composite images of the local-area-enhanced, natural-color data (the blue, green, and red wavelength bands) and color-infrared data (the green, red, and near-infrared wavelength bands). All image data were initially projected and maintained in Universal Transverse Mercator (UTM) map projection using the target area's local zone (42 for Ghazni2) and the WGS84 datum. The images for the Ghazni2 area are provided as embedded geotiff images, which can be read and used by most geographic information system (GIS) and image-processing software. The tiff world files (tfw) are provided, even though they are generally not needed for most software to read an embedded geotiff image.

Ghazni2 Mineral District↗

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Parwan mineral district in Afghanistan: Chapter CC in Local-area-enhanced, high-resolution natural-color and color-infrared satellite-image mosaics of mineral districts in Afghanistan

The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Defense Task Force for Business and Stability Operations, prepared databases for mineral-resource target areas in Afghanistan. The purpose of the databases is to (1) provide useful data to ground-survey crews for use in performing detailed assessments of the areas and (2) provide useful information to private investors who are considering investment in a particular area for development of its natural resources. The set of satellite-image mosaics provided in this Data Series (DS) is one such database. Although airborne digital color-infrared imagery was acquired for parts of Afghanistan in 2006, the image data have radiometric variations that preclude their use in creating a consistent image mosaic for geologic analysis. Consequently, image mosaics were created using ALOS (Advanced Land Observation Satellite; renamed Daichi) satellite images, whose radiometry has been well determined (Saunier, 2007a,b). This part of the DS consists of the locally enhanced ALOS image mosaics for the Parwan mineral district, which has gold and copper deposits. ALOS was launched on January 24, 2006, and provides multispectral images from the AVNIR (Advanced Visible and Near-Infrared Radiometer) sensor in blue (420–500 nanometer, nm), green (520–600 nm), red (610–690 nm), and near-infrared (760–890 nm) wavelength bands with an 8-bit dynamic range and a 10-meter (m) ground resolution. The satellite also provides a panchromatic band image from the PRISM (Panchromatic Remote-sensing Instrument for Stereo Mapping) sensor (520–770 nm) with the same dynamic range but a 2.5-m ground resolution. The image products in this DS incorporate copyrighted data provided by the Japan Aerospace Exploration Agency (©JAXA,2006, 2007), but the image processing has altered the original pixel structure and all image values of the JAXA ALOS data, such that original image values cannot be recreated from this DS. As such, the DS products match JAXA criteria for value added products, which are not copyrighted, according to the ALOS end-user license agreement. elevation angles (near summer solstice) and (2) the least cloud, cloud-shadow, and snow cover. The multispectral and panchromatic data were orthorectified with ALOS satellite ephemeris data, a process which is not as accurate as orthorectification using digital elevation models (DEMs); however, the ALOS processing center did not have a precise DEM. As a result, the multispectral and panchromatic image pairs were generally not well registered to the surface and not coregistered well enough to perform resolution enhancement on the multispectral data. Therefore, it was necessary to (1) register the 10-m AVNIR multispectral imagery to a well-controlled Landsat image base, (2) mosaic the individual multispectral images into a single image of the entire area of interest, (3) register each panchromatic image to the registered multispectral image base, and (4) mosaic the individual panchromatic images into a single image of the entire area of interest. The two image-registration steps were facilitated using an automated control-point algorithm developed by the USGS that allows image coregistration to within one picture element. Before rectification, the multispectral and panchromatic images were converted to radiance values and then to relative-reflectance values using the methods described in Davis (2006). Mosaicking the multispectral or panchromatic images started with the image with the highest sun-elevation angle and the least atmospheric scattering, which was treated as the standard image. The band-reflectance values of all other multispectral or panchromatic images within the area were sequentially adjusted to that of the standard image by determining band-reflectance correspondence between overlapping images using linear least-squares analysis. The resolution of the multispectral image mosaic was then increased to that of the panchromatic image mosaic using the SPARKLE logic, which is described in Davis (2006). Each of the four-band images within the resolution-enhanced image mosaic was individually subjected to a local-area histogram stretch algorithm (described in Davis, 2007), which stretches each band’s picture element based on the digital values of all picture elements within a 500-m radius. The final databases, which are provided in this DS, are three-band, color-composite images of the local-area-enhanced, natural-color data (the blue, green, and red wavelength bands) and color-infrared data (the green, red, and near-infrared wavelength bands). All image data were initially projected and maintained in Universal Transverse Mercator (UTM) map projection using the target area’s local zone (42 for Parwan) and the WGS84 datum. The final image mosaics were subdivided into two overlapping tiles or quadrants because of the large size of the target area. The two image tiles (or quadrants) for the North Bamyan area are provided as embedded geotiff images, which can be read and used by most geographic information system (GIS) and image-processing software. The tiff world files (tfw) are provided, even though they are generally not needed for most software to read an embedded geotiff image.

Parwan Mineral District↗

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Ghazni1 mineral district in Afghanistan

The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Defense Task Force for Business and Stability Operations, prepared databases for mineral-resource target areas in Afghanistan. The purpose of the databases is to (1) provide useful data to ground-survey crews for use in performing detailed assessments of the areas and (2) provide useful information to private investors who are considering investment in a particular area for development of its natural resources. The set of satellite-image mosaics provided in this Data Series (DS) is one such database. Although airborne digital color-infrared imagery was acquired for parts of Afghanistan in 2006, the image data have radiometric variations that preclude their use in creating a consistent image mosaic for geologic analysis. Consequently, image mosaics were created using ALOS (Advanced Land Observation Satellite; renamed Daichi) satellite images, whose radiometry has been well determined (Saunier, 2007a,b). This part of the DS consists of the locally enhanced ALOS image mosaics for the Ghazni1 mineral district, which has spectral reflectance anomalies indicative of clay, aluminum, gold, silver, mercury, and sulfur deposits. ALOS was launched on January 24, 2006, and provides multispectral images from the AVNIR (Advanced Visible and Near-Infrared Radiometer) sensor in blue (420-500 nanometer, nm), green (520-600 nm), red (610-690 nm), and near-infrared (760-890 nm) wavelength bands with an 8-bit dynamic range and a 10-meter (m) ground resolution. The satellite also provides a panchromatic band image from the PRISM (Panchromatic Remote-sensing Instrument for Stereo Mapping) sensor (520-770 nm) with the same dynamic range but a 2.5-m ground resolution. The image products in this DS incorporate copyrighted data provided by the Japan Aerospace Exploration Agency ((c)JAXA, 2008, 2009), but the image processing has altered the original pixel structure and all image values of the JAXA ALOS data, such that original image values cannot be recreated from this DS. As such, the DS products match JAXA criteria for value added products, which are not copyrighted, according to the ALOS end-user license agreement. The selection criteria for the satellite imagery used in our mosaics were images having (1) the highest solar-elevation angles (near summer solstice) and (2) the least cloud, cloud-shadow, and snow cover. The multispectral and panchromatic data were orthorectified with ALOS satellite ephemeris data, a process which is not as accurate as orthorectification using digital elevation models (DEMs); however, the ALOS processing center did not have a precise DEM. As a result, the multispectral and panchromatic image pairs were generally not well registered to the surface and not coregistered well enough to perform resolution enhancement on the multispectral data. Therefore, it was necessary to (1) register the 10-m AVNIR multispectral imagery to a well-controlled Landsat image base, (2) mosaic the individual multispectral images into a single image of the entire area of interest, (3) register each panchromatic image to the registered multispectral image base, and (4) mosaic the individual panchromatic images into a single image of the entire area of interest. The two image-registration steps were facilitated using an automated control-point algorithm developed by the USGS that allows image coregistration to within one picture element. Before rectification, the multispectral and panchromatic images were converted to radiance values and then to relative-reflectance values using the methods described in Davis (2006). Mosaicking the multispectral or panchromatic images started with the image with the highest sun-elevation angle and the least atmospheric scattering, which was treated as the standard image. The band-reflectance values of all other multispectral or panchromatic images within the area were sequentially adjusted to that of the standard image by determining band-reflectance correspondence between overlapping images using linear least-squares analysis. The resolution of the multispectral image mosaic was then increased to that of the panchromatic image mosaic using the SPARKLE logic, which is described in Davis (2006). Each of the four-band images within the resolution-enhanced image mosaic was individually subjected to a local-area histogram stretch algorithm (described in Davis, 2007), which stretches each band's picture element based on the digital values of all picture elements within a 500-m radius. The final databases, which are provided in this DS, are three-band, color-composite images of the local-area-enhanced, natural-color data (the blue, green, and red wavelength bands) and color-infrared data (the green, red, and near-infrared wavelength bands). All image data were initially projected and maintained in Universal Transverse Mercator (UTM) map projection using the target area's local zone (42 for Ghazni1) and the WGS84 datum. The images for the Ghazni1 area are provided as embedded geotiff images, which can be read and used by most geographic information system (GIS) and image-processing software. The tiff world files (tfw) are provided, even though they are generally not needed for most software to read an embedded geotiff image.

Ghazni1 Mineral District↗