USGS Science⌕ Search

SEARCH · USGS Science

Results for “Internal Report”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

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

Detrital zircon U-Pb reconnaissance of the Franciscan subduction complex in northwestern California

In northwestern California, the Franciscan subduction complex has been subdivided into seven major tectonostratigraphic units. We report U-Pb ages of &asymp;2400 detrital zircon grains from 26 sandstone samples from 5 of these units. Here, we tabulate each unit's interpreted predominant sediment source areas and depositional age range, ordered from the oldest to the youngest unit. (1) Yolla Bolly terrane: nearby Sierra Nevada batholith (SNB); ca. 118 to 98 Ma. Rare fossils had indicated that this unit was mostly 151-137 Ma, but it is mostly much younger. (2) Central Belt: SND; ca. 103 too 53 Ma (but poorly constrained), again mostly younger than previously thought. (3) Yager terrane: distant Idaho batholith (IB); ca. 52 to 50 Ma. Much of the Yager's detritus was shed during major core complex extension and erosion in Idaho that started 53 Ma. An eocene Princeton River-Princeton submarine canyon system transported this detritus to the Great Valley forearc basin and thence to the Franciscan trench. (4) Coastal terrane: mostly IB, &plusmn;SNB, &plusmn;nearby Cascade arc, &plusmn;Nevada Cenozoic ignimbrite belt; 52 to <32 Ma. (5) King Range terrane: dominated by IB and SNB zircons; parts 16-14 Ma based on microfossils. Overall, some Franciscan units are younger than previously thought, making them more compatible with models for the growth of subduction complexes by positive accretion. From ca. 118 to 70 Ma, Franciscan sediments were sourced mainly from the nearby Sierra Nevada region and were isolated from southwestern US and Mexican sources. From 53 to 49 Ma, the Franciscan was sourced from both Idaho and the Sierra Nevada. By 37-32 Ma, input from Idaho had ceased. The influx from Idaho probably reflects major tectonism in Idaho, Oregon, and Washington, plus development of a through-going Princeton River to California, rather than radical changes in the subduction system at the Franciscan trench itself.

California↗

A total intensity magnetic anomaly map of the Red Sea and its interpretation

All available total intensity magnetic field data have been used to prepare a total intensity magnetic anomaly map of the Red Sea at a scale of 1:2,000,000. These include airborne surveys at various elevations and numerous sea surveys. The data have been digitised at approximately one kilometre intervals and the International Geomagnetic Reference Field (IGRF epoch 1965.0) removed. The resulting values have been upward continued to a common height of 1.83 km (6000 ft), transferred to a base map and contoured at 100 nT intervals. The map is published in colour by the United States Geological Survey - Saudi Arabian Project. The map reveals three distinct zones: (1) large amplitude (approximately 800 nT), short wavelength (15 km peak to trough) anomalies over the axial trough which form lineations roughly parallel to the coastlines, (2) smaller amplitude (approximately 300 nT), longer wavelength (30 km) anomalies over the main trough and shelves which also form lineations parallel to the coastlines, and (3) very variable amplitude (100 to 400 nT), short wavelength (5 to 30 km) anomalies over the shield areas on both sides of the Red Sea. The large, linear anomalies terminate in the north at about 23 °N and in the south at about 15 N. The anomalies over the main trough and shelves converge south of approximately 15.5 N. Several NE trending magnetic features offset the linear, axial anomalies. Such features are confined to zones (1) and (2) and are thought to indicate leaky transform faults. The interpretation of the magnetic data is presented in three parts: the northern Red Sea (north of 24°N), the central Red Sea (18°N to 24°N) and the southern Red S ea (south of 18 N). For each part, there is a discussion of the extent of oceanic crust, the seafloor spreading history, the possible location of transform faults and the relationship of the magnetic anomalies to the deep holes. In Chapter 7, the results are brought together and their significance for Red Sea floor spreading, transforms and plate geometry discussed. The problems of the brine deeps are briefly discussed in Chapter 8.

Open-File Report↗

Geologic field trip guide to Mount Mazama and Crater Lake Caldera, Oregon

