USGS Science⌕ Search

SEARCH · USGS Science

Results for “Transactions of the American Fisheries Society”

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.

1,688 records · Page 64Linked to original sources
USGS2025

Identification of representative earthquakes for probabilistic tsunami hazard analysis (PTHA) using earthquake rupture forecasts and machine learning

As probabilistic tsunami hazard analysis (PTHA) focuses more on assessments for localized, populous regions, techniques are needed to identify a subsample of representative earthquake ruptures to make the computational requirements for producing high-resolution hazard maps tractable. Moreover, the greatest epistemic uncertainty in seismic PTHA is related to source characterization, which is often poorly defined and subjective. We address these two salient issues by applying streamlined earthquake rupture forecasts (ERFs), based on combinatorial optimization methods, to an unsupervised machine learning workflow for identifying representative ruptures. ERFs determine the optimal distribution of a millennia-scale sample of earthquakes by inverting the observed slip rate on major faults. We use two previously developed combinatorial optimization ERFs, integer programming and greedy sequential, to produce the optimal location of ruptures with seismic moments sampled from a regional Gutenberg–Richter magnitude–frequency distribution. These ruptures in turn are used to calculate peak nearshore tsunami amplitude, using computationally efficient tsunami Green's functions. An unsupervised machine learning workflow is then used to identify a small subsample of the earthquakes input to ERFs for onshore PTHA analysis. We eliminate epistemic uncertainty related to source distribution under traditional PTHA analysis; in its place, a quantifiable, less subjective and generally smaller uncertainty related to the input to ERFs is included. The Nankai subduction zone is used as a test case, where previous ERFs have been conducted. Results indicate that the locations of representative earthquakes are sensitive to choice of magnitude–area relation and to whether a minimum cumulative stress objective is imposed on the fault. In general, incorporating ERFs into PTHA provide a physically self-consistent method to incorporate fault slip information in determining representative earthquakes for onshore PTHA, eliminating a major source of epistemic uncertainty.

Nankai subduction zone↗
USGS2025

Toward a new framework to evaluate process-based model configurations and quantify data worth prior to calibration

Model criticism, discrimination, and selection methods often rely on calibrated model outputs. Because calibration can be computationally expensive, model criticism can first be undertaken by assessing model outputs obtained from limited prior parameter ensembles. However, such prior-based methods are often heuristic and do not formalize the notion of balancing model consistency with data and model complexity (i.e., model adequacy). We present a new framework to discriminate among candidate models prior to calibration that formalizes prior-to-calibration model adequacy into a metric to implicitly balance prior model output data coverage with model complexity represented by prior output (co)variance. The prior model adequacy metric “Mahalanobis distance deviation” quantifies the deviation of (a) the set of squared Mahalanobis distances of data from a prior model output distribution from (b) the set of squared Mahalanobis distances of data from their own distribution. A new data worth metric “discernment value” is also presented which quantifies the value of data for screening less-adequate models prior to calibration. Discernment value is calculated from the change in variance of a weighted average of prior model outputs from all candidate models due to less-adequate model outputs receiving lower weight. The framework is demonstrated using a one-dimensional groundwater flow model with eight possible configurations. A synthetic data network is used to test the framework. Results show the framework identifies the candidate models most similar to the true model used to create the synthetic data. Discernment values show variation in the value of different data types and locations for screening less-adequate models.

Water Resources Research↗
USGS2025

Mapping bedrock outcrops in the Sierra Nevada Mountains (California, USA) using machine learning

