Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Mining Engineering”

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.

168 records · Page 4Linked to original sources
USGS2026

Future water constraints on United States lithium mining under climate change

Lithium is necessary for low-carbon technologies that combat climate change, but lithium extraction is water-intensive. Changes in temperature and precipitation arising from climate change are altering water distribution, which could further strain supplies for new mines and industry, farms, and households. Here we explored how climate change, water use, and mining siting could impact lithium mining in the United States. We analyzed whether there would be sufficient water available to support the single existing and 22 proposed U.S. lithium mines at mid-century under four socioeconomic-climate scenarios and five climate models. Though dependent on socioeconomic-climate scenario, climate model, and lithium deposit type, available water supply in most subbasins would likely be unable to support new mines’ water demands, or even non-mining water demands from other sectors. Water scarcity could hinder the ability of the United States to produce enough lithium to meet domestic demand thereby necessitating higher imports.

conterminous United States↗
USGS2025

Estimating the hypothetical endowment of critical minerals and other commodities in porphyry copper mine waste in the Four Corners states, USA

Society is fundamentally dependent upon commodities that are used in end-use products for the aerospace, defense, energy, telecommunication, and transportation sectors, resulting in centuries of mining to supply these commodities and materials. Waste from these mining operations can remain on the landscape indefinitely, but there is a lack of national understanding of the distribution and scale of such waste features. The renewable energy transition will continue to increase demand for critical minerals and will result in increasing volumes of mine waste on the Earth’s surface. Reprocessing mine waste can reduce environmental risks and recover needed commodities to match growing demand for societal growth. Therefore, understanding the approximate abundance of commodities that may be available for recovery within mine waste features can be an important piece of domestic critical mineral supply.

Arizona, Colorado, New Mexico, Utah↗
USGS2025

A synthesis engine for constructing geologic maps of the United States

The geologic history of the United States is cataloged in thousands of geologic maps produced during many decades. However, the disparate nature of these individual maps makes it challenging to assess resources, research geologic histories, or characterize natural hazards holistically across the Nation. The U.S. House of Representatives 2020 appropriations bill for the U.S. Department of the Interior (H.R. 116-100) requires the U.S. Geological Survey to “bring together detailed national and continental-resolution [two-dimensional] and [three-dimensional] information produced throughout the Survey and by [F]ederal and [S]tate partners.” In response to this directive, this report presents a compilation and synthesis of geologic maps across the United States in the form of a relational database. The synthesis database includes thematic maps that synthesize the Nation’s geology, and retains the original input maps as well as linkages to standardized vocabularies to aid the discoverability of geologic information. Specifically, the synthesis database is targeted toward producing four National-resolution maps for the conterminous United States: Quaternary geology, the geology at the Earth’s surface, pre-Quaternary geology, and Precambrian geology. In addition, the synthesis database includes the infrastructure necessary to expand to additional resolutions in the future.

conterminous United States↗
USGS2026

Remote sensing enables basin-scale inventories of coal mine methane

Underground coal mines are important global sources of methane, but emission estimates are uncertain. We show that emission estimates for individual mines from aircraft remote-sensing surveys in the United States agree within 40% with direct measurements used for national emission reporting (IPCC Tier 3 estimate). Such direct measurements are unavailable in most countries, which rely on estimated emission factors (EFs) applied to coal-production rates. We find that EFs from IPCC Tier 1 and the Model for Calculating Coal Mine Methane (MC2M) methods overestimate U.S. emissions 3-fold due to incorrect dependence on mine depth. An IPCC Tier 2 method using measured basin-specific mine gas content agrees with direct emission measurements but does not account for gob well emissions and requires gas content data that are generally unavailable. We show that aircraft remote sensing for a small sample of mines can successfully estimate basin-specific EFs for ventilation shafts and gob wells, enabling estimates of basin- and national-scale emissions. We discuss how the method can be applied with satellite remote sensing to quantify coal emissions worldwide.

Alabama, Colorado, Kentucky, New Mexico, Ohio, Pen↗
USGS2025

Mine waste as a potential source of critical minerals and other commodities: Examples from the Four Corners states, USA

