Open Access

The Role of Energy in a Real Business Cycle Model with an Endogenous Capital Utilization Rate in a Government Sector: Lessons from Bulgaria (1999–2016)


Cite

Introduction and Motivation

The average labor productivity in Bulgaria in the period following the introduction of currency board (1997) is highly pro-cyclical.

This stylized fact is similar to the finding documented in Bils and Cho (1994) for the United States as well.

The classical explanation for this stylized fact, which is observed in many developed economies (Jorgenson and Griliches, 1967), is that the major inputs of production, labor, and capital are used more intensively during periods of expansion as compared to periods of recession. In order to quantitatively rationalize this phenomenon, and gain a deeper understanding of the transmission mechanism responsible for economic fluctuations, we introduce an endogenous utilization rate of physical capital stock into a relatively standard real business cycle model in a government sector in detail. We examine the quantitative importance of the variability in capital utilization and its relevance to generate plausible cyclical movements in aggregate variables. More specifically, we investigate whether allowing for cyclical capital utilization helps our augmented real business cycle model match the empirical business cycles in Bulgaria in the period after the introduction of the currency board arrangement. This period was chosen in our investigation due to the fact that the introduction of the hard exchange-rate peg achieved macroeconomic stability in Bulgaria, and thus the time series have good statistical properties.

The other novelty in this paper is the particular way capital utilization enters the model. We follow Finn (2000) by adopting the empirical regularity that capital utilization requires energy and argue for the importance of energy in the transmission of technological shocks.

Earlier studies, e.g. Kydland and Prescott (1988), Greenwood, Hercowitz, and Huffman (1988), and Finn(1995) have incorporated varying capital utilization in real-business-cycle models but without featuring the energy use channel. Still, they make use of the Keynes (1936) notion of ”user cost of capital” but in a neoclassical framework, where the changes in capital utilization affect the marginal efficiency of capital. We make use of that notion in the current paper as well.

In turn, there are two costs to the capital utilization decision that are at play in the current model: (i) a cost in terms of higher energy use and (ii) a cost in terms of a higher depreciation rate of physical capital stock. The first is a direct effect working through the production function, and following from fact that energy becomes a de facto factor of production through the link with utilization rate and capital stock. The second is an indirect channel, which is one of the novelties in this paper. This effect occurs due the presence of a depreciation cost of utilization and the linkage between it and the use of energy, which works indirectly through the accumulated stock of physical capital. We then use this artificial economy with endogenous capital utilization through energy use as a model in order to study the importance of energy price shocks on the main aggregate variables. In order to be able to draw plausible quantitative predictions, we calibrate the theoretical economy to approximate Bulgaria in the period 1999–2016. We find that a positive shock to energy prices is akin to a negative technological shock and propose an explanation for a technological shock.

The novelty, however, is that the transmission mechanism of energy price shocks stems from a relatively little explored relationship between energy usage and services provided by physical capital, described in this paper. Put differently, since energy enters the production function only because it is essential to the utilization of capital, the endogenous variations in utilization and energy use would be inter-related.

It comes as no surprise that unexpected changes in world energy prices are very important for an energy-intensive production in Bulgaria, a small open economy. Energy price hikes or drops can have important real effects on the Bulgarian economy as Bulgaria imports most of its energy inputs (oil and natural gas in particular) from the Russian Federation. Next, from the perspective of the Bulgarian economy, the price of the aggregate energy input is taken as given. Thus, the industry structure of the energy production is not of central importance and will be ignored in this paper. More specifically, we abstract away from the issue, as it is of limited relevance for the international transmission of how changes in the price of imported energy inputs affect the Bulgarian economy.

In another line of research, Rotemberg and Woodford (1996a) present a theory based on imperfect competition in the oil market to explain business cycle fluctuations in the US economy. Hamilton (1983, 1985, 1996) studies the effect of oil price on real output in the US.

Instead, what takes a central stage in this paper is the fact that energy prices directly affect the productivity and the profitability of all sectors in the economy, and thus aggregate output. Overall, the model with endogenous utilization rate through the use of energy performs better than earlier real business cycle models when data for Bulgaria is considered. In particular, consistent with the observed cyclical fluctuations in Bulgaria, total hours follow output movement. Nevertheless, as with the standard RBC model, the model with endogenous utilization rate of capital falls far short of generating wage variability as in data, and the wage rate in the model is very strongly pro-cyclical, while wages are acyclical in data.

