Public Debt Sustainability Under Uncertainty
An Invariant Set Approach

Contributor Notes

Author’s E-Mail Address:

The paper offers an approach to assessing the sustainability of public debt taking into account the effect of fiscal policy on output, as well as uncertainty in the model parameters and system dynamics. Uncertainty is specified in general terms, and the analysis is based on the notion of invariant sets. Examples are provided to illustrate how the method can be applied in practice.


The paper offers an approach to assessing the sustainability of public debt taking into account the effect of fiscal policy on output, as well as uncertainty in the model parameters and system dynamics. Uncertainty is specified in general terms, and the analysis is based on the notion of invariant sets. Examples are provided to illustrate how the method can be applied in practice.

1 Introduction

The question of what level of public debt a country can afford is central to fiscal policy. Vast literature is devoted to debt sustainability and an array of tools has been developed to guide assessment. A comprehensive review of previous work on the subject is beyond the scope of this paper and is not attempted here. We only mention a few studies to provide a context for the subsequent discussion and illustrate some of the difficulties that arise in the analysis.

Most of the earlier contributions have focused on ensuring solvency, interpreted as compliance with the government’s intertemporal budget constraint. For example, Blanchard et al. (1990) derived a simple sustainability criterion according to which the current level of debt should not exceed the present value of all future primary surpluses. This condition follows directly from the dynamic equation describing the evolution of debt and is relatively mild in the sense that it only requires the discounted value of future debt to tend to zero. Thus, debt trajectories that increase indefinitely as time goes to infinity are not ruled out as long as the rate of increase is lower than the interest-growth differential. Variants of this solvency condition have been widely used in the literature, e.g., EU Commission (2006) and Escolano (2010), among others.

The solvency approach to debt sustainability has been often criticized for not being particularly stringent as it allows for large primary deficits to be run initially, provided that future surpluses are high enough to offset them. However, there is no guarantee that the government would follow such a path; moreover, large corrections are costly both politically and economically. There is usually a limit up to which the government can increase taxes before collections start to decline due to the distortions introduced (Laffer curve effect). Similarly, expenditures cannot be reduced below a certain threshold consistent with the minimum level necessary for the smooth functioning of government. Therefore, considerations regarding the feasibility of adjustment are important and they have been reflected, for example, in IMF’s definition of debt sustainability (IMF, 2002).

For practical purposes, debt sustainability analysis (DSA) is often carried out by comparing a country’s debt level– typically measured as a ratio to GDP– to some indicative benchmark value (the so-called threshold approach). Benchmarks can be derived in various ways. The solvency condition itself defines such a limit if the future flow of revenues and non-interest expenditures, hence the primary balance, is constrained. Historical episodes of debt crises may also provide some guidance. Empirical work has been carried out, mostly based on signal-to-noise ratios, to identify the level of debt above which the likelihood of a crisis rises significantly (Baldacci et al., 2011).1 The IMF DSA framework is also underpinned by debt limits derived based on the signal approach.

Notwithstanding the difficulties associated with the identification of thresholds, the approach has gained a lot of traction and many countries have adopted fiscal rules containing explicit debt benchmarks. According to IMF (2016a), about 70 countries currently have fiscal frameworks that include caps on gross public debt. In a number of cases, these reflect rules established at the supranational level such as the Maastricht convergence criteria for the EU member states, for example. Since debt cannot be directly controlled by the government, often debt ceilings are complemented with rules about deficits or expenditures.

An alternative to the threshold approach is the sustainability definition proposed by Arrow et al. (2004), which, applied to indebtedness, can be interpreted as a requirement for the net worth of the entity of interest to be non-decreasing. Cuerpo et al. (2015) and IMF (2016b) apply this sustainability framework to households and non-financial corporations. Wyplosz (2011) discusses it in relation to government debt, pointing out that the concept can be made operational by requiring the debt to GDP ratio to be stationary, and since stationarity is difficult to ascertain, the ratio should follow a declining trend, allowing only for temporary increases.

