Search USGSSearch

SEARCH · Search USGS

Results for “Final Report”

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

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

At least 1,189 records · Page 66Linked to original sources

Population genetic structure of rare and endangered plants using molecular markers

This study was initiated to assess the levels of genetic diversity and differentiation in the remaining populations of Phyllostegia stachyoides and Melicope zahlbruckneri in Hawai`i Volcanoes National Park and determine the extent of gene flow to identify genetically distinct individuals or groups for conservation purposes. Thirty-six Amplified Fragment Length Polymorphic (AFLP) primer combinations generated a total of 3,242 polymorphic deoxyribonucleic acid (DNA) fragments in the P. stachyoides population with a percentage of polymorphic bands (PPB) ranging from 39.3 to 65.7% and 2,780 for the M. zahlbruckneri population with a PPB of 18.8 to 64.6%. Population differentiation (Fst) of AFLP loci between subpopulations of P. stachyoides was low (0.043) across populations. Analysis of molecular variance of P. stachyoides showed that 4% of the observed genetic differentiation occurred between populations in different k ī puka and 96% when individuals were pooled from all k ī puka. Moderate genetic diversity was detected within the M. zahlbruckneri population. Bayesian and multivariate analyses both classified the P. stachyoides and M. zahlbruckneri populations into genetic groups with considerable sub-structuring detected in the P. stachyoides population. The proportion of genetic differentiation among populations explained by geographical distance was estimated by Mantel tests. No spatial correlation was found between genetic and geographic distances in both populations. Finally, a moderate but significant gene flow that could be attributed to insect or bird-mediated dispersal of pollen across the different k ī puka was observed. The results of this study highlight the utility of a multi-allelic DNA-based marker in screening a large number of polymorphic loci in small and closely related endangered populations and revealed the presence of genetically unique groups of individuals in both M. zahlbruckneri and P. stachyoides populations. Based on these findings, approaches that can assist conservation efforts of these species are proposed.

Hawai`i

Adaptive harvest management for the Svalbard population of pink-footed geese: cooperator report

