1 Introduction

The Global Financial Crisis (GFC) and the COVID-19 pandemic have drawn renewed attention to shock transmissions in the economy across both production and financial networks. The economics literature has shown the possibility of “cascade effects”, whereby sectoral shocks propagate to direct trading partners and also to the aggregate economy. Acemoglu et al. (2012) show that input–output connections between sectors act as a potential propagation mechanism of idiosyncratic shocks throughout the entire economy.Footnote 1 Carvalho (2014) demonstrates the importance of the structure of a production network: in an economy characterized by a small number of production hubs supplying inputs to many different firms or sectors, microeconomic shocks tend to result in higher aggregate volatility. Moreover, Atalay (2017) estimates that sectoral shocks account for nearly two thirds of the volatility of aggregate output.Footnote 2

In the same vein, the finance literature has provided ample evidence of the connectedness between financial markets. The co-movements of stock returns and volatility between international stock markets have been documented by Hamao et al. (1990), King and Wadhwani (1990), Diebold and Yilmaz (2009), Zhu et al. (2018), Jawadi et al. (2019), inter alia. Other studies examine shock spillovers among related industries within specific sectors such as energy (Alli et al., 1994; Ewing et al., 2002) and financials (Elyasiani et al., 2007). Wang (2010) investigates the dynamics and causality among 30 US industries to identify leading-lagging industries. Neither of these studies, however, consider the underlying economic relationships as factors potentially driving the uncovered spillovers.

Our paper bridges these two lines of research. We quantify the shock spillovers of stock returns among all the economic sectors in the US economy. In this respect, we generalize—to the best of our knowledge, for the first time—the specific, intra-sectoral perspective so far offered by the financial studies. We then investigate how the intensity of the said financial spillovers depends on the strength of the intersectoral trade relationships, as captured by the set of input–output connections among all the sectors.

The first part of our analysis documents the existence and the intensity of financial shock spillovers among all the US non-governmental sectors. We employ a technique proposed by Diebold and Yilmaz (2012, 2014) to analyse the daily stock returns of 14 US economic sectors. The approach, which decomposes the forecast error variance under a generalized vector autoregressive (VAR) framework, allows us to evaluate the degree to which shocks to stock returns are transmitted: (1) between pairs of sectors; (2) from a sector to the rest of the market; and (3) from the rest of the market to a sector. The technique can also be used to compute the extent of total spillovers, which measures the overall degree of connectedness among all the sectors in the economy.Footnote 3 We augment the Diebold and Yilmaz (2012, 2014) approach with exogeneous variables and estimate a VARX system which accounts for macroeconomic factors, stripping the spillover measures from the common macroeconomic effects which are likely to affect all the sectors in the economy.Footnote 4 To gauge how shock spillovers propagate across sectors and evolve over time, we conduct both a static, full sample, and a dynamic, rolling-sample, analysis.

The second important contribution of this paper is to investigate the relationship between the intensity of shock spillovers, uncovered in the first stage, and the strength of the trade links between the sectors. We are motivated by recent empirical evidence which suggests that economic links explain contemporaneous correlations between firms’ stock returns and at the same time help predict future stock returns (Cohen & Frazzini, 2008; Menzly & Ozbas, 2010). To quantify the strength of the economic link between sectors, we follow Becker and Thomas (2011) and Ahern and Harford (2014) by using information from the Input–Output (IO) tables published by the US Bureau of Economic Analysis (BEA).Footnote 5

As a preview of our variance decomposition results, we observe significant shock transmission between US sectors. The total shock spillover among the 14 sectors in our sample is 85.8%. Shocks to a sector are found to account for up to 10.16% of the forecast error variances of its trading partners’ stock returns. Compared to the average sector’s own contribution of 14.2%, this cross variance share is considerably large. All shock spillover measures remain at high levels when the spillovers are measured both statically, using the full sample period, and dynamically, over the 200-day rolling windows.

The second part of our analysis confirms that business linkages have significant effects on shock transmissions between sectors: the closer the trading relationship in a sector pair, the stronger the degree of spillovers. Shocks to a sector’s important supplier account for a large fraction of the forecast error variance of the sector’s return. Additionally, we show that the total directional spillovers from and to a sector correlate with the number of close linkages between that sector and the others: a sector which is a main supplier/customer to many other sectors is likely to be a net shock transmitter in the network.

An interesting additional result is that the intersectoral shock transmission and its dependence on the trade links are affected by the market conditions. In particular, our analysis conducted on sub-samples representing different market conditions shows stronger shock transmissions between sectors during periods of stress such as the GFC and the COVID-19 pandemic.

We check for the robustness of our results by subjecting them to a series of sensitivity tests. Firstly, we find similar evidence of significant shock spillovers among sectors when we conduct our analysis on a sample including only non-financial sectors. Secondly, we show that the magnitude of the shock spillovers is insensitive to the order of the VARX model and the choice of the forecast horizon. Finally, the relationship between shock spillovers and the strength of the business linkages between sectors is robust to an alternative construction of the business linkage variables.

Our results provide important implications for investors and portfolio managers with exposure to specific sectors. Sector funds, such as agriculture or energy funds, are attractive to investors who want to focus their knowledge on a manageable set of companies with similar business operations. Several studies have reported the outperformance of sector fund strategies. Chen and Hackbarth (2020) use the data of active equity sector funds and show that a naïve equal-weighted portfolio provides reliable alpha and the outperformance is resilient in market downturns, with infrequent need of rebalancing. O’Neal (2000) finds that sector mutual funds have profitable trading strategies based on sector momentum effect. However, these sector-specific investment strategies entail greater total and systematic risk than index. Thus, for better risk management, our study suggests that practitioners should pay close attention to sectors with close business linkages because returns and volatility of their investments could be affected indirectly by shocks to the related sectors.

Understanding shock transmissions across sectors also has important policy implications. A policy change targeting a particular sector could potentially have considerable impacts on other sectors given their high degree of connectedness. Moreover, since our results show that sectors in the US economy are closely connected, large systemic idiosyncratic shocks hitting several sectors simultaneously could potentially have a devastating negative impact on the financial system especially during turbulent periods.

