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.

At least 397 records · Page 22Linked to original sources
USGS1970

Photogeologic and thermal infrared reconnaissance surveys of the Los Negritos-Ixtlan de los Hervores geothermal area, Michoacan, Mexico

New techniques, involving interpretation of panchromatic, ektachrome and ektachrome infrared aerographic photogaphs and thermographic infrared imagery recording emission from the earth's surface in middle and far infrared wavelengths (3–5μm and 8–14μm), are being introduced in geothermal investigations in Mexico to identify outstanding structural and geologic features in a rapid and economical manner. The object of this work is to evaluate the new airborne infrared techniques and equipment as a complement to the data obtained from panchromatic aerial photography. This project is part of the Mexican remote sensing program of natural resources carried out under the auspices of the Comision Nacional del Espacio Exterior and in which the Research Institute (Instituto de Investigaciones de la Industria Electrica) is actively participating. The present study was made cooperatively with the U.S. National Aeronautics and Space Administration and the U.S. Geological Survey. The Los Negritos-Ixtlan de los Hervores geothermal fields are located east of Lake Chapala at the intersection of the Sierra Madre occidental and the west-central segment of the neovolcanic axis of Mexico. The two principal zones of hydrothermal activity occur in a tectonic trench filled with lake sediments of the Quaternary intercalated with Quaternary and Holocene volcanic rocks and characterized by an intricate system of block-fault tectonics, part of the Chapala-Acambay tectonic system, along which there has been volcanic activity in modern time. Surface manifestations of geothermal activity consist of relatively high heat flow and hot springs, small geysers and small steam vents aligned along an E-W axis at Ixtlan, possibly at the intersection of major fault trends and mud volcanoes and hot pools aligned NE-SW at Los Negritos. More than 20 exit points of thermal waters are shown on infrared imagery to be aligned along an extension of the Ixtlan fault between Ixtlan and El Salitre. A narrow zone of hydrothermal alteration and deposition at the surface is identifiable on the infrared imagery of this area, closey related spatially to a resistivity low at depth. Extinct geothermal areas near El Salitre, Ixtlan, and farther west at San Gregorio are clearly delineated on both infrared images and infrared ektachrome photographs. Predawn infrared images also show high-angle fault zones suggesting the dominance of block tectonics in much of the area. Special image enhancement techniques applied to the original magnetic tape records will be required for more precise identification of warm ground zones and for a qualitative or semiquantitative estimate of ground radiance associated with anomalously high convective heat flow.

Geothermics↗
USGS2012

Mapping outdoor recreationists' perceived social values for ecosystem services at Hinchinbrook Island National Park, Australia

Coastal ecosystems are increasingly faced with human impacts. To better understand these changing conditions, biophysical and economic values of nature have been used to prioritize spatial planning efforts and ecosystem-based management of human activities. Less is known, however, about how to characterize and represent non-material values in decision-making. We collected on-site and mailback survey data ( n = 209), and analyzed these data using the Social Values for Ecosystem Services (SolVES) GIS application to incorporate measures of social value and natural resource conditions on Hinchinbrook Island National Park, Australia. Our objectives in this paper are to: 1) determine the spatial distribution and point density of social values for ecosystem services; 2) examine the relationship between social values and natural resource conditions; and 3) compare social value allocations between two subgroups of outdoor recreationists. Results suggest that high priority areas exist on Hinchinbrook's land and seascapes according to the multiple values assigned to places by outdoor recreationists engaged in consumptive (e.g., fishing) and non-consumptive (e.g., hiking) activities. We examine statistically significant spatial clustering across two subgroups of the survey population for three value types that reflect Recreation, Biological Diversity, and Aesthetic qualities. The relationship between the relative importance of social values for ecosystem services and spatially-defined ecological data is explored to guide management decision-making in the context of an island national park setting.

Hinchinbrook Island National Park↗
USGS2015

Biogeochemical aspects of uranium mineralization, mining, milling, and remediation

