Skip to: [Printable Version (PDF)] [Bibliography] [Footnotes]
Finance and Economics Discussion Series: 2015-094 Screen Reader version

Input Linkages and the Transmission of Shocks:
Firm-Level Evidence from the 2011 Tōhoku Earthquake*

Christoph E. Boehm
University of Michigan
Aaron Flaaen
Federal Reserve Board of Governors
Nitya Pandalai-Nayar
University of Michigan

October 22 ,2015

Abstract:

Using novel firm-level microdata and leveraging a natural experiment, this paper provides causal evidence for the role of trade and multinational firms in the cross-country transmission of shocks. Foreign multinational affiliates in the U.S. exhibit substantial intermediate input linkages with their source country. The scope for these linkages to generate cross-country spillovers in the domestic market depends on the elasticity of substitution with respect to other inputs. Using the 2011 Tōhoku earthquake as an exogenous shock, we estimate this elasticity for those firms most reliant on Japanese imported inputs: the U.S. affiliates of Japanese multinationals. These firms suffered large drops in U.S. output in the months following the shock, roughly one-for-one with the drop in imports and consistent with a Leontief relationship between imported and domestic inputs. Structural estimates of the production function for all firms with input linkages to Japan yield disaggregated production elasticities that are similarly low. Our results suggest that global supply chains are sufficiently rigid to play an important role in the cross-country transmission of shocks.


Keywords:Multinational Firms, International Business Cycles, Business Fluctuations, Elasticity of Substitution

JEL Classification: F44, F23, E32, F14

The spillover effects of trade and financial linkages has been a preeminent topic in international economics in recent decades. The large expansions in trade and foreign direct investment (FDI) in the past twenty years have generated much discussion on whether they increase volatility di Giovanni & Levchenko (2012), increase comovement Burstein et al. (2008); Frankel & Rose (1998) or lead to less diversified production and specialization Imbs (2004). Identifying the micro-foundations underlying the role of these linkages in the increased interdependence of national economies is challenging. Advanced economies are highly connected, and most variables influenced by any candidate mechanism are often correlated with other developments in the source and destination countries. There is often little in the way of exogenous variation to isolate any particular mechanism from a host of confounding factors. Moreover, the requisite data to examine these issues at the necessary detail and disaggregation have been, until recently, unavailable.

This paper provides empirical evidence for the cross-country transmission of shocks via the rigid production linkages of multinational firms. The principal mechanism at work is not new; the idea of input-output linkages as a key channel through which shocks propagate through the economy dates back to at least Leontief (1936) or Hirschmann (1958). Two advances in this paper permit a new quantitative evaluation of the nature and magnitude of these linkages. First, we utilize a novel dataset that, for the first time, links restricted U.S. Census Bureau microdata to firms' international ownership structure. This information permits a forensic focus on particular firms and their underlying behavior. Second, we utilize the March 2011 Tōhoku earthquake and tsunami as a natural experiment of a large and exogenous shock disrupting the production linkages originating from Japan.

We study the role of imported intermediate inputs in the transmission of this shock to the United States economy. Because disruptions to imports of final goods would be unlikely to affect U.S. production, we develop a new methodology for isolating firm-level imports of intermediate inputs. We show that the U.S. affiliates of Japanese multinationals are the most natural source of this transmission, due to their high exposure to imported intermediates from Japan. The scope for shocks to these imported inputs to pass through and affect the firm's U.S. production depends on how substitutable they are with inputs from alternative sources. In other words, the role of imported inputs in the transmission of shocks is governed by the elasticity of substitution with respect to domestic factors of production.

We estimate this elasticity using the relative magnitudes of high frequency input and output shipments in the months following the Tōhoku earthquake/tsunami. This proceeds in two steps. First, reduced form estimates corresponding to Japanese multinational affiliates on average show that output falls, without a lag, by a comparable magnitude to the drop in imports. These results suggest a near-zero elasticity of imported inputs. Second, we structurally estimate a firm-level production function that allows for substitution across different types of inputs. The structural estimation procedure we use is uniquely tailored to the experiment. In an initial period prior to the Tōhoku disruption, we infer information on the firm's productivity and optimal input mix. Then, applying this production function to the period of the disruption, we estimate the elasticity parameters based on how changes in the firm's input mix translate into changes in output.

This estimation strategy has a number of attractive features. Most importantly, it relies on very few assumptions. Direct estimation of the production function circumvents the many difficulties associated with specifying a firm's optimization problem in the period after the shock. Second, it yields transparent parameter identification. This is an advantage over traditional estimation strategies as it does not suffer from omitted variables and endogeneity concerns arising from correlated shocks. Third, it allows for the estimation across different subgroups of firms.

The structural estimates are broadly in agreement with the results from our reduced form exercise. For Japanese multinationals, the elasticity of substitution across material inputs is 0.2 and the elasticity between material inputs and a capital/labor aggregate is 0.03. For non-Japanese firms using inputs from Japan, the estimates of the elasticity of substitution across material inputs are somewhat higher, ranging from 0.42 to 0.62. While the high cost share and particularly low elasticity for Japanese affiliates explains their predominant contribution to the direct transmission of this shock to the U.S., the elasticity estimates for non-Japanese firms are still substantially lower than typical estimates used in the literature. We argue that the substantial share of intra-firm intermediate trade implies greater complementarities in aggregate trade than is currently recognized.

There are a number of important implications for such low values of the elasticity of substitution. This parameter appears in various forms in a wide span of models involving the exchange of goods across countries. As discussed by Backus et al. (1994) and Heathcote & Perri (2002) among others, this parameter is critically important for the behavior of these models and their ability to match key patterns of the data. Prior estimates of this parameter were based on highly aggregated data that suffered from concerns about endogeneity and issues of product composition.1 Reflecting the uncertainty of available estimates for this elasticity, it is a common practice to evaluate the behavior of these models along a wide range of parameter values.

It is well known that a low value for this parameter (interpreted as either substitution between imported and domestic goods in final consumption or as intermediates in production) improves the fit of standard IRBC models along several important dimensions. In particular, the elasticity of substitution plays a role in two highly robust failings of these models: i) a terms of trade that is not nearly as variable as the data, and ii) a consumption comovement that is significantly higher than that of output, whereas the data show the opposite relative ranking.2

To understand the relationship between the elasticity and comovement, it is helpful to recall that these models generate output comovement by inducing synchronization in factor supplies, a mechanism that by itself generally fails to produce the degree of comovement seen in the data. Complementarities among inputs together with heterogeneous input shocks will generate direct comovement in production, augmenting the output synchronization based on factor movements. Burstein et al. (2008) show that a low production elasticity of substitution between imported and domestic inputs reduces substitution following relative price movements, and thereby increases comovement. Johnson (2014) is a more recent and detailed treatment of this topic. It is also straightforward to see how a lower elasticity increases volatility in the terms of trade. When two inputs are highly complementary, deviations from the steady state mix are associated with large changes in their relative prices. In the words of (Heathcote & Perri, 2002, page 621): "greater complementarity is associated with a larger return to relative scarcity." The estimates in this paper have implications for the role of trade in firm-level and aggregate volatility. Other research has argued that firms can diversify risk arising from country specific shocks by importing (Caselli et al. (2014)) or that firms with complex production processes of several inputs are less volatile as each input matters less for production (Koren & Tenreyro (2013)). On the other hand, there is a well-established fact that complementarities and multi-stage processing can lead to the amplification of shocks as in Jones (2011) and Kremer (1993). We discuss the potential for measured amplification in our context in Section 4.

This paper also contributes to the empirical evidence on the role of individual firms in aggregate fluctuations, emanating from the work of Gabaix (2011). Other related evidence comes from di Giovanni et al. (2014), who use French micro-data to demonstrate that firm-level shocks contribute as much to aggregate volatility as sectoral and macroeconomic shocks combined. The so-called granularity of the economy is evident in our exercise; though the number of Japanese multinationals is small, they comprise a large share of total imports from Japan, and are arguably responsible for a measurable drop in U.S. industrial production following the Tōhoku earthquake (see Figure 3).

The strong complementarity across material inputs implies that non-Japanese input use falls nearly proportionately, thereby propagating the shock to other upstream (and downstream) firms in the U.S. economy and abroad. Many suppliers were thus indirectly exposed to the shock via linkages with Japanese affiliates that had i) high exposure to Japanese inputs and ii) a rigid production function with respect to other inputs. Network effects such as these can dramatically magnify the overall transmission of the shock (both across countries and within). And while such effects are commonly understood to exist, we provide unique empirical evidence of the mechanisms at work.

As is the case with most research based on an event-study, care should be taken in generalizing the results to other settings. As argued by Ruhl (2008), the elasticity of substitution is necessarily tied to the time horizon and nature of shocks to which it is applied. More generally, one might worry that the composition of Japanese trade or firms engaged in such trade is not representative of trade linkages more broadly. We believe our results are informative beyond the context of this particular episode for two reasons. First, the features of Japanese multinationals that are underlying the transmission of this shock are common to all foreign multinational affiliates in the U.S.3 Second, estimates corresponding to all firms in our sample also exhibit substantial complementarities, and as a whole these firms account for over 70 percent of U.S. manufacturing imports.

The next section describes the empirical strategy and data sources used in this paper. Section 2 presents reduced form evidence in support of a low production elasticity of imported inputs for Japanese multinational affiliates. In Section 3, we expand the scope of parameters we identify with a structural model of cross-country production linkages, and estimate the parameters of this model across several firm subgroups. Section 4 discusses the implications of these estimates, and details issues of aggregation and external validity. The final section concludes.


1 Empirical Strategy and Specification

This section outlines the empirical approach of using an event-study framework surrounding the 2011 Tōhoku event to estimate the production elasticity of imported inputs. We discuss the relevant details of this shock, document the aggregate effects, and then outline the empirical specification for the firm-level analysis.


1.1 Background

The Tōhoku earthquake and tsunami took place off the coast of Northeast Japan on March 11, 2011. It had a devastating impact on Japan, with estimates of almost twenty thousand dead or missing (Schnell & Weinstein (2012)) and substantial destruction of physical capital. The magnitude of the earthquake was recorded at 9.0 on the moment magnitude scale ( $$ \mathrm{M}_{\mathrm{w}}$$ ), making it the fourth largest earthquake event recorded in the modern era.4 Most of the damage and casualties were a result of the subsequent tsunami that inundated entire towns and coastal fishing villages. The effects of the tsunami were especially devastating in the Iwate, Miyagi, and Fukushima prefectures. The Japanese Meteorological Agency published estimates of wave heights as high as 7-9m (23-29ft), while the Port and Airport Research Institute (PARI) cite estimates of the maximum landfall height of between 7.9m and 13.3m (26-44ft).

Figure 1 shows the impact of the Tōhoku event on the Japanese economy. Japanese manufacturing production fell by roughly 15 percentage points between February and March 2011, and did not return to trend levels until July. Much of the decline in economic activity resulted from power outages that persisted for months following damage to several power plants - most notably the Fukushima nuclear reactor.5 Further, at least six Japanese ports (among them the Hachinohe, Sendai, Ishinomaki and Onahama) sustained damage and were out of operation for more than a month, delaying shipments to both foreign and domestic locations. It should be noted, however, that the largest Japanese ports (Yokohama, Tokyo, Kobe) which account for the majority of Japanese trade, re-opened only days after the event.

As expected, the economic impact of the event was reflected in international trade statistics, including exports to the United States. Figure 2 plots U.S. imports from Japan around the period of the Tōhoku event, with imports from the rest of the world for comparison. The large fall in imports occurs during the month of April 2011, reflecting the several weeks of transit time for container vessels to cross the Pacific Ocean. The magnitude of this drop in imports is roughly similar to that of Japanese manufacturing production: a 20 percentage point drop from March to April, with a recovery by July 2011.

More striking is the response of U.S. industrial production in the months following the event. Figure 3 demonstrates that there is indeed a drop in U.S. manufacturing production in the months following the Japanese earthquake. Although the magnitudes are much smaller -- roughly a one percentage point drop in total manufacturing and almost two percentage points in durable goods -- the existence of a measurable effect is clear.6

Though tragic, the Tōhoku event provides a rare glimpse into the cross-country spillovers following an exogenous supply shock. This natural experiment features many characteristics that are advantageous for this type of study. It was large and hence measurable, unexpected, and directly affected only one country. On the other hand, the short duration of the shock presents a challenge for measurement as it limits the available datasets with information at the required frequency.

1.2 Data

Several restricted-use Census Bureau datasets form the core of our firm-level analysis. The Longitudinal Business Database (LBD) collects the employment, payroll, and industry of all establishments operating in the United States, and is maintained and updated as described by Jarmin & Miranda (2002). Longitudinal linkages allow the researcher to follow the establishment over time, and the annual Company Organization Survey (COS) provides a mapping from establishments to firms. All of the analysis in this paper will be at the firm-level.

The Longitudinal Foreign Trade Transactions Database (LFTTD) links individual trade transactions to firms operating in the United States. Assembled by a collaboration between the U.S. Census Bureau and the U.S. Customs Bureau, the LFTTD contains information on the destination (or source) country, quantity and value shipped, the transport mode, and other details from point-of-trade administrative documents. Importantly for this study, the LFTTD includes import and export trade transactions at a daily frequency, which is easily aggregated to monthly-level trade flows. A number of important papers have utilized this resource, such as Bernard et al. (2007) and Bernard et al. (2006).

We utilize two novel extensions to this set of Census data products. First, a new link between two international corporate directories and the Business Register (BR) of the Census Bureau provides information on the international affiliates of firms operating in the United States. These directories allow us, for the first time, to identify those U.S. affiliates part of a foreign parent company, as well as those U.S. firms with affiliate operations abroad. This information is an important resource for identifying the characteristics of U.S. firms affected by the Tōhoku event. For information on these directories and the linking procedure, see Flaaen (2014) and Appendix B.1.

The second novel data resource is a system to classify firm-level import transactions as intermediate or final goods. Although intermediate input trade represents as much as two-thirds of total trade (see Johnson & Noguera (2012)), the LFTTD does not classify a trade transaction based on its intended use. To overcome this limitation, we use information on the products produced by U.S. establishments in a given industry to define a set of products intended for final sale for that industry.7 The remaining products are presumably used by establishments in that industry either as intermediate inputs or as capital investment. Details on this classification procedure are available in Appendix B.2. In the aggregate, this firm-level classification procedure yields estimates of the intermediate share of trade that are consistent with prior estimates: 64 percent of manufacturing imports are classified as "intermediates" in 2007.

In Appendix B.3.2 we outline methods to integrate geographic information on the severity of the earthquake to the Japanese locations of U.S.-based firms. Due to limitations in within-Japan supply-chain linkages, however, we condition our sample on this information for robustness purposes only.

The ideal dataset to evaluate the transmission of the Tōhoku event on U.S. firms would consist of high frequency information on production, material inputs, and trade, separated out by geographic and ownership criteria. Information from the LFTTD on import shipments is ideal for these purposes. Census data on production and domestic material usage, however, is limited. The Annual Survey of Manufacturers (ASM) contains such information, but at an annual frequency and only for a subset of manufacturing firms. Recognizing the challenges of high-frequency information on firms' U.S. production, we utilize a proxy based on the LFTTD -- namely the firm's exports of goods to North America (Canada and Mexico). The underlying assumption of this proxy is that all firms export a fixed fraction of their U.S. output to neighboring countries in each period. The advantage of this approach is the ability to capture the flow of goods at a specific point in time. There are few barriers to North American trade, and transport time is relatively short. Moreover, as documented in Flaaen (2014), exporting is a common feature of these firms, of which exports to North America is by far the largest component. The disadvantage of this approach is that it conditions on a positive trading relationship between firms in the U.S. and Canada/Mexico. We will assess the quality of this measure as a proxy for output in section 4.3.1.8

1.3 Basic Theory

Before moving to our firm-level analysis, it is useful to describe the basic theory underlying the features of firm-level production that we estimate. The transmission of shocks within a firm's production chain is governed by the flexibility of production with respect to input sourcing. Rather than model these complex networks directly, the literature typically summarizes this feature with the well-known elasticity of substitution within a C.E.S. production function. Our identification of this elasticity will rely on the relative impacts on output and imported inputs following the shock. Consider the C.E.S. production function

$$\displaystyle x=\left[ \left( 1-\mu \right) ^{\frac{1}{\psi }}\left[ F_{D}\right] ^{\frac{\psi -1}{\psi }}+\mu ^{\frac{1}{\psi }}\left[ IM\right] ^{\frac{\psi -1}{\psi}}\right] ^{\frac{\psi }{\psi -1}}$$ (1)

where output consists of combining a domestic bundle of factors FD (e.g. capital and labor) with a foreign imported input IM. The parameter $$ \mu $$ reflects the relative weight on the input IM in production, conditional on prices and a given elasticity value. Suppose the firm purchases its inputs in competitive markets with prices pD and pM, respectively, and sells its good at price px. Our approach in section 2 will be to estimate the parameter $$ \psi $$ governing the degree of substitution between these inputs, using information on the output elasticity with respect to imported inputs, $$ \frac{\partial \ln p_{x}x}{\partial \ln p_{M}M}$$ , in the months following the shock.

The first order conditions imply that

$$\displaystyle \frac{F_{D}^{\ast }}{IM^{\ast }}=\frac{1-\mu }{\mu }\left( \frac{p_{M}}{p_{D}}\right) ^{\psi },$$ (2)

where $$ F_{D}^{\ast }$$ and $$ IM^{\ast }$$ denote the optimal quantities of inputs. We will show the theoretical foundations underlying the intuitive result that a one-for-one drop in output with the fall in imported inputs implies an elasticity of zero. To do this, we make the following assumptions, all of which we will relax to some degree in the estimation framework in Section 3:
  1. Imported inputs shipments are disrupted, such that the firm receives a suboptimally low quantity of IM: $$ IM<IM^{\ast }$$ ;
  2. The firm is unable to adjust domestic inputs $$ F_{D}^{\ast }$$ or its price px after learning that it receives IM;
  3. The firm does not shut down.