Our study also contributes to the recent literature strand on the financial impacts of COVID-19 (see Corbet et al., 2020; Iqbal et al., 2021, among others). Akhtaruzzaman et al. (2021) report a significant increase in volatility and shock transmissions between stock markets during the COVID-19 outbreak. Sharif et al. (2020) employ a wavelet coherence approach and show that the pandemic is expected to have a long-term negative impact on the US economic uncertainty. A similar approach is used by Choi (2020) to investigate the effect of the COVID-19 on US sectors. These studies reveal that the economic policy uncertainty induced by COVID-19 has greater influence on all sectors volatility, compared to the effect of GFC. However, the impact of COVID-19 on the shock spillovers between sectors has not been examined. Thus, our work provides insights into sectoral shock transmissions during the pandemic, which are valuable for policy makers, regulators and practitioners.

Our paper is structured in three main parts. Part I investigates the degree of connectedness among the US economic sectors. Section 2 describes the generalized variance decomposition methodology, Sect. 3 discusses the sector stock return data and summary statistics, and Sect. 5 reports the empirical results of the shock transmissions between sectors. Part II extends the analysis to the relation between inter-sectoral shock spillovers and the strength of their business linkages (Sect. 5). Part III reports additional results. Section 6 examines the shock transmissions and the impact of business linkages in different market conditions. Section 7 collects our robustness checks. Finally, Sect. 8 concludes.

Part I: Shock Spillovers among the US Sectors

2 Methodology

To investigate the transmission of shocks, we follow Diebold and Yilmaz (2012, 2014) to use a generalized vector autoregression (VAR) system for \(N\) sectors in the economy, and measure the forecast error variance decompositions across sectors. To account for the common impact of macroeconomic factors on sector returns, we extend VAR to include exogenous variables (see also Alter & Beyer, 2014; Claeys & Vašíček, 2014). Specifically, consider a covariance stationary VARX(p) system of \(N\) economic sectors:

$$ \user2{y}_{t} = \sum\limits_{{k = 1}}^{p} {{\mathbf{\Phi }}_{k} } \user2{y}_{{t - k}} + \user2{b}^{\prime } \user2{x}_{t} + \user2{\varepsilon }_{t}, $$
(1)

where \({{\varvec{y}}}_{t}=\left({y}_{1t},{y}_{2t},\dots ,{y}_{Nt}\right){^{\prime}}\) is an \(N\times 1\) vector of sector excess returns at time \(t\), \({{\varvec{\varepsilon}}}_{t}\sim \left(0,{\varvec{\Sigma}}\right)\) is a vector of the error terms at time \(t\), which are independently and identically distributed with zero mean and variance–covariance matrix \({\varvec{\Sigma}}\), and \({{\varvec{x}}}_{t}\) is a vector of exogenous macroeconomic variables at time t. The corresponding moving average representation of the VAR(p) system in Eq. (1) is

$${{\varvec{y}}}_{t}=\sum_{k=0}^{\infty }{{\varvec{A}}}_{k}{{\varvec{\varepsilon}}}_{t-k},$$
(2)

where the \(N\times N\) coefficient matrix \({A}_{k}\) is calculated recursively:

\({{\varvec{A}}}_{k}={{\varvec{\Phi}}}_{1}{{\varvec{A}}}_{k-1}+{{\varvec{\Phi}}}_{2}{{\varvec{A}}}_{k-2}+\cdots +{{\varvec{\Phi}}}_{p}{{\varvec{A}}}_{k-p}\) for \(k>0\), \({{\varvec{A}}}_{k}\) is an identity matrix for \(k=0\), and \({{\varvec{A}}}_{k}=0\) for \(k<0\).

The moving average coefficients, known as variance decompositions, reveal the dynamics of the system. Variance decompositions show the fraction of the \(H\)-step-ahead error variance in forecasting \({y}_{i}\) accounted for by the shocks to \({y}_{j}\), \(\forall j\ne i\), for each \(i\) in the system. Since VAR innovations are generally contemporaneously correlated, the variance decompositions traditionally rely on identification schemes such as those based on Cholesky factorization to achieve orthogonality (Sims, 1980). Using the traditional approach, however, renders variance decompositions dependent on the ordering of the variables in the VAR system. Diebold and Yilmaz (2012, 2014) overcome the ordering-dependence problem by employing the generalized VAR framework proposed by Koop et al. (1996) and Pesaran and Shin (1998). Instead of orthogonalizing the shocks, the framework uses the historically observed distribution of the errors to account for the correlated shocks.

The coefficient matrix \({{\varvec{A}}}_{k}\) contains a large amount of useful information which can be extracted as follows. Firstly, using the elements in \({{\varvec{A}}}_{k}\), we calculate the cross variance shares, or pairwise directional spillovers, which are defined as the proportion of the \(H\)-step-ahead forecast error variance of sector \(i\)’s return attributable to the shocks to sector \(j\) (\(i,j=\mathrm{1,2},\dots ,N\); \(i\ne j\)). When \(i=j\), the measure is termed the own variance shares which are defined as the proportion of the forecast error variance of sector \(i\)’s return accounted for by its own shocks. According to the generalized VAR framework, the \(H\)-step-ahead forecast error variance decomposition is obtained as

$${\theta }_{ij}^{H}=\frac{{\sigma }_{jj}^{-1}\sum_{h=0}^{H-1}{\left({{\varvec{e}}}_{i}^{^{\prime}}{{\varvec{A}}}_{h}{\varvec{\Sigma}}{{\varvec{e}}}_{j}\right)}^{2}}{\sum_{h=0}^{H-1}\left({{\varvec{e}}}_{i}^{^{\prime}}{{\varvec{A}}}_{h}{\varvec{\Sigma}}{{{\varvec{A}}}_{h}^{^{\prime}}{\varvec{e}}}_{i}\right)}$$
(3)

where \({\varvec{\Sigma}}\) is the variance–covariance matrix of the shock vector in the VAR system, \({\sigma }_{jj}\) is the \({j}^{th}\) element of the diagonal of \({\varvec{\Sigma}}\) representing the standard deviation of the shock for the \({j}^{th}\) sector, \({{\varvec{e}}}_{i}\) is the selection vector with one as the \({i}^{th}\) element and zeros elsewhere, and \({{\varvec{A}}}_{h}\) is the coefficient matrix in the infinite moving average representation in Eq. (2), with \(h=0, 1, 2\dots , H-1\). Since shocks to the sectors in the system can be contemporaneously correlated, the sum of all the contributions to a sector’s variance (i.e., the sum of all the elements in each row of the variance decomposition table) can differ from one (i.e., \({\sum }_{j=1}^{N}{\theta }_{ij}^{H}\ne 1)\). For this reason, Diebold and Yilmaz (2012, 2014) normalize entries in the variance decomposition matrix by dividing the quantity in Eq. (3) by the row sum:

$${S}_{i\leftarrow j}^{H}={\widetilde{\theta }}_{ij}^{H}=\frac{{\theta }_{ij}^{H}}{\sum_{j=1}^{N}{\theta }_{ij}^{H}}$$
(4)

Note that \({\sum }_{j=1}^{N}{\widetilde{\theta }}_{ij}^{H}=1\) and \({\sum }_{i,j=1}^{N}{\widetilde{\theta }}_{ij}^{H}=N\) by construction. In our investigation, we refer to \({S}_{i\leftarrow j}^{H}\) as the pairwise spillover from sector \(j\) to sector \(i\) after \(H\) periods. For a system of \(N\) sectors, there are \({N}^{2}-N\) pairwise spillover measures in total.

Secondly, in addition to the measure of pairwise spillovers, we also calculate the total directional spillover from others to sector \(i\) as

(5)

and the total directional spillover to others from sector \(i\) as

(6)

The net directional spillovers from sector \(i\) to all other sectors can then be easily calculated from Eqs. (5) and (6) using

(7)

Finally, the total spillover between all sectors in the system can be obtain as

$${S}^{H}=\frac{1}{N}\sum_{\begin{array}{c}i,j=1\\ i\ne j\end{array}}^{N}{\widetilde{\theta }}_{ij}^{H}$$
(8)

which is the equal-weighted average of the sum of all the entries, excluding the own variance shares, in the variance decomposition matrix.

3 Data and summary statistics

To calculate sector excess returns, we use daily stock returns for all the stocks listed on four US major stock markets—NYSE, Nasdaq, Arca and Amex with share code 10 and 11—obtained from the CRSP database. We sort all stocks into 14 sectors using the NAICS codes. Our dataset spans the period starting from 3rd January 2005 to 31st December 2020, giving us a total of 4028 daily observations for each of the sectors.Footnote 6 Daily sector returns are calculated as the sum of the value-weighted returns of all the stocks in the sector, where the weights are the beginning-of-the-day market capitalization of the stocks. The sector excess return is computed as the difference between the sector return and the yield on the 3-month US Treasury bills, obtained from the Federal Reserve Bank of St. Louis (FRED) database.

To account for macroeconomic factors, we use the short-term risk-free rate, the term risk premium (the difference between the 10-year Treasury bond yield and the 3-month T-Bill yield), the credit risk premium (the difference between Moody’s Baa corporate bond yield and Moody’s Aaa corporate bond yield) and the foreign exchange rate index. We obtain the data on the 3-month T-bill yield, 10-year Treasury yield, Moody’s Baa and Aaa corporate bond yield, and the trade weighted USD indices against a broad group of major US trading partners from the Federal Reserve Bank of St. Louis (FRED) database. The Augmented Dickey-Fuller (ADF) test suggests the use of the first-differences of the yields and the foreign exchange rate index, and confirms the stationarity of all the 14 sector excess return series.

Table 1 presents the descriptive statistics of our data. The mean daily excess returns range from 0.033% (AGR and MNG) to 0.061% (INF). Among the 14 sectors, Mining (MNG) has the largest standard deviations. All the graphs of the sector excess returns in Fig. 1 show a common pattern of a significant increase in volatility following the GFC, especially after the collapse of Lehman Brothers in September 2008, and during the pandemic in 2020.

Table 1 Summary statistics of the sector returns
Fig. 1
figure 1figure 1

Sector stock returns. This figure illustrates the time series of 14 sector excess returns over the period 3 January 2005–31 December 2020

4 Shock spillovers among the US sectors

This section reports the results of the investigation of shock spillovers between sectors using the generalized variance decomposition method. We conduct both a static analysis for the full sample period and a rolling sample (dynamic) analysis of conditional connectedness.

4.1 Static analysis: full-sample spillovers

We first estimate the VARX model for the daily excess returns of a system of 14 sectors using the full sample from 2005 to 2020. Following Diebold and Yilmaz (2012, 2014), the forecast error variance across sectors is decomposed into parts attributable to the various variables in the system. The variance decomposition matrix, also known as the connectedness table, is shown in Table 2. The table reports the magnitude of the pairwise, directional and total spillovers, calculated from the first-order VARX at the 10-day-ahead forecast horizon.Footnote 7

Table 2 Full-sample Shock spillover

The pairwise spillovers between sectors are given in the off-diagonal elements of the \(14\times 14\) matrix shown in Table 2, while elements along the diagonal axis show the contribution of a sector to its own forecast variance. The own connectedness (diagonal elements) tends to account for the largest share of each sector’s forecast error variance, ranging from 10.94% Manufacturing (MFG) to 19.69% Agriculture, forestry, fishing, and hunting (AGR). The average cross variance share (off-diagonal elements) is 6.60%, showing that a sizable proportion of the forecast error variance of sector \(i\) comes from shocks to sector \(j\). The highest pairwise spillover, is 10.16%, from Manufacturing (MFG) to Information (INF). In contrast, the cross variance share of 3.35%, contributed by Mining (MNG) to Retail trade (RT), is the smallest.Footnote 8

The total directional spillovers (from others and to others) tend to be large pointing to strong connections among sectors. The total directional spillovers to a sector from all the other sectors (the row sum excluding the diagonal element) are reported in the rightmost column labelled “FROM others”. The total directional spillovers from a sector to all the others (the column sum) are reported along the bottom row labelled “TO others”.Footnote 9 According to the results presented in Table 2, the total directional spillover coming from other sectors varies between 80.31 and 89.06%, while the total directional spillover originating from a sector to the others ranges from 53.80 to 115.25%. Agriculture, forestry, fishing, and hunting (AGR) is found to be both the smallest shock transmitter and receiver. At the same time, Manufacturing (MFG), Professional and Business Services (PRO), and Wholesale Trade (WST) are the top three shock transmitters and receivers. These three sectors account for both the largest total directional spillovers to others (over 100%) and received from others (nearly 90%).