Natural uranium (U) occurs as a mixture of three radioactive isotopes: 238 U, 235 U, and 234 U. Only 235 U is fissionable and makes up about 0.7% of natural U, while 238 U is overwhelmingly the most abundant at greater than 99% of the total mass of U. Prior to the 1940s, U was predominantly used as a coloring agent, and U-bearing ores were mined mainly for their radium (Ra) and/or vanadium (V) content; the bulk of the U was discarded with the tailings (Finch et al., 1972). Once nuclear fission was discovered, the economic importance of U increased greatly. The mining and milling of U-bearing ores is the first step in the nuclear fuel cycle, and the contact of residual waste with natural water is a potential source of contamination of U and associated elements to the environment. Uranium is mined by three basic methods: surface (open pit), underground, and solution mining (in situ leaching or in situ recovery), depending on the deposit grade, size, location, geology and economic considerations (Abdelouas, 2006). Solid wastes at U mill tailings (UMT) sites can include both standard tailings (i.e., leached ore rock residues) and solids generated on site by waste treatment processes. The latter can include sludge or “mud” from neutralization of acidic mine/mill effluents, containing Fe and a range of coprecipitated constituents, or barium sulfate precipitates that selectively remove Ra (e.g., Carvalho et al., 2007). In this chapter, we review the hydrometallurgical processes by which U is extracted from ore, the biogeochemical processes that can affect the fate and transport of U and associated elements in the environment, and possible remediation strategies for site closure and aquifer restoration. This paper represents the fourth in a series of review papers from the U.S. Geological Survey (USGS) on geochemical aspects of UMT management that span more than three decades. The first paper (Landa, 1980) in this series is a primer on the nature of tailings and radionuclide mobilization from them. The second paper (Landa, 1999) includes coverage of research carried out under the U.S. Department of Energy’s Uranium Mill Tailings Remedial Action Program (UMTRA). The third paper (Landa, 2004) reflects the increased focus of researchers on biotic effects in UMT environs. This paper expands the focus to U mining, milling, and remedial actions, and includes extensive coverage of the increasingly important alkaline in situ recovery and groundwater restoration.

Applied Geochemistry↗
USGS2015

On formally integrating science and policy: walking the walk

The contribution of science to the development and implementation of policy is typically neither direct nor transparent. In 1995, the U.S. Fish and Wildlife Service (FWS) made a decision that was unprecedented in natural resource management, turning to an unused and unproven decision process to carry out trust responsibilities mandated by an international treaty. The decision process was adopted for the establishment of annual sport hunting regulations for the most economically important duck population in North America, the 6 to 11 million mallards Anas platyrhynchos breeding in the mid-continent region of north-central United States and central Canada. The key idea underlying the adopted decision process was to formally embed within it a scientific process designed to reduce uncertainty (learn) and thus make better decisions in the future. The scientific process entails use of models to develop predictions of competing hypotheses about system response to the selected action at each decision point. These prediction not only are used to select the optimal management action, but also are compared with the subsequent estimates of system state variables, providing evidence for modifying degrees of confidence in, and hence relative influence of, these models at the next decision point. Science and learning in one step are formally and directly incorporated into the next decision, contrasting with the usual ad hoc and indirect use of scientific results in policy development and decision-making. Application of this approach over the last 20 years has led to a substantial reduction in uncertainty, as well as to an increase in transparency and defensibility of annual decisions and a decrease in the contentiousness of the decision process. As resource managers are faced with increased uncertainty associated with various components of global change, this approach provides a roadmap for the future scientific management of natural resources.

Journal of Applied Ecology↗
USGS2016

Hydropower assessment of Bolivia—A multisource satellite data and hydrologic modeling approach

