USGS Science⌕ Search

SEARCH · USGS Science

Results for “Applied Energy”

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 469 records · Page 26Linked to original sources

Organic petrology and micro-spectroscopy of Tasmanites microfossils: Applications to kerogen transformations in the early oil window

The transformation of kerogen to hydrocarbons in the early stages of oil generation is critical for understanding the resource potential of liquid-rich shale plays. Organic petrology commonly is used for visual evaluation of type, quality, and thermal maturity of organic matter, but the relationship of visual petrographic changes to chemical transformations is not well characterized. To improve understanding of these processes, organic-walled microfossils of the unicellular green alga Tasmanites (composed of algaenan) in Upper Devonian Ohio Shale (Huron Member, Appalachian Basin) were analyzed by micro-spectroscopy techniques including micro-Fourier transform infrared (micro-FTIR), X-ray photoelectron (XPS), electron probe microanalysis (EPMA), and fluorescence. Immature to mid-oil window maturation sequences of core and outcrop samples with solid bitumen reflectance (BR) and vitrinite reflectance (VR) values ranging from 0.45 to 0.80 %Ro were used. Hydrous pyrolysis was applied to low-maturity (BR: 0.25–0.39 %Ro) Huron and time-correlative New Albany shale samples to create similar artificial maturation sequences for comparison. Micro-FTIR spectroscopy revealed a decrease in the CH 2 /CH 3 ratio with increasing maturity, indicating Tasmanites aliphatic chains become shorter and more branched. Oxygenated functional groups decreased relative to aliphatic stretching bands and increased aromaticity was noted at the highest maturities. In samples that were pyrolyzed for 72 h at temperatures of 300–320 °C (BR: 0.56–0.68 %Ro), Tasmanites showed similar trends, whereas at pyrolysis temperatures of 340 °C and higher (BR > 1.0 %Ro), Tasmanites was pseudomorphed by accumulations of solid bitumen, carbonate and sulfide. Replacement of Tasmanites by these phases in hydrous pyrolysis experiments ≥340 °C and its absence at higher maturities (peak oil, VR and BR ≥ 0.9 %Ro) in naturally matured samples, as documented in a previous study, implies that a large fraction of the algaenan component of original organic carbon is converted to petroleum during thermal maturation. XPS analysis indicated the molar proportion of aliphatic carbon increases with increasing thermal maturity, accompanied by decreases in oxygenated functional groups and olefinic carbon. EPMA of Tasmanites showed highest concentrations of S, with concentrations of redox-sensitive trace elements U, Mo, Ni and V generally at or below detection limits. Decrease in organic S with increasing thermal maturity may be related to cleavage of Tasmanites at C S linkages; however, this relationship was inconsistent and presence of adjacent or entrained nanoscale silicate or sulfide phases may impact measured trace element concentrations. Fluorescence microscopy and spectroscopy showed a red shift in spectral maxima and decreased emission intensities with increasing maturity, interpreted as due to non-radiative energy loss possibly because of increased aromaticity. Collectively, these results provide new insights into the in situ chemical transformations that accompany petrographic changes as oil-prone kerogen converts to petroleum with thermal advance from immature conditions into the mid-oil window.

Organic Geochemistry↗

Utilizing dimensional analysis with observed data to determine the significance of hydrodynamic solutions in coastal hydrology

In this paper, the authors present an analysis of the magnitude of the temporal and spatial acceleration (inertial) terms in the surface-water flow equations and determine the conditions under which these inertial terms have sufficient magnitude to be required in the computations. Data from two South Florida field sites are examined and the relative magnitudes of temporal acceleration, spatial acceleration, and the gravity and friction terms are compared. Parameters are derived by using dimensionless numbers and applied to quantify the significance of the hydrodynamic effects. The time series of the ratio of the inertial and gravity terms from field sites are presented and compared with both a simplified indicator parameter and a more complex parameter called the Hydrodynamic Significance Number (HSN). Two test-case models were developed by using the SWIFT2D hydrodynamic simulator to examine flow behavior with and without the inertial terms and compute the HSN. The first model represented one of the previously-mentioned field sites during gate operations of a structure-managed coastal canal. The second model was a synthetic test case illustrating the drainage of water down a sloped surface from an initial stage while under constant flow. The analyses indicate that the times of substantial hydrodynamic effects are sporadic but significant. The simplified indicator parameter correlates much better with the hydrodynamic effect magnitude for a constant width channel such as Miami Canal than at the non-uniform North River. Higher HSN values indicate flow situations where the inertial terms are large and need to be taken into account.

Florida↗

Comparative performance and trend of remotely sensed phenology and productivity metrics across the Western United States