For each of the 14 sectors, we calculate the net total directional spillovers (NET) as the difference between the total spillover from a sector “TO others” and the total spillover it receives “FROM others”. The NET values are reported along the last row of Table 2. It appears that Manufacturing (MFG), Professional and business services (PRO), and Wholesale Trade (WST) are the largest transmitters, receivers, and have the largest positive net total directional spillovers (26.20%, 24.08% and 16.42%, respectively). At the other extreme, Agriculture, forestry, fishing, and hunting (AGR), Mining (MNG), and Utilities (UTL) stand out as the three major net shock receivers in the network, as shown by their corresponding large negative net total directional spillovers (− 26.51%, − 21.54%, and − 20.11%, respectively). Figure 2 provides a graphical illustration of the total directional spillovers (TO, FROM, NET) by sector.

Fig. 2
figure 2

Directional spillovers in 2005–2020. This figure illustrates the total directional spillover from a sector TO others, the total directional spillover FROM others to a sector, and the net directional spillover of a sector. The measures are obtained from the generalized variance decomposition approach over the period 3 January 2005–31 December 2020

Finally, the total shock spillover among all the sectors is estimated to be 85.8%. Our estimate for the total connectedness among the US sectors is substantial but compares well with the 78% total connectedness index for 13 major US financial institutions in Diebold and Yilmaz (2014), and the 82% index of connectedness for nine equity markets in Diebold and Yilmaz (2015). Yarovaya et al. (2016) report a total shock transmission of 71% between 21 stock markets.Footnote 10 Taken together, the results reported in this section point to significant connectedness among the 14 US economic sectors.

4.2 Dynamic analysis: rolling-sample spillovers

The analysis in Sect. 4.1 employs information over the whole 2005–2020 period to calculate unconditional shock spillover measures among the 14 US sectors. While our previous findings reveal strong connectedness among the US sectors, they do not uncover the dynamics of the spillovers, i.e., how the connectedness evolves over time. During the 16-year period under examination, the US economy went through different market conditions and experienced a number of major shocks. The impacts of both the GFC and COVID-19 on the US stock market prompted the Federal Reserve to cut its policy rate to near zero and embark on large scale asset purchases to prop up the US economy. The COVID-19 pandemic also has profound negative impacts on the world economy as lockdown restrictions caused consumer demand on non-essential items to evaporate. It is possible that the nature of spillovers among the 14 sectors differed during these turbulent periods. We examine this conjecture by estimating the spillover measures over 200-day rolling-sample windows.Footnote 11

Figure 3 plots the time series of the total spillover measure over the rolling windows, alongside the volatility of the US stock market, measured as the standard deviation of the daily returns of the CRSP index over the corresponding 200-day windows. At first glance, the plot suggests that the total spillover between sectors remains considerably large at around 70–90% over the entire 2005–2020 period. It also reveals distinguished patterns across different market conditions. The first cycle spans the pre-crisis period 2005–2006 where the economic climate was relatively calm and inflation was low and stable, with the estimated total spillover of about 80% to 85%. The second episode reflects the large fluctuations during the financial crisis where the total spillover between sectors surges to nearly 90% at the end of 2008—the beginning of 2009 and market volatility peaks. Despite some ups and downs, the total spillover persists at this high level after the crisis and peaks once more in 2012, coinciding with the implementation of the third round of quantitative easing (QE3). The third cycle spans the second half of 2012–2016 when the total spillover reverses to its pre-crisis level of around 80–85%. The total spillover experiences a decline to just below 70% along with the lowest level of market volatility in 2017 and remains rather stably at around 70–80% during the calm market in 2018–2019. The last cycle observes a jump in the total spillover to nearly 90%, together with a high level of market volatility following the COVID-19 outbreak at the beginning of 2020. Both total spillover and market volatility decrease at the end of 2020 along with the vaccination rollout. Overall, we observe substantial total spillovers across the rolling windows. It is clear that the total spillovers among the 14 US sectors exhibit a time-varying nature, where shock spillovers tend to increase in turbulent market conditions.

Fig. 3
figure 3

Total spillover—rolling windows. This figure plots the total spillover between 14 sectors and market volatility over 200-day rolling sample windows. The horizon axis shows the ending dates of the rolling samples

For each of the 200-day rolling windows, we calculate the total directional spillovers from a sector “TO others”, “FROM others” to a sector, and the “NET” directional spillovers, which are computed as the difference between the first two spillover measures. Figure 4 plots the time series of these measures by sector.Footnote 12 For the purpose of presenting our results, we separate the 14 sectors into two groups. Figure 4.1 illustrates the results for the net shock transmitters, i.e., the sectors with (generally) positive NET total directional spillovers. Figure 4.2 refers to the net shock receivers, i.e., the sectors with (generally) negative NET total directional spillovers. While the “FROM” plots appear to be relatively smooth, the total spillovers from a sector “TO others” exhibit noticeably larger variation over time. This seems to suggest that shocks transmitted to a sector from other sectors wash out in aggregate as argued by Lucas (1977). Different patterns emerge when we compare the plots of the shock transmitters to those of the shock receivers. We observe that the “FROM” plots in Fig. 4.1 are smoother than the corresponding plots in Fig. 4.2. This implies that returns for sectors that are net shock receivers (e.g., Agriculture, forestry, fishing and hunting—AGR and Mining—MNG) are likely to be more susceptible to external shocks compared to returns for sectors that are net shock transmitters, such as Manufacturing (MFG).

Fig. 4
figure 4figure 4

Directional spillover—rolling windows. This figure plots the total directional spillovers from a sector to others, the total directional spillover from others to a sector and the net directional spillovers of 14 US sectors over 200-day rolling sample windows

The dynamics of the total directional spillovers also reflect how the COVID-19 pandemic has changed the way the world operates. Retail Trade (RT), which is usually a shock transmitter, experiences the largest decrease in shock transmission to other sectors in 2020, and becomes a shock receiver with the NET total directional spillover of almost -30%. In contrast, we observe an increase in the total directional spillovers from Transportation and warehousing (TPW) and Finance, insurance, real estate, rental and leasing (FIN) “TO others”. The Educational services, health care, and social assistance (EH) sector also switches from a shock receiver to a shock transmitter, with a significant increase in the total directional spillovers “TO others”. More detailed discussion of the pairwise spillovers between sectors during COVID-19 is provided in our sub-sample analysis in Sect. 6.

Part II: Business linkages and shock spillovers

5 The strength of the business linkages and shock spillovers between sectors