This study produced a geospatial database for use in a decision support system by the Bolivian authorities to investigate further development and investment potentials in sustainable hydropower in Bolivia. The study assessed theoretical hydropower of all 1-kilometer (km) stream segments in the country using multisource satellite data and a hydrologic modeling approach. With the assessment covering the 2 million square kilometer (km 2 ) region influencing Bolivia’s drainage network, the potential hydropower figures are based on theoretical yield assuming that the systems generating the power are 100 percent efficient. There are several factors to consider when determining the real-world or technical power potential of a hydropower system, and these factors can vary depending on local conditions. Since this assessment covers a large area, it was necessary to reduce these variables to the two that can be modeled consistently throughout the region, streamflow or discharge, and elevation drop or head. First, the Shuttle Radar Topography Mission high-resolution 30-meter (m) digital elevation model was used to identify stream segments with greater than 10 km 2 of upstream drainage. We applied several preconditioning processes to the 30-m digital elevation model to reduce errors and improve the accuracy of stream delineation and head height estimation. A total of 316,500 1-km stream segments were identified and used in this study to assess the total theoretical hydropower potential of Bolivia. Precipitation observations from a total of 463 stations obtained from the Bolivian Servicio Nacional de Meteorología e Hidrología (Bolivian National Meteorology and Hydrology Service) and the Brazilian Agência Nacional de Águas (Brazilian National Water Agency) were used to validate six different gridded precipitation estimates for Bolivia obtained from various sources. Validation results indicated that gridded precipitation estimates from the Tropical Rainfall Measuring Mission (TRMM) reanalysis product (3B43) had the highest accuracies. The coarse-resolution (25-km) TRMM data were disaggregated to 5-km pixels using climatology information obtained from the Climate Hazards Group Infrared Precipitation with Stations dataset. About a 17-percent bias was observed in the disaggregated TRMM estimates, which was corrected using the station observations. The bias-corrected, disaggregated TRMM precipitation estimate was used to compute stream discharge using a regionalization approach. In regionalization approach, required homogeneous regions for Bolivia were derived from precipitation patterns and topographic characteristics using a k -means clustering approach. Using the discharge and head height estimates for each 1-km stream segment, we computed hydropower potential for 316,490 stream segments within Bolivia and that share borders with Bolivia. The total theoretical hydropower potential (TTHP) of these stream segments was found to be 212 gigawatts (GW). Out of this total, 77.4 GW was within protected areas where hydropower projects cannot be developed; hence, the remaining total theoretical hydropower in Bolivia (outside the protected areas) was estimated as 135 GW. Nearly 1,000 1-km stream segments, however, were within the boundaries of existing hydropower projects. The TTHP of these stream segments was nearly 1.4 GW, so the residual TTHP of the streams in Bolivia was estimated as 133 GW. Care should be exercised to understand and interpret the TTHP identified in this study because all the stream segments identified and assessed in this study cannot be harnessed to their full capacity; furthermore, factors such as required environmental flows, efficiency, economics, and feasibility need to be considered to better identify a more real-world hydropower potential. If environmental flow requirements of 20–40 percent are considered, the total theoretical power available reduces by 60–80 percent. In addition, a 0.72 efficiency factor further reduces the estimation by another 28 percent. This study provides the base theoretical hydropower potential for Bolivia, the next step is to identify optimal hydropower plant locations and factor in the principles to appraise a real-world power potential in Bolivia.

Open-File Report↗
USGS2022

Assessing mineral supply concentration from different perspectives through a case study of zinc

Increasing demand for nonfuel mineral commodities has increased concerns regarding the reliability of their supplies. “Criticality” assessments over the past decade have attempted to capture this concern through a set of indicators, the most common of which quantifies the risk associated with market concentration by applying the Herfindahl-Hirschman Index (HHI) to the world production of a given commodity by country in a given year. Although this approach is useful, it inherently assumes that all of world production is available to the market and is thus potentially at risk. In this analysis, the HHI, as well as HHI weighted by country governance, is calculated for mined and refined zinc using the standard approach of using all world production data and comparing that to the HHI when using the best estimate of what is available to the world market. The results indicate that although the HHI of both mined and refined zinc world production has increased markedly over the past decade, the HHI for what is available to the market for mined and refined zinc has remained relatively constant and low, which is indicative of minimal supply risk. This is mainly owing to the fact that a large and increasing share of the world’s mined and refined zinc production comes from China, but that production supplies domestic consumption as well as small amounts of exports. As a result, the zinc materials that are available to the world market outside of China are produced by a relatively large and diverse set of countries. Although these analyses are specific to zinc, they are likely to be comparable for other commodities of which the largest producers are also the largest consumers and highlight the importance of examining different perspectives in criticality assessments.

Mineral Economics↗
USGS2016

Climate change implications for tropical islands: Interpolating and interpreting statistically downscaled GCM projections for management and planning