Vegetation phenology and productivity play a crucial role in surface energy balance, plant and animal distribution, and animal movement and habitat use and can be measured with remote sensing metrics including start of season (SOS), peak instantaneous rate of green-up date (PIRGd), peak of season (POS), end of season (EOS), and integrated vegetation indices. However, for most metrics, we do not yet understand the agreement of remotely sensed data products with near-surface observations. We also need summaries of changes over time, spatial distribution, variability, and consistency in remote sensing dataset metrics for vegetation timing and quality. We compare metrics from 10 leading remote sensing datasets against a network of PhenoCam near-surface cameras throughout the western United States from 2002 to 2014. Most phenology metrics representing a date (SOS, PIRGd, POS, and EOS), rather than a duration (length of spring, length of growing season), better agreed with near-surface metrics but results varied by dataset, metric, and land cover, with absolute value of mean bias ranging from 0.38 (PIRGd) to 37.92 days (EOS). Datasets had higher agreement with PhenoCam metrics in shrublands, grasslands, and deciduous forests than in evergreen forests. Phenology metrics had higher agreement than productivity metrics, aside from a few datasets in deciduous forests. Using two datasets covering the period 1982–2016 that best agreed with PhenoCam metrics, we analyzed changes over time to growing seasons. Both datasets exhibited substantial spatial heterogeneity in the direction of phenology trends. Variability of metrics increased over time in some areas, particularly in the Southwest. Approximately 60% of pixels had consistent trend direction between datasets for SOS, POS, and EOS, with the direction varying by location. In all ecoregions except Mediterranean California, EOS has become later. This study comprehensively compares remote sensing datasets across multiple growing season metrics and discusses considerations for applied users to inform their data choices. des indicators of vegetation timing and quality through metrics such as start of season (SOS), peak instantaneous rate of green-up date (PIRGd), peak of season (POS), end of season (EOS), and integrated vegetation indices. Few comparisons guide users in dataset selection, examine a large spatial extent, and include multiple metrics. This study compares metrics from 10 leading remote sensing datasets against a network of PhenoCam near-surface cameras throughout the Western United States from 2002-2014. Correlation (R2) and mean bias varied substantially by dataset, metric, and land cover. The closest association with PhenoCam measured phenology metrics represented a date (SOS, PIRGd, POS, and EOS) rather than a duration (length of spring, length of growing season), with R2 of individual datasets ranging from 0.03 (SOS) – 0.55 (PIRGd), and absolute value of mean bias ranging from 0.38 (PIRGd) – 37.92 days (EOS). Datasets had higher agreement with PhenoCam metrics in shrublands, grasslands, and deciduous/broadleaf forests than in evergreen forests. Productivity metrics agreed worse than phenology metrics, though some datasets showed high correlations in deciduous/broadleaf forests. Using the two datasets that agreed best with PhenoCam metrics and covered 1982-2016, we conducted a trend analysis to study changes to growing seasons. Trends in phenology exhibited substantial spatial heterogeneity in the direction of trend for both datasets. Variability of metrics increased over time in some areas, particularly in the Southwest. Approximately 60% of pixels had consistent trend direction (both earlier and later) for SOS, POS, and EOS. In all ecoregions except Mediterranean California EOS trended toward a later date. This study provides a comprehensive comparison of remote sensing datasets across many important phenology and productivity metrics and discusses considerations for users to make informed decisions about their data choices.

Arizona, California, Colorado, Idaho, Montana, New↗

Visualizing wading bird optimal foraging decisions with aggregation behaviors using individual-based modeling

Foragers on patchy landscapes must efficiently balance time between searching for and consuming resources to meet their daily energetic requirements. Spatial aggregation foraging behaviors may improve foraging efficiency by sharing information on locations of resource hotspots. Wading birds are an example of patch foragers that form colonial aggregations during the breeding season to obtain sufficient prey energy to sustain themselves and their offspring each day. Here, we describe a spatially-explicit simulation model of wading bird optimal foraging that represents information sharing through visual cues. The overall purpose of the model is to describe how wading bird daily foraging and reproductive success may change with alternative water control management practices that determine spatial availability of prey for wading birds on the landscape, throughout their breeding seasons. Wading birds are simulated as individuals that operate independently, sampling and selecting among patches based on a prey density tolerance threshold, but also use information from other birds to inform their selection decisions. Foraging success is evaluated against the fundamental objectives of (a) fulfilling daily energetic demands and (b) minimizing predation exposure, by tracking individual daily energetic intake and time spent foraging. In this way, the model approximates population level dynamics of wading bird aggregations that emerge through collective decision making of birds simulated at the lower individual level. Key results of this study suggest that aggregation behaviors may improve population-level foraging success rates, and the optimal settling threshold may modulate when resources become more scarce or difficult to find. Thus, the model addresses ecological theory on the advantages of foraging in groups versus independently. This technique is appropriate for evaluating wading bird populations that forage on patchy landscapes, such as seasonally-pulsed wetlands, wherever sufficient information is available to describe (1) foraging behavior (e.g., feeding rate, flight speeds, patch selection decisions), (2) key landscape characteristics, (3) spatial distributions of prey densities among foraging patches, and (4) changes in prey densities through time. The model was designed to predict qualitative, testable spatial patterns of wading bird foraging movements which can be compared with empirical observations and empirically-derived habitat suitability models. These techniques can also be applied to other bird species, such as shorebirds, or more generally to any species that transits between discrete foraging patches.

Ecological Modelling↗

Possible continuous-type (unconventional) gas accumulation in the Lower Silurian "Clinton" sands, Medina Group and Tuscarora Sandstone in the Appalachian Basin; a progress report of the 1995 project activities

INTRODUCTION: In the U.S. Geological Survey's (USGS) 1995 National Assessment of United States oil and gas resources (Gautier and others, 1995), the Appalachian basin was estimated to have, at a mean value, about 61 trillion cubic feet (TCF) of recoverable gas in sandstone and shale reservoirs of Paleozoic age. Approximately one-half of this gas resource is estimated to reside in a regionally extensive, continuous-type gas accumulation whose reservoirs consist of low-permeability sandstone of the Lower Silurian 'Clinton' sands and Medina Group (Gautier and others, 1995; Ryder, 1995). Recognizing the importance of this large regional gas accumulation for future energy considerations, the USGS initiated in January 1995 a multi-year study to evaluate the nature, distribution, and origin of natural gas in the 'Clinton' sands, Medina Group sandstones, and equivalent Tuscarora Sandstone. The project is part of a larger natural gas project, Continuous Gas Accumulations in Sandstones and Carbonates, coordinated in FY1995 by Ben E. Law and Jennie L. Ridgley, USGS, Denver. Approximately 2.6 man years were devoted to the Clinton/Medina project in FY1995. A continuous-type gas accumulation, referred to in the project, is a new term introduced by Schmoker (1995a) to identify those natural gas accumulations whose reservoirs are charged throughout with gas over a large area and whose entrapment does not involve a downdip gas-water contact. Gas in these accumulations is located downdip of the water column and, thus, is the reverse of conventional-type hydrocarbon accumulations. Commonly used industry terms that are more or less synonymous with continuous-type gas accumulations include basin- centered gas accumulation (Rose and others, 1984; Law and Spencer, 1993), tight (low-permeability) gas reservoir (Spencer, 1989; Law and others, 1989; Perry, 1994), and deep basin gas (Masters, 1979, 1984). The realization that undiscovered gas in Lower Silurian sandstone reservoirs of the Appalachian basin probably occurs in a continuous accumulation rather than in conventionally trapped, discrete accumulations represents a significant departure from the 1989 National Assessment (Mast and others, 1989; deWitt, 1993). In 1989, a direct assessment (field-size distributions required for play analysis were unavailable) of the Lower Silurian sandstone play gave, at a mean value, about 1.7 TCF of gas. The 1995 estimate (~30 TCF of gas) is so much greater than the 1989 estimate (~1.7 TCF of gas) because of the interpreted continuous nature of the accumulation and the assessment methodology applied. The methodology for continuous hydrocarbon accumulations assumes that the reservoirs in the accumulation are gas-saturated and takes into account: 1) estimated ultimate recovery (EUR) per well probability distributions, 2) optimum area that a well can drain (spacing), 3) number of untested drill sites having the appropriate spacing area, 4) success ratio of previously drilled holes, and 5) risk (Schmoker, 1995b). Davis (1984), Zagorski (1988, 1991), and Law and Spencer (1993) were among the first petroleum geologists to suggest that gas in the 'Clinton' sands and Medina Group sandstones was trapped in a basin-centered/deep basin accumulation. They recognized many of the earmarks of a basin-centered/deep basin accumulation such as low-permeability reservoirs, abnormally low formation pressure, coalesced gas fields, gas shows or production in most holes drilled, low water yields, and a general lack of structural control on entrapment. Ryder (1995) adopted this interpretation by defining four continuous-type gas plays (6728-6731) in the 'Clinton' sands-Medina Group interval (fig.1). Play 6728 (Clinton/Medina sandstone gas high potential) covers a 17,000 sq mi region of western New York, northwestern Pennsylvania, eastern Ohio, and a small part of westernmost West Virginia that is very favorable for future gas resources (fig.1). Also, this play includes a l

