EconBase
← Back to paper

Quantifying Demand Shocks in the Green and Digital Transition

Extracted main text — title through conclusion, appendix excluded. This is what our citation measures are computed over, published so the extraction can be checked by eye.

60,142 characters · 12 sections · 59 citation commands

Rendered from LaTeX for readability, not typeset faithfully. Citation keys are highlighted; maths is left as source; figures, tables and equation environments are summarised rather than reproduced; unrecognised commands are greyed out so nothing is silently dropped. Email addresses are removed.

Quantifying Demand Shocks in the Green and Digital Transition

\doublespacing {\noindentAbstract: We use web search data to construct monthly indexes of derived demand for cobalt, copper, and nickel, which are key inputs in technologies driving the energy and digital transitions. We incorporate these indexes into Structural Vector Autoregressive (SVAR) models of global metal markets and identify structural shocks using zero, sign, and magnitude restrictions. This approach disentangles supply shocks from several demand-side drivers of metal prices and isolates a transition demand (TD) shock linked to the diffusion of metal-intensive technologies. We find that TD shocks generate persistent price effects, especially for copper and nickel, whereas supply and metal-specific demand shocks are more immediate and less persistent.

\noindentKeywords: Google Trends; Metals; Structural VAR; Sign restrictions; Transition Demand Shock.

\pagenumbering{arabic} \doublespacing

Introduction

Raw materials such as cobalt, copper, and nickel play a central role in the energy and digital transitions, hereafter referred to as the twin transition. They are key inputs in technologies such as electric vehicles, renewable energy systems, and data centre infrastructure, and demand for these metals has grown rapidly in recent years. At the same time, their supply is geographically concentrated, making their prices increasingly sensitive to shifts in global demand and geopolitical developments. Understanding the drivers of metal prices is therefore essential for both macroeconomic analysis and policy.

Data on the production and consumption of most commodities -- including metals, agricultural products, and energy -- are typically released with long publication lags and, in some cases, are only available at an annual frequency. This limits their usefulness for real-time monitoring of price drivers and for short- to medium-term analysis of prices and demands BRTFore. The advent of big data has provided researchers with new data sources and methods to compensate for the lack of timely, high-frequency official statistics abraham2022big,varian2014big. Leading examples include the use of web search data choi2012predicting,FerraraSimoni,Moller, scanner data Ngscanner, satellite imagery donaldson2016view, and news content baker2016measuring,barbaglia2023forecasting.

In this paper, we construct novel monthly Derived Demand Indexes (DDIs) for cobalt, copper, and nickel using web search activity data from Google. We extract these indexes from a high-dimensional set of Google Trends series using targeted principal component methods. Our approach uses search queries related to products that utilize these metals as production inputs, with a particular emphasis on technologies driving the energy and digital transitions. Online search volumes provide a measure of consumers’ purchase intentions goel2010predicting,vosen2011forecasting. For instance, Moller show that web search activity captures people’s intentions to buy a house and thereby serves as a proxy for housing demand. This interpretation of web search activity is also related to the literature that uses it to proxy investor attention and its impact on asset prices andrei2015investor,da2011search,vozlyublennaia2014. We therefore build on this literature to interpret our indexes as proxies for derived demand.

The concept of derived demand is based on the idea that the demand for cobalt, copper, and nickel depends on their use as intermediate goods in the production of final goods such as batteries and smart devices lowinger1984product,verleger2011margin. In our framework the relevant source of final demand may change over time. This is consistent with the view that the “marginal market” -- that determines the price of a metal -- might change over time as the shares of products in the consumption bundle shift verleger2011margin. For example, the end-use shares of cobalt have shifted substantially because of cobalt's increased use in battery manufacturing graedel2022us,miatto2020rise; therefore, the battery market has become an increasingly important marginal market, whose dynamics affect the price of metals.\footnote{Several papers rely on the notion of derived demand to model energy and metal markets and to forecast their prices baumeister2018product,duarte2021commodity,ederington2019review,issler2014using,stuermer2017industrialization.}

The resulting DDIs are incorporated into monthly Structural Vector Autoregressive (SVAR) models of the global markets for cobalt, copper, and nickel. We identify structural shocks using a combination of zero, sign, and magnitude restrictions, which allows us to disentangle supply shocks from multiple demand-side shocks that drive the real prices of these metals. In particular, we isolate a transition demand (TD) shock linked to the uptake of metal-intensive technologies in the energy and digital transitions. Our framework enables a quantitative assessment of the relative contribution of each structural driver to price dynamics and highlights the growing macroeconomic relevance of TD shocks.

This paper contributes to the literature developing SVAR models for commodity markets, including those focused on energy BHaer,KilianAER2009 and metals romani2024understanding,BaumeisterSpecialFocus, as well as studies that aim to isolate shocks related to the energy transition boer2024energy,cross2023drivers. However, to the best of our knowledge, this paper is the first to exploit monthly indexes based on web search data to identify demand shocks related to the energy and digital transitions in SVAR models.

commentHowever, following boer2024energy, we might think of two distinct sets of demand shifters: aggregate demand shocks that related to the global business cycle and metal-specific demand shocks related to the energy transition or to changes in demand for above-ground inventories due to fears of future supply shortages. All these shocks are expected to increase the price of metals. Literature - measure of investor attention da2011search,vozlyublennaia2014 - house prices Moller - lab mkt: Borup empl Marinucci unemp - nowcasting choi2012predicting,FerraraSimoni niesert2020can international macroeconomic aggregates - vosen2011forecasting,woo2019forecasting private consumption - tourism consump Cheng2021 - Yu global oil consumption

Our empirical analysis delivers two main findings. First, aggregate demand shocks remain a major driver of metal prices, but transition demand (TD) shocks generate persistent price responses that build gradually and decay slowly, especially for copper and nickel. Second, the importance of TD shocks has increased in recent years, with particular strong evidence for copper and more moderate for nickel.

The rest of the paper is organized as follows. Section (ref) describes the construction of DDIs. Section (ref) outlines the identification strategy used in SVAR models. Section (ref) presents the results, and Section (ref) concludes.

Derived demand indexes based on web search data

Google Trend Data

The absence of reliable data on metal demand, together with the difficulty of tracking such quantities, motivates our construction of a specific index. Specifically, we use data from Google Trends (GT) data to construct the monthly derived demand indexes (DDIs) for cobalt, copper, and nickel related to the twin transition.\footnote{Google's share of the global desktop search market has consistently exceeded 85% over the period 2009--2023, making it the dominant internet search provider statcounter2026.} GT is an online tool that provides time series of search volumes for specific queries and geographical regions. Google does not disclose the total count of queries by keyword; instead, it provides an index that reflects the proportion of searches within a given geographic region and time period. The index is then scaled from 0 to 100, with 100 representing the “peak popularity” of a search term.

Given that commodity prices are determined at the global level, our analysis considers a set of GT series capturing global search activity sampled at monthly frequency from January 2004 to December 2023. We construct the $DDI$ using GT series for products and components that contain battery metals, such as lithium-ion batteries, as well as downstream products that use these batteries, such as smartphones and laptops. We consider three main product categories: smart devices (i.e., laptops, smartphones, tablets, and smart vacuums) and related batteries, electric vehicles (EVs), and solar photovoltaic panels (solar PV).

We therefore identify eleven primitive queries related to the purchase or replacement of these products and components: “buy an electric car”, “buy solar panels”, “buy lithium battery”, “buy a smartphone”, “buy a laptop”, “buy a tablet”, “replace phone battery”, “replace laptop battery”, “replace tablet battery”, “buy a powerbank”, and “buy robot vacuums”. These queries capture consumer demand for products related to the twin transition containing cobalt, copper, and nickel. When a keyword is entered into GT, the platform suggests up to 25 semantically related queries. We therefore expand our initial list by adding related queries for each primitive query. After manually removing unrelated queries, the resulting set contains 121 search terms. This set is further reduced by dropping low-volume queries, defined as GT series with more than 5% of their values equal to zero. This yields a total of 108 GT series, which are grouped into three main categories capturing the search activity for smart devices and related batteries (84 GT), EVs (9 GT), and solar PV (15 GT).

GT data are based on a random sample of Google search activity and are not perfectly reproducible over time. GT series change over time and vary depending on the Internet Protocol address used to extract the data. To address this issue, we retrieved each GT series twelve times, using three machines and collecting the data on different days, and computed the sample average across these samples Eichnauer,Medeiros.

Construction of the Derived Demand Indexes

Once the GT series are collected, we use Principal Component Analysis (PCA) methods to summarise the information contained in GT series into a single metal-specific factor, which we interpret as our DDIs. Because PCA methods are not sensitive to non-stationarity, the GT series must first be preprocessed Borup,choi2012predicting,Marinucci,FerraraSimoni,Moller,Yu.

We initially address the large volatility in the GT series by log-transforming each series as $gt_{i,t} = \ln\left(1+GT_{i,t}\right)$ and then apply the hamilton2018you filter to extract the cyclical component of $gt_{i,t}$, where $i$ denotes individual GT series.

comment\begin{equation} gt_{i,t+24} = \alpha_{i,0} + \alpha_{1,i} gt_{i,t} + \alpha_{2,i} gt_{i,t-1} + \ldots + \alpha_{i,12} gt_{i,t-11} + \nu_{i,t+24}. \end{equation} The intuition is that the residuals $\hat{\nu}_{i,t+24}$ capture the factors that explain why forecasts of the variable at the 24-month horizon might be incorrect, such as a recession impacting the economy. hamilton2018you shows that residuals from regression $gt_{i,t+24} = \alpha_{i,0} + \alpha_{1,i} gt_{i,t} + \alpha_{2,i} gt_{i,t-1} + \ldots + \alpha_{i,12} gt_{i,t-11} + \nu_{i,t+24}$ are a reasonable proxy of the cyclical component of a wide range of data generating processes and that such residuals are stationary, provided that the fourth differences of the underlying series are stationary.

As shown by hamilton2018you, this approach can be applied to both seasonally and non-seasonally adjusted data and is robust to the seasonal adjustment procedure. The method therefore allows us to remove seasonality and induce stationarity before PCA without taking a stance on the nature of trends or the form of seasonality characterizing GT series.\footnote{We also considered the pre-processing approach in Borup and Moller. This is based on a sequential testing method to distinguish between stationarity, deterministic trend and unit root behavior; however, it is not effective in our setting, as some of the GT still showed a trend after the procedure was applied.}

We then construct the DDIs using Targeted PCA (TPCA). This method, introduced by BaiNg, aims to filter out the noise from the GT data by reducing the number of variables fed into the PCA, thereby improving the estimation of the latent factor. TPCA proceeds in two steps: (i) the selection of GT series targeted to the variable of interest, (ii) the computation of a latent factor via standard PCA. In the first step, we select relevant GT series with a soft-thresholding rule based on the Elastic-Net of zou2005regularization.

Let $\widetilde{\mathbf{gt}}_t = \left(\hat{\nu}_{1,t},\hat{\nu}_{2,t},\ldots,\hat{\nu}_{G,t}\right)^\prime$ be the $(G\times1)$ vector containing the cyclical component of the GT series extracted with the method of hamilton2018you, with $G=108$. We use the Elastic Net to select a subset of query series $\widetilde{\mathbf{gt}}_{m,t}\in \widetilde{\mathbf{gt}}_t$ that are most relevant for explaining the cyclical component of the log real price of metal $m$, extracted using the Hamilton filter. Following Borup and Moller, we tune the Elastic Net parameters to assign equal weight to the LASSO and ridge penalty ($\alpha=0.5$), and to return a vector of ten GT series for each metal.\footnote{We select ten series as a standard choice in the Google Trend literature. In the Supplement (available from the authors upon request), we show that the results are robust when selecting the number of series with an information criterion.} This procedure yields three possibly overlapping sets of selected queries: $\widetilde{\mathbf{gt}}_{m,t}$ for $m=\text{cobalt (CO), copper (CU), and nickel (NI)}$. The metal-specific DDIs are then estimated as the first principal component of each of these vectors and are denoted as $DDI_{CO,t},DDI_{CU,t}$, and $DDI_{NI,t}$.

Validation of the DDI

Before using the DDIs in the SVAR analysis, we evaluate whether they capture transition-related demand for metals. We do so by comparing the indexes with metal prices and, in the Supplement, with sales data and available measures of metal consumption.

figure[figure omitted — 653 chars of source]

Figure (ref) plots, for each metal, the DDI (colored line) together with the cyclical component of the corresponding log real price (black line), where the latter is extracted using the hamilton2018you filter. In all cases, the two series exhibit a strong contemporaneous correlation, supporting the information content of the DDIs and the effectiveness of the TPCA method. Specifically, cobalt shows the strongest correlation (74.5%), followed by nickel (72.1%), while the correlation for copper is 49.6%. The top left panel illustrates how $DDI_{CO,t}$ reflects the sharp price surges in 2022, driven by soaring demand for metals used in EV production and a sluggish supply response. Similarly, the bottom panel shows that $DDI_{NI,t}$ captures the trough in nickel prices in 2008, associated with the global financial crisis.

In all panels, a clear discontinuity emerges during the COVID-19 pandemic. This is confirmed by the fact that excluding the COVID-19 period increases the correlation between the DDIs and prices for all three metals. The effect is particularly pronounced for copper, for which the correlation rises by about 13%. This pattern suggests that the relationship between search-based indicators and the underlying economic mechanisms temporarily weakened. During the pandemic, mobility restrictions and stay-at-home policies increased online search activity, thereby raising DDIs even as underlying metal demand weakened with the contraction in global economic activity. At the same time, metal prices declined because industrial demand collapsed, while inventories likely buffered short-run supply adjustments. Consequently, the usual co-movement between DDIs and metal prices temporarily weakened because search activity and realized demand moved in opposite directions. Additional validation results are reported in the Supplement (available from the authors upon request).

\FloatBarrier

comment\begin{table}[h!] \caption{Cross-correlations between $PPI_t$ and $DDI_{m,t+k}$} \scriptsize {3pt} \resizebox{\textwidth}{!}{ \begin{tabular}{l*{25}{S[table-format=1.3]}} \toprule & \multicolumn{12}{c}{$k<0$} & \multicolumn{1}{c}{$k=0$} & \multicolumn{12}{c}{$k>0$} \\ \cmidrule(lr){2-13} \cmidrule(lr){14-14} \cmidrule(lr){15-26} Metal & {$-12$} & {$-11$} & {$-10$} & {$-9$} & {$-8$} & {$-7$} & {$-6$} & {$-5$} & {$-4$} & {$-3$} & {$-2$} & {$-1$} & {$0$} & {$+1$} & {$+2$} & {$+3$} & {$+4$} & {$+5$} & {$+6$} & {$+7$} & {$+8$} & {$+9$} & {$+10$} & {$+11$} & {$+12$} \\ \midrule Cobalt & 0.498 & 0.513 & 0.526 & 0.539 & { 0.544} & 0.537 & 0.524 & 0.505 & 0.491 & 0.486 & 0.475 & 0.454 & 0.431 & 0.427 & 0.426 & 0.425 & 0.415 & 0.410 & 0.406 & 0.405 & 0.397 & 0.395 & 0.390 & 0.393 & 0.406 \\ Copper & 0.359 & 0.370 & 0.381 & 0.411 & 0.419 & 0.419 & 0.422 & { 0.427} & 0.424 & { 0.427} & 0.420 & 0.393 & 0.369 & 0.362 & 0.358 & 0.351 & 0.348 & 0.337 & 0.332 & 0.340 & 0.331 & 0.291 & 0.253 & 0.204 & 0.167 \\ Nickel & { 0.574} & 0.574 & 0.567 & 0.557 & 0.539 & 0.508 & 0.465 & 0.420 & 0.373 & 0.324 & 0.268 & 0.216 & 0.163 & 0.141 & 0.111 & 0.083 & 0.042 & 0.012 & -0.023 & -0.048 & -0.087 & -0.127 & -0.160 & -0.184 & -0.196 \\ \bottomrule \end{tabular} } \begin{minipage}{0.9\textwidth} \scriptsizeNotes: The table reports cross-correlations between the cyclical component of real prices ($c_{m,t}$) and $DDI_{m,t+k}$ for $k=-12,\ldots,+12$ over 12/2006--12/2023. Negative (positive) $k$ indicates that DDI leads (lags) prices. Bold values denote peak correlations. \end{minipage} \end{table} For all three minerals, the proxies deviate from the cyclical components of log real prices during the COVID-19 pandemic. This is perfectly explainable because in that moment this new source of demand increased because people forced to stay home bought a lot of devices (this peak is particularly noticeable for copper) . However, the price is mainly driven by the other uses of the minerals which were not hit by this positive effect. However, this should no be regarded as a flow. Indeed, the pandemic represented a particular moment. On one \begin{figure}[h!] \caption{Correlation between DDI and the cyclical component of Computer and Eletronics IP from 12/2006 to 12/2023.} \\ \caption*{\scriptsizeNotes: See notes to Figure (ref).} \end{figure} \begin{figure}[h!] \caption{Correlation between DDI and the cyclical component of Electrical Equipment IP from 12/2006 to 12/2023.} \\ \caption*{\scriptsizeNotes: See notes to Figure (ref).} \end{figure} \begin{figure}[h!] \caption{Correlation between DDI and the cyclical component of Battery Manufacturing IP from 12/2006 to 12/2023.} \\ \caption*{\scriptsizeNotes: See notes to Figure (ref).} \end{figure} \begin{table}[h!] \caption{Cross-correlations between $c_{m,t}$ and $DDI_{m,t+k}$} {4pt} \begin{tabular}{l*{13}{S[table-format=1.3]}} \toprule & \multicolumn{6}{c}{$k<0$} & \multicolumn{1}{c}{$k=0$} & \multicolumn{6}{c}{$k>0$} \\ \cmidrule(lr){2-7} \cmidrule(lr){8-8} \cmidrule(lr){9-14} Metal & {$-6$} & {$-5$} & {$-4$} & {$-3$} & {$-2$} & {$-1$} & {$0$} & {$+1$} & {$+2$} & {$+3$} & {$+4$} & {$+5$} & {$+6$} \\ \midrule Cobalt & 0.676 & 0.693 & 0.717 & 0.723 & 0.726 & 0.734 & { 0.745} & 0.738 & 0.726 & 0.712 & 0.696 & 0.684 & 0.673 \\ Copper & 0.539 & 0.549 & 0.553 & { 0.558} & 0.554 & 0.531 & 0.496 & 0.455 & 0.424 & 0.408 & 0.398 & 0.387 & 0.364 \\ Nickel & 0.589 & 0.615 & 0.642 & 0.661 & 0.681 & 0.708 & { 0.721} & 0.716 & 0.701 & 0.674 & 0.652 & 0.639 & 0.618 \\ \bottomrule \end{tabular} \begin{minipage}{0.9\textwidth} \scriptsizeNotes: The table reports cross-correlations between the cyclical component of real prices ($c_{m,t}$) and $DDI_{m,t+k}$ for $k=-6,\ldots,+6$ over 12/2006--12/2023. Negative (positive) $k$ indicates that DDI leads (lags) prices. Bold values denote peak correlations. \end{minipage} \end{table} \begin{table}[h!] \caption{Cross-correlations between $cons_{m,t}$ and $DDI_{m,t+k}$ (consumption)} \tiny {3pt} \begin{tabular}{l*{25}{S[table-format=1.3]}} \toprule & \multicolumn{12}{c}{$k<0$} & \multicolumn{1}{c}{$k=0$} & \multicolumn{12}{c}{$k>0$} \\ \cmidrule(lr){2-13} \cmidrule(lr){14-14} \cmidrule(lr){15-26} Metal & {$-12$} & {$-11$} & {$-10$} & {$-9$} & {$-8$} & {$-7$} & {$-6$} & {$-5$} & {$-4$} & {$-3$} & {$-2$} & {$-1$} & {$0$} & {$+1$} & {$+2$} & {$+3$} & {$+4$} & {$+5$} & {$+6$} & {$+7$} & {$+8$} & {$+9$} & {$+10$} & {$+11$} & {$+12$} \\ \midrule Copper & 0.208 & 0.233 & 0.235 & 0.252 & 0.240 & 0.251 & 0.236 & 0.309 & 0.296 & { 0.311} & 0.304 & 0.262 & 0.216 & 0.188 & 0.184 & 0.167 & 0.181 & 0.155 & 0.129 & 0.140 & 0.098 & 0.056 & 0.067 & -0.011 & -0.022 \\ Nickel & -0.067 & -0.077 & -0.082 & -0.055 & -0.002 & 0.040 & 0.080 & 0.113 & 0.125 & 0.141 & 0.160 & 0.151 & 0.163 & 0.163 & 0.178 & 0.203 & 0.233 & 0.260 & 0.272 & 0.272 & { 0.282} & 0.259 & 0.251 & 0.215 & 0.155 \\ \bottomrule \end{tabular} \begin{minipage}{0.9\textwidth} \scriptsizeNotes: The table reports cross-correlations between the cyclical component of consumption ($c_{m,t}$) and $DDI_{m,t+k}$ for $k=-12,\ldots,+12$. Negative (positive) $k$ indicates that DDI leads (lags) consumption. Bold values denote peak correlations for each metal. \end{minipage} \end{table} \begin{table}[h!] \caption{Cross-correlations between $DDI_{t+k}$ and copper: prices vs consumption.} {4pt} \begin{tabular}{l*{13}{S[table-format=1.3]}} \toprule & \multicolumn{6}{c}{$k<0$} & \multicolumn{1}{c}{$k=0$} & \multicolumn{6}{c}{$k>0$} \\ \cmidrule(lr){2-7} \cmidrule(lr){8-8} \cmidrule(lr){9-14} Series & {$-6$} & {$-5$} & {$-4$} & {$-3$} & {$-2$} & {$-1$} & {$0$} & {$+1$} & {$+2$} & {$+3$} & {$+4$} & {$+5$} & {$+6$} \\ \midrule Price & 0.539 & 0.549 & 0.553 & { 0.558} & 0.554 & 0.531 & 0.496 & 0.455 & 0.424 & 0.408 & 0.398 & 0.387 & 0.364 \\ Consumption & 0.236 & 0.309 & 0.296 & { 0.311} & 0.304 & 0.262 & 0.216 & 0.188 & 0.184 & 0.167 & 0.181 & 0.155 & 0.129 \\ \bottomrule \end{tabular} \begin{minipage}{0.9\textwidth} \scriptsize\textit{Notes:} The table reports cross-correlations between $DDI_{t+k}$ and the cyclical component of copper prices ($c_t$) and consumption ($cons_t$). Negative (positive) $k$ indicates that DDI leads (lags) the variable. Bold values denote peak correlations within each row. \end{minipage} \end{table} \begin{table}[h!] \caption{Cross-correlations between mean $DDI_{m,t+k}$ and mean of the cyclical component of mineral prices ($\text{for} \ m \in \{\text{CO},\text{CU},\text{NI}\}$).} {4pt} \begin{tabular}{l*{13}{S[table-format=1.3]}} \toprule & \multicolumn{6}{c}{$k<0$} & \multicolumn{1}{c}{$k=0$} & \multicolumn{6}{c}{$k>0$} \\ \cmidrule(lr){2-7} \cmidrule(lr){8-8} \cmidrule(lr){9-14} Series & {$-6$} & {$-5$} & {$-4$} & {$-3$} & {$-2$} & {$-1$} & {$0$} & {$+1$} & {$+2$} & {$+3$} & {$+4$} & {$+5$} & {$+6$} \\ \midrule Price & 0.714 & 0.738 & 0.752 & { 0.752} & 0.748 & 0.732 & 0.718 & 0.692 & 0.664 & 0.631 & 0.608 & 0.590 & 0.573 \\ \bottomrule \end{tabular} \begin{minipage}{0.9\textwidth} \scriptsize\textit{Notes:} The table reports cross-correlations between mean $DDI_{m,t+k}$ and mean of the cyclical component of mineral prices ($\text{for} \ m \in \{\text{CO},\text{CU},\text{NI}\}$). Negative (positive) $k$ indicates that DDI leads (lags) the variable. Bold values denote peak correlations within the row. \end{minipage} \end{table} \begin{table}[h!] \caption{Cross-correlations between mean $DDI_{m,t+k}$ and mean of the cyclical component of mineral prices and consumptions ($\text{for} \ m \in \{\text{CU},\text{NI}\}$).} {4pt} \begin{tabular}{l*{13}{S[table-format=1.3]}} \toprule & \multicolumn{6}{c}{$k<0$} & \multicolumn{1}{c}{$k=0$} & \multicolumn{6}{c}{$k>0$} \\ \cmidrule(lr){2-7} \cmidrule(lr){8-8} \cmidrule(lr){9-14} Series & {$-6$} & {$-5$} & {$-4$} & {$-3$} & {$-2$} & {$-1$} & {$0$} & {$+1$} & {$+2$} & {$+3$} & {$+4$} & {$+5$} & {$+6$} \\ \midrule Price & 0.610 & 0.640 & 0.660 & 0.673 & { 0.680} & 0.668 & 0.648 & 0.619 & 0.602 & 0.578 & 0.564 & 0.552 & 0.532 \\ Consumption & 0.181 & 0.243 & 0.242 & 0.240 & { 0.252} & 0.204 & 0.168 & 0.162 & 0.142 & 0.099 & 0.126 & 0.113 & 0.096 \\ \bottomrule \end{tabular} \begin{minipage}{0.9\textwidth} \scriptsize\textit{Notes:} The table reports cross-correlations between mean $DDI_{m,t+k}$ and mean of the cyclical component of mineral prices and consumptions ($\text{for} \ m \in \{\text{CU},\text{NI}\}$). Negative (positive) $k$ indicates that DDI leads (lags) the variable. Bold values denote peak correlations within each row. \end{minipage} \end{table}

Structural VAR models for global metal markets

We incorporate the DDIs into metal-specific SVAR models to identify supply- and demand-side drivers of the real price of each commodity:

equation[equation omitted — 312 chars of source]

where $n$ denotes the number of endogenous variables, $\mathbf{A}_{m,0}$ captures the contemporaneous relations among the model variables, while $\mathbf{c}_m$ and $\mathbf{A}_{m,j}$ are the vector of intercepts and the matrix of coefficients, respectively.

The vector of endogenous variables is defined as $\mathbf{y}_{m,t} = \left[\Delta prod_{m,t}^o, \Delta WIP_t,p_{m,t},DDI_{m,t} \right]^\prime$ for $m \in \{\text{CO},\text{CU},\text{NI}\}$, where $\Delta prod_{m,t}^o$ denotes the percentage growth rate of global metal ore production\footnote{As robustness checks, the production of refined metals is also considered. For instance, focusing on nickel, the distinction between ore and refined production reflects the difference between raw material extraction and processed output, with ore data capturing gross volumes of mined material and refined data representing the nickel content available for industrial use, thereby offering a more accurate measure of supply relevant to downstream economic activity.}, which is sourced from the WBMS Database. For cobalt, we use the sum of ore production for nickel and copper given that 99% of cobalt is mined as a byproduct of these metals cobaltinstitute2024cobalt. $\Delta WIP_t$ is the percentage growth rate of World Industrial Production BHaer. $p_{m,t}$ is the natural logarithm of the global real price of commodity $m$, obtained by deflating the nominal price using the US CPI from the Federal Reserve Economic Data (FRED), and expressed in constant 2023 US dollars. Specifically, for cobalt we use the London Metal Exchange (LME) spot price, whereas for copper and nickel we use the LME spot price adjusted for delivery to European ports (CIF). All price series are sourced from the International Monetary Fund (IMF) website. Finally, $DDI_{m,t}$ denotes the derived demand index introduced in Section (ref).

In the empirical analysis, we estimate all models with monthly data over the period from December 2006 to December 2023 yielding 205 observations. The VAR lag order is set to $p=3$ (other robustness checks using $p=6$ and $p=12$ yield qualitatively similar structural impulse responses). All the models are estimated using the Bayesian methodology detailed in AriasRubioWaggoner2018 and imposing a diffuse Normal–Inverse-Wishart prior for the reduced-form parameters. For further details about the estimation procedure we refer the reader to the Supplement (available from the authors upon request).

commentusing a flat Normal-Inverse-Wishart prior on the reduced-form parameters and identify structural shocks following the algorithm of .\footnote{Specifically, let $\mathbf{B}_m = \mathbf{A}_{m,0}^{-1}[\mathbf{c}_m, \mathbf{A}_{m,1}, \ldots, \mathbf{A}_{m,p}]$ denote the reduced-form counterpart of equation (ref) and $\boldsymbol{\Sigma}_m = (\mathbf{A}_{m,0}\mathbf{A}_{m,0}')^{-1}$ the reduced-form covariance matrix. We assume $p(\mathbf{B}_m, \boldsymbol{\Sigma}_m) \propto |\boldsymbol{\Sigma}_m|^{-(n+1)/2}$, which corresponds to a Normal-Inverse-Wishart prior with prior precision $N_0 = 0$ and degrees of freedom $v_0 = 0$. The posterior is then also Normal-Inverse-Wishart, with posterior precision $N_T = X_m'X_m$, posterior mean $\bar{\mathbf{B}}_{m,T} = \mathbf{B}_{m,OLS}$, posterior degrees of freedom $v_T = T$, and posterior scale matrix $S_T = \hat{\mathbf{U}}_m'\hat{\mathbf{U}}_m$, where $\hat{\mathbf{U}}_m$ are the OLS residuals. Further details are available in the Supplement.}

Identification

Starting from Equation (ref), we identify four structural shocks by combining: (i) a set of contemporaneous sign and zero restrictions imposed on the matrix of impact responses $\mathbf{A}_{m,0}^{-1}$; and (ii) restrictions on the magnitude of the response of real economic activity to what we label transition and aggregate demand shocks.

{Contemporaneous zero and sign restrictions.} As shown in Table (ref), we identify four structural shocks by enforcing a set of contemporaneous sign and zero restrictions. We normalize shocks so that they are all associated with an increase in the real price of metals (i.e., a negative supply and positive demand shocks).

table[table omitted — 871 chars of source]
itemize• Supply shock. A negative supply shock represents an exogenous disruption to the production process -- such as labor strikes or natural disasters affecting major mining operations -- that leads to a contraction in global metal production within the same month. This reduction in supply exerts upward pressure on real metal prices through reduced availability, as commonly assumed in the literature BRT26Copper,BaumeisterSpecialFocus,BamueisterPeersman2013,KilianMurphy2012. Importantly, we assume that the within-month response of $DDI_{m,t}$ is constrained to be zero. This assumption is motivated by two considerations: either by appealing to the notion that inventories smooth short-term fluctuations in final goods demand, or by considering that consumers’ intention to purchase products included in the DDI is unaffected on impact. • Aggregate demand (AD) shock. An AD shock is associated with fluctuations in the global business cycle and the resulting changes on the demand for industrial commodities. We follow the common assumption that an unanticipated aggregate demand expansion leads to an increase in global metal production, global economic activity, and the real price of metals BaumeisterSpecialFocus,BHaer,boer2024energy,KilianMurphy2012. In this case, the impact response of $DDI_{m,t}$ is left unrestricted. • Transition demand (TD) shock. This shock captures unexpected increases in demand for metals driven by their role as inputs in selected final goods associated with the energy and digital transition -- such as batteries, electric vehicles, solar panels, and smart devices. Conceptually, it represents a metal-specific demand shock arising from surges in demand for these goods. We assume that a positive TD shock leads to an immediate increase in real metal prices and $DDI_{m,t}$, while remaining agnostic about its contemporaneous effects on production and economic activity. • Metal-specific demand (MD) \textit{shock}. A MD shock captures unexpected changes in the demand for a specific metal that are not directly linked to the global business cycle or to the energy and digital transition. These shocks are a residual category encompassing innovations, such as precautionary demand shocks (e.g., anticipation of future supply disruptions), expectation shocks related to shifts in uncertainty, and unforeseen increases in demand arising from uses of the metal not captured by our DDIs. These shocks are assumed to have no contemporaneous impact on $DDI_{m,t}$, while generating an immediate increase in real metal prices and global production, and simultaneously inducing a contraction in global economic activity due to the cost-push effects of higher input prices boer2024energy,BaumeisterSpecialFocus,KilianMurphy2012.

{Dynamic “magnitude” restrictions.} Table (ref) shows that the identifying restrictions do not rule out cases in which the impact responses associated with AD and TD shocks share the same sign (see Column 2 and 3). In such situations, the identification of AD and TD shocks would fail. To avoid any potential confounding between AD and TD shocks, we impose an additional restriction: the response of real economic activity to an AD shock must be larger in magnitude than the response to a TD shock, both on impact and during the following three months. The idea of restricting the magnitude of a variable’s response to different shocks is also employed by cross2022role to disentangle precautionary and speculative oil demand shocks.\footnote{Related approaches include piffer, who propose restricting the shape of structural impulse responses, and amir2021identification, who develop a framework for enforcing ranking restrictions.}

Results

In this section, we report impulse response function, forecast error variance decomposition, and historical decomposition for the shocks identified in Section (ref).

Structural impulse response analysis

We compute structural impulse-response functions (IRFs) from 2,500 posterior draws that satisfy the identification restrictions described in Section (ref). We report impulse responses over a 24-month horizon to a one-standard-deviation shock and normalize the sign so that all shocks generate a positive price response (i.e., we consider a negative supply shock and positive demand shocks). We follow the Bayesian decision-theoretic approach of INOUE2022457, which conducts inference on the entire vector of impulse responses and therefore accounts for dependence across horizons and variables.\footnote{This differs from the conventional approach based on pointwise posterior medians and marginal credible intervals, which ignore this dependence and may imply response paths that are inconsistent with any feasible parameterization of the underlying VAR model. Our implementation uses an angular loss function, which emphasizes similarity in dynamic response paths and is invariant to the scaling of both the data and the normalization of structural shocks.} The Bayes estimator corresponds to the impulse-response vector that minimizes expected loss under the joint posterior distribution. Estimation uncertainty is summarized using joint credible sets defined as the subset of posterior draws with the lowest posterior risk whose cumulative probability equals the desired credibility level. Figures (ref)--(ref) present, for each metal, the impulse responses of ore production, WIP, and the real price to the four structural shocks. IRFs for variables entering the model in first differences are cumulated. Each panel reports the Bayes estimate (thick blue line), the posterior trajectories in the 68% joint credible set (thin gray lines), and, for comparison, conventional 68% pointwise credible intervals (dashed lines).

\noindentSupply and AD shocks. A natural starting point is the response to supply and aggregate demand (AD) shocks, whose properties are well-established in the commodity market literature. AD shocks generate the largest impact responses and remain among the most persistent shocks, with impact price responses (based on the Bayes estimate) of 6.50% for cobalt, 3.67% for copper, and 7.05% for nickel, reaching peak responses of 8.29%, 4.67%, and 9.29% within the first 4 months, in all cases. Consistent with their business-cycle interpretation, AD shocks generate a simultaneous expansion in both ore production and World IP. Supply shocks generate the second largest impact responses, with price responses of 5.54% for cobalt, 3.36% for copper, and 2.55% for nickel, accompanied by production contractions and negative World IP responses across all metals. Both patterns are consistent with results for the oil market BHaer and for metals bauer2024critical,boer2024energy.

figure[figure omitted — 836 chars of source]
figure[figure omitted — 356 chars of source]
figure[figure omitted — 355 chars of source]

\noindentTD shocks. TD shocks exhibit a pattern that differs from the other shocks in the model, and that is directly informative about the macroeconomic nature of the twin transition. Specifically, TD shocks have two distinctive features. First, they generate persistent and economically significant price effects that build up over time. On impact, price responses are moderate: 0.89% for cobalt, 2.19% for copper, and 1.45% for nickel. However, the response to TD shocks strengthens after impact: peak responses are substantially larger than the corresponding impact responses and are reached only after several months. In particular, the peaks are at 3.24% for copper at month 3 and 5.23% for nickel at month 6 later than the peaks associated with the other shocks. The ratio of peak to impact response is particularly striking for nickel, at 3.6, compared to 1.3 for AD shocks and 1.3 for supply shocks. For copper and nickel, TD impulse responses decline more slowly than those of supply shocks, despite their smaller initial effect. A different story arises for cobalt, where the evidence for a sustained price build-up following TD shocks is considerably weaker than for copper and nickel.

Second, and in sharp contrast to their price effects, TD shocks have mostly little effect on metal production or on WIP. Impact responses for ore production are 0.70%, 0.54%, and 0.80% for cobalt, copper, and nickel -- well below those of supply and AD shocks -- while WIP responses of -0.07%, 0.34%, and 0.19% are negligible and their credible intervals always include zero. These findings are consistent with the interpretation of TD shocks. Demand generated by the energy and digital transition for specific metal-intensive final goods -- batteries, electric vehicles, solar panels -- that is large enough to move metal prices persistently but does not constitute a sufficient share of global output to shift aggregate real activity on impact.

\noindentMD shocks. MD shocks are characterized by an absence of price amplification beyond impact: peak price responses occur at month 1 for all three metals, implying that the entire effect materializes on impact and then quickly disappears. Impact price responses are smaller than those of AD and supply shocks and, for copper and nickel, smaller than the peak TD responses documented above. This pattern is particularly notable given that MD shocks are associated with large production responses, combined with WIP contractions, consistent with a cost-push channel. The absence of persistent price dynamics distinguishes MD shocks sharply from TD shocks, which build up over several months, and confirms that the decomposition captures economically distinct phenomena rather than observationally equivalent demand shocks.

\noindentCross-metal heterogeneity. Two patterns of cross-metal heterogeneity are worth highlighting. The first concerns cobalt. Unlike copper and nickel, the cobalt TD shock generates a flat and uncertain response: the Bayes estimate peaks on impact and the joint credible set remains wide at all horizons. One possible explanation is cobalt's byproduct nature -- since no mine is dedicated exclusively to cobalt extraction, its production is tied to that of nickel and copper, which may make supply less responsive to metal-specific demand shifts and amplify the role of supply-side noise in the price signal. The second concerns the speed of price adjustment under TD shocks. The peak response for nickel occurs at month 6, twice as late as for copper (month 3), and the price profile remains elevated throughout the 24-month horizon. Given nickel's central role in battery cathodes for electric vehicles, this slower and more persistent adjustment is consistent with longer investment cycles and relatively inelastic supply in refined nickel markets. Further evidence on the relative contribution of each shock to price variability across horizons is provided in the next section through a forecast error variance decomposition.

comment\subsection{Structural impulse response analysis} \begin{figure}[t] \caption{Response of real prices to structural shocks.} \caption*{\scriptsize Notes: Impulse responses of cobalt, copper, and nickel real price to the four shocks. The plots are based on 2500 draws and show the 68% pointwise credible sets (light blue dashed lines), the Bayes estimators (blue lines) and the 68% joint credible sets (grey lines) under an angular loss function. The size of all shocks is one standard deviation. To express the responses in percentage points, all lines are multiplied by 100.} \end{figure}

Forecast Error Variance Decomposition

Table (ref) reports the forecast error variance decomposition (FEVD) of real metal prices across structural shocks, metals, and horizons. To examine how the importance of each shock has evolved as the twin transition has become increasingly important, we estimate the model separately over two subsamples of equal length: Subsample 1 (12/2006--05/2015) and Subsample 2 (06/2015--12/2023).\footnote{Each entry in Table (ref) represents the median contribution of each shock to the variance of the real price at a given horizon, normalized so that rows sum to one. Normally, the median across different draws would not sum to 1 because the median is non-additive. To interpret each element of the table as a percentage, we divide each entry by the sum of its row and multiply by 100.}

Consistent with the IRF analysis, AD shocks are the main driver of real metal prices, explaining more than half of total variance in most cases. Their dominance reflects major macroeconomic downturns in each subsample: the Global Financial Crisis in Subsample 1 and the COVID-19 pandemic in Subsample 2. Supply shocks account for about one-third of price variance on impact across metals and subsamples, but their contribution typically declines over longer horizons.

The FEVD also reveals two complementary patterns. First, for all metals and in both subsamples, the contribution of TD shocks increases with the forecast horizon. The same pattern is observed in the full sample case reported in the Supplement (available from the authors upon request). Second, for copper and nickel, the share of real price variance explained by TD shocks rises in the more recent subsample, consistent with the growing macroeconomic relevance of the twin transition. This effect is particularly pronounced for copper: in Subsample 2, TD shocks become the dominant driver of price variance at all horizons beyond impact, with their contribution increasing from 12.42% on impact to 59.52% at the 120-month horizon, compared to 25.34% in Subsample 1 at the same horizon.

Across metals and subsamples, MD shocks consistently account for the smallest share of price variance, rarely exceeding 10% at any horizon, with their contribution largest on impact. This stability across subsamples is itself informative: it implies that the redistribution of variance from Subsample 1 to Subsample 2 -- specifically, the rise of TD shocks at the expense of AD and supply shocks -- is not driven by residual demand variation captured by MD shocks. The quantitative importance of each shock across specific historical episodes is examined in the following subsection through historical decompositions.

table[table omitted — 2,023 chars of source]

Historical Decompositions

We rely on historical decompositions to attribute observed price movements to specific structural shocks over time. We focus on two price surge episodes shared across all three metal markets: the post--Global Financial Crisis period (January--June 2009) and the post-pandemic surge (December 2021--June 2022), which marked a turning point driven by rapid electric vehicle adoption and the accelerated deployment of renewable energy technologies IEAbattery,iea_ev_sales_2025.\footnote{According to iea_ev_sales_2025, total sales of EVs (including both plug-in hybrid electric vehicles and battery electric vehicles) more than tripled between 2020 and 2022, from around 3 million to 10.1 million units.}

The left panels of Figure (ref) show that the 2009 recovery was largely driven by AD shocks (light blue bar) across all three metals after March 2009, with supply shocks (blue bar) contributing to upward pressure for copper and nickel. TD shocks (orange bar) are generally small and in some cases negative; TD shares of total explained variation are 7.8% for copper and 9.3% for nickel.\footnote{To quantify the relative role of TD shocks in each episode, we compute their share of total explained price variation as $\sum_{t \in \mathcal{T}} |\hat{\psi}_{\mathrm{TD},t}| \,/\, \sum_{j} \sum_{t \in \mathcal{T}} |\hat{\psi}_{j,t}|$, where $\hat{\psi}_{j,t}$ is the contribution of shock $j$ at time $t$, $\mathcal{T}$ denotes the episode window, and absolute values prevent positive and negative contributions from canceling.}

figure[figure omitted — 201 chars of source]

Turning to the right panels, AD shocks remain the primary driver of the post-pandemic price surge for cobalt and nickel, while supply shocks also contribute positively after March 2022. For nickel, MD shocks (green bar) reflect heightened precautionary demand linked to Russia’s invasion of Ukraine, which triggered fears of supply losses from a key producer of high-grade nickel. Consistent with the definition of MD shocks, these concerns, rather than realized shortages, contributed to the surge in real nickel prices from March 2022 onwards, as market participants revised expectations and engaged in anticipatory stockpiling ReutersNickel.

Finally, TD shocks (orange bar) in the most recent subsample play a quantitatively significant important role for copper and, to a lesser extent, nickel, consistent with rapid EV adoption and the accelerated deployment of renewable energy technologies IEAbattery,iea_ev_sales_2025. The TD share for nickel rises from 9.3% to 14.2%, while for copper the increase is more pronounced, from 7.8% to 57%, making TD shocks the most important driver of copper price movements in the post-pandemic episode.

commentThis robustness check examines the effects of including in the model the percentage growth rate of inventories ($\Delta inv_{m,t}$), measured with LME Warehouse Stocks sourced from Refinitiv. In this case, a new shock—referred to as the Expectational Demand (ED) shock—is identified. All the restrictions detailed in Section (ref), including the magnitude, remain unchanged. The additional sign and zero restrictions related to inventories and the ED shock are summarized in Table (ref). In particular, a positive ED shock is assumed to share the same signs as the MD shock, except that it causes inventories to increase. The rationale is that an ED shock reflects precautionary inventory accumulation reflecting, for instance, the anticipation of future supply disruptions. This assumption follows an established approach in the literature and has been applied both to the oil market (KilianMurphy2012) and to metal markets (boer2024energy,BaumeisterSpecialFocus). \begin{table}[t] \caption{Sign and zero restrictions on the impact response matrix $\mathbf{A}_{m,0}^{-1}$} \begin{tabular}{c|ccccc} \hline & \multicolumn{5}{c}{shock} \\\cline{2-6} & (1) & (2) & (3) & (4) & (5) \\ variable & Supply shock & AD shock & ED shock & TD shock & MD shock \\\hline $\Delta prod_{m,t}^o$ & $-$ & $+$ & $+$ & $\ast$ & $+$ \\ $\Delta WIP_{t}$ & $-$ & $+$ & $-$ & $\ast$ & $-$ \\ $\Delta inv_{m,t}$ & $\ast$ & $\ast$ & $+$ & $\ast$ & $-$ \\ $p_{m,t}$ & $+$ & $+$ & $+$ & $+$ & $+$ \\ $DDI_{m,t}$ & $0$ & $\ast$ & $0$ & $+$ & $0$ \\\hline \end{tabular} \caption*{\scriptsizeNotes: “$\ast$” indicates that the element is left unrestricted. The signs follow the convention that all shocks are normalized to generate an increase in the real price of metals (i.e., a negative supply shock and positive demand shocks).} \end{table} \begin{itemize} • Lag order (6 (or 12?) instead of 3) • Minnesota prior instead of flat prior • Refined production instead of ores production (a part from cobalt) • Model that treats COVID in some way • Model that treats seasonality in some way • REA instead of WIP • Model estimated with external shock (Baumeister method) • PCA estimated on 15 GT series instead of 10 to address estimation uncertainty • For Haar prior problem: plot the IRF without the data and check if IRFs are similar. If so, the problem is significant. \end{itemize} We estimate an IV Local Projection using the TD shock of a metal as instrument for the respective price. Specifically we estimate: \begin{equation} \Delta \log(y_{t+h}) = \alpha^h + \beta^h p_{m,t} + \sum_{l=1}^{L} \phi^h_{x,l} X_{t-l} + \varepsilon_{t+h}, \quad for h = 0,\dots,24 . \end{equation} where $y_{t+h}$ is eihter: US consumer price index, US producer price index, US industrial production or EU consumer price index at time t + h, $\alpha^h$ is the constant, $p_{m,t}$ is the log of real copper price at time t (either in log levels or log first difference).The set of controls $X_{t-l}$ includes L = 12 lags of $y_t$ and $p_{m,t}$, as well as a global economic activity index from BHaer.

Conclusions

In this paper, we construct a novel monthly metal-specific derived demand indexes for cobalt, copper, and nickel using Google Trends web search data, and incorporate them into metal-specific SVAR models to identify supply- and demand-side shocks in global metal markets. Our key contribution is the identification of transition demand (TD) shocks associated with the demand for technologies crucial for the energy and digital transitions.

The SVAR model analysis shows that TD shocks generate persistent price effects that build gradually and decay slowly, especially for copper and nickel, whereas supply and metal-specific demand shocks tend to have more immediate and less persistent effects. Forecast error variance decompositions and historical decompositions analyses show that TD shocks have become more relevant in recent years, especially for copper and nickel, highlighting the growing macroeconomic relevance of transition-related demand.

Emerging sources of demand, including artificial intelligence and data centre infrastructure, may further strengthen the role of metals such as copper. As these forces are only partially captured in our indexes and have become more relevant in recent years, our estimates may understate the broader contribution of transition-related demand.

Overall, not all demand-driven price increases are alike: distinguishing among these shocks is essential for understanding commodity price dynamics in an economy undergoing rapid energy and digital transitions.