The rest of the paper is organized as follows. Section 2 describes the model framework and describes the decentralized competitive equilibrium system. Section 3 discusses the calibration procedure, and Section 4 presents the steady-state model solution. Sections 5 proceeds with the out-of-steady-state dynamics of model variables, and compares the simulated second moments of theoretical variables against their empirical counterparts. Section 6 concludes the paper.

Model Description

There is a representative household that derives utility out of consumption and leisure. The time available to households can be spent in productive use or as leisure. In addition, the household chooses optimally the rate at which capital stock is being utilized. The government taxes consumption spending and levies a common tax on all income, in order to finance wasteful purchases of government consumption goods, and government transfers. On the production side, there is a representative firm, which hires labor and utilized capital to produce a homogenous final good, which could be used for consumption, investment, government purchases, or energy consumption. Depreciation rate is endogenous and is a function of the endogenous capital utilization rate, and depends on the energy use.

Household’s Problem

There is a representative household, which maximizes its expected utility function, as in Finn (2000):

maxE0t=0βt{[ctψ(1ht)1ψ]1σ1σ}$$max\,{{E}_{0}}\sum\nolimits_{t=0}^{\infty }{{{\beta }^{t}}\left\{ \frac{{{\left[ c_{t}^{\psi }{{\left( 1-{{h}_{t}} \right)}^{1-\psi }} \right]}^{1-\sigma }}}{1-\sigma } \right\}}$$

where E0 denotes household’s expectations at period 0, ct denotes household’s private (non-energy) consumption in period t, ht are hours worked in period t, 0 < β < 1 is the discount factor, 0 < ψ < 1 is the relative weight that the household attaches to consumption and σ > 0 is the curvature of the utility function

This utility function is equivalent to a specification with a separable term containing government consumption, e.g. Baxter and King (1993).Since in this paper we focus on the exogenous (observed) policies, and the household takes government spending as given, the presence of such a term is irrelevant. For the sake of brevity, we skip this term in the utility representation above.

.

The household starts with an initial stock of physical capital k0> 0, and has to decide how much to add to it in the form of new investment, as well as the rate at which the stock of physical capital is being utilized. As a result, in every period, physical capital depreciates at an endogenous rate, which depends on the level of utilitization rate ut chosen by the household, so that 0 < δ(ut) < 1. Following Taubman and Wilkinson (1970), Greenwood, Hercowitz, and Huffman (1988), and Finn (1995, 2000), the functional form for the endogenous depreciation rate is as follows:

δ(ut)=ω0utω1ω1(0,1),$$\delta \left( {{u}_{t}} \right)={{\omega }_{0}}\frac{u_{t}^{{{\omega }_{1}}}}{{{\omega }_{1}}}\in \left( 0,1 \right),$$

where ω0> 0, ω1> 1. This depreciation function is consistent with Keynes’s (1936) notion of the “user cost of capital,” which argues that higher utilization causes faster depreciation, at an increasing rate, because of faster ”wear and tear” on the aggregate physical capital stock. In addition, as in Finn (2000), we assume that capital utilization requires energy, et.

This channel is missing from earlier studies, such as Taubman and Wilkinson (1970), Greenwood, Hercowitz, and Huffman (1988), and is one of the novelties of this paper.

More specifically, it will be postulated that energy spending complements the service flow from physical capital as follows:

etkt=a(ut)=v0utv1v1,$$\frac{{{e}_{t}}}{{{k}_{t}}}=a\left( {{u}_{t}} \right)={{v}_{0}}\frac{u_{t}^{{{v}_{1}}}}{{{v}_{1}}},$$

where ν0> 0, ν1> 1. The technical relationship function, a(.), the same as those developed by Finn (1995), postulates that energy is essential to the utilization of capital, with increases in utilization requiring more energy usage per unit of capital, at an increasing rate.

This modelling choice could be traced back to Jorgen and Grilliches (1967), who find that capital and electricity are complements in production. In addition, after some algebra, one can show that ut=(etkt)1v1(v1v0)1v0,$${{u}_{t}}={{\left( \frac{{{e}_{t}}}{{{k}_{t}}} \right)}^{\frac{1}{{{v}_{1}}}}}{{\left( \frac{{{v}_{1}}}{{{v}_{0}}} \right)}^{\frac{1}{{{v}_{0}}}}},$$