Crater Lake partly fills one of the most spectacular calderas of the world—an 8 by 10 kilometer (km) basin more than 1 km deep formed by collapse of the Mount Mazama volcano during a rapid series of explosive eruptions ~7,700 years ago. Having a maximum depth of 594 meters (m), Crater Lake is the deepest lake in the United States. Crater Lake National Park, dedicated in 1902, encompasses 645 square kilometers (km 2 ) of pristine forested and alpine terrain, including the lake itself, and virtually all of Mount Mazama. The geology of the area was first described in detail by Diller and Patton (1902) and later by Williams (1942), whose vivid account led to international recognition of Crater Lake as the classic collapse caldera. Because of excellent preservation and access, Mount Mazama, Crater Lake caldera, and the deposits formed by the climactic eruption constitute a natural laboratory for study of volcanic and magmatic processes. For example, the climactic ejecta are renowned among volcanologists as evidence for systematic compositional zonation within a subterranean magma chamber. Mount Mazama’s climactic eruption also is important as the source of the widespread Mazama ash, a useful Holocene stratigraphic marker throughout the Pacific Northwest United States, adjacent Canada, and offshore. A detailed bathymetric survey of the floor of Crater Lake in 2000 (Bacon and others, 2002) provides a unique record of postcaldera eruptions, the interplay between volcanism and filling of the lake, and sediment transport within this closed basin. Knowledge of the geology and eruptive history of the Mount Mazama edifice, enhanced by the caldera wall exposures, gives exceptional insight into how large volcanoes of magmatic arcs grow and evolve. In addition, many smaller volcanoes of the High Cascades beyond the limits of Mount Mazama provide information on the flux of mantle-derived magma through the region. General principles of magmatic and eruptive processes revealed by geologic research at Crater Lake have been incorporated not only in scientific investigations elsewhere, but also in the practical evaluation of local hazards (Bacon and others, 1997b) and geothermal resources (Bacon and Nathenson, 1996). The 1:24,000-scale geologic map of Mount Mazama and Crater Lake caldera (Bacon, 2008) is unusual because it portrays bedrock (outcrop), surficial, and lake floor geology. Caldera wall geology is depicted in detail on the accompanying geologic panoramas, and bedrock geology is shown in a 1:50,000-scale geologic map. This field guide supersedes earlier geology guides of Crater Lake (Bacon, 1987, 1989).

Oregon↗

Sedimentary carbonate-hosted giant Bayan Obo REE-Fe-Nb ore deposit of Inner Mongolia, China; a cornerstone example for giant polymetallic ore deposits of hydrothermal origin

Detailed, integrative field and laboratory studies of the textures, structures, chemical characteristics, and isotopically determined ages and signatures of mineralization of the Bayan Obo deposit provided evidence for the origin and characteristics favorable for its formation and parameters necessary for defining giant polymetallic deposits of hydrothermal origin. Bayan Obo is an epigenetic, metasomatic, hydrothermal rare earth element (REE)-Fe-Nb ore deposit that is hosted in the metasedimentary H8 dolostone marble of the Middle Proterozoic Bayan Obo Group. The metasedimentary sequence was deposited on the northern continental slope of the North China craton. The mine area is about 100 km south of the suture marking Caledonian subduction of the Mongolian oceanic plate from the north beneath the North China craton. The mineralogy of the deposit is very complex, consisting of more than 120 different minerals, some of which are epigenetic minerals introduced by hydrothermal solutions, and some of which are primary and secondary metamorphic minerals. The major REE minerals are monazite and bastnaesite, whereas magnetite and hematite are the dominant Fe-ore minerals, and columbite is the most abundant Nb mineral. Dolomite, alkali amphibole, fluorite, barite, aegirine augite, apatite, phlogopite, albite, and microcline are the most widespread gangue minerals. Three general types of ores occur at Bayan Obo: disseminated, banded, and massive ores. Broad zoning of these ore types occurs in the Main and East Orebodies. Disseminated ores are in the outermost zone, banded ores are in the intermediate zone, and massive ores are in the cores of the orebodies. On the basis of field relations, host rocks, textures, structures, and mineral assemblages, many varieties of these three types of ores have been recognized and mapped. Isotopic dating of monazite, bastnaesite, aeschynite, and metamorphic and metasomatic alkali amphiboles associated with the deposit provides constraints on the ages of mineralization and the history of the deposit. Textural relations, differences in chemical composition, and 232Th/208Pb internal isochron ages of monazite and bastnaesite samples indicate that many episodes of REE mineralization occurred at Bayan Obo, ranging from about 555 Ma to about 398 Ma. Initial 208Pb/204Pb ratios suggest different sources of REE's for different generations of REE minerals. Relative ages of Fe mineralization were deduced from textural relationships of Fe minerals with other, dated mineral phases in the deposit. Most Nb mineralization was in the area of the West Orebodies and resulted in disseminated ore. Aeschynite, an early stage of Nb mineralization (438+-25.1 Ma), occurs with huanghoite and alkali amphiboles in veins. The 40Ar/39Ar ages of amphiboles, as well as petrographic textures, were used to distinguish three periods of regional metamorphism in the Bayan Obo mine area: (1) Late Proterozoic, about 890 Ma, which recrystallized H8 carbonate to marble and crystallized lineated alkali amphiboles along foliation planes in the marble; (2) Caledonian, about 425-395 Ma, which resulted in metamorphic and metasomatic-metamorphic alkali amphiboles; and (3) Hercynian, about 300 Ma, based on biotite 40Ar/39Ar analyses from biotite schist and folded banded ores. The 40Ar/39Ar ages of metasomatic alkali amphiboles also place time constraints on the hydrothermal history of the ore deposit. Metasomatic amphiboles represent periods of intense hydrothermal activity, which began as early as 1.26 Ga; that date is based on the age of amphibole from a vein that crosscuts the H6 quartzite that underlies the H8 dolostone marble. Although much of the metasomatic amphibole formed during periods that overlapped the peak period of REE mineralization of banded ores, REE and alkali amphibole phases generally occur in different mineral assemblages or are of very different ages in the same assemblage and, therefore, may have been derived from