Accurate, high-resolution maps of bedrock outcrops can be valuable for applications such as models of land–atmosphere interactions, mineral assessments, ecosystem mapping, and hazard mapping. The increasing availability of high-resolution imagery can be coupled with machine learning techniques to improve regional bedrock outcrop maps. In the United States, the existing 30 m U.S. Geological Survey (USGS) National Land Cover Database (NLCD) tends to misestimate extents of barren land, which includes bedrock outcrops. This impacts many calculations beyond bedrock mapping, including soil carbon storage, hydrologic modeling, and erosion susceptibility. Here, we tested if a machine learning (ML) model could more accurately map exposed bedrock than NLCD across the entire Sierra Nevada Mountains (California, USA). The ML model was trained to identify pixels that are likely bedrock from 0.6 m imagery from the National Agriculture Imagery Program (NAIP). First, we labeled exposed bedrock at twenty sites covering more than 83 km 2 (0.13%) of the Sierra Nevada region. These labels were then used to train and test the model, which gave 83% precision and 78% recall, with a 90% overall accuracy of correctly predicting bedrock. We used the trained model to map bedrock outcrops across the entire Sierra Nevada region and compared the ML map with the NLCD map. At the twenty labeled sites, we found the NLCD barren land class, even though it includes more than just bedrock outcrops, accounted for only 41% and 40% of mapped bedrock from our labels and ML predictions, respectively. This substantial difference illustrates that ML bedrock models can have a role in improving land-cover maps, like NLCD, for a range of science applications.

California↗
USGS2025

lasertram: A Python library for time resolved analysis of laser ablation inductively coupled plasma mass spectrometry data

Laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) data has a wide variety of uses in the geosciences for in-situ chemical analysis of complex natural materials. Improvements to instrument capabilities and operating software have drastically reduced the time required to generate large volumes of data relative to previous methodologies. Raw data from LA-ICP-MS, however, is in counts per unit time (typically counts per second), not elemental concentrations and converting these count ratesto concentrations requires additional processing. For complex materials where the ablated volume may contain a range of material compositions, a moderate amount of user input is also required if appropriate concentrations are to be accurately calculated. In geologic materials such as glasses and minerals that potentially have numerous heterogeneities (e.g., microlites or other inclusions) within them, this is typically determiningwhether the total ablation signal should be filtered to remove these heterogeneities. This necessitates that the LA-ICP-MS data processing pipeline is one that is not automated, but is also designed to enable rapid and efficient processing of large volumes of data. Here we introduce , a Python library for the time resolved analysis of LA-ICP-MS data. We outline its mathematical theory, code structure, and provide an example of how it can be used to provide the time resolved analysis necessitated by LA-ICP-MS data of complex geologic materials. Throughout the pipeline we show how metadata and data are incrementally added to the objects created such that virtually any aspect of an experiment may be interrogated and its quality assessed. We also show, that when combined with other Python libraries for building graphical user interfaces, it can be utilized outside of a pure scripting environment. can be found at https://doi.org/10.5066/P1DZUR3Z

Applied Computing and Geosciences↗
USGS2024

Trimming the UCERF3-TD logic tree: Model order reduction for an earthquake rupture forecast considering loss exceedance

The Uniform California Earthquake Rupture Forecast version 3-Time Dependent depicts California’s seismic faults and their activity. Its logic tree has 5760 leaves. Considering 30 more model combinations related to ground motion produces 172,800 distinct models representing so-called epistemic uncertainties. To calculate risk to a portfolio of buildings, one also considers millions of earthquakes and spatially correlated ground-motion variability. We offer a tree-trimming technique that retains the probability distribution of portfolio loss and identifies the leading sources of uncertainty for further study. We applied it to a California statewide building portfolio and various levels of nonexceedance probability between one in 100 and one in 2500. We trimmed the logic tree from 172,800 leaves to as few as 15. The result: a supercomputer that would otherwise run 24 h to estimate the distribution of one-in-250-year loss can calculate it in moments with the reduced-order model. Others can use the reduced-order model to calculate risk to different California portfolios, and scientists can prioritize study to reduce the remaining epistemic uncertainty.

Earthquake Spectra↗
USGS2024

Self-guided decision support groundwater modelling with Python

The GMDSI tutorial notebooks repository provides learners with a comprehensive set of tutorials for self-guided training on decision-support groundwater modelling using Python-based tools. Although targeted at groundwater modelling, they are based around model-agnostic tools and readily transferable to other environmental modelling workflows. The tutorials are divided into three parts. The first covers fundamental theoretical concepts. These are intended as background reading for reference on an as-needed basis. Tutorials in the second part introduce learners to some of the core concepts parameter estimation in a groundwater modelling context, as well as providing a gentle introduction to the PEST, PEST++ and pyEMU software. Lastly, the third part demonstrates how to implement highly-parameterized applied decision-support modelling workflows. The tutorials aim to provide examples of both “how to use” the software as well as “how to think” about using the software. A key advantage to using notebooks in this context is that the workflows described run the same code as practitioners would run on a large-scale real- world application. Using a small synthetic model facilitates rapid progression through the workflow.