The growing demand for critical minerals and other mineral resources has raised concerns about possible supplies of these essential commodities. Mine waste is a potential source of these essential commodities. We compiled a geospatial database of publicly available data of the largest non-fuel mine waste features (>200,000 m 2 in areal extent) in the Four Corners states of the United States, where most of those features are from porphyry Cu deposits. The combined volume of those large porphyry Cu mine waste features is approximately 17 cubic kilometers, 60% of which is ore-related material such as tailings. Using publicly available data on density, grade, and previous recovery values, we estimate the contained endowments of Ag, As, Au, Bi, Cu, Mo, Re, S, Sb, Se, Te, and W. These estimates indicate endowments within ore-related mine wastes are collectively comparable to endowments of giant to supergiant deposits. If fully recovered, these commodities could meet current global demand from less than 1 year (Sb) to more than 200 years (Re), underscoring the enormous untapped resource potential of mine waste.

Arizona, Colorado, New Mexico, Utah↗
USGS2025

Water resources related to breccia pipe uranium mining in the Grand Canyon region

Introduction In the arid Grand Canyon region, water resources are limited to primarily the Colorado River and associated tributaries and to groundwater in the form of seeps and springs. Groundwater resources in the region supply water for human use and support diverse and rich ecosystems in the locations immediately surrounding the seeps and springs. Throughout the region, uranium resources occur and may interact with water resources in both mined and unmined uranium deposits. There is a need to better understand groundwater in the region and the effects from uranium mining in order to better manage the limited water resources in the area. This Fact Sheet summarizes results from U.S. Geological Survey studies that were conducted on this topic from 2012 to 2023.

Arizona↗
USGS2026

Vegetation cover and composition in environments surrounding uranium mines in the Grand Canyon ecosystem, Northern Arizona

Mining uranium from breccia-pipe deposits in the greater Grand Canyon region has occurred since the mid-1900s. However, possible ecosystem contamination with harmful levels of radionuclides may have occurred due to mining activities in the 21st century. In response, a 20-year Federal moratorium on new mining claims in the Grand Canyon watershed was initiated in 2012, to allow time to evaluate the potential effects of uranium exploration and mining on human health, wildlife, and water resources. This moratorium, nor the 2023 designation of the “Baaj Nwaavjo I’tah Kukveni–Ancestral Footprints of the Grand Canyon National Monument,” precludes operation or development of mining claims predating 2012. Vegetation is a core ecosystem component that may be affected by uranium mining (for instance, through uptake and storage of radionuclides from the air or soil) or may act as a vector of exposure to wildlife, livestock, and humans (for instance, via their consumption of contaminated plant tissues). To provide baseline information about the plant communities associated with uranium mines in the Grand Canyon region, the U.S. Geological Survey surveyed an approximately 200-meter-wide buffer surrounding four breccia-pipe deposits, each in a unique stage of mine development, and at one reference area (a livestock water tank) that underwent ground disturbance but contains no mineral deposits. We sectioned the buffer zones into 0.65–4.52 hectare plots, within which we (1) inventoried all plant species, (2) measured percent cover of plant species, plant functional groups, and ground surface types (dark cyanobacteria, lichen, moss, bedrock, rock, embedded litter, duff, plant bases, and bare soil) using line-point intercept, and (3) measured length and frequency of gaps between perennial plant canopies using canopy gap intercept. We found that plant composition at the mines and the reference area differed from one another but were all characteristic of expected regional vegetation patterns. We provide this data summary as potential baseline information for future research and management efforts.

Arizona↗
USGS2026

Alternative approaches to dryland reclamation enhance vegetation cover and soil stability at a former uranium mine

Introduction Conventional reclamation methods to recover dryland ecosystems after mineral extraction often have low success. Alternative reclamation techniques may help overcome the many challenges to establishing persistent native vegetation and stable soils. Objectives Compare conventional reclamation approaches that focus primarily on soil stabilization to alternative approaches that may improve native ecosystem restoration at the Kanab North Mine in northern Arizona. Methods In a fully crossed experiment, conventional seed mix (native and non-native species) and seeding methods (drill-seeding) were compared to alternative practices, including a native-only seed mix, connectivity modifiers (ConMods) with hand seeding, and a biocrust inoculation. Restoration metrics like plant establishment, species composition, soil surface cover and stability, biocrust establishment, and near-surface wind-driven sediment were monitored over the first 5 years post-reclamation. Results While all restoration metrics improved with time, compared to traditional methods, ConMods increased total, graminoid, and native plant species cover, and the native-only seed mix increased woody, native, and perennial species cover. Biocrust inoculation increased late-successional biocrust development, but biocrusts only persisted over 5 years in combination with ConMods. Soil stability increased in all treatments, but both ConMods and biocrust inoculation increased stability compared to traditional practices. Conclusion Overall, alternative approaches to reclamation tested here (ConMods, native seed, and biocrust inoculation) were more effective than traditional practices at increasing total plant cover, facilitating biocrust establishment, and promoting soil stability after uranium mine reclamation.