Bulletin↗

A rapid diagnostic test and mobile "lab in a suitcase" platform for detecting Ceratocystis spp. responsible for Rapid ‘Ōhi‘a Death

We describe a field compatible molecular diagnostic test for two new species of Ceratocystis that infect `ōhi`a (Metrosideros polymorpha) and cause the disease commonly known as Rapid `Ōhi`a Death. The diagnostic is based on amplification of a DNA locus within the internal transcribed spacer region that separates fungal 5.8S ribosomal genes. The assay uses forward and reverse primers, recombinase polymerase, and a fluorescent probe that allows isothermal (40oC) amplification and simultaneous quantification of a 115 base pair product with a battery operated fluorometer. DNA extractions are field compatible and can be done by heating wood drill shavings to 100oC in Instagene® solution containing Chelex® resin to bind potential amplification inhibitors. The initial heat treatment is followed by a short bead beating step with steel ball bearings and zirconium beads to release DNA. DNA is subsequently purified with a magnetic bead based extraction method that does not require silica columns or centrifugation. The assay is designed around a portable “lab-in-a-suitcase” platform that includes a portable fluorometer, miniature centrifuge, and heat block that operate off either 120V AC power sources or a 12 volt battery with a portable inverter, a magnetic rack designed for 1.5 ml tubes and magnetic bead DNA purification, pipettes and consumable reagents and tubes. The entire assay from DNA extraction to results can be performed in less than 90 minutes on up to six independent samples plus a positive and negative control. Sensitivity based on suspensions of Ceratocystis endoconidia (spores) that were added to wood shavings and processed under field conditions by Instagene® magnetic bead DNA extraction was up to 163 spores/mg wood for Species A and 55 spores/mg wood for Species B in 95% of replicates as determined by probit analysis. Sensitivity increased 5–10 fold to 19 spores/mg wood for Species A and 9 spores/mg wood for Species B when extractions were performed with a commercial, silica column based DNA purification kit. The test did not cross react with other common fungi that have been isolated from `ōhi`a.

Technical Report↗

GIS Representation of Coal-Bearing Areas in North, Central, and South America

Worldwide coal consumption and international coal trade are projected to increase in the next several decades (Energy Information Administration, 2007). A search of existing literature indicates that in the Western Hemisphere, coal resources are known to occur in about 30 countries. The need exists to be able to depict these areas in a digital format for use in Geographic Information System (GIS) applications at small scales (large areas) and in visual presentations. Existing surficial geology GIS layers of the appropriate geologic age have been used as an approximation to depict the extent of coal-bearing areas in North, Central, and South America, as well as Greenland (fig. 1). Global surficial geology GIS data were created by the U.S. Geological Survey (USGS) for use in world petroleum assessments (Hearn and others, 2003). These USGS publications served as the major sources for the selection and creation of polygons to represent coal-bearing areas. Additional publications and maps by various countries and agencies were also used as sources of coal locations. GIS geologic polygons were truncated where literature or hardcopy maps did not indicate the presence of coal. The depicted areas are not adequate for use in coal resource calculations, as they were not adjusted for geologic structure and do not include coal at depth. Additionally, some coal areas in Central America could not be represented by the mapped surficial geology and are shown only as points based on descriptions or depictions from scientific publications or available maps. The provided GIS files are intended to serve as a backdrop for display of coal information. Three attributes of the coal that are represented by the polygons or points include geologic age (or range of ages), published rank (or range of ranks), and information source (published sources for age, rank, or physical location, or GIS geology base).

Open-File Report↗

User's Guide for the MapImage Reprojection Software Package, Version 1.01

Scientists routinely accomplish small-scale geospatial modeling in the raster domain, using high-resolution datasets (such as 30-m data) for large parts of continents and low-resolution to high-resolution datasets for the entire globe. Recently, Usery and others (2003a) expanded on the previously limited empirical work with real geographic data by compiling and tabulating the accuracy of categorical areas in projected raster datasets of global extent. Geographers and applications programmers at the U.S. Geological Survey's (USGS) Mid-Continent Mapping Center (MCMC) undertook an effort to expand and evolve an internal USGS software package, MapImage, or mapimg, for raster map projection transformation (Usery and others, 2003a). Daniel R. Steinwand of Science Applications International Corporation, Earth Resources Observation Systems Data Center in Sioux Falls, S. Dak., originally developed mapimg for the USGS, basing it on the USGS's General Cartographic Transformation Package (GCTP). It operated as a command line program on the Unix operating system. Through efforts at MCMC, and in coordination with Mr. Steinwand, this program has been transformed from an application based on a command line into a software package based on a graphic user interface for Windows, Linux, and Unix machines. Usery and others (2003b) pointed out that many commercial software packages do not use exact projection equations and that even when exact projection equations are used, the software often results in error and sometimes does not complete the transformation for specific projections, at specific resampling resolutions, and for specific singularities. Direct implementation of point-to-point transformation with appropriate functions yields the variety of projections available in these software packages, but implementation with data other than points requires specific adaptation of the equations or prior preparation of the data to allow the transformation to succeed. Additional constraints apply to global raster data. It appears that some packages use the USGS's GCTP or similar point transformations without adaptation to the specific characteristics of raster data (Usery and others, 2003b). It is most common for programs to compute transformations of raster data in an inverse fashion. Such mapping can result in an erroneous position and replicate data or create pixels not in the original space. As Usery and others (2003a) indicated, mapimg performs a corresponding forward transformation to ensure the same location results from both methods. The primary benefit of this function is to mask cells outside the domain. MapImage 1.01 is now on the Web. You can download the User's Guide, source, and binaries from the following site: http://mcmcweb.er.usgs.gov/carto_research/projection/acc_proj_data.html