Given these assumptions, the following result holds:
Result 1   Under assumptions 1) to 3):
$$\displaystyle \frac{\partial \ln p_{x}x}{\partial \ln p_{M}IM}=\frac{1}{1+\left( \frac{IM^{\ast }}{IM}\right)^{\frac{\psi-1}{\psi }}\left( \frac{1-\mu }{\mu }\right) \left( \frac{p_{M}}{p_{D}}\right) ^{\psi -1}}\in \left( 0,1\right)$$ (3)

for any $$ \psi \in \left( 0,\infty \right) $$ .
Proof. See Appendix A.1 for details.

An immediate implication of this result is that the output elasticity is unity only when $$ \psi $$ approaches zero.9 In this case $$ \left( \frac{IM^{\ast }}{IM}\right) ^{\frac{\psi-1}{\psi }}\rightarrow 0$$ (recall that $$ IM<IM^{\ast }$$ ) and hence $$ \lim_{\psi \rightarrow 0}\frac{\partial \ln p_{x}x}{\partial \ln p_{M}IM}=1$$ . Hence, observing a one-for-one drop in the value of output with the value of imported inputs, we infer that $$ \psi $$ is close to zero. It is also easy to show that conditional on a value for $$ \psi \in \left( 0,\infty \right) $$ , the output elasticity in (3) is increasing in the parameter $$ \mu $$ . That is, conditional on a given drop in the imported input, a larger weight on this input leads to a larger percent response in output.

Our use of the natural experiment is critical for observing the effects of suboptimal input combinations $$ \left( F_{D}^{\ast },IM\right) $$ . To see this, suppose the firm could freely adjust FD after learning it will receive $$ IM<IM^{\ast }$$ . Then, it would choose FD such that $$ \frac{F_D}{IM}=\frac{F_D^{\ast}}{IM^{\ast}}$$ and the firm would contract one-for-one with the drop in imports. It is a well-known fact that constant returns to scale production functions in competitive environments lead to indeterminate firm size. This has the implication that:

$$\displaystyle \frac{\partial \ln \left( p_{x}x\right) }{\partial \ln \left( p_{M}IM\right) }=\frac{\partial \ln \left( p_{x}x\right) }{\partial \ln \left(p_{D}F_{D}\right) }=\frac{\partial \ln \left( p_{D}F_{D}\right) }{\partial \ln \left( p_{M}IM\right) }=1.$$ (4)

In this case it is not possible to learn anything about $$ \psi $$ from the joint behavior of output and the value of intermediate inputs. We provide evidence below that firms did not significantly adjust their domestic labor force following the disruption, so that a constant FD is indeed a reasonable assumption in this simple framework. To be sure, there are a number of alternative frameworks where such behavior would not hold. We discuss some of these in Appendix A.1, and show that the mapping $$ \lim_{\psi \rightarrow 0}\frac{\partial \ln p_{x}x}{\partial \ln p_{M}IM}=1$$ is more general.


2 Reduced Form Evidence


2.1 Framework

Our analysis of the production function (1) above demonstrates that a natural measure to evaluate the potential conduits of the Tōhoku shock to the United States would be the degree of reliance on Japanese imported inputs. This is best expressed as the firm's cost share of inputs from Japan, and can be constructed in a Census year by taking Japanese imported inputs and dividing by all other inputs (which includes production worker wages and salaries, the cost of materials, and the cost of new machinery expenditures). Exposure to Japanese imported inputs is heavily concentrated among Japanese affiliates. In the year 2007, which is the closest available Census year, this cost share was nearly 22% on average for Japanese affiliates (see Panel A of Table 1), compared to just 1% for other firms. For more detail on the heterogeneity across and within these firm groups, we construct a density estimate of such an exposure measure for the Japanese affiliates and non-Japanese multinationals. The results, shown in Figure 4, show little overlap between these distributions: there are few Japanese affiliates with low exposure to Japanese inputs, and few non-Japanese firms with substantial exposure.10

We now estimate the relative impacts on imported inputs and output for the Japanese affiliates as a group. To do this, we implement a dynamic treatment effects specification in which a firm is defined as being treated if it is owned by a Japanese parent company.11 The effect on these firms can be inferred from the differential impact of the variable of interest relative to a control group, which soaks up common seasonal patterns and other demand-driven factors in the U.S. market. While there are a number of competing methodologies for this type of estimation, we use normalized propensity score re-weighting due to the relatively favorable finite-sample properties as discussed in Busso et al. (2014), as well as for its transparent intuition. Consistent estimation of the average treatment effect on the treated requires the assumption of conditional independence: the treatment/control allocation is independent of potential outcomes conditional on a set of variables. As the average Japanese firm differs considerably from other firms in the data, we use other multinational firms - both US and non-Japanese foreign- as our baseline control group prior to reweighting. To compute the propensity scores for reweighting, we control for size and industry, which ensures the control group has a similar industrial composition and size distribution as our treated sample.12 Table 1 reports summary values for the sample, including statistics on the balancing procedure using the normalized propensity score.

The magnitude of the shock for a representative Japanese multinational is captured by the effect on total imported intermediate inputs at a monthly frequency.13 Including non-Japanese imported intermediates is important for applying the control group as a counterfactual, and the shares by source-country gives the necessary variation for identification: as shown in Table 1 the share of imported inputs from Japan is 70% of the total for Japanese firms and only 3.5% for non-Japanese multinationals. Let $$ V^{M}_{i,t}$$ be the value of intermediate imports of firm i in month t, after removing a firm-specific linear trend through March 2011. We fit the following regression:

$$\displaystyle V^{M}_{i,t} = \alpha_{i}+ \sum_{p=-4}^{9}\gamma_{p}E_{p}+ \sum_{p=-4}^{9}\beta_{p}E_{p}\mathrm{JPN}_{i,p}+ u_{i,t}$$ (5)

where $$ \alpha_{i}$$ are firm fixed-effects, $$ \gamma_{p}$$ are monthly fixed effects (with the indicator variables Ep corresponding to the calendar-months surrounding the event), and $$ u_{i,t}$$ is an error term. The baseline sample will consist of January 2009 to December 2011. We denote March 2011 as t=0.

The $$ \beta_{p}$$ coefficients are of primary interest. The $$ \mathrm{JPN}_{i,t}$$ is an indicator variable equal to one if the firm is owned by a Japanese parent company. Interacting these indicator variables with each month of the panel allows for a time-varying effect of Japanese ownership on a firm's overall intermediate input imports, particularly during and after the Tōhoku event. The $$ \beta_{p}$$ coefficients will estimate the differential effect of the Tōhoku event on Japanese multinational affiliates in the U.S., compared to the control group of non-Japanese firms. A useful interpretation of the $$ \{E_{p}\mathrm{JPN}_{i,p}\}$$ variables is as a set of instruments that captures the exogeneity of imports during these months, reflecting the source-country share of imports from Japan as evident in Table 1. To evaluate the differential impact on production for Japanese firms, we simply replace the dependent variable in equation (5) with the firm's North American exports, denoted $$ V^{NA}_{i,t}$$ .

It is important to highlight that equation (5) is in levels. There are several reasons for doing so, as opposed to using log differences or growth rates. First, allowing for the presence of zeros is important when the data are at a monthly frequency, particularly given the magnitude of the shock to imports for Japanese firms. The second reason is more conceptual. Because we are interested in calculating the average effect of these firms that represents (and can scale up to) the aggregate impact on the U.S. economy, it is appropriate to weight the firms based on their relative size. The levels specification does exactly this: the absolute deviations from trend will be greater for the bigger firms and hence will contribute disproportionately to the coefficient estimates.14 We discuss aggregation and weighted vs unweighted elasticity estimates in section 4.1.

In addition to the Conditional Independence Assumption highlighted earlier, the $$ \beta_{p}$$ coefficients are valid estimates of the mean effect for Japanese affiliates only in so far as the control group is not itself impacted by the shock. This Stable Unit Treatment Value Assumption (SUTVA) implies that general equilibrium effects or peer effects (e.g. strategic interaction) do not meaningfully effect the estimates. The share of imported inputs from Japan is low for the control group, and thus the shock is unlikely to have a measurable effect on imported inputs as a whole. We discuss the potential for strategic interaction in Appendix C.3.


2.2 Results: Total Manufacturing Sector

The top panel of Figure 5 plots the $$ \beta_{p}$$ coefficients from equation (5) for the months surrounding the Tōhoku event. Relative to the control group, there is a large drop in total intermediate input imports by Japanese firms in the months following the earthquake. The drop in intermediate inputs bottoms out at 4 million USD in t=3 (June 2011) and the point estimates do not return back to the pre-shock trend until month t=7 (October 2011).

More interesting are the results from panel B of Figure 5, which looks for evidence of the production/sales impact of this shock on Japanese firms via their North American exports. The differential time-path of N.A. exports also exhibits a substantial drop following the Tōhoku event, hitting a trough of 2 million USD below baseline in t=2 (May 2011). The standard errors, which are clustered at the firm level, are themselves interesting. As made clear via the 95-percent confidence bands on the point estimates of Figure 5, the standard errors increase dramatically in the months following the shock, a feature we interpret to reflect heterogeneous incidence and timing of the shocks (as well as the recoveries) for the Japanese multinationals.

To gain a sense of the average percentage drops of these two data series for Japanese multinationals as a group, we take the two plots of the differential dollar amounts from Figure 5 and divide by the average pre-shock level for these firms (see Table 1). The results, plotted jointly in Figure 6, show the fraction below pre-shock trend levels for these firms, on average. There is a remarkable correlation between these two series - whereby there is essentially a one-for-one drop in output for a given drop in intermediate imports. Using the mapping from Result 1, these reduced form results suggest a production function that is essentially Leontief in the imported input.

One potential concern with the interpretation of these results is separating out the intermediate input channel with other channels, such as a direct "productivity shock" affecting the U.S. operations of Japanese affiliates. Separating an ownership channel from an imported input channel is difficult due to lack of substantial overlap evident in Figure 4: few Japanese firms have low input exposure and few non-Japanese firms have high input exposure. In appendix C.7 we present results using a binary response model, as one attempt to disentangle the defining features of the import and output disruptions during this time.


3 Structural Estimation of Cross Country Input Linkages

The relative movements of imported inputs and output of Japanese multinational firms point to little substitutability of intermediate inputs. This section expands our analysis by structurally estimating the production function of firms affected by the Tōhoku shock. Unlike in the previous section, which used a set of instruments based on the differential import shares of Japanese intermediates, this estimation relies on leveraging the high degree of exogenous variation in Japanese inputs coming from the Tōhoku event, while also expanding the production function under study. This estimation serves multiple purposes. First, it is reassuring to find elasticities that are consistent with the heuristic evidence implied by our reduced-form results, when imposing a conventional production function framework. Second, by adding further structure, we can distinguish two elasticities: one between Japanese material inputs and other material inputs, and another between an aggregate bundle of material inputs and domestic capital/labor. Finally, by using an estimation procedure not relying on a control group we obtain separate estimates for Japanese and non-Japanese firms.

3.1 Framework

The estimation procedure will utilize information from two distinct periods: the six months preceding and the six months following the March 11 event. The pre-period, which we denote by $$ \tau-1$$ , yields information on the production function of the firm under profit-maximizing conditions. In the post-period, denoted $$ \tau$$ , we do not impose that the firm is optimizing over its input use, due to the fact that shipments from Japan are to some extent beyond the control of the firm. We assume that the firm's technology in any period t is given by the nested CES aggregate

$$\displaystyle x_{i,t}=\phi _{i}\left[ \mu_i^{\frac{1}{\zeta }}\left(K_{i,t}^{\alpha }L_{i,t}^{1-\alpha }\right) ^{\frac{\zeta -1}{\zeta }}+\left( 1-\mu_i \right) ^{\frac{1}{\zeta }}M_{i,t}^{\frac{\zeta -1}{\zeta }}\right] ^{\frac{\zeta }{\zeta -1}},$$ (6)

where
$$\displaystyle M_{i,t}=\left( \nu _{i}^{\frac{1}{\omega }}\left(m_{i,t}^{-J}\right) ^{\frac{\omega -1}{\omega }}+\left( 1-\nu _{i}\right) ^{\frac{1}{\omega }}\left( m_{i,t}^{J}\right) ^{\frac{\omega -1}{\omega }}\right) ^{\frac{\omega }{\omega -1}}.$$ (7)

In this production function $$ x_{i,t}$$ , $$ K_{i,t}$$ , and $$ L_{i,t}$$ denote the output, capital, and labor of firm i. The variable $$ M_{i,t}$$ denotes an aggregate of intermediate inputs of materials sourced from Japan ( $$ m_{i,t}^{J}$$ ) and materials sourced from all places other than Japan ( $$ m_{i,t}^{-J}$$ ), including domestic materials. We are interested in estimating $$ \omega$$ and $$ \zeta$$ , which parameterize the substitutability between Japanese and non-Japanese materials and that between the capital-labor aggregate and the aggregate of intermediate inputs. The parameters $$ \mu_i$$ and $$ \nu _{i}$$ are firm-specific weights and $$ \phi_{i}$$ parameterizes the firm's productivity, all of which we assume are constant over such a short time horizon. Further, we assume that the firm is monopolistically competitive and faces a CES demand function
$$\displaystyle p_{i,t}^{x}=\left( \frac{Y_{i,t}}{x_{i,t}}\right) ^{\frac{1}{\varepsilon}}.$$ (8)

As usual, $$ Y_{i,t}$$ is the bundle used or consumed downstream and serves as a demand shifter beyond the control of the firm.

3.1.1 Pre-Tsunami period

Period $$ \tau$$ corresponds to the period April-September 2011, and $$ \tau-1$$ the period September 2010 - February 2011. We exclude the month of March 2011. In period $$ \tau-1$$ the firm operates in a standard environment, choosing capital, labor, and materials to maximize

$$\displaystyle p_{i,\tau -1}^xx_{i,\tau -1} - w_{\tau -1}L_{i,\tau -1} - R_{\tau -1}K_{i,\tau -1} - p_{i,\tau -1}^{-J}m_{i,\tau -1}^{-J} - p_{i,\tau -1}^{J}m_{i,\tau -1}^{J}$$    

subject to (6), (7), and (8). The firm takes all factor prices as given. Material prices $$ p_{i,\tau -1}^{J}$$ and $$ p_{i,\tau -1}^{-J}$$ are firm-specific to indicate that different firms use different materials. It is straightforward to show that this optimization problem implies
$$\displaystyle K_{i,\tau-1}$$ $$\displaystyle = \frac{\alpha }{1-\alpha}\frac{w_{\tau-1}L_{i,\tau-1}}{R_{\tau-1}},$$ (9)
$$\displaystyle \nu_{i}$$ $$\displaystyle = \frac{\left(p_{i,\tau-1}^{-J}\right)^\omega m_{i,\tau-1}^{-J}}{\left(p_{i,\tau-1}^{J}\right)^\omega m_{i,\tau-1}^{J} + \left(p_{i,\tau-1}^{-J}\right)^\omega m_{i,\tau-1}^{-J}},$$ (10)

$$\displaystyle \mu_{i}$$ $$\displaystyle =\frac{\left( \left( \frac{R_{\tau-1}}{\alpha }\right) ^{\alpha }\left( \frac{w_{\tau-1}}{1-\alpha }\right) ^{1-\alpha }\right) ^{\zeta}K_{i,\tau-1}^{\alpha }L_{i,\tau-1}^{1-\alpha }} {\left( P_{i,\tau-1}^M\right) ^{\zeta}M_{i,\tau-1}+\left( \left( \frac{R_{\tau-1}}{\alpha }\right) ^{\alpha }\left( \frac{w_{\tau-1}}{1-\alpha }\right) ^{1-\alpha }\right) ^{\zeta}K_{i,\tau-1}^{\alpha }L_{i,\tau-1}^{1-\alpha }},$$ (11)

where
$$\displaystyle P_{i,\tau-1}^M=\left[ \nu _{i}\left( p_{i,\tau-1}^{-J}\right) ^{1-\omega }+\left(1-\nu _{i}\right) \left( p_{i,\tau-1}^{J}\right) ^{1-\omega }\right] ^{\frac{1}{1-\omega }}.$$    

We will use these relationships in the structural estimation that follows below.

3.1.2 Post-Tsunami period

At the beginning of period $$ \tau$$ many firms' production processes in Japan are disrupted. Obtaining the desired amount of shipments of materials from Japan may either be prohibitively expensive or simply impossible. Modeling firm behavior in this environment therefore requires modifications to the previous setup. One possibility is to assume that the quantity of materials that firms obtain from Japan is exogenous and that firms freely choose non-Japanese materials, capital and labor. This option is unattractive for two reasons. First, due to existing contracts it is unlikely that a firm is able to adjust the quantities of non-Japanese materials, capital, and labor without costs in such a short time frame. One remedy would be to add adjustment costs to the model. Although straightforward, this approach would require us to estimate additional parameters. Second, and more importantly, the materials sourced from Japan ($$ m_{i,t}^J$$ ) may not be exogenous for every firm. Some suppliers in Japan may have been unaffected by the earthquake and tsunami such that materials could be shipped as desired. Hence, using this approach would require us to distinguish between firms whose supply chains are disrupted and those whose are not. That is, we would have to classify firms based on an endogenous outcome.

For these reasons we prefer an alternative approach, namely to estimate the production function without specifying the full optimization problem. We only assume that in period $$ \tau$$ , firms operate the same technologies given by (6) and (7), and that no firm adjusts its capital stock such that $$ K_{i,\tau}=K_{i,\tau-1}$$ .15 Conditional on knowing the time-invariant features of the production function ($$ \phi_i$$ , $$ \mu_i$$ , $$ \nu_i$$ ), we next describe an estimation procedure that allows us to find the elasticity parameters most consistent with the observed input choices and output evident in the data.

3.2 Estimation

Recall that we use North American exports as a proxy for a firm's output $$ p_{i,t}^xx_{i,t}$$ , with the underlying assumption that the former is proportional to the latter. We continue here in the same spirit, though we now make this assumption explicit. Let $$ V_{i,t}^{NA}$$ be the value of North American exports at time t and define