Open-File Report↗

The effects of management practices on grassland birds — An introduction to North American grasslands and the practices used to manage grasslands and grassland birds

The Great Plains of North America is defined as the land mass that encompasses the entire central portion of the North American continent that, at the time of European settlement, was an unbroken expanse of primarily herbaceous vegetation. The Great Plains extend from central Saskatchewan and Alberta to central Mexico and from Indiana to the Rocky Mountains. The expanses of herbaceous vegetation are often referred to as native prairie or native grasslands. Native grasslands share the characteristics of a general uniformity in vegetation structure, dominance by grasses and forbs, a near absence of trees and shrubs, annual precipitation ranging from 25 to 100 centimeters, extreme intra-annual fluctuations in temperature and precipitation, and a flat to rolling topography over which fires can spread. To the west of the Great Plains lie the sagebrush communities of the Great Basin, which extend from British Columbia and Saskatchewan to northern Arizona and New Mexico and from the eastern slopes of the Sierra Nevada and Cascade mountain ranges to western South Dakota. Sagebrush communities share similar characteristics to native grasslands, but their location east of the Rocky Mountains creates a more moderating influence from prevailing westerly winds that affect timing of peak precipitation and growth form of dominant vegetation. Native grasslands and sagebrush communities harbor a diverse array of grassland, wetland, and woodland plant and animal communities that are uniquely adapted to the natural forces of the Great Plains and Great Basin, namely the interactive forces of climate, fire, and grazing. The arrival of European settlers to North America brought profound change to native grassland and sagebrush communities, including the establishment of permanent towns and cities, the proliferation of cropland-based agricultural systems, and the suppression of wildfires. The near extirpation of bison by the 1860s paved the way for dramatic changes in the dominant grazers and a shift in the disturbance patterns that historically influenced vegetation structure. The greatest threat to native grasslands and sagebrush communities in modern times is their loss due to conversion to rowcrop agriculture and to urbanization. Concomitant with habitat loss is a precipitous decline in populations of bird species that evolved with, and are uniquely adapted to, the native grassland and sagebrush habitats. Avian population trends are linked strongly to agricultural land use. Besides outright loss of suitable breeding habitat, agricultural practices affect birds through factors such as pesticide exposure, habitat fragmentation, shifts in predator community composition, and occurrence of brood parasites. Bird populations face other stressors, such as loss of habitat to and behavioral avoidance of urbanized areas, roads, and infrastructure associated with energy production. Despite the many anthropogenic changes to North American grassland and sagebrush communities, some bird species are adaptable and opportunistic in their habitat selection and now utilize one or more human-created habitats. Human-created habitats include pastures, hayfields, agricultural terraces, crop buffer strips, field borders, grassed waterways, fencerows, road rights-of-way, airports, reclaimed coal mines, and planted wildlife cover. Fields of seeded grasslands enrolled in Federal long-term set-aside programs, such as the Conservation Reserve Program in the United States and the Permanent Cover Program in Canada, provide important nesting habitat for grassland bird species. The array of habitats used by birds makes habitat and avian management a complex undertaking, and the scale (for example, local, regional, international) at which management actions can be implemented are such that a universal approach to managing grasslands for the conservation of the entire suite of bird species does not exist. Experienced land managers recognize that it is impossible to manage for all bird species simultaneously, and thus, prioritization is necessary towards those habitats or bird species that the manager or management agency ranks highest for a specific region or management unit. The primary tools available for management are burning, grazing, mowing, herbicide application, and idling, but before choosing a particular practice, a manager will want to consider issues of seasonality, intensity, and frequency. Despite the thousands of studies that are cited in this compendium, much remains unknown about the effects of management practices on bird species. The series of species accounts in this compendium review the current state of knowledge regarding management of grassland and sagebrush bird species and summarize information on the effects of management practices on individual species. The accounts do not give definitive statements on the effects of management practices for any particular species, primarily because there are very few replicated studies in which identical management practices have been applied in the same geographical area with consistent results, which are elements necessary to provide concrete recommendations for the management of a particular species in a particular area. Documentation of the effects of management treatments on individual species through statistically sound methods that incorporate multiple years and locations will further scientists’ and land managers’ knowledge far more than 1–2-year studies that are limited in scope as well as time, but studies of that scope and breadth are rare.

Professional Paper↗

Connectivity of Mojave Desert tortoise populations—Management implications for maintaining a viable recovery network