Another challenge to the traditional DSA is uncertainty. The solvency condition depends on the future dynamics of primary balances, interest rates and growth rate, none of which are known for sure. Celasun et al. (2007), building upon Garcia and Rigobon (2005), address this issue by developing a methodology which accounts for the risks surrounding the debt dynamics. Their approach is probabilistic and uses fan charts to illustrate the uncertainty around the central projection of debt. The proposed algorithm comprises an estimated fiscal reaction function, a vector autoregression (VAR) which models the non-fiscal determinants of debt dynamics (e.g., GDP growth rate, interest and exchange rates), and a debt equation. In stochastic simulations of the debt path, shocks are drawn from a distribution which has the same covariance matrix as the covariance matrix of the VAR residuals and forecasts of the non-fiscal variables are consistent with these shocks. In a similar vein, Cherif and Hasanov (2012) examine how macroeconomic shocks affect the debt dynamics in a VAR framework with debt feedback. A different approach is taken by Barnhill and Kopits (2003) who use Value at Risk techniques to estimate the risk-adjusted net worth of the government.

IMF’s revised DSA methodology (IMF, 2011 and 2013) addresses some of the shortcomings of earlier methods. In particular, the definition of sustainable public debt is broadened to include a number of conditions, notably that “the primary balance needed to at least stabilize debt under both the baseline and realistic shock scenarios is economically and politically feasible, such that the level of debt is consistent with an acceptably low rollover risk and with preserving potential growth at a satisfactory level.” 2 The IMF DSA provides a set of tools to assess the realism of macroeconomic projections, vulnerabilities arising from the structure of debt, sensitivity to macro-fiscal shocks, and impact from realization of contingent liabilities. Stochastic simulations and fan charts showing ranges for the possible evolution of the debt ratio are also part of the toolkit.

This paper contributes to the literature by offering an alternative framework for assessing debt sustainability which features the following characteristics.

First, the proposed approach explicitly takes into account the effect of fiscal policy on the output gap and makes it endogenous in the analysis. Since fiscal policy can affect demand, large fiscal contractions, especially over extended periods of time, may not be feasible or desirable, including due to possible hysteresis effects, whereby long and deep recessions destroy the economy’s productive potential. The IMF DSA also recommends that the potential impact of fiscal adjustment on growth and interest rates is factored in, but this impact is not fully endogenized.

Second, the methodology provides an estimate of a debt threshold which is entirely forward-looking and depends on assumptions about fiscal multipliers, interest and growth rates, as well as uncertainties around the future paths of these variables. Thus, the threshold is system-specific and does not rely on estimation and/or averaging based on past data as in other works. In addition, the method delivers a fiscal reaction function in the form of a linear feedback rule which prescribes how the primary balance should respond to debt and the output gap. Unlike most of the literature where fiscal reaction functions are estimated empirically based on past data, our rule is normative and is designed to stabilize the economy given the parameters of the model and system disturbances.

Third, uncertainty is inherent in the analysis– it affects the system not only through additive shocks to the dynamics, but also through key model parameters. Uncertainty is described in rather general terms, reflecting the view that a probabilistic representation may not be available (especially since the focus is on the future) or an outcome “on average” may not be desirable. Instead, the decision-maker would seek policies that perform well under a range of possibilities, for example, when fiscal multipliers fall in the interval [0.5,1.5], or when future interest rates could be anywhere between 3 and 7 percent. A similar type of bounds can be assumed for additive shocks, either based on historical experience or other relevant information. In this respect, the analysis is closely related to the literature on robustness (Hansen and Sargent, 2008). A robust approach can be justified based on model misspecification as decision makers seldom have complete knowledge of the underlying dynamics.

Fourth, the framework assumes an infinite time horizon. This is another important difference from the prevailing approaches to debt sustainability under uncertainty. Typically, DSAs fix a specific time frame and examine the debt path during this period. For example, IMF uses a 6-year horizon, consistent with the availability of World Economic Outlook projections. As argued by Wyplosz (2011), finite horizons are a major constraint since in theory the debt sustainability concept requires infinite time. In fact, the stochastic approach to debt sustainability works well only over relatively short intervals. For instance, in a fan chart presentation, uncertainty around the central debt path typically becomes very large after 5-10 periods, rendering the tool of little practical value for longer horizons.