All the results reported in Part I point to significant shock spillovers among the 14 US sectors. In Part II, we extend the error variance decomposition analysis by examining if the trade flows between pairs of sectors can explain the degree of return spillovers. Our hypothesis is that the financial connectedness observed in our data is driven by the commercial linkages between sectors, i.e., larger trade flows between sectors are associated with larger shock spillovers.

To formally test our hypothesis, we first construct variables which measure the strength of the trade linkages between each pair of sectors. Such measures based on data from the Input–Output (IO) accounts have been used in a number of finance studies. Menzly and Ozbas (2010), for instance, investigate the return spillovers between supplier and customer industries. They use the trade flows between industries, obtained from the IO tables, as weights to construct two separate portfolios of the representative supplier and the representative customer for each industry. Their results suggest that returns on an industry are affected by the lagged returns of its representative supplier and customer industries. Aobdia et al. (2014) also find significant return interdependence between an industry and its “related” industry, which is constructed as a portfolio of all its trading partner (both supplier and customer) industries. Ahern (2013) finds that return spillovers between industries are affected by the closeness of the relationships. Shocks to an industry return have immediate effects on its closely-related industries and delayed effects on its distant-connected industries. He also demonstrates that a central industry has greater exposure to market risk than a distant one. As a result, its stock returns comove more closely with the market returns. Ahern and Harford (2014) utilize the IO data to create a network of suppliers and customers and gauge the strength of the links between industries. They find that merger activities spill over the supplier-customer network in wave-like patterns.

Details of the construction of these variables and some descriptive statistics are discussed in Sect. 5.1 below. In what follows, we perform a series of cross-sectional regressions of the spillover measures, calculated from the full sample, on the measures of business linkages. Section 5.2 presents the results regarding the pairwise spillovers, while Sect. 5.3 analyzes how the importance and closeness of a sector to the others affect the total directional spillovers.

5.1 Business-linkage measures

We use the information from the IO accounts, provided by the Bureau of Economic Analysis (BEA), to measure the strength of economic linkages between pairs of sectors. These accounts report the value of commodities (goods and services) produced and traded among all the sectors in the US economy. BEA use information from the Economic Census to provide benchmark IO tables every five years (i.e., years ending 2 and 7) and estimated tables for the other years. In this paper, we use the IO accounts aggregated at the sector level, which consist of 15 sectors and 17 commodities. Our analysis includes the 14 non-governmental sectors which have been used in the variance decomposition analysis discussed above.

5.1.1 Measure construction

Our measures of economic linkages between pairs of sectors are calculated based on the information from the two main IO tables: the Make and the Use tables (for snapshots of these tables see Appendix). The Make table reports the value of each commodity produced by sectors. In this table, each row presents a sector while commodities are shown across columns. It is noteworthy that, although a commodity is mainly produced by a sector, the same commodity can also be supplied by other sectors. Equivalently, while a sector predominantly produces one commodity, it may also supply other commodities to other sectors. Thus, the total output of a sector, denoted by \({OUTPUT}_{i}\), is calculated as the sum of all the entries along the corresponding row, while the total output of a commodity produced by all the sectors is computed as the sum of all the entries in the corresponding column.

The Use table reports the value of each commodity purchased as an input for the production of other sectors (or consumed by final users). Elements along each row in this table show the commodity output while elements in each column present the sector input. Therefore, the total commodity output is computed as the sum of all the entries along the corresponding row while the total input value for a sector, denoted by \({INPUT}_{j}\), is calculated as the sum of all the commodity entries in the respective column. The sum of the total sector input and the total value added is the total sector output (\({OUTPUT}_{j}\)), presented in the last row of the Use table.

Following Ahern and Harford (2014) and Becker and Thomas (2011), we use the trading values of the commodities reported in the Make and the Use tables to construct the CUST and the SUPP matrices. These matrices quantify the importance of sectors as customers and suppliers to one another, respectively. Using the information in the Make table, we first calculate the subordinate SHARE matrix which shows the share of each sector’s commodity output in the total supply of each commodity. Specifically, the element in row \(i\), column \(c\) in the SHARE matrix, denoted \({SHARE}_{ic}\), is computed as:

$${SHARE}_{ic}=\frac{{Make}_{ic}}{{Total\,Supply}_{c}}$$
(9)

where \(i\) and \(c\) represent sector and commodity, respectively. The numerator \({Make}_{ic}\) shows the value of commodity \(c\) produced by sector \(i\) (element in row \(i\), column \(c\) of the Make table) while the denominator \({Total Supply}_{c}\) is the total supply of commodity \(c\), which is the total output of commodity \(c\) produced by all sectors in the economy (i.e., the sum of all entries in the commodity \(c\) column in the Make table) plus other components such as imports or changes in inventories.

Using elements from the SHARE matrix along with the information from the Use table, we construct the REVSHARE matrix, which contains the value of commodities transacted between sectors. The element in row \(i\), column \(j\) of the matrix, denoted by \({REVSHARE}_{ij}\), shows the total value of all the commodities which sector \(j\) purchases from sector \(i\) and is given by:

$${REVSHARE}_{ij}=\sum_{c=1}^{C}\left({SHARE}_{ic}\times {Use}_{cj}\right)$$
(10)

where \({SHARE}_{ic}\), defined previously in Eq. (9), is the proportion of commodity \(c\) produced by sector \(i\) (i.e., element in row \(i\), column \(c\) of the SHARE matrix) and \({Use}_{cj}\) (i.e., the element in row \(c\), column \(j\) of the Use table) presents the value of commodity \(c\) purchased by sector \(j\).Footnote 13 The commodity flows between sectors, obtained from the REVSHARE matrix, are the key ingredients in the construction of the CUST and the SUPP matrices, whose elements are computed as follows:

$${CUST}_{ij}=\frac{{REVSHARE}_{ij}}{{OUTPUT}_{i}}$$
(11)
$${SUPP}_{ij}=\frac{{REVSHARE}_{ij}}{{INPUT}_{j}}$$
(12)