Executive Summary The historic distribution of Mojave desert tortoises ( Gopherus agassizii ) was relatively continuous across the range, and the importance of tortoise habitat outside of designated tortoise conservation areas (TCAs) to recovery has long been recognized for its contributions to supporting gene flow between TCAs and to minimizing impacts and edge effects within TCAs. However, connectivity of Mojave desert tortoise populations has become a concern because of recent and proposed development of large tracts of desert tortoise habitat that cross, fragment, and surround designated conservation areas. This paper summarizes the underlying concepts and importance of connectivity for Mojave desert tortoise populations by reviewing current information on connectivity and providing information to managers for maintaining or enhancing desert tortoise population connectivity as they consider future proposals for development and management actions. Maintaining an ecological network for the Mojave desert tortoise, with a system of core habitats (TCAs) connected by linkages, is necessary to support demographically viable populations and long-term gene flow within and between TCAs. There are four points for wildlife and land-management agencies to consider when making decisions that could affect connectivity of Mojave desert tortoise populations (for example, in updating actions in resource management plans or amendments that could help maintain or restore functional connectivity in light of the latest information): Management of all desert tortoise habitat for persistence and connectivity . Desert tortoise populations continue to decline within most TCAs, and it is unlikely that trends are better in populations outside protected areas. Fragmentation exacerbates negative population trends by breaking large continuous populations into smaller isolated populations. Connectivity within large populations can enhance resilience to localized disturbances due to rescue by neighboring individuals. In contrast, smaller fragmented populations are resistant to rescue by their isolation and thus could suffer irreversible declines to extirpation from a variety of threats and stochastic events. Enhanced threat reduction to reverse declines within TCAs and to maintain occupied habitat in the surrounding matrix would help reduce the variability in population growth rates and improve the resilience of protected populations even while implementing efforts to improve connectivity. Each TCA has unique strengths and weaknesses regarding its ability to support minimum sustainable populations based on areal extent and its ability to support population increases based on landscape connection with adjacent populations. Considering how proposed projects (inside or outside of TCAs) affect connectivity and the ability of TCAs to support at least 5,000 adult tortoises (the numerical goal for each TCA) could help managers to maintain the resilience of TCAs to population declines. The same project, in an alternative location, could have very different impacts on local and regional populations. For example, within the habitat matrix surrounding TCAs, narrowly delineated corridors may not allow for natural population dynamics if they do not accommodate overlapping home ranges along most of their widths so that tortoises reside, grow, find mates, and produce offspring that can replace older tortoises. In addition, most habitat outside TCAs may receive more surface disturbance than habitat within TCAs. Therefore, managing the entire remaining matrix of desert tortoise habitat for permeability may be better than delineating fixed corridors. These concepts apply, especially given uncertainty about long-term condition of habitat, within and outside of TCAs under a changing climate. Ultimately, questions such as “ What are the critical linkages that need to be protected ?” could be better framed as “ How can we manage the remaining habitat matrix in ways that sustain ecological processes and habitat suitability for special status species ?” Land-management decisions made in the context of the latter question may be more conducive to maintenance of a functional ecological network. Limitations on landscape-level disturbance across habitat managed for the desert tortoise Clearly delineating habitat linkages and differentiating them from non-delineated areas by the uses that are permitted or prohibited within them by specific management guidelines can help achieve functional connectivity. Such guidelines would be most effective if they considered and accounted for all surface disturbances (for example, temporary disturbances such as fiberoptic lines or off-highway vehicle routes, right-of-ways, utility-scale solar development, urbanization) to the extent possible. A weighted framework that varies with the permanence or severity of the disturbance, and can be additive to quantify cumulative effects, could be useful (Xiong, 2020). For example, minor roads can alter tortoise movements independently of other features (Peaden and others, 2017; Hromada and others, 2020), but if the isolated dirt road is accompanied by a powerline that encourages raven predation (Xiong, 2020), then the two features together may be additive. Ignoring minor or temporary disturbance on the landscape could result in a cumulatively large impact that is not explicitly acknowledged (Goble, 2009); therefore, understanding and quantifying all surface disturbance on a given landscape is prudent. In California, the Bureau of Land Management established 0.1–1.0 percent caps on new surface-disturbance for TCAs and mapped linkages that address the issues described in number 1 of this list. Nevada, Utah, and Arizona currently do not have surface-disturbance limits. Limits comparable to those in the Desert Renewable Energy Conservation Plan (DRECP) would be 0.5 percent within TCAs and 1 percent within the linkages modeled by Averill-Murray and others (2013). Limits in some areas of California within the Desert Renewable Energy Conservation Plan, such as Ivanpah Valley, are more restrictive, at 0.1 percent. Continuity across the state line in Nevada could be achieved with comparable limits in the adjacent portion of Ivanpah Valley, as well as the Greater Trout Canyon Translocation Area and the Stump Springs Regional Augmentation Site. These more restrictive limits would help protect remaining habitat in the major interstate connectivity pathway through Ivanpah Valley and focal areas of population augmentation that provide additional population connectivity along the western flank of the Spring Mountains. In a recent study that analyzed 13 years of desert tortoise monitoring data, nearly all desert tortoise observations were at sites in which 5 percent or less of the surrounding landscape within 1 kilometer was disturbed (Carter and others, 2020a). To help maintain tortoise habitability and permeability across all other non-conservation-designated tortoise habitat, all surface disturbance could be limited to less than 5-percent development per square kilometer because the 5-percent threshold for development is the point at which tortoise occupation drops precipitously (Carter and others, 2020a). However, although individual desert tortoises were observed at development levels up to 5 percent, we do not know the fitness or reproductive characteristics of these individuals. This level of development also may not allow for long-term persistence of healthy populations that are of adequate size needed for demographic or functional connectivity; therefore, a conservative interpretation suggests that, ideally, development could be lower. Lower development levels would be particularly useful in areas within the upper 5th percentile of connectivity values modeled by Gray and others (2019). Reducing ancillary threats in places where connectivity is restricted to narrow strips of habitat, for example, narrow mountain passes or vegetated strips between solar development, could enhance the functionality of these vulnerable linkages. In such areas, maintaining multiple, redundant linkages could further enhance overall connectivity. Minimization of mortality from roads and maximization of passage under roads . Roads pose a significant threat to the long-term persistence of local tortoise populations, and roads of high traffic volume lead to severe population declines, which ultimately fragments populations farther away from the roads. Three points (a.–c.) pertain to reducing direct mortality of tortoises on the many paved roads that cross desert tortoise habitat and to maintaining a minimal level of permeability across these roads: Tortoise-exclusion fencing tied into culverts, underpasses, overpasses, or other passages below roads in desert tortoise habitat, would limit vehicular mortality of tortoises and provide opportunities for movement across the roads. Installation of shade structures on the habitat side of fences installed in areas with narrow population-depletion zones would limit overheating of tortoises that may pace the fence. Passages below highways could be maintained or retrofitted to ensure safe tortoise access, for example, by filling eroded drop-offs or modifying erosion-control features such as rip-rap to make them safer and more passable for tortoises. Wildlife management agencies could work with transportation departments to develop construction standards that are consistent with hydrologic/erosion management goals, while also incorporating a design and materials consistent with tortoise survival and passage and make the standards widely available. The process would be most effective if the status of passages was regularly monitored and built into management plans. Healthy tortoise populations along fenced highways could be supported by ensuring that land inside tortoise-exclusion fences is not so degraded that it leads to degradation of tortoise habitat outside the exclusion areas. For example, severe invasive plant infestations inside a highway exclusion could cause an increase of invasive plants outside the exclusion area and degrade habitat; therefore, invasive plants inside road rights of way could be mown or treated with herbicide to limit their spread into adjacent tortoise habitat and minimize the risk of these plants carrying wildfires into adjacent habitat. Adaptation of management based on new information . Future research will continue to build upon and refine models related to desert tortoise population connectivity and develop new ones. New models could consider landscape levels of development and be constructed such that they share common foundations to support future synthesis efforts. If model development was undertaken in partnership with entities that are responsible for management of desert tortoise habitat, it would facilitate incorporation of current and future modeling results into their land management decisions. There are specific topics that may be clarified with further evaluation: The effects of climate change on desert tortoise habitat, distribution, and population connectivity; The effects of large-scale fires, especially within repeatedly burned habitat, on desert tortoise distribution and population connectivity; The ability of solar energy facilities or similar developments to support tortoise movement and presence by leaving washes intact; leaving native vegetation intact whenever possible, or if not possible, mowing the site, allowing vegetation to re-sprout, and managing weeds; and allowing tortoises to occupy the sites; and The design and frequency of underpasses necessary to maintain functional demographic and genetic connectivity across linear features, like highways.