Arizona↗
USGS2026

Distinguishing natural from mining-related metal sources by including streambank groundwater data in a stream mass loading study

Distinguishing stream metal loading caused by mine features from that caused by natural background sources remains challenging, yet this distinction is essential for making effective remedial decisions at many legacy mine sites. We combine a stream tracer injection and synoptic sampling study with data from shallow near-stream groundwater wells to estimate left-bank versus right-bank metal loading contributions at the 100-m spatial scale. The study was performed in the East Mancos River, a mountain headwater stream in Colorado, USA. The dominant source of elevated stream metal concentrations could be either groundwater infiltration through right-bank Doyle Mine waste piles or natural acid rock drainage from hydrothermally altered bedrock located mainly on the left bank. For the five metals of concern (Cu, Al, Zn, Cd, and Mn), we find that 15 % of the load contributed by diffuse groundwater inputs in the section potentially influenced by Doyle mine waste originates from the right bank. This right-bank potential mine contribution equates to only 3 % of the total watershed instream load for these metals. Furthermore, apparent 3 H/ 3 He groundwater ages in segments contributing most of the right-bank metal loading are sufficiently old (9–12 yr) to suggest that infiltration through the waste piles, located only 140–180 m from the stream, is unlikely. Estimated potential Doyle mine loading contributions can therefore be considered maximum values. Study results thus indicate that Doyle mine waste piles are a minor source of metal loading under low-flow conditions, and streambank groundwater data can provide valuable additional information in stream mass loading studies.

Colorado↗
USGS2025

Controls on water quality below a reclaimed surface coal mine, southeastern Montana

Coal mining and reclamation can have a profound influence on hydrogeologic systems, with clear consequences for groundwater quality, yet their long-term influence on downgradient water quality over time following reclamation is less well documented. Geochemical trends were evaluated in water quality downgradient of a fully reclaimed landscape at the former Big Sky Mine in the Rosebud Creek watershed (southeastern Montana, USA), over a 3-year period (2020–2022), including bond release in 2022. Within 6 km downgradient from the reclaimed area, sulfate concentrations decreased from approximately 3500 to 1800 mg l −1 within the Miller Coulee alluvial aquifer. Major ions, δ 34 S SO4 values, and residence time tracers suggest that the observed decreases in sulfate concentration result from a combination of dilution by mixed-age inflows and incomplete transit of the high salinity plume from the mine boundary. Both bedrock and alluvial aquifers of the Rosebud Creek corridor contained contributions of millennia-old regional groundwater, which may serve to mitigate mine-derived high salinity waters. Rosebud Creek, which traverses the outflow zone of Miller Coulee in the study area, exhibited high sulfate concentrations during low flows and consistent downgradient increases in sulfate concentration. The possibility of plume dynamics in Miller Coulee suggests that the greatest water quality impacts may not yet have reached Rosebud Creek.

Montana↗
USGS2025

Temporal and spatial comparison of coal mine ventilation methane emissions and mitigation quantified using PRISMA satellite data and on-site measurements