The analysis in this paper is based on Lyapunov’s stability theory for dynamical systems and the related notion of set invariance. In a nutshell, to capture the effect of fiscal policy on growth, the standard debt equation is complemented with an equation describing the evolution of the output gap. Formally, in each period the two variables representing the state of the system – output gap and the debt ratio–are viewed as a point in the two-dimensional space. We define debt as sustainable if it does not increase indefinitely, i.e., if it remains bounded at all times. This definition can be operationalized as follows: find a set of points with the property that whenever the initial state falls into this set, it stays there forever. In other words, we are interested in all possible combinations of output gaps and debt levels for which it is guaranteed that debt will not embark on an explosive path when subjected to shocks of a given size, also taking into account uncertainty in the system parameters. Such sets, if they exist, will be called invariant. We note that while invariant set membership is a sufficient condition for debt to remain bounded, it is not necessary. Even if the initial point is outside the invariant set, or if the state is pushed out of it later by a large one-off shock, it would still be possible to stabilize the system. This is so because, in the case of linear dynamics the invariant sets are also attracting, so all trajectories that start outside the set tend to approach it as time tends to infinity (Khlebnikov et al., 2011). It cannot be guaranteed, however, that in the process of convergence some underlying constraints, e.g., on the size of the fiscal adjustment in any given period, will be satisfied.

Finally, it is important to stress that assessing the sustainability of public debt is a complex problem and the simple model presented here is not meant to substitute a full-fledged DSA. In particular, it cannot capture many of the qualitative elements involved in such assessments. Still, the invariant set method can be useful as a supplementary tool to more elaborate frameworks.

The rest of the paper is organized as follows. Section 2 presents the model and describes the methodology for debt sustainability assessment; section 3 illustrates how the methodology is applied to specific examples, and section 4 concludes. Appendix A contains all the technical material, including proofs of some of the results used in the main text.

2 Analytical framework

To study debt sustainability when the effect of fiscal policy on growth is explicitly taken into account, we consider the following system of equations:


where yt denotes the output gap at time t, defined as yt=Yt/YtP1, with Yt and YtP denoting actual GDP and potential GDP, respectively. The output gap is assumed to follow a first order autoregressive process. Without a fiscal action, and assuming that 0 < ρ < 1, the output gap would eventually close, with degree of persistence determined by the value of ρ. Although in practice the output gap is not directly observable, the level of current economic activity relative to its natural level is an important consideration in economic decision making, and various techniques have been developed to measure potential output and deviations from it.

A key assumption of the model is that the pace at which the output gap closes can be influenced by a change in the primary deficit ut. The effect of the fiscal action on the output gap depends on the size of the fiscal multiplier αt. An increase in the primary deficit relative to the baseline (here assumed to be zero), however, increases the debt level dt, defined as the ratio of public debt to potential GDP.3 Besides the primary deficit, the evolution of the debt ratio depends on the growth rate of potential output gt and the interest rate rt. Both output gap and public debt are subject to exogenous shocks captured by the additive terms wi,t, i = 1, 2 (scaled by the parameters δi, i = 1, 2). No probabilistic assumptions are made about the nature of the shocks; they are only assumed to belong to some compact set. It is convenient to work with ellipsoidal sets of the type Wt={wt:wtWt1wt1} given their analytical advantages and potential links to statistical inference. In particular, the shape matrix of the disturbance ellipsoid can be calibrated to an estimated covariance matrix of shocks.4 This is a generalization of the standard assumption in the engineering literature that ||w′twt|| ≤ 1.

Further to the additive uncertainty, there is a strong case for introducing uncertainty in the model parameters. Indeed, since debt sustainability assessments are forward-looking and the evolution of the debt ratio depends critically on future interest and growth rates that are not known precisely, it appears natural to allow these parameters to vary within a certain range, e.g., rt ∈ [rl, rh], gt ∈ [gl, gh], where subscripts “l” and “h” stand for “low” and “high”. Similarly, the fiscal multiplier is uncertain and most likely not constant, with higher values typically reported during economic downturns and when interest rates are low. While ideally these variables would be modeled as functions of the state of the system (e.g., higher debt triggers higher interest rates), such an approach is technically very challenging due to the nonlinearities involved.