Arizona, California, Nevada↗

Understanding the Avian-Impact Offset Method—A tutorial

Biodiversity offsetting, or compensatory mitigation, is increasingly being used in temperate grassland and wetland ecosystems to compensate for unavoidable environmental damage from anthropogenic disturbances such as energy development and road construction. Energy-extraction and -generation facilities continue to proliferate across the natural landscapes of the United States, yet mitigation tools to ameliorate the negative behavioral effects on wildlife from these types of facilities are rarely implemented. Scientists from the U.S. Geological Survey conducted a 10-year before-after-control-impact (commonly referred to as BACI) study that evaluated the displacement effects of wind facilities on breeding grassland birds. The study determined behavioral avoidance for 7 of 9 species. This research is notable because of its design, geographical scope, and duration, which allowed for the determination of immediate, short-term effects; delayed or sustained effects; and discrete distances at which effects occurred. In addition, the U.S. Fish and Wildlife Service and Ducks Unlimited conducted a 3-year concurrent-year paired-reference study to determine behavioral avoidance for five species of dabbling ducks. By quantifying displacement rate from these two studies, U.S. Geological Survey and U.S. Fish and Wildlife Service scientists developed the Avian-Impact Offset Method (AIOM) to quantify and compensate for loss in value of breeding habitat. The AIOM converts the biological value (that is, number of bird pairs) lost by way of avoidance and estimates the site-specific number of hectares of grasslands and number of wetlands needed to compensate for displaced pairs of grassland birds and waterfowl. By converting biological value to traditional units of measure in which land is described and purchased or sold, the AIOM lends itself readily to the delivery of offsetting measures such as easement protections and restoration projects. The AIOM tool is applicable to wind, solar, oil, gas, and transportation infrastructure. This tutorial was designed to increase awareness of the AIOM and to promote its proper application. The tutorial is divided into four sections, each of which explains a discrete topic concerning aspects of behavioral displacement. The first section provides geographical and biological context, and the second section describes the field and statistical methods and results. The third section provides step-by-step instructions for applying the AIOM to several scenarios involving grassland birds or waterfowl at wind or oil facilities. The fourth section describes decision-support tools created to implement the AIOM. The appendices provide the actual field protocols constituting the methods for the research, provide detailed results by species and wind facility for that research, and provide detailed instructions for downloading and applying the decision-support tools.

Open-File Report↗

Application of Raman spectroscopy as thermal maturity probe in shale petroleum systems: Insights from natural and artificial maturation series

Raman spectroscopy was studied as a thermal maturity probe in a series of Upper Devonian Ohio Shale samples from the Appalachian Basin spanning from immature to dry gas conditions. Raman spectroscopy also was applied to samples spanning a similar thermal range created from 72-h hydrous pyrolysis (HP) experiments of the Ohio Shale at temperatures from 300 to 360 °C and isothermal HP experiments lasting up to 100 days of similar Devonian–Mississippian New Albany Shale. Raman spectra were treated by automated evaluation software based on iterative and simultaneous modeling of signal and baseline functions to decrease subjectivity. Spectra show robust correlation to measured solid bitumen reflectance (BR o ) values and were therefore used to construct logarithmic regression relationships for calculation of BR o equivalent values. Raman spectra show considerable differences between natural samples and HP residues with similar measured BR o values, indicating as-yet undetermined differences in carbon chemistry. We speculate this result may be due to differences in the sampling interactions of Raman vs reflectance measurements, and the incomplete nature of maturation reactions in the time-limited hydrous pyrolysis residues. Samples used in this study are similar in organic assemblage (dominantly solid bitumen) to other commonly exploited North American shale petroleum systems, i.e., Bakken, Barnett, Duvernay, Fayetteville, and Woodford shales. Therefore, results presented herein may be broadly applicable to other important shale plays. However, caution is suggested and Raman spectroscopy as a thermal probe may need individual calibration in each shale play due to differences in solid bitumen carbon chemistry.