Journal of Open Source Education↗
USGS2026

Nonstationary demographic state-space models using unreplicated counts for species undergoing environmental stressors

A fundamental task in ecological statistics is to estimate abundance and growth rate distributions from wildlife monitoring data to inform conservation management. Modeling time series of wildlife populations presents a number of challenges from both statistical and ecological perspectives, including discreteness; lack of replication; nonstationarity; and observation, demographic, and other phenomenological processes. Nonstationary dynamics are often exhibited by populations undergoing environmental stressors. Models must account for these characteristics to produce reliable estimates of abundance and trends, yet estimation can be challenging with unreplicated data. We propose nonstationary demographic state-space models using unreplicated counts for populations undergoing environmental stressors. A reduced growth rate model matches the complexity of the unreplicated count data, and a fecundity bound on growth rate distributions allows the separation of processes affecting growth rates like environmental stressors from those affecting abundance external to growth rates like migration. NDSSMs allow for the embedding of nonstationary model components, and we explore the use of changepoints, volatility clustering, and migration processes. We apply the proposed nonstationary models in case studies of herons affected by predator/competitor reestablishment and three bat species affected by a fungal pathogen causing white-nose syndrome. Nonstationary models outperform stationary models and generalized linear mixed effects models according to model scoring and visual inspection of predictions, and provide estimates more consistent with published values. Incorporating migration improves model fit universally, even with approximate one-way immigration, most likely because populations are extirpated, recolonized, and increase multiple-fold over the upper bound set by species fecundity. In addition, estimates of the timing and severity of the environmental stressor differed for models with migration. Including nonstationary and demographic components in a fecundity-bounded growth rate model improves inference and benefits interpretability of hyperparameters. In turn, this adjusts uncertainties in predictions of abundance and growth rates over time, providing the ingredients needed for informed conservation analysis and for directing future monitoring of at-risk species.

Journal of Agricultural, Biological and Environmen↗
USGS2025

Three-dimensional seismic velocity models for the San Francisco Bay region, California from joint body-wave and surface-wave tomography validated by waveform simulation

A high-resolution seismic velocity model for both P and S waves is essential for the San Francisco Bay (SFB) region to accurately simulate earthquake ground motions and assess seismic hazards in this high-risk region. We present new three-dimensional P-wave and S-wave velocity (Vp and Vs) models of the SFB region developed using joint tomographic inversion with a rich set of earthquake and active-source body-wave data and ambient-noise surface-wave data. We designed a dynamic balance strategy to effectively balance the contributions of body-wave and surface-wave data during each iteration of the inversion. Owing to the complementary resolution of body-wave and surface-wave data, our new models are significantly improved compared to previous tomographic models. Major geologic features in the SFB region, including low-velocity basins, high-velocity Franciscan Complex, granite, and gabbro rocks, and across-fault velocity contrasts, are imaged in our Vp and Vs models. Compared to the geology-based U.S. Geological Survey velocity model, our new model has better defined structures for some important areas that are of great concern for seismic hazard assessment, such as the Sacramento-San Joaquin Delta, and the 2014 magnitude 6.0 South Napa earthquake area. We validated our model by simulating waveforms for a set of moderate-magnitude earthquakes. The good performance of our model in waveform simulation indicates its potential for improving ground motion simulation and seismic hazard assessment in the SFB region.

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

The Sand AbrasioN Device for Aeolian Research (SANDAR): A new experimental device for investigating how wind transport affects sand on Earth and Mars

