USGS Science⌕ Search

SEARCH · USGS Science

Results for “Explorer”

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

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

At least 1,459 records · Page 81Linked to original sources

Water quality and hydrology of the Silver River Watershed, Baraga County, Michigan, 2005-08

The Silver River Watershed comprises about 69 square miles and drains part of northeastern Baraga County, Michigan. For generations, tribal members of the Keweenaw Bay Indian Community have hunted and fished in the watershed. Tribal government and members of Keweenaw Bay Indian Community are concerned about the effect of any development within the watershed, which is rural, isolated, and lightly populated. For decades, the area has been explored for various minerals. Since 2004, several mineral-exploration firms have been actively investigating areas within the watershed; property acquisition, road construction, and subsurface drilling have taken place close to tributary streams of the Silver River. The U.S. Geological Survey, in cooperation with Keweenaw Bay Indian Community, conducted a multi-year water-resources investigation of the Silver River Watershed during 2005-08. Methods of investigation included analyses of streamflow, water-quality sampling, and ecology at eight discrete sites located throughout the watershed. In addition, three continuous-record streamgages located within the watershed provided stage, discharge, specific conductance, and water-temperature data on an hourly basis. Water quality of the Silver River Watershed is typical of many streams in undeveloped areas of Upper Michigan. Concentrations of most analytes typically were low, although several exceeded applicable surface-water-quality standards. Seven samples had concentrations of copper that exceeded the Michigan Department of Environmental Quality standards for wildlife, and one sample had concentrations of cyanide that exceeded the same standards. Concentrations of total mercury at all eight sampling sites exceeded the Great Lakes Basin water-quality standard, but the ratio of methylmercury to total mercury was similar to the 5 to 10 percent found in most natural waters. Concentrations of arsenic and chromium in bed sediments were near the threshold-effect concentration. A qualitative ecological assessment of fishes and macroinvertebrates showed that intolerant salmonids were present at most sampled sites, and macroinvertebrate communities were indicative of near-excellent or excellent conditions at all eight sites. This baseline information will aid in an ongoing monitoring effort designed to protect the water resources of the

Scientific Investigations Report↗

Arc-related porphyry molybdenum deposit model

This report provides a descriptive model for arc-related porphyry molybdenum deposits. Presented within are geological, geochemical, and mineralogical characteristics that differentiate this deposit type from porphyry copper and alkali-feldspar rhyolite-granite porphyry molybdenum deposits. The U.S. Geological Survey's effort to update existing mineral deposit models spurred this research, which is intended to supplement previously published models for this deposit type that help guide mineral-resource and mineral-environmental assessments. Arc-related porphyry molybdenum deposits are a substantial resource for molybdenum metal and may have anomalous concentrations of tungsten. The deposits contain low-grade ore (0.03-0.22 percent molybdenum) as molybdenite, but are large-tonnage, making them amenable to bulk mining open-pit techniques. The mineralizing system usually has fluorine contents of less than 0.1 percent. The cogenetic intrusion is a differentiated calc-alkaline granitoid, typically granodiorite to quartz monzonite in composition, with low rubidium and niobium, and moderate to high strontium concentrations. Metals and hydrothermal fluids are sourced from these intrusions, with an additional meteoric fluid component contributing to peripheral alteration but not adding more metal. The lithology of the surrounding country rocks is not important to the formation of these deposits, but a surrounding carbonate unit may be altered to skarn that contains economic mineralization. The creation of contact-metamorphosed hornfels adjacent to the intrusion is common. Formation of arc-related porphyry molybdenum deposits typically occurs within a continental arc environment related to arc-continent or continent-continent collision and subduction. Few deposits are found in an island arc setting. Most classified arc-related porphyry molybdenum deposits are located in the western cordillera of North America, notably in British Columbia and Alaska. Hydrothermal alteration provides a key component to the identification of a deposit. Alteration usually is zoned from a core of potassic plus/minus silicic alteration outwards through phyllic to propylitic alteration. Argillic alteration may be irregular in shape and will overprint earlier hydrothermal alteration. Exploration should be limited to magmatic arc belts that have been unroofed and eroded to levels of a few kilometers depth. Important geological vectors toward areas of higher grade mineralization include intensity of hydrothermal alteration, veining, and faulting. Anomalous levels of molybdenum, tungsten, copper, lead, or zinc in soils, tills, stream sediments, and drainage waters may indicate the presence of an arc-related porphyry molybdenum deposit. Geophysical exploration techniques have been met with minimal success because of the overall low concentration of associated sulfide and oxide minerals. Geoenvironmental concerns are generally low because of low volumes of sulfide minerals. Most deposits are marginally acid-generating to non-acid-generating with drainage waters being near-neutral pH because of the acid generating potential of pyrite being partially buffered by late-stage calcite-bearing veins. The low ore content results in a waste:ore ratio of nearly 1:1 and large tailings piles from the open-pit method of mining.

Scientific Investigations Report↗

Mineral-deposit model for lithium-cesium-tantalum pegmatites

