USGS Science⌕ Search

SEARCH · USGS Science

Results for “Applied Economics”

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.

377 records · Page 21Linked to original sources

Depositional setting, petrology and chemistry of Permian coals from the Paraná Basin: 2. South Santa Catarina Coalfield, Brazil

In Brazil economically important coal deposits occur in the southern part of the Paraná Basin, where coal seams occur in the Permian Rio Bonito Formation, with major coal development in the states of Rio Grande de Sul and Santa Catarina. The current paper presents results on sequence stratigraphic interpretation of the coal-bearing strata , and petrological and geochemical coal seam characterization from the South Santa Catarina Coalfield, Paraná Basin. In terms of sequence stratigraphic interpretation the precursor mires of the Santa Catarina coal seams formed in an estuarine-barrier shoreface depositional environment , with major peat accumulation in a high stand systems tract (Pre-Bonito and Bonito seams), a lowstand systems tract (Ponta Alta seam, seam A, seam B) and a transgressive systems tract (Irapuá, Barro Branco and Treviso seams). Seam thicknesses range from 1.70 to 2.39 m, but high proportions of impure coal (coaly shale and shaley coal), carbonaceous shale and partings reduce the net coal thickness significantly. Coal lithoypes are variable, with banded coal predominant in the Barro Branco seam, and banded dull and dull coal predominantly in Bonito and Irapuá seams, respectively. Results from petrographic analyses indicate a vitrinite reflectance range from 0.76 to 1.63 %Rrandom (HVB A to LVB coal). Maceral group distribution varies significantly, with the Barro Branco seam having the highest vitrinite content (mean 67.5 vol%), whereas the Irapuá seam has the highest inertinite content (33.8 vol%). Liptinite mean values range from 7.8 vol% (Barro Branco seam) to 22.5 vol% (Irapuá seam). Results from proximate analyses indicate for the three seams high ash yields (50.2 – 64.2 wt.%). Considering the International Classification of in-Seam Coals, all samples are in fact classified as carbonaceous rocks (> 50 wt.% ash). Sulfur contents range from 3.4 to 7.7 wt.%, of which the major part occurs as pyritic sulfur. Results of X-ray diffraction indicate the predominance of quartz and kaolinite (also pyrite). Gypsum , gibbsite , jarosite and calcite were also identified in some samples. Feldspar was noted but is rare. The major element distribution in the three seams (coal basis) is dominated by SiO 2 (31.3 wt.%, mean value), Al 2 O 3 (14.5 wt.%, mean value) and Fe 2 O 3 (6.9 wt.%, mean value). Considering the concentrations of trace elements that are of potential environmental hazards the Barro Branco, Bonito and Irapuá seams (coal base) are significantly enriched in Co (15.7 ppm), Cr (54.5 ppm), Li (59.3 ppm), Mn (150.4 ppm), Pb (58.0 ppm) and V (99.6 ppm), when compared to average trace elements contents reported for U. S. coals. Hierarchical cluster analysis identified, based on similarity levels, three groups of major elements and seven groups of trace elements. Applying discriminant analyses using trace and major element distribution, it could be demonstrated that the three seams from Santa Catarina show distinct populations in the discriminant analyses plots, and also differ from the coals of Rio Grande do Sul analyzed in a previous study.

Santa Catarina↗

A practical primer on geostatistics

Introduction The Challenge —Most geological phenomena are extraordinarily complex in their interrelationships and vast in their geographical extension. Ordinarily, engineers and geoscientists are faced with corporate or scientific requirements to properly prepare geological models with measurements involving a small fraction of the entire area or volume of interest. Exact description of a system such as an oil reservoir is neither feasible nor economically possible. The results are necessarily uncertain. Note that the uncertainty is not an intrinsic property of the systems; it is the result of incomplete knowledge by the observer. The Aim of Geostatistics —The main objective of geostatistics is the characterization of spatial systems that are incompletely known, systems that are common in geology. A key difference from classical statistics is that geostatistics uses the sampling location of every measurement. Unless the measurements show spatial correlation, the application of geostatistics is pointless. Ordinarily the need for additional knowledge goes beyond a few points, which explains the display of results graphically as fishnet plots, block diagrams, and maps. Geostatistical Methods —Geostatistics is a collection of numerical techniques for the characterization of spatial attributes using primarily two tools: probabilistic models, which are used for spatial data in a manner similar to the way in which time-series analysis characterizes temporal data, or pattern recognition techniques. The probabilistic models are used as a way to handle uncertainty in results away from sampling locations, making a radical departure from alternative approaches like inverse distance estimation methods. Differences with Time Series —On dealing with time-series analysis, users frequently concentrate their attention on extrapolations for making forecasts. Although users of geostatistics may be interested in extrapolation, the methods work at their best interpolating. This simple difference has significant methodological implications. Historical Remarks —As a discipline, geostatistics was firmly established in the 1960s by the French engineer Georges Matheron, who was interested in the appraisal of ore reserves in mining. Geostatistics did not develop overnight. Like other disciplines, it has built on previous results, many of which were formulated with different objectives in various fields. Pioneers —Seminal ideas conceptually related to what today we call geostatistics or spatial statistics are found in the work of several pioneers, including: 1940s: A.N. Kolmogorov in turbulent flow and N. Wiener in stochastic processing; 1950s: D. Krige in mining; 1960s: B. Mathern in forestry and L.S. Gandin in meteorology Calculations —Serious applications of geostatistics require the use of digital computers. Although for most geostatistical techniques rudimentary implementation from scratch is fairly straightforward, coding programs from scratch is recommended only as part of a practice that may help users to gain a better grasp of the formulations. Software —For professional work, the reader should employ software packages that have been thoroughly tested to handle any sampling scheme, that run as efficiently as possible, and that offer graphic capabilities for the analysis and display of results. This primer employs primarily the package Stanford Geomodeling Software (SGeMS) - recently developed at the Energy Resources Engineering Department at Stanford University - as a way to show how to obtain results practically. This applied side of the primer should not be interpreted as the notes being a manual for the use of SGeMS. The main objective of the primer is to help the reader gain an understanding of the fundamental concepts and tools in geostatistics. Organization of the Primer —The chapters of greatest importance are those covering kriging and simulation. All other materials are peripheral and are included for better comprehension of these main geostatistical modeling tools. The choice of kriging versus simulation is often a big puzzle to the uninitiated, let alone the different variants of both of them. Chapters 14, 18, and 19 are intended to shed light on those subjects. The critical aspect of assessing and modeling spatial correlation is covered in chapter 7. Chapters 2 and 3 review relevant concepts in classical statistics. Course Objectives —This course offers stochastic solutions to common problems in the characterization of complex geological systems. At the end of the course, participants should have: an understanding of the theoretical foundations of geostatistics; a good grasp of its possibilities and limitations; and reasonable familiarity with the SGeMS software, thus opening the possibility of practically applying geostatistics.