Emission monitoring at the facility level (bottom-up, BU) is key for accurate reporting of coal mine methane (CMM) emissions. Recent advances in aerial and satellite observations (top-down, TD) indicate that these methods have the potential to support CMM emissions monitoring and reporting of, as well as track the effectiveness of, mitigation actions. That said, studies have shown discrepancies between BU and TD estimations. Performing TD monitoring with concurrent BU measurements at the same mines may help address the observed discrepancies and improve quantification methods to narrow the gap between BU and TD data. This paper presents a comparison of concurrent BU-TD methane emission fluxes monitored from two ventilation shafts at a southwestern Virginia, USA, longwall mine to complement the existing body of studies on satellite-based monitoring of coal mines by incorporating continuous ground-based monitoring with concurrent TD monitoring of methane emission fluxes to address the gap and provide valuable insights into temporal emission patterns. The shafts were monitored on multiple dates between 2020 and 2023. BU monitoring was performed at the exhaust fans, while TD quantifications used PRISMA hyperspectral satellite data and two different wind reanalysis datasets (i.e., GEOS-FP and ERA5). This mine also offered a unique opportunity for BU-TD comparisons before and after ventilation air methane (VAM) oxidizer operation, which was installed at one of the shafts in 2022. The results showed that TD-estimated mean fluxes were generally lower than BU data, which were attributed to quantification difficulties associated with the low albedo caused by heavy vegetation and the terrain of the area. However, despite the discrepancies in mean emissions, both the interquartile range and the data range of the distributions generally overlapped, and the estimates correctly showed the emission trends.

Virginia↗
USGS2025

Site-specific, extended ShakeMaps for earthquake engineering applications

The U.S. Geological Survey (USGS) routinely produces ShakeMaps of shaking intensity across the globe. Due to practical constraints, the number of response spectral periods was limited to three standard periods (0.3, 1.0, and 3.0 sec). We have recently developed the tools that are necessary to expand this functionality to include 22 periods (matching the current U.S. National Seismic Hazard Model periods) as well as the orientation-independent components (e.g., “RotD50”). We refer to ShakeMap products that include these extensions as “extended ShakeMaps.” The added level of complexity motivated us to also develop a user-friendly tool called the “ShakeMap Sampling Tool” (SST) that gives all the estimated shaking metrics for a specific location (or list of locations). Additionally, we develop a web application where users can input locations of interest and view/download the SST results. We further familiarize users with the concept of “Composite ShakeMaps.” For earthquakes sequences such as a mainshock and larger foreshocks and aftershocks, this provides a map of the maximum value of each shaking metric, which is useful for overall loss estimates, the full extent of ground failure triggering potential, and a better portrayal of the repeated shaking levels at a given point for a series of earthquakes. Such a site-specific shaking history facilitates earthquake forensics at building or infrastructure sites for which damage may be of concern, as described in the Disproportionate Damage Earthquake trigger specified in the IEBC (2018, Section 405.2.2) and in developing ATC-145 guidelines (Guidelines for Post-Earthquake Assessment, Repair, and Retrofit of Buildings). The composite ShakeMap can be combined with the SST for a variety of earthquake-hazard applications, such as systematically inferring triggering shaking estimates at specific sites of geotechnical interest for landsliding, liquefaction, and lateral-spreading hazards.

Kahramanmaraş↗
USGS2024

Geophysical identification of potential groundwater – surface water interactions in the Bonita Peak Mining District, San Juan County, Colorado

This report includes the results and interpretation of geophysical investigations to aid in the location, identification, and conceptual site model (CSM) development of the near surface groundwater and surface water hydrology of the Upper Cement Creek and California Gulch of the Animas Rivers within the Bonita Peak Mining Districts (BPMD). Mine influenced water (MIW) from numerous mine adits within the BPMD has altered the geochemistry and ecology of these ground and surface waters. Many studies of varying scales have investigated the geology, hydrogeology, geochemistry, and aqueous chemistry of these adits and impacted areas. Field investigations using non-invasive geophysical methods were conducted in the summer of 2019 to aid in identifying potential groundwater inputs to these streams and further the CSM of the near surface within and along these streams. Two field deployments were completed in August and September of 2019. In August 2019, Phase 1 was completed using multifrequency electromagnetic induction (EMI) and thermal infrared imaging (TIR) to measure the bulk electrical conductivity, magnetic susceptibility, and the temperature response at known and suspected groundwater seeps into the Upper Cement Creek and California Gulch stretches of these alpine streams. Results and interpretation from Phase 1 were used to co-locate stretches of these streams for further detailed fiber optic distributed temperature system (FODTS) monitoring of the temperature dynamics for submeter identification of potential groundwater seeps during Phase 2. Phase 2 included additional EMI measurements, FODTS measurements, and repeated some Phase 1 measurements for improved data quality. All these data were used for the combined qualitative interpretation identifying seep locations. The data from both Phase 1 and Phase 2 are available in the USGS data release at: https://www.sciencebase.gov/catalog/item/5f18a30082cef313ed84879b. These data include detailed near-surface electrical conductivity and magnetic susceptibility mapping, discrete TIR imaging, and several kilometers of FODTS data. These data were georeferenced and used to identify potential groundwater seeps into these streams. The results indicate a complex near surface geophysical and hydrogeologic regime. Seeps are classified and spatially correlated to the EMI, TIR, and geochemical sampling results. This report is a thorough compilation of the geophysical field efforts in the summer of 2019, the data processing, and interpretations of these data. These data alone should not be used for final interpretations of the complex hydrogeology within these alpine systems. Rather, these data should be combined with other known data from various investigations to further refine the CSM and guide decisions for future investigations, watershed protective measures, and other site decisions. Geographical information system (GIS) layers of these results are available to project managers for such integration and combined interpretations.