Lithium-cesium-tantalum (LCT) pegmatites comprise a compositionally defined subset of granitic pegmatites. The major minerals are quartz, potassium feldspar, albite, and muscovite; typical accessory minerals include biotite, garnet, tourmaline, and apatite. The principal lithium ore minerals are spodumene, petalite, and lepidolite; cesium mostly comes from pollucite; and tantalum mostly comes from columbite-tantalite. Tin ore as cassiterite and beryllium ore as beryl also occur in LCT pegmatites, as do a number of gemstones and high-value museum specimens of rare minerals. Individual crystals in LCT pegmatites can be enormous: the largest spodumene was 14 meters long, the largest beryl was 18 meters long, and the largest potassium feldspar was 49 meters long. Lithium-cesium-tantalum pegmatites account for about one-fourth of the world’s lithium production, most of the tantalum production, and all of the cesium production. Giant deposits include Tanco in Canada, Greenbushes in Australia, and Bikita in Zimbabwe. The largest lithium pegmatite in the United States, at King’s Mountain, North Carolina, is no longer being mined although large reserves of lithium remain. Depending on size and attitude of the pegmatite, a variety of mining techniques are used, including artisanal surface mining, open-pit surface mining, small underground workings, and large underground operations using room-and-pillar design. In favorable circumstances, what would otherwise be gangue minerals (quartz, potassium feldspar, albite, and muscovite) can be mined along with lithium and (or) tantalum as coproducts. Most LCT pegmatites are hosted in metamorphosed supracrustal rocks in the upper greenschist to lower amphibolite facies. Lithium-cesium-tantalum pegmatite intrusions generally are emplaced late during orogeny, with emplacement being controlled by pre-existing structures. Typically, they crop out near evolved, peraluminous granites and leucogranites from which they are inferred to be derived by fractional crystallization. In cases where a parental granite pluton is not exposed, one is inferred to lie at depth. Lithium-cesium-tantalum LCT pegmatite melts are enriched in fluxing components including H2O, F, P, and B, which depress the solidus temperature, lower the density, and increase rates of ionic diffusion. This, in turn, enables pegmatites to form thin dikes and massive crystals despite having a felsic composition and temperatures that are significantly lower than ordinary granitic melts. Lithium-cesium-tantalum pegmatites crystallized at remarkably low temperatures (about 350–550 °C) in a remarkably short time (days to years). Lithium-cesium-tantalum pegmatites form in orogenic hinterlands as products of plate convergence. Most formed during collisional orogeny (for example, Kings Mountain district, North Carolina). Specific causes of LCT pegmatite-related magmatism could include: ordinary arc processes; over thickening of continental crust during collision or subduction; slab breakoff during or after collision; slab delamination before, during, or after collision; and late collisional extensional collapse and consequent decompression melting. Lithium-cesium-tantalum pegmatite deposits are present in all continents including Antarctica and in rocks spanning 3 billion years of Earth history. The global age distribution of LCT pegmatites is similar to those of common pegmatites, orogenic granites, and detrital zircons. Peak times of LCT pegmatite genesis at about 2640, 1800, 960, 485, and 310 Ma (million years before present) correspond to times of collisional orogeny and supercontinent assembly. Between these pulses were long intervals when few or no LCT pegmatites formed. These minima overlap with supercontinent tenures at ca. 2450–2225, 1625–1000, 875–725, and 250–200 Ma. Exploration and assessment for LCT pegmatites are guided by a number of observations. In frontier areas where exploration has been minimal at best, the key first-order criteria are an orogenic hinterland setting, appropriate regional metamorphic grades, and the presence of evolved granites and common granitic pegmatites. New LCT pegmatites are most likely to be found near known deposits. Pegmatites tend to show a regional mineralogical and geochemical zoning pattern with respect to the inferred parental granite, with the greatest enrichment in the more distal pegmatites. Mineral-chemical trends in common pegmatites that can point toward an evolved LCT pegmatite include: increasing rubidium in potassium feldspar, increasing lithium in white mica, increasing manganese in garnet, and increasing tantalum and manganese in columbite-tantalite. Most LCT pegmatite bodies show a distinctive internal zonation featuring four zones: border, wall, intermediate (where lithium, cesium, and tantalum are generally concentrated), and core. This zonation is expressed both in cross section and map view; thus, what may appear to be a common pegmatite may instead be the edge of a mineralized body. Neither lithium-cesium-tantalum pegmatites nor their parental granites are likely to cause serious environmental concerns. Soils and country rock surrounding a LCT pegmatite, as well as waste from mining operations, may be enriched in characteristic elements relative to global average soil and bedrock values. These elements may include lithium, cesium, tantalum, beryllium, boron, fluorine, phosphorus, manganese, gallium, rubidium, niobium, tin, and hafnium. Among this suite of elements, however, the only ones that might present a concern for environmental health are beryllium and fluorine, which are included in the U.S. Environmental Protection Agency drinking-water regulations with maximum contaminant levels of 4 micrograms per liter and 4 milligrams per liter, respectively.

Scientific Investigations Report↗

Descriptive models for epithermal gold-silver deposits

Epithermal gold-silver deposits are vein, stockwork, disseminated, and replacement deposits that are mined primarily for their gold and silver contents; some deposits also contain substantial resources of lead, zinc, copper, and (or) mercury. These deposits form in the uppermost parts of the crust, at depths less than about 1,500 meters below the water table, and at temperatures below about 300 °C. Most epithermal gold-silver deposits are genetically related to hydrothermal systems associated with subaerial volcanism and intrusion of calc-alkaline magmas along convergent plate margins. These deposits formed throughout most of geologic time, although most known deposits are Cenozoic, which reflects preferential preservation of these shallowly formed deposits in tectonically unstable regions. Epithermal gold-silver deposits range in size from tens of thousands to greater than 1 billion metric tons of ore and have gold contents of 0.1 to greater than 30 grams per metric ton and silver contents of less than 1 to several thousand grams per metric ton. Historically, these deposits have been an important source of gold and silver and are estimated to contain about 8 percent of global gold. The wide range of tonnage-grade characteristics makes epithermal gold-silver deposits an attractive target for small and large exploration and mining companies. This report constitutes a new descriptive model for epithermal gold-silver deposits. It summarizes characteristics of known deposits, including their geological, geophysical, geochemical, and geoenvironmental aspects. Models concerning the genesis of epithermal gold-silver deposits are discussed. The application of descriptive and genetic aspects of the model to mineral exploration and resource assessment of undiscovered deposits is described. Finally, areas where additional research is needed to better understand the genesis of these deposits are identified. An extensive summary table outlining the characteristics of about 100 epithermal gold-silver deposits is included as an appendix; this summary table includes most of the world’s largest epithermal gold-silver deposits, and many smaller, well-studied deposits.

Scientific Investigations Report↗

Geology and undiscovered resource assessment of the potash-bearing, Middle Devonian (Givetian), Prairie Evaporite, Elk Point Basin, Canada and United States