On Earth and Mars, aeolian transport causes sand grains to become abraded, resulting in mineralogic and textural changes. Understanding how sands evolve, or mature, with transport via experimental studies is important for understanding the origins, geologic history, and cycling of sediments, as well as dust production. Previous experimental works have used a variety of methods to simulate aeolian transport in the laboratory, but practical limitations and similitude concerns have limited such research. Here, we present and validate the Sand AbrasioN Device for Aeolian Research (SANDAR), a modified air mill that uses pressurized air to circulate sand around a small abrasion chamber, simulating the effects of aeolian transport. This device is re-circulating to simulate long-distance transport, and it allows for repeated analyses of well-constrained sediment samples, revealing their evolution over time. It is compatible with the grain sizes (74–500 μm) and grain impact velocities (∼0.6–3.7 m/s) typically expected for natural aeolian environments, and is also adaptable for diverse applications simulating different wind conditions. We show that the SANDAR achieves similitude of kinetic energy with respect to saltating sand on both Earth and Mars. SEM and optical microscope imaging reveal that the SANDAR produces microtextures on the surfaces of sand grains similar to those found with natural aeolian transport, demonstrating that it effectively simulates the mechanical effects of aeolian processes. Thus, the SANDAR is a valid tool for use in experimental research to improve our understanding of sedimentary processes across the Solar System.

Aeolian Research↗
USGS2025

Persistent deep long-period seismicity near the Lassen Volcanic Center

Deep long-period (DLP) earthquakes have been observed at many volcanic settings around the world and linked to the magmatic processes that drive volcanic unrest. At the Lassen Volcanic Center (LVC) of the Cascade arc, limited detection of DLP activity hinders classification of anomalous seismic behavior and its relationship to the LVC magmatic system. This study uses a template matching approach with seismic data from temporary nodal and permanent stations to detect and locate DLP earthquakes near the LVC between 2017 and 2024. Within the DLP catalog of 611 events, a transition occurs from scattered (2017–2020) to oscillatory (2021–2024) occurrence rates. During the oscillatory period, regional earthquakes with large amplitude velocity waveforms observed near the LVC are associated with abrupt changes in DLP occurrence rates. Continued monitoring of DLP activity has the potential to better define the processes that drive volcanic unrest at the LVC in the future.

California↗
USGS2024

Updated three-dimensional temperature maps for the Great Basin, USA

As part of the periodic update of the geothermal energy assessments for the USA (e.g., last update by Williams and others, 2008), a new three-dimensional temperature map has been constructed for the Great Basin, USA. Williams and DeAngelo (2011) identified uncertainty in estimates of conductive heat flow near land surface as the largest contributor to uncertainty in previously published temperature maps. The new temperature maps incorporate new conductive heat flow estimates developed by DeAngelo and others (2023). Predicted temperatures at depth are compared with representative measurements (for conductively dominated conditions), showing good agreement under relatively simple uniform conditions. Inputs included radiogenic heat production for all layers of 1.89 μW/m 3 , effective bulk thermal conductivity of 2.7 W/m/°C for all rocks underlying sedimentary basins, and a previously published (Williams and DeAngelo, 2011) empirically driven estimate of increasing thermal conductivity with depth in sedimentary sequences. The resulting three-dimensional temperature model is published in a USGS data release associated with this manuscript (Burns and others, 2023).

Arizona, California, Idaho, Nevada, Oregon, Utah↗
USGS2026

Cultural benefits knowledge and social learning in Elwha River dam removal and ecosystem decision-making

In a retrospective case study of the seminal Elwha River dam removal in Washington State, USA, we explore pathways through which ecosystem services (ES) knowledge can inform environmental decision-making. We specifically investigate how knowledge of cultural benefits linked to the Elwha River ecosystem was made available. Through interviews and historical document review we identify diverse forms of cultural benefits knowledge (i.e., how individuals or groups know the contributions of an ecosystem to their identities, experiences, and capabilities), including both knowledge in the form of practices – embodied or enacted forms of knowledge, and knowledge in the form of products – quantitative or written documentation. Recognizing and learning from the diverse forms in which benefits knowledges are conveyed can support more accurate and equitable consideration of the diverse values, well-beings, and human-nature relationships linked to ecosystems. With this case study, we illuminate pathways for meaningful consideration of multiple knowledge forms, and the plural ecosystem values they convey, with relevance for environmental decision processes in the Puget Sound and beyond.