Colorado↗
USGS2026

Simulated ground motion dataset in the Azores Plateau, Portugal, on rock and soil sites

Building on a previously developed bedrock dataset, this study extends the Azores Plateau ground motion simulations to include soil-amplified records and introduces a comprehensive validation framework. Soil amplification is modeled using one-dimensional soil profiles. A stochastic source-based approach is employed to generate the dataset, incorporating randomization of input-model parameters to account for the aleatory uncertainty in seismic activity. The accuracy of the dataset is verified through a comprehensive validation framework, showing that the randomization effectively captures variance and inter-period correlation observed in records. This work provides a robust dataset for advancing seismic hazard and risk assessment in the Azores Plateau.

central and eastern Azores islands↗
USGS2025

Lithium from magma to mine in an early Yellowstone hotspot caldera

Renewable energy technologies rely on the extraction of metals not historically in high demand, such as lithium (Li), for which ore deposit models are incompletely understood. One of the world’s largest Li deposits is hosted in lake sediments of the 16.4 Ma McDermitt caldera, which formed during the early stages of Yellowstone hotspot volcanism in the western United States. Eruptive and posteruptive mobility of Li are major challenges in elucidating deposit formation. Melt inclusions preserved in quartz crystals provide a means to assess pre-eruptive magmatic Li contents. Concentrations of Li determined by ion microprobe for melt inclusions in a McDermitt rhyolite lava are 400−1350 ppm, compared to 20−70 ppm Li in matrix rhyolite glasses. Synthesis with melt inclusion data for eight additional calderas demonstrates a recurrence of Li-rich rhyolitic magmas (200−2000 ppm Li) in the western part of the Yellowstone hotspot track. However, unlike the multicyclic caldera complexes with overlapping fault networks that may have compromised Li retention, the McDermitt caldera remained a closed hydrologic system throughout its evolution. Modeling indicates 100 km 3 of resurgent magma could yield 25−150 Mt Li in a magmatic fluid and supports accumulation of Li-rich magmatic fluid in a closed intracaldera lake, followed by evaporative concentration and sequestration of Li within clay minerals to generate the McDermitt deposit.

Idaho, Nevada, Oregon, Wyoming↗
USGS2025

Aeromagnetic and magnetotelluric imaging of west-central Idaho and the Stibnite-Yellow Pine mining district: A regional to district perspective

Aeromagnetic and magnetotelluric (MT) data are used to better understand the geology and mineral resources near the Stibnite-Yellow Pine mining district in central Idaho. The reduced-to-pole (RTP) transformation of regional-scale aeromagnetic data shows that allochthonous island-arc rocks west of the Salmon River suture are significantly more magnetic than the Laurentian continental rocks east of the suture and that the granitoids of the Idaho batholith have moderate to low magnetization in both early, metaluminous, and late, peraluminous phases. Application of tilt derivative to aeromagnetic data highlights major crustal-scale structures. The 5-km upward continued magnetic data indicate island-arc rocks have deep magnetic sources. The 110-km-long MT profile images resistivity structure to depths around 30 km. At shallow depths, resistivity corresponds to mapped geologic units, with moderate resistivities underlying volcanic and roof-pendant metasedimentary rocks and moderate to high resistivities occurring beneath the Idaho batholith. Crustal-scale moderate resistivities beneath the suture image the results of tectonomagmatic processes that accompanied suturing and translating allochthonous terranes. Low resistivity values beneath and fringing the batholith are derived from metasedimentary rocks that may have served as a melt source and reductant during melt generation and provided metals during later ore formation. In the Stibnite-Yellow Pine mining district, a high-resolution aeromagnetic compilation is shown to correlate with mapped lithologies and mineral deposit-related structures. The RTP transform distinguishes magnetic and nonmagnetic granitoid phases of the Idaho batholith. The tilt derivative highlights metasedimentary rocks, some of which are favorable ore hosts. The Meadow Creek fault hosts the Stibnite and Hangar Flats deposits and is imaged as a magnetic low due to hydrothermal alteration. Reconstructions of magnetic anomaly offsets and orebodies indicate around 3 km of post-95 Ma dextral separation, with some or all of the offset inferred to postdate the main Au mineralization episode (61–66 Ma).