In matrix notation, system (1)- (2) can be written as:


For the matrices At and Bt we assume that Atco{A1, A2} and Btco{B1, B2} (here co denotes convex hull), with


The autoregressive parameter ρ is kept constant for simplicity. With these assumptions, it is enough to consider the convex hull of only two A matrices, rather than all four possible combinations of growth and interest rates; the intermediate cases (1 + rl) / (1 + gl) and (1 + rh)/(1 + gh) can be obtained as convex combinations of A1 and A2. The proposed framework remains applicable when an interval for ρ is specified as well but this would entail adding more constraints to the ensuing optimization problem (see below). Note that formulation (3) is more general than (1)- (2) in that it allows for the possibility of shocks to the output gap equation to affect also debt and vice versa, if matrix D is non-diagonal.

The main question of interest is under what combinations of growth, interest rates, fiscal multipliers, and shocks debt will remain sustainable. This, according to our definition, is equivalent to finding the largest set of points E (if it exists) that is invariant with respect to the dynamics described by (3). Invariance is understood in the following sense: for any initial x0E, the state xt is guaranteed to stay in E for all t > 0. In the context of debt sustainability, the existence of such a set implies that debt will not grow indefinitely even in the presence of uncertainties as specified above. As in the case of additive disturbances, we restrict the class of invariant sets to ellipsoids. This choice is motivated by two main reasons: first, any non-empty compact convex set can be reasonably well approximated by an ellipsoid, and second, ellipsoidal sets arise naturally from the application of Lyapunov’s stability theory to linear systems. The link can perhaps be best understood through a simple example in continuous time.

Consider the following system of ordinary differential equations defined on a domain Ω ⊂ R2:


and suppose that the system has an equilibrium at zero, i.e., f(0)=0. We want to know under what conditions the origin is a stable equilibrium in the sense that every trajectory that starts close to zero at t0 remains so for all t > t0. A famous result by A. Lyapunov states that if there exists a continuously differentiable function V: Ω → R, such that V(0) = 0, V(x) > 0, ∀x ∈ Ω\{0} and V˙(x)0x ∈ Ω, then x = 0 is stable. The last condition is equivalent to requiring that ∇V(x). f(x) ≤ 0, i.e., the gradient of V(x) and the velocity vector f(x) do not form an acute angle, implying that on the boundary of the sublevel set V(x) ≤ a, all velocity vectors point to the interior of the set when the angle is obtuse (Figure 1). Thus, within the invariant set a point can wander around, but once it gets close to the border it will move back.

Figure 1:
Figure 1:

Example of set invariance.

Citation: IMF Working Papers 2017, 057; 10.5089/9781475586251.001.A001

In physical systems, the Lyapunov function typically has the meaning of total energy and the stability condition requires that total energy decreases along every trajectory. In an economics context, one could interpret V(x) as the value function of an infinite-horizon linear-quadratic control problem whereby the decision-maker minimizes a loss function that penalizes deviations of the state from equilibrium. With this interpretation, the definition of sustainability adopted here is closely related to the criterion of Arrow et al. (2004). This criterion is derived from an optimization problem where a representative agent maximizes intertemporal utility. Sustainability is achieved when the value function associated with this problem is non-decreasing with time. Since in our setup we can think of the Lyapunov function as the value function of a minimization problem, it is natural to require that the loss function be non-increasing.

In a discrete-time framework, the condition V˙(x)0 is replaced with V(xt+1)V(xt) ≤ 0. Further, it is well known that for linear systems the Lyapunov function has the form:


where Q is some positive definite matrix. Note that the sublevel sets of V(x) are ellipsoids, so there is a close link between Lyapunov functions and the invariant ellipsoids introduced earlier. Specifically, since Q is positive definite, the set E={xt:xtP1xt1} where P−1 = Q defines an ellipsoid. In order for this ellipsoid to be invariant, for any xtE we need xt+1E as well.