Open-File Report↗

Seismological analyses of the 2010 March 11, Pichilemu, Chile Mw 7.0 and Mw 6.9 coastal intraplate earthquakes

On 2010 March 11, a sequence of large, shallow continental crust earthquakes shook central Chile. Two normal faulting events with magnitudes around M w 7.0 and M w 6.9 occurred just 15 min apart, located near the town of Pichilemu. These kinds of large intraplate, inland crustal earthquakes are rare above the Chilean subduction zone, and it is important to better understand their relationship with the 2010 February 27, M w 8.8, Maule earthquake, which ruptured the adjacent megathrust plate boundary. We present a broad seismological analysis of these earthquakes by using both teleseismic and regional data. We compute seismic moment tensors for both events via a W-phase inversion, and test sensitivities to various inversion parameters in order to assess the stability of the solutions. The first event, at 14 hr 39 min GMT, is well constrained, displaying a fault plane with strike of N145°E, and a preferred dip angle of 55°SW, consistent with the trend of aftershock locations and other published results. Teleseismic finite-fault inversions for this event show a large slip zone along the southern part of the fault, correlating well with the reported spatial density of aftershocks. The second earthquake (14 hr 55 min GMT) appears to have ruptured a fault branching southward from the previous ruptured fault, within the hanging wall of the first event. Modelling seismograms at regional to teleseismic distances (Δ > 10°) is quite challenging because the observed seismic wave fields of both events overlap, increasing apparent complexity for the second earthquake. We perform both point- and extended-source inversions at regional and teleseismic distances, assessing model sensitivities resulting from variations in fault orientation, dimension, and hypocentre location. Results show that the focal mechanism for the second event features a steeper dip angle and a strike rotated slightly clockwise with respect to the previous event. This kind of geological fault configuration, with secondary rupture in the hanging wall of a large normal fault, is commonly observed in extensional geological regimes. We propose that both earthquakes form part of a typical normal fault diverging splay, where the secondary fault connects to the main fault at depth. To ascertain more information on the spatial and temporal details of slip for both events, we gathered near-fault seismological and geodetic data. Through forward modelling of near-fault synthetic seismograms we build a kinematic k −2 earthquake source model with spatially distributed slip on the fault that, to first-order, explains both coseismic static displacement GPS vectors and short-period seismometer observations at the closest sites. As expected, the results for the first event agree with the focal mechanism derived from teleseismic modelling, with a magnitude M w 6.97. Similarly, near-fault modelling for the second event suggests rupture along a normal fault, M w 6.90, characterized by a steeper dip angle (dip = 74°) and a strike clockwise rotated (strike = 155°) with respect to the previous event.

Pichilemu↗

Dynamically downscaled climate simulations over North America: Methods, evaluation, and supporting documentation for users

We have completed an array of high-resolution simulations of present and future climate over Western North America (WNA) and Eastern North America (ENA) by dynamically downscaling global climate simulations using a regional climate model, RegCM3. The simulations are intended to provide long time series of internally consistent surface and atmospheric variables for use in climate-related research. In addition to providing high-resolution weather and climate data for the past, present, and future, we have developed an integrated data flow and methodology for processing, summarizing, viewing, and delivering the climate datasets to a wide range of potential users. Our simulations were run over 50- and 15-kilometer model grids in an attempt to capture more of the climatic detail associated with processes such as topographic forcing than can be captured by general circulation models (GCMs). The simulations were run using output from four GCMs. All simulations span the present (for example, 1968-1999), common periods of the future (2040-2069), and two simulations continuously cover 2010-2099. The trace gas concentrations in our simulations were the same as those of the GCMs: the IPCC 20th century time series for 1968-1999 and the A2 time series for simulations of the future. We demonstrate that RegCM3 is capable of producing present day annual and seasonal climatologies of air temperature and precipitation that are in good agreement with observations. Important features of the high-resolution climatology of temperature, precipitation, snow water equivalent (SWE), and soil moisture are consistently reproduced in all model runs over WNA and ENA. The simulations provide a potential range of future climate change for selected decades and display common patterns of the direction and magnitude of changes. As expected, there are some model to model differences that limit interpretability and give rise to uncertainties. Here, we provide background information about the GCMs and the RegCM3, a basic evaluation of the model output and examples of simulated future climate. We also provide information needed to access the web applications for visualizing and downloading the data, and give complete metadata that describe the variables in the datasets.