The law of motion for physical capital is then

kt+1=it+(1δ(ut))kt,$${{k}_{t+1}}={{i}_{t}}+\left( 1-\delta \left( {{u}_{t}} \right) \right){{k}_{t}},$$

and the real interest rate is rt, hence the before-tax effective (utilized) physical capital income of the household

in period t equals rtutkt. In addition to capital income, the household can generate labor income. Hours supplied to the representative firm are rewarded at the hourly wage rate of wt, so pre-tax labor income equals wtht. Lastly, the household owns the firm in the economy and has a legal claim on all the firm’s profit, πt.

Next, the household’s problem can be now simplified to

max{ct,ut,ht,kt+1}t=0E0t=0βt{[ctψ(1ht)1ψ]1σ1σ},$$\underset{\left\{ {{c}_{t}},{{u}_{t}},{{h}_{t}},{{k}_{t+1}} \right\}_{t=0}^{\infty }}{\mathop{\max }}\,{{E}_{0}}\sum\nolimits_{t=0}^{\infty }{{{\beta }_{t}}\left\{ \frac{{{\left[ c_{t}^{\psi }{{\left( 1-{{h}_{t}} \right)}^{1-\psi }} \right]}^{1-\sigma }}}{1-\sigma } \right\}},$$

s.t.

(1+τc)ct+kt+1(1δ(ut))kt+ptet=(1τy)[ωtht+rtutkt]+gtt+πt,$$\begin{array}{*{35}{l}}\left( 1+{{\tau }^{c}} \right){{c}_{t}}+{{k}_{t+1}}-\left( 1-\delta \left( {{u}_{t}} \right) \right){{k}_{t}}+{{p}_{t}}{{e}_{t}}= \\\left( 1-{{\tau }^{y}} \right)\left[ {{\omega }_{t}}{{h}_{t}}+{{r}_{t}}{{u}_{t}}{{k}_{t}} \right]+g_{t}^{t}+{{\pi }_{{{t},}}} \\\end{array}$$

where τc is the tax on consumption, τy is the proportional income tax rate (0 < τc, τy < 1), levied on both labor and capital income, pt is the relative (to the aggregate consumption price index) energy price, et denotes energy use in period t, and gtt$g_{t}^{t}$denotes government transfers. The household takes the two tax rates cy}, government spending categories,{gc,gt}t=0,$\left\{ {{g}^{c}},{{g}^{t}} \right\}_{t={0}}^{\infty },$profit {πt}t=0,$\left\{ {{\pi }_{t}} \right\}_{t={0}}^{\infty },$the realized technology process {At}t=0,$\left\{ {{A}_{t}} \right\}_{t={0}}^{\infty },$prices {ptwtrt}t=0,$\left\{ {{p}_{{{t}'}}}{{w}_{{{t}'}}}{{r}_{t}} \right\}_{t={0}}^{\infty },$and chooses {ct,ht,ut,kt+1}t=0$\left\{ {{c}_{t}},{{h}_{t}},{{u}_{t}},{{k}_{t+1}} \right\}_{t=0}^{\infty }$to maximize its utility subject to the budget constraint.

Note that by choosing kt + 1 the household is implicitly setting investment it optimally. Similarly, by choosing the utilization rate, the household is determining the time-varying depreciation rate. Lastly, by choosing the level of physical capital and the rate of capital utilization, that determines optimally the level of energy use.

The constraint optimization problem generates the following optimality conditions:

ct:[ctψ(1ht)1ψ]σψctψ1=λt(1τc)$${{c}_{t}}:{{\left[ c_{t}^{\psi }{{\left( 1-{{h}_{t}} \right)}^{1-\psi }} \right]}^{-\sigma }}\psi c_{t}^{\psi -1}={{\lambda }_{t}}\left( 1-{{\tau }^{c}} \right)$$ht:[ctψ(1ht)1ψ]σ(1ψ)(1ht)ψ=λt(1τy)wt$${{h}_{t}}:{{\left[ c_{t}^{\psi }{{\left( 1-{{h}_{t}} \right)}^{1-\psi }} \right]}^{-\sigma }}\left( 1-\psi \right){{\left( 1-{{h}_{t}} \right)}^{-\psi }}={{\lambda }_{t}}\left( 1-{{\tau }^{y}} \right){{w}_{t}}$$ut:δ(ut)+pta(ut)=ω0utω11+ptv0utv11=(1τy)rt$$\begin{align}& {{u}_{t}}:{{{{\delta }'}}^{\left( {{u}_{t}} \right)}}+{{p}_{t}}{{{{a}'}}^{\left( {{u}_{t}} \right)}}={{\omega }_{0}}u_{t}^{{{\omega }_{1}}-1}+{{p}_{t}}{{v}_{0}}u_{t}^{{{v}_{1}}-1} \\ & \,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,=\left( 1-{{\tau }^{y}} \right){{r}_{t}} \\ \end{align}$$kt+1:λt=βEtλt+1[1+(1τy)rt+1ut+1δ(ut+1)pt+1a(ut+1)]$$\begin{align}& {{k}_{t+1}}:{{\lambda }_{t}}=\beta {{E}_{t}}{{\lambda }_{t+1}}\left[ 1+\left( 1-{{\tau }^{y}} \right){{r}_{t+1}}{{u}_{t+1}} \right. \\ & \,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\left. -\delta \left( {{u}_{t+1}} \right)-{{p}_{t+1}}a\left( {{u}_{t+1}} \right) \right] \\ \end{align}$$limtβtλtkt+1=0$$\underset{t\to \infty }{\mathop{\lim }}\,{{\beta }^{t}}{{\lambda }_{t}}{{k}_{t+1}}=0$$

where λt is the Lagrangean multiplier attached to household’s budget constraint in period t.

The interpretation of the first-order conditions above is as follows: the first one states that for each household, the marginal utility of consumption equals the marginal utility of wealth, corrected for the consumption tax rate. The second equation states that when choosing labor supply optimally, at the margin, each hour spent by the household working for the firm should balance the benefit from doing so in terms of additional income generated, and the cost measured in terms of lower utility of leisure. The third equation describes the optimal utilization rate, which requires that the change in the depreciation rate, or the marginal cost in terms of an increased depreciation rate resulting from utilizing capital at a higher rate, together with the marginal cost in terms of additional energy used in the capital utilization, equal the after tax return on utilized capital. In other words, the marginal benefit resulting from physical capital services should balance with the user cost of capital at the margin. The fourth equation is the so-called Euler condition, which describes how the household chooses to allocate physical capital over time. The last condition is called the “transversality condition” (TVC): it states that at the end of the horizon, the value of physical capital should be zero.

Firm Problem

There is a representative firm in the economy, which produces a homogeneous product. The price of output is normalized to unity. The production technology is Cobb-Douglas and uses both utilized (effective) physical capital, utkt, and labor hours, ht, to maximize static profit

t=At(utkt)αht1αrtutktwtht,$$\prod{_{t}={{A}_{t}}{{\left( {{u}_{t}}{{k}_{t}} \right)}^{\alpha }}h_{t}^{1-\alpha }-{{r}_{t}}{{u}_{t}}{{k}_{t}}-{{w}_{t}}{{h}_{t}},}$$

where At denotes the level of technology in period t. Since the firm rents the capital from households, the problem of the firm is a sequence of static profit maximizing problems. In equilibrium, there are no profits, and each input is priced according to its marginal product, i.e.:

utkt:αytutkt=rt$${{u}_{t}}{{k}_{t}}:\alpha \frac{{{y}_{t}}}{{{u}_{t}}{{k}_{t}}}={{r}_{t}}$$ht:(1α)ytht=wt$${{h}_{t}}:\left( 1-\alpha \right)\frac{{{y}_{t}}}{{{h}_{t}}}={{w}_{t}}$$

In addition, using the link between energy, capital, and utilization we can express output as follows:

yt=Atht1α[kt(11v1)et1v1(v1v0)1v1]α$${{y}_{t}}={{A}_{t}}h_{t}^{1-\alpha }{{\left[ k_{t}^{\left( 1-\frac{1}{{{v}_{1}}} \right)}e_{t}^{\frac{1}{{{v}_{1}}}}{{\left( \frac{{{v}_{1}}}{{{v}_{0}}} \right)}^{\frac{1}{{{v}_{1}}}}} \right]}^{\alpha }}$$

The equation specifies output as a function of labor, capital, and energy, showing the direct effect of energy on output.