The U.S. Geological Survey (USGS) assessed undiscovered potash resources in the Elk Point Basin in Canada and the United States as part of a global mineral resource assessment. The Elk Point Basin is a large, Middle Devonian (Givetian) intracratonic evaporite basin covering approximately 1,200,000 square kilometers (km 2 ) and filled mainly with marine evaporite and minor clastic sedimentary rocks that contain stratabound potash-bearing salt. The potash-bearing salt is concentrated in four stratigraphic members (Patience Lake, Belle Plaine, White Bear, and Esterhazy) in the upper 100 meters (m) of the Prairie Evaporite and are separated by beds of halite (NaCl) that contain lesser—presently non-economic—amounts of sylvite (KCl) and carnallite (KMgCl 3 ·6H 2 O). The principal ore-bearing salt contains mainly sylvite. Four permissive tracts were defined that permit the presence of undiscovered stratabound potash (both sylvite- and carnallite-bearing salt) using geological criteria. Permissive tracts are defined by the spatial extent of each stratigraphic member that is at least 1 m thick, are less than 3 kilometers (km) from the surface, contain at least 4 percent equivalent potassium oxide (K 2 O), and contain the currently known resources. The permissive tracts include known potash deposits and potash occurrences as wells or mines not in production and show where undiscovered potash resources may be present. Well data are used to define the extent, thickness, average K 2 O equivalent grades, and volumes of each member. Data were supplied by the Saskatchewan Geological Survey or were obtained from published National Instrument (NI) 43-101 technical reports and other published reports, such as annual 10-K reports or news releases. The Elk Point Basin is the world’s largest source of potash, producing 23.0 million metric tons (Mt) of potassium chloride (KCl) (the equivalent of about 14.4 Mt of K 2 O) in 2018. In terms of global importance, the Elk Point Basin may contain 40 to greater than 50 percent of the world’s potash resources. Since 1962, potash companies have mined more than 1.5 trillion metric tons of ore containing 605 Mt of KCl (the equivalent of about 380 Mt of K 2 O). The total value of the ore produced through 2018 is on the order of $70 trillion (CAD). Potash is currently produced from eight conventional and three underground solution mines at depths ranging from 900 m to nearly 1,800 m. Estimates of the amount of potash in the Elk Point Basin vary considerably and the data and methods used in those estimations are not well documented. Known potash resources are approximately 99 billion metric tons (Bt) of ore containing 22 Bt of K 2 O equivalent. As a result of new mine openings and increased production capacity at existing mines, the total production capacity of mines in the Elk Point Basin has increased significantly (to about 32.8 Mt of KCl or 22.8 Mt of K 2 O equivalent per year). Additional production capacity of about 31 Mt of KCl (or 17 Mt of K 2 O equivalent) per year could be realized over the next decade if several current (as of 2019) exploration and development projects reach production status. Stratabound potash-bearing salt of the Prairie Evaporite presently underlies a total area of about 188,000 km 2 and has a total volume of about 2,690 cubic kilometers (km 3 ). Post-depositional solution processes considerably modified the mineralogy and presence of the potash-bearing salt. These changes had a profound effect on the volume and grade of potash resources that remained in the Prairie Evaporite and are a major consideration of exploration and mining operations as well as in this assessment of undiscovered potash resources. This USGS assessment includes the locations and possible amounts of undiscovered potash resources in the Prairie Evaporite. Volumes for each stratigraphic member were computed using member thicknesses and areal extent modified by actual, estimated geologic loss owing to salt dissolution and extraction ratios, as well as estimated distribution of carnallite and sylvite. Both sylvite- and carnallite-bearing salts were assessed for potash in this study. The assessment uses modern published grade and tonnage data. The amount of undiscovered potash is estimated by using Monte Carlo simulations to combine volume estimates of the potash-bearing members with probability distributions for average grade and bulk density. Mean potash grades (expressed as percentage of K 2 O equivalent) calculated using drill core analyses are 17.76 for the Patience Lake Member, 15.98 for the Belle Plaine Member, 10.66 for the White Bear Member, and 15.30 for the Esterhazy Member. Geologic losses reported as extraction ratios during mining may range from 27.5 to 41.6 percent and are dependent on mining method and local geologic conditions. The assessment determined that mean estimated undiscovered K 2 O equivalent resources for the Patience Lake, Belle Plaine, White Bear, and Esterhazy Members are 340, 220, 34, and 190 Bt, respectively, and estimated a total mean of 790 Bt for the entire Prairie Evaporite above a depth of 3 km. The total mineralized rock tonnage is estimated to be about 5,000 Bt. Most of the assessed potash is located within Saskatchewan with lesser amounts in Alberta and Manitoba as well as Montana and North Dakota within the United States. Although carnallite is mined for potash in Europe, it has historically been avoided in mining plans for potash-producing companies in Saskatchewan because of mining, processing, and grade considerations. Carnallite-rich salt is locally present in concentrations and volumes that could be a significant resource of magnesium chloride (MgCl 2 ) obtained as a byproduct of processing the carnallite for potash. Previously estimated reserves (not NI 43-101 compliant) of mineralized material from 1955 to 2019 are 695 Mt at 22.1 percent MgCl 2 . The total amount of K 2 O equivalent as carnallite was estimated during this USGS assessment to be about 120 Bt (or 180 Bt KCl). With uncertainties in defining the areal extent of carnallite in each of the potash-bearing members, the amount of MgCl 2 as carnallite in the Elk Point Basin could be approximately 180 Bt.

Elk Point basin↗

Economic filters for evaluating porphyry copper deposit resource assessments using grade-tonnage deposit models, with examples from the U.S. Geological Survey global mineral resource assessment: Chapter H in Global mineral resource assessment

An analysis of the amount and location of undiscovered mineral resources that are likely to be economically recoverable is important for assessing the long-term adequacy and availability of mineral supplies. This requires an economic evaluation of estimates of undiscovered resources generated by traditional resource assessments (Singer and Menzie, 2010). In this study, simplified engineering cost models were used to estimate the economic fraction of resources contained in undiscovered porphyry copper deposits, predicted in a global assessment of copper resources. The cost models of Camm (1991) were updated with a cost index to reflect increases in mining and milling costs since 1989. The updated cost models were used to perform an economic analysis of undiscovered resources estimated in porphyry copper deposits in six tracts located in North America. The assessment estimated undiscovered porphyry copper deposits within 1 kilometer of the land surface in three depth intervals. The results of the updated engineering cost model analysis for open-pit porphyry copper deposits are in agreement with the grade-tonnage boundary defining positive economic returns for copper deposits developed between 1989 and 2008. This correspondence demonstrates that the updated engineering cost equations are performing well and appear to be appropriate to evaluate the economic status of open-pit porphyry copper mines under current, and potentially future, economic conditions. Economic filters based on these simplified engineering cost models provide a method for estimating potential tonnages of undiscovered metals that may be economic in individual assessment areas. One implication of the economic filter results for undiscovered copper resources is that global copper supply will continue to be dominated by production from a small number of giant deposits. This domination of resource supply by a small number of producers may increase in the future, because an increasing proportion of new deposit discoveries are likely to occur in remote areas and be concealed deep beneath covering rock and sediments. Extensive mineral exploration activity will be required to meet future resource demand, because these deposits will be harder to find and more costly to mine than near-surface deposits located in more accessible areas. Relatively few of the new deposit discoveries in these high-cost settings will have sufficient tonnage and grade characteristics to assure positive economic returns on development and exploration costs.

Scientific Investigations Report↗

Porphyry copper assessment of Europe, exclusive of the Fennoscandian Shield: Chapter K in Global mineral resource assessment