If the matrix Q (or equivalently P) were known, we would only need to verify whether the current state falls inside the set E; if it does, then we can ascertain that the debt trajectories will remain bounded even in case of realization of the worst combination of parameters and shocks, as long as they are within the specified limits. All that it takes to stabilize the system is to adjust the primary deficit in accordance with a linear rule which reacts to the output gap and the debt ratio, i.e., to follow the rule ut = Kxt, where K in this case is a (1 × 2) matrix. The problem, however, is that the matrix Q is not known in advance and there can be many such matrices or none. In the latter case, the problem is not feasible, implying that debt may not be sustainable, and in the former case, we need a criterion to choose among the possible options. It appears reasonable to aim for the largest invariant set since typically policy-makers and markets are interested in the highest debt ratio that can be stabilized and eventually reduced. Thus, a natural criterion is to seek the invariant ellipsoid with the largest volume. Since the volume of an ellipsoid is proportional to the determinant of its shape matrix, this leads to a constrained maximization problem of the following kind:


subject to


Problem (4)-(7) is not standard in the sense that the variable of interest is not a scalar or vector but a matrix, and this matrix must possess certain properties, notably to be symmetric and positive definite, and to satisfy additional constraints arising from the nature of the problem. Such are constraints (5)-(6) which represent linear matric inequalities or LMIs (see Boyd et al., 1994).

The first LMI ensures that the ellipsoidal set defined by the shape matrix P−1 is invariant (see Appendix A). The second LMI is necessary because of economic considerations; it arises from the definition of the output gap. Clearly, neither actual, nor potential GDP can be negative, so as a minimum we have to impose the requirement y ≥ −1. More generally, it may be desirable to constrain the size of the output gap to some value ymax, i.e., to have ||yt|| ≤ ymax for all t. To incorporate this condition (which has to be symmetric in order to be able to use LMI techniques), it is convenient to introduce an auxiliary equation of the form


where zt isa vector that depends linearly on the state xt, and to establish a more general result based on the constraint ||zt|| ≤ zmax (see Proposition 2 in Appendix A). The constraint related to the output gap definition is a special case of inequality (6) with zmax = 1 and C being a (2 × 2) matrix with an entry of 1 in the (1,1) position and zeroes elsewhere.

In the above problem, Y is an auxiliary matrix which is related to matrix K appearing in the feedback law and τ is a free parameter taking values in (0,1). 5

In order for problem (4)-(7) to have a solution, the LMI constraints must be feasible. Even for simple problems, feasibility is often difficult to establish explicitly, so numerical methods are used. While a complete analytical characterization of feasibility in our case is challenging as well, it can be shown that if δ12w111τ>ymax2, i.e., when the shock on the output gap equation is large relative to the maximum admissible output gap, LMI (6) is not feasible and the problem does not admit a solution (Appendix A).

In summary, the invariant set approach to debt sustainability reduces the problem of stabilizing debt and output in the presence of uncertainty to that of finding a matrix P with certain properties, as described above. The existence of such a matrix ensures that the debt ratio will eventually converge, and moreover, if at any point in time the state of the system (represented by the combination of debt level and output gap) falls within the ellipsoidal set determined by P, all future states will remain within this set, as long as the shocks do not exceed the specified bounds. If in a given period there is a large one-off shock that takes the state outside the invariant set, the system would still be stabilizable; however, the required adjustment will likely be larger, and it will take longer to move to the equilibrium. On the other hand, if no such matrix exists, then the system cannot be stabilized and the debt ratio would likely grow indefinitely. An important feature of the approach is that whenever debt is sustainable it provides a fiscal rule – a linear function of the state variables which guarantees stabilization even if the worst admissible combination of shocks occurs.

Below we illustrate how this analytical framework can be used to determine the stability regions for specific values of the parameters involved. The numerical solutions have been obtained using CVX, a package for specifying and solving convex programs (Grant and Boyd, 2008 and 2013) and the figures have been generated with the aid of the Ellipsoidal Toolbox by A. Kurzhanski and P. Varaiya (2007).