Open-File Report↗

Numbers and presence of guarding dogs affect wolf and leopard predation on livestock in northeastern Iran

Livestock predation can pose socio-economic impacts on rural livelihoods and is the main cause of retaliatory killings of carnivores in many countries. Therefore, appropriate interventions to reduce livestock predation, lower conflict and promote coexistence are needed. Livestock guarding dogs have been traditionally used to reduce predation, yet details regarding the use of dogs, especially the number of dogs per herd effectively required, are rarely studied. In this study, we assessed how the number and presence of guarding dogs in a herd can reduce livestock losses to leopard and wolf in corrals at night and on grazing grounds in day-time. Using systematic interview surveys (2016-2019), we documented sheep/goat losses per attack (predation rates) from 139 shepherds across 32 villages around Golestan National Park, Iran. We analysed the effects of the number of dogs, presence of dogs, presence of shepherds, seasons, corral quality, livestock number, dog size, distance to villages and distance to reserve on predation rates using generalized linear models. For the leopard model, dog presence significantly decreased ( β = –1.80, 95% confidence interval –2.61 to –0.81) predation rates during day-time to 1.41 individuals per attack. For wolf attacks in corrals at night, predation rates significantly decreased ( β = –0.29, –0.54 to –0.04) with increasing dog numbers. Also, shepherd presence ( β = –0.56, –1.10 to –0.10) and herd size (β = –0.36, –0.60 to –0.12) significantly reduced predation rates. In the wolf day-time model, shepherd presence significantly decreased ( β = –0.93, –1.74 to –0.10) predation rates. Our study suggests that (1) using dogs can reduce, but not eliminate, predation by leopards during day-time; (2) with every additional dog, predation rates by wolves in corrals at night are likely to decrease on average by 25.2%; and (3) the presence of shepherds in corrals at night and during day-time can reduce predation rates.

Azizabad No-Hunting Area, Golestan National Park↗

Technique for the extraction and partial chemical analysis of fluid-filled inclusions from minerals

A method has been developed for the extraction and limited chemical analysis of the materials in solution in the fluid from the very minute fluid-filled inclusions such as commonly occur in whitish or milky quartz. The method may also be applied, with some reservations, to a variety of other minerals. As the amounts of substances in the fluid are small compared with possible contaminants, great care is needed in sample purification and cleaning; an electrolytic method has been found to be the only satisfactory final cleaning step. Following this, the inclusions are opened by ball milling of the cleaned sample, with deionized water, in an alumina ball mill using alumina grinding media. The ions present in the resultant slurry are separated from the ground quartz by electrodialysis and analyzed. Other methods, such as the decrepitation of a sample in an absorption train, are used to estimate the amount of H2O and CO2 in the inclusions. The most significant part of the analytical work has been to determine the ratios between the alkali metals. The materials in solution in the fluid-filled inclusions from 11 samples of quartz have been analyzed by the ball milling-electrodialysis method. Although there are large differences between samples, the average weights of the alkali metal ions found, for all samples, in milligrams per kilogram of quartz, are: Li+-0.92, Na+-99, K+-133, Rb+-0.45, Cs+-0.38 (atomic ratios, in the same sequence: 0.03/1.00/0.79/0.001/0.0007). In addition, 27 to 193 milligrams of CI", and 5 to 140 milligrams of SOr were found, per kilogram of quartz. H2O and COa were determined on only one sample. Six of the samples were from gold-quartz veins in the Grass Valley district, California. In these six the Na+/K+ ratios were all very similar, but the amounts of Li+, Rb+ and Cs+ found varied greatly. Although there are serious limitations to this and to all other techniques developed, it is felt that the results presented are sufficiently encouraging to warrant further study and possible application to specific geologic problems, such as the identification of epochs of quartz deposition

Economic Geology↗

Relative influence of human harvest, carnivores, and weather on adult female elk survival across western North America

Well-informed management of harvested species requires understanding how changing ecological conditions affect demography and population dynamics, information that is lacking for many species. We have limited understanding of the relative influence of carnivores, harvest, weather and forage availability on elk Cervus elaphus demography, despite the ecological and economic importance of this species. We assessed adult female survival, a key vital rate for population dynamics, from 2746 radio-collared elk in 45 populations across western North America that experience wide variation in carnivore assemblage, harvest, weather and habitat conditions. Proportional hazard analysis revealed that 'baseline' (i.e. not related to human factors) mortality was higher with very high winter precipitation, particularly in populations sympatric with wolves Canis lupus. Mortality may increase via nutritional stress and heightened vulnerability to predation in snowy winters. Baseline mortality was unrelated to puma Puma concolor presence, forest cover or summer forage productivity. Cause-specific mortality analyses showed that wolves and all carnivore species combined had additive effects on baseline elk mortality, but only reduced survival by <2%. When human factors were included, ‘total’ adult mortality was solely related to harvest; the influence of native carnivores was compensatory. Annual total mortality rates were lowest in populations sympatric with both pumas and wolves because managers reduced female harvest in areas with abundant or diverse carnivores. Mortality from native carnivores peaked in late winter and early spring, while harvest-induced mortality peaked in autumn. The strong peak in harvest-induced mortality during the autumn hunting season decreased as the number of native carnivore species increased. Synthesis and applications. Elevated baseline adult female elk mortality from wolves in years with high winter precipitation could affect elk abundance as winters across the western US become drier and wolves recolonize portions of the region. In the absence of human harvest, wolves had additive, although limited, effects on mortality. However, human harvest, and its apparent use by managers to offset predation, primarily controls overall variation in adult female mortality. Altering harvest quotas is thus a strong tool for offsetting impacts of carnivore recolonization and shifting weather patterns on elk across western North America.

Colorado;Idaho;Montana;Utah;Washington;Wyoming↗

Harvesting wildlife affected by climate change: a modelling and management approach for polar bears