The U.S. Geological Survey (USGS) collaborated with European geologists to assess resources in porphyry copper deposits in Europe, exclusive of Scandinavia (Sweden, Denmark, Norway, and Finland) and Russia. Porphyry copper deposits in Europe are Paleozoic and Late Cretaceous to Miocene in age. A number of the 31 known Phanerozoic deposits contain more than 1 million metric tons of contained copper, including the Majdanpek deposit, Serbia; Assarel, Bulgaria; Skouries, Greece; and Rosia Poeni, Romania. Five geographic areas were delineated as permissive tracts for post-Paleozoic porphyry copper deposits. Two additional tracts were delineated to show the extent of permissive igneous rocks associated with porphyry copper mineralization related to the Paleozoic Caledonian and Variscan orogenies. The tracts are based on mapped and inferred subsurface distributions of igneous rocks of specific age ranges that define areas where the occurrence of porphyry copper deposits within 1 kilometer of the Earth’s surface is possible. These tracts range in area from about 4,000 to 93,000 square kilometers. Although maps at a variety of different scales were used in the assessment, the final tract boundaries are intended for use at a scale of 1:1,000,000. The post-Paleozoic deposits in Europe all formed in conjunction with the tectonic evolution of southern Europe as the former Tethyan Ocean closed by convergence of the African and Arabian Plates with Europe, accompanied by accretion of microcontinents to the southern Eurasian Plate and development and demise of magmatic arcs and ocean basins. Many of the deposits formed in extensional or post-collisional settings; these tectonic environments are increasingly being recognized as environments where porphyry copper deposits occur. Probabilistic estimates of undiscovered porphyry copper deposits were made for four Phanerozoic permissive tracts; the other tracts are discussed qualitatively. Assessment participants estimated numbers of undiscovered deposits at different levels of confidence for the four tracts. These estimates were then combined with grade and tonnage models using Monte Carlo simulation to generate probabilistic estimates of amounts of in-place undiscovered resources. Additional resources that may be present in extensions of known deposits were not evaluated. Assessment results are reported in tables and graphs as expected amounts of metal and rock in undiscovered deposits at different quantile levels, as well as the arithmetic mean for each commodity for each tract. This assessment estimated a mean of 14 undiscovered porphyry copper deposits within the four permissive tracts for which estimates were made. On the basis of global grade and tonnage models, mean (arithmetic) estimated resources that could be associated with undiscovered deposits are about 46 million metric tons of copper and about 2,600 metric tons of gold, as well as byproduct molybdenum and silver. Reliable reported identified resources for the 27 deposits in the assessed areas total about 44 million metric tons of copper and about 2,300 metric tons of gold. Exploration for gold-rich porphyry systems is ongoing in some parts of historical copper mining districts in central Europe and in northwesternmost (European) Turkey. Political and social conflicts, environmental concerns associated with historical mining, and the global economic situation have had negative effects on exploration, development, and mining in Europe for many years. The assessment includes an overview with summary tables. Detailed descriptions of each tract, including the rationales for delineation and assessment, are given in appendixes A–G. Appendix H describes a geographic information system (GIS) that includes tract boundaries and point locations of known porphyry copper deposits and significant prospects.

Scientific Investigations Report↗

Porphyry copper assessment of East and Southeast Asia: Philippines, Taiwan (Republic of China), Republic of Korea (South Korea), and Japan

The U.S. Geological Survey collaborated with member countries of the Coordinating Committee for Geoscience Programmes in East and Southeast Asia (CCOP) on an assessment of the porphyry copper resources of East and Southeast Asia as part of a global mineral resource assessment. The assessment covers the Philippines in Southeast Asia, and the Republic of Korea (South Korea), Taiwan (Province of China), and Japan in East Asia. The Philippines host world class porphyry copper deposits, such as the Tampakan and Atlas deposits. No porphyry copper deposits have been discovered in the Republic of Korea (South Korea), Taiwan (Province of China), or Japan. Thirteen geographic areas were delineated as tracts that are permissive for porphyry copper deposits in the assessed area. Individual tracts range from about 3,000 to 100,000 square kilometers in area. Permissive tracts are delineated on the basis of mapped distributions of igneous rocks of specific age ranges that define subduction-related magmatic arcs or magmatic belts that might contain porphyry copper deposits. Most of these magmatic arcs are subduction related, although some porphyry deposits and prospects are present in back-arc or poorly understood tectonic settings. Maps at various scales were used in the compilation; however, the final tract boundaries are intended for use at a scale of 1:1,000,000. Numbers of undiscovered deposits were estimated at different levels of confidence for 10 permissive tracts in the Philippines including one area that extends to eastern Taiwan (Republic of China); permissive tracts in South Korea and Japan are discussed qualitatively. Estimates of numbers of undiscovered deposits were combined with grade and tonnage models using Monte Carlo simulation to estimate amounts of undiscovered resources. Grades and tonnages of known porphyry copper deposits in the study area were compared with global grade and tonnage models to determine the appropriate model for simulation of undiscovered resources. Most of the known deposits are best described as copper-gold subtypes of porphyry copper deposits. For some permissive tracts, a general porphyry copper-gold-molybdenum model was used. Thirty-eight porphyry copper deposits are known in the Philippines; the mean number of undiscovered deposits was estimated to be 28. Mean (arithmetic) resources that could be associated with the undiscovered deposits are 90 million metric tons of copper and 5,800 metric tons of gold, as well as byproduct molybdenum and silver. Additional resources that could be discovered in extensions to known deposits were not evaluated. Assessment results, presented in tables and graphs, indicate expected amounts of total contained metal and mineralized rock in undiscovered deposits at different quantile levels, as well as the arithmetic mean for each tract. The Philippines have a long history of porphyry exploration cycles and mine development, interrupted at times by political and social unrest, environmental concerns, and natural disasters. Changes in mining laws within the region and the recent high price of gold on the world market have prompted renewed interest in porphyry copper deposits in the region. South Korea and Japan have been thoroughly explored for many types of mineral deposits. Available data suggest that the permissive rocks in South Korea typically are too deeply eroded to preserve porphyry copper deposits. Porphyry copper systems may be present in Japan, but are likely to lie at depths greater than the 1 kilometer from the surface protocol adopted for this study. Descriptions of the geologic basis for delineating each tract, the data used, the geologic criteria and rationale for the assessment, and results of the assessment are included in appendixes along with the description of a geographic information system (GIS) that includes tract boundaries, known porphyry copper deposits and significant prospects, and assessment results.

Scientific Investigations Report↗

Platinum-group elements in southern Africa: mineral inventory and an assessment of undiscovered mineral resources

The platinum-group elements, platinum, palladium, rhodium, ruthenium, iridium, and osmium, possess unique physical and chemical characteristics that make them indispensable to modern technology and industry. However, mineral deposits that are the main sources of these elements occur only in three countries in the world, raising concerns about potential disruption in mineral supply. Using information in the public domain, mineral resource and reserve information has been compiled for mafic and ultramafic rocks in South Africa and Zimbabwe that host most of the world’s platinum-group element resources. As of 2012, exploration and mining companies have delineated more than 20 billion metric tons of mineralized rock containing 42,000 metric tons of platinum, 29,000 metric tons of palladium, and 5,200 metric tons of rhodium, primarily in mafic and ultramafic intrusions of the Bushveld Complex and the Great Dyke, in southern Africa. Additional mineralized rock is likely to occur in extensions to the well-explored and characterized volumes of mineralized rock. Underexplored extensions of stratabound platinum-group element (PGE) deposits in the Bushveld Complex in South Africa may contain 65,000 metric tons of platinum, palladium, and rhodium to a depth of 3 km. Rocks enriched in PGE, which occur near the contact of the Bushveld Complex with older Transvaal Supergroup sedimentary rocks, may contain 1,100 metric tons of platinum and 1,370 metric tons of palladium (mean estimate to a depth of 1 km). A stratabound platinum-group element deposit in the Great Dyke in Zimbabwe may contain 6,900 metric tons of undiscovered platinum, palladium, and rhodium. By comparison, the global net demand for PGE in 2012 was approximately 460 metric tons. Since the 1920s, mining has recovered 7,200 and 107 metric tons of platinum-group elements from the Bushveld Complex and the Great Dyke, respectively. The large layered intrusions in southern Africa—the Bushveld Complex and the Great Dyke—are now and will continue to be a major source of the world’s supply of PGE. Mining will not deplete the identified mineral resources and reserves or potential undiscovered mineral resources for many decades; however, in the near-term, PGE supply could be affected by social, environmental, political, and economic factors.