As shown in Eqs. (11) and (12), \({CUST}_{ij}\) (the element in row \(i\), column \(j\) in the CUST matrix) shows the fraction of sector \(i\)’s revenue accounted for by sector \(j\) and is obtained by dividing \({REVSHARE}_{ij}\) by the total output value of sector i \({(OUTPUT}_{i}\) in the Make table). \({SUPP}_{ij}\) (the element in row \(i\), column \(j\) in the SUPP matrix) shows the fraction of sector \(j\)’s total inputs purchased from sector \(i\) and is calculated by dividing \({REVSHARE}_{ij}\) by the total value of sector \(j\)’s inputs (\({INPUT}_{j})\).Footnote 14 Thus, when considering a pair of sectors \(i\) and \(j\), we can compute a total of four business linkage measures: \({CUST}_{ji}\), \({SUPP}_{ij}\), \({CUST}_{ij}\), and \({SUPP}_{ji}\); they quantify the importance of the customer role of \(i\) to \(j\), the supplier role of \(i\) to \(j\), the customer role of \(j\) to \(i\), and the supplier role of \(j\) to \(i\), respectively. These variables will be employed as regressors in the cross-sectional regressions, where the regressants are the spillover measures from the variance decomposition matrix. The coefficients on these measures represent the impact of the trading relationship on the shock spillovers between sectors.

5.1.2 Summary statistics

We construct the measures of economic linkages between pairs of sectors using the IO tables for each year and use in our analysis the average values over the sample period.Footnote 15 Table 3 presents some summary statistics. For each sector, the columns labelled Average report the mean SUPP (CUST) values of a sector relative to all its partners over the sample period; they represent the average role of a sector as a supplier to (customer of) the other 13 sectors. We identify whether a sector serves as a main supplier or a main customer and report in the table the number of trading partners for which the SUPP and the CUST variables are larger than the thresholds of 1%, 5% and 10% (see also Horvath, 1998). A graphical representation of the supplier and customer network between sectors, using the average SUPP (CUST) values over the sample period, is provided in Fig. 5. The direction of an arrow shows the role of a sector as a supplier (customer) to its trading partner. The thickness and the shade of the arrow represent the strength of the relationship, based on the 1%, 5%, and 10% thresholds.

Table 3 Summary statistics of CUST and SUPP variables
Fig. 5
figure 5

SUPP and CUST network. This figure shows the supplier and customer network of 14 sectors, using the average values of the CUST and SUPP variables during the whole sample period. The direction of the arrow shows the role of a sector in the trading relationship (i.e., a sector as a supplier / customer of the other sector). The thickness and colour of the arrow show the strength of the business linkage. A thick black arrow represents a relationship variable greater than 10%. A medium dark grey arrow represents a relationship variable between 5–10%. A thin grey arrow represents a relationship variable between 1–5%. The thinnest light grey arrow represents a relationship variable smaller than 1%

Considering the supplier role of a sector, according to Table 3 and Fig. 5, three sectors stand out: Professional and business services (PRO), Manufacturing (MFG), and Finance, insurance, real estate, rental and leasing (FIN) supply, on average, 9.75%, 9.54% and 9.29% of inputs to other sectors, respectively. Their dominant position in the supplier network is also confirmed by the number of sectors they serve as main suppliers: each of these three sectors provides over 1% of the inputs of each of the other 13 sectors, more than 5% of the inputs of 10–12 sectors, and more than 10% of the inputs of 4–7 sectors.

The finding that a few sectors play a disproportionately important role as input suppliers to others supports the notion of the first-order interconnections as described by Acemoglu et al. (2012). In their theoretical model, shocks to a dominant sector are shown to generate strong spillovers to direct customer sectors, leading to large aggregate volatility. Our earlier variance decomposition results are in line with Acemoglu et al. (2012). According to the estimated measures of directional spillovers, Professional and business services (PRO), Manufacturing (MFG) and Finance, insurance, real estate, rental and leasing (FIN) are found to be net shock transmitters. Importantly, the estimated net spillovers for these three sectors (PRO, MFG and FIN) are relatively large compared to the values for the other net shock transmitters.

On the other side of the coin, Educational services, health care, and social assistance (EH), Agriculture, forestry, fishing, and hunting (AGR), and Other services, except government (OS) are the three sectors that account for the smallest proportion of other sectors’ inputs, supplying only 0.28%, 0.45% and 0.81% of inputs to other sectors, respectively. While the latter two (AGR and OS) are among the sectors with the smallest outputs in the economy, Educational services, health care, and social assistance (EH) is the second largest supplier to final users, and has the largest proportion of its revenue accounted for by personal consumption. Consistent with the findings reported in Table 3 and Fig. 5, the results for the variance decomposition matrix shown in Table 2 show that these three industries are net shock receivers.

In what regards the customer role, the descriptive statistics in Table 3 designate Manufacturing (MFG) as the largest customer sector, accounting on average for 12.89% of the other sectors’ revenue. This sector contributes to more than 1% of the sales of 11 (out of 13) sectors, and to more than 10% of the revenues of 5 other sectors, as also illustrated by the number of thick arrows in Fig. 5. Once again, the importance of MFG in the supplier-customer network is consistent with its output magnitude and key role in a modern economy. Finance, insurance, real estate, rental and leasing (FIN) ranks as the second largest customer, with the average CUST value of 4.43%. Mining (MNG) and Agriculture, forestry, fishing, and hunting (AGR) are the smallest customers. A relative comparison of the average values of the CUST and SUPP variables together with the numbers of main customers and suppliers show that US sectors tend to have slightly more diversified customers than suppliers.

Table 3 also reports the number of trading partners with which a sector has strong linkages—either in terms of CUST or SUPP—at the 1%, 5% and 10% thresholds, respectively. These numbers show the closeness of a sector to other sectors in the network regardless of its role. Manufacturing (MFG), Wholesale trade (WST), Professional and business services (PRO) and Finance, insurance, real estate, rental and leasing (FIN) are the central sectors in the economy. Each of these sectors has strong trade linkages (higher than 1%) with all their trading partners (i.e., 13 other sectors). Agriculture, forestry, fishing, and hunting (AGR) is the least-connected sector; it is closely related to two other sectors at the 1% threshold, and has no close linkage at the 5% and 10% thresholds.

The last two columns of Table 3 report information about sector size: the number of firms and the number of establishments. This data is obtained from the US Census Bureau’s Statistics of US Businesses (SUSB).Footnote 16 Utilities (UTL) is the smallest sector, with 6,003 firms and 17,759 establishments, whereas Professional and business services (PRO) is the largest sector, with 1,141,406 firms and 1,315,263 establishments.

5.2 Pairwise spillovers and business linkages