Open-File Report↗

The cyclic and fractal seismic series preceding an mb 4.8 earthquake on 1980 February 14 near the Virgin Islands

Seismic activity in the 10 months preceding the 1980 February 14, m b 4.8 earthquake in the Virgin Islands, reported on by Frankel in 1982, consisted of four principal cycles. Each cycle began with a relatively large event or series of closely spaced events, and the duration of the cycles progressively shortened by a factor of about 3/4. Had this regular shortening of the cycles been recognized prior to the earthquake, the time of the next episode of seismicity (the main shock) might have been closely estimated 41 days in advance. That this event could be much larger than the previous events is indicated from time-to-failure analysis of the accelerating rise in released seismic energy, using a non-linear time- and slip-predictable foreshock model. Examination of the timing of all events in the sequence shows an even higher degree of order. Rates of seismicity, measured by consecutive interevent times, when plotted on an iteration diagram of a rate versus the succeeding rate, form a triangular circulating trajectory. The trajectory becomes an ascending helix if extended in a third dimension, time. This construction reveals additional and precise relations among the time intervals between times of relatively high or relatively low rates of seismic activity, including period halving and doubling. The set of 666 time intervals between all possible pairs of the 37 recorded events appears to be a fractal; the set of time points that define the intervals has a finite, non-integer correlation dimension of 0.70. In contrast, the average correlation dimension of 50 random sequences of 37 events is significantly higher, close to 1.0. In a similar analysis, the set of distances between pairs of epicentres has a fractal correlation dimension of 1.52. Well-defined cycles, numerous precise ratios among time intervals, and a non-random temporal fractal dimension suggest that the seismic series is not a random process, but rather the product of a deterministic dynamic system.

Geophysical Journal International↗

User's Guide for Mixed-Size Sediment Transport Model for Networks of One-Dimensional Open Channels

This user's guide describes a mathematical model for predicting the transport of mixed sizes of sediment by flow in networks of one-dimensional open channels. The simulation package is useful for general sediment routing problems, prediction of erosion and deposition following dam removal, and scour in channels at road embankment crossings or other artificial structures. The model treats input hydrographs as stepwise steady-state, and the flow computation algorithm automatically switches between sub- and supercritical flow as dictated by channel geometry and discharge. A variety of boundary conditions including weirs and rating curves may be applied both external and internal to the flow network. The model may be used to compute flow around islands and through multiple openings in embankments, but the network must be 'simple' in the sense that the flow directions in all channels can be specified before simulation commences. The location and shape of channel banks are user specified, and all bedelevation changes take place between these banks and above a user-specified bedrock elevation. Computation of sediment-transport emphasizes the sand-size range (0.0625-2.0 millimeter) but the user may select any desired range of particle diameters including silt and finer (<0.0625 millimeter). As part of data input, the user may set the original bed-sediment composition of any number of layers of known thickness. The model computes the time evolution of total transport and the size composition of bed- and suspended-load sand through any cross section of interest. It also tracks bed -surface elevation and size composition. The model is written in the FORTRAN programming language for implementation on personal computers using the WINDOWS operating system and, along with certain graphical output display capability, is accessed from a graphical user interface (GUI). The GUI provides a framework for selecting input files and parameters of a number of components of the sediment-transport process. There are no restrictions in the use of the model as to numbers of channels, channel junctions, cross sections per channel, or points defining the cross sections. Following completion of the simulation computations, the GUI accommodates display of longitudinal plots of either bed elevation and size composition, or of transport rate and size composition of the various components, for individual channels and selected times during the simulation period. For individual cross sections, the GUI also allows display of time series of transport rate and size composition of the various components and of bed elevation and size composition.

Water-Resources Investigations Report↗

Uncertainty and extreme events in future climate and hydrologic projections for the Pacific Northwest: providing a basis for vulnerability and core/corridor assessments

The purpose of this project was to (1) provide an internally-consistent set of downscaled projections across the Western U.S., (2) include information about projection uncertainty, and (3) assess projected changes of hydrologic extremes. These objectives were designed to address decision support needs for climate adaptation and resource management actions. Specifically, understanding of uncertainty in climate projections – in particular for extreme events – is currently a key scientific and management barrier to adaptation planning and vulnerability assessment. The new dataset fills in the Northwest domain to cover a key gap in the previous dataset, adds additional projections (both from other global climate models and a comparison with dynamical downscaling) and includes an assessment of changes to flow and soil moisture extremes. This new information can be used to assess variations in impacts across the landscape, uncertainty in projections, and how these differ as a function of region, variable, and time period. In this project, existing University of Washington Climate Impacts Group (UW CIG) products were extended to develop a comprehensive data archive that accounts (in a reigorous and physically based way) for climate model uncertainty in future climate and hydrologic scenarios. These products can be used to determine likely impacts on vegetation and aquatic habitat in the Pacific Northwest (PNW) region, including WA, OR, ID, northwest MT to the continental divide, northern CA, NV, UT, and the Columbia Basin portion of western WY New data series and summaries produced for this project include: 1) extreme statistics for surface hydrology (e.g. frequency of soil moisture and summer water deficit) and streamflow (e.g. the 100-year flood, extreme 7-day low flows with a 10-year recurrence interval); 2) snowpack vulnerability as indicated by the ratio of April 1 snow water to cool-season precipitation; and, 3) uncertainty analyses for multiple climate scenarios.