Energy & Fuels↗

3-D geological modeling for numerical flow simulation studies of gas hydrate reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. This study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-step gridding were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a 3-D trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

Alaska↗

Evapotranspiration over spatially extensive plant communities in the Big Cypress National Preserve, southern Florida, 2007-2010

Evapotranspiration (ET) was quantified over plant communities within the Big Cypress National Preserve (BCNP) using the eddy covariance method for a period of 3 years from October 2007 to September 2010. Plant communities selected for study included Pine Upland, Wet Prairie, Marsh, Cypress Swamp, and Dwarf Cypress. These plant communities are spatially extensive in southern Florida, and thus, the ET measurements described herein can be applied to other humid subtropical locations such as the Everglades. The 3-year mean annual ET was about 1,000, 1,050, 1,100, 930, and 900 mm (millimeters) at the Dwarf Cypress, Wet Prairie, Cypress Swamp, Pine Upland, and Marsh sites, respectively. Spatial differences in annual ET were considerable due to the recovery of the Marsh site from extensive forest fire and drought conditions. Temporal variability in annual ET was relatively small at sites that were well watered (Dwarf Cypress, Wet Prairie, Cypress Swamp, Pine Upland) over the 3-year study. In other words, locations that were well watered appeared to have similar annual ET rates. Diurnal water-level variability was observed in response to ET and was less at flooded sites than at dry sites. For example, surface-water levels declined about 1.5 mm in response to ET at the flooded Cypress Swamp site during July 18-22, 2008 and declined about 10 mm in response to ET at the flooded Dwarf Cypress site from April 18-27, 2008. Specific yield was computed using ET estimates and diurnal water-level variability measured at the dry Pine Upland site as a check on the accuracy of the eddy covariance method. Water levels repeatedly dropped about 15 mm on average in response to ET at the Pine Upland site from April 27 to May 4, 2008. ET was about 3 mm on each of these days, resulting in a reasonable estimate for specific yield of 0.2 for the sandy soils at the Pine Upland site. Monthly ET estimates exhibited seasonal variation. ET was generally greatest between March to October when solar radiation was relatively large, and least from November to February when solar radiation was small. Monthly ET was greatest in the spring and summer at the Cypress Swamp site, reaching rates as large as 140 mm. The large ET rates at this site coincide with the most active period of cypress growth during late spring and early summer. Cypress trees begin to senesce in late summer reducing transpiration. Net radiation and available energy explained most of the variability in ET observed at all five sites. Mean annual and monthly net radiation varied among the sites in response to cloud cover and the albedo of the land surface and plant community. Net radiation was greatest at the Cypress Swamp site, averaging about 130 W/m 2 (watts per square meter) during the 3-year study. Net radiation was generally less at the Dwarf Cypress site, averaging about 115 W/m 2 over 3 years. The Dwarf Cypress site apparently has the largest albedo, which likely is due to the sparse canopy and a highly reflective, calcareous, periphyton-covered land surface. Furthermore, mean annual net radiation was least in the first year of the study, which likely was due to greater cloud cover during a relatively wet year. In contrast, net radiation was greatest in the second year of the study, which likely was due to less cloud cover during a relatively dry year.

Scientific Investigations Report↗

Remote sensing-based actual evapotranspiration assessment in a data-scarce area of Brazil: A case study of the Urucuia Aquifer System

The large groundwater reserves of the Urucuia Aquifer System (UAS) enabled agricultural development and economic growth in the western Bahia State, in northeastern Brazil. Over the last several years, concern has grown around the aquifer’s diminishing water levels, and water balance (WB) studies are in demand. Considering the lack of measured actual evapotranspiration (ET a ), a major component of the water cycle, this work uses the Operational Simplified Surface Energy Balance (SSEBop) model to estimate ET a , and compares it to basin-scale estimates from the Soil Moisture Accounting Procedure (SMAP) monthly model and from an annual WB closure method, based on gridded meteorological data and the Gravity Recovery and Climate Experiment (GRACE) product. Additionally, a comparative assessment of different versions of the SSEBop parameterization was performed. Moderate Resolution Imaging Spectroradiometer (MODIS) imagery was used to implement eight different versions of the SSEBop algorithm over the UAS between 2000 and 2013. SSEBop and SMAP ET a yielded similar seasonal patterns, with correlation coefficient (r) up to 0.65, mean difference (MD) of 0.8 mm/month and mean absolute difference (MAD) of 18.5 mm/month. Comparison of SSEBop annual ET a estimates to annual SMAP and WB closure estimates yielded low MD (12.1 and −7.3 mm/year, respectively) and MAD (82.5 and 82.8 mm/year, respectively), but also low r values (0.00 and 0.37, respectively). The comparison of the different SSEBop versions indicated the need to incorporate a calibration step of the aerodynamic heat resistance (r ah ) parameter. SSEBop results were also used for land cover and drought monitoring. Analysis indicates that agriculture, associated with an increasing trend of atmospheric evaporative demand, is responsible for the decrease in groundwater levels and streamflow in the studied time period.

Urucuia Aquifer System↗

Spatial variability in biodegradation rates as evidenced by methane production from an aquifer

Accurate predictions of carbon and energy cycling rates in the environment depend on sampling frequencies and on the spatial variability associated with biological activities. We examined the variability associated with anaerobic biodegradation rates at two sites in an alluvial sand aquifer polluted by municipal landfill leachate. In situ rates of methane production were measured for almost a year, using anaerobic wells installed at two sites. Methane production ranged from 0 to 560 μmol · m -2 · day -1 at one site (A), while a range of 0 to 120,000 μmol · m -2 · day -1 was measured at site B. The mean and standard deviations associated with methane production at site A were 17 and 57 μmol · m -2 · day -1 , respectively. The comparable summary statistics for site B were 2,000 and 9,900 μmol · m -2 · day -1 . The coefficients of variation at sites A and B were 340 and 490%, respectively. Despite these differences, the two sites had similar seasonal trends, with the maximal rate of methane production occurring in summer. However, the relative variability associated with the seasonal rates changed very little. Our results suggest that (i) two spatially distinct sites exist in the aquifer, (ii) methanogenesis is a highly variable process, (iii) the coefficient of variation varied little with the rate of methane production, and (iv) in situ anaerobic biodegradation rates are lognormally distributed.