3 Application

As an application of invariant set methods to debt sustainability, we consider three distinct scenarios depending on the interest-growth differential. Scenario (1) assumes that future growth rates are consistently higher than interest rates, that is gl > rh; Scenario (2) makes the opposite assumption, namely that rl > gh and finally, Scenario (3) considers an intermediate case where gl < rl < gh < rh. In all simulations, we assume that the fiscal multiplier α takes values in the interval [0.3, 1] and the autoregressive parameter ρ is equal to 0.7, values consistent with the typical findings in the literature.6 The shape matrix of the disturbance ellipsoid W−1 is diagonal with entries 0.012 and 0.022 which correspond to persistent shocks to the output gap of up to 1 percent of potential GDP per year and to debt of up to 2 percent of potential GDP per year. These values are purely illustrative and do not draw on any particular empirical work. Further, the matrix D is the identity matrix and for simplicity we fix the free parameter τ = 0.95. In principle, we cannot optimize with respect to τ but we can do a grid search to find the value that maximizes the objective. The results presented below, however, are not very sensitive to the choice of this parameter. In all simulations, the initial condition is set at (—0.02; 0.6), that is, the economy starts with a 2 percent negative output gap and 60 percent debt to potential GDP ratio. For the graphical presentation of the state trajectory, we take A = 0.5A1 + 0.5A2, B = 0.5B1 + 0.5 B2 and the disturbances are generated randomly from uniform distributions on [—0.01,0.01] for w1 and [−0.02,0.02] for w2, respectively.

3.1 Scenario 1: Growth rates higher than interest rates

In this scenario, we assume that future interest rates are consistently lower than growth rates: rl < rh < gl < gh. Specifically, r ∈ [0.03,0.045] and g ∈ [0.05,0.06]. The maximal invariant set and the simulated primary deficit path are shown in Figure 2. As one would expect, under these favorable interest and growth assumptions, debt sustainability is not an issue; initial debt can be 8 times larger than potential GDP, and the state would still remain in the invariant set. Although Scenario 1 may correspond to the economic realities that many countries currently face in an environment of historically low interest rates, it is unlikely that this situation will be sustained in the long run. Therefore, from a policy point of view, the case where the growth rate is always higher than the interest rate is perhaps the least interesting one.

Figure 2:
Figure 2:

Scenario 1: Invariant set for rl < rh < gl < gh.

Citation: IMF Working Papers 2017, 057; 10.5089/9781475586251.001.A001

3.2 Scenario 2: Interest rates higher than growth rates

Scenario 2 assumes that gl < gh < rl < rh. In this scenario, r ∈ [0.05,0.07], g ∈ [0.03,0.04] and it corresponds to a dynamically efficient economy. As one would expect, maintaining debt sustainability is more challenging in this case. Figure 3 confirms this; the invariance bounds for debt are more than twice as narrow as in the previous scenario. Moreover, the implied adjustment of the primary balance in the first period of 12 percent of potential GDP is unrealistically large. It is underpinned by the high (in absolute terms) coefficient on debt in the fiscal reaction function, which for this example has the form F = − 0.64y − 0.22d. For comparison, the coefficient on debt in fiscal reaction functions typically found in the empirical literature ranges between 0.02 and 0.10 (Berti et al., 2016).

Figure 3:
Figure 3:

Scenario 2a: Invariant set for gl < gh < rl < rh.

Citation: IMF Working Papers 2017, 057; 10.5089/9781475586251.001.A001

Given that adjustment of the magnitude suggested by the model is not feasible in most cases, introducing an additional constraint to limit the size of the primary balance change in each period seems warranted: || ut|| ≤ umax. 7 This gives rise to another LMI:


The proof that this LMI is equivalent to the constraint on the control variable follows the same steps as Proposition 2 in Annex A and can be found in Nazin et al. (2007), for example. In line with the literature (see IMF, 2016a), we choose the maximum adjustment to be 4 percent of potential GDP per year in the subsequent simulations.