Note that if the depreciation rate is held constant, then the transmission of energy price shocks is restrained only to the effect of energy input on output through the production function channel. However, when depreciation rate is endogenous and depends on the utilization of physical capital, and then in turn through it on the use of energy, then energy has an additional indirecteffect on output, which operates through the capital stock. As we show later in the paper, the combination of those direct and indirect effects produces important difference in the dynamics of model variables

Government

In the model setup, the government is levying taxes on labor and capital income, as well as consumption, in order to finance spending on wasteful government purchases, and government transfers. The government budget constraint is as follows:

gtc+gtt=τcct+τy[wtht+rtutkt]$$g_{t}^{c}+g_{t}^{t}={{\tau }^{c}}{{c}_{t}}+{{\tau }^{y}}\left[ {{w}_{t}}{{h}_{t}}+{{r}_{t}}{{u}_{t}}{{k}_{t}} \right]$$

Tax rates and government consumption-to-output ratio would be chosen to match the average share in data, and government transfers would be determined residually in each period so that the government budget is always balanced.

Dynamic Competitive Equilibrium (DCE)

For a given process followed by technology {A}t=0$\left\{ A \right\}_{t=0}^{\infty }$average tax rates initial capital stock, k0 the decentralized dynamic competitive equilibrium is a list of sequences {ct,it,kt,ut,et,ht}t=0$\left\{ {{c}_{t}},{{i}_{t}},{{k}_{t}},{{u}_{t}},{{e}_{t}},{{h}_{t}} \right\}_{t=0}^{\infty }$for the household, a sequence of government purchases and transfers {gtc,gtt}t=0,$\left\{ g_{t}^{c},g_{t}^{t} \right\}_{t=0}^{\infty },$and input prices {pt,wt,rt}t=0$\left\{ {{p}_{t}},{{w}_{t}},{{r}_{t}} \right\}_{t=0}^{\infty }$that (i) the household maximizes its utility function subject to its budget constraint; (ii) the representative firm maximizes profit; (iii) government budget is balanced in each period; and (iv) all markets clear.

Data and Model Calibration

To characterize business cycle fluctuations with an endogenous depreciation rate in Bulgaria, we will focus on the period following the introduction of the currency board (1999–2016). Quarterly data on output, consumption, and investment was collected from National Statistical Institute (2017), while the real interest rate is taken from Bulgarian National Bank Statistical Database (2017). The calibration strategy described in this section follows a long-established tradition in modern macroeconomics: first, as in Vasilev (2016), the discount factor, β= 0.982, is set to match the steady-state capital-to-output ratio in Bulgaria, k/y = 13.964, in the steady-state Euler equation. The labor share parameter, 1 − α = 0.571, is obtained as in Vasilev (2017d), and equals the average value of labor income in aggregate output over the period 1999–2016. This value is slightly higher as compared to other studies on developed economies, due to the overaccumulation of physical capital, which was part of the ideology of the totalitarian regime that was in place until 1989. Next, the average income tax rate was set to τy = 0.1. This is the average effective tax rate on income between 1999 and 2007, when Bulgaria used progressive income taxation, and equal to the proportional income tax rate introduced as of 2008. Similarly, the tax rate on consumption is set to its value over the period, τc = 0.2.

In terms of parameters characterizing the household’s preferences, the curvature of the utility function is set to σ = 2, as in Hansen and Singleton (1983). Note that this parameter does not enter steady-state computation, and only affects cyclical fluctuations. Next, the relative weight attached to the utility out of leisure in the household’s utility function, ψ, is calibrated to match that in steady-state consumers would supply one-third of their time endowment to working. This is in line with the estimates for Bulgaria (Vasilev 2017a) as well over the period studied. In sum, the steady-state depreciation rate of physical capital in Bulgaria, δ = 0.013, was taken from Vasilev (2016). It was estimated as the average quarterly depreciation rate over the period 1999–2014. In addition, the steady-state capital utilization rate is normalized to unity, thus ω0 = 0.013. The curvature parameter, ω1, does not enter the steady state, and only matters for cyclical fluctuations. As in Finn (2000), we set ω1 = 1.25. Next, the scale parameter ν0 set to average value of energy-to-physical capital ratio, e/k. Again, the curvature parameter of the energy-utilization function, ν1, does not enter the steady state, and only matters for cyclical fluctuations. As in Finn (2000), we set ν1 = 1.61.