The potential ecological and economic effects of climate change for tropical islands were studied using output from 12 statistically downscaled general circulation models (GCMs) taking Puerto Rico as a test case. Two model selection/model averaging strategies were used: the average of all available GCMs and the average of the models that are able to reproduce the observed large-scale dynamics that control precipitation over the Caribbean. Five island-wide and multidecadal averages of daily precipitation and temperature were estimated by way of a climatology-informed interpolation of the site-specific downscaled climate model output. Annual cooling degree-days (CDD) were calculated as a proxy index for air-conditioning energy demand, and two measures of annual no-rainfall days were used as drought indices. Holdridge life zone classification was used to map the possible ecological effects of climate change. Precipitation is predicted to decline in both model ensembles, but the decrease was more severe in the “regionally consistent” models. The precipitation declines cause gradual and linear increases in drought intensity and extremes. The warming from the 1960–90 period to the 2071–99 period was 4.6°–9°C depending on the global emission scenarios and location. This warming may cause increases in CDD, and consequently increasing energy demands. Life zones may shift from wetter to drier zones with the possibility of losing most, if not all, of the subtropical rain forests and extinction risks to rain forest specialists or obligates.

Puerto Rico↗
USGS1964

Preliminary results of hydrogeologic investigations in the valley of the Humboldt River near Winnemucca, Nevada

Most of the ground water of economic importance and nearly all the ground water closely associated with the flow o# the Humboldt River in the. 40-mile reach near Winnemucca, Nev., are in unconsolidated sedimentary deposits. These deposits range in age from Pliocene to Recent and range in character from coarse poorly sorted fanglomerate to lacustrine strata of clay, silt, sand, and gravel. The most permeable deposit consists of sand and gravel of Lake Lahontan age--the so-called medial gravel unit--which is underlain and overlain by fairly impermeable silt and clay also of Lake Lahontan age. The ultimate source of nearly all the water in the study area is precpitation within the drainage basin of the Humboldt River. Much of this water reaches the study, area as flow or underflow of the Humboldt River and as underflow from other valleys tributary to the study area. Little if any flow from the tributary streams in the study area usually reaches the Humboldt River. Most of the tributary streamflow within the study area evaporates or is transpired by vegetation, but a part percolates downward through unconsolidated deposits of the alluvial fans flanking the mountains and move downgradient as ground-water underflow toward the Humboldt River. Areas that contribute significant amounts of ground-water underflow to. the valley of the Humboldt River within the study area are (1) the valley of the Humboldt River upstream from the study area, (2) the Pole Creek-Rock Creek area, (3) Paradise Valley, and (4) Grass Valley and the northwestern slope of the Sonoma Range. The total average underflow from these areas in the period 1949-61 was about 14,000-19,000 acre-feet per year. Much of this underflow discharged into the Humboldt River within the study area and constituted a large part of the base flow of the river. Streamflow in the Humboldt River increases substantially in the early spring, principally because of runoff to the river in the reaches upstream from the study area. The resulting increase of the stage of the river causes the river to lose large amounts of water by infiltration to the ground-water reservoir in the study area. In addition, there is much recharge to the ground-water reservoir in the spring and early summer as a result of seepage losses from irrigation ditches and the downward percolation of some of the excess water applied for irrigation. The average net increase of ground water in storage in the deposits beneath and adjacent to the flood plain of the Humboldt River during the spring and early summer is about 10,000 acre-feet.

Nevada↗
USGS1967

Silver-bearing black calcite in western mining districts

The name black calcite has been applied from time to time to a dark gray to black variety of calcite or aragonite in metalliferous ore deposits in the Western States. Most of the material shows curved, roughly rhombic cleavage. The color is due to dispersed minute grains of one or more oxides of manganese, so dispersed that they rarely show any relation to cleavage surfaces or outward crystal forms.Recent study shows that the black calcite from numerous metalliferous deposits in the Western States contains small amounts of several metals, especially zinc, lead, and silver; in places, the crude black calcite contains as much as one percent silver and the black residue remaining after solution of the calcite as much as 1,500 ounces of silver to the ton. Black calcite is recognized in some deposits of Late Cretaceous to early Tertiary age where the host is Paleozoic carbonate rocks; thus far, this variety contains the larger amounts of silver. It is also present in some vein deposits where the host is generally volcanic rocks of middle Tertiary age; these are known in many districts in the southwest and even though the lead and zinc content is noteworthy, the silver content rarely exceeds one ounce to the ton (8).By studying polished sections of black calcite from the Aurora mine, White Pine County, Nevada, three silver-bearing manganate minerals have been recognized; (1) an argentian chalcophanite containing 7.5 percent silver (Ag 2 O), to which, in a following paper, the name aurorite is given; (2) argentian todorokite containing 3.9 percent silver (Ag 2 O), and (3) hydrous silver-bearing lead manganese oxide, containing 1.18 percent silver (Ag 2 O).