Arizona, California, Colorado, Idaho, Montana, Nev↗

Geology of the Andover Granite and surrounding rocks, Massachusetts

Field and petrographic studies of the Andover Granite and surrounding rocks have afforded an opportunity for an explanation of its emplacement and crystallization. The investigation has contributed secondarily to an understanding of the geologic history of southeastern New England, particularly as it is revealed in the Lawrence, Wilmington, South Groveland, and Reading quadrangles of Massachusetts. The Andover Granite and Sharpners Pond Tonalite together comprise up to 90 percent of the Acadian(?) subalkaline intrusive series cropping out within the area of study. The subalkaline series locally invades a sequence of early to middle Paleozoic and possibly Precambrian metasedimentary and metavolcanic rocks. Much of the subalkaline series and most of the Andover Granite is confined between two prominent east-northeast trending faults or fault systems. The northern fault separates the mildly metamorphosed Middle Silurian(?) Merrimack Group on the north from a highly metamorphosed and thoroughly intruded Ordovician(?) sequence on the south. The southern 'boundary '' fault is a major structural discontinuity characterized by penetrative, diffuse shearing over a zone one-half mile or more in width. The magmatic nature of the Andover Granite is demonstrated by: (1) sharply crosscutting relationships with surrounding rocks; (2) the occurrence of tabular-shaped xenoliths whose long directions parallel the foliation within the granite and whose internal foliation trends at a high angle to that of the granite; (3) continuity with the clearly intrusive Sharpners Pond Tonalite; (4) the compositional uniformity of the granite as contrasted with the compositional diversity of the rocks it invades; (5) its modal and normative correspondence with (a) calculated norms of salic extrusives and (b) that of the ternary (granite) minimum for the system NaAlSi3O8-KAlSi3O8-SiO2. Orogenic granites, as represented by the Andover, contrast with post-orogenic granites, represented locally by the Peabody Granite, in their phase composition and texture. Unlike the Peabody, the Andover Granite is thought to have been thoroughly recrystallized through the unmixing of initially homogeneous phases with the concomitant development of extremely intricate, allotriomorphic textures. Textural relationships between potassium and plagioclase feldspars and among quartz and the two feldspars, suggest that the Andover Granite has evolved through exsolution of a single hypersolvus feldspar (or two coexisting subsolvus feldspars of only slightly disparate compositions) into discrete grains of plagioclase and potassium feldspar, much along the lines proposed by Tuttle (1952). A hypothesis is proposed for the origin of myrmekite whereby it is evolved indirectly through exsolution of a homogeneous, hypersolvus, calcalkali feldspar in the presence of a silica reservoir. Where the An 'molecule' is contained in the primary mix crystal, exsolution into potassium and plagioclase feldspar phases normally requires a paired exchange between Ca-Al and K-Si. Should the silicon requirements of the developing potassium feldspar be met by the matrix silica reservoir, the concomitantly evolving plagioclase may become stoichiometrically enriched in silicon and ultimately develop into myrmekite. Discrete unmixing of pure alkali feldspar proceeds through simple alkali ion exchange; ternary compostions high in An are more apt to fall initially in the two-feldspar field, thereby reducing the unmixing potential. General restriction of myrmekite to plagioclase of calcic albite to oligoclase composition is explained accordingly.

Open-File Report↗

Environmental considerations related to mining of nonfuel minerals

Throughout most of human history, environmental stewardship during mining has not been a priority partly because of the lack of applicable laws and regulations and partly because of ignorance about the effects that mining can have on the environment. In the United States, the National Environmental Policy Act of 1969, in conjunction with related laws, codified a more modern approach to mining, including the responsibility for environmental stewardship, and provided a framework for incorporating environmental protection into mine planning. Today, similar frameworks are in place in the other developed countries of the world, and international mining companies generally follow similar procedures wherever they work in the world. The regulatory guidance has fostered an international effort among all stakeholders to identify best practices for environmental stewardship. The modern approach to mining using best practices involves the following: (a) establishment of a pre-mining baseline from which to monitor environmental effects during mining and help establish geologically reasonable closure goals; (b) identification of environmental risks related to mining through standardized approaches; and (c) formulation of an environmental closure plan before the start of mining. A key aspect of identifying the environmental risks and mitigating those risks is understanding how the risks vary from one deposit type to another—a concept that forms the basis for geoenvironmental mineral-deposit models. Accompanying the quest for best practices is the goal of making mining sustainable into the future. Sustainable mine development is generally considered to be development that meets the needs of the present generation without compromising the ability of future generations to meet their own needs. The concept extends beyond the availability of nonrenewable mineral commodities and includes the environmental and social effects of mine development. Global population growth, meanwhile, has decreased the percentage of inhabitable land available to support society’s material needs. Presently, the land area available to supply the mineral resources, energy resources, water, food, shelter, and waste disposal needs of all Earth’s inhabitants is estimated to be 135 square meters per person. Continued global population growth will only increase the challenges of sustainable mining. Current trends in mining are also expected to lead to new environmental challenges in the future, among which are mine-waste management issues related to mining larger deposits for lower ore grade; water-management issues related to both the mining of larger deposits and the changes in precipitation brought about by climate change; and greenhouse gas issues related to reducing the carbon footprint of larger, more energy-intensive mining operations.