Southern Africa↗

Porphyry copper assessment of northeast Asia: Far East Russia and northeasternmost China: Chapter W in Global mineral resource assessment

The U.S. Geological Survey assesses resources (mineral, energy, water, environmental, and biologic) at regional, national, and global scales to provide science in support of land management and decision making. Mineral resource assessments provide a synthesis of available information about where mineral deposits are known and suspected to be in the Earth’s crust, which commodities may be present, and estimates of amounts of resources in undiscovered deposits. A probabilistic mineral resource assessment of undiscovered resources associated with porphyry copper deposits in northeast Asia—composed mainly of Far East Russia and a small part of northeasternmost China—was completed as part of a global mineral resource assessment. Porphyry copper deposits are the main source of copper globally. Russia is an important source of copper, consistently ranking as 6th, 7th, or 8th in world production since 2000, and ranked 7th in 2014. Most of this production has been from magmatic copper-nickel-platinum-group element, volcanogenic massive sulfide, and sediment-hosted copper deposit types. The purpose of the assessment was to (1) compile a database of known deposits and significant prospects, (2) delineate permissive areas (tracts) for undiscovered porphyry copper deposits that may be present in the upper kilometer of the Earth’s crust, and (3) provide probabilistic estimates of amounts of copper (Cu), molybdenum (Mo), gold (Au), and silver (Ag) that could be contained in undiscovered porphyry copper deposits in the tracts. The assessment was completed by the U.S. Geological Survey in collaboration with geologists from the Russian Academy of Sciences and industry consultants. The database of known deposits, significant prospects, and prospects includes an inventory of mineral resources in two known porphyry copper deposits, as well as key characteristics derived from available exploration reports for 70 significant porphyry copper prospects and 86 other prospects. Resource and exploration and development activity are updated with information current through February 2013.

Scientific Investigations Report↗

Modeling uncertainty in coal resource assessments, with an application to a central area of the Gillette coal field, Wyoming

Standards for the public disclosure of mineral resources and reserves do not require the use of any specific methodology when it comes to estimating the reliability of the resources. Unbeknownst to most intended recipients of resource appraisals, such freedom commonly results in subjective opinions or estimations based on suboptimal approaches, such as use of distance methods. This report presents the results of a study of the third of three coal deposits in which drilling density has been increased one order of magnitude in three stages. Applying geostatistical simulation, the densest dataset was used to check the results obtained by modeling the sparser drillings. We have come up with two summary displays of results based on the same simulations, which individually and combined provide a better assessment of uncertainty than traditional qualitative resource classifications: (a) a display of cell 90 percent confidence interval versus cumulative cell tonnage, and (b) a histogram of total resources. The first graph allows classification of data into any number of bins with dividers to be decided by the assessor on the basis of a discriminating variable that is statistically accepted as a measure of uncertainty, thereby improving the quality and flexibility of the modeling. The second display expands the scope of the modeling by providing a quantitative measure of uncertainty for total tonnage, which is a fundamental concern for stockholders, geologists, and decision makers. Our approach allows us to correctly model uncertainty issues not possible to predict with distance methods, such as (a) different levels of uncertainty for individual beds with the same pattern and density of drill holes, (b) different local degrees of reduction of uncertainty with drilling densification reflecting fluctuation in the complexity of the geology, (c) average reduction in uncertainty at a disproportionately lesser rate than the reduction in area per drill hole, (d) the proportional effect of higher uncertainty in areas of higher tonnages, despite a regular drilling pattern, (e) the possibility of a local increase in uncertainty despite drilling densification to reflect a more complex geology as the deposit is known in more detail, and (f) for exactly the same drilling pattern, tonnage per individual beds with different uncertainty than the aggregated tonnage. These results should be considered realistic improvements over distance methods used for quantitative classification of uncertainty in coal resource, such as U.S. Geological Survey Circular 891. 1 The approach should be a welcome addition to the toolkit of Competent Persons preparing public disclosures according to international mineral codes such as those promoted by the Combined Reserves International Reporting Standards Committee 2 and the Joint Ore Reserve Committee. 3 1 Wood, G.H., Jr., Kehn, T.M., Carter, M.D., and Culbertson, W.C., 1983, Coal resources classification system of the U.S. Geological Survey: U.S. Geological Survey Circular 891, 65 p. 2 CRIRSCO (Combined Reserves International Reporting Standards Committee), 2013, International reporting template for the public reporting of exploration results, mineral resources and mineral reserves: Accessed February 2014 at http://www.crirsco.com/crirsco_template_may2013.pdf. 3 JORC (Joint Ore Reserves Committee), 2012, Australasian code for reporting of exploration results, mineral resources and ore reserves: Accessed September 2014 at http://www.jorc.org/docs/jorc_code2012.pdf.

Wyoming↗

Application of a Weighted Regression Model for Reporting Nutrient and Sediment Concentrations, Fluxes, and Trends in Concentration and Flux for the Chesapeake Bay Nontidal Water-Quality Monitoring Network, Results Through Water Year 2012