The conservation of many wildlife species requires understanding the demographic effects of climate change, including interactions between climate change and harvest, which can provide cultural, nutritional or economic value to humans. We present a demographic model that is based on the polar bear Ursus maritimus life cycle and includes density-dependent relationships linking vital rates to environmental carrying capacity ( K ). Using this model, we develop a state-dependent management framework to calculate a harvest level that (i) maintains a population above its maximum net productivity level (MNPL; the population size that produces the greatest net increment in abundance) relative to a changing K , and (ii) has a limited negative effect on population persistence. Our density-dependent relationships suggest that MNPL for polar bears occurs at approximately 0·69 (95% CI = 0·63–0·74) of K . Population growth rate at MNPL was approximately 0·82 (95% CI = 0·79–0·84) of the maximum intrinsic growth rate, suggesting relatively strong compensation for human-caused mortality. Our findings indicate that it is possible to minimize the demographic risks of harvest under climate change, including the risk that harvest will accelerate population declines driven by loss of the polar bear's sea-ice habitat. This requires that (i) the harvest rate – which could be 0 in some situations – accounts for a population's intrinsic growth rate, (ii) the harvest rate accounts for the quality of population data (e.g. lower harvest when uncertainty is large), and (iii) the harvest level is obtained by multiplying the harvest rate by an updated estimate of population size. Environmental variability, the sex and age of removed animals and risk tolerance can also affect the harvest rate. Synthesis and applications . We present a coupled modelling and management approach for wildlife that accounts for climate change and can be used to balance trade-offs among multiple conservation goals. In our example application to polar bears experiencing sea-ice loss, the goals are to maintain population viability while providing continued opportunities for subsistence harvest. Our approach may be relevant to other species for which near-term management is focused on human factors that directly influence population dynamics within the broader context of climate-induced habitat degradation.

Journal of Applied Ecology↗

The HayWired Earthquake Scenario

Foreword The 1906 Great San Francisco earthquake (magnitude 7.8) and the 1989 Loma Prieta earthquake (magnitude 6.9) each motivated residents of the San Francisco Bay region to build countermeasures to earthquakes into the fabric of the region. Since Loma Prieta, bay-region communities, governments, and utilities have invested tens of billions of dollars in seismic upgrades and retrofits and replacements of older buildings and infrastructure. Innovation and state-of-the-art engineering, informed by science, including novel seismic-hazard assessments, have been applied to the challenge of increasing seismic resilience throughout the bay region. However, as long as people live and work in seismically vulnerable buildings or rely on seismically vulnerable transportation and utilities, more work remains to be done. With that in mind, the U.S. Geological Survey (USGS) and its partners developed the HayWired scenario as a tool to enable further actions that can change the outcome when the next major earthquake strikes. By illuminating the likely impacts to the present-day built environment, well-constructed scenarios can and have spurred officials and citizens to take steps that change the outcomes the scenario describes, whether used to guide more realistic response and recovery exercises or to launch mitigation measures that will reduce future risk. The HayWired scenario is the latest in a series of like-minded efforts to bring a special focus onto the impacts that could occur when the Hayward Fault again ruptures through the east side of the San Francisco Bay region as it last did in 1868. Cities in the east bay along the Richmond, Oakland, and Fremont corridor would be hit hardest by earthquake ground shaking, surface fault rupture, aftershocks, and fault afterslip, but the impacts would reach throughout the bay region and far beyond. The HayWired scenario name reflects our increased reliance on the Internet and telecommunications and also alludes to the interconnectedness of infrastructure, society, and our economy. How would this earthquake scenario, striking close to Silicon Valley, impact our interconnected world in ways and at a scale we have not experienced in any previous domestic earthquake? The area of present-day Contra Costa, Alameda, and Santa Clara Counties contended with a magnitude-6.8 earthquake in 1868 on the Hayward Fault. Although sparsely populated then, about 30 people were killed and extensive property damage resulted. The question of what an earthquake like that would do today has been examined before and is now revisited in the HayWired scenario. Scientists have documented a series of prehistoric earthquakes on the Hayward Fault and are confident that the threat of a future earthquake, like that modeled in the HayWired scenario, is real and could happen at any time. The team assembled to build this scenario has brought innovative new approaches to examining the natural hazards, impacts, and consequences of such an event. Such an earthquake would also be accompanied by widespread liquefaction and landslides, which are treated in greater detail than ever before. The team also considers how the now-prototype ShakeAlert earthquake early warning system could provide useful public alerts and automatic actions. Scientific Investigations Report 2017–5013 and accompanying data releases are the products of an effort led by the USGS, but this body of work was created through the combined efforts of a large team including partners who have come together to form the HayWired Coalition (see chapter A). Use of the HayWired scenario has already begun. More than a full year of intensive partner engagement, beginning in April 2017, is being directed toward producing the most in-depth look ever at the impacts and consequences of a large earthquake on the Hayward Fault. With the HayWired scenario, our hope is to encourage and support the active ongoing engagement of the entire community of the San Francisco Bay region by providing the scientific, engineering, and economic and social science inputs for use in exercises and planning well into the future. As HayWired volumes are published, they will be made available at https://doi.org/10.3133/sir20175013 .

California↗

Oyster model inventory: Identifying critical data and modeling approaches to support restoration of oyster reefs in coastal U.S. Gulf of Mexico waters