Finally, the processes followed by TFP processes and energy prices are estimated from the detrended series by running an AR(1) regression and saving the residuals. Tab. 1 summarizes the values of all model parameters used in the paper.

Model Parameters

ParameterValueDescriptionMethod
β0.982Discount factorCalibrated
α0.429Capital shareData average
1 − α0.571Labor shareCalibrated
ψ0.873Relative weight attached to consumptionCalibrated
σ2.000Curvature parameter, utility functionSet
δ0.013Depreciation rate on physical capitalData average
ω00.013Scale parameter, depreciation functionCalibrated
ω11.250Curvature parameter, depreciation functionSet
ν00.0143Scale parameter, energy utilization functionData average
ν11.610Curvature parameter, energy utilization functionSet
τy0.100Average tax rate on incomeData average
τc0.200VAT/consumption tax rateData average
A0.604Steady-state value of TFP processCalibarated
p1.000Steady-state energy price levelCalibrated
ρa0.701AR(1) persistence coefficient, TFP processEstimated
ρp0.980AR(1) persistence coefficient, energy price processEstimated
σa0.044st. error, TFP processEstimated
σp0.013st. error, energy processEstimated
Steady State

Once the values of model parameters were obtained, the steady-state equilibrium system solved, the “big ratios” can be compared to their averages in Bulgarian data. The results are reported in Tab. 2. The steady-state level of output was normalized to unity (hence the level of technology A differs from 1, which is usually the normalization done in other studies), which greatly simplified the computations. Next, the model matches consumption-to-output and government purchases ratios by construction. The investment ratios are also closely approximated, despite the closed-economy assumption and the absence of foreign trade sector. The shares of income are also identical to those in data, which is an artifact of the assumptions imposed on functional form of the aggregate production function.

Data Averages and Long-run Solution

VariableDescriptionDataModel
ySteady-state outputN/A1.000
c/y(non-energy) Consumption-to-output ratio0.6240.624
i/yInvestment-to-output ratio0.2010.175
pe/yEnergy consumption-to-output ratio0.1510.151
gt/yGovernment transfers-to-output ratio0.2200.149
wh/yLabor income-to-output ratio0.5710.571
ruk/yCapital income-to-output ratio0.4290.429
hShare of time spent working0.3330.333
r¯After-tax net return on capital0.0140.016

The after-tax return, where r¯ = (1 − τy )r − δ is also relatively well-captured by the model. Lastly, given the absence of debt, and the fact that transfers were chosen residually to balance the government budget constraint, the result along this dimension is understandably not so close to the average ratio in data.

Out of Steady-State Model Dynamics

Since the model does not have an analytical solution for the equilibrium behavior of variables outside their steady-state values, we need to solve the model numerically. This is done by log-linearizing the original equilibrium (non-linear) system of equations around the steady state. This transformation produces a first-order system of stochastic difference equations. First, we study the dynamic behavior of model variables to an isolated shock to the total factor productivity process, and then we fully simulate the model to compare how the second moments of the model perform when compared against their empirical counterparts.

Impulse Response Analysis

This subsection documents the impulse responses of model variables to a 1% surprise innovation to technology, as well as an unexpected 1% change in energy prices.

This price is to be interpreter as an aggregate category, comprosing electricity, coal, natural gas, and petroleum.

The impulse response functions (IRFs) are presented in Fig. 1 and 2.

Fig. 1

Impulse responses to a 1% surprise innovation in technology.

Fig. 2

Impulse responses to a 1% surprise innovation in energy price.

Impulse Responses to Technology Shocks

As a result of the one-time unexpected positive shock to total factor productivity, output increases upon impact. This expands the availability of resources in the economy, so the use of output consumption, investment, energy use and government consumption also increase contemporaneously. At the same time, the increase in productivity increases the after-tax return on the two factors of production, labor, and capital. The representative households then respond to the incentives contained in prices and start accumulating capital, and supplies more hours worked. In turn, the increase in capital input feeds back in output through the production function and that further adds to the positive effect of the technology shock. Lastly, the utilization rate increases as well, following the increase in the return on capital, but this also increases the endogenous depreciation rate. In the labor market, the wage rate increases, and the household increases its hours worked. In turn, the increase in total hours further increases output, again indirectly.