Washington↗
USGS2026

GST-1: A high-resolution global sediment thickness model

Global Sediment Thickness 1 (GST-1) is a high-resolution sedimentary thickness model calculated on a 0.125° x 0.125° grid. It modifies the sediment thickness of the 1° x 1° Earth Crustal Model 1 (ECM1) by means of 3D inversions of free air gravity anomalies. GST-1 is calculated by performing structural inversions on high-density contrasts across two crustal boundaries: the sediment – basement interface and the crystalline crust – upper mantle interface. The inversions are calculated in each of ten overlapping 3D models that span the globe. These ten models are merged to obtain the GST-1 global model, providing an eight-fold increase in lateral spatial resolution in comparison with ECM1 and CRUST 1.0. Our sediment thickness model exploits the nearly continuous sampling of gravity data when compared to the irregular, sparse sampling of seismic refraction data. Sediment thickness values in GST-1 are in excellent agreement with independently derived cross sections from well-studied sedimentary basins, and within expected resolution limits of seismic refraction data. GST-1 offers a robust, high resolution global model of sedimentary thickness to support studies of sedimentary basins.

Tectonophysics↗
USGS2024

Testing megathrust rupture models using tsunami deposits

The 26 January 1700 CE Cascadia subduction zone earthquake ruptured much of the plate boundary and generated a tsunami that deposited sand in coastal marshes from northern California to Vancouver Island. Although the depositional record of tsunami inundation is extensive in some of these marshes, few sites have been investigated in enough detail to map the inland extent of sand deposition and depict variability in tsunami deposit thickness and grain size. We collected 129 cores in marshes of the Salmon River estuary in Oregon and reanalyzed 114 core logs from a 1987–88 study that mapped the inland extent of circa 1700 CE sandy tsunami deposits. The ca. 1700 CE tsunami deposit in the Salmon River estuary is easily recognized in cores ≤1 m deep in which a buried marsh peat is overlain by a well sorted sand bed with a sharp lower contact that thins and fines inland. We use tsunami deposit data and models of sandy tsunami sediment transport (using Delft3D-FLOW) to test 15 rupture models that could represent a ca. 1700 CE earthquake. At least 12–16 m of slip offshore of the Salmon River, which results in 0.8–1.0 m of coastal coseismic subsidence, is required to match the ca. 1700 CE sand deposit's inland extent, which is consistent with models of heterogeneous megathrust slip in ca. 1700 CE. Our methods of detailed tsunami deposit mapping, combined with sediment transport modeling, can be used to test models of megathrust ruptures and their tsunamis to potentially improve earthquake and tsunami hazard assessments.

JGR Earth Surface↗
USGS2012

Development of vegetation based soil quality indices for mineralized terrane in arid and semi-arid regions

Soil quality indices (SQIs) are often management driven and attempt to describe key relationships between above- and below-ground parameters. In terrestrial systems, indices that were initially developed and modified for agroecosystems have been applied to non-agricultural systems in increasing number. We develop an SQI in arid and semi-arid ecosystems of the Western US impacted by different types of geologic mineralization using the relationship between vegetation community parameters and soil abiotic and biotic properties. We analyze these relations in soils associated with three different mineralization types: podiform chromite, Cu/Mo porphyry, and acid-sulfate gold vein systems at four different sites in California and Nevada. Soil samples were collected from undisturbed soils in both mineralized and nearby unmineralized substrates as well as from waste rock and tailings. Aboveground net primary productivity (ANPP), canopy cover and shrub density were measured for the vegetative communities. Minimum data sets were developed based on correlations between the soil and vegetation parameters, refined using principal components analysis, scored using non-linear functions, and combined into an overall SQI. The indices are comprised of one or two microbial parameters and three to six abiotic parameters, the latter consisting of nutrients and metals. Given the preliminary development of this approach, the parameters and combinations to arrive at an SQI for a given site cannot at this time be correlated or compared with that of another site. This SQI approach provides a means of quantifying disturbed ecosystem recovery resulting from mining, and could be applied to other disturbances in a way that readily distills the information for potential use by land managers. However, severely disturbed areas with little to no aboveground biomass, such as unreclaimed tailings, have likely crossed an ecological threshold that precludes the use of this type of monitoring tool.