Executive Summary Along the coast of the U.S. Gulf of Mexico, the eastern oyster ( Crassostrea virginica ) plays important ecological and economic roles. Commercial landings from this region account for more than 50 percent of all U.S. landings; these oyster reefs also provide varied ecosystem services, including nursery habitat for many fish and macroinvertebrate species, shoreline protection, and water-quality maintenance. Declining trends in both total oyster production and functional reef area across this region have spurred investment in restoration of oyster resources, with specific calls for restoration projects to develop a network of reefs and identify broodstock and sanctuary reef restoration sites. Decision making related to restoration and establishment of a network of oyster reefs in the Gulf of Mexico requires information on both the environment and the effects of the environment on the oyster life cycle (including larval movement, survival, oyster recruitment, reproduction, growth, and mortality). Here, we examined the current state of data and model development in this region with the goal of providing an overview of oyster modeling approaches and an inventory of available data and existing oyster models. This report is meant to provide an overview to managers for understanding existing efforts and identify a path forward to most efficiently inform oyster resource management and restoration planning in moving from a single reef management approach to a reef network management approach. Numerous models related to some aspect of the oyster life cycle have been built, calibrated, and validated for various Gulf of Mexico estuaries over the last few decades (over 30 models identified). These models, which could inform site restoration, can be classified into four approaches: (1) oyster Habitat Suitability Index (HSI) models; (2) larval transport models; (3) on-reef oyster models that may include oyster growth, mortality and reproduction, and substrate persistence; and (4) coupled larval transport on-reef metapopulation models that simulate the entire oyster life cycle. The data requirements, model complexity and assumptions, and transferability vary by approach. Specifically, some approaches may offer greater accessibility, flexibility, and transferability spatially or temporally, with minimal data input, but only provide broad information to support site selection. In contrast, other approaches may require significant site-specific data for their construction and validation but may provide more accurate and location-specific data to support site selection for broodstock reefs. Regardless of modeling approach used, data on environmental drivers, such as salinity, water temperature, or water flow impacting oyster metabolism and movement, are required at appropriate spatial and temporal scales. While numerous data collection platforms, environmental models, and research products exist within Gulf of Mexico estuaries to provide important environmental data to use as drivers in the oyster models, significant variability in temporal and spatial coverage of the data, and variation in the availability of future condition models, exists across estuaries. This variation influences the spatial and temporal scales at which oyster models may be developed and impacts the calibration and validation of the oyster models within a given estuary, affecting its potential ability to address specific management or restoration questions. While multiple modeling approaches exist for informing site selection of broodstock or sanctuary oyster reefs, the development, calibration, and validation of a single modeling platform presents the most efficient, transferable, and useful tool for managers across the Gulf of Mexico. The development of a single modeling platform would involve using standardized input variables, governing equations, and assumptions for the modeled oyster processes and outputs, and for standardized calibration and validation procedures that could be applied within each estuary. The differences among estuary applications would require substituting only estuary-specific environmental data, and calibrating and validating the modeling approach with local oyster data. Two modeling approaches likely to be useful include (1) development of a general geospatial HSI modeling framework that could be applied consistently across estuaries and (2) a mechanistic coupled larval transport on-reef metapopulation model requiring only estuarine specific calibration and hydrodynamic models. Both approaches benefit from existing work across multiple Gulf of Mexico estuaries and could provide valuable support for oyster restoration, but may differ in their ability to address specific questions related to oyster restoration. HSI models specifically guide restoration practitioners in determining suitable habitat based on available data. The HSI approach, while currently more widely used and accessible, requires more development of larval suitability and larval input and output components in order to inform reef connectivity. A metapopulation approach considering the full oyster life cycle that simulates both on-reef oyster growth, mortality, reproduction, substrate persistence, and larval transport (ideally with larval growth and mortality) would provide the greatest detail and level of understanding but requires significant up-front investment. The larval oyster model and on-reef oyster model are usually developed independently for systems, although the two approaches can be coupled to represent the entire oyster life cycle in order to characterize and assess a reef metapopulation. This approach may be less accessible and much more data-intensive, however, and it requires some expertise to run and apply to inform oyster resource management. Ultimately, the development of single modeling platforms for each of these approaches would provide flexible tools applicable across all Gulf of Mexico oyster supporting estuaries. By using a single platform for model development, testing, calibrating and validating, and evaluation of modeled future scenarios, oyster restoration scientists and managers would not only be able to examine different scenario outcomes within a single estuary, but could also have comparable modeled results to evaluate potential outcomes, across estuaries and regions, that are not confounded by varying modeled data inputs, governing equations, assumptions, or user judgement.

Alabama, Florida, Louisiana, Mississippi, Texas↗

Multiscale guidance and tools for implementing a landscape approach to resource management in the Bureau of Land Management

The Bureau of Land Management (BLM) is implementing a landscape approach to resource management (hereafter, landscape approach) to more effectively work with partners and understand the effects of management decisions. A landscape approach is a set of concepts and principles used to guide resource management when multiple stakeholders are involved and goals include diverse and sustainable social, environmental, and economic outcomes. Core principles of a landscape approach include seeking meaningful participation of diverse stakeholders, considering diverse resource values in multifunctional landscapes, acknowledging the tradeoffs needed to meet diverse objectives in the context of sustainable resource management, and addressing the complexity of social and ecological processes by embracing interdisciplinarity and considering multiple and broad spatial and temporal perspectives. In chapter 1, we outline the overall goal of this report: to provide a conceptual foundation and framework for implementing a landscape approach to resource management in the BLM, focusing on the role of multiscale natural resource monitoring and assessment information. In chapter 2, we describe a landscape approach to resource management. BLM actions taken to implement a landscape approach include a major effort to compile broad-scale data on natural resource status and condition across much of the west. These broadscale data now provide a regional context for interpreting monitoring data collected at individual sites and informing decisions made for local projects. We also illustrate the utility of using multiscale data to understand potential effects of different resource management decisions, define relevant terms in landscape ecology, and identify spatial scales at which planning and management decisions may be evaluated. In chapter 3, we describe how the BLM Rapid Ecoregional Assessment program and Assessment, Inventory and Monitoring program may be integrated to provide the multiscale monitoring data needed to inform a landscape approach. We propose six core, broad-scale indicators of natural resource status and condition: the amount, spatial distribution, patch size and connectivity of ecosystems and wildlife habitats, and the pattern of existing development across the landscape. Additional supplemental broad-scale indicators may include fire return intervals, distributions of invasive species, and vulnerability of ecosystems to a changing climate. Landscape intactness is an additional derived indicator that is calculated from one or more of the core and supplemental broad-scale indicators. We then outline a process for assessing broad-scale indicators that is consistent with the overall BLM Assessment, Inventory, and Monitoring process, facilitating development of a multiscale natural resource monitoring program. Finally, we describe how broad-scale indicators of natural resource status and condition may guide field monitoring implemented through the BLM Assessment, Inventory and Monitoring program and help address complex management questions. In chapter 4, we consider the specific question of assessing the ecological integrity of rangelands across the western United States. We first define ecological integrity and its relation to land health. We then suggest that a combination of six local-scale indicators collected through field sampling at individual sites and five complementary broad-scale indicators together provide information on the composition, structure, and function of rangelands. The terrestrial monitoring indicators collected at the level of individual field sites are the amount of bare ground, vegetation composition (including invasive plants and plants of management concern), vegetation height, and the proportion of the soil surface in large intercanopy gaps. The broad-scale indicators are vegetation amount, distribution, patch size, connectivity, and productivity, along with the pattern of terrestrial development. Our suggested approach to quantifying ecological integrity focuses specifically on informing management of public lands for multiple resource uses, and illustrates how existing data collected through BLM monitoring and assessment programs may be used together to provide multiscale information on land condition across broad extents. In chapter 5, we develop a method for quantifying landscape intactness and apply this method to the western United States. Our multiscale index of landscape intactness is designed to be defensible, decomposable, and easy to understand. The foundation of the multiscale index of landscape intactness is the surface disturbance footprint of anthropogenic development, including energy and urban development, roads and railroads, cultivated croplands, surface mines and quarries, and energy transmission lines and pipelines. The index represents a gradient of anthropogenic influence as represented by development summarized at two spatial scales of analysis: 2.5 and 20 kilometers. We provide several example applications of the index, illustrating how these data may inform natural resource decisions at the spatial extent of BLM field and district offices, states, ecoregions, and the western United States. We find that 19.2 percent of lands managed by the BLM across the 17 western states of the conterminous United States had the highest landscape intactness. The largest intact areas occur on public lands at high elevations or in the Great Basin. We believe the frameworks, processes, and analyses provided in this report will improve the ability of the BLM to identify and evaluate potential direct and indirect effects of management actions (such as habitat restoration and renewable energy development), and assist the BLM in further implementing a landscape approach to resource management.