Nevada↗
USGS1992

Applying the DRASTIC model: A review of county-scale maps

The potential for contamination of ground water has become an issue of great concern to citizens and government alike, especially in the last decade. Numerous methods for assessing ground water vulnerability to contamination have been proposed and used. These methods include field-scale deterministic models that predict the rate of migration and fate of specific chemicals, and regional models that attempt to show general trends of ground-water vulnerability to contamination. To address regional ground-water vulnerability, a model may focus on characteristics of the soil, the geologic materials at and above the unconfined water table, or on the larger part of the geologic framework containing aquifers. A model's focus significantly affects the characteristics and, therefore, the utility of maps produced by the model. For example, a model focusing on the soil or vadose zone cannot adequately characterize the contamination potential of confined aquifers. The most commonly-used method for regional assessment of ground-water vulnerability is called "DRASTIC," an acronym for the seven factors that comprise the model. The U.S. Environmental Protection Agency (EPA) funded the development of the DRASTIC model, and used it as an assessment tool in their National Pesticide Survey (USEPA, 1992). In that Survey and in other assessments of contamination potential from application of agricultural chemicals, the pesticide version of DRASTIC is used; it differs from the standard DRASTIC model only in the degree of weighting applied to the seven factors. The DRASTIC model was developed by a committee of technical advisors who used the consensus approach to specify the relative significance of each factor. DRASTIC was designed primarily as a regional tool for prioritization, or screening, to indicate those areas (of at least 100 acres in size) which are generally more sensitive to contamination and, therefore, in need of more detailed mapping and evaluation or monitoring. The committee also intended DRASTIC to support decisions on allocation of scarce monitoring or remediation resources, and to serve as an educational tool. A detailed user's manual for the DRASTIC model (Aller and others, 1987) provides descriptions of common hydrogeologic settings across the United States and expected values for the seven factors; with this information and a specified "weight" or multiplier applied to each factor, DRASTIC scores can be calculated for each setting, and a map of these DRASTIC scores generated. The authors of this model intended it to be sufficiently objective and straightforward that it could be effectively used by persons with but a rudimentary knowledge of hydrogeologic principles. By adjusting the factor values for actual or interpreted local conditions a more knowledgeable person should be able to produce a somewhat more realistic and, therefore, reliable map of ground water vulnerability than by using the expected values for a given hydrogeologic setting. Other models for evaluating regional ground-water vulnerability have been developed; of importance to this report are models developed by the Illinois State Geological Survey (ISGS), the Wisconsin Department of Natural Resources (WDNR), and a joint effort between the U.S. Geological Survey (USGS) and the ISGS. All of these models, as well as DRASTIC, can only estimate the relative potential for contamination, either with a numerical scheme or a hierarchy of contamination potential map units. The model developed by WDNR (Schmidt, 1987) is based on a numerical scheme. It was used to generate a statewide contamination susceptibility map of Wisconsin (Schmidt and Kessler, 1987); Schmidt (1987) defines ground water contamination susceptibility as the "ease with which water (and, presumably, accompanying contaminants) at the surface can reach the water table." The model is not limited to aquifers, but rather considers all unconfined ground water whether in sandy surface aquifers or low-permeability glacial till. Ground water in confined aquifers is not addressed. Contamination susceptibility is estimated by a factor-weighting scheme similar in concept to DRASTIC. The factors (soil texture, surficial deposits, depth to the water table, bedrock type, and glacial drift thickness) are assigned an arbitrary value based on perceived importance. These numbers are then weighted and summed to produce the contamination susceptibility score. This model is compared and contrasted with other models elsewhere in this report. In contrast, the ISGS method avoids the use of a numerical system to rate the relative potential for contamination, and instead orders the map units in a hierarchy from relatively low to relatively high contamination potential. Also in contrast to the two models described above, the ISGS model relies solely on the textural character of the geologic framework to a specified depth. For example, a statewide map (Berg and others, 1984, scale 1:500,000) addresses land burial of wastes, and shows the contamination potential of aquifers within geologic units in the upper 50 feet. Another ISGS map (Keefer and Berg, 1990) addresses the contamination of major, economically important aquifers, and therefore evaluates contamination potential to greater depths (to greater than 300 feet). In general, the ISGS contamination potential maps evaluate aquifers, and do not evaluate potential for contamination of ground water at the unconfined water table in geologic units that are not aquifers (for example, in glacial till or finegrained lake sediments). With the cooperation of the ISGS, I have developed a model (Seller and Berg, in press) for the regional assessment of aquifer contamination potential that is based on ISGS techniques, adapted to a broader map area where detailed information may be unavailable. This model was used to generate a map of aquifer contamination potential for an area encompassing parts of five states near southern Lake Michigan and Lake Erie; the model and map are currently being refined and evaluated.