As Figure 4 suggests, when the constraint on the primary deficit is taken into account, the invariant set shrinks further, so that the initial condition is no longer inside it. Still, eventually the state moves toward equilibrium (reflecting the attracting property of the invariant ellipsoid), but outside of the set the constraint on control is violated (right panel). Once the point gets into the invariant set all constraints are satisfied. The calculated feedback rule prescribes zero response to the output gap and deficit reduction of 0.08 percent of potential GDP for each percentage point of the debt ratio. If we relax slightly the constraint on the primary balance from 4 percent to 5 percent, the initial point falls into the stability region and all constraints are satisfied. A similar result is obtained if we reduce the shock to the debt equation from 2 percent to below 1 percent of potential GDP. If, on the other hand, the size of the admissible permanent shocks is doubled (i.e., matrix W−1 has diagonal elements 0.022 and 0.042), then an invariant set no longer exists and debt cannot be stabilized.

Figure 4:
Figure 4:

Scenario 2b: Invariant set for gl < gh < rl < rh, constrained control.

Citation: IMF Working Papers 2017, 057; 10.5089/9781475586251.001.A001

As an alternative to imposing an explicit limit on the improvement of the primary balance, one can opt for a tighter constraint on the output gap, such that at no point in time actual GDP is allowed to deviate from potential by more than a specified amount. Figure 5 presents the case when the maximum output gap is set at 15 percent. Under this assumption, the invariant set becomes smaller and the required adjustment to move the state to the stability region is again very large in the initial period (the fiscal reaction function is F = −0.61x1 − 0.30x2), so constraining the output gap does not appear to be a viable alternative to limiting the size of adjustment.

Figure 5:
Figure 5:

Scenario 2c: Invariant set for gl < gh < rl < rh, constrained output.

Citation: IMF Working Papers 2017, 057; 10.5089/9781475586251.001.A001

3.3 Scenario 3: Intermediate case

Finally, Scenario 3 considers the intermediate, and perhaps most realistic case where in some periods interest rates are higher than potential GDP growth rates and in others they are lower (gl < rl < gh < rh). Nominal growth under this scenario is assumed to be between 4 and 6 percent and interest rates between 5 and 7 percent. The results are shown on Figure 6. Again, the adjustment in the initial period is quite substantial, which calls for imposing constraints on the change in the primary balance. The maximal invariant set under the 4 percent constraint is presented on Figure 7.

Figure 6:
Figure 6:

Scenario 3a: Invariant set for gl < gh < rl < rh.

Citation: IMF Working Papers 2017, 057; 10.5089/9781475586251.001.A001

Figure 7:
Figure 7:

Scenario 3b: Invariant set for gl < rl < gh < rh, constrained control.

Citation: IMF Working Papers 2017, 057; 10.5089/9781475586251.001.A001

4 Conclusion

This paper discusses how invariant set techniques can be used to aid debt sustainability assessments. The method provides an estimate of a debt threshold and a fiscal reaction function which is guaranteed to stabilize debt taking into account bounded parameter uncertainty and exogenous shocks. In the case of linear systems, there is a close link between invariant ellipsoids and Lyapunov functions, which in turn can be related to the debt sustainability criterion proposed by Arrow et al. (2004).

Application of the invariant set approach leads to an optimization problem which involves a system of linear matrix inequalities. This problem can be efficiently solved using numerical methods and the paper provides several examples which differ in terms of key model parameters. Simulations suggest that debt sustainability issues arise in scenarios where future interest rates are consistently higher than potential growth rates and the size of the primary balance adjustment is constrained. In those scenarios, the resulting invariant sets are relatively tight; debt can still be stabilized (if shocks are not large) but at the cost of violating the constraints on adjustment.

Although the analysis in the paper is restricted to a simple system of two equations, the method has a rather general applicability. It can be used for an arbitrary linear system, e.g., an estimated VAR. Besides debt sustainability, invariant sets could be useful in addressing problems of stabilization in other areas, such as monetary policy. Deriving analytical criteria for the existence of invariant sets, however, remains a challenge.

Public Debt Sustainability Under Uncertainty: An Invariant Set Approach
Author: Rossen Rozenov