This document describes progress to date on the development of a harvest‐management strategy for maintaining pink‐footed goose abundance near their target level by providing for sustainable harvests in Norway and Denmark. Many goose populations in western Europe have increased dramatically in recent decades. The Svalbard population of pink‐footed geese (Anser brachyrhynchus) is a good example, increasing from about 10 thousand individuals in the early 1960’s to roughly 80 thousand today. Although these geese are a highly valued resource, the growing numbers of geese are causing agricultural conflicts in wintering and staging areas. The African‐Eurasian Waterbird Agreement (AEWA; http://www.unep‐aewa.org/) calls for means to manage populations which cause conflicts with certain human economic activities. We compiled relevant demographic and weather data and specified an annual‐cycle model for pink-footed geese that reconciles the different dates of monitoring activities and the timing of harvest-management decisions. We then developed dynamic models for survival and reproductive processes and parameterized them using available data. By combining varying hypotheses about survival and reproduction, we developed a suite of nine models that represent a wide range of possibilities concerning the extent to which demographic rates are density dependent or independent, and the extent to which spring temperatures are important. These nine models varied significantly in their predictions of the harvest required to stabilize current population size, ranging from a low of about 500 to a high of about 17 thousand. For comparison, the harvest in Norway and Denmark was about 11 thousand in 2011 and the population increased from 70 to 80 thousand. We relied on the passive form of adaptive management in formulating a harvest strategy. In passive adaptive management, alternative population models and their associated weights of evidence are explicitly considered in the development of an optimal harvest strategy. Unlike active adaptive management, however, there is no explicit consideration of how harvest management actions could reduce uncertainty as to the most appropriate model of population dynamics. In optimizing a harvest strategy, we assumed equal probabilities for all nine models and assumed relatively course control over harvest. We used a management objective that seeks to maximize sustainable harvest, but avoids harvest decisions that are expected to result in a subsequent population size different than the population goal of 60 thousand. Optimal harvest strategies were calculated using stochastic dynamic programming, and Monte Carlo simulations were used to investigate expected strategy performance. The optimal passive adaptive‐management strategy is expected to maintain mean population size near 60 thousand, regardless of the most appropriate model. However, mean harvest rates and harvests varied substantially depending on the most appropriate model of population dynamics. With an average number of days above freezing in May in Svalbard, optimal harvest rates (i.e., the proportion of the population to be harvested in autumn) increase rapidly once there are more than about 50 thousand birds in the population. Generally, optimal harvests were on the order of 10 – 20 thousand for population sizes > 60 thousand, and 0 – 5 thousand for population sizes < 60 thousand. For the observations of young of 15.4 thousand and adults of 54.6 thousand in autumn 2010, and 10 days above freezing in May 2011 (a relatively warm spring compared to the average of about 7), the optimal harvest rate in autumn of 2011 would have been 0.16, or a harvest of about 14 thousand. Based on the optimal strategy, hunting‐season closures would be required as the number of adults in the autumn population falls below about 52 thousand, regardless of the number of young in the population. As the number of adults and young decrease, the number of warm days in May required to keep the hunting season open increases. We also investigated the ability of the optimal strategy to stabilize the population at around 60 thousand birds, assuming varying values of the maximum harvest rate that could be implemented. Harvest strategies that contained a maximum harvest rate of 0.16 (equivalent to a harvest of about 17 thousand) were effective at stabilizing the population at 60 thousand within 4‐5 years, regardless of climate scenario. Harvest strategies with a maximum harvest rate of 0.12 (harvest ≈ 13 thousand) were also able to stabilize the population near 60 thousand, although it took more time. Harvest strategies with a maximum harvest rate of 0.08 (harvest ≈ 8 thousand) were unsuccessful at stabilizing the population at 60 thousand. Continued monitoring of the pink‐footed goose population on an annual basis is critical to an informed harvest management strategy. At a minimum, the ground census in November should be continued to determine population size and proportion of young. Continued estimates of harvest from Norway and Denmark are also necessary to help judge the credibility of the alternative population models. However, an adaptive management process that relies on periodic updating of model weights will depend on acquiring either estimates of the realized harvest rate of adults or the age composition of the harvest. We also recommend that a census conducted during spring migration be operationalized, and that estimates of survival based on mark‐recapture data be updated. Finally, the International Working Group has expressed a desire to adopt a three‐year cycle of decision making related to the regulation of pink‐footed goose harvests. The idea is that once a target harvest level is adopted, it would remain in place for three years, after which time population status would be assessed and a potentially new management action chosen. We have developed a preliminary framework to implement a three‐year cycle using stochastic dynamic programming, and we hope to have it fully operational later this year . We note, however, that application of this 3‐year framework will still require annual resource monitoring and assessments to facilitate learning, and to allow managers the opportunity to respond to any unforeseen change in resource conditions.

Report

Perched water tables on hillsides in western Oregon: II. Preferential downslope movement of water and anions

Perched water tables on hillsides in western Oregon potentially provide a means by which pollutants from agricultural and domestic sources may enter surface waters and consequently degrade the quality of these waters. This paper reports the results of experiments which were carried out to investigate the flow of solutes and water from buried line sources in and above perched water tables on three different hillsides in western Oregon. In the saturated soil and rock mantles of two hillsides, evidence suggests that water and solutes flow preferentially through large continuous voids. The evidence includes the very high maximum rate of anion movement downslope, the failure of pore velocities calculated from Darcy's Law to predict the rate and direction of anion movement, the short time before high levels of anions are detected downslope, and the spatial variability in the maximum rates of anion movement over the depth and width of sites. On the third hillside, tracer anion movement and soil water pressure measurements indicate that preferential flow occurs through large continuous pores during heavy rainfall while the surrounding soil and rock mantle remains unsaturated. Finally, experiments with Rodamine WT dye show that a variety of continuous voids of different origins and sizes conduct water and solutes preferentially in the soil and upper regions of rock on these hillsides.

Oregon

Terrestrial Ecosystems - Topographic Moisture Potential of the Conterminous United States

As part of an effort to map terrestrial ecosystems, the U.S. Geological Survey has generated topographic moisture potential classes to be used in creating maps depicting standardized, terrestrial ecosystem models for the conterminous United States, using an ecosystems classification developed by NatureServe. A biophysical stratification approach, developed for South America and now being implemented globally, was used to model the ecosystem distributions. Substrate moisture regimes strongly influence the differentiation and distribution of terrestrial ecosystems, and therefore topographic moisture potential is one of the key input layers in this biophysical stratification. The method used to produce these topographic moisture potential classes was based on the derivation of ground moisture potential using a combination of computed topographic characteristics (CTI, slope, and aspect) and mapped National Wetland Inventory (NWI) boundaries. This method does not use climate or soil attributes to calculate relative topographic moisture potential since these characteristics are incorporated into the ecosystem model though other input layers. All of the topographic data used for this assessment were derived from the USGS 30-meter National Elevation Dataset (NED ) including the National Compound Topographic Index (CTI). The CTI index is a topographically derived measure of slope for a raster cell and the contributing area from upstream raster cells, and thus expresses potential for water flow to a point. In other words CTI data are 'a quantification of the position of a site in the local landscape', where the lowest values indicate ridges and the highest values indicate stream channels, lakes and ponds. These CTI values were compared to independent estimates of water accumulation by obtaining geospatial data from a number of sample locations representing two types of NWI boundaries: freshwater emergent wetlands and freshwater forested/shrub wetlands. Where these shorelines (the interface between the NWI wetlands and adjacent land) occurred, the CTI values were extracted and a histogram of their statistical distributions was calculated. Based on an evaluation of these histograms, CTI thresholds were developed to separate periodically saturated or flooded land, mesic uplands (moderately moist), and uplands. After the range of CTI values for these three different substrate moisture regimes was determined, the CTI values were grouped into three initial topographic moisture potential classes. As a final step in the generation of this national data layer, the uplands classification was subdivided into either very dry uplands or dry uplands. Very dry uplands were defined as uplands with relatively steep, south-facing slopes, and identification of this class was based on the slope and aspect datasets derived from the NED. The remaining uplands that did not meet these additional criteria were simply re-classified as dry uplands. The final National Topographic Moisture Potential dataset for the conterminous United States contains four classes: periodically saturated or flooded land (CTI = 18.5), mesic uplands (12 =< CTI < 18.5), dry uplands (CTI < 12), and very dry uplands (CTI < 12, Slope > 24 degrees and 91 degrees =< Aspect =< 314 degrees). This map shows a smoothed and generalized image of the four topographic moisture potential classes. Additional information about this map and any of the data developed for the ecosystems modeling of the conterminous United States is available online at http://rmgsc.cr.usgs.gov/ecosystems/.

Scientific Investigations Map

Reconstructing the Quaternary depositional history using geologic mapping and three-dimensional modeling of the subsurface near Fort Morgan, northeastern Colorado

Centered on Fort Morgan, Colorado, this study is intended to build from previous work by adding a three-dimensional (3D) view of the subsurface to better understand the depositional history of Quaternary deposits. A 1:100,000 scale geologic map was made by combining previous geologic maps, regional soil maps, and recent field investigations. In addition to the geologic mapping, drill hole lithologic data from water wells and oil and gas exploration were compiled and lithologic units simplified to best represent the stratigraphy of the Quaternary deposits. From these subsurface data, a 3D subsurface model was constructed, trimmed at the surface by a digital elevation model, and a bedrock surface foundation gridded from drill hole data was added. The surface of the 3D model was then compared visually to the surficial geologic map. Cross sections were constructed from the 3D model and compared to site-specific drilling that was done as part of this project. Finally, the model was examined in detail to reconstruct the depositional history of the subsurface alluvial and eolian units. Alluvial and fluvial drainage basins exposed in the subsurface have a greater areal extent than the present-day narrow drainages. Older eolian sand in the subsurface tends to be interbedded with loess indicating coeval deposition. Holocene sand, both eroded from bedrock exposed at the surface north of the study area and reworked from the South Platte River, buries most of the interbedded older sand and loess.

Colorado

Guidelines for monitoring and adaptively managing restoration of Chinook salmon ( Oncorhynchus tshawytscha ) and steelhead ( O. mykiss ) on the Elwha River

As of January, 2014, the removal of the Elwha and Glines Canyon dams on the Elwha River, Washington, represents the largest dam decommissioning to date in the United States. Dam removal is the single largest step in meeting the goals of the Elwha River Ecosystem and Fisheries Restoration Act of 1992 (The Elwha Act) &mdash; full restoration of the Elwha River ecosystem and its native anadromous fisheries (Section 3(a)). However, there is uncertainty about project outcomes with regards to salmon populations, as well as what the &lsquo;best&rsquo; management strategy is to fully restore each salmon stock. This uncertainty is due to the magnitude of the action, the large volumes of sediment expected to be released during dam removal, and the duration of the sediment impact period following dam removal. Our task is further complicated by the depleted state of the native salmonid populations remaining in the Elwha, including four federally listed species. This situation lends itself to a monitoring and adaptive management approach to resource management, which allows for flexibility in decision-making processes in the face of uncertain outcomes. The Elwha Monitoring and Adaptive Management (EMAM) guidelines presented in this document provide a framework for developing goals that define project success and for monitoring project implementation and responses, focused upon two federally listed salmon species &mdash; Puget Sound Chinook salmon (Oncorhynchus tshawytscha) and Puget Sound steelhead (O. mykiss). The framework also should serve as a guide to help managers adaptively manage fish restoration actions during and following dam removal. The document is organized into seven sections, including an introduction (Section 1), a description of the adaptive management approach (Section 2), suggested modifications to the existing restoration strategy developed in previous Elwha River restoration documents (section 3), specific descriptions of an adaptive management framework, including establishment of goals, performance indicators, and potential adaptive management responses to monitoring information (section 4), monitoring tools and methods for use in evaluating performance and project outcomes (section 5), and brief sections on data record keeping and reporting (Section 6) and an estimated budget (section 7). The purpose of the EMAM guidelines is to propose (1) refinement of existing goals established in previous documents (e.g., Ward et al. (2008), U.S. Department of the Interior, Department of Commerce, and Lower Elwha S&rsquo;Klallam Tribe (1994)); (2) an adaptive management framework, (3) specific trigger values for relevant performance indicators that guide the adaptive management approach, (4) a specific monitoring strategy for evaluating outcomes of restoration activities; (5) a data management strategy, (6) information needed for adjusting goals when observations indicate conditions are different from anticipated. When taken together, our proposed adaptive management guidelines rely upon setting goals and objectives for each species of interest, which are monitored by relevant performance indicators and measurable trigger values that define success within each phase of the project. The guidelines themselves are arranged in a hierarchy for each species of interest. The levels of this hierarchy are goals, objectives, performance indicators, decision rules, triggers, and decisions (i.e., management/policy response). The monitoring and adaptive management approach provided is based on monitoring several categories of performance indicators, each containing associated &lsquo;trigger&rsquo; values which, when met, alters restoration activities (e.g., hatchery releases and/or strategies) through four successive restoration phases. Performance indicators proposed in these EMAM guidelines are based upon Viable Salmon Population (VSP) metrics, including abundance, productivity, distribution, and diversity (McElhany et al. 2000). Trigger values for each performance indicator are developed for four different restoration phases: Preservation, Recolonization, Local Adaptation, and Viable Natural Population. These biologically-based phases each have a set of objectives that are based on resource management scenarios, including the dam removal project itself, which change largely based on the level of active management required and the degree, if any, of resource utilization. Thus, details of prescribed management actions for each phase are based upon different needs specific to that phase. The creation of biologically-based phases is one of the major differences between our proposed EMAM guidelines and previously presented plans for Elwha River Restoration Project management. Changed largely in response to the recommendations of the most recent of three Hatchery Scientific Review Group project reviews (HSRG 2012), the goal-oriented phases replaced the previous system of temporal changes centered around the decommissioning of the dams (i.e., before, during, and after dam removal). By focusing on outcomes associated with rebuilding salmon populations instead of an engineering schedule, the guidelines are more amenable to an adaptive management framework and the ability for management actions to influence outcomes, particularly in the periods during and following dam removal. Trigger values for each performance indicator were generally developed using existing data from the Elwha River watershed, the Puget Sound region, or other Pacific Northwest rivers (i.e., elsewhere in Washington State, Oregon, British Columbia) modified to be relevant for Chinook salmon and steelhead recovery in the Elwha River. By meeting all of the trigger value levels for all performance indicators for a set amount of time within a management phase, the guidelines call for moving to the next phase. This next phase has a new set of trigger values for the same performance indicators. For example, upon moving from the Preservation phase to the Recolonization phase, the trigger value for intrinsic potential increases. Intrinsic potential is a pre-defined estimate of the total extent of available habitat within a watershed for adult and juvenile fish, specific to the target species and is therefore a performance indicator of spatial distribution. By the final Viable Natural Population phase, the entire intrinsic potential of the watershed is being occupied by the species of interest. For those cases when a performance indicator is not exceeding the target value for a particular phase after a certain time period, the trigger values provided in this document, as well as a series of exogenous variables, are explored that may help explain why the performance indicator is not being met. These exogenous variables include variables that are not part of the suite of performance indicators, such as hatchery production, harvest, habitat, and ecosystem indicators. In these cases where the program is stuck in a particular recovery phase, the situation could be caused by the selection of inappropriate trigger values or unforeseen environmental conditions. If the former, adaptive management would call for existing monitoring data to be used for modifying trigger values to an appropriate level. If one of the exogenous variables is found to be preventing the program moving to the next phase, then appropriate changes to management would be advised. For each performance indicator and many of the exogenous variables, a set of monitoring tools were proposed. Data standards were also proposed for data generated by each monitoring tool. Data management, record keeping, and reporting of monitoring and adaptive management activities and results are also outlined. Management of data from the focused monitoring program and documenting the outcomes of trigger value evaluations and associated decisions from the adaptive management approach are key components of the EMAM guidelines. Without a clear history of data generated and adaptive management decisions taken by managers, the ability to learn through adaptive management breaks down. In addition to the long time period involved, another complication is the fact that the data will likely be collected by different federal and state agencies, tribal staff, and others. Having a system of reporting developed should help alleviate potential problems. The restoration of the migration route to spawning and rearing habitats upstream of the former Glines Canyon Dam represents a great opportunity for salmon on the Olympic Peninsula. By removing two aging structures, it will be possible for all 5 species of salmon and steelhead to return to wild stretches of the Elwha River and major floodplain habitat characterized by multiple channels, as well as significant portions of numerous tributaries. Measuring the progress of restoration, from the perspective of both salmon populations and the ecosystem upon which they depend, is a great test for a collaborative team of scientists. The normally challenging conditions of working in a steep gradient, high velocity wilderness river are exacerbated by the release of millions of cubic yards of sediment that had accumulated in the reservoirs. After the first two years of the dam decommissioning process, this release has changed the ecology of the river, estuary, and nearshore habitats downstream of the dams. Our goal in developing the guidelines described is to provide a roadmap for tracking what hopefully will become a successful outcome. If successfully implemented, this information should prove useful as others begin planning for the removal, alteration, or reconstruction of dams throughout North America and elsewhere, an inevitable outcome of an aging dam infrastructure.

Washington

Southern sea otter (Enhydra lutris nereis) population biology at Big Sur and Monterey, California --Investigating the consequences of resource abundance and anthropogenic stressors for sea otter recovery

The range of the southern sea otter ( Enhydra lutris nereis ) spans most of the central California coast from Half Moon Bay to Gaviota. Some coastal areas within this range are heavily developed and highly impacted by humans, while other areas are wild and largely pristine. Determining the relative importance of food resource abundance, environmental conditions, and anthropogenic increases in pathogens and pollutants to population change in sea otters is critical to understanding limitations to population growth. To investigate the causal links between the sluggish population growth of sea otters in central California and factors that could be driving variation in survival and reproduction, we designed a study to compare two distinct subpopulations—one in an area of low human impact (Big Sur) and one in an area of high human impact (Monterey). Between 2008 and 2011, the U.S. Geological Survey and collaborators conducted a telemetry-based study of sea otters at these two locations. The results of this study were not consistent with the hypothesis that sea otters adjacent to human population centers (Monterey) experience higher exposure to pollutants and pathogens than those in lower impacted areas (Big Sur). In fact, based on serological analysis, female sea otters from Big Sur showed higher exposure rates to Toxoplasma gondii than did female otters from Monterey, while domoic acid exposure appeared to be similar at both sites. Gene expression (specifically transcription) analysis did not indicate any consistent differences between the two populations that would have suggested a response to pathogen or toxin exposure, although there were temporal changes in gene transcription for sea otters at Big Sur following potential exposure to run-off from wildfires that occurred during the study. Together, these metrics suggest that variation in exposure to environmental stressors occurred, but patterns were not clearly attributable to differences in human population densities or land-use patterns. When compared to Monterey, sea otters in Big Sur spent more time feeding, had a higher degree of dietary specialization, were in poorer body condition, and had lower survival rates (both pups and adults). Together, these metrics suggest that otters at Big Sur had greater nutritional stress, consistent with lower per-capita resource abundance. Overall, study results indicate that density-dependent population regulation, mediated by per-capita resource abundance, is the most significant factor currently limiting population growth in the center part of the range. Additionally, spatial and temporal variation in environmental and anthropogenic stressors also can affect sea otter health, although patterns of variation are complex and are not simply a function of proximity to human populations. We also found that exposure to environmental stressors (either natural or anthropogenic in origin) often is associated with resource limitation. Finally, our results indicate that sea otter populations are structured at relatively small spatial scales, and the processes that regulate population abundance (including density-dependent resource abundance) also occur at these smaller, more local scales.

California

Conceptualization and analysis of ground-water flow system in the Coastal Plain of Virginia and adjacent parts of Maryland and North Carolina

The ground-water flow system in the Coastal Plain of Virginia and adjacent parts of Maryland and North Carolina consists of a water table aquifer and an underlying sequence of confined aquifers and intervening confining units composed of unconsolidated sand and clay. A digital flow model was developed to enhance knowledge of the behavior of the ground-water flow system in response to its development. Ten pumping periods covering 90 yr of withdrawal simulated the history of ground-water development. Simulated potentiometric-surface maps for 1980 show lowered water levels and the development of coalescing cones of depression around the cities of Franklin, Suffolk, and Williamsburg and the town of West Point, all in Virginia. The largest simulated decline in water level, about 210 ft was near Franklin. Water budgets indicate that over the period of simulation (1891-1980): (1) pumpage from the model area increased by about 105 Mgal/d; (2) lateral boundary outflow increased by about 5 Mgal/d; (3) ground-water flow to streams and coastal water decreased by about 107.5 Mgal/d; (4) lateral boundary inflow increased by about 0.7 Mgal/d, and (5) water released from aquifer storage increased by about 1.6 Mgal/d. Simulated rates of recharge into the confined aquifer system at the end of the final pumping period (1980) varied up to 3.8 in/yr. and simulated rates of discharge out of the confined system varied up to 2.2 in/yr. Results of simulations show an increase of about 110 Mgal/d into the confined system from the unconfined system over the period of simulation. This increase in flow into the confined system affected local discharge of ground water to streams and regional discharge to coastal water. Lowering the storage coefficient of the aquifer had a minimal effect simulated water levels, whereas increasing the storage coefficient had a much more significant effect.

Maryland, North Carolina, Virginia

To burn or not to burn Oriental bittersweet: A fire manager’s conundrum

Oriental bittersweet ( Celastrus orbiculatus ) is a highly invasive liana (woody vine) that occurs throughout the Eastern United States. This twining plant can blanket and girdle adjacent vegetation, affecting succession and damaging trees. In areas where prescribed fire is a management tool, the response of Oriental bittersweet to fire needs to be quantified, rather than relying on anecdotal evidence. Currently, in areas already infested with this species, there are no strategies for prioritizing pre- or post-fire treatments on Oriental bittersweet. This largely results from a lack of understanding of the nature of post-fire resprouting by this species. Sprouting of bittersweet can at least double with fire and sprouts appear to respond to fire with an increase in growth rate (Pavlovic and Young pers. obs.). Beyond this basic need to understand the interaction between fire and Oriental bittersweet resprouting, we need to investigate how fire may interact with light, soil moisture, litter and other environmental factors to either increase or decrease abundance of this species. Finally, it is unknown how fire regimes influence the distribution of Oriental bittersweet on the landscape; thus we need to model the distribution of Oriental bittersweet in a fire impacted landscape. If we determine through our research that fire enhances the spread of this species, modification of fire suppression tactics and potential fire exclusion zones may be necessary. Thus we will be able to provide land managers throughout the Eastern US with data-driven decision support tools for more successful management of this species in fire dependent and invaded areas.

Indiana

An interactive code (NETPATH) for modeling NET geochemical reactions along a flow PATH, version 2.0

NETPATH is an interactive Fortran 77 computer program used to interpret net geochemical mass-balance reactions between an initial and final water along a hydrologic flow path. Alternatively, NETPATH computes the mixing proportions of two to five initial waters and net geochemical reactions that can account for the observed composition of a final water. The program utilizes previously defined chemical and isotopic data for waters from a hydrochemical system. For a set of mineral and (or) gas phases hypothesized to be the reactive phases in the system, NETPATH calculates the mass transfers in every possible combination of the selected phases that accounts for the observed changes in the selected chemical and (or) isotopic compositions observed along the flow path. The calculations are of use in interpreting geochemical reactions, mixing proportions, evaporation and (or) dilution of waters, and mineral mass transfer in the chemical and isotopic evolution of natural and environmental waters. Rayleigh distillation calculations are applied to each mass-balance model that satisfies the constraints to predict carbon, sulfur, nitrogen, and strontium isotopic compositions at the end point, including radiocarbon dating. DB is an interactive Fortran 77 computer program used to enter analytical data into NETPATH, and calculate the distribution of species in aqueous solution. This report describes the types of problems that can be solved, the methods used to solve problems, and the features available in the program to facilitate these solutions. Examples are presented to demonstrate most of the applications and features of NETPATH. The codes DB and NETPATH can be executed in the UNIX or DOS1 environment. This report replaces U.S. Geological Survey Water-Resources Investigations Report 91-4078, by Plummer and others, which described the original release of NETPATH, version 1.0 (dated December, 1991), and documents revisions and enhancements that are included in version 2.0. 1 The use of trade, brand or product names in this report is for identification purposes only and does not constitute endorsement by the U.S. Geological Survey.

Water-Resources Investigations Report

Surficial geology

Surficial materials are those at or near the Earth's surface. They constitute, by far, the largest and most used part of the ground around us. Areas not covered by surficial deposits--bare bedrock--form probably less than 5 percent of our land surface. Most surficial deposits are composed of poorly consolidated clay, silt, sand, or gravel-sized particles that are produced chiefly by erosion and are transported by and finally deposited by water, wind, or ice, but are also partly produced by the in-situ weathering of bedrock. The major genetic categories are shown on the map by various colors; the principal compositional types in each category are indicated by patterns. Colors, patterns, and letter symbols are given in the explanation which follows this summary description. Surficial deposits have characteristics important to our environment--water-bearing properties, mineral resources, and suitability as a natural foundation for buildings. They are also susceptible to flooding, erosion, ground subsidence, land slides, and earthquakes.

National Atlas of the United States of America

Water resources of the Mississippi Headwaters Watershed, North-central Minnesota

The Mississippi Headwaters watershed is a 7,608 square mile area in north-central Minnesota which includes all land drained by the Mississippi River above the Crow Wing River. From its source in Lake Itasca, 1,460 feet above mean sea level, the Mississippi River follows a semicircular 376 mile course to where it leaves the watershed at an altitude of 1,150 feet. The origin of the river is in glacial moraine. The river flows through end moraines and rolling till plains, across an extensive outwash plain occupied by large reservoir lakes, and finally across an extensive marshy plain which is the bed of an ancient glacial lake. Glacial deposits in the watershed include till, lenses of sand and gravel in till, outwash deposits of sand and gravel, and lake deposits of fine sand, silt, and clay. Beneath the glacial drift is an uneven surface of Precambrian bedrock. Bedrock consists of igneous intrusives and metamorphosed sedimentary formations. Water resources in the area are abundant. Lakes and streams occupy about 8 percent of the surface of the area. Ground-water supplies are available from glacial drift and from bedrock.

Minnesota

Geologic maps of the Descartes region of the Moon Apollo 16 pre-mission maps

The evolution of the Moon’s s urface shows three successive stages (Hartmann, 1970). The earliest , from the beginning of the decipherable record to the Imbrian P eriod, was dominated by frequent large -scale impacts which produced overlapping and interfering craters and g iant multi-ringed basin with ejecta blankets. The second stage , extending through most of the Imbrian to early Eratosthenian Periods , was characterized by volcanism, producing basalt flows of the maria an vari ed materials on the terrae . The final and longest stage , extending through the greater part of the Eratosthe n ian Period and the whole of the Copernican , was one of reduced geologic activity. Major modifications of the geologic structure by large impacts of volcanism occurring only locally, although changes by small-scale impact and mass wasting affected the entire surface. The De s cartes region, in the central highlands of the Moon about 500 km southeast of the center of the lunar disk, is dominated by features of the middle st a ge and contains a n exceptionally favorable site for their study.

IMAP

Evaluation of the ground-water resources of parts of Lancaster and Berks Counties, Pennsylvania

Secondary openings in bedrock are the avenues for virtually all ground-water flow in a 626-sqare-mile area in Lancaster and Berks Counties, Pennsylvania. The number, size, and interconnection of secondary openings are functions of lithology, depth, and topography. Ground water actively circulates to depths of 150 to 300 feet below land surface. Total average annual ground-water recharge for the area is 388 million gallons per day, most of which discharges to streams from local, unconfined flow systems. A digital ground-water flow model was developed to simulate unconfined flow under several different recharge and withdrawal scenarios. On the basis of lithologic and hydrologic differences, the modeled area was sub-divided into 22 hydrogeologic units. A finite-difference grid with rectangular blocks, each 2,015 by 2,332 feet, was used. The model was calibrated under steady-state and transient conditions. The steady-state calibration was used to determine hydraulic conductivities and stream leakage coefficients and the transient calibration was used to determine specific yields. The 22 hydrogeologic units fall into four general lithologies: Carbonate rocks, metamorphic rocks, Paleozoic sedimentary rocks, and Triassic sedimentary rocks. Average hydraulic conductivity ranges from about 8.8 feet per day in carbonate units to about .5 feet per day in metamorphic units. The Stonehenge Formation (limestone) has the greatest average hydraulic conductivity--85.2 feet per day in carbonate units to about 0.11 feet per day in the greatest gaining-strem leakage coefficient--16.81 feet per day. Specific yield ranges from 0.06 to 0.09 in carbonate units, and is 0.02 to 0.015, and 0.012 in metamorphic, Paleozoic sedimentary, and Triassic sedimentary units, respectively. Transient simulations were made to determine the effects of four different combinations of natural and artificial stresses. Natural aquifer conditions (no ground-water withdrawals) and actual aquifer conditions (current ground-water withdrawals) were simulated for two years under normal seasonal and hypothetical drought (60-percent reduction in winter-spring recharge) conditions. In October, 6 months after the hypothetical drought, simulated declines in water-table altitude due to the drought occurred everywhere and ranged from a median of 3.6 feet in carbonate units to 8.7 feet in carbonate units. Simulated base flows for five major streams were reduced by 33 to 51 percent during the hypothetical drought. Also in October, maximum simulated declines in water-table altitude due to ground-water withdrawls ranged from 33 feet in carbonate units to 79 feet in Triassic sedimentary units. Simulated base flows for five major streams were reduced by the amount of ground water withdrawn. Finally, again in October, maximum simulated declines in water-table altitude due to the combination of hypothetical drought and ground-water withdrawls ranged from 38 feet in carbonate units to 109 feet in Triassic sedimentary units. Due to aquifer dewatering, simulated declines were as much as 24 feet greater than the sum of the separate simulated declines that were caused by hypothetical drought and ground-water withdrawals. Some of the greatest simulated declines were in well fields, operated by three municipalities that experienced water-supply problems during the 1980-81 drought.

Water-Resources Investigations Report

Evidence for Paleocene dinosaurs in the Ojo Alamo Sandstone, San Juan Basin, New Mexico

Dinosaur fossils from the Ojo Alamo Sandstone, San Juan Basin, New Mexico, have been a puzzle to geologists for decades because other stratigraphic and paleontologic evidence indicates that this formation is Paleocene in age. The Paleocene age of the Ojo Alamo is based on the following evidence: (1) the lowermost part of the Nacimiento Formation, which conformably overlies and intertongues with the Ojo Alamo Sandstone in the Ojo Alamo type area in the southern part of the basin, contains Paleocene (Puercan) mammalian fossils; (2) plant megafossils collected from the Ojo Alamo indicate a Tertiary age for this rock unit; (3) pollen-productive rock samples collected from the Ojo Alamo at several localities in the basin have, without exception, yielded Paleocene-age palynomorph assemblages, and conversely, no Cretaceous pollen has ever been reported from Ojo Alamo rock samples. We report here the first discovery of Paleocene palynomorphs from Ojo Alamo Sandstone outcrops in the northern San Juan Basin. The pollen-producing rock samples were collected from a coaly, carbonaceous shale bed 12-m above the base of the Ojo Alamo and 3-m below the level of a large hadrosaur femur (1310-mm long and weighing about 130 kg) collected from this site. The bone was recovered from a cobble-to-boulder conglomeratic sandstone but shows no evidence of abrasion by erosion indicating almost instant burial after death of the animal and early mineralization of the femur. It is our opinion that this bone was not reworked from underlying Cretaceous strata. The data from the San Juan River site indicate that dinosaurs did not die out at the end of the Cretaceous in the San Juan Basin, but lived on for as much as a few hundred thousand years into the Paleocene before finally becoming extinct.

New Mexico

Petroleum systems of the Malay Basin Province, Malaysia

The offshore Malay Basin province is a Tertiary oil and gas province composed of a complex of half grabens that were filled by lacustrine shales and continental clastics.These deposits were overlain by clastics of a large delta system that covered the basin.Delta progradation was interupted by transgressions of the South China Sea to the southeast, which finally flooded the basin to form the Gulf of Thailand.Oil and gas from the Oligocene to Miocene lacustrine shales and Miocene deltaic coals is trapped primarily in anticlines formed by inversion of the half grabens during the late Miocene.Hydrocarbon reserves that have been discovered amount to 12 billion barrels of oil equivalent.The U.S. Geological Survey assessment of the estimated quantities of conventional oil, gas and condensate that have the potential to be added to reserves by the year 2025 for this province is 6.3 billion barrels of oil equivalent (BBOE) (U. S. Geological Survey World Energy Assessment Team, 2000).

Basin Province

Status of a reconnaissance field study of the Susitna basin, 2011

The Alaska Division of Geological & Geophysical Surveys (DGGS) and Alaska Division of Oil and Gas (DOG), in collaboration with the U.S. Geological Survey (USGS) performed reconnaissance field studies for ten days in late June 2011, in the Susitna basin, directly north of Cook Inlet, south-central Alaska (fig. 1). The purpose of our investigation was to reconnoiter outcrops in the basin and along its periphery to gather new information towards understanding the basin formation history and stratigraphy. This reconnaissance data represents the first step toward better understanding the basin’s hydrocarbon potential, a key component of DGGS’s multi-year In- State Gas Program. This program is focused on collecting baseline geologic information from potential frontier gas basins to encourage new exploration to help, in part, reduce the high cost of energy in rural Alaska. Our work represents the first season of this three-year project. Preliminary results from year two, a companion project within the Nenana and Tanana basins in interior Alaska, are described by Wartes and others (2013). DGGS plans to return to the Susitna basin for follow-up fieldwork during the third and final year of the program. The motivation for developing a better understanding of the Susitna basin stems from the recognition that the Susitna basin shares similar age coal-bearing strata with the adjacent, petroliferous Cook Inlet forearc basin (Barnes, 1966; Reed and Nelson, 1980) and with exhumed strata in the Matanuska Valley forearc basin (Trop and others, 2003) (figs. 1 and 2). Cook Inlet basin has eight producing oil fields, more than 25 producing gas fields, and likely contains many additional undiscovered oil and gas accumulations (LePain and others, in press). Most of the Cook Inlet gas is of microbial origin and apparently was sourced from abundant coalbeds of primarily Miocene age in the Tyonek, Beluga, and Sterling Formations (Claypool and others, 1980; Magoon, 1994). If the biogenic gas model for Cook Inlet is applicable to the Susitna basin, then the latter may be a viable source for Alaska Railbelt and rural energy needs. This brief overview report summarizes the reconnaissance field data collected in the Susitna basin during the first summer of the program. As the data are developed, this report will be followed by interpretive technical reports addressing the stratigraphy, reservoir quality, coal quality and gas potential, hydrocarbon seal integrity, subsurface structure, and uplift history of the basin and sub-basin margins.

Alaska

A probabilistic assessment of tephra-fall hazards at Hanford, Washington, from a future eruption of Mount St. Helens

Hanford, Washington (USA) is the construction site of a multi-billion-dollar high-level nuclear waste treatment facility. This site lies 200 kilometers (km) east of Mount St. Helens (MSH), the most active volcano in the contiguous United States. Tephra from a future MSH eruption could pose a hazard to the air intake and filtration systems at this plant. In this report, we present a probabilistic estimate of the amount of tephra that could fall, and the concentrations of airborne ash that could occur at the Hanford Site during a future eruption. Mount St. Helens has produced four large explosive eruptions in approximately the past 500 years, suggesting that its annual probability of eruption ( P 1 ) is roughly 4/500=0.008. Assuming that a large eruption occurs, we calculate the probability ( P 3|1 ) of a given fall deposit thickness or airborne concentration at Hanford by running about 10,000 simulations of ash-producing eruptions using the atmospheric transport model Ash3d. In each simulation, we calculate the pattern of tephra dispersal, deposit thickness at Hanford, and airborne ash concentration at ground level. As input for each simulation, we choose meteorological conditions from a randomly chosen time in the historical record between 1980 and 2010, using data from the European Centre for Medium-Range Weather Forecasting (ECMWF) Reanalysis (ERA) Interim model. The volume (dense-rock equivalent) of each simulated eruption is randomly chosen from a uniform probability distribution on a log scale from the range of magma volumes (0.008–2.3 cubic kilometers [km 3 ]) estimated for late Holocene eruptions at MSH. Plume heights and durations of each eruption are chosen using empirical correlations between volume, height, and eruption rate, which account for the fact that larger eruptions have higher plumes and last longer. We construct summary tables of final deposit thickness ( T ), maximum ground-level airborne concentration ( C max ), and average ground-level airborne concentration ( C avg ) during tephra-fall for each run. Each table is sorted and ranked by decreasing value of T , C max , or C avg . Conditional probabilities ( P 3|1 ) are derived by dividing rank by n+1, where n is the total number of successful runs. For example, a deposit thickness of 5.10 centimeters (cm) from run 446 is ranked 123 of 9,785 successful runs, yielding P 3|1 =123/9,786=0.01257. Its annual probability is P = P 1 · P 3|1 =0.008×0.01257=0.000101. By interpolation, the deposit thickness ( T 10k ) having an annual probability of 1 in 10,000 ( P = 0.0001) is 5.11 cm. Analogous concentration values are C max,10k =3,819 and C avg,10k =1,513 milligrams per cubic meter (mg/m 3 ), respectively. Independent calculations using the known mass accumulation rate of the deposit (=0.001–0.006 kilograms per square meter per second [kg/m 2 /s]), aggregate fall velocities ( u =0.3–0.8 meters per second [m/s]), and the simple formula , yield similar results, although highly variable fall velocities add significant uncertainty. This formula implies that deposit accumulation rates of millimeters (mm) to greater than 1 cm per hour, which are not uncommon during heavy ash fall, are associated with airborne concentrations of 10 2 –10 3 milligrams per cubic meter (mg/m 3 ). These concentrations are much higher than published measurements (10 -3 –10 1 mg/m 3 ), which record only suspended particles sampled in sheltered areas. During heavy ashfall, most fine ash falls as aggregates. Whether such aggregates will be ingested into air ducts will depend on the aggregate size and fall rate, the fragility of the aggregates, the air duct geometry, intake velocity, and other factors.

Washington