In the Chesapeake Bay watershed, estimated fluxes of nutrients and sediment from the bay’s nontidal tributaries into the estuary are the foundation of decision making to meet reductions prescribed by the Chesapeake Bay Total Maximum Daily Load (TMDL) and are often the basis for refining scientific understanding of the watershed-scale processes that influence the delivery of these constituents to the bay. Two regression-based flux and trend estimation models, ESTIMATOR and Weighted Regressions on Time, Discharge, and Season (WRTDS), were compared using data from 80 watersheds in the Chesapeake Bay Nontidal Water-Quality Monitoring Network (CBNTN). The watersheds range in size from 62 to 70,189 square kilometers and record lengths range from 6 to 28 years. ESTIMATOR is a constant-parameter model that estimates trends only in concentration; WRTDS uses variable parameters estimated with weighted regression, and estimates trends in both concentration and flux. WRTDS had greater explanatory power than ESTIMATOR, with the greatest degree of improvement evident for records longer than 25 years (30 stations; improvement in median model R 2 = 0.06 for total nitrogen, 0.08 for total phosphorus, and 0.05 for sediment) and the least degree of improvement for records of less than 10 years, for which the two models performed nearly equally. Flux bias statistics were comparable or lower (more favorable) for WRTDS for any record length; for 30 stations with records longer than 25 years, the greatest degree of improvement was evident for sediment (decrease of 0.17 in median statistic) and total phosphorus (decrease of 0.05). The overall between-station pattern in concentration trend direction and magnitude for all constituents was roughly similar for both models. A detailed case study revealed that trends in concentration estimated by WRTDS can operationally be viewed as a less-constrained equivalent to trends in concentration estimated by ESTIMATOR. Estimates of annual mean flow-adjusted (ESTIMATOR) and flow-normalized (WRTDS) concentration for years initially constituting the end of a water-quality record showed a similar degree of variability as data for additional years were incrementally added and the initial estimates “aged.” On the basis of the results of this broad comparison of the two models, the U.S. Geological Survey is adopting WRTDS as the primary model for estimating constituent fluxes and trends throughout the CBNTN. Nutrient and sediment flux and trend estimates, based on WRTDS, are summarized narratively and tabulated in appendixes for all stations for which fluxes or trends were reported through water year 2012. WRTDS also was used to explore the sensitivity of flux and trend estimates to three data-quality issues common in many large-scale monitoring networks and evident in some of the CBNTN records. The potential effects of inconsistency in annual sampling effort and inconsistency in storm sampling effort were explored by way of a subsampling experiment using eight of the most densely sampled long-term (1985–2012) stations in the CBNTN as baseline datasets. From each dataset, a set of 10 “design guideline” subsamples was selected, consisting of 12 monthly samples and 8 targeted storm samples per year. The selection was conducted in a manner that preserved the overall intensity of storm sampling in the baseline data. These 10 subsamples were further manipulated to create “heterogeneous” subsamples by removing storm samples prior to 2003. The maximum relative difference between flow-normalized flux estimated in a single year from any of the 10 design guideline subsamples and values estimated in the corresponding year from baseline data was smallest for dissolved inorganic nitrogen (median of 8 stations = 6 percent of baseline estimate), but more appreciable for total phosphorus and sediment (medians of 22 and 32 percent, respectively). The maximum relative difference between flow-normalized flux estimated from from the 10 heterogeneous subsamples and values estimated in the corresponding year from baseline data was more pronounced, with medians for 8 stations of 15, 30, and 53 percent of the corresponding baseline estimates for dissolved inorganic nitrogen, total phosphorus, and sediment, respectively. The worst-case maximum relative differences between flow-normalize flux estimated in a single year from the 10 heterogeneous subsamples and values estimated in the corresponding year from baseline data were 25 percent for dissolved inorganic nitrogen, 37 percent for total phosphorus, and 250 percent for sediment. The results for the heterogeneous subsamples indicate that changes in storm sampling frequency can result in appreciable distortion of estimated trends in flow-normalized flux, especially for total phosphorus and sediment. Trend lines estimated from heterogeneous subsamples tended to converge with the trend lines estimated from baseline data after 2003. In contrast, 2003–12 trends based on subsamples truncated by discarding all data prior to the induced heterogeneity in 2003 showed appreciable biases and differences in slope, relative to the corresponding 2003–12 segment of the trend computed from the design guideline subsamples. Overall, the results indicate that for particulate constituents, load and trend estimates computed using long-term records recently converted to CBNTN design guideline sampling protocols will be most reliable if the trend is computed using the entire record, but reported only for the period that design guideline sampling protocols were followed. Inconsistencies related to changing laboratory methods were also examined via two manipulative experiments. In the first experiment, increasing and decreasing “stair-step” patterns of changes in censoring level, overall representing a factor-of-five change in the laboratory reporting limit, were artificially imposed on a 27-year record with no censoring and a period-of-record concentration trend of –68.4 percent. Trends estimated on the basis of the manipulated records were broadly similar to the original trend (–63.6 percent for decreasing censoring levels and –70.3 percent for increasing censoring levels), lending a degree of confidence that the survival regression routines upon which WRTDS is based are generally robust to data censoring. The second experiment considered an abrupt disappearance of low-concentration observations of total phosphorus, associated with a laboratory method change and not reflected through censoring, near the middle of a 28-year record. By process of elimination, an upward shift in the estimated flow-normalize concentration trend line around the same time was identified as a likely artifact resulting from the laboratory method change, although a contemporaneous change in watershed processes cannot be ruled out. Decisions as to how to treat records with potential sampling protocol or laboratory methods-related artifacts should be made on a case-by-case basis, and trend results should be appropriately qualified.

Maryland, New York, Pennsylvania, Virginia↗

Field trip guide to Mount St. Helens, Washington—Recent and ancient volcaniclastic processes and deposits

This field guide explores volcanic effusions, sediments, and landforms at Mount St. Helens in Washington. A detailed synopsis outlines the eruptive history of Mount St. Helens from about 300,000 years ago through 1980 and beyond. The five days in the field include about 28 stops and 12 potential stops. Exposures in valleys surrounding Mount St. Helens reveal records of diverse Pleistocene and Holocene processes including debris avalanche, lahar, huge water wave on a nearby lake, pyroclastic density currents (surge and flow), tephra fall, lava flow, the growth of domes, and Pleistocene glaciation. Many of the stops explore effects of the several catastrophes that constituted the 18 May 1980 eruption and made Mount St. Helens famous.

Washington↗

Potential effects of energy development on environmental resources of the Williston Basin in Montana, North Dakota, and South Dakota—Executive summary

Executive Summary The Williston Basin, which includes parts of Montana, North Dakota, and South Dakota in the United States and parts of Manitoba and Saskatchewan in Canada, has been explored as a potential source of energy resources since the early 20th century; however, commercially viable petroleum drilling and recovery began in earnest in the 1950s. When oil prices rose in the mid-1980s, the number of wells also increased and then subsequently declined. Interest in the Williston Basin increased again in the mid-2000s with the application of new drilling technology. Since then, development has increased rather quickly. Most of this new development has been facilitated by advances in horizontal drilling and hydraulic fracturing technologies. The North Dakota Department of Mineral Resources reported an increase of more than 10,000 producing wells between 2000 and the spring of 2016. In total, 84 percent of those 10,000 wells target the Bakken Formation, which is now home to one of the Nation’s largest energy booms. Current estimates suggest that exploration and drilling activities are expected to continue for the next 20 to 50 years; however, future activity will likely ebb and flow in response to energy prices. Although most energy has been developed on non-Federal property, more than 2,000 wells were started on federally managed lands in the three States that contain the Williston Basin between 2004 and 2015, though these numbers do not reflect whether or not these wells targeted the Bakken Formation. Executive Order no. 13604 (March 22, 2012) directs Federal agencies to improve the timeliness of the permitting process for extracting publically owned minerals, while minimizing negative environmental effects. This means that Federal agencies need information about how energy development may affect other resources they are tasked with managing. One example of where information about potential effects of development may be useful is the Bureau of Land Management’s permitting process. Permits may include stipulations or special conditions that limit unforeseen negative consequences or ameliorate potential conflicts of future development. Federal agencies also need to coordinate permitting actions to ensure that development complies with existing regulations (for example, the Endangered Species Act [16 U.S.C. § 1531 et seq.] or the National Environmental Protection Act [42 U.S.C. § 4321 et seq.]) without unnecessarily restricting or delaying development. Part of this coordination entails agreeing on the information that will be used to assess the potential effects of energy development, which should also improve efficiency of the permitting process. Within the Williston Basin, a group of Federal agencies called the Bakken Federal Executive Group is developing coordination strategies for numerous energy-related issues on Federal lands. This report was developed in cooperation with the Bureau of Land Management to provide them with the best available scientific information to support documentation of potential effects on resources that Federal agencies manage. This report summarizes information about the effects of energy development on air, water, and biological resources within the U.S. part of the Williston Basin. The topics discussed in the report were based on a prioritized list of information needs elicited from the Bakken Federal Executive Group. The list was developed using a process known as structured decision making or decision analysis. This process began with an initial scoping workshop to determine the range of decisions made by those involved directly in managing energy development and resources on public land. U.S. Geological Survey staff then developed a simple quantitative ranking tool to assess which information needs were of greatest importance to those decisions.