$$\displaystyle \kappa _{i}=\frac{V_{i,\tau-1}^{NA}}{p_{i,\tau-1}^{x}x_{i,\tau-1}}.$$ (12)

In words, $$ \kappa _{i}$$ is the fraction of firm i's shipments exported to Canada and Mexico in the six months preceding the tsunami. We next make two assumptions that allow us to construct an estimation equation. First, we assume that a relationship analogous to (12) continues to hold in period $$ \tau$$ , except for a log-additive error $$ u_{i,\tau}$$ . That is,
$$\displaystyle \ln V_{i,\tau}^{NA}= \ln \kappa _{i}p_{i,\tau}^{x}x_{i,\tau} + u_{i,\tau}.$$ (13)

The second assumption is that $$ \mathbb{E}\left[u_{i,\tau} \vert X_{i} \right]=0$$ where Xi is a vector of all right-hand-side variables. Setting the conditional mean of $$ u_{i,\tau}$$ to zero is a standard exogeneity assumption requiring that, loosely speaking, the error is uncorrelated with all right-hand-side variables. It rules out, for example, that in response to a fall in Japanese intermediate imports firms export a fraction of their shipments to Canada and Mexico that systematically differs from $$ \kappa _{i}$$ . We provide evidence in section 4.3.1 that demonstrates that this is a reasonable assumption.

Using equation (6) we can rewrite (13) as

$$\displaystyle \ln \left( V_{i,\tau}^{NA}\right) =\ln \left( \kappa _{i}\phi _{i}\right) + \ln \left(p_{i,\tau}^x \left[ \mu_i ^{\frac{1}{\zeta }} \left( K_{i,\tau}^{\alpha }L_{i,\tau}^{1-\alpha }\right)^{\frac{\zeta -1}{\zeta }}+\left( 1-\mu_i \right) ^{\frac{1}{\zeta }}\left(M_{i,\tau}\right)^{\frac{\zeta -1}{\zeta }}\right] ^{\frac{\zeta }{\zeta -1}}\right) + u_{i,\tau}.$$ (14)

Values for $$ \nu_i$$ and $$ \mu_i$$ are obtained from equations (10) and (11).16 Using (12), the intercept can be constructed from the previous period
$$\displaystyle \kappa_{i}\phi_{i}=\frac{V_{i,\tau-1}^{NA}}{p_{i,\tau-1}^x\left[ \mu_i^{\frac{1}{\zeta }}\left(K_{i,\tau-1}^{\alpha }L_{i,\tau-1}^{1-\alpha }\right)^{\frac{\zeta -1}{\zeta }} + \left( 1-\mu_i \right)^{\frac{1}{\zeta }}\left(M_{i,\tau-1}\right)^{\frac{\zeta -1}{\zeta }}\right]^{\frac{\zeta }{\zeta -1}}}.$$    

Notice that $$ \kappa _{i}$$ and $$ \phi_{i}$$ are not separately identified. Under standard assumptions, we can consistently estimate equation (14) using, e.g., nonlinear least squares. The only parameters to calibrate are the rental rate of capital $$ R_\tau$$ and the capital share in the capital/labor aggregate $$ \alpha$$ . We estimate the two elasticities, $$ \zeta$$ and $$ \omega$$ . Notice that $$ \omega$$ appears in the intermediate aggregate $$ M_{i,\tau}$$ as shown in equation (7). The estimates ( $$ \hat{\zeta}$$ , $$ \hat{\omega}$$ ) solve $$ \displaystyle \min_{\{\zeta, \omega\}}\sum\limits_{i=1}^N (u_{i,\tau})^2$$ .

Why do we restrict the sample to the year surrounding the Tōhoku event? To understand this, recall that a principal difficulty of estimating production functions lies in unobserved inputs and productivity. Since both are unobserved by the econometrician, they are absorbed into the error term. However, because they are known to the firm, other input choices depend on them. Hence, right-hand-side variables and the error term will generally be correlated, rendering estimates inconsistent.17

By restricting the sample period to a single 12-month interval, the assumption of constant firm productivity seems appropriate. If productivity is constant, it cannot be correlated with the error term, thereby ruling out one of the concerns.18 The fact that the Tōhoku event was an unexpected shock negates much of the concern about endogeneity arising from unobserved inputs. To see why, consider the case when the firm anticipates an input disruption in a future period. Firm adjustment of unobserved inputs in expectation of this shock will impact input choices - leading to an endogeneity problem where inputs are correlated with the shock. Put simply, the unexpected nature of the Tōhoku event works towards equalizing the information sets between the econometrician and the firm because factor choices are not affected prior to the shock being realized.19

Before turning to the data we briefly discuss the intuition of parameter identification. Unlike other approaches to estimating elasticities of substitution (e.g. Feenstra et al. (2014)), our method does not rely on the response of relative values to a change in relative prices.20 In fact, in an econometric sense, our approach treats all inputs as independent variables.

A simple example illustrates how the parameters are identified. Consider the production function (6) and suppose that, for a particular firm, the initial period yields a value of $$ (1-\mu)=0.4$$ . The elasticity $$ \zeta$$ determines how deviations from this measure of the optimal input mix between the intermediate aggregate $$ M_{i,\tau}$$ and the capital labor aggregate translate into measured output. Thus, if we observe comparatively fewer intermediates $$ M_{i,\tau}$$ , reflecting a different mix of inputs than that given by 0.4, we obtain an elasticity estimate for $$ \zeta$$ that best matches the response in output. Because the estimates for $$ \mu $$ , $$ \nu$$ , and $$ \kappa_i\phi_i$$ are themselves functions of the elasticities, this procedure must iterate across the parameter space to find the estimate most consistent with the data. Similar reasoning applies for the identification of the $$ \omega$$ elasticity based on relative movements in Japanese materials, non-Japanese materials, and output. The estimates we obtain are the best fit across the firms in each sample.

3.3 Connecting Model and Data

Estimation of the model requires data on employment, Japanese and non-Japanese material inputs, exports to North America, and output prices for periods $$ \tau-1$$ and $$ \tau$$ . Since data on firm-specific capital stocks are hard to obtain and likely noisy, we use equation (9) to construct it from firm payroll and a semi-annual rental rate of 7 percent for period $$ \tau-1$$ .21 Recall that the capital stock is not adjusted over this time horizon so that $$ K_{i,\tau}=K_{i,\tau-1}$$ . The parameter $$ \alpha$$ is calibrated to $$ 1/3$$ .22

Quarterly employment information comes from the Business Register, which we adjust to reflect the average value over the 6 month periods we study, as they do not align with the quarters defined within a calendar year.23 As discussed in earlier sections, the LFTTD contains firm-level data of Japanese imports and North American exports. For non-Japanese material inputs, we would ideally combine the non-Japanese imported materials with information on domestic material usage for these firms. As information on domestic material inputs is not available in Census data at this frequency, we utilize information on the total material expenditures from the Census of Manufacturers (CM) to construct a firm-level scaling factor to gross up non-Japanese intermediate imports. Put differently, we impute non-Japanese material inputs from non-Japanese input imports. For each firm, we construct the scaling factor as

$$\displaystyle \frac{P^M_iM_i - p^J_im^J_{i}}{p^{-J}_im^{-J}_{i}}$$ (15)

from the latest CM year. Because the closest available CM year is 2007 in our data, there is some concern about missing or outdated information for this factor. We mitigate this by using industry-specific means for missing values, and winsorizing large outliers at the 90th/10th percentiles.

Regarding information on prices, the LFTTD records the value and quantity of each trade transaction (at the HS10 level), and thus it is possible to construct the associated price, or "unit-value" of each shipment directly.24 Aggregating up these shipments into a firm-month observation is complicated, of course, by the differing quantity units. Lacking any better alternative, we simply average the transaction prices using the dollar value of each transaction as weights.

Finally, we restrict the sample of firms to those that have regular imports from Japan and non-Japan over the periods we study, as well as regular North American exports.25 While this substantially limits the number of firms in each sample, the shares of trade represented by these firms in each category remains very high (see Table 2).

We obtain standard errors using bootstrap methods, which also allow us to account for the uncertainty implied by the imputation of non-Japanese material inputs. We draw randomly with replacement from our set of firms to construct 5000 distinct bootstrap samples. For each of these samples, the non-Japanese materials share is imputed as described above prior to estimation.

3.4 Summary of Results

The results of the estimation are shown in Table 2. The elasticity between material inputs for Japanese affiliates is 0.2, while the elasticity between the aggregate material input and capital/labor is 0.03. Together, these estimates are indeed consistent with the reduced-form evidence for the ($$ \psi $$ ) elasticity from section 2.2. The relative magnitudes are also intuitive: while Japanese imported inputs are strong complements with other material inputs -- consistent with the high share of intra-firm transactions comprising this trade -- there is even less scope for substitution between material inputs and domestic capital/labor.

The estimation procedure also allows us to estimate these elasticities for two samples of non-Japanese firms: non-Japanese multinationals and non-multinational firms. While the estimates for the $$ \zeta$$ elasticity are indeed very close for these other samples, the elasticity estimates corresponding to material inputs are higher, at 0.6 and 0.4 respectively. The lower share of intra-firm imports from Japan for the non-Japanese multinationals aligns with the argument that this type of trade is the key source of non-substitutability in the short-run. On the other hand, the low estimates for non-multinational firms, which have essentially zero intra-firm imports, may point to other mechanisms at work beyond the role of intra-firm trade. More generally, however, the estimates for these parameters are all significantly lower than what is commonly assumed (typically unity or higher) in the literature.

Although the number of firms included in this estimation is small (550 firms in total across the three subgroups), they account for a large share of economic activity in the United States. Looking at their combined share of total trade, these firms account for over 80% of Japanese intermediate imports, 68% of non-Japanese intermediate imports, and well over 50% of North American exports. Such high concentration of trade among relatively few firms is consistent with other studies using this data (see Bernard et al. (2007)).


4 Discussion

The structural estimates of the model are broadly in agreement with the evidence in section 2.2: imported inputs are strong complements with other inputs in production. The rigidity of the production function for multinational firms in particular is likely due to i) the high degree of intra-firm trade in what is presumably highly specialized inputs, and ii) a high-degree of supplier concentration. Our results have a number of important implications for how we should think about multinational firms, as well as aggregate topics such as volatility and business cycle co-movement.


4.1 Aggregation

Before relating our estimates to macroeconomic topics, it is necessary to discuss aggregation. Indeed, in any study utilizing micro-level estimates to inform macro-level objects of interest, the details of aggregation and heterogeneity are of critical importance. Work by Imbs & Méjean (2015) argues that imposing homogeneity across sectors when estimating consumption elasticities can be overly restrictive, creating a heterogeneity bias which can be quantitatively large. In our case one could discuss aggregation along various dimensions: across products, industries, firms, and so on. We examine the effects of product-level aggregation in Appendix C.2.

A primary concern is how to translate the results from the firm-level subsamples into estimates that would pertain to macro-oriented models. As a first step, the final column in Table 2 shows the elasticity estimates when aggregating across all firms in the sample. The results are consistent with the estimates by subgroup, suggesting substantial complementarities across inputs. All estimates in Table 2, however, correspond to the average across firms in each group, and do not take into account heterogeneity in firm size within the groups. It is relatively straightforward to modify our estimation procedure to weight firms according to their relative size.26 We report the results from this modified estimation in Panel B of Table 3. When comparing the results to those in Table 2, it is evident that the weighted estimates are not substantially different than the unweighted estimates. Although the samples of firms comprising these estimates do not amount to the total manufacturing sector of the United States, they do account for the considerable majority of U.S. trade.

4.2 Implications

The rigid production networks of foreign-owned multinationals will have direct consequences on destination and host economies. Previous literature has hypothesized that input linkages could generate business-cycle comovement, but supportive empirical evidence has been difficult to find. This paper can be seen as a first step in establishing empirical evidence for a causal relationship between trade, multinational firms, and business cycle comovement. In a companion paper (Boehm et al. (2014)), we evaluate the quantitative importance of such complementarities of imported inputs by multinational affiliates. When separately accounting for intermediate input trade by multinationals and traditional trade in final goods, the model distinguishes between the production elasticity of imported inputs and the traditional "Armington" elasticity used to bundle together international goods for consumption. The complementarities in import linkages by multinationals increases value-added comovement in the model by 11 percentage points relative to a benchmark without such firms.

This model shares similarities with several other existing models, particularly Burstein et al. (2008). A key advantage of Boehm et al. (2014), however, is a tight link to Census data for matching other features of multinationals and trade. Johnson (2014) also looks at the role of vertical linkages on comovement, but applies greater input-output structure on the model. Such features will generate increases in value-added comovement in his model, the magnitude of which becomes significant only when the elasticity of substitution among inputs is sufficiently low. Other work also identifies multinationals as a key source of the transmission of shocks: Cravino & Levchenko (2015) demonstrates that foreign multinational affiliates can account for about 10 percent of aggregate productivity shocks.27

The low value for $$ \omega$$ indicates the presence of spillovers beyond the immediate effect from Japan. That is, imports from non-Japanese locations are lower as a result of the shock in Japan28, and we would presume this applies to suppliers within the United States as well. Specifically, upstream suppliers (in countries other than Japan as well as within the U.S.) were affected indirectly via their exposure to those firms with direct exposure to Japanese inputs, combined with the rigidity of their production with respect to those inputs. Downstream suppliers that rely on the inputs from the disrupted firms would likewise be adversely affected. The presence of such spillovers combined with the large network of input linkages can indeed magnify the total effect of the transmission of the shock to the U.S. market. Such effects are also evident in a related paper, Carvalho et al. (2014), which finds large spillovers in both upstream and downstream firms in Japan following the 2011 earthquake.

Another branch of literature on the diversification of risk has studied whether firms using complex production structures with several intermediates could be less volatile (Koren & Tenreyro (2013)). Kurz & Senses (2013) establish that firms with substantial imports and exports have lower employment volatility than domestic firms in the medium to long term, which they attribute partly to the diversification of risk.29 The key result in this paper points to a possibly overlooked fact: the extent of the benefits from diversification depends heavily on the substitutability of inputs. Conditional on a given number of inputs used in production, a firm will likely experience greater volatility if each input is key to the production process and inputs are subject to heterogeneous shocks.30 Conceptually, an increase in the use of imported inputs should not be viewed necessarily as diversification. A fragmentation of production can lead to an increased supply chain risk that is an important counterweight to whatever efficiencies such complex input sourcing might afford, particularly when the production elasticities are low.

The rigid production networks of multinational firms also influences our understanding of why firms segment production across countries. In a related paper, Flaaen (2014) shows that despite the presence of substantial and complex import linkages with the source country (consistent with a vertical framework of FDI), the motive for multinational production appears to be to serve the domestic market (consistent with the horizontal framework of FDI). The result could be called "horizontal FDI with production sharing."31 Evidence of strong complementarities in this production sharing, however, presents a puzzle. Why does the firm replicate only select portions of the supply chain, considering the penalties for disruptions and mismatched inputs are so great? It is perhaps the case that the segments of the production chain that remain in the source country have a location-specific component that is not easily transferable when the firm moves production abroad.32 A greater understanding of these sourcing decisions is an area for future research.

4.3 Robustness


4.3.1 Mis-measurement of Firm Production

A natural concern with our analysis is the use of N.A. exports as a proxy for firm-level production. Perhaps it is the case that export shipments fall disproportionately more than domestic shipments following a shock to production. If this were true, the N.A. exports would indeed be a poor proxy for production, and its usefulness in evaluating a production elasticity substantially compromised.

To evaluate this concern, we narrow our study to the automotive sector, which has data on production, sales, and inventory at a monthly frequency. Using the Ward's electronic databank, we obtain plant-level information on production, and model-line information on inventory and sales.33 The specification is identical to equation (5), where the dependent variable is now Qjit: production of plant j of firm i in month t. The Japanese multinational firms are, in this case, those automakers with plants located in North America but whose parent company is headquartered in Japan.34

Figure 7 shows the results, where we once again divide by pre-shock levels to gain a sense of the percentage effects of these changes. Relative to their U.S. counterparts, Japanese automakers in the United States experienced large drops in production following the Tōhoku event. Production bottomed out in May of 2011 -- two months after the event -- at almost 60 percent below trend.35 The point estimates return to a level near zero in September of 2011, implying that the shock affected production for nearly 6 months.36 We interpret these results to be largely supportive of the results obtained using the exports-based proxy for production. The percentage drops in the two series are remarkably similar: a trough of 59% at t=2 in the automotive data vs 53% at t=2 using the proxy. We conclude that, at least for this exercise, the proxy appears to be providing valuable information on a firm's U.S. production behavior.37

4.3.2 Other Robustness Work

We discuss a number of other results and extensions to this work in the online appendix. We explore evidence for input inventories prior to the shock, evidence for product-level heterogeneity, movements in other domestic inputs, import prices movements, and exports behavior back to Japan.

4.4 External Validity

Finally, we discuss the external validity of this result. The exogenous variation we use to identify this elasticity is tied to a particular event in time, making generalization subject to some caveats. On the other hand, there are few, if any, estimates of this parameter in the existing literature. The critical question is whether the mechanisms underlying the elasticity estimates are operative beyond the circumstances surrounding this event study.

The pattern of strong intermediate input linkages with the source country is not restricted to Japanese affiliates only. As shown in Flaaen (2014), over 45 percent of the imports for all foreign multinational affiliates are sourced from the country of the parent firm. The cost share of imported intermediates from the source country is 0.12 for all foreign affiliates, which is lower than the 0.22 for Japanese affiliates but still much larger than the representative importing firm in the United States. The cost share of all imported inputs is actually quite close: 35 percent for Japanese affiliates vs 32 percent for all foreign affiliates.