Over time, as capital is being accumulated, its after-tax marginal product starts to decrease, which lowers the households’ incentives to save. As a result, physical capital stock eventually returns to its steady-state and exhibits a hump-shaped dynamics over its transition path. The rest of the model variables return to their old steady states in a monotone fashion as the effect of the one-time surprise innovation in technology dies out.

Impulse Responses to Unanticipated Energy Prices

As a result of an unexpected one-time increase in the price of the aggregate energy input, illustrated in Fig. 2, the consumption of energy decreases, while its substitute, the non-energy private consumption, increases. Due to the relative scarcity of energy, illustrated in the increased valuation of energy, capital utilization rate increases. In turn, the time-varying endogenous depreciation rate increases as well, which in turn decreases capital accumulation and investment. As a result of lower capital availability, real interest rate goes up.

Next, from the complementarity of capital and labor in the Cobb-Douglas production function, hours fall as well, and the wage rate in the economy increases. Interestingly, aggregate output falls as well upon impact of the energy price shock, so an increase in energy prices is akin to a negative productivity shock, as energy could be expressed as a direct input in the production function. Next, as government consumption and transfers follow private out-put, both fall as well. Over time, all variables return to their steady state, but the negative effects from one-time unexpected increases in energy prices has a long-term negative effect on the economy.

Simulation and Moment-Matching

As in Vasilev (2017b), we will now simulate the model 10,000 times for the length of the data horizon. Both empirical and model simulated data are detrended using the Hodrick-Prescott (1980) filter. Tab. 3 summarizes the second moments of data (relative volatilities to output, and contemporaneous correlations with output) versus the same moments computed from the model-simulated data at quarterly frequency.

The model-predicted 95% confidence intervals are available request.

To minimize the sample error, the simulated moments are averaged out over the computer-generated draws. As in Vasilev (2016, 2017b, 2017c), the model matches quite well the absolute volatility of output and investment. By construction, government consumption in the model varies as much as output. However, the model in this paper underestimates the variability in consumption, due to the presence of energy consumption, which

Business Cycle Moments

σyDataModel
σcy0.550.14
σiy1.771.97
σgy1.211.00
σhy0.630.63
σwy0.830.39
σy/hy0.860.39
σuy3.220.63
σwh1.321.61
corr(c, y)0.850.47
corr(i, y)0.610.75
corr(g, y)0.311.00
corr(h, y)0.490.96
corr(w, y)-0.010.97
corr(u, y)-0.47-0.96
corr(h, y/h)-0.140.92

acts as a substitute for non-energy consumption. Still, the model is qualitatively consistent with the stylized fact that consumption generally varies less than output, while investment is more volatile than output.

With respect to the labor market variables, the variability of employment predicted by the model is almost identical to that in data, but the variability of wages in the model is much lower than that in data. This is yet another confirmation that the perfectly competitive assumption, e.g. Vasilev (2009), does not describe very well the dynamics of labor market variables. In addition, as in Vasilev (2017b, 2017c), the model fails in matching unemployment volatility, which in this model varies as much as the employment rate.

The reason behind this mismatch could be driven by several possible explanatory factors: the fact that the model misses the “out-of the-labor-force” segment, as well as the significant emigration to the older EU member states.

Next, in terms of contemporaneous correlations, the model systematically over-predicts the pro-cyclicality of the main aggregate variables such as consumption, investment, and government consumption.

This, however, is a common limitation of this class of models. However, along the labor market dimension, the contemporaneous correlation of employment with output, and unemployment with output, is relatively well-matched. With respect to wages, the model predicts strong cyclicality, while wages in data are acyclical. This shortcoming is well-known in the literature and an artifact of the wage being equal to the labor productivity in the model.

In the next subsection, as in Vasilev (2016), we investigate the dynamic correlation between labor market variables at different leads and lags, thus evaluating how well the model matches the phase dynamics among variables. In addition, the autocorrelation functions (ACFs) of empirical data, obtained from an unrestricted VAR(1) are put under scrutiny and compared and contrasted to the simulated counterparts generated from the model.

Auto- and Cross-Correlation

This subsection discusses the auto-(ACFs) and cross-correlation functions (CCFs) of the major model variables. The coefficients such as empirical ACFs and CCFs at different leads and lags are presented in Tab. 4 against the averaged simulated AFCs and CCFs. Following Canova (2007), this is used as a goodness-of-fit measure.

Autocorrelations for Bulgarian data and the model economy