Open-File Report↗

Cadmium isotope fractionation during coal combustion: Insights from two U.S. coal-fired power plants

Coal combustion, one of the principal energy sources of electricity in the United States, produces over 100 million tons of coal combustion products (CCPs) per year in the U.S. The reuse and disposal of CCPs has the potential to release toxic trace elements, including cadmium (Cd), into the environment. In this study, we investigated CCPs, including bottom ash (BA), economizer fly ash (EFA), and fly ash (FA), as well as feed coal (FC) and pulverized coal (PC) collected from two U.S. coal-fired power plants in New Mexico and Ohio with different coal supplies. The New Mexico plant uses high volatile C bituminous, low-sulfur coals mined from the San Juan Basin (Cretaceous Fruitland Formation) and the Ohio plant uses high volatile A bituminous, high-sulfur central Appalachian Basin coals (Upper Pennsylvanian Monongahela Formation). Mineralogical and elemental analysis showed that these CCP samples consist of ∼70% amorphous Al-Si-rich glasses and ∼30% mineral phases of quartz (SiO 2 ) and mullite (Ai 6 Si 2 O 13 ). The Cd isotope compositions (δ 114 Cd, normalized to NIST Cd standard 3108) of FA and EFA samples (ranging from −0.51 to +0.47‰) are distinctively heavier than those of BA samples (−0.75 to −0.52‰) in both power plants. We interpret this Cd isotope difference as a result of Cd condensation from the gas phase during flue gas cooling, instead of evaporation of Cd phase during coal combustion. Cd condensation is the main process to generate the isotopically heavy Cd signatures that preferentially partition on the fine FA particles. We also investigated Cd isotope compositions in different leachate products from a series of batch-leaching experiments with these CCPs, using diluted acetic acid, hydroxyl ammonium chloride, hydrogen peroxide followed by ammonium acetate, and 5% nitric acid, as a possible means to identify CCP-released Cd in the environment. Unusually and significantly heavier Cd isotope compositions were observed in each leachate of FA samples (+1.10 to +7.09‰), which fall far outside from the range of Cd isotope ratios observed in natural soils and rocks, but less so for the EFA samples (−0.43 to +1.18‰). Such an observation is consistent with the interpretation that isotopically heavy Cd preferentially partitions on the fine FA particles after coal combustion and is readily to be released during these leaching experiments. This study demonstrates that high-temperature coal combustion can lead to a very large degree of fractionation of Cd isotopes that can be used as a unique tracer for identifying anthropogenic metal inputs in the environment. The major Cd isotope fractionation process occurs as the Cd gas phase condenses on fine FA particles during the flue gas cooling stage after coal combustion.

New Mexico, Ohio↗

Calcite U-Pb dating unravels the age and hydrothermal history of the giant Shuiyindong Carlin-type gold deposit in the Golden Triangle, south China

The ages of Carlin-type gold deposits in the Golden Triangle of South China have long been questioned due to the general lack of minerals unequivocally linked to gold deposition that can be precisely dated using conventional radiogenic isotope techniques. Recent advances in U-Pb methods show that calcite can be used to constrain the ages of hydrothermal processes, but few studies have been applied to ore deposits. Herein, we show that this approach can be used to constrain the timing of hydrothermal activity that generated and overprinted the giant Shuiyindong Carlin-type gold deposit in the Golden Triangle. Three stages of calcite (Cal-1, Cal-2, and Cal-3) have been recognized in this deposit based on crosscutting relationships, cathodoluminescence colors, and chemical (U, Pb, and rare earth element [REE]) and isotope (C, O, Sr) compositions. Cal-1 is texturally associated with ore-stage jasperoid and disseminated Au-bearing arsenian pyrite in hydrothermally altered carbonate rocks, which suggests it is synmineralization. Cal-2 fills open spaces and has a distinct orange cathodoluminescence, suggesting that it precipitated during a second fluid pulse. Cal-1 and Cal-2 have similar carbonate rock-buffered chemical and isotopic compositions. Cal-3 occurs in veins that often contain realgar and/or orpiment and are chemically (low U, Pb, and REE) and isotopically (higher δ 13 C, lower δ 18 O and Sr i values) distinct from Cal-1 and Cal-2, suggesting that it formed from a third fluid. U-Pb isotope analyses, by laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) for U-rich Cal-1 and Cal-2 and by LA-multicollector (MC)-ICP-MS for U-poor Cal-3, yield well-defined age constraints of 204.3 to 202.6, 191.9, and 139.3 to 137.1 Ma for Cal-1, Cal-2, and Cal-3, respectively. These new ages suggest that the Shuiyindong gold deposit formed in the late Triassic and was overprinted by hydrothermal events in the early Jurassic and early Cretaceous. Given the association of Cal-3 with orpiment and realgar, and previous geochronologic studies of several other major gold deposits in the Golden Triangle, we infer that the latest stage of calcite may be associated with an early Cretaceous regional gold metallogenic event. Combined with existing isotopic ages in the region, these new ages lead us to propose that Carlin-type gold deposits in the Golden Triangle formed during two metallogenic episodes in extensional settings, associated with the late Triassic Indochina orogeny and early Cretaceous paleo-Pacific plate subduction. This study shows that the calcite U-Pb method can be used to constrain the timing of Carlin-type gold deposits and successive hydrothermal events.

Golden Triangle, south China↗

Anthropogenic influences on the input and biogeochemical cycling of nutrients and mercury in Great Salt Lake, Utah, USA