Montana, North Dakota, South Dakota↗

Bend-scale geomorphic classification and assessment of the Lower Missouri River from Sioux City, Iowa, to the Mississippi River for application to pallid sturgeon management

Management actions intended to increase growth and survival of pallid sturgeon ( Scaphirhynchus albus ) age-0 larvae on the Lower Missouri River require a comprehensive understanding of the geomorphic habitat template of the river. The study described here had two objectives relating to where channel-reconfiguration projects should be located to optimize effectiveness. The first objective was to develop a bend-scale (that is, at the scale of individual bends, defined as “cross-over to cross-over”) geomorphic classification of the Lower Missouri River to help in the design of monitoring and evaluation of such projects. The second objective was to explore whether geomorphic variables could provide insight into varying capacities of bends to intercept drifting larvae. The bend-scale classification was based on geomorphic and engineering variables for 257 bends from Sioux City, Iowa, to the confluence with the Mississippi River near St. Louis, Missouri. We used k-means clustering to identify groupings of bends that shared the same characteristics. Separate 3-, 4-, and 6-cluster classifications were developed and mapped. The three classifications are nested in a hierarchical structure. We also explored capacities of bends to intercept larvae through evaluation of linear models that predicted persistent sand area or catch per unit effort (CPUE) of age-0 sturgeon as a function of the same geomorphic variables used in the classification. All highly ranked models that predict persistent sand area contained mean channel width and standard deviation of channel width as significant variables. Some top-ranked models also included contributions of channel sinuosity and density of navigation structures. The sand-area prediction models have r-squared values of 0.648–0.674. In contrast, the highest-ranking CPUE models have r-squared values of 0.011–0.170, indicating much more uncertainty for the biological response variable. Whereas the persistent sand model documents that physical processes of transport and accumulation are systematic and predictable, the poor performance of the CPUE models indicate that additional processes will need to be considered to predict biological transport and accumulation.

Iowa, Missouri↗

Assessing water quality from highway runoff at selected sites in North Carolina with the Stochastic Empirical Loading and Dilution Model (SELDM)

In 2015, the U.S. Geological Survey (USGS) entered into a cooperative agreement with the North Carolina Department of Transportation (NCDOT) to develop a North Carolina-enhanced variation of the national Stochastic Empirical Loading and Dilution Model (SELDM) with available North Carolina-specific streamflow and water-quality data and to demonstrate use of the model by documenting selected simulation scenarios. The USGS developed the national SELDM in cooperation with the Federal Highway Administration (FHWA) to provide the tools and techniques necessary for performing stormwater-quality simulations. SELDM uses a stochastic mass-balance approach to estimate combinations of flows, concentrations, and loads of stormwater constituents from the site of interest (often a highway catchment; nonhighway areas, such as a large impervious area at a shopping center complex, also can be used) and the basin upstream from the stormwater outfall to assess the risk for adverse effects of runoff. SELDM also can be used to simulate the effectiveness of volume reduction, hydrograph extension, and water-quality concentration reductions by stormwater best management practices (BMPs), which are designed to help mitigate the effects of runoff on receiving water bodies. Some of the statistical inputs needed for the North Carolina-enhanced SELDM were either calculated or augmented using local or regional data from North Carolina. Streamflow statistics used by SELDM were determined for 266 streamgages across North Carolina on the basis of data available through the 2015 water year. Recession ratio statistics used for triangular hydrographs were also developed for 30 streamgages across the State. The NCDOT identified previous research reports on highway-runoff and BMP studies in North Carolina for review of potential data addition to the national FHWA Highway-Runoff Database (HRDB). Following USGS review of these data, a total of 25,087 event mean concentration values and 1,140 storm events for 39 highway-runoff sites and 195 analytes were uploaded to the national HRDB from six North Carolina highway-runoff research reports and a recent USGS bridge deck runoff study. Using data for 27 streamgages in North Carolina, a total of 57 water-quality transport curves were developed for seven constituents for use in simulating water-quality conditions in the upstream basin. Performance data for three BMPs (bioretention, grass strip or swale, and wetland channel) from NCDOT research data were incorporated into the North Carolina-enhanced SELDM for volume-reduction statistics, including the effectiveness of treating four water-quality constituents (total suspended solids, total nitrogen, total phosphorus, nitrate plus nitrite) and turbidity. Simulations using the North Carolina-enhanced SELDM are presented for two hypothetical upstream basins in the Piedmont ecoregion and one hypothetical highway site to demonstrate how simulations can be used to provide risk-based information about potential effects of stormwater runoff on downstream water quality and the potential for mitigating those risks by using BMPs. The first group of simulations explores the stochastic variability in dilution factors (the ratio of the highway runoff to the total downstream stormflow) for a hypothetical Piedmont rural creek having drainage areas ranging from 1 to 100 square miles. The second group of simulations examines dilution factors based on variations in precipitation, streamflow, and recession ratios for two hypothetical Piedmont upstream basins (rural and urban) where the drainage area was held constant at 25 square miles. These simulations indicate the sensitivity of results to variations in each of the three variables. The third group of simulations examines the effects of varied concentrations in the upstream basin on water-quality conditions downstream from the highway crossing. Variations in upstream water-quality conditions for three constituents (suspended sediment concentration, total nitrogen, and total phosphorus) are based on water-quality transport curves selected from among the 57 curves developed as part of this study to represent low-, medium-, and high-concentration statistics. Simulations completed for this third group also examine the potential effects of grass swale and bioretention BMP treatment on total nitrogen and total phosphorus concentrations in highway runoff. The BMP performance data from the NCDOT research reports were applied in this group of simulations. The stochastic mass-balance approach used in SELDM analyses and simulations provides a strong tool for engineers and water-resource managers to use in exploring a wide range of possible hydrologic and water-quality inputs and their effects on downstream water quality. The results of this study can not only aid engineers and managers in planning for potential adverse effects of runoff at site-specific locations, they can also help the USGS and other Federal and State agencies with oversight responsibilities in stormwater-quality issues to continue gathering data on potential water-quality effects in receiving streams.