For each pair of sectors \(i\) and \(j\), we obtain two directional spillover measures showing the fraction of the forecast error variance of sector \(i\) which is due to shocks to sector \(j\), denoted by \({S}_{i\leftarrow j}^{H}\), and the fraction of the forecast error variance of sector \(j\) which is due to shocks to sector \(i\), denoted by \({S}_{j\leftarrow i}^{H}\). These measures are the \({ij}^{th}\) and the \({ji}^{th}\) elements in the variance decomposition matrix shown in Table 2 and are employed as the dependent variable in the cross-sectional regression model specified in Eq. (13) below.

$${S}_{i\leftarrow j}^{H}={\alpha }_{0}+{\alpha }_{1}{CUST}_{ij}+{\alpha }_{2}{SUPP}_{ji}+{\alpha }_{3}{CUST}_{ji}+{\alpha }_{4}{SUPP}_{ij}+{\alpha }_{5}{Size}_{j}+{\alpha }_{6}{Size}_{i}+{u}_{ij}$$
(13)

where \({S}_{i\leftarrow j}^{H}\) is the pairwise spillover from sector \(j\) to sector \(i\), showing the proportion of the \(H\)-step-ahead forecast error variances of sector \(i\) accounted for by shocks to sector \(j\). The first four regressors—\({{CUST}_{ij}, { SUPP}_{ji}, CUST}_{ji},\) and \({SUPP}_{ij}\)—are the business linkage variables constructed from the IO accounts. They quantify the importance of a sector as a supplier or a customer for its trading partner. \({SUPP}_{ij}\) shows the supplier role of sector \(i\) to sector \(j\), while \({CUST}_{ij}\) shows the customer role of sector j to sector \(i\).\({SUPP}_{ji}\) and \({CUST}_{ji}\) are defined similarly. In line with Ahern (2013), we control also for the impact of the sector size on shock spillovers. We expect a larger sector, with a large number of constituent firms, to have stronger influence on other sectors in the economy. We construct 91 trading pairs from the 14 sectors in our sample, giving us 182 pairwise directional spillover observations. Note that because the value for the dependent variable is obtained from the variance decomposition matrix, we follow Lewis and Linzer (2005) and Weiß et al. (2014) to estimate and report the heteroskedasticity-consistent standard errors as discussed in White (1980).

Table 4 reports the estimates of Eq. (13). We find that the degree of pairwise spillovers from sector \(j\) to sector \(i\) is positively related to the supplier share of sector \(j\) in the total input of sector \(i\). In other words, shocks to an important supplier sector contribute strongly to its partner’s forecast error variance. This finding highlights the importance of the supplier role for intersectoral shock spillovers. It may be attributed (partly) to the high aggregation at the sectoral level, which reduces considerably the possibility of substituting input suppliers. Our results are in line with the findings in Barrot and Sauvagnat (2016) who document the negative impact on a firm’s sales growth following input disruption (caused by a natural disaster) to the firm’s supplier.

Table 4 Pairwise spillovers and business linkages in 2005–2020

The size of the trading sectors appears to positively correlate with the pairwise shock spillovers as the coefficients associated with the number of firms (column 1) and the number of establishments (column 2) are positive and statistically significant. This suggests that a sector which has a large number of firms or establishments, with an extended trading network with other sectors, tends to easily transmit and receive shocks to/ from other sectors. This finding is in line with Herskovic et al. (2020), which confirms the size impacts on volatility transmission at firm level.

5.3 Total directional spillovers and business linkages

We now investigate whether total directional spillovers are affected by the role a sector plays as either a supplier or a customer. Table 5 reports the results for the total directional spillover from a sector TO others (columns 1–3), FROM others to a sector (columns 4–6), and the NET total directional spillover (columns 7–9). The explanatory variable is the number of trading partners with close business linkages between a sector and the others regardless of its role (i.e., whether supplier or customer to the other sectors), as reported in the last three columns of Table 3 at various thresholds. The results reported in Table 5 reveal a positive significant relationship between the total directional spillovers and the number of close economic linkages at all three thresholds: 1%, 5% and 10%. Our findings are in line with Carvalho (2014), who demonstrates that hub sectors, which are close to the majority of other sectors, comove more with others. Ahern (2013) also indicates that shocks to sectors which are central in the network are more likely to contribute more to the aggregate system-wide volatility. Given the small sample size, however, these findings should be interpreted with care.

Table 5 Total directional spillovers and business linkages in 2005–2020

Part III: Extensions and robustness tests

6 Market conditions

Our 16-year sample period spans different market conditions. In this section, we examine the impact of economic relationship on shock transmissions during the sub-periods: the pre-crisis 2005–2006, the financial crisis 2007–2008, the post-crisis 2009–2019, and the COVID-19 pandemic 2020. Specifically, we fit the VARX model for 14 sectors in each sub-period and calculate the spillover measures using the generalized forecast variance approach. We then regress these measures on the average business linkage variables calculated for each sub-period. Panels A-D of Table 6 show the error variance decomposition matrices computed using data for the sub-periods. We immediately notice that the total spillovers increase from 82.32% in the pre-crisis period to 88% in the crisis and 84.31% in the post-crisis period. The lowest total directional spillover from a sector TO others is 38.37% (Mining—MNG) in Panel A, which is considerably smaller than 61.74% (Mining—MNG) in Panel B, and 47.8% (Agriculture, forestry, fishing, and hunting—AGR) in Panel C. The lowest total directional spillover FROM others to a sector also increases from 69.83% (Mining—MNG) in the pre-crisis period to 83.62% (Mining—MNG) and 76.81% (Agriculture, forestry, fishing, and hunting—AGR) in the crisis and post-crisis period, respectively. In terms of pairwise spillover, the smallest cross-variance share is 1.23% from Mining (MNG) to Retail trade (RT) in the sub-periods 2005–2006, 2.96% from Mining (MNG) to Finance, insurance, real estate, rental and leasing (FIN) in 2007–2008, and 3.42% from Agriculture, forestry, fishing, and hunting (AGR) to Mining (MNG) and to Utility (UTL) in the sub-period 2009–2019. Most net spillovers declined relative to their values prior to the financial crisis. This is true for both the largest net shock receivers and transmitters. Overall, our findings confirm stronger shock spillovers between sectors during and after the financial crisis.

Table 6 Shock spillover—separate samples