A related concern is whether the estimates for Japanese affiliates are driven solely by the automotive sector. The ideal check would be to run industry-by-industry subgroup estimates for the elasticities, thereby generating heterogeneity that could be assessed relative to expectations. Unfortunately, the small number of firms applicable for this analysis, combined with disclosure requirements associated with Census Bureau data usage, prevents this degree of detail. Instead, we address this concern by splitting the sample into a motor vehicle and non-motor vehicle subsample. We do this for the Japanese multinationals as well as the total sample. The results for these four subsamples are reported in Panel C of Table 3. We find that the low elasticity estimates are not driven exclusively by firms in the motor-vehicle sector.

When viewed in light of the substantial fraction of intra-firm imports comprising multinational affiliate trade, the low elasticity of substitution should not come as a surprise. One would not expect close substitutes for the sort of specialized products reflecting firm-specific knowledge that likely comprises this trade. Moreover, such a low estimate for an elasticity of this nature is not without precedent.38 Using different methodologies, recent work by Atalay (2014) highlights strong complementarities between intermediate inputs, using industry-level data for the United States.39

Any elasticity estimate is tied to the time-horizon to which it corresponds. Ruhl (2008) emphasizes the difference between elasticities implied by responses to temporary vs permanent shocks. Larger values are calculated for an elasticity following a permanent shock, owing in part to firm responses along the extensive margin. In our context, we estimate the elasticity subject to a short-lived shock where the structure of the supply chain is plausibly fixed and extensive margin movements of supplier relationships would not apply. For this reason the elasticity parameters ( $$ \omega, \zeta$$ ) should likely generalize to other contexts of this horizon and for shocks of this general duration. Even for a long-lived shock, the estimated elasticities would remain relevant while the firm makes changes to its network of suppliers. Evaluating whether there is evidence for long-term supply-chain reorganization following the Tōhoku event is an area of ongoing work.

5 Conclusions

Using a novel firm-level dataset to analyze firm behavior surrounding a large exogenous shock, this paper reveals the mechanisms underlying cross-country spillovers. We find complementarities in the global production networks of Japanese affiliates, such that the U.S. output of these firms declined dramatically following the Tōhoku earthquake, roughly in line with the decline in imported inputs. The elasticity of substitution between imported and domestic inputs that would best match this behavior is very low - nearly that implied by a Leontief production function. The reliance on intra-firm imports by multinational affiliates from their source country is the most plausible explanation for such strong complementarities in production. Structural estimates of disaggregated elasticities are similarly low, and imply spillovers to upstream and downstream firms in the U.S. and abroad. The large impacts to Japanese affiliates together with the propagation to other U.S. firms explains the large transmission of the shock to the U.S. economy in the aggregate.

These elasticities play a critical role in the way international trade impacts both source and destination economies. Such complementarities between domestic and foreign goods have been shown to improve the ability of leading theoretical models to fit key moments of the data. We emphasize here the distinction between substitutability between domestic and foreign final goods (a "consumption" elasticity of substitution, or the so-called Armington elasticity) and substitutability between domestic and foreign intermediate goods (a "production" elasticity of substitution). In a companion paper (Boehm et al. (2014)), we document the behavior of a model with such complementarities in imported intermediates, and discuss how these elasticity parameters interact. Calibrating this model to the share of multinational affiliate trade in intermediates yields an increase in value-added comovement of 11 p.p.

Such rigid production networks will also play a role in aggregate volatility, productivity growth and dispersion, and the international ownership structure of production. The novel datasets described in this paper may help to shed light on these and other areas of research in the future.


Bibliography

D. A. Ackerberg, et al. (2006).
`Structural Identification of Production Functions'.
Working paper, University of Michigan.
E. Atalay (2014).
`How Important Are Sectoral Shocks?'.
Working Paper.
D. K. Backus, et al. (1994).
`Dynamics of the Trade Balance and the Terms of Trade: The J-Curve?'.
American Economic Review 84(1):84-103.
D. K. Backus, et al. (1995).
`International Business Cycles: Theory and Evidence'.
In T. F. Cooley (ed.), Frontiers of Business Cycle Research. Princeton University Press, Princeton.
A. Bernard, et al. (2007).
`Firms in International Trade'.
Journal of Economic Perspectives 21(3):105-130.
A. Bernard, et al. (2006).
`Transfer Pricing by U.S.-Based Multinational Firms'.
Working paper, Yale School of Management.
C. Boehm, et al. (2014).
`Complementarities in Multinational Production and Business Cycle Dynamics'.
Working paper, University of Michigan.
A. Burstein, et al. (2008).
`International Trade, Production Sharing and the Transmission of Business Cycles'.
Journal of Monetary Economics 44(4):775-795.
M. Busso, et al. (2014).
`New Evidence on the Finite Sample Properties of Propensity Score Matching and Reweighting Estimators'.
Review of Economics and Statistics forthcoming.
V. M. Carvalho, et al. (2014).
`Supply Chain Disruptions: Evidence from the Great East Japan Earthquake'.
Working paper.
F. Caselli, et al. (2014).
`Diversification through Trade'.
Working paper.
J. Cravino & A. Levchenko (2015).
`Multinational Firms and International Business Cycle Transmission'.
Working paper, University of Michigan.
J. di Giovanni & A. Levchenko (2012).
`Country Size, International Trade, and Aggregate Fluctuations in Granular Economies'.
Journal of Political Economy 120(6):1083-1132.
J. di Giovanni, et al. (2014).
`Firms, Destinations, and Aggregate Fluctuations'.
Econometrica 82(4):1303-1340.
R. C. Feenstra, et al. (2014).
`In Search of the Armington Elasticity'.
Working Paper 20063, NBER.
A. Flaaen (2013).
`Constructing a Bridge Between the LexisNexis Directory of Corporate Affiliations and the U.S. Business Register'.
Working paper, University of Michigan.
A. Flaaen (2014).
`Multinational Firms in Context'.
Working paper, University of Michigan.
J. A. Frankel & A. K. Rose (1998).
`The Endogeneity of the Optimum Currency Area Criteria'.
Economic Journal 108(449):1009-1025.
X. Gabaix (2011).
`The Granular Origins of Aggregate Fluctuations'.
Econometrica 79:733-772.
P. Goldberg, et al. (2010).
`Imported Intermediate Inputs and Domestic Product Growth: Evidence from India'.
Quarterly Journal of Economics 125(4):1727-1767.
L. Halpern, et al. (2015).
`Imported Inputs and Productivity'.
American Economic Review forthcoming.
J. Heathcote & F. Perri (2002).
`Financial Autarky and International Business Cycles'.
Journal of Monetary Economics 49:601-627.
A. Hirschmann (1958).
`Interdependence and Industrialization'.
In The Strategy of Economic Development. Yale University Press, New Haven.
J. Imbs (2004).
`Trade, Finance, Specialization, and Synchronization'.
Review of Economics and Statistics 86(3):723-734.
J. Imbs & I. Méjean (2015).
`Elasticity Optimism'.
American Economic Journal: Macroeconomics forthcoming.
R. Jarmin & J. Miranda (2002).
`The Longitudinal Business Database'.
mimeo, available at http://www.vrdc.cornell.edu/info7470/2007/Readings/jarmin-miranda-2002.pdf.
R. C. Johnson (2014).
`Trade in Intermediate Inputs and Business Cycle Comovement'.
American Economic Journal: Macroeconomics 6(4):39-83.
R. C. Johnson & G. Noguera (2012).
`Accounting for Intermediates: Production Sharing and Trade in Value Added'.
Journal of International Economics 86(2):224-236.
C. I. Jones (2011).
`Intermediate Goods and Weak Links in the Theory of Economic Development'.
American Economic Journal: Macroeconomics 3(2):1-28.
W. Keller & S. Yeaple (2013).
`The Gravity of Knowledge'.
American Economic Review 103(4):1414-1444.
M. Koren & S. Tenreyro (2013).
`Technological Diversification'.
American Economic Review 103(1):378-414.
M. Kremer (1993).
`The O-Ring Theory of Economic Development'.
Quarterly Journal of Economics 108(3):551-575.
P. Krishna & A. Levchenko (2014).
`Comparative Advantage, Complexity, and Volatility'.
Journal of Economic Behavior and Organization 94:314-329.
C. J. Kurz & M. Z. Senses (2013).
`Importing, Exporting and Firm-Level Employment Volatility'.
Working Paper 2013-44, Board of Governors of the Federal Reserve System.
W. Leontief (1936).
`Quantitative Input and Output Relations in the Economic Systems of the United States'.
The Review of Economics and Statistics 18(3):105-125.
J. Peek & E. S. Rosengren (1997).
`The International Transmission of Financial Shocks: The Case of Japan'.
American Economic Review 87(4):495-505.
J. Peek & E. S. Rosengren (2000).
`Collateral Damage: Effects of the Japanese Bank Crisis on Real Activity in the United States'.
American Economic Review 90(1):30-45.
J. Pierce & P. Schott (2012).
`A Concordance Between U.S. Harmonized System Codes and SIC/NAICS Product Classes and Industries'.
Journal of Economic and Social Measurement 37(1-2):61-96.
V. A. Ramey (1989).
`Inventories as Factors of Production and Economic Fluctuations'.
American Economic Review 79(3):338-354.
N. Ramondo, et al. (2015).
`Intra-firm Trade and Vetical Fragmentation in U.S. Multinational Corporations'.
Working paper, Columbia Business School.
K. J. Ruhl (2008).
`The International Elasticity Puzzle'.
Working paper, NYU Stern School of Business.
M. Schnell & D. E. Weinstein (2012).
`Evaluating the Economic Response to Japan's Earthquake'.
Working Paper 12-P-003, RIETI.
R. Stern, et al. (1976).
`Price Elasticities in International Trade'.
In The Strategy of Economic Development. Macmillan, London.
D. J. Wald, et al. (2006).
`ShakeMap Manual: Technical Manual, User's Guide and Software Guide'.
Tech. rep., U.S. Geological Survey.
E. Walker (1997).
`The Census Bureau Business Register: Basic Features and Quality Issues'.
paper presented at the Joint Statistical Meetings, Anaheim, CA.
N. Wasi & A. Flaaen (2014).
`Record Linkage Using Stata: Pre-Processing, Linking and Reviewing Utilities'.
Working paper, University of Michigan.


Table 1: Summary Statistics

Source: LFTTD, DCA, and UBP as explained in the text.

Panel A data are for year 2007. Panel B reports the baseline average values of N.A. exports and intermediate input imports, as well as the characteristics of that trade, for the two groups of firms: Japanese affiliates and other multinational firms. The statistics are calculated in the three months prior to the Tōhoku earthquake: Dec. 2010, Jan. 2011, and Feb 2011. The control group of other multinational firms has been re-weighted using the normalized propensity score, from a specification including the level of N.A. exports, int imports, and industry dummies. The final three columns report balancing tests of the equality of the means between the treated and control group.

Panel A: Cost Share Of Imported Inputs

  Japanese Firms Non Multinationals
from Japan 21.8 1.0
from all countries 35.0 17.5

Panel B: Treatment Effects Sample Details

  Japanese Firms Non Multinationals Balancing Tests: t Balancing Tests: $$ p>\lvert t \rvert$$ % Reduct $$ \lvert \mathrm{bias} \rvert$$
N.A. Exports 3,504,894 3,413,058 0.38 0.706 79.1
N.A. Exports share intra-firm 72.0 52.2      
Intermediate Input Imports 8,075,893 7,596,761 0.87 0.384 88
Intermediate Input Imports share from Japan 70.0 3.5      
Intermediate Input Imports share intra-firm 86.0 21.7      
Industry (Avg) - - 0.009 0.965 97.8




Table 2: Firm-Level Estimation: Results and Sample Details

Source: CM, LFTTD, DCA, and UBP as explained in the text.

This table reports the results from the firm-level estimation detailed in section 3. Panel A outlines the parameters that are calibrated prior to estimation. The top two rows of Panel B reports the point estimates of the elasticities, and the corresponding 95 percent confidence intervals using a bootstrapping procedure. (See Appendix C.8 for more details on the measurement of dispersion for these estimates.) Rows 3 and 4 report other estimates related to the calculated production functions. The final rows of Panel B describe features of the estimation samples.

Panel A: Calibration

Parameter Value
Rt 0.07
$$ \alpha$$ $$ 1/3$$

Panel B: Estimation Results

  Japanese Multinationals Non-Japanese Multinationals Non-Multinationals All Firms
$$ \omega$$ 0.201 0.624 0.423 0.552
$$ \omega$$ [0.02 0.43] [0.16 0.69] [0.26 0.58] [0.21 0.62 ]
$$ \zeta$$ 0.032 0.038 0.032 0.037
$$ \zeta$$ [0.030 0.673] [0.035 0.508] [0.029 1.68] [0.034 0.038]
Weight on K/L Aggregate ($$ \bar{\mu}$$ ) 0.223 0.514 0.278 0.409
Weight on JPN Materials ( $$ 1-\bar{\nu}$$ ) 0.173 0.044 0.147 0.096
Number of Firms 105 304 141 550
Share of Total Trade: JPN int imports 0.60 0.23 0.03 0.86
Share of Total Trade: Non-JPN int imports 0.02 0.66 0.01 0.69
Share of Total Trade: N.A. exports 0.08 0.47 0.01 0.56




Table 3: Firm-Level Estimation: Other Results

Source: CM, LFTTD, DCA, and UBP as explained in the text.

This table reports additional estimation results. Panel B recalculates the results from Table 2 using a vector of weights to assign larger firms a greater share in the estimation. Panel C divides the samples based on the motor vehicle industry.

Panel A: Calibration

Parameter Value
Rt 0.07
$$ \alpha$$ $$ 1/3$$

Panel B: Estimation Results (Weighted)

  Japanese Multinationals Non-Japanese Multinationals Non-Multinationals All Firms
$$ \omega$$ 0.157 0.611 0.543 0.606
$$ \omega$$ [0.02 0.40] [0.30 1.23] [0.305 0.57] [0.28 0.70]
$$ \zeta$$ 0.241 0.038 0.032 0.037
$$ \zeta$$ [0.03 0.884] [0.034 0.51] [0.029 0.55] [0.034 0.038 ]
Number of Firms 105 304 141 550

Panel C: Estimation Results: MV Sector

  Japanese Mult. Motor Vehicles Japanese Mult. Non-Motor Vehicles All Firms Motor Vehicles All Firms Non-Motor Vehicles
$$ \omega$$ 0.311 0.094 0.414 0.574
$$ \omega$$ [0.019 0.398] [0.016 0.59] [0.27 0.60] [0.16 0.66]
$$ \zeta$$ 0.032 0.071 0.037 0.037
$$ \zeta$$ [0.030 0.48] [0.028 1.27] [0.031 0.64] [0.033 0.037]
Number of Firms 35 70 100 450



Figure 1: Index of Japanese Industrial Production: Manufacturing Jul.2010 - Jan.2012
Source: Japanese Ministry of Economy, Trade, and Industry (METI). The series are logged, HP-Filtered, after seasonally adjusting.

See link for accessible version.

Accessible version



Figure 2: U.S. Imports from Japan and Rest of World, Jul.2010 - Jan.2012
Source: U.S. Census Bureau (FT900: U.S. International Trade in Goods and Services). The series are logged, HP-Filtered, after seasonally adjusting.

See link for accessible version.

Accessible version



Figure 3: U.S. Industrial Production: Manufacturing and Durable Goods
Source: Federal Reserve Board, Industrial Production and Capacity Utilization - G.17 Series . Series is Seasonally Adjusted.

See link for accessible version.

Accessible version



Figure 4: Density of Firm-Level Exposure to Japanese Imported Inputs: By Firm Type
Source: LFTTD-DCA-UBP as explained in text. The estimates correspond to year 2010. This figure displays density estimates of the log exposure measure to Japanese imported inputs, separately for Japanese affiliates and non-Japanese multinational firms. The measure is defined as the ratio of Japanese imported inputs to total imported inputs plus U.S. salaries and wages. Estimates at either tail are suppressed for confidentiality purposes.

See link for accessible version.

Accessible version



Figure 5: Dynamic Treatment Effects: Japanese Firms
Source: LFTTD-DCA-UBP as explained in text. These figures report the intermediate imports and North American exports of the U.S. affiliates of Japanese firms relative to a control group of other multinational firms. The values are coefficient estimates taken from an interaction of a Japanese-firm dummy with a monthly dummy - additional baseline monthly dummies remove seasonal effects. See equation 5 in the text. Standard errors are clustered at the firm level.

See link for accessible version.

Accessible version



Figure 6: Relative Imported Inputs and Output (Proxy) of Japanese Firms: Fraction of Pre-Shock Level
Source: LFTTD-DCA-UBP as explained in text. This figure reports the intermediate imports and output proxy (North American exports) of the U.S. affiliates of Japanese firms relative to a control group of other multinational firms. The values are percent changes from the pre-shock level of each series, defined as the average of the months December 2010, January 2011, and February 2011.

See link for accessible version.

Accessible version



Figure 7: Assessing the Output Proxy Using Monthly Automotive Production
Source: Ward's Automotive Database. This figure reports the production levels of Japanese auto plants relative to a control group of non-Japanese auto plants. The values are percent changes from a pre-shock level, defined as the average of the months December 2010, January 2011, and February 2011. See equation A16 in the text. For purposes of comparison, we also include the equivalent measure corresponding to total manufacturing of Japanese affiliates using the output proxy from Census data (from Figure 6). The Japanese automakers are Honda, Mazda, Mitsubishi, Nissan, Toyota, and Subaru. For the sake of clarity, we suppress the standard errors for the automotive series, though there are 4 months with below zero production based on a 95 percent confidence interval. See Appendix C.10 for more details.

See link for accessible version.

Accessible version



A. Basic Theory Appendix


A.1 Proof of Result 1

Suppose that the firm solves

$$\displaystyle \max p_{x}x-p_{D}F_{D}-p_{M}IM $$
subject to

$$\displaystyle x=\left[ \left( 1-\mu \right) ^{\frac{1}{\psi }}\left[ F_{D}\right] ^{\frac{\psi -1}{\psi }}+\mu ^{\frac{1}{\psi }}\left[ IM\right] ^{\frac{\psi -1}{\psi}}\right] ^{\frac{\psi }{\psi -1}}$$
and

