..... The examples of sunk costs include the costs for the initial marketing
International Trade Dynamics with Sunk Costs and Productivity Shocks⇤ Daisuke Fujii† The University of Chicago Job Market Paper‡ November 25, 2013
Abstract This paper offers a unified framework to analyze both short- and long-run trade dynamics in a consistent manner. It explains “the international elasticity puzzle”, a low trade elasticity against temporary shocks and a high trade elasticity against permanent shocks, studied in Ruhl (2008). The model in this paper extends the idea of export sunk costs and uncertainty to more general productivity shock processes by embedding the classical theory of “export hysteresis” into a continuous-time trade model with heterogeneous firms, and considers the effects of both aggregate and idiosyncratic productivity shocks. A sharp analytical characterization of the equilibrium elucidates the microfoundations of trade dynamics linking a static trade model with heterogeneous firms and an international macroeconomic model. Due to the sunk costs and uncertainty, firms do not change their export status against small temporary shocks. Aggregate productivity shocks and export sunk costs explain the elasticity puzzle because of the different adjustments on the extensive margin. If the productivity shocks are idiosyncratic, the economy is in a steady state, with individual firms moving around within a stationary distribution of productivities. Export hysteresis gives rise to a region of firm productivity where both exporters and non-exporters coexist given the same current productivity. The full model incorporates both types of shocks and offers realistic microfoundations of trade dynamics including simultaneous export entry and exit, an evolving productivity density of exporters, and the sluggish trade response to aggregate shocks.
⇤
I am grateful to Samuel Kortum and Nancy Stokey for their advice and constant encouragements. I would also like to thank Fernando Alvarez, Costas Arkolakis, David Burk, Lorenzo Caliendo, Thomas Chaney, Robert Lucas, Robert Shimer, Serginio Sylvain and the seminar participants of International Trade Lunch at the University of Chicago and Yale University, and Capital Theory Workshop at the University of Chicago for helpful comments. All remaining errors are my own. † Email:
[email protected] ‡ Latest version available at the author’s website: https://sites.google.com/site/fujii0622/research
1
1
Introduction
This paper develops a unified theoretical framework that captures both “export hysteresis,” the reluctance of firms to enter or exit from export markets in response to cyclical trade shocks, and the “international elasticity puzzle,” a low trade elasticity against temporary shocks and a high trade elasticity against long-lived shocks. It incorporates export fixed costs and sunk costs, as well as productivity shocks, in a continuous-time framework with heterogeneous firms. Both aggregate and idiosyncratic productivity shocks are allowed. To develop intuition, two specialized versions of the model, each with only one type of shock, are considered. Then the full model is studied. A sharp analytical characterization of the equilibrium elucidates the microfoundations of trade dynamics linking a static trade model with heterogeneous firms and an international macroeconomic model. The elasticity of exports with respect to the terms of trade is estimated to be low in international real business cycle (IRBC) models, typically ranging between 0.5 and 2. On the other hand, in static general equilibrium models of trade that compare steady states, the elasticity with respect to trade libralizations is estimated to be quite high, typically ranging between 5 and 10. Ruhl (2008) recognizes this discrepancy, and coins the term “the international elasticity puzzle”. He reconciles the contradictory findings by building a model in which business cycles are caused by temporary productivity shocks, and exporting entails an up-front sunk cost. Due to the sunk cost, few firms enter the export market in response to temporary shocks while a larger number of firms start exporting when there is a tariff reduction, which is considered to be permanent. The different adjustments on the extensive margin explain the elasticity puzzle. This paper extends that model of export sunk costs and uncertainty to more general productivity shock processes and incorporates rich microfoundations for trade dynamics. The key economic mechanism is export hysteresis caused by sunk costs and uncertainty. Many empirical studies (Roberts and Tybout (1997), Bernard and Jensen (2004), and Das, Roberts and Tybout (2007)) find the existence of substantial up-front sunk costs for exporting and a high persistence of export participation. In an environment where an action of discrete choice requires a fixed cost and the future profit is uncertain, “hysteresis” behavior is prevalent. Hysteresis is defined as the failure of an effect to revert itself when its underlying cause is removed. The export hysteresis was first investigated by Baldwin (1988) and many others. Dixit (1989) analyzes the optimal entry and exit decisions of a firm when the output price follows a random walk and the action requires a sunk cost. A firm waits to enter the export market until the output price rises sufficiently high, higher than the amortized sunk cost. After entry, the firm will not exit if the output price goes back to the original level. The firm continues to export until the output price drops to a sufficiently low level, lower than its continuation fixed cost, implying that the firm might run a temporary negative profit. The gap between the entry and exit triggers is called “the band of inaction”, which causes the export hysteresis behavior. The band of inaction is studied extensively in the literature of investment under uncertainty pioneered by Dixit and Pindyck (1994). The model in this paper 2
embeds the classical theory of export hysteresis into a continuous-time trade model to examine the effects of both aggregate and idiosyncratic productivity shocks. By employing the continuous time setup and techniques developed in stochastic optimal control literature, the model generates sharp characterizations of equilibrium variables, and the export dynamics is described as a solution to a partial differential equation. Hence, the comparative statics or dynamics can be seen in a transparent way. The main contribution of this paper is to present a framework of trade dynamics, which is consistent with micro and macro empirical regularities at both short- and long-run time horizons. The aggregate productivity shocks and export sunk costs are the essential ingredients to reconcile the international elasticity puzzle. Introducing idiosyncratic productivity shocks that follow a Brownian motion leads to the endogenous Pareto distribution of cross-sectional productivity, and gives rise to a range of firm productivity where both exporters and non-exporters coexist given the same current productivity. The full model captures the effects of both types of shocks and offers realistic microfoundations of trade dynamics such as simultaneous export entry and exit, an evolving productivity density of exporters, and the sluggish trade response to aggregate shocks. To build intuition, the model is presented step by step. After laying out the basic setup and explaining the mechanism of export hysteresis, I first focus on the case with aggregate productivity shocks only. The idiosyncratic firm productivity is fixed at a Pareto distribution. The aggregate shock hits all firms proportionally. At any date, exporters and non-exporters are separated by a cutoff firm, which evolves stochastically. If the export profit of the cutoff firm is in the inaction region, there is no firm entry nor exit, and the cutoff productivity does not change. When the export profit of the cutoff firm hits an upper bound, less productive firms enter the export market keeping the cutoff profit constant. A similar result holds when the cutoff profit hits the lower bound. Due to the inaction region, the extensive margin (set of exporters) of trade does not respond to a small temporary shock, so the trade elasticity only reflects the intensive margin in the short run. When a series of favorable shocks hits a country, less productive firms start entering and the extensive margin becomes operative. This results in a large trade elasticity in the long run. Therefore, the aggregate shock and sunk costs alone explain the international elasticity puzzle. Ruhl (2008) reviews the empirical evidence on the impact of the extensive margin adjustments in the long run in response to trade liberalizations. Kehoe and Ruhl (2013) also report that the extensive margin accounts for an important part of trade growth following tariff reductions whereas it does not change in the absence of trade policy or structural change. After characterizing the equilibrium properties of the model with aggregate shocks only, I turn to the case where productivity shocks are idiosyncratic. A firm-level productivity evolves as a geometric Brownian motion which gives rise to the double-Pareto distribution of the cross-sectional productivity as in Luttmer (2007). The idiosyncratic shocks make the export profit stochastic, and hence, we observe the export hysteresis at the firm level. Nonetheless, the aggregate variables and exporter distribution are constant in the stationary equilibrium. A range of firm productivity 3
in which both exporters and non-exporters coexist given the same productivity level emerges. The model explains the smooth increase in export participation as firm size becomes larger. The exporters in the “mixed region” are the ones whose productivity hit the entry margin, but dropped to a lower level. The non-exporters are those whose productivity has not reached the entry trigger yet since the last time they exited from the export market. This version of the model is essentially the same as the one in Impullitti et al. (2013), but simpler because I abstract from free entry and sunk costs for the domestic market. This allows a sharp characterization of the stationary equilibrium and clear comparative statics. The full model considers the effects of both types of shocks. Both aggregate and idiosyncratic shocks follow a Brownian motion. Firms derive the optimal strategies for export entry and exit, which characterizes the inaction region. The distribution of relative firm productivity is Pareto, but the exporter distribution evolves stochastically. The exporter dynamics can be characterized as the solution to a partial differential equation, which enables us to analyze the dynamics in response to any aggregate shock paths without simulations. The full model exhibits simultaneous export entry and exit, “the mixed region” of exporters and non-exporters, and the sluggish transition dynamics of exporter’s distribution in response to aggregate shocks. It offers an alternative explanation to the J-curve relationship between the terms of trade and exports as in Backus et al. (1994). Rather than the time to build capital, the time to “build” the set of exporters explains the J-curve. Besides Ruhl (2008), the papers by Alessandria and Choi (2007, 2012) are closely related to this research. They develop an open-economy dynamic stochastic general equilibrium (DSGE) model with heterogeneous firms, export sunk costs and both idiosyncratic and aggregate productivity shocks to examine the aggregate effects of trade liberalizations. Though it can incorporate more general features, the model is constructed in discrete time, and does not provide closed-form expressions for the firm distribution nor optimal entry/exit strategies. Ghironi and Melitz (2005) also construct a two-country DSGE model with heterogeneous firms to study the exporter dynamics and persistent deviations from PPP, but the model does not have sunk costs, and therefore, there is no hysteresis. This paper aims to elucidate the mechanism of the exporter dynamics in more transparent way using a continuos-time model. In addition to Dixit (1989) and Dixit and Pindyck (1994), many analytical techniques are borrowed from Stokey (2009). The self-contained book covers the necessary mathematical tools for stochastic control models with fixed costs, and contains rich economic applications. As she suggests, those modeling tools and concepts may have a wide applicability in many economic fields. The model developed in this paper can be applied to other situations in which heterogeneous agents face sunk costs and uncertainty. The paper proceeds as follows. The next section lays out the basic structure of the economy with heterogeneous firms. Section 3 describes a firm’s dynamic problem of entry and exit, and analyzes the properties of the inaction region for a general diffusion process. Sections 4 characterizes the export dynamics with aggregate shocks only, which explains the international elasticity puzzle. Section 5 characterizes the stationary equilibrium in which only idiosyncratic shocks hit firms. Section 6 4
presents the full model which incorporates both aggregate and idiosyncratic productivity shocks. Section 7 calibrates the full model to match the U.S. data and performs numerical analyses. Section 8 concludes.
2
Setup
Time is continuous and indexed by t. Consumers have “love for variety”, and their preferences are modeled by a CES utility function as in many other trade models. A continuum of firms with measure M produces the differentiated varieties. There is a measure L of consumers, and each consumer supplies one unit of labor inelastically. Labor is the only factor of production. There are two countries, home and foreign. Foreign variables are denoted by an asterisk. To understand the different effects of aggregate and idiosyncratic productivity shocks, I first focus on the case with only aggregate shocks up through Section 4. Aggregate productivity shocks hit all firms proportionally. Due to the export sunk costs, firms solve not only a static profit maximization problem but a dynamic entry and exit problem for exporting.
2.1
Demand
Within a country, all consumers are identical and have the same homothetic risk-neutral preference. In what follows, I will focus on variables in the home country. Foreign variables are described analogously with an asterisk. At any date t, the representative consumer maximizes the following expected discounted utility U t = Et
ˆ
1
t
e
⇢(s t) (1 ⇠) ⇠ qs Cs ds
,
where ⇢ is the preference discounting rate, and qt is the consumption of freely traded homogeneous good. The homogeneous good serves as a numeraire with its price normalized to one. The composite good Ct is produced by combining differentiated goods according to the following CES aggregator Ct =
ˆ
ct (!)
✓ 1 ✓
d!
(✓✓1)
,
⌦t
where ct (!) is the consumption of a variety !. The elasticity of substitution across varieties is denoted by ✓ and common for both countries . The set of available goods ⌦t will be characterized in a recursive general equilibrium evolving over time. Total expenditure in the home country (home GDP) is denoted by Yt . The demand for the homogenous good is qt = (1
5
⇠) Yt
and the demand for a differentiated variety ! at time t is ct (!) =
✓
pt (!) Pt
◆
✓
Ct =
✓
pt (!) Pt
◆
✓
⇠Yt , Pt
where the consumption-based price index is defined as Pt =
ˆ
pt (!)
1 ✓
d!
1 1 ✓
⌦t
There is no borrowing and saving decisions by households because all consumers are identical and trade is assumed to be balanced. Thus, total expenditure must equal the sum of labor income and total profits of the home firms at any time t, giving the following budget constraint wL + ⇧t = Yt where ⇧t is the total profits (both from domestic and foreign markets) of the home firms in the differentiated good sector.
2.2
Trade Costs
There are three types of costs pertaining to export. The standard “iceberg” transportation cost of exporting goods from home to foreign is denoted by ⌧ . If one unit of good is shipped, only a fraction 1 ⌧
units arrives. In addition to this variable cost, a firm must pay a fixed cost f regardless of the
quantity it exports. This fixed cost must be paid at any given time for a firm to remain as an exporter. There is no fixed cost for the domestic sales, so firms never exit from the domestic market. The third type of trade costs is a sunk cost for export entry denoted by . When a firm starts exporting, it must pay this sunk cost, which is assumed to be larger than the fixed cost: f < . This is a one-time sunk cost so that the firm does not have to pay it again if it continues to export. However, if a firm exits from exporting and wishes to resume exporting at a later date, it must pay again. The examples of sunk costs include the costs for the initial marketing research and establishing distribution networks in the foreign market. All trade costs are the same across firms within a country. The fixed and sunk export costs f and are measured as a portion of the destination’s GDP. Thus, in the unit of numeraire, the fixed and sunk export costs are f Yt⇤ and Yt⇤ respectively. This is similar to the usage of foreign labor for trade costs if there is no freely traded homogeneous good and wage rates are determined in a general equilibrium. The proportionality of those costs to the foreign GDP simplifies the firms’ export entry and exit problem. Also, they are assumed to be iceberg costs, meaning that the incurred costs don’t go to labor income. This assumption simplifies the equilibrium conditions, but does not change the nature of the derived results.
6
The three types of trade costs have distinct effects on trade. The per-unit cost ⌧ affects the intensive margin, that is, how much each exporter exports. A permanent reduction in ⌧ also induces more firms to enter the export market affecting the extensive margin (number of exporters) in the long run. The fixed cost f is responsible for the separation of exporters and non-exporters. Since f is assumed to be the same for all firms, only efficient firms who are productive enough to recoup the fixed cost can export. The sunk cost controls the sluggishness of trade adjustment on the extensive margin in a dynamic horizon. If is large, the set of exporters is very persistent. The three different roles played by each trade cost will become clear in Section 4.
2.3
Production
The homogenous and differentiated goods are produced under constant returns to scale. One unit of labor in the home country produces w units of the homogeneous good. Since the homogeneous good is freely traded and its price is normalized to be 1, the wage in the home country is w if this good is produced. I only consider the equilibrium in which both countries produce some of the numeraire. This allows the countries to differ both in size L and in productivity w as in Chaney (2008). There is a continuum of differentiated good producers in both countries. Firms are differentiated by their idiosyncratic productivity
which is constant over time (the stochastic evolution of
will
be incorporated in Section 5 and 6). Also, a country is hit by an aggregate productivity shock Zt at any given time. This aggregate shock hits only the differentiated good sector. Hence, the economy-wide productivity is given by w whereas Zt is the relative productivity of the differentiated good sector, which follows a diffusion process specified in the next section. Firms have the following linear production function y (Zt ; ) = Zt l, where y is the output and l is the labor used. Each firm maximizes its expected discounted sum of profits. The firm’s optimization problem has two components: a static monopolist’s profit maximization problem in each period and a dynamic discrete choice problem of export entry and exit. Let pD,t ( ) and pX,t ( ) denote the domestic and export prices respectively of a home firm with productivity . Hereafter, I will use subscripts D and X to denote domestic and export variables. The prices charged by firm
are ✓
◆
w ✓ 1 Zt ✓ ◆ ✓ ⌧w pX,t ( ) = ⌧ pD,t ( ) = ✓ 1 Zt
pD,t ( ) =
✓
7
The variable profits (gross of fixed export cost) are ✓ ◆ ✓ ◆ ⇠ pD,t ( ) (1 ✓) Z t Pt ✓ 1 Yt = a1 Yt ✓ Pt w 8 ✓ ◆(1 ✓) ✓ ◆ > Zt Pt⇤ ✓ 1 ⇤ ⇤ < ⇠ pX,t ( ) Yt = a1 Yt if firm exports ✓ Pt⇤ ⌧w ⇡X,t ( ) = > : 0 otherwise
⇡D,t ( ) =
where a1 =
⇠ ✓
⇣
✓
✓ 1
⌘1
✓
. The net export profit is ⇡X,t ( )
f Yt⇤
= a1
⇣
Zt Pt⇤ ⌧w
⌘✓
1
(1)
(2)
f Yt⇤ if firm
exports. Due to the sunk export cost and productivity uncertainty, the decision of export entry and exit will be optimally determined by the firm’s dynamic discrete choice problem.
2.4
Productivity Distribution
For both countries, the relative productivity distribution is assumed to be Pareto with a shape parameter
and a lower bound
min
> 0, so its cumulative distribution and density functions are ✓
G( ) = 1
min ( +1)
◆
g( ) = min
for
2[
min , 1).
Assumption 1.
The following assumption ensures that aggregate variables are finite. >✓
1
At any given time, there is a cutoff firm which separates exporters and non-exporters. Denote the cutoff by ¯t . Any firm with > ¯t exports at time t and the rest sell only in the domestic market. ¯t is a key variable in the recursive equilibrium and evolves over time depending on the history of aggregate shocks. Define the following average productivities ˜= ˜X,t =
ˆ "
1
1
✓ 1
dG ( )
1 ✓ 1
=⌫
min
1 G ¯t
ˆ
1 ¯t
✓ 1
dG ( )
min
#
1 ✓ 1
= ⌫ ¯t
⇣ ⌘ 1 ✓ 1 where ⌫ = . When productivity follows a Pareto distribution, the average productivity (✓ 1) ˜ is proportional to the lower bound. This result reflects the scale free property of the Pareto distribution, and simplifies the expressions of aggregate variables. The measure of firms in the home country is fixed at M . There is no firm entry or exit for the 8
domestic market, so M is constant over time.1 Then the measure of exporting firms is ⇥
G ¯t
MX,t = M 1
2.5
⇤
✓ ¯ ◆ t
=M
min
Aggregate Variables and Equilibrium
In terms of aggregate implications, this model is isomorphic to one where M firms of productivity ˜ produce only for the domestic market and MX,t firms of productivity ˜X,t export to the foreign market as demonstrated in Melitz (2003). Let p˜D,t and p˜X,t be the average domestic and foreign prices of the home country. Using the previous results, we have p˜D,t p˜X,t
⇣ ⌘ ✓ ✓ ◆ w = pD,t ˜ = ✓ 1 Zt ⌫ min ⇣ ⌘ ✓ ✓ ◆ w⌧ = pX,t ˜X,t = ✓ 1 Zt ⌫ ¯t
The home and foreign price indices can be written as2 "
¯t ⇤ ✓ 1
Pt = a 2 M ⇤ "
Pt⇤ = a2 M ¯t
where a2 = Notice that
⇣
✓
✓ 1 Pt⇤ is
⌘
✓ 1
✓
✓
w⇤ ⌧ ⇤ Zt⇤
w⌧ Zt
◆1
◆1 ✓
✓
+M(
+ M⇤ (
min )
min )
✓ 1
✓ 1
✓
✓ w⇤ Zt⇤
w Zt
◆1
◆1
✓
✓
#
#
1 1 ✓
(3)
1 1 ✓
(4)
1 ✓ min
. The home export profits are determined by the price relative to Pt⇤ . decreasing in both Zt and Z ⇤ , and increasing in ¯t . A favorable technological ⌫
t
shock in either country reduces the price level and increases welfare. Conditional on Zt and Zt⇤ , high ¯t implies that the price of imported goods (from foreign’s perspective) is higher, so the price index will be high as well. ⇥
In the current model, aggregate productivity shocks Zt = [Zt , Zt⇤ ] and exporter cutoffs ¯ t = ⇤ ¯t , ¯⇤ are state variables. Unlike the Chaney-Melitz type trade models, the cutoff productivity t 1
Alternatively, we can incorporate free entry with fixed cost and exogenous firm death shock to allow firm entry. In that case, any aggregate productivity shocks are absorbed by the adjustment in M , and the average total (domestic + foreign) profit is constant over time. To explicitly illustrate the different adjustments of intensive and extensive margins against different kinds of shocks, I shut down the adjustment on M . 2 Using the expressions for the average prices, the home price index can be written as Pt =
h
M p˜1D,t✓
+
⇤ MX,t
1 ✓ p˜⇤X,t
i
1
1
✓
"
= M
✓✓
✓ ✓
1
◆
w Zt ⌫ min
9
◆1
✓
+M
⇤
✓ ¯⇤ ◆ t min
✓✓
✓ ✓
1
◆
w⇤ ⌧ ⇤ ⇤ Zt ⌫ min
◆1
✓
#
1
1
✓
¯t cannot be pinned down by the structural parameters and destination’s characteristics. Due to the sunk export cost and export hysteresis behavior, we need a complete history of the shocks Z to obtain ¯ t . Given Zt and ¯ t , the price indices are calculated according to (3) and (4). The average net profit in the home country is MX,t [˜ ⇡X,t f Yt⇤ ] M ⇣ ⌘ = ⇡X,t ˜X,t are the average profits of all firms and exporters ⇡ ˜t = ⇡ ˜D,t +
⇣ ⌘ where ⇡ ˜D,t = ⇡D,t ˜ and ⇡ ˜X,t
respectively. The total home expenditure is
Yt = wL + M ⇡ ˜t In this model, ⇡ ˜t acts like wage in standard general equilibrium models. The average profit ⇡ ˜t affects both income and expenditure of a country and the vector [˜ ⇡t , ⇡ ˜t⇤ ] is the fixed point of the two excess demand functions. The derivations of the equilibrium average profits ⇡ ˜t and ⇡ ˜t⇤ are relegated to Appendix A.1. Trade is not necessarily balanced in the differentiated good sector, but I assume any imbalances in the differentiated good sector is made up for in the homogeneous good sector so that aggregate trade is balanced.
3
Export Hysteresis
3.1
Firms’s Dynamic Problem
Since there is no fixed or sunk entry costs for the domestic market, all firms sell in the domestic market. For the foreign market, they solve a dynamic optimal entry and exit problem due to the export sunk cost and profit uncertainty. Let ⇡ ˆX,t = be written as
⇡ ˆX,t
✓
⇡X,t Yt⇤
Zt⇤ ¯ , t, Zt
be the export profit normalized by the foreign GDP. Using (2) and (4), it can
◆
= a3 ⇥
✓ 1
| {z } fixed, idiosyncratic
0
B ⇥@ |
M
⇣
¯t
⌘✓
1 1
⇣
Zt⇤ Zt
w + M⇤ w⇤ ⌧ min {z stochastic, aggregate
⌘✓
1 C
1A
(5)
}
where a3 is a constant.3 If we take the log of (5), it is clear that the stochastic movement of the export profit is the same for all firms. Only the level is different due to the idiosyncratic productivity heterogeneity. Notice that the normalized export profit depends on the ratio of the productivity shocks 3
a3 =
Zt⇤ Zt .
If aggregate productivities rise by the same magnitude in both countries leaving
a1 ✓ a1 ( min )✓ 1 2
=
⇠ ✓
⇣
(✓ 1)
⌘
(1 ✓) min
10
Zt⇤ Zt
constant, export profit is not affected. In that situation, the export price pX,t will decrease due to the home productivity shock, but the foreign price index will also decrease by the same proportion due to the foreign shock, which offsets the extra profit. Also, the home exporter cutoff ¯t is the aggregate state variable which evolves over time depending on the history of productivity shocks as shown below. Since > ✓ 1, when ¯t increases, ⇡X,t ( ) increases. For exporters, an increase in ¯t implies a higher foreign price index implying milder competition there. In the dynamic entry and exit problem, firms take the stochastic evolution of ⇡ ˆX as exogenous, and compare the expected value of being an exporter to the sunk costs to determine the optimal entry and exit thresholds. Conditional on exporting, each firm follows the monopolist’s profit maximization rule described in Section 2.3. Given that result, they solve a dynamic stochastic control problem of export entry and exit. Recall that the fixed and sunk export costs f Yt⇤ and Yt⇤ are proportional to the foreign GDP. Normalizing the export profit and the associated costs by Yt⇤ , we can analyze the stochastic evolution of ⇡ ˆX,t in relation to f and to derive the optimal entry and exit policies. The proportionality of export costs to Yt⇤ prove useful in this analysis since any effect of Yt⇤ on the export profit is offset by p
( )
the same change in export costs, and hence, we can focus on the effect of the relative price X,t on Pt⇤ Zt⇤ ¯ trade. Because ⇡ ˆX,t in (5) only depends on , and , it is convenient to define the log of relative productivity as
Zt
zt = ln
✓
Zt Zt⇤
◆
Then, the state variables of ⇡ ˆX,t zt , ¯t , are zt and ¯t . In this section, it is assumed that the effect of ¯t on export profit is negligible so that the only stochastic variable of ⇡ ˆX,t (zt , ) is zt . This assumption corresponds to the case of a small home country discussed in Section 4. The general case where the export profit is a function of both zt and ¯t is discussed in Appendix A.3. Firms maximize the expected discounted sum of export profits net of trade costs by optimally choosing the timings of entry and exit. Due to the aggregate shock, the export profit follows a diffusion process. This is a classical problem of investment under uncertainty (e.g. Dixit (1989) and Dixit and Pindyck (1994)). An exporter and a non-exporter can be viewed as assets that are call options on each other, and the optimal strategy is characterized by a pair of threshold productivity levels. The export sunk cost and uncertainty coming from aggregate productivity shocks give rise to export hysteresis studied by Baldwin (1988) and Baldwin and Krugman (1989).
3.2
Hamilton-Jacobi-Bellman (HJB) Equations
Hereafter, I will use ⇡ (z) to denote the normalized export profit ⇡ ˆX,t (z, ) to simplify notations. The subscript
is dropped since the structure of the problem is identical for all firms. Yet, it should
be understood that all value functions and derived triggers depend on
11
. Consider the following
general diffusion process for z dz = µ (z) dt + (z) dW where W is a standard Brownian motion (Wiener process) and µ (z) and
(z) are the infinitesimal
drift and variance parameters respectively. They are time invariant, but can depend on the current status of z. Let V0 (z) to be the value function for a non-exporter when the current aggregate shock is z, and similarly V1 (z) for an exporter. By Ito’s lemma, we obtain 1 dV0 = V00 (z) dz + V000 (z) 2 Therefore,
0
E (dV0 ) = µ (z) V0 (z) +
1 2
2
2
(z) dt
00
(z) V0 (z) dt
The asset equilibrium condition gives the following Hamilton-Jacobi-Bellman (HJB) equations for V0 and V1 4 ⇢V0 (z) = µ (z) V00 (z) +
1 2
2
(z) V000 (z)
1 2 ⇢V1 (z) = ⇡ (z) f + µ (z) V10 (z) + (z) V100 (z) | {z } | 2 {z } dividend capital gain
(6) (7)
The interpretation of the above HJB equations is straightforward. The return on the asset (LHS) must equal the sum of dividend and the capital gain (RHS). For non-exporters, there is no dividend profit flow, but V0 is still positive due to the option value of entering the export market in the future. For exporters, the dividend is the net export profit ⇡ (z)
f and the capital gain is the option value
of exiting from export market. Both equations are second-order linear ODE, but (6) is homogeneous whereas (7) is nonhomogeneous. The ODE’s (6) and (7) characterize the general solutions for V0 and V1 . To pin down the particular solutions, we need boundary conditions. For a second-order linear ODE, two boundary conditions are required to characterize the particular solution. Since we have two of those ODE’s, we need four boundary conditions in total. The optimal solution to this stochastic dynamic control problem is characterized by two threshold levels of z: zL and zH with zL < zH . zH is the threshold at which a non-exporter becomes an exporter, and zL is the threshold at which an exporter becomes a non-exporter. Since zL < zH , this is an inaction region. If the current aggregate shock z is in this range [zL , zH ], the optimal policy is to remain the status quo; exporters should continue exporting and non-exporters should not start exporting. Thus, the current aggregate productivity z is not a sufficient statistic to determine if a particular firm is an exporter or not if z falls in the inaction band for the firm. The definition of zL 4
See Chapter 3 of Stokey (2009) for the detailed treatment of HJB equations. Basically, HJB equations can be derived by considering a small change t, applying E (dV ) and taking the limit t ! 0.
12
and zH suggests the following value-matching conditions5 (8)
V0 (zL ) = V1 (zL ) V0 (zH ) = V1 (zH )
(9)
These boundary conditions pin down the form of V0 and V1 . For optimality, we need the following smooth pasting conditions6 V00 (zL ) = V10 (zL )
(10)
V00 (zH )
(11)
=
V10 (zH )
Above four conditions (8) - (11) characterize the particular solutions for V0 and V1 , and the optimal thresholds zL and zH . Because the two HJB equations (6) and (7) share the same homogeneous part, it is convenient to define the following function as in Dixit (1989) U (z) = V1 (z)
V0 (z)
Then, we can collapse the two HJB equations into f + µ (z) U 0 (z) +
1 2
(z) U 00 (z)
(12)
U (zL ) = 0, U (zH ) = , U 0 (zL ) = 0, U 0 (zH ) = 0
(13)
⇢U (z) = ⇡ (z)
2
The value-matching and smooth pasting conditions become
Any two of the above boundary conditions determine the particular solution for U (z) and the other two conditions are required to solve zL and zH . Therefore, for any general diffusion process of z with µ (z) and
(z), above expressions (12) and (13) determine zL and zH . An example of V0 , V1
and U is illustrated in Figure 1. Because V0 (z) only consists of the option value of entering the export market, it must have limz!0 V0 (z) = 0 since when z is very small, the probability of rising to zH in any given finite time is almost zero. Also V1 (z) consists of the profit flow which is increasing 5
The non-exporter value function can be expressed as ⇢ 0 1 V0 (z) = max V1 (z) , V0 (z) + µ (z) V0 (z) + 2
2
(z) V
00
(z)
⇢V0 (z) dt
for the region of z where being a non-exporter is optimal. At zH , the firm is indifferent between being an exporter and non-exporter implying V0 (zH ) = V1 (zH ) . A similar argument holds for the case of zL . 6 See Chapter 4 of Dixit and Pindyck (1994) and Stokey (2009) for more detailed treatments of this condition. If the smooth pasting conditions are not satisfied, two value functions are connected at a kink. It can be shown that the optimal switching strategies are no longer valid due to the local convexity around the kink.
13
in z and the option value of exiting which is decreasing in z. V1 (z) exhibits a U-shape: when z is small, the option value dominates, and when z is large, the profit flow dominates. When z follows a Brownian motion, closed-form solutions for V0 (z) and V1 (z) are available, but in general, one must solve the ODEs numerically.
3.3
Band of Inaction
Because zL < zH , the band of inaction [zL , zH ] causes the export hysteresis behavior. We can convert the band of inaction of z to that of the export profit, [⇡L , ⇡H ]. Consider a firm whose export profit ⇡X is in the inaction region [⇡L , ⇡H ]. The firm starts exporting when the export profit ⇡ goes above ⇡H but does not stop exporting when ⇡ comes back to the original level. When ⇡ decreases to low enough, which is ⇡L , the firm will exit. Dixit (1989) analytically characterizes the band of inaction, and performs comparative statics when ⇡ follows a geometric Brownian motion. He demonstrates that ⇡L < f and f + ⇢ < ⇡H . The inaction region is illustrated in Figure 2. The standard Marshallian trigger profit levels are f and f + ⇢. When there is no uncertainty, a firm should start exporting when ⇡X rises to f + ⇢, which is the amortized full cost of exporting. Similarly, an exporter should exit when its export profit drops to f . The presence of uncertainty expands the Marshallian inaction band [f, f + ⇢]. Thus, it is possible that an exporter receives negative profit when ⇡ falls in [f, ⇡L ] since the firm hopes ⇡ to rise in the near future. Similarly, the profit level of f + ⇢ is not enough to attract a non-exporter to become an exporter. Due to the uncertainty, ⇡ must rise more than f + ⇢ to attract a non-exporter. As
2 ⇡
! 0, ⇡H ! f + ⇢
and ⇡L ! f . Also, as becomes smaller, the inaction region [⇡L , ⇡H ] shrinks but doesn’t vanish. In essence, sunk cost creates the inaction band and the uncertainty parameter
2
amplifies it.
Export hysteresis is the key ingredient to derive the results in the following sections.
4
Model with Aggregate Shocks
This section analyzes the trade dynamics in which the source of uncertainty only comes from the aggregate productivity shocks. I focus on the case of small home country to elucidate the workings of export hysteresis on the international elasticity puzzle. The analysis of the general case is presented in Appendix A.3 because the full model also uses the small home country assumption. The following assumption is placed in this section. Assumption 2. The home mass of firms is small: M ⇡ 0. With this assumption, the share of home firms in the foreign price index is negligible so that ¯t does not affect Pt⇤ nor ⇡X,t . The normalized export profit (5) can be expressed as ⇡ ˆX,t (z, ) = S ( ) e(✓
14
1)z
,
(14)
value V0 HzL
30
V1 HzL
25
k
20 15 10 5
zL 0.2
zH 0.4
0.6
0.8
1.0
z
(a) V0 (z) and V1 (z)
value 6
5
4
U@zD
3
2
1
zL 0.2
zH 0.4
0.6
0.8
(b) U (z)
Figure 1: An example of V0 (z) , V1 (z) and U (z)
15
1.0
z
dont' export
inaction region
π_L
export
f
f+ρκ
π_H
negative profits π
Figure 2: Entry and exit decisions where S ( ) = a4
✓ 1 .7
The scaling constant S ( ) does not change over time. Because the state
variable of ⇡ ˆX,t is z alone, we can apply the results of the previous section to derive the inaction region. Another ingredient for the problem is the diffusion process of z, which is described below.
4.1
Band of Inaction: Brownian Motion
Suppose z follows a Brownian motion with the following law of motion dz = µdt + dW where W is a Wiener process (standard Brownian motion). The infinitesimal drift and variance parameters µ (z) and
(z) don’t depend on the current level of z. The HJB equations for exporters
and non-exporters become ⇢V0 (z) = µV00 (z) + ⇢V1 (z) = Se(✓
1 2
2
V000 (z)
f + µV10 (z) +
1)z
(15) 1 2
2
V100 (z)
(16)
In this case, we obtain closed-form solutions for V0 and V1 . Lemma 1. The value functions that solve the HJB equations (15) and (16) take the following form V0 (z) = Ae↵z V1 (z) = Be 7
a4 =
a3 M⇤
⇣
w⇤ ⌧w
⌘✓
1
(
min )
(✓ 1)
16
(17) z
+
⇡ (z) ⌘
f ⇢
(18)
where ↵ = =
µ+
p p
µ
µ2 + 2⇢ 2
µ2 + 2⇢
2
⌘ = ⇢
2
⇡ (z) = Se(✓
2
(✓
2
2
1)2
µ (✓
1)
1)z
and A and B are constants determined by associated boundary conditions.
8
Proof. See Appendix A.2. The non-exporter’s value function V0 (z) is valid in the interval ( 1, zH ) and the exporter’s
value function V1 (z) is valid in (zL , 1). The U (z) function becomes U (z) = Be
z
Ae↵z +
⇡ (z) ⌘
f , ⇢
and the associated boundary conditions provide the following system of four equations which characterizes A, B, zL and zH . 2
F1 (A, B, zL , zH )
6 6 F2 (A, B, zL , zH ) F =6 6 F (A, B, z , z ) L H 4 3 F4 (A, B, zL , zH )
3
2
Ae↵zL
Be
zL
7 6 7 6 Ae↵zH Be zH 7=6 7 6 ↵Ae↵zL + Be zL 5 4 ↵Ae↵zH + Be zH
Se(✓
1)zL
⌘
Se(✓
(✓ (✓
f ⇢
3
2
0
3
7 6 7 6 7 + 7 7=6 0 7 (✓ 1)z 6 7 1) Se ⌘ L 7 5 4 0 5 (✓ 1)z 0 1) Se ⌘ H
1)zH
⌘
+
+
f ⇢
(19)
We do not obtain the explicit analytical expressions for zL and zH , but we can do comparative statics analytically using the implicit function theorem. Figure 1 plots the value functions when z follows a Brownian motion with the parameters summarized in Table 1. The exit and entry triggers are zL = 0.311 and zH = 0.98, that imply ⇡L = 0.638 and ⇡H = 1.74. We confirm that 0.638 = ⇡L < f = 1 and 1.35 = f + ⇢ < ⇡H = 1.74. Thus, it is possible that an exporter receives negative profit when ⇡ falls in [f, ⇡L ] hoping ⇡ to rise in the near future. Similarly, the profit level 8 The second part of equation (18) has a nice interpretation. Since z follows a Brownian motion with µ and 2 , the aggregate productivity ZZ⇤ = ez follows a geometric Brownian motion with a drift µ + 12 2 and variance 2 . The additional 12 2 term comes from the convexity of exponential function combined with Ito’s lemma. Then, ⇡ (z) follows a geometric Brownian motion with a drift µ⇡ = (✓ 1) µ + 12 (✓ 1)2 2 and variance ⇡2 = (✓ 1)2 2 . So the following relationship holds ˆ 1 ⇡ f ⇡ f = =E e ⇢t (⇡t f ) dt ⌘ ⇢ ⇢ µ⇡ ⇢ 0
Thus, the second part of (18) is the expected present value of exporting if the firm exports forever given the initial export profit ⇡. This is explained in Dixit (1989). The first part of (18) Be s corresponds to the option value of exiting from the export market optimally. Similarly, V0 (s) = Ae↵s is the option value of entering the export market.
17
value
µ 0.01
2
0.01
✓
1 1.5
⇢ 0.07
S 0.4
f 1
5
Table 1: Baseline parameters for the Brownian motion of f + ⇢ is not enough to attract a non-exporter to become an exporter. Due to the uncertainty, ⇡ must rise more than f + ⇢ to attract the non-exporter. Again, it should be noted that each firm
derives different inaction regions of the aggregate
shock [zL, , zH, ] because the scaling constant S = a4
✓ 1
is different. Let ⇡L, = ⇡X, (zL, ) and
⇡H, = ⇡X, (zH, ) be the lower and upper threshold levels of export profit for firm . With the Brownian motion assumption, we have the following lemma. Lemma 2. When the aggregate shock z follows a Brownian motion, the entry and exit trigger profits ⇡L, and ⇡H, are homogeneous of degree one in f and jointly. When z follows a Brownian motion, the export profit follows a geometric Brownian motion. The above result can be found in Dixit (1989) and Dixit and Pindyck (1994). Since ⇡X , f and are measured in the same unit, any rescaling by f and does not change the nature of the problem when ⇡X follows a geometric Brownian motion.
4.2
Evolution of the Cutoff Firm
Each firm
solves the above problem to determine the optimal threshold levels of the aggregate
shock: zL, and zH, . The following lemma can be proved. Lemma 3. When z follows a Brownian motion, the threshold levels of export profit ⇡L = ⇡X, (zL, ) and ⇡H = ⇡X, (zH, ) do not depend on . Even though each firm derives different threshold productivity levels, the threshold export profit levels are the same for all firms. This is evident by considering the parallel movements of ln ⇡ ˆX, of different ’s as pointed out in 3.1. We can also see the result by considering the movement of the export profit. When z follows a Brownian motion with µ and geometric Brownian motion with a drift µ⇡ = (✓
1) µ+ 12
(✓
1)
2
2, 2
the export profit ⇡ follows a
and variance
2 ⇡
= (✓
1)2
2.
Then, the HJB equation (12) can be written in terms of ⇡ ⇢U (⇡) = ⇡
f + µ⇡ ⇡U 0 (⇡) +
1 2
2 2 00 ⇡ ⇡ U (⇡)
All firms face the same problem and derive the common thresholds ⇡L and ⇡H . Because firm’s idiosyncratic productivity
is constant over time and the aggregate shock affects the export profit of all firms proportionally, there is a cutoff firm ¯t which separates exporters and non-exporters at any given time. Firms with > ¯t export and firms with < ¯t are non-exporters. 18
Denote the export profit of the cutoff firm by ⇡ ¯t = ⇡X ¯t . The aggregate shock influences all firm’s export profit including that of the cutoff firm ⇡ ¯t . Whenever the extensive margin of trade (set of exporters) changes, the entry or exit must start from the cutoff firm ¯t . Figure 3 illustrates the evolution of ⇡ ¯ and ¯ in this environment. The aggregate productivity shock z (the top line plot of Figure 3) starts rising at t1 which raises ⇡ ¯ (the middle plot) as well. At t2 , ⇡ ¯ hits the upper ¯ threshold ⇡H . At this point, firms below enters the export market, which implies a decrease in ¯ (the bottom plot). When ⇡ ¯ is at ⇡H , any further positive shocks will be absorbed by the adjustment in ¯ keeping ⇡ ¯ constant at ⇡H . At t3 , z starts falling, and hence, ⇡ ¯ falls as well. Yet, this decline in ⇡ ¯ does not affect ¯ because ⇡ ¯ is still above ⇡L which is the exit trigger. So the current export profit is not low enough to induce the threshold firm to exit. The aggregate shock z bumps up between t4 and t5 but this does not change ¯ either because ⇡ ¯ does not hit the upper threshold. From t5 , z and ⇡ ¯ start falling again, and ⇡ ¯ hits the lower threshold ⇡L at t6 . Exporters start exiting from firm ¯ first. The negative shock continues until t7 , so ¯ rises between t6 and t7 making ⇡ ¯ constant at ⇡L . At t7 , the shock turns to positive, and ¯ stops rising.
4.3
Short Run vs. Long Run
The interaction between ⇡ ¯ and ¯ along with the inaction band [⇡L , ⇡H ] generates the main results in this section. The inaction band emerges from the sunk export cost and uncertainty. As long as the cutoff profit ⇡ ¯ remains within the inaction band, the set of exporters does not change. When the cutoff profit is at one of the thresholds, any further shocks are completely absorbed by the updating of the cutoff firm ¯. The current model generates a low trade elasticity in the short run and a high elasticity in the long run as summarized in the next proposition. Proposition 1. (Short run vs. Long run) Let ✏xz be the trade elasticity with respect to the aggregate shock
Z Z⇤ .
With Assumptions 1 and 2, ✏xz =
(
✓
1 if ⇡ ¯ 2 (⇡L , ⇡H ) (short run)
if ⇡ ¯ = ⇡L or ⇡H (long run)
From Assumption 1, the long-run trade elasticity
is larger than the short-run elasticity ✓
1.
I will use ✏XY to denote the elasticity of X with respect to Y throughout the paper. Denote the export sales of firm
and the aggregate export sales (normalized by Yt⇤ ) by xt ( ) and Xt
respectively. Using the results in Section 2, it is informative to express the aggregate export as Xt =
⇣ ⌘ MX,t ⇥ xt ˜X,t = a4 ⇥ ¯t | {z } | {z } set of exporters average export
19
⇥ e(✓
1)z ¯✓ 1 t
t1
t2
t3
t4
t5
t6
Time !!!!!!
Figure 3: Dynamics of ln ⇡ ¯ and ln ¯
20
t7
Figure 4: The cutoff exporter ¯ and export profit ⇡ ¯ !
where a4 is a! constant. Taking the log of the above expression yields ln Xt = ln a4 + (✓
1) zt
(
(✓
1)) ln ¯t
The trade elasticity with respect to z is ✏xz = (✓
1) + (✓
1
) ✏ ¯z
where ✏ ¯z is the elasticity of ¯ with respect to z. When ⇡ ¯ 2 (⇡L , ⇡H ), ✏ ¯z = 0 and when ⇡ ¯ is on ✓ 1 Z one of the boundaries, ¯ is being updated to keep ⇡ ¯ / ¯ Z⇤ constant which implies ✏ ¯z = 1 confirming Proposition 1. Because the elasticities of ⇡ ¯ and ¯ with respect to z are ✓ 1 and 1 respectively, when there is a linear movement in z, ⇡ ¯ and ¯ also have linear movements with slopes multiplied by ✓
1 and
1 as shown in Figure 3. The log of the cutoff profit ln ⇡ ¯ follows a regulated
Brownian motion with the two-sided regulators ln ⇡L and ln ⇡H . The long-run stationary probability density of x = ln ⇡ ¯ is (x) = where
=
2(✓ 1)µ (✓ 1)2 2
e
xH
e
xL
e
x
and xi = ln ⇡i for i = L, H (see Chapter 3 of Dixit and Pindyck (1994) for
derivation). If µ = 0, the long-run average of ln ⇡ ¯ is
ln ⇡L +ln ⇡H . 2
In most of the time, the cutoff profit ⇡ ¯ is in the inaction band, so only the intensive margin (average export sales of the current exporters) adjust to the movement of z (short run). When 21
either positive or negative shocks keep hitting the economy, the extensive margin (change in ¯) will be operative, which generates a larger trade elasticity. When = 0, the inaction band vanishes (⇡L = ⇡H ) and the result of Chaney (2008) is recovered. Because there is no sunk cost, the extensive margin always adjusts to a change in z generating the trade elasticity .
5
Model with Idiosyncratic Shocks
This section presents a trade model in which there is no aggregate shock, but idiosyncratic productivity evolves as a Brownian motion. The growth process of firm productivities is borrowed from Gabaix (1999) and Luttmer (2007). I abstract from the balanced growth path analysis, but incorporate export sunk costs in the model. Due to the sunk costs and idiosyncratic productivity uncertainty, each firm exhibits the export hysteresis. Yet, a stationary equilibrium emerges in which the productivity distribution of all firms and exporters are both modeled by a double Pareto distribution as in Luttmer (2007). The same idea is investigated in Impullitti et al. (2013), but their analysis is restricted for the case of two symmetric countries. The current model can be extended to asymmetric multi-country case.
5.1
Productivity Shocks and Stationary Distribution
The basic setup is the same as in Section 2, but the Pareto distribution of productivity endogenously emerges from the random growth of idiosyncratic productivities. The log of idiosyncratic productivity
is assumed to evolve as a Brownian motion. Define s = ln . The law of motion of s
is (20)
ds = µdt + dW
where W is a standard Brownian motion. The entry point of productivity is fixed at s⇤ . When new firms enter the economy, they all start from productivity s⇤ . After that, each idea follows a diffusion process as specified in (20). There is an exogenous death shock rate
which hits all ideas uniformly.
Let a be the age of a firm (time elapsed after birth) and f (a, s) be the density of s of age a. The following Kolmogorov Forward Equation (henceforth KFE) describes the evolution of productivity density f (a, s)9 @f (a, s) = @a
f (a, s)
µ
@f (a, s) 1 + @s 2
2@
2f
(a, s) @s2
(21)
The first term on the RHS reflects the reduction of density due to the exogenous death shock. The second term reflects the outflow of the density due to the drift of s. The third term reflects the inflow of density due to randomness. The solution to the above second-order PDE is the following 9
See Chapter 3 of Dixit and Pindyck (1994) for more details. The KFE is “the balance equation”, which links the densities of evolving particles at different points in time.
22
discounted normal density with mean s⇤ + µa and variance a a
f (a, s) = e
p
1
e
2
2⇡a
2
[s (s⇤ +µa)]2
(22)
2a 2
The derivation of the above solution can be found in Harrison (1985) and Luttmer (2007). The intuition for the above density is clear. All ideas starts from s⇤ and follow a Brownian motion. Since the drift and variance are proportional to the elapsed time a, the distribution of productivities of age cohort a has a mean s⇤ + µa and variance a
2.
Due to the death shock , the density of age
cohort a must be discounted by a. This implies that old ideas have a small weight since most of their cohort have died out up until time a. The first term e
a
reflects this discounting. Conditional
on surviving, they have larger variance. Let f (s) be the stationary probability density of s. This is the marginal density of f (a, s) with ´1 respect to s: f (s) = 0 f (a, s) da. By the definition of stationarity, we have @f@a(s) = 0. Then, we obtain the following KFE which describes f (s) 0=
f (s)
µ
@f (s) 1 + @s 2
2@
2f
(s) @s2
(23)
This is a homogeneous second-order linear ODE (not PDE). The general solution to this ODE takes the following form 1s
f (s) = C1 e where
1
and
2
2 2
µ
Quadratic formula gives 1
2 1, 2
(24)
are the two solutions to the following characteristic equation 1 2
so that
2s
+ C2 e
> 0. Since
> 0,
1
>
2µ 2
= =
µ+
=0
p µ2 + 2
2
2
p µ2 + 2
µ
and
(25)
2
2
2
> 0. The constants C1 and C2 are determined by
additional conditions. The stationary probability density f (s) must satisfy the following conditions as well (this is from Arkolakis (2011)) lim f (s) = 0
s! 1
ˆ µ [f (s⇤ )
f (s⇤ +)] +
s⇤
f (s) ds + 1
1 2
2
⇥
ˆ
1
f (s) ds = 1
s⇤
f 0 (s⇤ ) 23
0 8s 2 ( 1, 1)
f (s)
f 0 (s⇤ +)
⇤
=
(26) (27) (28) (29)
The first condition is the lower boundary condition. For the shock process of s,
1 is the absorbing
barrier. The second and third conditions ensure that f (s) is a probability density. The last condition
emerges by integrating the KFE (23). This states that the net inflows to the density are equal to the net outflows. The particular solution of (24) which satisfies the above four conditions is given by f (s) =
8 > < > :
1 2 1+ 1 2 1
+
e
1 (s
s⇤ )
for s < s⇤
2 (s
s⇤ )
for s
2
e 2
(30)
s⇤
where derivations can be found in Reed (2001) and Arkolakis (2011).10 This is a skewed (asymmetric) Laplace distribution with a mode at s⇤. For the convergence of the stationary distribution, it is assumed that
> 0. The resulting stationary distribution of
g( ) =
where
⇤
⇤
8 > > > < > > > :
1 2 1+
2(
1 2 1
+
2
(
1
1
⇤)
1
2
1
⇤)
2
for
= es is the following11
all firms with
H ],
0 and s > sL
where 1 (a, s|sH ) = p
a
"
N
✓
s
sH p
µa a
◆
µ(sH sL ) 2 /2
e
N
✓
s + sH
2s p L a
µa
◆#
and N (·) a standard normal density. More details about this solution can be found in Harrison (1985). As a tends to zero, the first term in the square brackets approaches to a point mass as sH and the second term becomes zero since s+sH > 2sL . Also, if sL tends to
1, the above expression
becomes equation (22) derived earlier. The size density of exporters f1 (s) is given by f1 (s) =
ˆ
1
1
0
where
1
and
2
min e(
1 2
f1 (a, s) da =
+
2
1 + 2 )(s
⇥
sL ) , e(
1 (sH
e
sL )
1 + 2 )(sH
⇤
2 (s
1 e
sL )
1
sL )
for s > sL
were computed in subsection 5.1. Since the entry barrier sH is a point mass, above
expression is also a probability density. Thus,
f1 (s) =
8 > > >
> > :
2
⇥
h ⇤ e 1
1 e
1 (sH
sL )
1 (s
sL )
2 (s
e
sL )
e( 1 + 2 )(sH sL ) 1 ⇥ ⇤ e e 1 (sH sL ) 1 1+ 2 1 2
i
for s 2 (sL , sH )
2 (s
sL )
for s > sH
This was introduced in Luttmer (2007) and used in Impullitti et al. (2012) as well. The corresponding probability density for
= es is
g1 ( ) =
where (
2
L
= esL and
H
8 > > > > > > > > < > > > > > > > > :
1 2 1
+
2
2⇣
⌘
6 L 4 ⇣ ⌘ +
1
H L
1 2 1
2
3
1+ 2
2⇣
H L
6 4 ⇣
17 5
1 ⌘
H L
(
1
= esH . All firms with
for
2
L
1+ 2
⌘
2 +1)
1 >
3 17 5 H
2(
L,
H)
(32) (
2 +1) 2
for
>
H
L
export and their Pareto tail index is
+ 1) which is equal to the tail index of the overall stationary distribution. This confirms the
scale free property of the Pareto distribution. Figure 5 depicts the overall stationary density of firms and exporter’s density. From (31), we can easily show that the share of firms with 1 1+ 2
⇤
is
0.6
0.5
0.4
0.3
0.2
Non-exporters 0.1
Exporters
*
fH
f
0.0
1
fL
2
3
4
5
Figure 5: Stationary Density of Firms As in Meliz (2003), define the following average productivity for exporters ˜=
ˆ
1
✓ 1
dG1 ( )
1 ✓ 1
L
We have the following result. Lemma 4. The average productivity for exporters ˜ is given by the following expression 2
where k1 =
˜=6 4
1
1) and k2 =
+ (✓
2
✓
1 2
k1 k2
(✓
◆
⇣
⇣
H L H L
3✓11
⌘ k1 ⌘
17 5 1
1
(33)
L
1).
Proof. See Appendix A.4. By using L’Hopital’s Rule, we can confirm that lim
L! H
˜=
h
2 2
(✓ 1)
i
1 ✓ 1
H
which is also
obtained in Ghironi and Melitz (2005). For further analyses, it is useful to define q = lemma 2, we can see that ⇡X ( implies that q =
H L
L)
and ⇡X (
H)
H L
. Using
are homogeneous of degree one in f and . This
is determined by the ratio
f.
The average productivity of exporters (33)
summarizes relevant information for all aggregate variables. Though not pursued in this paper, the results in this section can be extended to asymmetric multi-country general equilibrium of trade extending Chaney (2008)13 . 13
Let i and j be the subscripts for the source and destination countries. If we assume that
27
ij fij
is common for all
6
Full Model
This section presents a unified model which incorporates both aggregate and idiosyncratic productivity shocks. The basic setup is the same as before, but now there are two types of productivity shocks st (idiosyncratic) and zt (aggregate) that follow a Brownian motion. Assumptions in this section are summarized below. 1. (small home country): The home mass of firms relative to foreign is small: 2. (aggregate shock): The log of relative productivity z = ln
Z Z⇤
motion dz =
M M⇤
⇡ 0.
follows a drift-less Brownian
z dW
3. (idiosyncratic shock): The log of idiosyncratic productivity s = ln ˆ evolves as a Brownian motion ds = µs dt +
s dW
4. (independence): Idiosyncratic shocks are i.i.d, and independent from aggregate shocks. 5. (parameters):
>✓
1,
1)2
> (✓
2 s
2
Again, the small country assumption simplifies the diffusion process of the export profit, and is a good approximation for most of the countries. Because there is no fixed or sunk entry costs for domestic sales, firms never exit from the domestic market unless they are hit by the death shock. The aggregate shock does not affect the stationary distribution of the relative firm productivities. The stationary density of all firms is given by equation (31). The entry rate of new firms is fixed at M to keep the mass of operating firms constant. New firms are born at productivity s⇤ as nonexporters. The fixed and sunk export costs are proportional to the foreign GDP, and the diffusion process of the normalized export profit ⇡ ˆX is analyzed as in Section 4. Unlike Section 5, the density of the cross-sectional productivities of exporters is not constant due to the aggregate fluctuations. The dynamics of the exporter’s density depends on the optimal entry and exit strategy of firms and the aggregate shock path, which is analyzed below. country pairs but allow variation in fij , then ⇣ ⌘qij is equal but fij varies across countries. Using (33), we can define the average price for exporters as p˜ij = p ˜ij and the expression of the price index of country j becomes Pj =
✓ ✓
1
N X
Lkj
⌫kj (wk ⌧kj )
1 ✓
k=1
where ⌫kj =
Mk
k2,k
2,k
✓
+ 2,k k1,k k2,k
1,k
Other comparative statics can be conducted as well.
28
◆
k
qkj1,k 1,k + 2,k
qkj
!11✓
1 1
(34)
6.1
Band of Inaction
By the small country assumption, a firm’s export profit can be written as ⇡ ˆX (zt , st ) = Se(✓
1)(zt +st )
where S is a scaling constant. The aggregate and idiosyncratic productivity shocks enter the export profit function symmetrically. Because of these two shocks, the export profit also follows a diffusion process. The associated HJB equation is a second-order PDE (not ODE): rU (z, s) = ⇡ (z, s)
f
µs U s +
1 2
2 s Uss
+
1 2
2 z Uzz
(35)
where r = ⇢ + . In the zs plane, two curves separate the domain into three regions of different optimal strategies: export, inaction, and do not export. We can also express the threshold levels of s as functions of z: sL (z) and sH (z). The curves sL (z) and sH (z) characterize the lower and upper thresholds of the export hysteresis. The following boundary conditions must be satisfied U (z, sL (z)) = 0, U (z, sH (z)) = , Uz (z, sL (z)) = Us (z, sL (z)) = Uz (z, sH (z)) = Us (z, sH (z)) = 0 (36) This is a free boundary problem of a partial differential
equation14 .
In general, the problem is quite
difficult to analyze because not only the surface function, but the boundary conditions are also unknown. Nevertheless, when both shocks follow a Brownian motion, we can reduce the problem to an ODE with an analytical solution. Proposition 2. Given the assumptions in this section, the value function U (z, s) has the following form U (z, s) = Be
(z+s)
Ae↵(z+s) +
⇡ (z + s) ⌘
f , r
where where p
µ2s + 2r ( z2 + s2 ) ( 2 + s2 ) p z µs µ2s + 2r ( z2 + s2 ) = ( z2 + s2 ) 2 2 z + s ⌘ = r (✓ 1)2 µs (✓ 2
↵ =
µs +
1)
and A and B are constants determined by associated boundary conditions. 14 An example of this type of problem is called Stefan problem in physics, which analyzes the heat evolution of ice passing to water.
29
s 1.2 1.0 10
UHz,sL
0.8
Export
5 0.2
0 -5
0.6 0.4
0.0
0.0
0.5
s
0.0
-0.2
z
Inaction region
0.2
-0.2
1.0
(a) Plot of U (z, s)
Do not export 0.0
0.2
0.4
0.6
0.8
1.0
1.2
z
(b) Contour lines of U (z, s) = 0 and U (z, s) =
Figure 6: Value function U (z, s) and its contour sets Proof. See Appendix A.5. The result is intuitive. A firm’s export profit depends on the combined productivity v = z + s, which follows a Brownian motion with drift µv = µs and variance
2 v
=
2 z
+
2 s
since z and s are
independent. Then, we can solve the lower and upper thresholds vL and vH as in subsection 4.1, that characterize the two curves sL (z) = vL
z and sH (z) = vH
z. Along these lines, the boundary
conditions are satisfied. Note that the value function U (z, s) are symmetric about the line z = s. The two curves that separate the regions of optimal strategies are both linear with a slope of -1. An important result is that the distance between sH (z) and sL (z) are constant at vH
vL for any
z. This result is used below to transform the KFE. The graph of U (z, s) and the contour lines at U (z, s) = 0 and U (z, s) = are shown in Figure 6.
6.2
Dynamics of the Exporter’s Density
The above result gives the optimal entry and exit strategies for firms given the aggregate shock z. Consider a continuous path of z (t), which determines the paths of sL (z (t)) and sH (z (t)). At any date t, all firms with s > sH (z (t)) are exporters whereas firms with s < sL (z (t)) are non-exporters. Write l (t) = sL (z (t)) and h (t) = sH (z (t)). We can characterize the dynamics of the exporters’ density in the inaction range [l (t) , h (t)] as follows. Let m (s, t) be the density of exporters s in the inaction region [l (t) , h (t)] at time t. Because the productivities follow a Brownian motion, m (s, t) must satisfy the following KFE as in Section 5 @m (s, t) = @t
m (s, t)
µs
@m (s, t) 1 + @s 2
2@ s
2 m (s, t)
@s2
(37)
Let m0 (s) be the initial distribution of exporters in the inaction range [l (0) , h (0)] . The boundary 30
conditions are m (s, 0) = m0 (s) m (l (t) , t) = 0 1 2
m (h (t) , t) = 1
+
e
2 (h(t)
s⇤ )
2
The first condition specifies the initial distribution. The second condition reflects that l (t) is an absorbing barrier for exporters at time t. The third condition smoothly matches the density of exporters at h (t), where the expression is derived in (30). I only consider the aggregate shock paths by which h (t) does not go below s⇤15 . This is a moving boundary problem for a p.d.e. The domain changes over time, which makes the problem difficult to analyze16 . In our case, the problem can be reduced to a simpler one using the fact that the length of the inaction range is h (t) is constant. Let d = h (t)
l (t)
l (t). Consider the following transformation of the density function;
n (s, t) = m (s + l (t) , t). The domain of the function n (s, t) is fixed in s 2 [0, d]. The new function n (s, t) must satisfy the following p.d.e @n (s, t) = @t
n (s, t)
z 0 (t)
@n (s, t) @s
µs
@n (s, t) 1 + @s 2
2@ s
2 n (s, t)
@s2
(38)
The boundary conditions become n (s, 0) = n0 (s) n (0, t) = 0 n (d, t) = b (t) where n0 (s) = m0 (s
l (0)) and b (t) =
1 2 1+ 2
e
2 (h(t)
s⇤ ) .
We transformed the moving boundary
problem to the one with fixed domains. The new function n (s, t) can be numerically solved in Mathematica. Note the new term z 0 (t) @n(s,t) @s . The growth rate of the aggregate shock amplifies the effect of µs in the KFE because the entry and exit decisions are based on the export profit ⇡ (z + s) in which z and s are perfect substitutes. When the aggregate shock path is linear in time, we obtain the analytical expression of n (s, t) (See Appendix A.6). 15 If this happens, new firms are all exporters. The calculations of the exporter dynamics and its stationary density are complicated. 16 In physics, this type of problem appears in the analysis of the heat equation with moving insulating barriers. Various numerical, and some analytical solutions can be found in the literature.
31
Description
Parameter
Value
preference discount rate death shock rate elasticity of substitution variable trade cost sunk export cost relative to the fixed cost drift parameter (idiosyncratic) standard deviation (idiosyncratic) standard deviation (aggregate)
⇢
0.05 0.09 2 0.3 20.5 -0.0163 0.081 0.036
✓ ⌧ f
µs s z
Table 2: Parameters for calibration
7
Calibration and Numerical Analysis
7.1
Parameters
This section calibrates the full model to match the relevant U.S. data on the firm size distribution and exporter dynamics. With the calibrated model, I perform several comparative dynamics. The preference discounting rate ⇢ is set to 0.05, which implies the annual interest rate of 5 percent. For the benchmark case, ⌧ = 1.3 to match the results in Obstfeld and Rogoff (2000). I choose 2 for the elasticity of substitution between varieties ✓ to reflect the low trade elasticity in the short run. Following Impullitti et al. (2012), I set the death shock rate
= 0.09. This is based on
the average annual death rate for manufacturing firms in the 1998-2004 period from the 2004 U.S. Census data. The aggregate trade elasticity which augments the extensive margin adjustment is governed by the upper Pareto tail index which is
2.
I target the estimate from Eaton and Kortum (2002),
= 8.28. I borrow the drift p parameter of the idiosyncratic shocks µs = 1.63% from µs + µ2s +2 s2 Arkolakis (2011). Since 2 = , this gives s = 8.1%. As in Ruhl (2008), z is set 2 2
s
to match the volatility of the logged output in the U.S. for the period of 1950-2000, which gives z
= 3.6%. The fixed and sunk costs for exporting are calculated to match two moments of the U.S. exporter
distribution. I normalize s⇤ = 0 and the scaling constant for the export profit to be S = 1. Although the model does not have a stationary equilibrium, I set z (0) = 0 and interpret this as a benchmark. The moments of the U.S. firm distribution are matched to their counterparts in the stationary equilibrium with the absence of the aggregate shock. Bernard and Jensen (2004b) find that firms enter and exit from the export market at different productivity levels using U.S. data for the period of 1983-1992. According to their statistics, the exit trigger productivity is two-thirds of the entry trigger productivity on average. The implied inaction band is sH sL = 0.41. Impullitti et al. (2012) calibrate the ratio of the sunk and fixed export costs as the current 17
model17 .
I target sH
sL = 1 and obtain
f
f
= 91, which leads to sH
sL = 1.64 in
= 20.5. Next, Bernard et al. (2003) report
From lemma 2, we know that the inaction band is determined by the ratio of the sunk and fixed export costs.
32
that 21 percent of U.S. manufacturing firms are exporters in 1992. Given sL and sH , the share of exporters is18
which gives sH
7.2
⇥ e 2 sH e m1 = e 1 (sH sL ) = 0.184 and sL = 0.816.
1 (sH
sL ) 2 (sH
e
⇤ 1
sL )
,
Trade Dynamics
Given an aggregate shock path and initial exporter’s distribution, the full model has the ability to generate nonlinear export dynamics. The dynamics of the exporter distribution is characterized by the p.d.e (38), which enables us to compute aggregate export for any t. Two scenarios are investigated below with the implications on the short- and long-run trade elasticities. Those aggregate shock paths are particular realizations of z, which all firms expect to follow a Brownian motion. For both scenarios, the economy starts from the stationary distribution, and the KFE (38) is numerically solved to characterize n (s, t), which is converted back to m (s, t). The corresponding density function of = es , m ˆ ( , t), is then computed to obtain ˜t (the average productivity of exporters), which summarizes the model’s implications for aggregate export. Log of export (normalized by Yt⇤ ) is ln Xt = ln
"ˆ
where a5 is a constant and ˜t =
1
x( )m ˆ ( , t) d
L (t)
h´
1
L (t)
trade adjustments are represented by (✓
✓ 1m ˆ(
#
, t) d
= a5 + (✓ i
1) zt and (✓
h i 1) ln ˜t + zt ,
1 ✓ 1
. The intensive and extensive margins of 1) ln ˜t respectively. The log of export in
the initial period ln X0 is normalized to be 0. 7.2.1
Regime Switching
Consider the following regime-switching scenario. The aggregate shock z is initially 0, jumps by 10% within a quarter, and stays at the new level for 8 years as shown in Figure 7a. In the initial quarter, the total export increases by 31% implying the trade elasticity of 3.1. This number is a composite of the intensive margin elasticity of ✓ of 2.1. Let { {
L,1 ,
and
H,1 }
H,1