Professional Paper↗

Geology and mineral deposits of an area in the Departments of Antioquia and Caldas (Subzone IIB), Colombia

The Inventario Minero National (IMN), a four-year cooperative geologic mapping and mineral resources appraisal project, was accomplished under an agreement between the Republic of Colombia and the U. S. Agency for International Development from 1964 through 1969. Subzone IIB, consisting essentially of the east half of Zone comprises nearly 20,000 km2 principally in the Department of Antioquia but including also small parts of the Departments of Caldas and Tolima. The rocks in IIB range from Precambrian to Holocene. Precambrian feldspar-quartz gneiss occupies a mosaic of fault-bounded blocks intruded by igneous rocks between the Oto fault and the Rio Magdalena. Paleozoic rocks are extensive, and include lightly metamorphosed graptolite-bearing Ordovician shale at Cristalina, and a major suite of graphitic quartz-mica schist, feldspathic and aluminous gneiss, quartzite, marble, amphibolite, and other rocks. Syntectonic intrusive gneiss included many of the older rocks during a late Paleozoic(?) orogeny, which was accompanied by Abukuma-type metamorphosing from lowermost greenschist to upper amphibolite facies. A Jurassic diorite pluton bounded by faults cuts volcanic rocks of unknown age east of the Otu fault. Cretaceous rocks are major units. Middle Cretaceous carbonaceous shale, sandstone, graywacke, conglomerate, and volcanic rocks are locally prominent. The Antioquian batholith (quartz diorite) of Late Cretaceous age cuts the middle Cretaceous and older rocks. A belt of Tertiary nonmarine clastic sedimentary rocks crops out along the Magdalena Valley. Patches of Tertiary alluvium are locally preserved in the mountains. Quaternary alluvium, much of it auriferous, is widespread in modern stream valleys. Structurally IIB constitutes part of a vast complex synclinorium intruded concordantly by syntectonic catazonal or mesozonal felsic plutons, and by the later epizonal post-tectonic Antioquian batholith. Previously unrecognized major wrench faults are outstanding structural features of IIB. Some are traceable for several hundred kilometers and probably have displacements measurable in kilometers, although only the Palestina fault, with right-lateral displacement of 27.7 km, is accurately documented. Correlations of rocks mapped in IIB with those of outlying areas including neighboring IIA are discussed.

Open-File Report↗

Hydrology of the shallow aquifer and uppermost semiconfined aquifer near El Paso, Texas

The availability of fresh ground water in El Paso and adjacent areas that is needed to meet increased demand for water supply concerns local, State, and Federal agencies. The Hueco bolson is the principal aquifer in the El Paso area. Starting in the early 1900s and continuing to the 1950s, most of the municipal and industrial water supply in El Paso was pumped from the Hueco bolson aquifer from wells in and near the Rio Grande Valley and the international border. The Rio Grande is the principal surface-water feature in the El Paso area, and a major source of recharge to the shallow aquifer (Rio Grande alluvium) within the study area is leakage of flow from the Rio Grande. The shallow aquifer and the underlying Hueco bolson aquifer are in general hydrologic connection, but wells that penetrate these aquifers have different water levels and water quality. The configuration (slope) of the predevelopment (pre- 1903) water table indicates that ground water in the Hueco bolson flowed south from the Texas-New Mexico State line toward the El Paso Valley and the Rio Grande. There the water moved upward from the Hueco bolson fill through the alluvium and discharged as seepage at the land surface or through evapotranspiration on the river flood plain. Subsequent (post-January 1903) development of the Hueco bolson has caused water levels to decline in both the bolson aquifer and the shallow aquifer, and these declines have reversed the original vertical hydraulic gradient and now cause vertical leak- 1 Hydrologist, U.S. Geological Survey. 2 Geologist, El Paso Water Utilities-Public Service Board. age downward. Water-level declines in the shallow aquifer have in turn induced increasing amounts of leakage from the Rio Grande and from irrigation canals. The reversal from upward to downward in vertical hydraulic gradient between the Rio Grande alluvium and the underlying Hueco bolson aquifer has induced shallow water in the alluvium to move downward into the deeper aquifer. The introduction of water from the alluvium probably has led to a gradual water-quality deterioration of ground water in the Hueco bolson aquifer. The extent of any deterioration is a major concern because the dissolved solids concentration in water from some wells is approaching 1,000 milligrams per liter and already has exceeded this limit in other wells.