Despite the ecological and economic importance of Great Salt Lake (GSL), little is known about the input and biogeochemical cycling of nutrients and trace elements in the lake. In response to increasing public concern regarding anthropogenic inputs to the GSL ecosystem, the US Geological Survey (USGS) and US Fish and Wildlife Service (USFWS) initiated coordinated studies to quantify and evaluate the significance of nutrient and Hg inputs into GSL. A 6??? decrease in ??15N observed in brine shrimp (Artemia franciscana) samples collected from GSL during summer time periods is likely due to the consumption of cyanobacteria produced in freshwater bays entering the lake. Supporting data collected from the outflow of Farmington Bay indicates decreasing trends in ??15N in particulate organic matter (POM) during the mid-summer time period, reflective of increasing proportions of cyanobacteria in algae exported to GSL on a seasonal basis. The C:N molar ratio of POM in outflow from Farmington Bay decreases during the summer period, supportive of the increased activity of N fixation indicated by decreasing ??15N in brine shrimp and POM. Although N fixation is only taking place in the relatively freshwater inflows to GSL, data indicate that influx of fresh water influences large areas of the lake. Separation of GSL into two distinct hydrologic and geochemical systems from the construction of a railroad causeway in the late 1950s has created a persistent and widespread anoxic layer in the southern part of GSL. This anoxic layer, referred to as the deep brine layer (DBL), has high rates of SO42 - reduction, likely increasing the Hg methylation capacity. High concentrations of methyl mercury (CH3Hg) (median concentration = 24 ng/L) were observed in the DBL with a significant proportion (31-60%) of total Hg in the CH3Hg form. Hydroacoustic and sediment-trap evidence indicate that turbulence introduced by internal waves generated during sustained wind events can temporarily mix the elevated CH3Hg concentrations in the DBL with the more biologically active upper brine layer (UBL). Brine shrimp collected during the summer/fall time periods contained elevated Hg concentrations (median concentration = 0.34 mg/kg, dry weight (dw)) relative to samples collected during the spring (median concentration < 0.2 mg/kg, dw). Higher Hg in brine shrimp during the summer and fall may reflect the higher proportion of adult brine shrimp during this time period, resulting in an increased time for bioaccumulation of Hg. Eared grebes (Podiceps nigricollis) consume brine shrimp from GSL during the fall molting period. Median Hg concentrations in eared grebe livers increased by almost three times during the 3-5 month fall molting period. Selected duck species utilizing GSL have consistently exceeded the US Environmental Protection Agency (USEPA) screening level for Hg (0.3 mg/kg Hg wet weight), resulting in the issuance of warnings against unlimited human consumption of breast muscle tissue.

Applied Geochemistry↗

The origin of jasperoid in limestone

The name jasperoid has been applied to rocks that consist mainly of silica and that have formed by replacement. This paper considers only those jasperoids formed by replacement of limestone . Major problems involved in the origin of such jasperoid include: source of the silica; nature of solutions that dissolve, transport, and precipitate silica; and the mechanism of replacement of limestone by silica. The answers to these problems are of practical as well as scientific interest because many jasperoid bodies are closely related to mineralization. Silica may be derived from: juvenile silica of magmatic origin ; silica leached from underlying rocks by hydrothermal solutions; silica locally derived from enclosing rocks by circulating solutions; and silica carried downward in ground water from the weathering of overlying rocks. The nature and the concentration of other substances in solutions influence, in a complex manner, the ability of these solutions to dissolve, transport, and precipitate silica. Nevertheless, the following generalizations can be made. The solubility and rate of solution of silica in water at moderate pressure increase slowly with temperature up to 200 degrees C.; from 200 degrees to 360 degrees C. they increase rapidly; above 360 degrees C. solubility is pressure dependent, increasing steadily at high pressure and decreasing slightly at moderate pressure due to the formation of a vapor phase. The pH has slight effect on the ionic solubility of silica in the range from pH 1 to ph 9 at low temperature. The effect of other components on the solubility of silica is probably subordinate to that of temperature above 200 degrees C., but becomes increasingly important as the temperature falls below that point. Most jasperoid bodies form by both replacement and silica deposition in voids, with replacement dominant during the early phase, and precipitation dominant during later phases. Replacement of limestone by silica is favored by relatively low temperature acid solutions, and the presence of CO 2 . As limestone dissolves, Ca ions are released to promote the precipitation of colloidal silica. Acid solutions then diffuse through this gelatinous film to continue dissolving limestone behind it; the Ca ions diffusing outward cause the precipitation of more colloidal silica at the solution-gel interface. As the gel mass ages, it shrinks, hardens, and ruptures. More silica is then deposited in the fractures. Eventually the gel crystallizes to a dense mass of aphanitic quartz and chalcedony, with shrinkage cracks and vugs filled or coated with younger coarse-grained quartz and other minerals that have been deposited directly from solution. The theory that relatively low temperature favors the formation of jasperoid replacement bodies in carbonate rocks, and high temperature inhibits their formation, offers an explanation for the gap that is observed in some districts between contact metasomatic lime silicates and siliceous replacement of limestone . This gap is characterized by the lack of any reaction between limestone and silica-bearing solutions moving through it.

Economic Geology↗

Cascadia Margin cold seeps: Subduction zone fluids, gas hydrates, and chemosynthetic habitats