$$\displaystyle p_{x}=\left( \frac{Y}{x}\right) ^{\frac{1}{\varepsilon }} $$
The first order conditions are

$$\displaystyle \left( 1-\frac{1}{\varepsilon }\right) \left( Y\right) ^{\frac{1}{\varepsilon }}\left( x\right) ^{\frac{1}{\psi }-\frac{1}{\varepsilon }}\left( 1-\mu \right) ^{\frac{1}{\psi }}\left[ F_{D}\right]^{-\frac{1}{\psi }}=p_{D} $$

$$\displaystyle \left( 1-\frac{1}{\varepsilon }\right) \left( Y\right) ^{\frac{1}{\varepsilon }}\left( x\right)^{\frac{1}{\psi }-\frac{1}{\varepsilon }}\mu ^{\frac{1}{\psi }}\left[ IM \right] ^{-\frac{1}{\psi }}=p_{M} $$
Dividing one by the other gives

$$\displaystyle \frac{F_{D}^{\ast }}{IM^{\ast }}=\frac{1-\mu }{\mu }\left( \frac{p_{M}}{p_{D}}\right) ^{\psi }. $$
The same equation can be obtained under perfect competition.

Now take the production function and multiply it by px

$$\displaystyle p_{x}x=p_{x}\left[ \left( 1-\mu \right) ^{\frac{1}{\psi }}\left[ F_{D}\right]^{\frac{\psi -1}{\psi }}+\left( p_{M}\right) ^{-\frac{\psi -1}{\psi }}\mu ^{\frac{1}{\psi }}\left[ p_{M}IM\right] ^{\frac{\psi -1}{\psi}}\right] ^{\frac{\psi }{\psi -1}} $$
Taking logs gives
$$\displaystyle \ln \left( p_{x}x\right) =$$ $$\displaystyle \frac{\psi }{\psi -1}\ln \left( p_{x}\left[\left( 1-\mu \right) ^{\frac{1}{\psi }}\left[ F_{D}\right] ^{\frac{\psi -1}{\psi }}+\left( p_{M}\right) ^{-\frac{\psi -1}{\psi }}\mu^{\frac{1}{\psi }}\left[ p_{M}IM\right] ^{\frac{\psi -1}{\psi }}\right] \right)$$ (A1)
$$\displaystyle =$$ $$\displaystyle \frac{\psi }{\psi -1}\ln \left( p_{x}\left[ \left( 1-\mu \right) ^{\frac{1}{\psi }}\exp \left( \frac{\psi -1}{\psi }\ln \left[ F_{D}\right] \right)+\left( p_{M}\right) ^{-\frac{\psi -1}{\psi }}\mu ^{\frac{1}{\psi }}\exp\left( \frac{\psi -1}{\psi }\ln \left[ p_{M}IM\right] \right) \right] \right)$$ (A2)

Before differentiating, recall the assumption that the firm takes prices pM as given and that it cannot change px after learning about the shock. Then

$$\displaystyle \frac{\partial \ln p_{x}x}{\partial \ln p_{M}M} =$$ $$\displaystyle \frac{\psi }{\psi -1}\frac{p_{x}\left( P_{M}\right) ^{-\frac{\psi -1}{\psi }}\mu ^{\frac{1}{\psi } }\exp \left( \frac{\psi -1}{\psi }\ln \left[ p_{M}IM\right] \right) \frac{\psi -1}{\psi }}{p_{x}\left[ \left( 1-\mu \right) ^{\frac{1}{\psi }}\exp \left( \frac{\psi -1}{\psi }\ln \left[ F_{D}\right] \right) +\left(p_{M}\right) ^{-\frac{\psi -1}{\psi }}\mu ^{\frac{1}{\psi }}\exp \left( \frac{\psi -1}{\psi }\ln \left[ p_{M}IM\right] \right) \right] }$$ (A3)
$$\displaystyle =$$ $$\displaystyle \frac{1}{1+\left( \frac{1-\mu }{\mu }\right) ^{\frac{1}{\psi }}\left[ \frac{F_{D}}{IM}\right] ^{\frac{\psi -1}{\psi }}}$$ (A4)

We evaluate this elasticity at

$$\displaystyle \frac{F_{D}^{\ast }}{IM}=\frac{IM^{\ast }}{IM}\frac{1-\mu }{\mu }\left( \frac{ p_{M}}{p_{D}}\right) ^{\psi } $$
so that

$$\displaystyle \frac{\partial \ln p_{x}x}{\partial \ln p_{M}IM}=\frac{1}{1+\left( \frac{IM^{\ast }}{IM}\right) ^{\frac{\psi-1}{\psi }}\frac{1-\mu }{\mu }\left( \frac{p_{M}}{p_{D}}\right) ^{\psi -1}} $$

A.2 On Flexibility in Domestic Inputs

Under the assumption of perfect competition, the first order conditions are:

$$\displaystyle x\left( 1-\mu \right) =\left( p_{D}\right) ^{\psi }F_{D} $$

$$\displaystyle x\mu =\left( p_{M}\right) ^{\psi }IM $$
If the firm takes prices px, pM, and pD as given, the following elasticities are immediate:

$$\displaystyle \frac{\partial \ln \left( p_{x}x\right) }{\partial \ln \left( p_{D}F_{D}\right) }=\frac{\partial \ln \left( p_{x}x\right) }{\partial \ln \left( p_{M}M\right) }=\frac{\partial \ln \left( p_{D}F_{D}\right) }{\partial \ln \left( p_{M}M\right) }=1. $$
The above equations demonstrate that a constant returns to scale production function combined with these assumptions on market structure imply that the output elasticity will equal one for all values of the elasticity of substitution. For this reason, we require some assumptions limiting the flexibility of domestic inputs following the import disruption.

Below we show an alternative way of understanding the interaction of competitive factor markets, changes in domestic inputs, and the mapping of the output elasticity into parameter values for the elasticity of substitution. Consider the total derivative of $$ \ln(x)$$ :

$$\displaystyle d\ln x = \frac{\partial \ln x}{\partial IM}d \ln IM + \frac{\partial \ln x}{\partial F}d \ln F$$ (A5)

$$\displaystyle d\ln x = \frac{ \mu^{\frac{1}{\psi}} (IM)^{\frac{\psi-1}{\psi}}d\ln IM}{( 1-\mu)^{\frac{1}{\psi}} \left[ F_D\right]^{\frac{\psi -1}{\psi }}+ \mu^{\frac{1}{\psi}} \left[IM\right]^{\frac{\psi -1}{\psi }}} + \frac{(1- \mu)^{\frac{1}{\psi}} (F_D)^{\frac{\psi-1}{\psi}} d \ln F_D }{( 1-\mu)^{\frac{1}{\psi}} \left[ F_D\right]^{\frac{\psi -1}{\psi }}+ \mu^{\frac{1}{\psi}}\left[IM\right]^{\frac{\psi -1}{\psi }}}$$ (A6)

Dividing by $$ d\ln IM$$ yields:

$$\displaystyle \frac{d\ln x}{d\ln IM} = \frac{ \mu^{\frac{1}{\psi}} (IM)^{\frac{\psi-1}{\psi}}}{( 1-\mu)^{\frac{1}{\psi}} \left[ F_D\right]^{\frac{\psi -1}{\psi }}+ \mu^{\frac{1}{\psi}} \left[IM\right]^{\frac{\psi -1}{\psi }}} + \frac{(1- \mu)^{\frac{1}{\psi}} (F_D)^{\frac{\psi-1}{\psi}}}{( 1-\mu)^{\frac{1}{\psi}} \left[ F_D\right]^{\frac{\psi -1}{\psi }}+ \mu^{\frac{1}{\psi}}\left[IM\right]^{\frac{\psi -1}{\psi }}} \frac{d \ln F_D }{d \ln IM}$$    

Now, as before, combining the first order conditions from the profit maximization problem, we have:

$$\displaystyle \frac{F_D(\cdot)}{IM} = \frac{1-\mu}{\mu}\left(\frac{p_D}{p_M}\right)^{-\psi}$$ (A7)

Log-differentiating this expression:
$$\displaystyle d \ln\left(\frac{F_D}{IM}\right)$$ $$\displaystyle = -\psi d \ln \left(\frac{p_D}{p_M}\right)$$    
$$\displaystyle d \ln F_D - d \ln IM$$ $$\displaystyle = -\psi d \ln \left(\frac{p_D}{p_M}\right)$$    
$$\displaystyle \frac{d \ln F_D}{d \ln IM}$$ $$\displaystyle = 1 - \psi \frac{d \ln \left(\frac{p_D}{p_M}\right)}{d \ln IM}$$ (A8)

Finally, we have:

$$\displaystyle \frac{d\ln x}{d\ln IM} = \frac{ \mu^{\frac{1}{\psi}} (IM)^{\frac{\psi-1}{\psi}}}{( 1-\mu)^{\frac{1}{\psi}} \left[ F_D\right]^{\frac{\psi -1}{\psi }}+ \mu^{\frac{1}{\psi}} \left[IM\right]^{\frac{\psi -1}{\psi }}} + \frac{(1- \mu)^{\frac{1}{\psi}} (F_D)^{\frac{\psi-1}{\psi}}\left[1 - \psi \frac{d \ln \left(\frac{p_D}{p_M}\right)}{d \ln IM}\right]}{( 1-\mu)^{\frac{1}{\psi}} \left[ F_D\right]^{\frac{\psi -1}{\psi }}+ \mu^{\frac{1}{\psi}}\left[IM\right]^{\frac{\psi -1}{\psi }}}$$ (A9)

Thus, if there is no change in the relative input price following the disruption in IM of the firm: $$ \frac{d \ln \left(\frac{p_M}{p_D}\right)}{d \ln IM}=0$$ , then the output elasticity will be equal to one regardless of the value of $$ \psi $$ . On the other hand, any assumptions that yield a non-zero change in the relative input prices will then yield the result that $$ \frac{d\ln x}{d\ln IM}=1$$ provided $$ \psi \rightarrow 0$$ .


B. Data Appendix


B.1 Matching Corporate Directories to the Business Register

The discussion below is an abbreviated form of the full technical note (see Flaaen (2013)) documenting the bridge between the DCA and the Business Register.

B.1.1 Directories of International Corporate Structure

The LexisNexis Directory of Corporate Affiliations (DCA) is the primary source of information on the ownership and locations of U.S. and foreign affiliates. The DCA describes the organization and hierarchy of public and private firms, and consists of three separate databases: U.S. Public Companies, U.S. Private Companies, and International - those parent companies with headquarters located outside the United States. The U.S. Public database contains all firms traded on the major U.S. exchanges, as well as major firms traded on smaller U.S. exchanges. To be included in the U.S. Private database, a firm must demonstrate revenues in excess of $1 million, 300 or more employees, or substantial assets. Those firms included in the International database, which include both public and private companies, generally have revenues greater than $10 million. Each database contains information on all parent company subsidiaries, regardless of the location of the subsidiary in relation to the parent.

The second source used to identify multinational firms comes from Uniworld Business Publications (UBP). This company has produced periodic volumes documenting the locations and international scope of i) American firms operating in foreign countries; and ii) foreign firms with operations in the United States. Although only published biennially, these directories benefit from a focus on multinational firms, and from no sales threshold for inclusion.

Because there exist no common identifiers between these directories and Census Bureau data infrastructure, we rely on probabilistic name and address matching -- so-called "fuzzy merging" -- to link the directories to the Census data infrastructure.

B.1.2 Background on Name and Address Matching

Matching two data records based on name and address information is necessarily an imperfect exercise. Issues such as abbreviations, misspellings, alternate spellings, and alternate name conventions rule out an exact merging procedure, leaving the researcher with probabilistic string matching algorithms that evaluate the "closeness" of match -- given by a score or rank -- between the two character strings in question. Due to the large computing requirements of these algorithms, it is common to use so-called "blocker" variables to restrict the search samples within each dataset. A "blocker" variable must match exactly, and as a result this implies the need for a high degree of conformity between these variables in the two datasets. In the context of name and address matching, the most common "blocker" variables are the state and city of the establishment.

The matching procedure uses a set of record linking utilities described in Wasi & Flaaen (2014). This program uses a bigram string comparator algorithm on multiple variables with differing user-specified weights.40 This way the researcher can apply, for example, a larger weight on a near name match than on a perfect zip code match. Hence, the "match score" for this program can be interpreted as a weighted average of each variable's percentage of bigram character matches.

B.1.3 The Unit of Matching

The primary unit of observation in the DCA, UBP, and BR datasets is the business establishment. Hence, the primary unit of matching is the establishment, and not the firm. However, there are a number of important challenges with an establishment-to-establishment link. First, the DCA (UBP) and BR may occasionally have differing definitions of the establishment. One dataset may separate out several operating groups within the same firm address (i.e. JP Morgan - Derivatives, and JP Morgan - Emerging Markets), while another may group these activities together by their common address. Second, the name associated with a particular establishment can at times reflect the subsidiary name, location, or activity (i.e. Alabama plant, processing division, etc), and at times reflect the parent company name. Recognizing these challenges, the primary goal of the matching will be to assign each DCA (UBP) establishment to the most appropriate business location of the parent firm identified in the BR. As such, the primary matching variables will be the establishment name, along with geographic indicators of street, city, zip code, and state.

B.1.4 The Matching Process: An Overview

The danger associated with probabilistic name and address procedures is the potential for false-positive matches. Thus, there is an inherent tension for the researcher between a broad search criteria that seeks to maximize the number of true matches and a narrow and exacting criteria that eliminates false-positive matches. The matching approach used here is conservative in the sense that the methodology will favor criteria that limit the potential for false positives at the potential expense of slightly higher match rates. As such, the procedure generally requires a match score exceeding 95 percent, except in those cases where ancillary evidence provides increased confidence in the match.41

This matching proceeds in an iterative fashion, in which a series of matching procedures are applied with decreasingly restrictive sets of matching requirements. In other words, the initial matching attempt uses the most stringent standards possible, after which the non-matching records proceed to a further matching iteration, often with less stringent standards. In each iteration, the matching records are assigned a flag that indicates the standard associated with the match.

See Table A1 for a summary of the establishment-level match rate statistics by year and type of firm. Table A2 lists the corresponding information for the Uniworld data.

B.1.5 Construction of Multinational Indicators

The DCA data allows for the construction of variables indicating the multinational status of the U.S.-based establishment. If the parent firm contains addresses outside of the United States, but is headquartered within the U.S., we designate this establishment as part of a U.S. multinational firm. If the parent firm is headquartered outside of the United States, we designate this establishment as part of a Foreign multinational firm. We also retain the nationality of parent firm.42

There can be a number of issues when translating the DCA-based indicators through the DCA-BR bridge for use within the Census Bureau data architecture. First, there may be disagreements between the DCA and Census on what constitutes a firm, such that an establishment matches may report differing multinational indicators for the same Census-identified firm. Second, such an issue might also arise due to joint-ventures. Finally, incorrect matches may also affect the degree to which establishment matches agree when aggregated to a firm definition. To address these issues, we apply the following rules when using the DCA-based multinational indicators and aggregating to the (Census-based) firm level. There are three potential cases:43

Potential 1: A Census-identified firm in which two or more establishments match to different foreign-country parent firms

  1. Collapse the Census-identified firm employment based on the establishment-parent firm link by country of foreign ownership
  2. Calculate the firm employment share of each establishment match
  3. If one particular link of country of foreign ownership yields an employment share above 0.75, apply that link to all establishments within the firm.
  4. If one particular link of country of foreign ownership yields an employment share above 0.5 and total firm employment is below 10,000, then apply that link to all establishments within the firm.
  5. All other cases require manual review.

Potential 2: A Census-identified firm in which one establishment is matched to a foreign-country parent firm, and another establishment is matched to a U.S. multinational firm.

  1. Collapse the Census-identified firm employment based on the establishment-parent firm link by type of DCA link (Foreign vs U.S. Multinational)
  2. Calculate the firm employment share of each establishment match
  3. If one particular type of link yields an employment share above 0.75, apply that link to all establishments within the firm.
  4. If one particular type of link yields an employment share above 0.5 and total firm employment is below 10,000, then apply that link to all establishments within the firm.
  5. All other cases require manual review.

Potential 3: A Census-identified firm in which one establishment is matched to a non-multinational firm, and another establishment is matched to a foreign-country parent firm (or U.S. multinational firm).

Apply same steps as in Potential 2.


B.2 Classifying Firm-Level Trade

The firm-level data on imports available in the LFTTD does not contain information on the intended use of the goods.44 Disentangling whether an imported product is used as an intermediate input for further processing -- rather than for final sale in the U.S. -- has important implications for the nature of FDI, and the role of imported goods in the transmission of shocks. Fortunately, the Census Bureau data contains other information that can be used to distinguish intermediate input imports from final goods imports. Creating lists of the principal products produced by firms in a given detailed industry in the United States should indicate the types of products that, when imported, should be classified as a "final" good - that is, intended for final sale without further processing. The products imported outside of this set, then, would be classified as intermediate goods.45 Such product-level production data exists as part of the "Products" trailer file of the Census of Manufacturers. As detailed in Pierce & Schott (2012) (see page 11), combining import, export, and production information at a product-level is useful for just such a purpose.

B.2.1 Creating a NAICS-Based set of Final/Intermediate Products

As part of the quinquennial Census of Manufacturers (CM), the Census Bureau surveys establishments on their total shipments broken down into a set of NAICS-based (6 digit) product categories. Each establishment is given a form particular to its industry with a list of pre-specified products, with additional space to record other product shipments not included in the form. The resulting product trailer file to the CM allows the researcher to understand the principal products produced at each manufacturing establishment during a census year.

There are several data issues that must be addressed before using the CM-Products file to infer information about the relative value of product-level shipments by a particular firm. First, the trailer file contains product-codes that are used to "balance" the aggregated product-level value of shipments with the total value of shipments reported on the base CM survey form. We drop these product codes from the dataset. Second, there are often codes that do not correspond to any official 7-digit product code identified by Census. (These are typically products that are self-identified by the firm but do not match any of the pre-specified products identified for that industry by Census.) Rather than ignoring the value of shipments corresponding to these codes, we attempt to match at a more aggregated level. Specifically, we iteratively try to find a product code match at the 6, 5, and 4 digit product code level, and use the existing set of 7-digit matches as weights to allocate the product value among the 7-digit product codes encompassing the more aggregated level.