Open-File Report↗
USGS2021

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↗
USGS2016

A study of the 2015 Mw 8.3 Illapel earthquake and tsunami: Numerical and analytical approaches

The September 16, 2015 Illapel, Chile earthquake triggered a large tsunami, causing both economic losses and fatalities. To study the coastal effects of this earthquake, and to understand how such hazards might be accurately modeled in the future, different finite fault models of the Illapel rupture are used to define the initial condition for tsunami simulation. The numerical code Non-hydrostatic Evolution of Ocean WAVEs (NEOWAVE) is employed to model the tsunami evolution through the Pacific Ocean. Because only a short time is available for emergency response, and since the earthquake and tsunami sources are close to the coast, gaining a rapid understanding of the near-field run-up behavior is highly relevant to Chile. Therefore, an analytical solution of the 2 ? 1 D shallow water wave equations is considered. With this solution, we show that we can quickly estimate the run-up distribution along the coastline, to first order. After the earthquake and tsunami, field observations were measured in the surrounded coastal region, where the tsunami resulted in significant run-up. First, we compare the analytical and numerical solutions to test the accuracy of the analytical approach and the field observations, implying the analytic approach can accurately model tsunami run-up after an earthquake, without sacrificing the time necessary for a full numerical inversion. Then, we compare both with field run-up measurements. We observe the consistency between the two approaches. To complete the analysis, a tsunami source inversion is performed using run-up field measurements only. These inversion results are compared with seismic models, and are shown to capture the broad-scale details of those models, without the necessity of the detailed data sets they invert.

Illapel↗
USGS2010

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↗
USGS2018

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↗
USGS1988

Review of mechanisms, methods, and theory for determining recharge to shallow aquifers in North Dakota