Priority Geographic Area: The outer continental shelf and upper continental slope from Canada/U.S. border offshore Washington State to the Mendocino Fracture Zone (Northern California), entirely within the U.S. Exclusive Economic Zone (EEZ), from the outermost shelf to at least 2000 m water depth (Figure 1). Description of Priority Area: Since 2015, over a thousand water column gas plumes originating at seafloor gas seeps have been discovered landward of the Cascadia deformation front (e.g., Embley et al., 2016; Johnson et al., 2015, 2019; Merle and Embley, 2016; NA-95 Cruise Report, 2018; Riedel et al., 2018), adding to those that had long been known on Hydrate Ridge (e.g., Heeschen et al., 2003; Tréhu et al., 2004). The recently-discovered seeps stretch from offshore Vancouver Island to the Mendocino Fracture Zone and from the outer shelf to ~2000 m water depth, occurring both landward and seaward of the nominal limit for gas hydrate stability zone on the upper continental slope (Figure 1). Hundreds of seeps likely remain undiscovered. Water column imaging is incomplete both within the target geographic area and farther seaward, between the 2000 m isobath and the deformation front, which is the subject of an imaging study described in a white paper by Watt et al. The recently-discovered Cascadia Margin cold seeps partially overlap an important active margin gas hydrate province (Spence et al., 2001; Tréhu et al., 2003, 2004), as well as an area where sediments on the North American plate are folded and faulted and affected by fluids generated in the subduction complex beneath the Cascadia forearc (e.g., Saffer and Tobin, 2011). Several Ocean Drilling Program expeditions have focused on hydrate systems offshore Vancouver and Oregon (e.g., Riedel et al., 2009; Tréhu et al., 2004) and on the connection between the shallow and deep hydrogeologic systems. Cabled observatories now continuously monitor physical, chemical, and venting processes on south Hydrate Ridge (OOI; e.g., Philip et al., 2016a) and offshore Vancouver Island (NEPTUNE; e.g. Römer et al., 2016). Outside of these well-studied gas hydrate areas, a subset of the recently-discovered Cascadia seeps, including some that we visited with R/V Falkor in 2019 (e.g., https://schmidtocean.org/cruise/methane-seeps-at-edge-of-hydrate-stability/), also likely emit methane associated with shallow subseafloor gas hydrate systems. Other seeps are delivering not only methane, but also deep-derived gases (Baumberger et al., 2018, 2020) to the seafloor. Many Cascadia Margin seeps have also been recognized at water depths too shallow (e.g., 175 m) to be connected to gas hydrate dynamics. These seeps are postulated to be emitting gas and fluids that originated deep in accretionary wedge before migrating up normal faults generated during forearc extension associated with large earthquakes (Johnson et al., 2019). Only a small fraction of the recently discovered U.S. Cascadia Margin water column gas plumes has so far been verified by ROVs (Hercules from E/V Nautilus in 2016 and 2018; SuBastian from R/V Falkor in 2018 and 2019) to correspond to seafloor seeps. Careful scientific mapping, investigation, and sampling at the seeps have also been limited (e.g., Baumberger et al., 2018, 2020; Merle and Embley, 2016; Seabrook et al., 2018; Greinert et al. 2019). This white paper focuses on expanding exploration of already-identified U.S. Cascadia Margin cold seeps through a multipronged and multidisciplinary discovery program that could be accomplished with a variety of NOAA assets. The goals of the proposed exploration activities are to develop high-resolution maps of seep fields from deep ocean vehicles; to verify (and sample) seafloor gas emissions at the locations of water column plumes for compositional and isotopic studies; to map, sample, and conduct analyses on chemosynthetic communities and deep-sea coral habitats near seep sites to document species distributions and habitats as a function of depth and latitude along the margin; to collect seep geologic samples that can constrain the timing of methane emissions through geochronology; and to record environmental data (e.g., CTD) near the seafloor and in the water column above the seeps. Seafloor mapping using shipboard systems (multibeam/backscatter) would be needed to characterize seafloor features near seep sites. Water column imaging (EK60/80 and/or multibeam WCD data) conducted before and after seafloor explorations would capture active methane plumes and constrain temporal variations in seep emissions (e.g., Kannberg et al., 2013; Philip et al., 2016a, 2016b), which are known to vary on time scales as rapid as tidal cycles on this margin (e.g., Römer et al., 2016). What are the characterization and data needs in this area? Check all that apply: __x_ Biology, Geology, Physical Oceanography, Chemistry ___ Marine Archaeology ___ Other Provide a list or brief description of the data needed within this area, from your perspective: 1. Water column backscatter to image active gas plumes 2. High-resolution multibeam bathymetry, seafloor backscatter, and shallow sub-bottom imaging 3. Visual characterization and ground truthing of potential seeps, including high-resolution mapping and photography from near-seafloor vehicles; collection of seep-associated species, corals, sediments, authigenic carbonates, gases, and seawater Describe relevance to national security, conservation, and/or the economy: The Cascadia margin seeps provide significant ecosystem services, including habitat for commercially important fishes and support for diversity along the continental margin. Methane seeps are also biological hotspots for krill, plankton, and crustaceans, which in turn sustain higher trophic levels (e.g., whales). Methane-derived authigenic carbonates serve as a hard substrate for deep-sea corals and sponges on millennial time scales. The studies proposed here will elucidate the relationship among seep environments, deep-sea corals, sponges, fisheries, and other organisms and provide new insight into subduction zone and hydrate-associated fluids in this important seismogenic zone. The studies address fishery management concerns and inform future conservation of sensitive species (e.g., deep-sea corals) and benthic habitats. From your perspective, what makes this area unique? The Cascadia Margin seeps are a critical component of the leaky margin that stretches from Baja California to the Aleutian Arc along the Pacific coastline of North America. Cold seeps have been intensely studied on the Gulf of Mexico and U.S. Atlantic passive margins with a focus on chemosynthetic communities, deep-sea corals, and leakage of microbially-generated and/or thermogenic hydrocarbons; however, the recently-discovered Cascadia Margin seeps, as well as active margin seep systems in general, remain more poorly characterized. Such seeps not only contribute to the ocean carbon cycle (e.g., Pohlman et al., 2011), thereby fueling the base of the food chain in these settings, but also emit subduction zone fluids that provide clues about processes within the seismogenic zone and the accretionary complex. The Cascadia seeps area allows both biological (e.g., benthic habitats, coral distributions) and physical processes (e.g., generation of subduction zone fluids) to be studied along both depth (perpendicular to the deformation front) and latitudinal gradients.

California, Oregon, Washington↗

Scenario earthquake hazards for the Long Valley Caldera-Mono Lake area, east-central California (ver. 2.0, January 2018)

As part of the U.S. Geological Survey’s (USGS) multi-hazards project in the Long Valley Caldera-Mono Lake area, the California Geological Survey (CGS) developed several earthquake scenarios and evaluated potential seismic hazards, including ground shaking, surface fault rupture, liquefaction, and landslide hazards associated with these earthquake scenarios. The results of these analyses can be useful in estimating the extent of potential damage and economic losses because of potential earthquakes and also for preparing emergency response plans. The Long Valley Caldera-Mono Lake area has numerous active faults. Five of these faults or fault zones are considered capable of producing magnitude ≥6.7 earthquakes according to the Uniform California Earthquake Rupture Forecast, Version 2 (UCERF 2) developed by the 2007 Working Group on California Earthquake Probabilities (WGCEP) and the USGS National Seismic Hazard Mapping Program. These five faults are the Fish Slough, Hartley Springs, Hilton Creek, Mono Lake, and Round Valley Faults. CGS developed earthquake scenarios for these five faults in the study area and for the White Mountains Fault Zone to the east of the study area. In this report, an earthquake scenario is intended to depict the potential consequences of significant earthquakes. A scenario earthquake is not necessarily the largest or most damaging earthquake possible on a recognized fault. Rather it is both large enough and likely enough that emergency planners should consider it in regional emergency response plans. In particular, the ground motion predicted for a given scenario earthquake does not represent a full probabilistic hazard assessment, and thus it does not provide the basis for hazard zoning and earthquake-resistant building design. Earthquake scenarios presented here are based on fault geometry and activity data developed by the WGCEP, and are consistent with the 2008 Update of the United States National Seismic Hazard Maps (NSHM). Alternatives to the NSHM scenario were developed for the Hilton Creek and Hartley Springs Faults to account for different opinions in how far these two faults extend into Long Valley Caldera. For each scenario, ground motions were calculated using the current standard practice: the deterministic seismic hazard analysis program developed by Art Frankel of USGS and three Next Generation Ground Motion Attenuation (NGA) models. Ground motion calculations incorporated the potential amplification of seismic shaking by near-surface soils defined by a map of the average shear wave velocity in the uppermost 30 m ( V S30 ) developed by CGS. In addition to ground shaking and shaking-related ground failure such as liquefaction and earthquake induced landslides, earthquakes cause surface rupture displacement, which can lead to severe damage of buildings and lifelines. For each earthquake scenario, potential surface fault displacements are estimated using deterministic and probabilistic approaches. Liquefaction occurs when saturated sediments lose their strength because of ground shaking. Zones of potential liquefaction are mapped by incorporating areas where loose sandy sediments, shallow groundwater, and strong earthquake shaking coincide in the earthquake scenario. The process for defining zones of potential landslide and rockfall incorporates rock strength, surface slope, and existing landslides, with ground motions caused by the scenario earthquake. Each scenario is illustrated with maps of seismic shaking potential and fault displacement, liquefaction, and landslide potential. Seismic shaking is depicted by the distribution of shaking intensity, peak ground acceleration, and 1.0-second spectral acceleration. One-second spectral acceleration correlates well with structural damage to surface facilities. Acceleration greater than 0.2 g is often associated with strong ground shaking and may cause moderate to heavy damage. The extent of strong shaking is influenced by subsurface fault dip and near surface materials. Strong shaking is more widespread in the hanging wall regions of a normal fault. Larger ground motions also occur where young alluvial sediments amplify the shaking. Both of these effects can lead to strong shaking that extends farther from the fault on the valley side than on the hill side. The effect of fault rupture displacements may be localized along the surface trace of the mapped earthquake fault if fault geometry is simple and the fault traces are accurately located. However, surface displacement hazards can spread over a few hundred meters to a few kilometers if the earthquake fault has numerous splays or branches, such as the Hilton Creek Fault. Faulting displacements are estimated to be about 1 meter along normal faults in the study area and close to 2 meters along the White Mountains Fault Zone. All scenarios show the possibility of widespread ground failure. Liquefaction damage would likely occur in the areas of higher ground shaking near the faults where there are sandy/silty sediments and the depth to groundwater is 6.1 meters (20 feet) or less. Generally, this means damage is most common near lakes and streams in the areas of strongest shaking. Landslide potential exists throughout the study region. All steep slopes (>30 degrees) present a potential hazard at any level of shaking. Lesser slopes may have landslides within the areas of the higher ground shaking. The landslide hazard zones also are likely sources for snow avalanches during winter months and for large boulders that can be shaken loose and roll hundreds of feet down hill, which happened during the 1980 Mammoth Lakes earthquakes. Whereas methodologies used in estimating ground shaking, liquefaction, and landslides are well developed and have been applied in published hazard maps; methodologies used in estimating surface fault displacement are still being developed. Therefore, this report provides a more in-depth and detailed discussion of methodologies used for deterministic and probabilistic fault displacement hazard analyses for this project.

California↗

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↗

Quantifying resilience

The biosphere is under unprecedented pressure, reflected in rapid changes in our global ecological, social, technological and economic systems. In many cases, ecological and social systems can adapt to these changes over time, but when a critical threshold is surpassed, a system under stress can undergo catastrophic change and reorganize into a different state. The concept of resilience, introduced more than 40 years ago in the ecological sciences, captures the behaviour of systems that can occur in alternative states. The original definition of resilience forwarded by Holling ( 1973 ) is still the most useful. It defines resilience as the amount of disturbance that a system can withstand before it shifts into an alternative stable state. The idea of alternative stable states has clear and profound implications for ecological management. Coral reefs, for example, are high-diversity systems that provide key ecosystem services such as fisheries and coastal protection. Human impacts are causing significant, ongoing reef degradation, and many reefs have shifted from coral- to algal-dominated states in response to anthropogenic pressures such as elevated water temperatures and overfishing. Understanding and differentiating between the factors that help maintain reefs in coral-dominated states vs. those that facilitate a shift to an undesired algal-dominated state is a critical step towards sound management and conservation of these, and other, important social&ndash;ecological systems. Resilience has gained popularity among both academicians and laypeople, as a term meant to describe a systems&rsquo; ability to withstand disturbance. Resilience has become a buzzword in the last decade, as shown by its increasing appearance in calls for research proposals and scientific citation data bases. The term resilience has in many cases lost the clarity of the original definition and in fact is frequently used in a manner in direct opposition to the original definition. Many current uses of the concept are loose and incorrect. The term is becoming increasingly used in a normative sense (Brand & Jax 2007 ), as if resilience were a desirable quality of systems. However, even systems in highly undesirable states, such as macro-algae dominated reefs, or city cores in poverty traps, may be highly resilient, which is to say they withstand attempts to transform them into different (desirable) states. Operationalizing the concept of resilience for application and management has been difficult. Misuse of the term can have significant negative impacts, because resilience is being used to help guide responses to natural disasters and to assess the sustainability of ecosystems and urban systems and has been driving international research priorities. Resilience has been argued to be a basic emergent property of systems, a process or a rate. We focus on the original concept as described by Holling, which is that of an emergent system property; when a system is in a desirable state and managers wish to enhance resilience, or when the system is in an undesirable state and managers wish to erode resilience and foster a transformation to an alternative state. Fostering or eroding resilience is a process. When a system is perturbed but resilience is not exceeded, then the recovery can be measured as a rate. Several frameworks to operationalize resilience have been proposed. A decade ago, a special feature focused on quantifying resilience was published in the journal Ecosystems (Carpenter, Westley & Turner 2005 ). The approach there was towards identifying surrogates of resilience, but few of the papers proposed quantifiable metrics. Consequently, many ecological resilience frameworks remain vague and difficult to quantify, a problem that this special feature aims to address. However, considerable progress has been made during the last decade (e.g. Pope, Allen & Angeler 2014 ). Although some argue that resilience is best kept as an unquantifiable, vague concept (Quinlan et al . 2016 ), to be useful for managers, there must be concrete guidance regarding how and what to manage and how to measure success (Garmestani, Allen & Benson 2013 ; Spears et al . 2015 ). Ideas such as &lsquo;resilience thinking&rsquo; have utility in helping stakeholders conceptualize their systems, but provide little guidance on how to make resilience useful for ecosystem management, other than suggesting an ambiguous, Goldilocks approach of being just right (e.g. diverse, but not too diverse; connected, but not too connected). Here, we clarify some prominent resilience terms and concepts, introduce and synthesize the papers in this special feature on quantifying resilience and identify core unanswered questions related to resilience.

Journal of Applied Ecology↗