Panel D reveals the changes in sectoral shock spillovers due to the COVID-19 outbreak in 2020, which is consistent with our findings in the rolling window analysis. As compared to the spillovers in Panel C, the greatest change in the total directional spillovers, hence the greatest changes in the net spillovers, is a decrease of 28.41% in Retail Trade (RT). A closer look at the pairwise spillover shows that the largest changes in pairwise spillovers are also related to the Retail Trade (RT) sector: an increase of 4.79% in the shock spillover that Retail Trade (RT) receives from Information (INF) and a decrease of 4.78% in the shock spillover from Retail Trade (RT) to Arts, entertainment, recreation, accommodation, and food services (ES). The fact that Retail trade (RT) changes from a shock transmitter to a shock receiver in the system reflects how COVID-19 affects business operations. Retail businesses have been hardly hit by lockdown and social distancing measures implemented in many countries, including the US. Thus, it is understandable why Retail trade (RT) becomes less influential and more prone to shock spillovers. In the same vein, Arts, entertainment, recreation, accommodation, and food services (ES) also suffers a decline in shock transmission, with a decrease of 1.29% of pairwise spillover to others on average. Finally, due to the pandemic, pairwise shock spillovers from Educational services, health care, and social assistance (EH) to every other sector increase by an average of 2.06%, resulting in an increase of 26.72% in the total directional spillover TO others, changing it from a shock receiver to a shock transmitter.

Tables 7 and 8 report the results of the extended analysis in different market conditions.Footnote 17 All coefficients representing the impacts of the supplier/customer roles on intersectoral shock transmissions preserve their signs and significance across the three sub-periods. The magnitude of the estimates on the trade linkage variables weakens slightly during the crisis period, but recovers in the post-crisis period. Similar to Ahern (2013), we note that the economic linkages are not the only channel for shock transmissions between sectors. In addition to their connectedness in the production network, sectors connect to each other via complicated links of their constituent firms, in terms of various relationships such as ownership and geographic closeness. Thus, in periods of stress, shock spillovers may be more strongly affected by other relationships than their economic linkages. However, even though the impacts of the business linkages weaken in the crisis period, our results clearly demonstrate that they invariably matter for intersectoral shock transmissions in all market conditions.

Table 7 Pairwise spillovers and business linkages—separate samples
Table 8 Total directional spillovers and business linkages—separate samples

7 Robustness checks

This section collects a series of sensitivity tests to check for the robustness of our results.

Non-financial sectors. Regardless of its actual trading relationship with the other sectors, the financial sector is more sensitive to shocks from others and more likely to transmit significant shocks to other sectors. To ensure that our results are not driven by it, we drop the Finance, insurance, real estate, rental, and leasing (FIN) sector from our sample. We confirm the importance of the business linkages on volatility spillovers between non-financial sectors (results available upon request).

Forecast horizon and VARX order. We check the sensitivity of our results to the choice of forecast horizon and VARX order. We estimate the model in Eq. (13) using the spillover measures obtained from the VARX specification of orders from one to four as well as the forecast horizons of 1 to 5 days and 20 days. We obtain qualitatively similar results (available upon request) to those reported in Tables 4 and 5.

Alternative business-linkage variables. Finally, we check whether our cross-sectional regression results are robust to the choice of the IO tables. In addition to using the average values of the business linkage variables over the full sample period, we use the variables calculated from the IO tables of every year from 2005 to 2019. The unreported results are qualitatively similar to our main findings, confirming that the impact of business linkages on shock spillovers between sectors persists regardless of the choice of year for the IO tables.

8 Conclusions

We investigate the degree of financial connectedness as well as the characteristics of shock transmissions among the 14 US economic sectors using a VARX system based on the error variance decomposition approach proposed by Diebold and Yilmaz (2012, 2014). The total shock spillover among all sectors is found to be considerably large (85.8%), pointing to significant financial connectedness among the sectors. We find that, on average, shocks to a sector’s return can explain 6.6% of the forecast error variance of its trading partner. Manufacturing (MFG), Professional and business services (PRO), and Wholesale Trade (WST) are the top three net shock transmitters, while Agriculture, forestry, fishing, and hunting (AGR) and Mining (MNG) are the smallest shock transmitters. Our dynamic analysis using 200-day rolling windows and the analysis on subsamples show evidence of stronger intersectoral shocks spillovers in turbulent markets during the 2007–2008 Financial Crisis and the COVID-19 pandemic. We observe an increase in the total directional shock spillovers from Educational services, health care, and social assistance (EH) to other sectors, while Retail trade (RT) becomes more sensitive to shock spillovers, changing from a shock transmitter to a shock receiver. This evidence reflect the impact of COVID-19 on different sectors of the economy.

We also examine the relevance of supplier-customer linkages on the stock return spillovers across sectors. We uncover a strong relation between supplier-customer linkages and shock transmissions between US sectors: shock spillovers from a sector to its partner are significant and can be explained by the supplier role of the sector. In addition, we find that the total directional spillovers from / to a sector are affected by the closeness of the sector to others and that the number of close trading partners of a sector with other sectors is positively related to the net shock spillover which the sector transmits to others.

The robustness of our results is confirmed under different settings. We obtain qualitatively similar results when the analysis excludes the financial sector from the sample. Our results are robust to the use of business linkage variables calculated using information from the IO tables for different years. They are also insensitive to the order of the VARX model and the choice of forecast horizon.

Our findings are useful for investors whose portfolios concentrate in some specific sectors or industries. While sector funds are more manageable and shown to deliver outperformance, they entail greater total and systematic risk. Thus, our evidence of sectoral shock transmission taking into account business linkages provides practical implications for risk management. Investors can observe the shocks to the sector’s major suppliers to better predict the volatility of their positions. Moreover, our findings are useful for policy makers and regulators, since understanding how shocks to a sector can affect others, ultimately resulting in aggregate fluctuation, is crucial, especially in periods of stress such as the financial crisis or the COVID-19 pandemic.

Our work is not without limitations. Our first stage results would likely be affected if investors take various business horizons into consideration. One fruitful direction for future research thus relates to the rich time–frequency dynamics of volatility connectedness in the spirit of Barunik and Krehlik (2018). We have not pursued this avenue as the approach does not match the frequency of the supplier-customer measures employed in our second stage analysis. We therefore leave this unexplored area to future researchers.