We now discuss how this file can be used to assemble a set of NAICS product codes that are the predominant output (final goods) for a given NAICS industry. Let xpij denote the shipments of product p by establishment i in industry j during a census year. Then the total output of product p in industry j can be written as:

$$\displaystyle X_{pj} = \sum_{i=1}^{I_j} x_{pij} ,$$
where Ij is the number of firms in industry j. Total output of industry j is then:

$$\displaystyle X_{j} = \sum_{p=1}^{P_j} X_{pj} .$$
The share of industry output accounted for by a given product p is therefore:

$$\displaystyle S_{pj} = \frac{X_{pj}}{X_j}. $$

One might argue that the set of final goods products for a given industry should be defined as the set of products where $$ S_{pj}>0$$ . That is, a product is designated as a "final good" for that industry if any establishment recorded positive shipments of the product. The obvious disadvantage of employing such a zero threshold is that small degrees of within-industry heterogeneity will have oversized effects on the classification.

Acknowledging this concern, we set an exogenous threshold level W such that any p in a given j with $$ S_{pj}>W$$ is classified as a final good product for that industry. The upper portion of Table A3 documents the number of final goods products and the share of intermediate input imports based on several candidate threshold levels. The issues of a zero threshold are quite clear in the table; a small but positive threshold value (0.1) will have a large effect on the number of products designated as final goods. This shows indirectly that there are a large number of products produced by establishments in a given industry, but a much smaller number that comprise the bulk of total value.

There are several advantages to using the CM-Products file rather than using an input-output table.46 First, within a given CM year, the classification can be done at the firm or establishment level rather than aggregating to a particular industry. This reflects the fact that the same imported product may be used as an input by one firm and sold to consumers as a final product by another. Second, the CM-Products file is one of the principal data inputs into making the input-output tables, and thus represents more finely detailed information. Related to this point, the input-output tables are produced with a significant delay - the most recent available for the U.S. is for year 2002. Third, the input-output tables for the U.S. are based on BEA industry classifications, which imply an additional concordance (see below) to map into the NAICS-based industries present in the Census data.

We now turn to the procedure to map firm-level trade into intermediate and final goods using the industry-level product classifications calculated above.

B.2.2 Mapping HS Trade Transactions to the Product Classification

The LFTTD classifies products according to the U.S. Harmonized Codes (HS), which must be concorded to the NAICS-based product system in order to utilize the classification scheme from the CM-Products file. Thankfully, a recent concordance created by Pierce & Schott (2012) can be used to map the firm-HS codes present in the LFTTD data with the firm-NAICS product codes present in the CM-Products data.

A challenge of this strategy is that the LFTTD exists at a firm-level, while the most natural construction of the industry-level classification scheme is by establishment. More concretely, for multi-unit, multi-industry firms, the LFTTD is unable to decompose an import shipment into the precise establishment-industry of its U.S. destination. 47 While recognizing the caution that should be used in this regard, we adopt the approach that is commonly used in such circumstances: the industry of the firm is defined as that industry encompassing the largest employment share.

Once the firm-level trade data is in the same product classification as the industry-level filter created from the CM-Products file, all that is left is to match the trade data with the filter by NAICS industry. Thus, letting Mij denote total imports from a firm i (firm i is classified as being in industry j), we can then categorize the firm's trade according to:

$$\displaystyle \left. \begin{array}{l} M^{\mathrm{int}}_{ij} = \sum\limits_{p\notin P_j} M_{ipj} \\ \\ M^{\mathrm{fin}}_{ij} = \sum\limits_{p\in P_j} M_{ipj} \end{array} \right \}$$   where$$\displaystyle \qquad P_j= \left\{p \ \vert \ S_{pj} \geq W\right\}.$$ (A10)

The bottom section of Table A3 shows some summary statistics of the intermediate share of trade according to this classification system, by several values of the product-threshold W. There are at least two important takeaways from these numbers. First, the share of intermediates in total imports is roughly what is reported in the literature using IO Tables. Second, the share of total trade occupied by intermediate products is not particularly sensitive to the exogenous threshold level. While there is a small increase in the share when raising the threshold from 0 to 0.1 (about 3 percentage points), the number is essentially unchanged when raising it further to 0.2.


B.3 Sample Selection

B.3.1 Constructing the Baseline Dataset

This section will discuss the steps taken to construct the sample used in section 2.1.

Beginning with the raw files of the LFTTD export/import data, we drop any transactions with missing firm identifiers, and those pertaining to trade with U.S. territories. Next, we merge the LFTTD files with the HS-NAICS6 product concordance from Pierce & Schott (2012); if there is no corresponding NAICS6 code for a particular HS code, then we set NAICS6 equal to XXXXXX. We then aggregate up to the level of Firm-Country-Month-NAICS6, and then create extracts according to three sets of destinations/sources: Japan, Non-Japan, and North America (Canada and Mexico). Then, assigning each firm to an LBD-based industry (see below), we run the NAICS-based trade codes through the intermediate/final goods filter discussed in Appendix B.2. The firms' monthly trade can then be split into intermediate and final goods components. We repeat this step for years 2009, 2010, and 2011.

Using the Longitudinal Business Database, we drop inactive, ghost/deleted establishments, and establishments that are not in-scope for the Economic Census. To create the sample of manufacturing firms in the U.S., we first create a firm industry code defined as the industry encompassing the largest share of firm employment. We then drop non-manufacturing firms. Next, we merge the LBD for each year with the DCA-Bridge (see section B.1) containing multinational indicators. We then apply the rules specified above for clarifying disagreements with the DCA-based multinational indicators. After creating monthly copies of each firm, we merge by firm-month to the trade data. Missing information of trade data is altered to represent zeros. We repeat these steps for years 2009-2011, and then append the files together. Firms that do not exist in all three years are dropped from the sample.


B.3.2 GIS Mapping of Earthquake Intensity Measures to Affiliate Locations

As part of the Earthquake Hazards Program, the U.S. Geological Survey produces data and map products of the ground motion and shaking intensity following major earthquakes. The preferred measure to reflect the perceived shaking and damage distribution is the estimated "Modified Mercalli Intensity (MMI)" which is based on a relation of survey response and measured peak acceleration and velocity amplitudes. The USGS extends the raw data from geologic measurement stations and predicts values on a much finer grid using standard seismological inferences and interpolation methods. The result is a dense grid of MMI values covering the broad region affected by the seismic event. For more information on this methodology, see Wald et al. (2006).

To utilize this information, we take all Japanese addresses from the DCA/Uniworld directories that correspond to any U.S. operation via an ownership link. We geocode these addresses into latitude/longitude coordinates using the Google Geocoding API, and then compute the inverse distance-weighted mean of the relevant seismic intensity measure based on a 10km radius surrounding a given establishment. The firm identifiers within the corporate directories allow us to create firm-specific measures (average and maximum values, by manufacturing/non-manufacturing), which can then be brought into the baseline Census dataset via the bridges discussed in appendix B.1.


Table A1: DCA Match Statistics: 2007-2011

Total

  # of DCA Establishments Matched to B.R. Percent Matched
2007 112,346 81,656 0.73
2008 111,935 81,535 0.73
2009 111,953 81,112 0.72
2010 111,998 79,661 0.71
2011 113,334 79,516 0.70

U.S. Multinationals

  # of DCA Establishments Matched to B.R. Percent Matched
2007 22,500 16,396 0.73
2008 23,090 16,910 0.73
2009 22,076 16,085 0.73
2010 21,667 15,785 0.73
2011 21,721 15,557 0.72

Foreign Multinationals

  # of DCA Establishments Matched to B.R. Percent Matched
2007 10,331 7,555 0.73
2008 9,351 6,880 0.74
2009 11,142 8,193 0.74
2010 11,308 8,181 0.72
2011 11,619 8,357 0.72



Table A2: Uniworld Match Statistics: 2006-2011

U.S. multinationals include only the establishment identified as the U.S. headquarters.

Foreign Multinationals

  # of Uniworld Establishments Matched to B.R. Percent Matched
2006 3,495 2,590 0.74
2008 3,683 2,818 0.76
2011 6,188 4,017 0.65

U.S. Multinationals

  # of Uniworld Establishments Matched to B.R. Percent Matched
2007 4,043 3,236 0.80
2009 4,293 3,422 0.80



Table A3: Appendix Table Comparing the Results from Threshold Values W

Threshold Values

Number of Final Good Products per Industry

  W = 0 W=0.1 W=0.2
Median 19 1 1
Mean 25 1.52 1.14
Min 1 1 0
Max 154 6 3

Implied Share of Intermediate Inputs

  W = 0 W=0.1 W=0.2
Imports 60.9 63.90 63.97
Exports 52.0 54.96 55.04


Figure A1: Geographic Distribution of Earthquake Intensity and Affiliate Locations
Source: USGS and DCA/Uniworld Directories This figure plots the geographic distribution of the Tōhoku earthquake, based on recorded measurements taken directly after the event. The "Modified Mercalli Intensity" (MMI) scale is constructed based on a relation of survey response and measured peak acceleration and velocity amplitudes from prior major seismic events. Each dot corresponds to a geocoded Japanese affiliate location corresponding to a firm with U.S. operations. For more details, see Appendix B.3.2.

See link for accessible version.

Accessible version


C. Appendix: Other Results

C.1 Intermediate Input Inventories

Inventories are another obvious feature that should influence the relationship between input shipments, production, and the elasticity of substitution. In particular, inventories of intermediate inputs allow the firm to absorb unforeseen shocks to input deliveries without an impact on the production process.48 As it relates to the production elasticity, however, the presence of these inventories should serve to diminish or delay the production impact, thereby increasing the elasticity relative to what it would be without such inventories.

In fact, it is striking the extent to which we do not see any evidence for the role of intermediate input inventories in the production impacts of Figure 5 (Panel B) or Figure 7. The effect on production appears to be almost immediate, indicating that the stock of inventories of imported intermediates is low (less than one month's supply) for these firms.

We obtain a rough sense of the degree of inventory holdings from the Census of Manufacturers micro-data. Combining information on the beginning period stock of materials inventories with the annual usage of materials, we calculate the average monthly supply of inventories for each firm.49 Table A4 calculates the production-weighted averages over a select set of firm groups.50 We see that on average, Japanese multinationals hold a little over 3-weeks supply of intermediate inputs as inventory. This is slightly less than non-multinational firms, a fact that aligns with the oft-cited "lean" production processes made famous by Japanese firms in previous decades. Though these data are for the year 2007, there is little reason to believe these relative magnitudes have changed substantially over a period of a few years. For completeness, Table A4 also reports the corresponding estimates for output inventories.51

Low inventory holdings combined with an inelastic production function suggests that firms are willing to tolerate some degree of expected volatility in their production. Either the costs of holding inventories or diversifying sources of supply are sufficiently high, or firms believe the probability of disruption is low. In either case, these lean production strategies carry a greater potential for the propagation of shocks across countries, perhaps affecting firms with limited knowledge of their indirect exposure through complicated production chains.


C.2 Multi-Products and Sub-Optimal Mix

In the frameworks used in sections 2.1 and 3, we consider the aggregate bundles of imported intermediates, abstracting away from product-level detail. In reality, the firms in our dataset often import many distinct intermediate inputs from Japan. The structure of a CES production function implies that if each of these within-country inputs was non-substitutable with one another (a further, nested Leontief structure), the production impact of a disruption in the supply of just one input could be amplified relative to the value of that input.52 We evaluate this possibility below.

This is particularly true given the heterogeneous impact of the Tōhoku event across Japan (see Figure A1). This could translate into considerable dispersion in the impact on the products imported by a particular U.S. firm or Japanese affiliate. With product-level shocks, considering the effect on the aggregate import bundle amounts to assuming either 1) perfect substitutability among products, or 2) that the firm maintains an optimal within-country product mix at all times.

To be concrete, it may be more accurate to view the Mt in equation (1) as a further C.E.S. function of multiple products. Thus, we can define the proper measurement of this variable as

$$\displaystyle V^M_{i,t} = P^M_{i,t}\left(\sum_{n=1}^N \eta_n^{\frac{1}{\chi}} (m^J_{n,i,t})^{\frac{\chi-1}{\chi}} \right)^{\frac{\chi}{\chi-1}},$$ (A11)

where now $$ V^M_{i,t}$$ is the value based on a combination of N distinct products, with weights $$ \eta_n$$ and elasticity $$ \chi$$ .

Product-level heterogeneity in the production impact of the shock combined with imperfect coordination among input suppliers implies that the aggregate (measured) import bundle for a particular firm may turn out to be suboptimal. In this case, we are measuring $$ \widehat{V^M}_{i,t} = \sum_{n=1}^N (p^m_{n,t}m_{n,t}) \geq V^M_{i,t}$$ . And the lower the elasticity of substitution among products, the more severe the disconnect between the measured imports and the "effective" imports -- that which is actually useful in downstream production.

A suboptimal product mix indicates that measured imports ( $$ \widehat{V^M}_{i,t}$$ ) are greater than the effective imports ($$ V^M_{i,t}$$ ). As a result the measured output response to the import shock will be larger than otherwise, resulting in a downward "bias" in the elasticity estimates from section 2.1 and 3.53 Such an effect is decreasing in the product-level elasticity parameter $$ \chi$$ , as complementarity itself is the driving force between differences in $$ \widehat{V^M}_{i,t}$$ and $$ V^M_{i,t}$$ . In addition, the effect is also increasing in the degree of deviation from the optimal product mix.

Does this exert a quantitatively large effect on our point estimates? Given the emphasis on low inventories and lean production processes in downstream operations, one might expect that across-product adjustment would take place before sending the inputs abroad. To analyze this empirically, we analyze whether there are significant deviations in the product composition of Japanese imports during the months following the Tōhoku event. To do this, we construct a measure of the distance of a firm's import bundle from a benchmark, which we will interpret to be the optimal bundle. Let $$ t=s^*$$ be such a benchmark date. Then, using the product-level information in the LFTTD data we construct for each firm i, the share of total imports from Japan for a given product code n. Defining this share to be $$ s_{n,i,t}$$ , we then construct the average product-level distance from optimum $$ DO_{i,t}$$ as:

$$\displaystyle DO_{i,t} = \frac{1}{N^i}\sum_{n=1}^{N^i} \left(\lvert s_{n,i,t} - s_{n,i,s^*} \rvert \right)$$ (A12)

where Ni is the total number of products imported by firm i. We define the period $$ s^*$$ to be the months of April-June of 2010, and then evaluate DOi at a monthly frequency, with particular interest in the months following the Tōhoku event. While there may be natural movements in the bundle of products imported from Japan, evidence for substantial coordination failure in product composition or heterogeneity in product-level shocks would come from any abnormal jumps in this index in the months of the disruption. One can calculate this at various levels of product aggregation (i.e. HS4, HS6, HS8, HS10), though we report results using the HS6 level.54

The results of this exercise are shown in Figure A3. We plot the average DOi across Japanese firms for each month (the figure shows a 3-month moving average) during the period 2009-2011. Mechanically, this measure should be relatively close to zero in the months consisting of the benchmark (April-June 2010). While there is a secular rise in this measure on either side of this benchmark period, there do not appear to be any large jumps in the months directly following the Tōhoku event. More interesting, perhaps, are the considerably larger values for this measure during early 2009, which might reflect the effects of the trade collapse associated with the Great Recession. We interpret Figure A3 as evidence that the potential for suboptimal mix across products from Japan does not pose a serious problem to our measurement in previous sections.


C.3 Strategic Behavior

Another possibility that could affect the interpretation of the results from Figure 6 might be strategic behavior, particularly on the part of the competitors of Japanese firms in the United States. These firms could raise production or prices following the negative supply shock affecting their competitors, which would serve to bias downward the $$ \beta_{p}$$ coefficients from the equation with $$ X^{NA}_{i,t}$$ as the dependent variable.55 To evaluate this possibility, we turn to the Ward's automotive data and consider the production of non-Japanese automakers in the months directly following the Tōhoku event. Figure A4 plots the relative production of these firms, using time-series variation only. There appears to be no quantitatively meaningful responses in the months following March 2011. This should not come as a surprise given capacity constraints and utilization adjustment costs, particularly given the short time horizon of this shock.


C.4 Effects on Unit Values (Prices) of Trade

Traditionally, estimating the elasticity of substitution is accomplished via price and quantity data for products over extended periods of time. For the short horizon we consider in this paper, there are several reasons why prices may not have the scope to adjust. Many supplier relationships negotiate prices for longer periods of time than one or two months. Second, and perhaps more importantly, Table 1 demonstrates that the large majority of imported intermediate inputs are intra-firm. The observed prices of these transactions are transfer-prices (within firm) and not likely to change reflecting any short-term disturbance.

The LFTTD contains information on quantities as well as values for each trade transaction, recorded at a highly disaggregated product definition (HS 10 digit). This allows for the construction of unit values (prices) for each firm-product-month observation, which allows for an analysis of price movements surrounding the Tōhoku event.

The majority of the data construction is identical to that in section B.3, however there are a number of modifications. First, we drop all transactions with missing or imputed quantities in the LFTTD, and then aggregate to the Firm-HS10-month frequency, separately for each type of trade transaction: 1) Related-Party imports from Japan; 2) Non Related-Party imports from Japan; 3) Related-Party exports to Canada/Mexico; and 4) Non Related-Party exports to Canada/Mexico. Next, we select only those firms identified as manufacturing in the LBD. We keep the related-party and arms-length transactions separate as one may expect these prices to behave differently following a shock. As above, we keep only manufacturing firms, append the annual files together, and then select only those firms identified as a multinational in either 2009, 2010, or 2011.