Applied and Environmental Microbiology↗

Applied geochemistry, geology and mineralogy of the northernmost Carlin trend, Nevada

Investigations in the northernmost Carlin trend were undertaken to advance understanding of the geochemical signatures and genesis of precious metal deposits in the trend . Two fundamental geologic relationships near the trend significantly affect regional geochemical distributions: a remarkably intact lower Paleozoic stratigraphic sequence of siliceous rocks in the upper plate of the middle Paleozoic Roberts Mountains thrust, and the widespread repetition of rocks high in the upper plate during late Paleozoic thrusting that thickens the cover above mineralized rock in the lower plate. A compilation of previously published chemical analyses of 440 stream sediment samples and 115 rocks from two 7 1/2-minute quadrangles, as well as new chemical analyses of approximately 1,000 drill core samples in a 1,514 m (4,970 ft) hole through the Rodeo Creek deposit were used to construct three-dimensional element distribution models that highlight metal zonation in the mineralized systems. The Rodeo Creek deposit comprises deep Ag base-metal ± Au-mineralized rock below the Roberts Mountains thrust and contains an unusually high Ag/Au ratio greater than 30. Stacked geochemical halos related to the deposit are confined to the lower plate of the Roberts Mountains thrust and include two horizons of Hg, Cu, and Zn anomalies-as much as 180 m above the deposit-that mostly result from mercurian sphalerite. Extremely subtle indications of mineralization in the upper plate of the Roberts Mountains thrust above the deposit include arsenopyrite overgrowths on small pyrite crystals in 50- to 75-μm-wide clay-carbonate veinlets that lack alteration halos, arsenical rims on small disseminated crystal of recrystallized diagenetic pyrite, and partial replacement of diagenetic pyrite by tennantite. Some of these minerals contain anomalously high Au. However, these As-(Au)-bearing rocks most likely represent another locus of largely untested mineralized rock rather than distal halos related to either the Rodeo Creek or the nearby Dee and Storm gold deposits. Application of micromineralogic techniques helped to identify mineral assemblages that are specific to mineralization and provided an empirical foundation for interpretations of geochemical halos in the Carlin trend . District-scale geochemical patterns of several elements in stream sediments and surface rocks coincide with the northernmost Carlin trend and can be used to explore for Carlin -type deposits. Concentrations of elevated As and Sb in stream sediments (as much as 54 ppm As) have northwest-elongate lobate patterns that clearly outline the trend across a width of approximately 4 km. Arsenic contents of exposed rocks (as much as 90 ppm As) strongly correlate with As contents of derivative stream sediments, and rock contents of Sb show a somewhat lesser but nonetheless strong and similar correspondence. Factor analysis of stream-sediment data shows that those factor scores that are correlated with As, Sb, Au, and Pb also are high along the trend and suggest that mineralized rocks may be present. Although As was not detected by scanning electron microscope-energy dispersive spectrometer (SEM-EDS) studies in heavy mineral concentrates of high-As stream sediments in the Carlin trend , X-ray absorption near-edge spectra (XANES) of selected light fractions of stream sediment samples indicate that Al-bearing phases, such as gibbsite, amorphous Al oxyhydroxides, or aluminosilicate clay minerals host most of the As(V). The best fit, visually and in terms of the lowest residual, was obtained by a model compound of As(V) sorbed to gibbsite. Thus, most As in stream sediments derived from altered rock within the Carlin trend apparently is contained in light fractions. The geochemical character of young, unconsolidated, postmineral deposits that cover mineralized rocks on the Carlin trend partly results from mineralized sources along the trend . Concentration of As in the Miocene Carlin Formation shows an exceptionally well developed progressive increase to about 30 ppm As as altered rock surrounding the trend is approached. Mineralized and/or altered rock fragments probably have been shed directly into the sedimentary basin of the Carlin Formation, and migration of As, now fixed as As(V), also may have occurred in the supergene environment after material was recycled out of the Carlin Formation and into present-day gulleys.

Nevada↗

The aerosphere as a network connector of organisms and their diseases

Aeroecological processes, especially powered flight of animals, can rapidly connect biological communities across the globe. This can have profound consequences for evolutionary diversification, energy and nutrient transfers, and the spread of infectious diseases. The latter is of particular consequence for human populations, since migratory birds are known to host diseases which have a history of transmission into domestic poultry or even jumping to human hosts. In this chapter, we present a scenario under which a highly pathogenic avian influenza (HPAI) strain enters North America from East Asia via post-molting waterfowl migration. We use an agent-based model (ABM) to simulate the movement and disease transmission among 10 6 generalized waterfowl agents originating from ten molting locations in eastern Siberia, with the HPAI seeded in only ~10 2 agents at one of these locations. Our ABM tracked the disease dynamics across a very large grid of sites as well as individual agents, allowing us to examine the spatiotemporal patterns of change in virulence of the HPAI infection as well as waterfowl host susceptibility to the disease. We concurrently simulated a 12-station disease monitoring network in the northwest USA and Canada in order to assess the potential efficacy of these sites to detect and confirm the arrival of HPAI. Our findings indicated that HPAI spread was initially facilitated but eventually subdued by the migration of host agents. Yet, during the 90-day simulation, selective pressures appeared to have distilled the HPAI strain to its most virulent form (i.e., through natural selection), which was counterbalanced by the host susceptibility being conversely reduced (i.e., through genetic predisposition and acquired immunity). The monitoring network demonstrated wide variation in the utility of sites; some were clearly better at providing early warnings of HPAI arrival, while sites further from the disease origin exposed the selective dynamics which slowed the spread of the disease albeit with the result of passing highly virulent strains into southern wintering locales (where human impacts are more likely). Though the ABM presented had generalized waterfowl migration and HPAI disease dynamics, this exercise demonstrates the power of such simulations to examine the extremely large and complex processes which comprise aeroecology. We offer insights into how such models could be further parameterized to represent HPAI transmission risks as well as how ABMs could be applied to other aeroecological questions pertaining to individual-based connectivity.