North Carolina↗

Trends in streamflow and concentrations and flux of nutrients and total suspended solids in the Upper White River at Muncie, near Nora, and near Centerton, Indiana

The U.S. Geological Survey (USGS), in cooperation with The Nature Conservancy, completed a study to estimate and assess trends in streamflow and annual mean concentrations and flux of nutrients (nitrate plus nitrite, total Kjeldahl nitrogen, and total phosphorus) and total suspended solids at three USGS streamgages (hereafter referred to as “study gages”) on the Upper White River at Muncie (USGS station 03347000), near Nora (USGS station 03351000), and near Centerton (USGS station 03354000), Indiana. Water-quality data used in the analyses were collected by several agencies between calendar years 1991 and 2017, and streamflow (discharge) data were collected by the USGS. For most of the water-quality constituents, there were suitable data to facilitate an analysis of the 26-year period extending from calendar years 1991 to 2017 (water years 1992 to 2017); however, shorter analytical periods were necessary for total Kjeldahl nitrogen for the study gages at Muncie and near Centerton and for total suspended solids for the study gage near Centerton. Temporal trends in streamflows at the study gages for the period extending from water years 1978 to 2017 were assessed using Exploration and Graphics for RivEr Trends (EGRET) and Mann-Kendall and Pettitt tests. With just one exception, the annual maximum and mean daily streamflows and the annual minimum 7-day mean streamflows at the study gages demonstrated upward trends (increasing streamflows) in the EGRET analyses. The exception was the annual 7-day minimum streamflow at the study gage near Nora, which indicated no trend. Mann-Kendall tests also indicated that the average trend for the annual maximum daily, annual mean daily, and annual 7-day minimum streamflow statistics between water years 1978 and 2017 was upward at each of the study gages; however, only the trends in the annual mean daily streamflows at the study gage at Muncie and the annual maximum daily streamflows at the study gages near Nora and near Centerton were statistically significant at a 0.05 probability level. The Pettitt tests indicated that a statistically significant step trend (abrupt change) in annual mean daily streamflows occurred at each of the study gages around water year 2001. The seasonal distributions of total suspended solids, total phosphorus, nitrate plus nitrite, and total Kjeldahl nitrogen concentrations at the study gages were evaluated to identify patterns and other distinguishing characteristics by examining boxplots of concentrations as a function of month of the year. Seasonal distributions of nitrate plus nitrite concentrations and total suspended solids concentrations differed from each other but were generally similar among the three study gages for a given constituent. Median concentrations of nitrate plus nitrite were highest during the January–June months, whereas median concentrations of total suspended solids were highest during June and July. Seasonal distributions of total phosphorus concentrations were similar at the study gages near Nora and near Centerton, but the seasonal distribution was noticeably different at the study gage at Muncie, which had monthly median concentrations that were substantially lower than at the two downstream study gages (near Nora and near Centerton). The seasonal distribution of total Kjeldahl nitrogen concentrations differed in pattern among the three study gages; however, in general, some of the higher monthly median total Kjeldahl nitrogen concentrations at each study gage were associated with the late spring and summer periods. The Weighted Regressions on Time, Discharge, and Season (WRTDS) method implemented in EGRET was used to estimate water-year annual mean daily concentrations and flux of nutrients and total suspended solids, as well as estimates of concentrations and flux that were “normalized” to remove the effect of year-to-year variation in streamflow. The approximate coefficients of determination for the WRTDS regression models ranged from a high of 0.82 for total phosphorus for the study gage near Centerton to a low of 0.19 for nitrate plus nitrite for the study gage near Nora. Loads and yields of total suspended solids, total phosphorus, nitrate plus nitrite, and total Kjeldahl nitrogen were estimated for analytical periods consisting of the longest periods of concurrent record at the three study gages. Loads of each of the constituents increased sequentially from the most upstream study gage to the most downstream study gage; however, the same was not true for yields. The highest yields of total suspended solids, total phosphorus, and total Kjeldahl nitrogen occurred at the most upstream study gage (at Muncie); however, the highest yield of nitrate plus nitrite occurred at the most downstream study gage (near Centerton). WRTDS bootstrap tests were used to assess the magnitude, direction, and likelihood of changes in annual flow-normalized mean daily concentrations and flux of total suspended solids, total phosphorus, nitrate plus nitrite, and total Kjeldahl nitrogen at the study gages between water years 1997 and 2017. Changes in flow-normalized concentrations and flux of the constituents between water years 1997 and 2017 were mostly downward (decreasing). The exceptions were likely to highly likely upward (increasing) changes in (1) flow-normalized annual mean daily concentration and annual flux for total suspended solids and total phosphorus at the study gage at Muncie, (2) flow-normalized annual mean daily total phosphorus concentration at the study gage near Centerton, (3) flow-normalized annual flux of total phosphorus at the study gage near Centerton, and (4) flow-normalized annual mean daily nitrate plus nitrite concentration at the study gage near Centerton. Although an upward change in flow-normalized nitrate plus nitrite concentrations was likely at the study gage near Centerton, flow-normalized annual flux of nitrate plus nitrite at that study gage was determined to have a highly likely downward change. EGRET and Exploration and Graphics for RivEr Trends Confidence Intervals (EGRETci) analyses can be used to improve our understanding of how concentrations and flux change as functions of time and streamflow, as well as provide information on how the relations between streamflow and constituent concentrations have changed within the calendar year between any 2 years included in the analyses. Examples of those uses, illustrating changes between calendar years 1992 and 2017, were given for total suspended solids concentrations at the study gage near Nora and for nitrate plus nitrite concentrations at the study gage near Centerton.

Indiana↗

Grade and tonnage model for tungsten skarn deposits—2020 update

This report presents an updated grade and tonnage model for tungsten skarn deposits. As a critical component of the U.S. Geological Survey’s three-part form of quantitative mineral resource assessment, robust grade and tonnage models are essential to transforming mineral resource assessments into effective tools for decision makers. Using the best data available at the time of publication, this represents the first attempt in nearly 30 years to capture current mineral inventory and cumulative production data for worldwide tungsten skarn deposits. The accuracy of modern assessments of undiscovered tungsten skarn resources is highly influenced by the use of current data on the distribution of the grades and tonnages of well-explored tungsten skarn deposits. Primary factors affecting the changes to these distributions in the model presented here compared with those of previous models are the inclusion of important deposits, especially those in China that had been omitted in previous models; expanded mineral inventories resulting from increased exploration; and changes to international reporting standards. These factors have resulted in dramatic increases in average ore tonnage and slight decreases in the average grade of tungsten skarn deposits compared with previous models. Large increases in contained metal are observed among many of the individual deposits incorporated within this model that were also included in previous tungsten skarn grade and tonnage models. This report also provides recommendations for input parameters related to grade and tonnage models to use with software tools designed to facilitate the three-part form of quantitative mineral resource assessments.

Scientific Investigations Report↗