At the product level, there is little reason to suspect trends or seasonal variation over this short of a time period. Moreover, there is no concern here about accounting for zeros in the data. As such we take a firm j's imports (exports) of product p in month t, and run the following specification in logs ( m$$ _{p,j,t} = log($$M$$ _{p,j,t})$$ :

m$$\displaystyle _{p,j,t} = \alpha_{pj} + \sum_{i=-19}^{9}\gamma_{i}E_{i}+ \sum_{i=-19}^{9}\beta_{i}E_{i}D_{j,i}+ u_{j,t}$$ (A13)

where $$ \alpha_{pj}$$ are firmXproduct fixed-effects, $$ \gamma_{i}$$ are monthly fixed effects (with the dummy variable $$ E_{i}'s$$ corresponding to each calendar month), and $$ u_{j,t}$$ are random effects. The variables $$ D_{j,t}$$ are dummy variables equal to one if the firm is owned by a Japanese parent company.

A qualitative version of the results is shown in Table A7. The results confirm that there are few significant price movements on import or export transactions for either Japanese or non-Japanese multinationals surrounding the Tōhoku event.


C.5 Domestic Inputs

It is also possible to evaluate the response of domestic inputs directly, using the limited information we have on quarterly firm-level employment and payroll information, taken from the Census Bureau's Business Register (BR).56

The Standard Statistical Establishment List (SSEL) contains quarterly employment and payroll information for all employers (with some small exceptions) in the U.S. economy. This list is held separately as a single-unit(SSEL-SU) and multi-unit (SSEL-MU) file. The Report of Organization Survey (ROS) asks firms to list the establishments which report under a particular EIN, and this information is then recorded to the firm identifier on the Multi-Unit File. To build a quarterly employment series at the firm-level, we link the EIN variables on the SU file with the firm-identifier linked with each EIN on the MU file. In principle, the four quarters of payroll listed on the SSEL is combined by Census to create an annual payroll figure for each establishment, which is the value recorded in the LBD. Similarly, the employment variable corresponding to the 1st quarter (week of March 12) from the SSEL is that used by the LBD.

Once we merge the SSEL-based data with quarterly employment and payroll to the LBD for a particular year, we conduct a series of reviews to ensure that the annual payroll (and 1st quarter employment) roughly align. Any establishments with disagreements between the SSEL-based payroll and LBD-based payroll such that the ratio was greater than 2 or less than 0.6 were dropped.

After these modifications were made, the remainder of the data construction was similar to that in section B.3. We merge multinational indicators from the DCA, drop non-manufacturing firms, append the 2009, 2010, and 2011 files together, and keep only those firms that exist in each year. Using the same set of firms as a control group as specified in section 2.1, we run the following regression:

$$\displaystyle \Delta$$   emp$$\displaystyle _{j,t} = \sum_{i=-3}^{3}\gamma_{i}E_{i}+ \sum_{i=-3}^{3}\beta_{i}E_{i}D_{j,i}+ u_{j,t}$$ (A14)

where $$ \Delta$$   emp$$ _{j,t} \equiv ln($$emp$$ _{j,t}/$$   emp$$ _{j,t-4})$$ , where emp$$ _{j,t}$$ indicates employment at firm j in quarter t. We also re-run the equation specified in equation A14 using payroll $$ pay_{j,t}$$ as the dependent variable (where $$ \Delta$$   pay$$ _{j,t} \equiv ln($$pay$$ _{j,t}/$$   pay$$ _{j,t-4})$$ . The qualitative results are shown in Table A6. We find no significant effects on either employment or payroll for Japanese firms in the quarter(s) following the shock. Of course, there are a number of reasons -- principally labor adjustment costs -- why one would expect little, if any, impacts on employment following this short-lived shock. Press releases dispatched by the Japanese automakers during this time indicated that no layoffs would occur. Rather, the firms indicated that they would use the production stoppages for employee skill and safety training.


C.6 Alternate Specifications for Treatment Effects Regressions

Our results from section 2.2 are based on a sample including all Japanese multinationals in manufacturing, and therefore uses a levels specification to allow for zeros in the firm-month observations. Because larger firms exhibit greater absolute deviations from trend, this roughly amounts to weighting firms based on size, such that the results correspond to a representative firm based on the aggregate effect of the group.

To see this, and to explore how the levels specification influences our interpretation, we repeat the analysis on a subset of the firms for which we can view the percentage changes directly. Specifically, we drop any firms with zeros in any month for intermediate imports or N.A. exports during the sample, and then take logs and HP-filter each series to obtain percentage deviations from trend for each firm.57 The results of this exercise are shown in Panel A of Figure A5. We suppress standard errors for the sake of clarity; the drops are significant at the 95% level for between 2-4 months following the shock. If we rerun these regressions while also weighting according to the pre-shock size of firms, we obtain a picture that looks much closer to Figure 6, see Panel B of Figure A5.

These results indicate that the larger firms appear to be affected the most from this shock. This could be partly a result of our proxy being less effective for smaller firms that may not engage in consistent exports to North America.


C.7 Probit Model of Import/Output Disruptions

We specify a simple probit model to understand the relative importance of various firm-level characteristics in the import and output declines following the tsunami. The model is

$$\displaystyle Pr(X^D_{ik}=1) =\Phi\left[ \beta_1\mathrm{JPN}_{ik} + \beta_2Exposed_{ik} + \beta_3MMI_{ik} + \beta_4Port_{ik} + \gamma_k\right]$$ (A15)

where the dependent variable ($$ X^D_{ik}$$ ) is an indicator equal to one if the N.A. exports of firm i in industry k are on average 20% below trend during the five months following the Tōhoku event. The independent variables are also indicators: $$ \mathrm{JPN_{ik}}$$ , for affiliates of Japanese multinationals; Exposedik, for firms with an exposure to Japanese inputs above 0.05 of total material; MMIik for firms with an elevated MMI value pertaining to their average Japanese manufacturing locations; and Portik for firms that typically rely on imports via ports damaged by the tsunami.58 The $$ \gamma_k$$ term allows for industry-specific intercepts. To evaluate the determinants of an input disruption from Japan, we replace the dependent variable with $$ J^D_{ik}$$ , an indicator for a drop in Japanese imported inputs of 20% relative to trend.

Panel A of Table A5 evaluates firm characteristics predicting a drop in U.S. output ($$ X^D_{ik}$$ ), as measured by our proxy. The columns (1)-(4) show the results from different specifications with various combinations of the covariates in equation (A15). Both Japanese ownership and high exposure to Japanese inputs significantly increase the probability of an output disruption, as expected. In columns (3) and (4), we demonstrate that Japanese ownership is substantially more indicative of an output decline than high input exposure alone. In Panel B, we replace the dependent variable with the binary measure of a drop in Japanese intermediate inputs ($$ J^D_i$$ ). The results from these regressions indicate, unsurprisingly, that high exposure to Japanese imports are highly predictive of a subsequent disruption following the Tōhoku event. Apart from their exposure to imports from Japan, the Japanese affiliates are no more likely to suffer a disruption to these imports (see column 8).59 While the results from Table A5 are somewhat inconclusive, they nevertheless point to unique features of the production function of Japanese affiliates that yields direct pass-through of Japanese shocks to the U.S. economy. Our estimation procedure that follows should help to clarify this point further.


C.8 Bootstrapping Standard Errors

We use bootstrapping methods to compute measures of the dispersion of our point estimates. Using random sampling with replacement within each group of firms, we create 5000 new artificial samples and re-run the estimation procedure. The standard deviation of the point estimates across these bootstrap samples is shown in Table 2. To gain a more complete picture of the dispersion, we create density estimates for each sample of firms across the parameter space for the elasticities. These densities are shown in Figure A7.


C.9 Effects on U.S. Exports to Japan

Another dimension of the transmission of the Tōhoku shock to the United States is U.S. exports back to Japan. To the extent that firms in the U.S. receive inputs from Japan for processing and re-shipment back to Japan, one might expect the U.S. exports to Japan may fall following the Tōhoku event. On the other hand, U.S. firms may have increased shipments to Japan following the shock in order to offsett what were large production and supply shortages within Japan. To evaluate this, we re-run the specification in equation (5) but replace $$ V^{M}_{i,t}$$ , the value of intermediate imports of firm i in month t, with $$ V^{JEXP}_{i,t}$$ , the value of Japanese exports of firm i in month t. The results are shown in Figure A6. As is clear from the figure, we do not see strong evidence to support either hypothesis regarding this particular trade flow, at least as it pertains to Japanese multinationals in particular.


C.10 Ward's Automotive Data

Ward's electronic databank offers a variety of data products for the global automotive industry at a monthly frequency. We obtain Japanese production (by model), North American production (by plant and model), U.S. inventory (by model), and North American sales (by model) all for the period January 2000 to December 2012. The inventory and sales data also contain the country of origin, so one can separate out these variables based on whether a particular model was imported vs domestically-produced. The series cover the universe of the assembly operations of finished cars and light trucks. Unfortunately, there is no information on input shipments.

For the plant-level analysis of production, the base sample consists of 167 plants active at some point during 2000-2012. We remove plants that were not continuously in operation during the period 2009-2012, and combine several plants that are recorded separately in the data, but are in effect the same plant. After these modifications, the sample reduces to 62 plants, 22 of which are owned by a Japanese parent. The average monthly production in the three months preceding the shock is 12,904 for Japanese plants, and 14,903 for Non-Japanese plants. The specification is identical to that in section 2.1:

$$\displaystyle Q_{i,t} = \alpha_0 + \alpha_i + \sum_{p=-14}^{9}\gamma_pE_p + \sum_{p=-14}^{9}\beta_p E_p\mathrm{JPN}_{i,p}+ u_{i,t}$$ (A16)

where here the variable $$ Q_{i,t}$$ is auto production by plant i in month t, after removing a plant-specific trend though March 2011. Because these plants can be tracked with some confidence back in time, it is reasonable here to remove seasonality directly, rather than assume a shared seasonal component between the treated and control groups as in section 2.2. We use the X12-ARIMA model, provided by the National Bank of Belgium, and apply it to each series before correcting for trend. The results for the Japanese plants are mostly similar, as shown in table A8.


Table A4: Summary Statistics: Inventories by Firm Type

Source: CM, LFTTD, DCA, and UBP as explained in the text. The data are for year 2007. This table reports the average monthly supply of inventories [(usage/12)/beginning period inventory stock] for materials and output.

  Japanese Multinationals Non Multinationals
Inputs 0.83 1.08
Output 0.31 0.45



Figure A2: Evidence of Potential Spillovers: Relative Non-Japanese Imported Inputs of Japanese Firms: Fraction of Pre-Shock Level
Source: LFTTD-DCA-UBP as explained in text. This figure reports the non-Japanese intermediate importsof the U.S. affiliates of Japanese firms relative to a control group of other multinational firms. The values are percent changes from the pre-shock level of each series, defined as the average of the months December 2010, January 2011, and February 2011. Standard errors are suppressed to report the series in percentage terms. The drops are significant at the 95 percent level for months 0, 1, 2, 3, and 5 following the earthquake.

See link for accessible version.

Accessible version


Figure A3: Japanese Products: Average Distance from Benchmark Cost Shares: JPN Multinationals
Source: LFTTD-DCA-UBP as explain in the text Underlying this figure is the calculation of the average total (absolute) deviations from a benchmark measure of a firm's cost shares across input products from Japan. See equation A12 in the text. The figure reports the mean across the Japanese multinationals used in the section 3.

See link for accessible version.

Accessible version


Figure A4: Automotive Production, Inventory, Sales by Firm Type, Distributed Lag Model
Source: Ward's Automotive Database This figure reports North American production, and U.S. sales and inventory data according to firm type: Japanese and non-Japanese firms. The values are coefficient estimates taken from a distributed lag model, exploiting time-series variation only. The underlying series have been seasonally adjusted, logged, and HP-Filtered Standard errors are suppressed in the interests of clarity. The Japanese automakers are Honda, Mazda, Mitsubishi, Nissan, Toyota, and Subaru.

See link for accessible version.

Accessible version


Figure A5: Relative Inputs and Output (Proxy) of Japanese Firms (Reduced Sample) Logged, HP-Filtered
Source: LFTTD-DCA-UBP as explained in text. These figures report the relative percentage deviations from trend of Japanese affiliates relative to a control group of other multinational firms. The values are coefficient estimates taken from an interaction of a Japanese-firm dummy with a monthly dummy - additional baseline monthly dummies remove seasonal effects. These results reflect a reduced sample with no firm-month zeros in imported inputs or N.A. exports. The data is logged, and HP-filtered using a monthly smoothing parameter.

See link for accessible version.

Accessible version


Figure A6: Dynamic Treatment Effects: Relative Japanese Exports of Japanese Firms
Source: LFTTD-DCA-UBP as explained in text. These figures report the Japanese exports of the U.S. affiliates of Japanese firms relative to a control group of other multinational firms. The values are coefficient estimates taken from an interaction of a Japanese-firm dummy with a monthly dummy - additional baseline monthly dummies remove seasonal effects. Standard errors are clustered at the firm level.

See link for accessible version.

Accessible version



Table A5: Predicting Japanese Import and U.S. Output Disruption by Firm Characteristics

*** p < 0.01 , **p <0.05 , * p <0.1
Source: LFTTD, DCA, UBP, and USGS as explained in the text. This table reports the results of a probit model prediction of JPN import and N.A. exports (output) disruption based on firm characteristics. See section 2.1 for a definition of the variables.

  (1) Panel A: Disruption to U.S. Output (proxy), $$ X_i^D=1$$ (2) Panel A: Disruption to U.S. Output (proxy), $$ X_i^D=1$$ (3) Panel A: Disruption to U.S. Output (proxy), $$ X_i^D=1$$ (4) Panel A: Disruption to U.S. Output (proxy), $$ X_i^D=1$$ (5) Panel B: Disruption to Japanese Imports, $$ J_i^D=1$$ (6) Panel B: Disruption to Japanese Imports, $$ J_i^D=1$$ (7) Panel B: Disruption to Japanese Imports, $$ J_i^D=1$$ (8) Panel B: Disruption to Japanese Imports, $$ J_i^D=1$$
Japan 0.443***   0.352*** 0.347** 0.707***   0.310*** 0.686***
Japan (0.0921)   (0.117) (0.152) (0.0917)   (0.115) (0.150)
Exposed   0.351*** 0.145 0.140   0.814*** 0.636*** 0.991***
Exposed   (0.0886) (0.112) (0.149)   (0.0880) (0.110) (0.144)
JPN*Exp       -0.00771       -0.848***
JPN*Exp       (0.228)       (0.222)
MMI -0.176*** -0.121* -0.178*** -0.178*** 0.346*** 0.389*** 0.341*** 0.306***
MMI (0.0676) (0.0646) (0.0676) (0.0683) (0.0691) (0.0667) (0.0694) (0.0704)
Ports -0.174 -0.144 -0.197   0.248 0.217 0.168 0.174
Ports (0.224) (0.225) (0.226)   (0.211) (0.212) (0.213) (0.213)
Constant -0.674 -0.674 -0.674 -0.674 -4.672 -4.672 -4.672 -4.668
Constant (0.681) (0.681) (0.681) (0.681) (85.78) (85.78) (85.78) (85.00)
Industry Dummies Yes Yes Yes Yes Yes Yes Yes Yes
Observations 2451 2451 2451 2451 2451 2451 2451 2451


Figure A7: Density Estimates of Elasticities Across Bootstrap Samples
Source: LFTTD-DCA-UBP as explained in text.

See link for accessible version.

See link for accessible version.

Accessible version



Table A6: Dynamic Treatment Effects: Quarterly Employment/Payroll Surrounding Tōhoku Event

Log 4-Quarter Difference

Source: SSEL and DCA as explained in the text.

Robust standard errors (clustered at the firmXProduct level) pertaining to each sign coefficient are indicated by: *** p<0.01, ** p<0.05, * p<0.1. This table reports qualitative features of firm employment and firm payroll in the quarters surrounding the Tōhoku earthquake and tsunami. The first set of coefficients correspond to quarter dummies, whereas the second set (JPNx) correspond to the interaction of a Japanese firm dummy with quarter dummies. See equation A14 in the text. The dependent variable is the four-quarter log difference of employment (payroll).

Independent Variables (1) Employment (2) Payroll
Q2_2010 (t=-3) pos*** pos***
Q3_2010 (t=-2) pos*** pos***
Q4_2010 (t=-1) pos*** pos***
Q1_2011 (t=0) pos*** pos***
Q2_2011 (t=1) pos*** pos***
Q3_2011 (t=2) pos*** pos***
Q4_2011 (t=3) pos*** pos***
JPNxQ2_2010 (t=-3) neg neg
JPNxQ3_2010 (t=-2) neg neg
JPNxQ4_2010 (t=-1) neg neg
JPNxQ1_2011 (t=0) neg neg
JPNxQ2_2011 (t=1) neg neg
JPNxQ3_2011 (t=2) neg neg
JPNxQ4_2011 (t=3) neg pos
constant neg*** neg***
Firm Fixed Effects Yes Yes
Observations    
R-squared    



Table A7: Dynamic Treatment Effects: Unit Values of Trade Surrounding Tōhoku Event

Source: LFTTD, DCA, and UBP as explained in the text.

Robust standard errors (clustered at the firmXProduct level) pertaining to each sign coefficient are indicated by: *** p<0.01, ** p<0.05, * p<0.1 .

This table reports qualitative features of the unit values of trade surrounding the 2011 Tōhoku earthquake and tsunami. The first set of coefficients correspond to monthly dummies, whereas the second set (JPNx) correspond to the interaction of a Japanese firm dummy with monthly dummies. See equation A13 in the text.

Log Unit-Value

Independent Variables (1) JPN Imports Related Party (2) JPN Imports Non-Related Party (3) N.A. Exports Related Party (4) N.A. Exports Non-Related Party
Sep 2010 (t=-6) neg** pos pos* pos
Oct 2010 (t=-5) pos neg pos** pos
Nov 2010 (t=-4) pos pos pos** pos
Dec 2010 (t=-3) pos neg pos pos
Jan 2011 (t=-2) neg pos neg pos
Feb 2011 (t=-1) pos neg pos** pos
Mar 2011 (t=0) neg pos pos pos
Apr 2011 (t=1) pos pos pos pos
May 2011 (t=2) neg pos neg pos**
Jun 2011 (t=3) pos** neg pos** neg
Jul 2011 (t=4) neg neg pos neg
Aug 2011 (t=5) pos pos neg pos
Sep 2011 (t=6) pos pos pos pos**
Oct 2011 (t=7) neg neg pos pos
Nov 2011 (t=8) pos neg pos neg
Dec 2011 (t=9) neg pos pos** pos
JPNxSep 2010 (t=-6) pos** neg* neg** neg
JPNxOct 2010 (t=-5) neg* pos pos pos
JPNxNov 2010 (t=-4) neg pos neg neg
JPNxDec 2010 (t=-3) neg neg* pos pos
JPNxJan 2011 (t=-2) pos neg neg neg
JPNxFeb 2011 (t=-1) neg pos pos pos**
JPNxMar 2011 (t=0) pos pos neg neg
JPNxApr 2011 (t=1) neg pos neg neg
JPNxMay 2011 (t=2) pos neg pos neg
JPNxJun 2011 (t=3) neg pos* neg neg
JPNxJul 2011 (t=4) pos neg pos neg
JPNxAug 2011 (t=5) neg* neg* neg pos
JPNxSep 2011 (t=6) neg neg neg neg
JPNxOct 2011 (t=7) pos neg neg neg
JPNxNov 2011 (t=8) neg neg neg pos
JPNxDec 2011 (t=9) neg neg pos neg
constant pos neg neg neg
FirmXProduct Fixed Effect Yes Yes Yes Yes
Observations        
R-Squared        



Table A8: Dynamic Treatment Effects: N.A. Automotive Production

Source: Ward's Automotive Yearbook

Robust standard errors (clustered at the plant level) in parentheses. *** p<0.01, ** p<0.05, * p<0.1 .

VARIABLES (1) Prod (2) Prod
Nov 2010 (t=-4) 91.06 17.78
Nov 2010 (t=-4) (649.9) (608.8)
JPN x Nov 2010 (t=-4) -195.8 -341.7
JPN x Nov 2010 (t=-4) (841.9) (799.2)
Dec 2010 (t=-3) -1,973*** 310.3
Dec 2010 (t=-3) (467.5) (497.5)
JPN x Dec 2010 (t=-3) -385.0 -408.3
JPN x Dec 2010 (t=-3) (736.5) (706.4)
Jan 2011 (t=-2) -611.5 1,083*
Jan 2011 (t=-2) (637.3) (618.7)
JPN x Jan 2011 (t=-2) 781.0 -1,092
JPN x Jan 2011 (t=-2) (792.1) (804.6)
Feb 2011 (t=-1) (401.9) (394.7)
JPN x Feb 2011 (t=-1) -1,142 -1,210*
JPN x Feb 2011 (t=-1) (696.2) (666.8)
Mar 2011 (t=0) 4,356*** 1,483***
Mar 2011 (t=0) (524.9) (389.1)
JPN x Mar 2011 (t=0) -3,515*** -2,592***
JPN x Mar 2011 (t=0) (812.0) (842.7)
Apr 2011 (t=1) -216.2 305.5
Apr 2011 (t=1) (707.7) (620.4)
JPN x Apr 2011 (t=1) -6,239*** -6,099***
JPN x Apr 2011 (t=1) (1,303) (1,282)
May 2011 (t=2) 1,584*** 799.1
May 2011 (t=2) (525.4) (511.3)
JPN x May 2011 (t=2) -7,244*** -6,625***
JPN x May 2011 (t=2) (1,651) (1,740)
Jun 2011 (t=3) 1,366** -499.3
Jun 2011 (t=3) (623.6) (594.9)
JPN x Jun 2011 (t=3) -4,564*** -3,423**
JPN x Jun 2011 (t=3) (1,248) (1,320)
Jul 2011 (t=4) -4,512*** 123.3
Jul 2011 (t=4) (878.4) (606.2)
JPN x Jul 2011 (t=4) -2,143 -3,723***
JPN x Jul 2011 (t=4) (1,430) (1,045)
Aug 2011 (t=5) 685.6 -1,323**
Aug 2011 (t=5) (744.0) (648.1)
JPN x Aug 2011 (t=5) -1,275 -1,108
JPN x Aug 2011 (t=5) (970.8) (1,012)
Sep 2011 (t=6) -836.5 -1,895***
Sep 2011 (t=6) (663.7) (641.5)
JPN x Sep 2011 (t=6) -359.4 40.37
JPN x Sep 2011 (t=6) (930.7) (959.8)
Oct 2011 (t=7) -338.0 -1,434**
Oct 2011 (t=7) (662.3) (632.4)
JPN x Oct 2011 (t=7) 93.27 -265.4
JPN x Oct 2011 (t=7) (885.6) (785.8)
Nov 2011 (t=8) -1,393** -1,443**
Nov 2011 (t=8) (582.8) (601.2)
JPN x Nov 2011 (t=8) -1,318 -2,059*
JPN x Nov 2011 (t=8) (1,159) (1,183)
Dec 2011 (t=9) -4,511*** -1,619**
Dec 2011 (t=9) (774.4) (655.5)
JPN x Dec 2011 (t=9) 759.1 24.95
JPN x Dec 2011 (t=9) (1,105) (803.9)
Constant -1,535*** -1,683***
Constant (89.30) (91.95)
Plant Fixed Effects Yes Yes
Remove Plant-Specific Pre-Shock Trend Yes Yes
Remove Seasonal Component No Yes
Observations 2,976 2,976
R-squared 0.260 0.272




Footnotes

* Boehm and Pandalai-Nayar: University of Michigan Department of Economics Flaaen: Federal Reserve Board of Governors, 20th St. and Constitution Ave, NW, Washington DC 20551, aaron.b.flaaen@frb.gov. We would like to thank Andrei Levchenko, Kyle Handley, Matthew Shapiro, and Linda Tesar for valuable comments, suggestions, and support. We also thank our discussants Rob Johnson, Eduardo Morales, Joan Monras, and Sebnem Kalemli-Ozcan, as well as seminar participants at the IMF, Michigan, Stanford, FRB, Colorado, Washington, BU, Dallas Fed, FREIT-EITI, Barcelona GSE-SI, NBER SI-ITM, SED, FREIT -RMET and Johns Hopkins-SAIS for helpful comments and suggestions. This research has received research grants from the Michigan Institute for Teaching and Research in Economics (MITRE), for which we are very thankful. Support for this research at the Michigan RDC from NSF (ITR-0427889) is also gratefully acknowledged. Any opinions and conclusions expressed herein are those of the authors and do not necessarily represent the views of the U.S. Census Bureau, the Board of Governors, or its research staff. All results have been reviewed to ensure no confidential information is disclosed. Return to Text
1. For a very useful compendium of this research from this era, see Stern et al. (1976). Return to Text
2. Due to the robust nature of these shortcomings, Backus et al. (1995) refer to them as the "price anomaly" and "quantity anomaly" respectively. Return to Text
3. Intra-firm trade accounts for a large majority of the trade of Japanese affiliates. More generally, the intra-firm share of imported intermediates for all foreign affiliates in the U.S. is 71 percent. Return to Text
4. Since 1900, the three earthquakes of greater recorded magnitude are: the 1960 Great Chilean earthquake (magnitude 9.5), the 1964 Good Friday earthquake in Prince William Sound, Alaska (magnitude 9.2); and the 2004 Sumatra-Andaman earthquake (magnitude 9.2). Return to Text
5. For precautionary reasons, all nuclear power plants were immediately shut down following the earthquake, and remained largely offline until 2015 or later. Because the electricity infrastructure exists on two separate grids (a 60Hz to the south and west, and 50Hz to the north and east), the reduction in power supply in Northeast Japan was not easily remedied, and power outages persisted for months. Return to Text
6. At the level of total U.S. GDP, both Deutsche Bank and Goldman Sachs revised 2nd quarter U.S. estimates down by 50 basis points explicitly due to the events in Japan. Return to Text
7. Note that products intended for final sale for a given industry may still be used as intermediates for other firms in a different industry. Alternatively, such "final goods" can be sold directly to consumers for ultimate consumption. Return to Text
8. Another consideration with the use of this proxy is whether it more accurately reflects production or sales, as the two are distinct in the presence of output inventories. In our case, this depends on whether the inventories are held in the U.S. or Canada/Mexico. Without further evidence, we interpret the proxy to be capturing some mix between production and sales. The structural estimation in section 3 will allow for such a mix. Return to Text
9. There is a second case which we do not examine, where $$ \psi\rightarrow \infty$$ and $$ p_{M}<p_{D}$$ and thus the firm only uses IM. We discard this scenario because such a firm would not show up in our data (i.e. this implies zero U.S. employment). Return to Text
10. The exposure measure used in Figure 4 is from 2010 and does not include the cost of domestic material usage. Return to Text
11. We could have also used a threshold of Japanese input usage for the classification of treatment status. Doing so yields estimates that are very similar, which is due to the patterns evident in Figure 4. We have also tried conditioning on our geographic information (i.e. the firm-level Japanese MMI index) in defining treatment status. The results are largely unchanged from those we report here, and for the sake of clarity we report results pertaining to the full sample. Return to Text
12. Using the predicted values (p) from the first stage regression, the inverse probability weights are $$ \frac{1}{1-p}$$ for the control group and $$ \frac{1}{p}$$ for the treated group. To normalize the weights such that the treated firms have weights equal to one, we then multiply each set of weights by p. Return to Text
13. We consider Japanese and non-Japanese intermediate imports separately in section 3. Return to Text
14. See Appendix C.6 for more discussion, as well as results obtained using other specifications. Importantly, in a reduced sample abstracting from zeros, a weighted regression using percentage changes directly yields estimates that are very close to those presented here. Return to Text
15. We have explored relaxing this assumption to allow our measure of K to vary across {$$ \tau$$ ,$$ \tau-1$$ }. The results are quantitatively similar. Return to Text
16. After constructing $$ \mu_i$$ according to equation (11) we average by industry to reduce the level of noise. Return to Text
17. This problem is discussed in greater detail in, for example, Ackerberg et al. (2006). Return to Text
18. Of course, the size and exogeneity of the shock also helps with this concern: any idiosyncratic productivity movements during this time are surely subsumed by the earthquake/tsunami. Return to Text
19. An unobserved input that could remain operative in our case is that of factor utilization. Since the scope for substantial adjustment along this dimension seems quite limited, we remain confident that our estimates would be robust to the inclusion of this missing ingredient. Return to Text
20. We observe little systematic variation in prices (see appendix C.4), and so this approach seems more appropriate in this setting. Return to Text
21. This comes from assuming a real interest rate of 4 percent, an annual depreciation rate of 10 percent, and then adjusting for a semi-annual frequency. The estimates are insensitive to alternative values of the rental rate. Return to Text
22. In principle it is possible to construct a firm-specific value for $$ \alpha$$ , using value-added information available in a census year. We are currently exploring the feasibility of this option. Return to Text
23. Specifically: $$ L_{\tau-1} = \frac{1}{6}Emp_{2010Q3} + \frac{1}{2}Emp_{2010Q4} + \frac{1}{3}Emp_{2011Q1} $$ and $$ L_{\tau} = \frac{1}{2}Emp_{2011Q2} + \frac{1}{2}Emp_{2011Q3}$$ . Return to Text
24. Those transactions with missing or imputed quantity information are dropped. Future efforts will evaluate whether it is possible to recover the quantity values from prior transaction details. Return to Text
25. Specifically, we drop any firm that has more than 3 months of zeros for any of these values, in period $$ \tau-1$$ or $$ \tau$$ . Return to Text
26. Since the appropriate measure of size in our context is output, we follow our convention and use the relative amounts of North American exports in the period before the shock as the weights. Return to Text
27. Of course, shocks can be passed through to affiliates through other means as well. See Peek & Rosengren (1997) and Peek & Rosengren (2000) for the case of U.S. affiliates from Japan. Return to Text
28. See Appendix Figure A2, which replicates Figures 5 (Panel A) and 6, but only for non-Japanese imports. Return to Text
29. An interesting result from Kurz & Senses (2013) is that firms that only import are actually more volatile than the domestic-only benchmark. Return to Text
30. Krishna & Levchenko (2014) outline theoretical results showing that for a given elasticity value (in their case, Leontief), volatility in output per worker should be actually decreasing in the number of inputs used. An interesting extension would be to look at results for total output. Return to Text
31. Ramondo et al. (2015) is another example arguing for a more nuanced framework for MP. Return to Text
32. The model of knowledge sharing in Keller & Yeaple (2013) is one attempt to analyze the dynamics between such transfers being accomplished in embodied (intra-firm trade) or disembodied (direct communication) form. Alternatively, domestic content requirements may provide incentives to produce inputs in one location over another. Return to Text
33. Appendix C.10 details further features of this data and explains how the sample was constructed. Return to Text
34. These firms are Honda, Mitsubishi, Nissan, Toyota, and Subaru. Return to Text
35. The average monthly plant-level production at these firms during December 2010 through February 2011 was about 12,200 units a month. The magnitude of the drop in May was -7200 units. Return to Text
36. We describe additional results on the behavior of inventories, sales, and production in Japan in Appendix C.3. Return to Text
37. In addition, one might be concerned that the N.A. exports series may be contaminated with Japanese imports whose country of ultimate destination is Canada/Mexico (a.k.a "in-transit shipments" - imports to Canada/Mexico via U.S.). These shipments should not be picked up in the reporting systems underlying the LFTTD. According to section 30.2(d)(1) of the U.S. Code of Federal Regulations, "In-transit shipments of goods from one foreign country to another where such goods do not enter the consumption channels of the United States are excluded from filing the Electronic Export Information (EEI)." Additionally, the Army Corps of Engineers has suspended the requirement to file the Form 7513, Shippers Export Declaration (SED) for In-transit Goods leaving the United States via vessel. Return to Text
38. Recent work finds that material inputs from foreign countries are imperfectly substitutable with domestic inputs for Hungary (Halpern et al. (2015)) and India (Goldberg et al. (2010)) Return to Text
39. The point estimate for the elasticity of substitution among intermediate inputs from Atalay (2014) is 0.03. Return to Text
40. The term bigram refers to two consecutive characters within a string (the word bigram contains 5 possible bigrams: "bi", "ig", "gr", "ra", and "am"). The program is a modified version of an existing string comparator algorithm by Michael Blasnik, and assigns a score for each variable between the two datasets based on the percentage of matching bigrams. See Flaaen (2013) or Wasi & Flaaen (2014) for more information. Return to Text
41. The primary sources of such ancillary evidence are clerical review of the matches, and additional parent identifier matching evidence. Return to Text
42. The multinational status of firms from the UBP directories are more straightforward. Return to Text
43. Some of these cases also apply to the UBP-BR bridge. Return to Text
44. This is one advantage of the survey data on multinational firms available from the Bureau of Economic Analysis. There are, however, a number of critical disadvantages of this data source, as outlined in Flaaen (2014). Return to Text
45. To be more precise, this set will include a combination of intermediate and capital goods. Return to Text
46. Another option is to use the CM-Materials file, the flip side of the CM-Products file. Unfortunately, the CM-Materials file contains significantly more problematic product codes than the Products file, and so concording to the trade data is considerably more difficult. Return to Text
47. It is worth pointing out that the most obvious way that this would materialize is by vertical integration of the firm in its U.S. operations. Provided that the industry designation of the firm pertains to its most downstream operations, then this is would not serve to bias the firms' classification of imported goods, as the upstream products are not actually "final" goods for that firm. Return to Text
48. The existence of final good inventories, on the other hand, makes a distinction between the production and sales of a particular product. Here, the presence of final good inventories implies that the firm can continue to sell from existing inventory stocks even while production is temporarily affected. Return to Text
49. Unfortunately, the CM data does not report imported materials inventory separately. Return to Text
50. These numbers are broadly similar, though somewhat lower than other estimates in the literature. See Ramey (1989) for one example. Return to Text
51. At first glance, the average monthly supply of these output inventories looks surprisingly low. On the other hand, it is probably the case that inventories are held jointly by the manufacturer and wholesale/retail establishments Thus, considering the inventories of manufacturers alone could potentially under-represent the "true" level of output inventories available for smoothing out production disturbances. Return to Text
52. This point has been made in somewhat differing contexts, by Kremer (1993) and Jones (2011). Return to Text
53. Because the source of this downward pressure on the estimate for $$ \psi $$ (or $$ \omega$$ ) is itself a very low product-level elasticity, it is unclear whether this should be considered a bias in the traditional sense. Return to Text
54. The level of aggregation we use attempts to balance concerns along two dimensions. With less product aggregation (i.e. HS10 level), one might be concerned with the inherent lumpiness of product-level firm imports. More product aggregation, on the other hand, could mask important product differences within a particular product grouping. Return to Text
55. Specifically, in equation (5) the $$ \gamma_{p}$$ 's would be higher than would be expected without the shock, and hence the $$ \beta_{p}$$ 's artificially low. Return to Text
56. The BR itself receives quarterly payroll and employment information for business and organizational employers from the IRS: Form 941, the Employer's Quarterly Federal Tax Return. For more information on the BR (formerly the SSEL), see Walker (1997). Return to Text
57. We re-weight the control group as described in section 2.1. Return to Text
58. Specifically, the $$ MMI_{ik}=1$$ if the average Japanese manufacturing establishment corresponding to a U.S. firm is above the median (roughly an MMI of 5.2) of all firms with Japanese manufacturing locations. The affected ports are: Onahama, Hitachi, Kashima, Haramachi, Shiogama, Sendai, Shimizu, Ishinomaki, Hashinohe, Miya Ko, Kamaishi, Ofunato, and Kessennuma. Return to Text
59. The combined effect of the coefficients on Japan and JPN*Exp is -0.16, and not significant. Return to Text

This version is optimized for use by screen readers. Descriptions for all mathematical expressions are provided in LaTex format. A printable pdf version is available. Return to Text