Texas↗

Sixth International Limnogeology Congress: abstract volume, Reno, Nevada, June 15-19, 2015

Limnogeology is the study of modern lakes and lake deposits in the geologic record. Limnogeologists have been active since the 1800s, but interest in limnogeology became prevalent in the early 1990s when it became clear that lake deposits contain continental environmental and climate records. A society that is focused on limnogeology would allow greater communication and access to research on these important subjects and contribute to providing sound science used to understand rapid global changes in our modern world; thus, the International Association of Limnogeology was founded in 1995 at the first International Limnogeology Congress (ILIC) held in Copenhagen, Denmark. The Sixth International Limnogeology Congress (ILIC6) was held in Reno, Nevada, from June 15&ndash;19, 2015. The ILIC meetings have been held every 4 years since the first meeting in1995 and were subsequently convened in Brest, France (1999), Tucson, Arizona, USA (2003), Barcelona, Spain (2007), and Konstanz, Germany (2011). The Congress in Reno, USA marks the second time the Congress has been held in the United States and more than 150 scientists from every part of the world participated. About one-half of the participants were from North America, together with scientists from Europe, South America, Asia, Africa, Australia, and New Zealand. The format of the Reno Congress followed the format originated at the Tucson Congress (ILIC3), which is unusual for scientific meetings. Nine keynote speakers spread throughout the Congress gave 1-hour talks, with the rest of the time available for viewing posters that were presented by the bulk of the participants. Keynote presentations were diverse and showed the breadth of research that is being done in lake systems worldwide. The abstracts of the keynote speakers and about 140 poster presentations are included in this volume. These posters cover a variety of limnologic, paleolimnologic, and limnogeologic topics including contaminant histories of lakes, the role of groundwater in lake processes, the formation of minerals in lake sediments, terminal lakes, how lakes reveal climate changes and paleohydrologic processes, the impact of volcanic emissions on lakes, as well as the biologic and chemical evolution of lake systems. The U.S. Geological Survey has sponsored each ILIC that has been held in the United States because of the importance of understanding paleoclimate and contaminant histories of lakes, two main themes of the Congress. This volume provides a permanent record of the wide variety of studies that are being conducted in modern lakes and ancient lake deposits worldwide, and it provides a stepping stone for any one desiring further discussion of the work that was presented at ILIC6.

Open-File Report↗

State of the Earth’s cryosphere at the beginning of the 21st century: Glaciers, global snow cover, floating ice, and permafrost and periglacial environments

This chapter is the tenth in a series of 11 book-length chapters, collectively referred to as &ldquo;this volume,&rdquo; in the series U.S. Geological Survey Professional Paper 1386, Satellite Image Atlas of Glaciers of the World. In the other 10 chapters, each of which concerns a specific glacierized region of Earth, the authors used remotely sensed images, primarily from the Landsat 1, 2, and 3 series of spacecraft, in order to analyze that glacierized region and to monitor changes in its glaciers. Landsat images, acquired primarily during the period 1972 through 1981, were used by an international team of glaciologists and other scientists to study the various glacierized regions and (or) to discuss related glaciological topics. In each glacierized region, the present distribution of glaciers within its geographic area is compared, wherever possible, with historical information about their past areal extent. The atlas provides an accurate regional inventory of the areal extent of glacier ice on our planet during the 1970s as part of an expanding international scientific effort to measure global environmental change on the Earth&rsquo;s surface. However, this chapter differs from the other 10 in its discussion of observed changes in all four elements of the Earth&rsquo;s cryosphere (glaciers, snow cover, floating ice, and permafrost) in the context of documented changes in all components of the Earth System. Human impact on the planet at the beginning of the 21st century is pervasive. The focus of Chapter A is on changes in the cryosphere and the importance of long-term monitoring by a variety of sensors carried on Earth-orbiting satellites or by a ground-based network of observatories in the case of permafrost. The chapter consists of five parts. The first part provides an introduction to the Earth System, including the interrelationships of the geosphere (cryosphere, hydrosphere, lithosphere, and atmosphere), the biosphere, climate processes, biogeochemical cycles, and the critically important hydrologic cycle, in which glacier ice is the second largest reservoir of water after the oceans. The second part assesses the state of glaciers in all of the glacierized regions of the planet, primarily as drawn in the other 10 chapters. It includes sections on ice cores and the climate record they contain, volumetric changes in glaciers, harnessing spaceborne sensors to measure changes in glaciers, and related topics. The third part summarizes trends in global snow cover. The fourth part summarizes long-term changes in area and thickness of floating ice, including polar sea ice and freshwater (lake and river) ice. The fifth part assesses the loss of permafrost and changes in periglacial environments at high latitudes and high altitudes.

Professional Paper↗