California↗
USGS2023

Hidden system identification: Basin modeling as a tool for examining sedimentary geothermal resource potential

Three-dimensional (3D) geologic and temperature models have been developed for the onshore U.S. Gulf Coast. The results from these models identify areas of moderate- to high-temperature (90°-150°C and >150°C; respectively) geothermal resources at depths <6 km. This modeling study addresses the fundamental challenge of predicting where opportune temperature and lithology coincide. Unlike traditional geothermal systems with surface expressions of hydrothermal circulation (e.g., hot springs, fumaroles, sinter), sedimentary geothermal systems (SGS) are generally hidden. Historically, simplified efforts to predict subsurface temperatures in sedimentary basins have focused on linear temperature extrapolation that does not consider the variable thermal properties of different lithologies or lithologic changes with depth (e.g., compaction, lithification). Therefore, the need to understand basin architecture and predict temperatures in 3D within SGS is paramount to identifying geothermal resources and determining economic feasibility. Basin modeling software has long been used to characterize the subsurface conditions of sedimentary basins, including temperature, in the pursuit of finding hydrocarbons. This tool can also be adapted to evaluate the potential of geothermal resources in a sedimentary basin by predicting the confluence of desirable temperatures and reservoir lithologies. In this work, PetroMod basin modeling software was used to create a regional geologic model of the onshore U.S. Gulf Coast, covering over 500,000 km 2 calibrated to temperature data from wells. Inputs include structural surfaces from commercial databases, lithology information derived from published literature, and corrected bottom-hole temperatures (BHT) from over 6,000 wells. The resulting 3D geologic model can be used to predict temperatures throughout the basin. Maps were exported showing the depth, depositional unit, and reservoir lithology at which temperatures of 90°C and 150°C were reached, revealing over 400,000 km 2 of moderate- to high-temperature resources at depths <6 km. These maps function as a first-order screening tool to identify areas where low-, moderate-, or high-grade resource potential may exist, based on temperature and if optimal reservoir lithologies or depositional units of interest are present. Depending on the success criteria of a project, the same maps can be exported for any isotherm or incorporate other 1407 Gardner and Birdwell subsurface properties. The methodology employed in this work can be applied in any sedimentary basin with available subsurface data. Further calibration incorporating other data, including pressure and porosity, can expand the utility of basin modeling for geothermal evaluations. Basin modeling is a powerful but underutilized tool for identifying prospective geothermal resources in sedimentary basins.

Conference Paper↗
USGS2025

Understanding market sensitivity: Estimation of supply and demand elasticities for non-fuel minerals

In today’s rapidly changing economic landscape, understanding market responsiveness to price changes and the factors influencing commodity prices has become increasingly relevant. Price elasticities serve as indicators of how variations in market conditions affect supply and demand, providing insights into the sensitivity of commodity markets to price fluctuations. This paper presents a comprehensive analysis of price elasticities of supply and demand for 74 non-fuel mineral commodities including precious metals, base metals, minor metals, and industrial minerals that are utilized across various industries. We employ various econometric techniques, including fixed effects models for panel data and two-stage dynamic ordinary least squares (2S-DOLS) alongside autoregressive distributed lag (ARDL) models for time series analysis, to derive robust estimates of price elasticities. Our findings reveal variability in elasticities among different commodities and indicate that all studied mineral commodities exhibit price inelastic supply and demand in the short run, which we define as one year for the purposes of our analysis, given that the data is all annual. This research provides original estimates of price responsiveness for a wide range of commodities that have not been previously addressed in the literature, thereby enhancing the understanding of market dynamics in the mineral sector. Given that price elasticities can be influenced by factors such as market structure, technological advancements, mining costs, and industry-specific demand drivers, we use variables that serve as proxies for these factors.

Mineral Economics↗