Idaho↗
USGS2025

Gas emissions from the Sulphur Bank Mercury Mine hydrothermal system, Clear Lake volcanic field, California

The Sulphur Bank Mercury Mine (SBMM) hydrothermal system offers insights into active degassing processes in the Clear Lake volcanic field (CLVF), a high-threat region based on its record of Holocene eruptions and proximity to populated areas. Here we present chemical and isotopic analyses of gas samples collected between 2015 and 2023, along with the first comprehensive CO 2 flux survey of the SBMM area conducted in 2023. Sampled gases are CO 2 - and CH 4 -rich (≥84 and 6 mol% in dry gas, respectively) with high mantle-derived helium contributions ( 3 He/ 4 He = 6.54–7.86 R C /R A ). Carbon isotopic compositions of CO 2 (δ 13 C = −10.0 to −9.5 ‰) and CH 4 (δ 13 C = −35.8 ‰) indicate mixed sources, with significant contributions from metamorphism of organic-rich Franciscan Complex rocks hosting the hydrothermal system. Modeling of gas compositions shows that scrubbing by interaction with air-saturated groundwater strongly influences observed compositional variability. From our CO₂ flux measurements, we estimate the deeply derived CO 2 emission rate from the SBMM hydrothermal area (0.2 km 2 ) at 240 t d −1 , comparable to many quiescently degassing volcanoes worldwide. We also provide a first-order estimate of CH 4 emissions at approximately 0.5 t d −1 . Our findings establish crucial baseline data for future volcanic monitoring efforts, enhancing detection capabilities for potential changes in this active hydrothermal system. This work contributes to the broader understanding of volatile contributions from volcanic and metamorphic sources to the global carbon budget, while highlighting the strong influence of bedrock geology on gas compositions in the CLVF.

California↗
USGS2025

Modeling byproduct and coproduct mine production and mineral substitution using multidimensional supply curves: Application to the Cu-Co-Ni system and beyond

Rapid demand growth is expected for many metals used in the energy transition. Many of these metals are byproducts of other commodities. Byproduct production’s price response is tied to host mineral economics, complicating its supply dynamics. Moreover, many of these metals are used in applications where the material properties desired are difficult to substitute; effectively, limiting how quickly demand can adapt to changes in commodity price. Previous work has demonstrated the interconnectivity of jointly produced mineral commodities from the supply side, where the copper–cobalt–nickel system was used and demand was assumed independent across commodities. Studies to understand byproduct-coproduct market interconnectivity on the demand side are limited, while studies on the interconnectivity of supply and demand simultaneously are even more so. We propose a modification to the multicommodity supply curve method to enable inter-commodity effects on demand simultaneous with supply. In batteries, high cobaltprices may push consumers to transition to high-nickel chemistries, causing the nickel demand surface to decrease with nickel price but increase with cobaltprice, creating a two-dimensional demand surface. Below cross-price elasticities of 0.05, inter-commodity effects were found to be negligible, potentially permitting exclusion of these effects for many commodities. This additional demand curve complexity introduces potential computation challenges alongside the capacity to model many interrelated commodity systems such as rare earth elements, ferroalloys, country-oriented subsidies or restrictions, and bifurcated sustainable metals markets. By presenting the work done on multicommodity supply surfaces to date and potential new directions, this work aims to catalyze the next round of innovative approaches to modeling jointly produced commodities.

Conference Paper↗