Book chapter↗

Resonance of a fluid-driven crack: Radiation properties and implications for the source of long-period events and harmonic tremor

A dynamic source model is presented, in which a three-dimensional crack containing a viscous compressible fluid is excited into resonance by an impulsive pressure transient applied over a small area ΔS of the crack surface. The crack excitation depends critically on two dimensionless parameters called the crack stiffness, C = ( b /μ)( L / d ), and viscous damping loss, F = (12η L )/(ρ ƒ d 2 α), where b is the bulk modulus, η is the viscosity, ρ ƒ is the density of the fluid, μ is the rigidity, α is the compressional velocity of the solid, L is the crack length, and d is the crack thickness. The first parameter characterizes the ability of the crack to vibrate and shapes the spectral signature of the source, and the second quantifies the effect of fluid viscosity on the duration of resonance. Resonance is sustained by a very slow wave trapped in the fluid-filled crack. This guided wave, called the crack wave, is similar to the tube wave propagating in a fluid-filled borehole; it is inversely dispersive, showing a phase velocity that decreases with increasing wavelength, and its wave speed is always lower than the acoustic velocity of the fluid, decreasing rapidly as the crack stiffness increases. The source spectrum shows many sharp peaks characterizing the individual modes of vibration of the crack; the variation of spectral shape, both in the number and width of peaks, is surprisingly complex, reflecting the interference between the lateral and longitudinal modes of resonance, as well as nodes for these modes. The far-field spectrum is marked by narrow-band dominant and subdominant peaks that reflect the interaction of the various source modes. The frequency of the dominant spectral peak radiated by the source is independent of the radiation direction. The frequency, bandwidth, and spacing of the resonant peaks are strongly dependent on the crack stiffness, larger values of the stiffness factor shifting these peaks to lower frequencies and decreasing their bandwidth. The excitation of a particular mode depends on the position of the trigger and on the extent of the crack surface affected by the pressure transient. Fluid viscosity decreases the amplitudes of the main spectral peaks, smears out the finer structure of the spectrum, and greatly reduces the duration of the radiated signal. The energy loss by radiation is stronger for high frequencies, producing a seismic signature that is marked by a high-frequency content near the onset of the signal and dominated by a longer-period component of much longer duration in the signal coda. Such signature is in harmony with those displayed by long-period events observed on active volcanoes and in hydrofracture experiments. The very low velocity which is possible in a crack with high stiffness ( C ≥ 100) also provides an attractive explanation for very long period tremor, such as type 2 tremor at Aso volcano, Japan, without the requirement of an unrealistically large magma container. The standing wave pattern set up on the crack surface by the sustained resonance in the fluid is observable in the near field of the crack, suggesting that the location and extent of the source may be estimated from the mapping of the pattern of nodes and antinodes seen in its vicinity. According to the model, the long-period event and harmonic tremor share the same source but differ in the boundary conditions for fluid flow and in the triggering mechanism setting up the resonance of the source, the former being viewed as the impulse response of the tremor generating system and the latter representing the excitation due to more complex forcing functions.

Journal of Geophysical Research Solid Earth↗

Stimulation of methane generation from nonproductive coal by addition of nutrients or a microbial consortium

Biogenic formation of methane from coal is of great interest as an underexploited source of clean energy. The goal of some coal bed producers is to extend coal bed methane productivity and to utilize hydrocarbon wastes such as coal slurry to generate new methane. However, the process and factors controlling the process, and thus ways to stimulate it, are poorly understood. Subbituminous coal from a nonproductive well in south Texas was stimulated to produce methane in microcosms when the native population was supplemented with nutrients (biostimulation) or when nutrients and a consortium of bacteria and methanogens enriched from wetland sediment were added (bioaugmentation). The native population enriched by nutrient addition included Pseudomonas spp., Veillonellaceae , and Methanosarcina barkeri . The bioaugmented microcosm generated methane more rapidly and to a higher concentration than the biostimulated microcosm. Dissolved organics, including long-chain fatty acids, single-ring aromatics, and long-chain alkanes accumulated in the first 39 days of the bioaugmented microcosm and were then degraded, accompanied by generation of methane. The bioaugmented microcosm was dominated by Geobacter sp., and most of the methane generation was associated with growth of Methanosaeta concilii . The ability of the bioaugmentation culture to produce methane from coal intermediates was confirmed in incubations of culture with representative organic compounds. This study indicates that methane production could be stimulated at the nonproductive field site and that low microbial biomass may be limiting in situ methane generation. In addition, the microcosm study suggests that the pathway for generating methane from coal involves complex microbial partnerships.

Texas↗

Aquatic pollution increases use of terrestrial prey subsidies by stream fish

Stream food webs are connected with their riparian zones through cross-ecosystem movements of energy and nutrients. The use and impact of terrestrial subsidies on aquatic consumers is determined in part by in situ biomass of aquatic prey. Thus, stressors such as aquatic pollutants that greatly reduce aquatic secondary production could increase the need for and reliance of stream consumers on terrestrial resource subsidies. To test this hypothesis, we surveyed stream fish, their diets and resource availability in 16 subalpine streams over a regional gradient of trace metals known to strongly impact aquatic insect communities (i.e. fish prey) in the Colorado Rocky Mountains, USA. Fish increased their reliance on terrestrial insect prey as stream metals increased. Relative biomass of terrestrial insects in stomach contents of Brook and Brown Trout increased with respect to aquatic insect biomass and total stomach contents. Drifting insect biomass showed a declining trend for aquatic, but not terrestrial insects, over the metal gradient. Trout densities were unrelated to metal concentrations in streams where we found fish. Synthesis and applications . Our results indicate that diets of aquatic consumers can become more terrestrial as aquatic stressors that limit in situ food production increase and that these subsidies may compensate for loss of aquatic resources. This work implies an important connection between preserving aquatic–terrestrial linkages and management of fish populations in stressed watersheds. Specifically, intact riparian zones and aquatic–terrestrial linkages are likely to be important for maintaining trout production in streams with moderate metal contamination.

Colorado↗