Effective management of ground-water resources requires knowledge of all components of the water budget for the aquifer of interest. Efforts to simulate ground-water flow prior to development and the effects of proposed pumping in several of North Dakota's shallow glacial aquifers have been hindered by the lack of reliable estimates of ground-water recharge. This study was done to (1) review the methods that have been used to measure recharge, (2) review the theory of unsaturated flow and the methods for characterizing the physical properties of unsaturated media, (3) consider the relative merits of a rigorous data-intensive approach versus an estimation approach to the study of recharge, and (4) review past and current agronomic research in North Dakota for applicability of the research and the data generated to the study of recharge. Direct, quantitative techniques for evaluating recharge are rarely applied. The theory for computing fluxes in unsaturated media is well established and numerous physics-based models that effectively implement the theory are available, but the data required for the models generally are lacking. Many parametric approaches have been developed to avoid the large data requirements of the physics-based approaches for analyzing flow in the unsaturated zone. However, the parametric approaches normally include fitting coefficients that must be calibrated for every study site, thereby detracting from the general utility of the parametric approach. The functional relation of matric potential to moisture content is required for physics-based soil-water models, whether analytic or numeric. Laboratory methods to determine these relations are tedious, costly, and may not give results representative of the soils as they occur in the field. Many models have been proposed to estimate the moisture-characteristic curve and hydraulic-conductivity function from basic soil properties, but none yield results that are universally satisfactory. In situ methods, because they require minimal disturbance of the soil profile and may be used repeatedly on the same soil mass, have become the preferred means for acquiring physical data, especially hydraulic conductivity. Hydro logic investigations, except for recent studies of hazardous-waste disposal sites, rarely have included physical characterizations of unsaturated media. Any of four phenomena could hinder attempts to simulate unsaturated flow in settings typical of North Dakota; variability of soil properties, hysteresis, frozen ground, and macropore development. The spatial and temporal variability of soil properties probably is the greatest complicating phenomenon and must be dealt with by detailed characterization of the properties. Hysteresis can detract from the accuracy of flow calculations for some soils under certain conditions but, for the present, our scant knowledge of soil physical properties is a greater hindrance to reliable soi1-water mode 1 ing than is the hysteresis phenomenon. A1 though seasona1ly frozen ground undoubtedly affects hydrologic processes in North Dakota, much more research is needed before meaningful quantitative treatment is possible. Finally, macropores can influence soil-water movement significantly, but macropore development may not be common on the intensively farmed, coarse-textured soils that typically overlie North Dakota's glacial aquifers. Lysimetry currently is the only reliable means of analyzing macropore flow. The soil-related research that has been conducted in North Dakota to date (1983) provides little of the type of information required to estimate ground-water recharge. Useful data could be developed by systematically evaluating the hydraulic characteristics of the prominent soil types overlying North Dakota's shallow glacial aquifers. These data would be required to enable use of a physics-based approach to estimating recharge. The size of the aquifer under study, its economic value, and the resources available for data collection should be considered when choosing between parametric or physics-based methods.

Water-Resources Investigations Report↗
USGS2022

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↗
USGS1958

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↗
USGS1985

Electrical geophysical investigations of massive sulfide deposits and their host rocks, West Shasta copper-zinc district

The West Shasta copper-zinc district, Shasta County, California, contains many volcanogenic sulfide deposits within Middle Devonian rhyolites that have not been highly metamorphosed. The district was selected by the U.S. Geological Survey for intensive geological, geochemical, and geophysical study under the Development of Assessment Techniques (DAT) project because accessible exposures have been created by erosion and mining. This report describes the geophysical methods applied to characterize the electrical properties of selected West Shasta massive sulfide deposits and their host rocks, at both small (less than 25 ft) and large (greater than 25 ft) scales. The electrical techniques used galvanic (spectral induced polarization--SIP) and induction (very low frequency--VLF, slingram, and time domain electro-magnetics--TDEM) methods.In situ spectral induced polarization measurements were carried out to determine whether or not conductive anomalies in the district could be differentiated by their polarization signatures. The sulfide, in situ, induced polarization-phase spectral signatures (the induced polarization effect as a function of frequency) have much less character and lack the distinctive shape reported for other massive sulfide deposits; however, they do have some identifiable massive sulfide traits, such as low resistivity and variable polarizability. The nondescript sulfide spectral signature is attributed to the poor development of polarization processes due to a high percentage of resistive, nonpolarizable gangue minerals, lack of pore space, and limited electrolytic fluids. Large-scale spectral induced polarization measurements over the Hornet orebody have a greater polarization than the in situ measurements. This observation, in addition to the fact that much of the Hornet sulfide body has been removed by previous mining activity, suggests that the dominant polarization processes occur at the ground-water-sulfide interface.Combined use of induction techniques, which have different depths of penetration, were used to locate conductive anomalies and determine their shape and depth. All the induction surveys over the Hornet orebody detected the conductive tabular-shaped massive pyritic sulfide deposit hosted in resistive rhyolite. Shallow penetrating induction methods near the Keystone mine detected a conductive fault zone where a block of shale has been downfaulted into volcanic rock. Integrated interpretation of deeper penetrating induction data over this conductive fault zone indicates that parts of the shale are also conductive, demonstrating that the integrated use of several induction methods provides better conductor definition than a single method.

Economic Geology↗