k
MethodStatistic0123
Datacorr(ut, ut − k)1.0000.7650.5520.553
Modelcorr(ut, ut − k)1.0000.9550.9010.837
(s.e.)(0.000)(0.027)(0.051)(0.073)
Datacorr(nt, ntk)1.0000.4840.0090.352
Modelcorr(nt, nt − k)1.0000.9550.9010.837
(s.e.)(0.000)(0.027)(0.051)(0.074)
Datacorr(yt, yt − k)1.0000.8100.6630.479
Modelcorr(yt, yt − k)1.0000.9550.9010.836
(s.e.)(0.000)(0.027)(0.050)(0.073)
Datacorr(at, at − k)1.0000.7020.4490.277
Modelcorr(at, at − k)1.0000.9550.9000.836
(s.e.)(0.000)(0.026)(0.050)(0.072)
Datacorr(ct, ct − k)1.0000.9710.9520.913
Modelcorr(ct, ct − k)1.0000.9550.9030.845
(s.e.)(0.000)(0.026)(0.050)(0.073)
Datacorr(it, it − k)1.0000.8100.7220.594
Modelcorr(it, it − k)1.0000.9540.9010.841
(s.e.)(0.000)(0.026)(0.050)(0.073)
Datacorr(wt, wt − k)1.0000.7600.7830.554
Modelcorr(wt, wt − k)1.0000.9200.9000.836
(s.e.)(0.000)(0.026)(0.050)(0.073)

As seen from Tab. 4, the model compares relatively well vis-a-vis data. Empirical ACFs for output and investment are slightly outside the confidence band predicted by the model, while the ACFs for total factor productivity and household consumption are well-approximated by the model. The persistence of labor market variables are also relatively well-described by the model dynamics. Overall, the model with energy-utilization channel generates too much persistence in output and both employment and unemployment and is subject to the criticism in Nelson and Plosser (1992), Cogley and Nason (1995) and Rotemberg and Woodford (1996b), who argue that the RBC class of models do not have a strong internal propagation mechanism besides the strong persistence in the TFP process. In those models, e.g. Vasilev (2009), and in the current one, labor market is modelled in the Walrasian market-clearing spirit, and output and unemployment persistence is low.

Next, as seen from Tab. 5, over the business cycle, increase in data labor productivity leads to employment. The model, however, cannot account for this fact.

As in the standard RBC model a technology shock can be regarded as a factor shifting the labor demand curve, while holding the labor supply curve constant. Therefore, the effect between employment and labor productivity is only a contemporaneous one.

Dynamic correlations for Bulgarian data and the model economy

k
MethodStatistic-3-2-10123
Datacorr(nt, (y/n)t − k)-0.342-0.363-0.187-0.1440.4750.4700.346
Modelcorr(nt, (y/n)t − k)0.1230.1950.2920.9180.2880.2210.171
(s.e.)(0.320)(0.286)(0.250)(0.146)(0.243)(0.281)(0.317)
Datacorr(nt, wt − k)0.3550.4520.4470.328-0.040-0.390-0.57
Modelcorr(nt, wt − k)0.1230.1950.2920.9180.2880.2210.171
(s.e.)(0.320)(0.286)(0.250)(0.146)(0.243)(0.281)(0.317)
Conclusions

We introduce a pro-cyclical endogenous utilization rate of physical capital stock into a real business cycle model augmented with a detailed government sector. We calibrate the model to Bulgarian data for the period following the introduction of the currency board arrangement (1999–2016). We investigate the quantitative importance of the endogenous depreciation rate, and the capital utilization mechanism working through the use of energy for cyclical fluctuations in Bulgaria. In particular, a positive shock to energy prices in the model works like a negative technological shock. Allowing for variations in factor utilization and the presence of energy as a factor of production improves the model performance against data, and in addition this extended setup dominates the standard RBC model framework with constant depreciation and a fixed utilization rate of physical capital (e.g., Vasilev, 2009).

Still, the failure of the model along the labor market dimension – the high pro-cyclicality of wages and the low variability of the price of labor relative to that observed in data – suggest that the setup should depart from the perfectly competitive paradigm. As a suggestion for future research, the model should focus on the important frictions in the labor market, which forms almost two-thirds of total income (and much quantitatively much more important than the share of capital and energy), and extend the model along the lines of Vasilev (2016, 2017b, 2017c).