Next Article in Journal
Development Prospects of Tourist Passenger Shipping in the Polish Part of the Vistula Lagoon
Previous Article in Journal
Evaluating Polarity Trend Amidst the Coronavirus Crisis in Peoples’ Attitudes toward the Vaccination Drive
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Feasibility and Economic Impacts of the Energy Transition

1
iMMC (Institute of Mechanics, Materials and Civil Engineering), Université Catholique de Louvain, 1348 Louvain-la-Neuve, Belgium
2
Université Lille, CNRS, IESEG School of Management, UMR 9221-LEM-Lille Économie Management and IRES, Université Catholique de Louvain, 1348 Louvain-la-Neuve, Belgium
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(10), 5345; https://0-doi-org.brum.beds.ac.uk/10.3390/su13105345
Submission received: 4 March 2021 / Revised: 30 April 2021 / Accepted: 4 May 2021 / Published: 11 May 2021

Abstract

:
There is currently no consensus regarding whether or not renewable energies are capable of supplying all of our energy needs in the near future. To shed new light on this controversy, this paper develops a methodology articulating a macroeconomic model with two sectors (energy and non-energy) and an energy model that is able to calculate the maximum potentials of solar and wind energy. The results show that, in a business-as-usual context, a complete energy transition on a global scale is unachievable before the end of the century. The reason lies in the increasing capital needs of the energy sector, which slows, if not stops, economic growth and the energy transition. A complete transition can be achieved by 2070 provided that (i) energy demand is kept under control at its current level, (ii) a sufficient rate of capital growth is sustained (above its historical level), and (iii) substantial progress is made in terms of energy efficiency. However, this strategy requires a significant increase in the savings rate, with a negative impact on consumption, which ends up stagnating at the end of the transition.

1. Introduction

In recent years, a controversy has developed as to whether, within a few decades, renewable energies could meet all of the world energy needs, without jeopardising economic growth. Many articles have been published in the literature, but there is clearly no consensus regarding the answer to this question. The dominant position answers positively [1,2,3,4,5,6,7,8]: an energy transition (ET) to a fully renewable energy system is possible within a few decades while ensuring economic growth (the interested reader will find further contributions to the dominant position in the extensive literature reviews made by Moriarty and Honnery [9] and Capellán-Pérez et al. [10]). Such a transition would be made possible by (i) increased energy efficiency, (ii) technical progress increasing the quality of the means of energy production (e.g., solar panels), and (iii) a reduction in the manufacturing costs of these means (in particular, by exploiting the economies of scale allowed by mass production).
There are distinctions within this dominant position. Some studies propose scenarios where the ET is accompanied by an increase in energy consumption [7]. The development of renewable resources should then not only replace fossil resources, but also allow for further energy growth. Other studies present scenarios, which, on the contrary, are based on the possibility of an absolute decoupling between the production of goods and services and energy consumption. In this context, economic growth would be accompanied by a decrease or stagnation of energy consumption [4,11]. This decrease would be made possible thanks to technical progress, but also to more sober behaviours in terms of energy use.
The position that is presented in the previous paragraphs is dominant, but not unanimous. Various recent contributions [9,10,12,13,14], which we shall refer to by the expression “critical literature”, have highlighted that the ET could result in an exacerbation of the capital needs of the energy sector, with a possible crowding out effect on the rest of the economy regarding the allocation of investment. This effect would negatively affect economic growth and, if strong enough, the ET could be accompanied by a prolonged phase of (imposed) economic decline.
Some of these contributions are based on models that highlight the close relationship between the energy sector and the rest of the economy (the dominant position is also based on models, particularly of the IAM type. For a critique of the latter, see [10]). It is interesting to note that they are based on different, sometimes even opposing, methodologies and assumptions. Some models, such as GEMBA [15,16] and MEDEAS [10], belong to systems theory. Others belong to the economic theory of growth [17,18]. Unlike the former, they assume market equilibrium and explicitly model the behaviour of economic agents. Only the EENGM model of Režný and Bureš [19] was found to combine the two methodologies. This article also provides an interesting review of a set of energy-economic models. For another review of this type of models, see Rye and Jackson [20].
In all of these contributions, the crowding out effect mentioned above is associated with a decrease in energy efficiency (known as the acronym Energy Return On Investment (EROI)), that affects both non-renewable and renewable energy technologies. On the one hand, the EROI of fossil resources would decrease as they become increasingly difficult to extract. On the other hand, the EROI of renewable energy technologies appears to be (significantly) lower than the historical levels that are achieved by the EROI of fossil-based energy production. On top of that, despite technical progress, the EROI of renewables could decrease in the future due to a location effect linked to the lower "quality” of remaining available sites for additional installations to harvest these resources.
This location effect is particularly well highlighted by Dupont et al. [21,22] at the global level for wind and solar energy. As stated by Moriarty and Honnery [9], these two sources will play the main role in the upcoming ET.
The aim of this paper is to couple the macroeconomic model that was developed by Fagnart and Germain [23] and Dupont et al. [24], accounting for the crowding out effect, with the energy models of Dupont et al. [21,22], accounting for the location effect. The macroeconomic model, an input–output type model with two sectors (energy and non-energy), has the advantage of being relatively high-level (as compared to other much more detailed models) and to lead to transparent analytical results linking economic growth to the penetration rate of renewable energies. The energy models are able to calculate the maximum potentials of solar and wind energy on a global scale, from a very detailed description of the local potentials.
The articulation of the previous models allows us to answer the following questions: given a future trajectory of the energy demand, hypotheses on the evolution of the technical parameters of the economy and on the behaviour of economic agents in terms of consumption and savings, (i) is a complete ET possible and (ii), if so, at what pace? As a corollary, it is possible to determine the conditions that ensure that the ET is completed over a given period of time.
The structure of the article is as follows. Section 2.1 presents the macroeconomic model and EROI formula, Section 2.2 the energy model, and Section 2.3 and Section 2.4 the articulation between these two models. Section 2.5 presents the calibration of the macroeconomic model for a reference year, based on [24], and Section 2.6 formulates the hypotheses regarding the future evolution of some exogenous parameters. Section 3 presents the results of the various simulations of the ET, based on assumptions regarding (i) the behaviour of economic agents in terms of consumption and savings and (ii) the potential and pace of technical progress achievable in the coming decades. The end of the section discusses the main assumptions and their impacts on the simulation results.

2. Materials and Methods

2.1. The Macroeconomic Model

The macroeconomic model (MM) was extensively presented in a previous publication [24]. In the following subsections we recall its main characteristics for the understanding of the following work.

2.1.1. Main Features

The model (as illustrated in Figure 1) is that of a closed economy, with two sectors: (i) the energy sector (ES), which produces energy for the rest of the economy, and (ii) the non-energy sector (FS), which produces a general-purpose good used for intermediate and final consumptions and for investment. The output of the ES is distributed between energy intermediate consumptions (EIC) of the two sectors and deliveries for final consumption, final referring to private (household) as well as public (government) consumption. Both of the sectors use three factors of production: capital, energy, and goods (in the following, we will refer to the intermediate consumption of goods that are produced by the non-energy sector as non-energy intermediate consumption (NEIC)). The capital stock of the economy is shared between the two sectors according to their respective needs.

2.1.2. Mathematical Formulation

The energy and non-energy sectors are identified by the indices e and f, respectively. The model is dynamic in the sense that the variables are discrete functions of time. However, time-dependency is omitted in the equations to simplify the notations.
The sectoral outputs and factor consumptions are linked by the following relationships:
Y i = K i v i = E i q i = X i x i
where:
  • Y i , K i , E i , and X i refer to the production, capital stock, EIC, and NEIC of sector i, respectively. Energy production and consumption are measured in units of energy ( u e ); good production, NEIC, and capital stock in units of good produced by the FS ( u g ).
  • v i , q i , and x i refer to the capital, EIC, and NEIC intensities of sector i, respectively, i.e., the quantities of these inputs per unit of output. These intensities are likely to change over time as technologies evolve.
The total capital stock is shared between the two sectors:
K = K e + K f
Capital accumulation is described by:
Δ K = I δ K = I net
where I and I net refer to the gross and net investments, respectively. δ refers to the average depreciation rate of capital ( 0 < δ < 1 ). It is the weighted average of the sectoral depreciation rates:
δ = δ e K e + δ f K f K
Further, the growth rate of a given variable x is defined as g x = Δ x x . In particular, the capital growth rate is defined as:
g K = Δ K K
The production of the FS is divided into sectoral NEIC, deliveries for final consumption, and investment:
Y f = X e + X f + C f + I
The output of the ES is divided between sectoral EIC and final energy consumption (that is defined as private and public consumption of energy, similarly to the final consumption of goods):
Y e = E e + E f + C e
The total consumption is defined as the sum of final consumptions (of goods and energy):
C = C f + p C e
The price of the production good is chosen as the numeraire and is equal to 1. With p being the real price of energy (in u g / u e ), the sectoral gross added values (AV i ) are calculated as:
AV f = Y f X f p E f = [ 1 x f p q f ] Y f
AV e = p [ Y e E e ] X e = [ p ( 1 q e ) x e ] Y e
By definition, the gross domestic product (GDP) is the sum of the sectoral added values, i.e.:
GDP = AV f + AV e = [ 1 x f p q f ] Y f + [ p ( 1 q e ) x e ] Y e
The second equality derives from Equations (9) and (10) and it splits the GDP in terms of sectoral deliveries to the final demand. From Equations (6)–(11), one can easily rewrite Equation (11) to express the GDP as the sum of total consumption and investment:
GDP = C + I
The sectoral net added values (NV i ) are obtained from the gross sectoral added values by subtracting the sectoral consumptions of capital (understood as the depreciation of capital in the two sectors):
NV i = AV i δ i K i = AV i δ i v i Y i , i = e , f
Given Equations (9) and (10), one obtains:
NV f = [ 1 w f p q f ] Y f
NV e = [ p ( 1 q e ) w e ] Y e
where, by definition:
w i = x i + δ i v i , i = e , f
The technical parameters w i represent the consumption of goods of the two sectors per unit produced, in the form of NEIC and capital consumption.
By definition, the net domestic product (NDP) is the sum of the sectoral net added values. Given the previous equations, one obtains:
NDP = NV f + NV e = [ 1 w f p q f ] Y f + [ p ( 1 q e ) w e ] Y e
From Equations (3) and (12), one also obtains:
NDP = C + I net
NDP accounts for the useful output of the economy, which is used for total consumption and economic growth.
The macroeconomic savings rate is defined by:
s = I GDP
It represents the share of GDP invested in capital production, with the remaining going to final consumption. It follows that:
C = [ 1 s ] GDP
The share of energy in total consumption is defined by:
η = p C e C
Parameters s and η concisely describe savings and consumption choices.
A variable of interest for energy-economy models is the EROI and, in this paper, we will look at the point-of-use EROI. At the societal level, it corresponds to the ratio between the production of the ES that leaves the ES (i.e., the energy delivered, = Y e E e ) and all of its inputs expressed as energy flows (EIC E e , NEIC X e and capital δ e K e ). In Appendix A, it is shown that the point-of-use EROI at the societal level can be calculated, as follows:
EROI = Y e E e E e + E X e + E δ e K e = Y e ( 1 q e ) E e + q f w e Y e = 1 q e q e + q f w e
E X e and E δ e K e being the energy content of the goods flows X e and δ e K e .

2.2. The Energy Model

The energy model (EM) decomposes the ES into two subsectors: the renewable and non-renewable energy subsectors. The renewable energy subsector is based on the methodology that was developed by Dupont et al. [21,22]. It relies on the postulate, as advocated by Moriarty and Honnery [9], that wind and solar energy will represent the two main sources of energy in a 100% renewable world. This is a position shared by the International Energy Agence (IEA) and the U.S. Energy Information Administration (EIA) [25,26]. Indeed, although biomass is currently the leading source of renewable energy worldwide, it mainly concerns domestic use by the most deprived populations. Its intensive development for energy use is a source of serious controversy due to the potential competition with other essential uses (mainly food supply, but also biodiversity conservation, CO 2 capture, etc.). Regarding other renewable energy sources, either their potential for future development is limited (hydropower), or their technical potential is relatively low (wave and tidal energies), or their development does not seem as promising as expected (geothermal energy). These considerations are shared by [4], since, in their projection for 2050, biomass is absent, hydropower represents 4% of the total production, and the sum of geothermal, tidal, and wave energies represents just over 1% of the total.
The aim of the EM is to quantify the potential of these two energies at a global level.
A spatial discretisation of the world is at the root of the EM. The spatial resolution is 0.75° × 0.75° (4000 km 2 on average), which corresponds to a total of 114,000 terrestrial sites. For each terrestrial site, the area that is available for the implementation of wind or solar power plants is estimated by combining different databases, including, among others, the type of land use and the local slope of the land (which is a limiting factor for solar power plants), or the water depth and distance from the coast (for offshore wind). Based on the local resources (wind speed distribution and solar irradiation), the theoretical annual production of each technology considered is then calculated for each site.
The technologies considered are the following. Regarding wind energy, onshore and offshore wind up to 1000 m depth are taken into account. Among solar technologies, the ones taken into account are: (i) photovoltaic centralised power plants and residential installations (on rooftops), with 2 types of panels of different efficiencies and technologies, and (ii) 3 types of concentrated solar power plants (parabolic power plants with or without thermal storage, and solar towers with thermal storage). The wind and solar thermal technologies are optimised to be best adapted to the local site conditions. Wind turbines are optimised in terms of (i) their design to be best adapted to local climatic conditions and (ii) their spacing to reduce losses due to the wake effect (i.e., the interactions between successive rows of wind turbines). The design of thermal power plants (i.e., the size of the heliostat field in relation to the nominal power of the installation) is optimised to be best adapted to local sunlight conditions.
At each site, solar technologies (photovoltaic and thermal) compete for land use and the model selects the one with the best EROI. However, this is not the case between wind and solar power plants, as a given area can be occupied simultaneously by both types of installations.
Once the technological choices have been made at the level of each site, the EM calculates, per site and for each of the technologies used on the site, (i) the production potential, (ii) the direct and indirect energy inputs that are required per unit of energy produced and, by relating the two previous quantities, and (iii) the site EROI. On this basis, by ranking all combinations of sites and technologies in descending order of EROI, the EM provides the direct and indirect energy inputs E X r e and E K r e required producing a given aggregated amount of renewable energy Y r e (it is shown in Appendix B that the EROI ranking criterion leads to the same results as a ranking by profitability).
Thus, the EM provides two curves E X r e ( Y r e ) and E K r e ( Y r e ) over the interval [ 0 , Y r e m a x ] , where Y r e m a x designates the maximum production potential of wind and solar energy (potential being calculated by the EM, expressed in units of energy per year). These curves correspond to a discrete set of points, with each value corresponding to the potential of a given technology in a given geographical location. These two curves are used by the MM to calculate the technical parameters of the ES.

2.3. Articulation between the Two Models

The energy flows E X r e and E K r e are the energy contents (in units of energy u e ) of X r e and K r e , which refer, respectively, to the NEIC and capital consumption of the renewable energy subsector (in units of goods produced by the FS u g ). These quantities are related by:
X r e = E X r e q f
K r e = E K r e q f
where q f is the energy intensity (i.e., the energy content) of a unit of good that is produced by the FS.
It is then possible to calculate the technical coefficients of the renewable energy subsector, depending on its production:
q r e ( Y r e ) = E r e ( Y r e ) Y r e x r e ( Y r e ) = X r e ( Y r e ) Y r e v r e ( Y r e ) = K r e ( Y r e ) Y r e
Two assumptions are taken regarding the direct energy inputs E r e . First, internal energy inputs are not accounted for in the MM (see [24]). As wind and solar facilities only consume operationnal electricity that is directly produced by the installation, the direct energy input related to the production of renewable electricity is assumed to be zero. However, E r e should account for the losses in electricty that are associated with transport and distribution. According to the IEA statistics (https://www.iea.org/data-and-statistics, last acces on 10 May 2021) and MEDEAS [27], they represent about 9.5% of the electricity thta is consumed worldwide, i.e., 8.7% of the electricity produced. It gives that q r e = E r e Y r e = 8.7 % . The increase in distribution losses with renewable penetration (as described in the MEDEAS documentation) is not taken into account in the model, i.e., q r e remains constant during the energy transition.
The intensities x r e and v r e as a function of Y r e are shown in Figure 2. The increasing rate of these intensities highlights the location effect mentioned in the introduction. Because of the lower “quality” of the sites that can receive the additional installations exploiting wind and solar resources, increasing capital and NEIC are needed to produce one unit of energy.
The renewable energy penetration rate is defined by the ratio ϰ = Y r e Y e . If Y e is given, then the total renewable ( Y r e = ϰ Y e ) and non-renewable ( Y n r e = ( 1 ϰ ) Y e ) energy productions can be deduced. The technical parameters q r e , x r e , v r e can also be expressed as functions of ϰ .
The parameters of the non-renewable energy subsector are determined by the following ratios:
q n r e = E n r e Y n r e x n r e = X n r e Y n r e v n r e = K n r e Y n r e
If they are known, then the technical parameters of the ES are then determined as an average, weighted by ϰ , of the technical parameters of the two subsectors:
q e ( ϰ ) = ϰ q r e ( ϰ ) + ( 1 ϰ ) q n r e
x e ( ϰ ) = ϰ x r e ( ϰ ) + ( 1 ϰ ) x n r e
v e ( ϰ ) = ϰ v r e ( ϰ ) + ( 1 ϰ ) v n r e

2.4. Allowable Space and Transition Curve

2.4.1. Allowable Space

The concept of allowable space (AS) emerges from the attempt to determine an exhaustive set of trajectories satisfying a given criterion. In the present case, the retained criterion is simply that the net sectoral added values are both positive:
NV i 0 , i = e , f
This criterion reflects the fact that both sectors should contribute positively to the useful output of the economy (as measured by the NDP).
From this criterion, an allowable interval for the energy price p can be derived (Appendix C):
I p p ̲ = w e 1 q e < p < 1 w f q f = p ¯
From Equations (13) and (14), we know that NV f decreases and NV e increases with p. Beyond p ¯ , the price of energy is so high that NV f is negative. Below p ̲ , the price of energy is so low that NV e is negative.
Under assumptions regarding the economic choices in terms of consumption and investment, I p can be transformed into an interval for k, the capital/energy production ratio of the economy (Appendix D):
k = K Y e
I k = [ k ̲ ( x f , v f , q f , x e , v e , q e , z ) , k ¯ ( x f , v f , q f , x e , v e , q e , z ) ]
where z is a vector of behavioural parameters, which contains either s and η or s and g K , where g K refers to the capital growth rate. As was the case for the interval bounds for p, beyond k ¯ , NV f becomes negative (p is too high) and below k ̲ , NV e becomes negative (p is too low).
Among the technical parameters on which I k depends, v e , q e , and x e are functions of ϰ according to Equations (21)–(23), while the others do not vary with ϰ . Therefore, we can write Equation (27) as:
I k = [ k ̲ ( ϰ ) , k ¯ ( ϰ ) ]
k ̲ ( ϰ ) and k ¯ ( ϰ ) are the top and bottom envelopes of the allowable space (AS), defined as all the ( ϰ , k ) couples that satisfy the criterion from Equation (24). The envelopes k ¯ ( ϰ ) and k ̲ ( ϰ ) are both increasing functions of ϰ , partly because renewable energies are more NEIC and capital intensive than non-renewable energies. The equations determining these envelopes are specified in Appendices Appendix D.2 and Appendix D.4, depending on the assumed economic choices in terms of consumption and investment.
An example of AS is illustrated by the space between the lower ( k ̲ ( ϰ ) ) and upper ( k ¯ ( ϰ ) ) curves shown in Figure 3.

2.4.2. Transition Curve

A tansition curve (TC) is a function Φ such that:
k = Φ ( ϰ ) , ϰ [ 0 , 1 ]
derivable, increasing and contained in the AS. Consequently, a TC respects the following inequalities:
k ̲ ( ϰ ) Φ ( ϰ ) k ¯ ( ϰ ) , ϰ [ 0 , 1 ]
An example of TC is shown by the dotted curve presented in Figure 3. It is entirely located between the two envelopes and, therefore, respects the previous inequalities.

2.4.3. The Reduced Dynamic Model

The reduced dynamic model (RDM) is obtained from the MM, while integrating the notion of TC. It makes it possible to fully characterize the trajectory of an economy during its ET and, more precisely, to numerically determine whether (i) the ET is feasible or not in a finite time and (ii) if it is feasible, the time that is needed to complete it. Thanks to the MM, it is then possible to see the impact of the ET on all economic variables.
Two versions of the RDM are considered, which correspond to Exercises 1 and 2 described below, and that depend on different assumptions in terms of consumption and investment:
  • Exercise 1: s and η fixed.
  • Exercise 2: g K and η fixed.
Exercise 1 assumes a constant savings rate and a constant share of energy expenditure in total consumption. In this case, the capital growth rate g K is endogenous. Exercise 2, instead, assumes that the economy has a given capital growth rate g K to achieve at constant η , and it results in an endogeneous savings rate. Exercises 1 and 2 are based on the common assumption that η = p C e / C is constant. As part of this work, in order to test the sensitivity of the results to this hypothesis, a third exercise was analysed where s is constant and final good and energy consumption are assumed to be strictly complementary (in other words, the ratio C e / C f is constant). When compared to Exercise 1, in that case the substitution between final good and energy is less easy at the level of consumption choices. Exercise 1 and this third exercise are in fact polar cases in the sense that, for a given budget, they are respectively based on the assumption that the agents maximise either a Cobb-Douglas or a Leontief utility function (whose elasticities of substitution are, respectively, unitary and zero).
The results from Exercise 3 were most of the time favourable to the ET. For some TCs and in the framework of the sustainable development scenario, it can even be completed in a few decades. However, these results are questionable, as they imply a sharp rise of the price of energy during the ET, making the constant C e / C f hypothesis unlikely.
Exercise 1—s and η Fixed
Within the context of Exercise 1, it is possible, thanks to the MM, to express g K as a function of k, ϰ , s and the other parameters (Appendix D.3):
g K ( k , ϰ ) = s 1 η [ 1 s ] [ 1 x f v f [ 1 x f ] v e ( ϰ ) v f + x e ( ϰ ) k ] δ ( ϰ , k )
where
δ ( ϰ , k ) = δ f + [ δ e δ f ] v e ( ϰ ) k
The first term of the right-hand side member of Equation (30) is an increasing function of k, η , and s (because η < 1 ), and a decreasing one of ϰ . As the second term, δ ( ϰ , k ) , is also an increasing function of k and a decreasing one of ϰ , g K follows the same trends. It follows that, if g K is positive, K also increases, and, thus, for a given Y e , k also increases. Ceteris paribus, there exists a positive feedback loop between g K and k. However, this loop does not result in explosive growth, since, if k tends to , g K tends to a positive constant. On the other hand, g K is negatively affected by the ET. As ϰ increases, so do x e and v e , resulting in a lower share of the FS output being available for its own investment. This exacerbation of the needs of the ES (mainly in terms of capital) to the detriment of the FS, mathematically translated by the fact that g K ϰ < 0 , is the manifestation of the crowding out effect mentioned in the introduction.
Exercise 2— g K and η fixed
Within the context of Exercise 2, it is possible, thanks to the MM, to express s as a function of k, ϰ , g K , and the other parameters (Appendix D.5):
s ( k , ϰ ) = [ 1 η ] [ g K + δ ( ϰ , k ) ] 1 x f v f η [ g K + δ ( ϰ , k ) ] [ 1 x f v f v e ( ϰ ) + x e ( ϰ ) ] 1 k
where δ ( ϰ , k ) is defined by Equation (31). s is a decreasing function of k and an increasing function of g K and ϰ , provided that δ f is not too much higher than δ e . For a constant Y e , an economy with more capital produces more goods, which allows a lower relative savings rate to achieve the investment needed to sustain the desired growth rate g K . On the contrary, a higher renewable penetration rate diverts capital away from the FS, which results in less goods being produced and the need of more relative savings to reach g K .
The RDM for Exercise 1 consists of Equations (28)–(30), (21)–(23), and the identity defining g K (Equation (5)). The RDM for Exercise 2 consists of the same equations, except that Equation (32) replaces Equation (30).
Although the time dependency has not been explicited in the formulas, the RDM is a (discrete) dynamic model. The solution procedure of the RDM works, as follows. Given the evolutions over time of the energy demand Y e , of the technological parameters (in particular, q f ), and of the TC, the RDM is used to determine the trajectory of ϰ ( t ) and k ( t ) during the ET and the time that is needed to complete it (i.e., to reach ϰ = 1 ). Knowing the capital stock K ( t 1 ) , the capital growth rate g K ( t 1 ) is either calculated via Equation (30) (Exercise 1) or imposed (Exercise 2), and the capital stock K ( t ) is calculated via (5). Y e ( t ) is exogenously determined and k ( t ) can be deduced via Equation (26). The penetration rate of renewable energy ϰ ( t ) is then determined by inverting Equation (29). Finally, thanks to Equations (21)–(23), the values of the ES parameters ( v e ( t ) , q e ( t ) , and x e ( t ) ) are calculated. Once the previous variables have been determined, all of the other economic variables (sectoral production, consumption, GDP, etc.) can be calculated via the MM.
The first time periods of an example of ET in the AS are illustrated by the black dots shown in Figure 3, in the particular case where all exogenous technical parameters are constant (the envelopes of the AS and the CT are, therefore, not time-dependent). It is important to clearly distinguish the TC (dotted curve) from the trajectory representing the ET (the points). Mathematically, the TC is continuous while the trajectory is a discrete set. If the TC is postulated and, therefore, exogenous, the pace of the ET along the curve is endogenous and it is this pace that determines the ET actually observed.
Finally, it should be noted that there are two contradictory forces influencing the ET. On the one hand, the positive loop between k and g K noted above favours the growth of k and, thus, the ET. On the other hand, the ET results in an increase of x e and v e , which have a negative effect on growth. If an increase of k was to be accompanied by a too strong increase in these parameters, then g K could become negative, leading to further decreases of k and g K , which would block the transition. Therefore, for the ET to have a chance to continue, it must be done with enough available capital.

2.5. Calibration of the Macroeconomic Model

In [24], the variables of the MM are calibrated on the basis of global observations or estimates for a so-called reference year. A sensitivity analysis shows that the calibration is not very sensitive to variations in these parameters.
Subsequently, the units of good ( u g ) and energy ( u e ) of the theoretical model are translated into real units: respectively, 2010 US dollars (US$ 2010) and Joules (J).
Observed data
These are (i) IEA data for final energy consumption (https://www.iea.org/data-and-statistics, last access on 10 May 2021), and (ii) World Bank statistics for economic data: GDP and savings rate (https://data.worldbank.org: Gross savings (% of GDP) (s) and Gross added value at basic prices (GVA) (constant 2010 US$)) (GDP), last access on 10 May 2021).
Estimated data
Five model parameters have to be estimated in order to allow for the calibration of the whole model. These parameters are:
  • Sectoral capital depreciation rates. Based on the literature, the lifetime of capital in the ES is estimated at 25 years, and that in the FS at 20 years. This gives δ e = 1 T e = 0.04 and δ f = 1 T f = 0.05 .
  • The share of total capital devoted to the ES, as measured by the ratio μ = K e K . The calibration of the GEMBA model [16] provides the value μ =6.5%.
  • The share of the ES in the GDP, as measured by the ratio α = AV e GDP . Based on the OECD input-output tables [28], an estimate of 6% is found for the year 2015 (the last year for which the tables are currently available). This estimate is consistent with other estimates found in the literature [29,30,31].
  • The share of energy expenditure in the total household and state budget measured by η (Equation (19)). On the basis of the OECD input-output tables, an estimate of 4% is found for 2015.
Calibrated data
The observed and estimated data allow for the calibration of the model parameters and variables, for the chosen reference year (2018). The technical parameters of the two sectors and the EROI are listed in Table 1.
Knowing the technical parameters of the ES in 2018 (Table 1), the corresponding renewable penetration rate (according to the IEA statistics, wind and solar accounted for 2% of the total primary energy supply in 2018, i.e., ϰ ( 2018 ) = 0.02 ), and the evolution of the technical parameters of the renewable subsector thanks to the EM, the technical parameters of the non-renewable subsector (i.e., q n r e , x n r e , v n r e ) are found via Equations (21)–(23), see Table 2).
Figure 4 presents the evolution of the technical parameters of the ES as a function of the penetration rate ϰ , in the particular case where Y e is also set at its 2018 value (thus, for ϰ = 1 , Y r e = Y e ( 2018 ) ). ϰ evolves between 0 (all energy production is provided by non-renewable energies) and 1 (all of the energy production is provided by wind and solar technologies).
  • Figure 4a shows that the capital intensity v e increases considerably with the penetration rate ϰ : v e ( 1 ) / v e ( 0 ) > 45 , which is consistent with the common assumption that renewable energies are much more capital intensive than non-renewable energies.
  • Figure 4b shows that x e increases slightly with the ET: x e ( 1 ) / x e ( 0 ) = 1.3 . The NEIC requirements are not much higher for renewable energies than for non-renewable energies.
  • Because q n r e = 9.8 % and q r e = 8.7 % , the energy intensity q e slightly decreases with ϰ (Figure 4c).
Historical values of s and g K
Within the framework of the simulations, either the value of the savings rate s or the capital growth rate g K will be set at its historical level. To do this, the value averaged over the last five years was used. We find, for the period 2014–2018, for s, an average value of 25.05% and, for g K , an average value of 2.27%.

2.6. Future Evolution of the Parameters

Future changes in the technical parameters of the two sectors and energy demand Y e are based on the following assumptions:
  • The evolution of the technical parameters of the ES with the ET is described by Equations (21)–(23), where the technical parameters of the non-renewable energy subsector are assumed to be constant and equal to their calibrated value in 2018 (see Table 2). Assuming that the parameters of the non-renewable energy subsector will not change in the future can be considered to be optimistic. Given the maturity of the technologies, technological progress seems to be limited and not able to counterbalance the depletion of the resources. It is therefore likely that in the future there will be an increase in the energy and capital intensities of the non-renewable energy subsector.
  • Based on past values, the parameters x f and v f are assumed to be constant. Indeed, the MM calibration performed over the last 30 years shows that these parameters have not changed since 1990 (see Appendix E).
  • The changes in Y e and q f are based on the IEA’s "World Energy Outlook” 2020 report [25]. There are two scenarios: (i) the business-as-usual scenario (BAU), where changes in Y e and q f are extrapolated on the basis of the last 30 years (by 2040, Y e increases by 45% and q f decreases by 15%) and (ii) the sustainable development scenario (SD), which assumes a constant energy demand Y e and a rapidly decreasing energy intensity q f in the coming decades (by 2040, q f would decrease by 40%, which corresponds to a rate of decline never observed in the past) (see Appendix F for further details). In the long term, both of the scenarios have the same technical progress potential in the FS ( q f is asymptotically reduced by 50%), the difference between the two scenarios being the pace at which this limit is reached. A third scenario based on the IEA’s New Policies scenario was also tested, but did not lead to fundamentally different results from those of the BAU scenario.

3. Results

Thereafter, we consider TCs in the form of weighted averages of the AS envelopes, i.e., curves of the form:
Φ ( ϰ ) = [ 1 θ ] k ̲ ( ϰ ) + θ k ¯ ( ϰ )
where (i) θ is the weighting parameter and (ii) k ̲ ( ϰ ) and k ¯ ( ϰ ) are the envelopes of the AS defined by Equations (A9) (Exercise 1) and (A13) (Exercise 2). As the energy intensity of the FS q f decreases over time (see Section 2.6), the envelopes of the AS and the TC are both functions of t. θ is exogenous, constant, or a function of t, chosen so as to always respect inequalities 0 θ 1 . The formulation from Equation (33) then ensures that the TC remains within the AS in each time period.

3.1. Exercise 1—s and η Given

3.1.1. Impact of θ

The purpose of this subsection is to scan the entire AS and measure how the results change when the TC move through this space. A simple way to do this is to consider all of the TCs with a constant weighting parameter θ .
The graphs shown in Figure 5 provide, for the two scenarios: (i) ϰ ( 2100 ) , the renewable penetration rate in 2100, (ii) ϰ l , the asymptotical value of ϰ , (iii) if it is completed in a finite time, the year of the end of the ET, T f , and (iv) if any, the residual capital growth rate at the end of the transition, g K l . When the transition is completed in a finite time ϰ l = 1, if not, ϰ l is obtained when ϰ no longer varies, because g K tends towards 0. There is a residual capital growth rate when the transition is completed in a finite time. However, this residue is very low and never exceeds 0.06%.
Business-as-usual scenario (red curves)
  • Figure 5a gives the value of the penetration rate ϰ in 2100. The curve is slightly increasing, but ϰ ( 2100 ) is approaching 30% at best.
  • Figure 5b shows that the ET is never completed in a finite time. Nevertheless, we see that ϰ l is approaching 1 for θ = 1 , and ϰ l is above 65% for all values of θ . The non-completion of the ET is due to the fact that the growth rate of the capitak stock g K decreases down to zero.
  • As mentioned above, the residual growth rate of the capitak stock g K l is 0 for all the simulations. The economy tends towards a stationary state with no growth.
Sustainable development scenario (blue curves)
  • Figure 5a shows that the renewable penetration rate ϰ in 2100 is significantly higher for SD than for BAU.
  • Figure 5b shows that, if θ < 0.86 , no scenario allows for a complete ET in a finite time.
  • Figure 5c shows that for θ 0.86 the duration of the ET decreases with θ (from 248 years for θ = 0.86 to 131 years for θ = 1 ). The minimum value of T f is reached for θ = 1 (a complete transition by 2149).
  • Figure 5d shows that the residual capital growth rate g k l is 0 for θ < 0.86 and sightly increasing for θ 0.86 . At best, it reaches 0.056%.
The previous results only concern TCs with a constant weighting rate. However, numerically, the following conjecture is verified:
Consider a TC Φ v defined by (33), whose weighting parameter is variable and described by the sequence { θ t : t 0 } . Let θ ¯ ( 1 ) be the maximum of this sequence. Let Φ θ ¯ be the TC whose weighting parameter is constant and equal to θ ¯ . Subsequently, T f ( Φ θ ) > T f ( Φ θ ¯ ) .
The interest of this result is to provide a minimum bound for the duration of an ET based on a TC with variable weighting parameter, from the values of T f calculated for TCs with constant weighting parameter (the values given by Figure 5c).
As a corollary of this conjecture, the duration of any ET will be greater than or equal to T f ( Φ 1 ) , where Φ 1 designates the TC characterised by θ = 1 (i.e., greater than 104 years for SD). As a reminder, simulations with θ = 1 would correspond to an economy where the net added value of the FS is zero. The price of energy is so high that all of the added value of the economy comes from the ES.

3.1.2. Impact of the Technical Progress Rate

Whether for BAU or SD, the technical progress potential is assumed to be very high, since the energy intensity of the FS is asymptotically divided by 2. The difference between the two scenarios lies in the potential already realised in 2040 (30% for BAU, 80% for SD).
Further increasing the potential realised in 2040 under SD does not fundamentally change the results of the previous subsection. For the record, the evolution of q f is governed by Equation (A15). If the rate of change r q f is increased sharply to the point that the potential of technical progress is realised in a few years, the end of the ET (as measured by T f ) is slightly advanced, but remains beyond 2100.
Consequently, it appears from the previous simulations that, unless assuming unreasonable progress in terms of energy efficiency, it is not possible to achieve a complete ET before the end of the century, even when assuming energy demand control (as required by the SD scenario).

3.2. Exercise 2— η and g K Fixed

Previous simulations have shown how difficult it would be to achieve a complete ET before the end of the century if consumption and savings choices remain unchanged. The aim of this section is to study an alternative strategy, aiming at supporting capital growth through an increase in the savings rate.

3.2.1. Impact of θ

As in Section 3.1.1, the following simulations are based on a constant weighting parameter θ . Table 3 shows the values of T f , s, and ϰ (in 2050, 2100, and T f ) obtained for Exercise 2 for θ = 0 and θ = 1 . The changes in the variables between these extreme values are monotonous. The capital growth rate g K is fixed at its average value over the last five years: 2.27%.
Business-as-usual scenario
The table shows that:
  • The ET is completed in a finite time for θ = 0 and θ = 1 and, thus, for all values of θ [ 0 , 1 ] . T f is slightly increasing with θ , but it varies little ( T f goes from 2129 for θ = 0 to 2133 for θ = 1 ). In Exercise 1, an increase of θ stimulated capital growth and, thereby, accelerated the ET. In Exercise 2, this effect is not present, since g K is fixed.
  • The ET is far from being completed by the end of the century. ϰ ( 2100 ) varies between 46% for θ = 0 and 17% for θ = 1 .
  • As for T f , the observed values of s at the end of the ET evolve little with θ (from 44% for θ = 0 to 41% for θ = 1).
In summary, this analysis reveals that all of the TCs included in the AS lead to a complete ET in about 115 years, with a significant increase in the savings rate.
Sustainable development scenario
Previous results have shown that, under the BAU scenario, maintaining the capital growth rate at its historical level through a gradual increase in the savings rate would allow a complete ET in about 115 years. The end of the transition does not occur until the next century.
However, in the context of Exercise 1, we have seen that the SD scenario made it possible to significantly reduce the duration of the ET when compared to BAU (see Section 3.1.1). Intuition then suggests that the combination of the hypotheses of Exercise 2 and the SD scenario should enable a complete ET before the end of the century.
In order to verify this, the previous analysis made for BAU is repeated for SD. Table 3 shows the values obtained. We observe that:
  • The previous intuition is confirmed, since the completion of the ET is obtained after a period of 67 to 72 years, i.e., a value of T f around 2090.
  • Maintaining the capital growth rate g K during the ET requires a significant increase in the savings rate s. However, this increase is less than for BAU, as s reaches between 37 and 40% at the end of the ET, depending on the value of θ .
From the above, it appears that the SD scenario, combined with the support of capital growth, allows for completing the ET before 2100.

3.2.2. SD: g K Required for T f = 2070

With the previous simulations, in the best case scenario, we are still beyond certain target dates retained in the literature in the fight against climate change. Indeed, according to a recent IPCC report [32], limiting the global temperature increase to 2 °C requires the achievement of carbon neutrality by 2070 (and a reduction of 25% in emissions by 2030 as compared to their 2010 level). The latter key date (2070) is also chosen by the IEA to achieve global carbon neutrality as part of the sustainable development scenario [33].
The previous point showed that maintaining the capital growth rate at its historical level (i.e., 2.27%, its average value over the last five years period) within the framework of SD made it possible to complete the ET in about seventy years. Intuition suggests that increasing the capital growth rate should lead to an even shorter transition.
The following simulation results are based on the objective of an end of the ET in 2070, and they aim at determining the growth rate g K necessary to reach this goal, as a function of θ . Table 4 shows that:
  • g K increases slightly with θ , from 2.91% to 3.12%.
  • The savings rate that was observed at the end of the ET decreases slightly with θ , while still remaining above 40%.
    Note: similar simulations can be carried out within the framework of BAU. To complete the transition by 2070, a capital growth rate of around 4.1% is required, depending on θ . To maintain this growth rate, a savings rate of more than 55% is then necessary (s at the end of the ET decreases from 58.6% for θ = 0 to 56.2% for θ = 1 ).

3.2.3. Temporal Analysis

Previous results show that an ET completed before the end of the century requires a significant increase in the savings rate, highlighting that the transition has a profound effect on the economy. Using the MM, the following analysis aims to shed light on the behaviour of economic variables during the transition of the ES.
The observations below correspond to an ET that was completed in 2070 with a capital growth rate g K equal to 3%, under the SD scenario. The TC is chosen to ensure a gradual transition without delay, i.e., a constant weighting rate θ = 0.418 . TCs that are associated with lower θ values lead to a “strong” start of the transition followed by a slowdown. TCs that are associated with higher θ values lead to a delay before the start of the ET.
Figure 6 and Figure 7 show that:
  • Maintaining the growth rate g K during the ET requires a gradual increase in the savings rate s (Figure 6b). This increase is significant since s increases from 25 to about 40% at the end of the ET. This figure is similar to those of Režný and Bureš [19] and of Fagnart et al. [17].
  • The transition is characterised by a strong increase in the share of capital devoted to the ES ( μ = K e K ), which exceeds 40% at the end of the trajectory, as shown in Figure 6c. This result is similar to those of Dale [16] and Fagnart et al. [17].
  • Unsurprisingly, the ET results in a significant increase in the real price of energy p (Figure 6d). Despite the decrease of q f due to technical progress, the growth of the economy results in an increased demand for energy, which pushes up its price. At the end of the ET, the price p reaches six times its initial value. This results contradicts the one of Jacobson [4], claiming that the price of energy will not increase with the ET.
  • The rise in the price of energy translates into a significant increase in the share of the ES gross added value in GDP ( α = AV e GDP ), as shown in Figure 6e. The increase in the share of the ES net added value in NDP (as measured by α n = NV e GDP ) is more moderate.
  • Figure 6f illustrates the evolution of the FS production ( Y f ), as well as how it is distributed between final consumption C f , investment I, and intermediate consumption deliveries to the two sectors X ( = X e + X f ) . While final output remains growing during the transition, there is a decline in the growth of consumption of goods C f , which, like C, is almost constant at the end of the ET. In other words, because of the increase in the savings rate, the allocation of Y f changes to the detriment of C f and to the benefit of I (the share of the intermediate consumption deliveries remaining almost constant).
  • The increase in s translates into a slowdown in aggregate consumption C with respect to GDP. C becomes quasi-stationary at the end of the ET, while the economy is still growing (Figure 7a).
  • Figure 7b illustrates how the ES production Y e is split between self-consumption E e , deliveries to the FS E f , and final energy consumption C e . It can be observed that the share of E f in Y e increases significantly to the point of taking up most of the energy production at the end of the ET.
  • The ET is accompanied by a strong increase in v e , the ES capital intensity. As the share of capital in the economy that is devoted to this sector continues to grow, the capital intensity of the economy, which is defined as v = K GDP , increases (Figure 7c).
  • In Appendix G, we show that the GDP growth is decomposed into the growth of the capital stock and the decline in the capital intensity of the economy (defined as g v = Δ v v ): g g K g v . Given the previous observation (v is increasing), the GDP growth rate declines despite the constancy of g K , as confirmed by Figure 7d.
  • Figure 7e illustrates the decrease in the ratios GDP / Y T (where Y T = Y f + p Y e denotes the total output of the economy, including all intermediary consumptions) and NDP / GDP during the ET. The decrease in the first ratio indicates that the ET results in a faster growth of the intermediary consumptions as compared to GDP (the end-use consumption). The decrease of the second ratio is linked to the increase in the share of capital consumption δ K in GDP, which is induced by the increase in capital intensity v (the increase of v outweighs the decrease of δ , the latter being due to the fact that the life span of capital is higher in the ES than in the FS). The share of useful production (i.e., for consumption or growth) in total production tends to decrease with the ET.
  • Figure 7f shows that the societal point-of-use EROI only slightly decreases during the ET. Several opposing forces are at work: (i) the ET results in an increase of v e and x e , which lowers the EROI, (ii) the decrease of q e increases the EROI, and (iii) the technical progress behind the efficiency gains in the FS increases the EROI.
In conclusion, a strategy aiming at supporting the growth of capital through an increase in the savings rate makes it possible to reach the desired duration for the ET. The price to pay is a slowdown in the growth of total consumption, which ends up stagnating at the end of the transition.

4. Discussion

The simulation results confirm the main findings of the models of the critical literature that are cited in the introduction [9,10,12,13], in particular those presenting the evolution of certain macroeconomic variables in the context of the energy transition [17,19]. It is observed that, if completed in time, the energy transition would be accompanied by a considerable increase in the economy’s savings rate, exceeding 40% in all of the scenarios, as well as in the share of capital that is devoted to the energy sector, rising from less than 10% to more than 40%. At unchanged consumption behaviours (i.e., Exercise 1), the energy transition is almost never completed (except for the SD scenario and a few values of θ ) because the economy tends towards a stationary state.
In the model of Fagnart et al. [17], the transition is fully completed after a certain period of time, but it is accompanied (with a few exceptions) by a decrease in GDP. This is also the case with the MEDEAS model, where the economy decreases from 2050 onwards in most of the scenarios envisaged [10]. When compared to these general equilibrium models, we are in a partial equilibrium context where the energy demand does not depend on the future evolution of the economy, which might explain why we do not observe a decrease in GDP.
It is important to stress that the obtained results can be considered to be optimistic. Indeed, the underlying assumptions can be considered quite "conservative”; in particular:
  • At the beginning of the simulation, the technical progress potential lowering the FS energy intensity is assumed to be significant, since it is (asymptotically) divided by 2. For the BAU scenario, 55% of the gains are obtained in 50 years, and 80% in 100 years. For the SD scenario, most of the gains (97%) are assumed to be obtained before 50 years.
  • Production factor intensities ( q n r e , x n r e , and v n r e ) of the non-renewable energy subsector are assumed to be constant. However, given the maturity of technologies, technological progress does not seem likely to continue to offset the gradual depletion of resources and the fact that their extraction is becoming increasingly costly in terms of energy and capital. Therefore, a future increase in these intensities seems to be plausible. The same goes for the energy intensity of the renewable subsector ( q r e ), which is expected to increase with renewable penetration, as modeled in MEDEAS [27].
  • Capital is assumed to be perfectly malleable, in the sense that it can move between the non-energy and energy sectors on the one hand, and between the two subsectors of the ES on the other.
  • The site selection criterion is a criterion of cost-effectiveness, whether this is measured in energy or economic terms, as indicated in Section 2.3. However, many other factors come into play when setting up a new renewable energy park. This means that, in practice, less profitable sites could be selected before more profitable ones.
Assumptions 1 and 2 explain a large part of the relatively small decrease in EROI during the ET when compared to other models, while, at the same time, the growth rate of the economy declines. Equations (15) and (20) show that the EROI would decrease more if (i) q f decreased less as a result of lower technical progress potential, (ii) v e and x e increased even more as a result of the increase in x n r e and v n r e , and (iii) q e increased as a result of the increase of q r e .
Hypothesis 3 allows for reducing a capital stock faster than what its “natural” depreciation rate would allow. It is not applicable if (i) the stock is increasing or (ii) it decreases less quickly than its depreciation rate. Within the framework of the simulations studied, cause (i) is verified for the FS and the renewable energy subsector because their respective stocks are always growing. Cause (ii) applies to the non-renewable energy subsector during much of the ET, but no longer when its end is approaching. Hypothesis 3 then becomes operative and everything happens “as if” thermal power plants were transformed into solar panels or wind turbines, which is, of course, unrealistic.
Removing any of the hypotheses 1 to 4 would complicate the ET and result in its lengthening.

5. Conclusions

This paper develops a modeling scheme articulating a macroeconomic input-output model with two sectors (energy and non-energy) and an energy mode that is able to calculate the maximum potentials of solar and wind energy in order to study the feasibility and pace of the energy transition at the global level. This methodology allows to take two effects into account: a crowding out effect of the non-energy sector in terms of capital allocation, and a location effect due to the declining “quality” of the sites that can receive the additional wind and solar installations.
The main findings obtained thanks to this paper are as follows. At an unchanged savings rate, the transition to a 100% renewable energy system seems unachievable before the end of the century. Controlling energy demand at its current level will drastically reduce the duration of the transition, but it is not enough. The reason lies under the higher capital needs of the renewable energy technologies, resulting in the crowding out of the FS in terms of capital allocation. This effect weighs ever more negatively on growth, slowing down and then blocking the transition in most of the scenarios considered, as the economy tends towards a stationary state.
The strategy aiming at supporting capital growth through an increase in the savings rate allows to counteract the crowding out effect and to reduce significantly the duration of the ET. However, if the demand for energy continues to grow significantly (an increase of 150% by 2100 is expected in the business-as-usual scenario), the transition is only completed after the end of the century. The energy demand should be kept at its current level and the growth rate of capital should be maintained at a value above its historical level in order to allow for a complete transition by 2070. However, because of the significant increase in the savings rate, total consumption ends up stagnating at the end of the transition (while the economy is still growing). On top of that, despite the maintenance of the capital growth rate, the GDP growth rate decreases as a result of the increase in the capital intensity of the economy (as measured by the capital/GDP ratio).
The technical progress affecting the energy intensity of the FS proves to be a fundamental determinant of whether or not the ET is feasible in a finite time. Within the framework of the strategy of the previous point, this intensity must decrease by at least 50% by the end of the ET in 2070.
These results confirm those of the critical literature reviewed in the introduction, in particular concerning the impact of the crowding out effect on growth and the need to control energy demand in order to be able to complete the ET in a sufficiently short period of time.
In order to increase the relevance of the model and to refine the results, several developments are possible. Modelling the non-renewable energy subsector in the same way as what was done for the renewable energy subsector and removing the hypothesis of perfect malleability of capital would contribute to greater realism. Endogenising (i) consumption and investment behaviours and (ii) energy demand would enrich the interactions between the economy and the ET. Taking into account the labour factor next to the capital factor while requiring that their respective remunerations be positive would make it possible to refine the space of allowable trajectories.
Finally it would be interesting to add the cumulative CO 2 emissions to the simulations. Indeed, although we manage, under certain conditions, to complete the transition by 2070, there is no guarantee that the cumulative quantities of CO 2 emitted do not exceed the critical threshold necessary to limit the increase in global temperature.

Author Contributions

All authors have contributed to the different stages. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations and notations are used in this manuscript:
ETenergy transition
EROIEnergy Return On Investment
MMmacroeconomic model
EMenergy model
ESenergy sector
FSnon-energy sector
EICenergy intermediate consumption
NEICnon-energy intermediate consumption
u e units of energy
u g units of goods produced by the FS
AV i sectoral gross added values
NV i sectoral net added values
GDPgross domestic product
NDPnet domestic product
ASallowable space
TCtransition curve
RDMreduced dynamic model
BAUbusiness-as-usual scenario
SDsustainable development scenario
q i sectoral energy intermediate consumption (EIC) intensities
v i sectoral capital intensities
x i sectoral non-energy intermediate consumption (NEIC) intensities

Appendix A. The Societal EROI

Point-of-use EROI is defined as follows:
EROI = Energy delivered Direct + Indirect energy inputs to deliver the energy
In order to produce an amount of energy Y e , the ES consumes E e as EIC (i.e., direct energy consumption) (in u e ), X e as NEIC, and uses a stock of capital K e .
Since the EROI only accounts for energy flows, it is necessary to estimate the embodied energy included in the indirect energy inputs, i.e., the intermediate consumption flow X e and capital stock K e . As annual flows are considered in the model, we take into account the capital depreciation in the denominator of the EROI. The capital flow in the denominator of the EROI is estimated as the energy content of δ e K e , i.e., all the energy that was required to produce, install, maintain and dismantle the fraction of the ES capital stock consumed during the period under consideration. This estimate is therefore based on the energy content of the share of the investment (at time t) used to replace the capital that has become obsolete (at time t). An alternative approach would have been to estimate the embodied energy in the capital stock based on the energy content of past investments. Since q f decreases over time, this approach would lead to a higher estimate of the embodied energy contained in the flow measuring capital depreciation.
Let E X e and E δ e K e (both in u e ) be the energy embodied in the flows measuring NEIC X e and capital depreciation δ e K e . Since (i) the production of the FS is assumed to be homogeneous and (ii) it takes q f unit of energy to produce one unit of good, this embodied energy flow can be estimated by:
E X e + E δ e K e = q f ( X e + δ e K e ) = q f ( x e + δ e v e ) Y e = q f w e Y e
given Equations (1) and (15). The energy delivered (i.e., the amount that reaches the FS plus the amount used for final consumption) is defined as the gross energy produced Y e , minus the direct energy inputs E e . The EROI then becomes:
EROI = Y e E e E e + E X e + E δ e K e = Y e ( 1 q e ) E e + q f ( x e + δ e v e ) Y e = 1 q e q e + q f w e
One finds Equation (20).

Appendix B. Site Selection Criteria

In the EM, it is considered that the sites with the best EROI will be chosen first. For this approach to be economically consistent, this criterion must be the same as the criterion of selecting the sites with the highest profits.
For the renewable technologies of interest here (wind and solar), (i) the energy intensity of an installation is the same for all sites considered, (ii) the capital required K e and the maintenance costs included in X e are fixed costs, which are not site-dependent.
In view of the above, the profits related to sites a and b for a given technology are written:
Π a = p ( Y e , a E e , a ) X e δ e K e
Π b = p ( Y e , b E e , b ) X e δ e K e
Profits per unit of energy produced are then written:
π a = p ( 1 q e , a ) X e + δ e K e Y e , a = p ( 1 q e , a ) w e , a
π b = p ( 1 q e , b ) X e + δ e K e Y e , b = p ( 1 q e , b ) w e , b
If Π a > Π b , then necessarily Y e , a > Y e , b , which implies w e , a < w e , b (while q e , a = q e , b ).
Now, the point-of-use EROIs of the two sites are written respectively (from Equation (20)):
EROI a = 1 q e , a q e , a + q f w e , a
EROI b = 1 q e , b q e , b + q f w e , b
It follows that Π a > Π b EROI a > EROI b .

Appendix C. Interval for p

Criterion defined in Equation (24) requires that the net sectoral added values are both positive. With Equations (9) and (13), criterion NV f > 0 translates into the following inequality:
1 p q f x f δ f v f = 1 p q f w f > 0
With Equations (10) and (14), criterion NV e > 0 translates into the following inequality:
p ( 1 q e ) x e δ e v e = p ( 1 q e ) w e > 0
which implies 1 w f > p q f and p ( 1 q e ) > w e , hence the following allowable interval for p:
I p w e 1 q e < p < 1 w f q f
I p can be transformed into an interval for k , but this will depend on assumptions made about the parameters, in particular the behavioural parameters s (or g K ) and η .

Appendix D. Calculation of the Allowable Spaces

Thereafter, the allowable range for the ratio k = K Y e is established for the 2 exercises, as well as the equation that links ϰ and k with g K (Exercise 1) or s (Exercise 2).

Appendix D.1. Preliminaries

For the developments that will follow, let us define D f = C f + I , the final demand to the FS.
Equations (6) and (7) give the following relationships for final demands D f and C e :
D f = Y f X e X f = [ 1 x f ] Y f x e Y e
C e = Y e E e E f = [ 1 q e ] Y e q f Y f
By the expenditure approach of GDP, we have:
GDP = D f + p C e = D f + η C = D f + η [ 1 s ] GDP
Hence the following relationship between GDP, D f , η and s (useful for Exercise 1):
GDP = D f 1 η [ 1 s ]
And, as GDP = C + I , the following relationship between C, D f , I and η (useful for Exercise 2):
C + I = D f + η C C = D f I 1 η

Appendix D.2. Exercise 1—k Interval

With Equations (18), (19), (A4)–(A6), p is expressed as a function of the technical parameters and the outputs Y i :
p = η C C e = η [ 1 s ] GDP C e = η [ 1 s ] 1 η [ 1 s ] D f C e = η [ 1 s ] 1 η ] 1 s ] [ 1 x f ] Y f x e Y e [ 1 q e ] Y e q f Y f
Let y = Y f Y e , then the last equality becomes:
p = η [ 1 s ] 1 η [ 1 s ] [ 1 x f ] y x e 1 q e q f y
As a result, and with η = η [ 1 s ] 1 η [ 1 s ] , A = w e 1 q e , and B = 1 w f q f , Equation (A3) becomes:
A < η [ 1 x f ] y x e 1 q e q f y < B
A [ 1 q e q f y ] < η [ [ 1 x f ] y x e ] < B [ 1 q e q f y ]
A [ 1 q e ] + η x e < [ A q f + η [ 1 x f ] ] y [ B q f + η [ 1 x f ] ] y < B [ 1 q e ] + η x e
A [ 1 q e ] + η x e A q f + η [ 1 x f ] < y < B [ 1 q e ] + η x e B q f + η [ 1 x f ]
I y y ̲ = w e + η x e w e q f 1 q e + η [ 1 x f ] < y < [ 1 w f ] [ 1 q e ] q f + η x e 1 w f + η [ 1 x f ] = y ¯
To I y corresponds an interval for k. Indeed, k and y are linked by:
k = K Y e = K f + K e Y e = v f Y f + v e Y e Y e = v y y + v e
Therefore, we have:
k ̲ ( ϰ ) = v y w e ( ϰ ) + η x e ( ϰ ) w e ( ϰ ) q f 1 q e ( ϰ ) + η [ 1 x f ] + v e ( ϰ ) k ¯ ( ϰ ) = v y [ 1 w f ] [ 1 q e ( ϰ ) ] q f + η x e ( ϰ ) 1 w f + η [ 1 x f ] + v e ( ϰ )
These bounds are increasing functions of ϰ (since x e ( ϰ ) , v e ( ϰ ) and thus w e ( ϰ ) > 0 , and q e ( ϰ ) < 0 ), and decreasing functions of q f .
Indeed, the EM results show that solar and wind power are more intensive in goods produced by the FS than the technologies they replace. It should also be noted that the lower the energy intensity of the FS q f , the higher the envelopes are in the { ϰ , k } space.

Appendix D.3. Exercise 1—Capital Growth Rate

Given Equations (3) and (17), the capital growth rate is calculated as follows:
g K = Δ K K = s GDP K δ
Yet Equation (A6) gives:
s GDP K = s 1 η [ 1 s ] D f K
Substituting D f by Equation (A4), one finds:
s GDP K = s 1 η [ 1 s ] [ 1 x f ] Y f x e Y e K = s 1 η [ 1 s ] [ 1 x f ] K K e v f x e Y e K = s 1 η [ 1 s ] [ 1 x f ] k v e v f x e k = 1 v f s 1 η [ 1 s ] [ [ 1 x f ] [ 1 v e k ] x e v f k ]
Given Equation (4) and (A11), Equation (A10) rewrites:
g K = 1 v f s 1 η [ 1 s ] [ [ 1 x f ] [ 1 v e k ] x e v f k ] δ ( k ) = 1 v f s 1 η [ 1 s ] [ 1 x f [ 1 x f ] v e + x e v f k ] δ f + [ δ f δ e ] v e k
The first term is an increasing function of k, and δ ( k ) is an increasing function of k if δ f δ e , and negative otherwise. It follows that g K is an increasing function of k if δ f δ e is not too high.

Appendix D.4. Exercise 2—k Interval

Given Equations (A7) and (19) (which gives p = η C / C e ):
p = η 1 η D f I C e
Equation (A10) gives an expression for I ( I = [ g K + δ ( k ) ] K ). Combined with Equations (A7) and (A5) the previous equation becomes:
p = η 1 η [ 1 x f ] Y f x e Y e [ g K + δ ( k ) ] K [ 1 q e ] Y e q f Y f = η 1 η [ 1 x f ] K v e Y e v f x e Y e [ g K + δ ( k ) ] K [ 1 q e ] Y e q f K v e Y e v f = η 1 η [ 1 x f ] k v e v f x e [ g K + δ f [ δ f δ e ] v e k ] k 1 q e q f k v e v f = η 1 η [ 1 x f v f [ g K + δ f ] ] k + [ δ f δ e ] v e [ 1 x f ] v e v f x e 1 q e + q f v e v f q f v f k
With notations η = η 1 η , A = 1 η w e 1 q e , B = 1 η 1 w f q f , Equation (A3) becomes:
A < [ 1 x f v f [ g K + δ f ] ] k + [ δ f δ e ] v e [ 1 x f ] v e v f x e 1 q e + q f v e v f q f v f k < B
Both the numerator and the denominator are positive, so the inequality rewrites as follows:
A ( 1 q e + q f v e v f q f v f k ) < ( 1 x f v f ( g K + δ f ) ) k + ( δ f δ e ) v e ( 1 x f ) v e v f x e < B ( 1 q e + q f v e v f q f v f k )
Which gives two inequalities:
A ( 1 q e + q f v e v f ) ( δ f δ e ) v e + ( 1 x f ) v e v f + x e < ( A q f v f + 1 x f v f ( g K + δ f ) ) k
( B q f v f + 1 x f v f ( g K + δ f ) ) k < B ( 1 q e + q f v e v f ) ( δ f δ e ) v e + ( 1 x f ) v e v f + x e
If δ f < δ e , or if δ f δ e > 0 is not too high, the left-hand side of the first inequality is positive, which implies that the other members of the 2 inequalities are also positive (because A < B ). Hence the following interval is found for k:
A ( 1 q e + q f v e v f ) ( δ f δ e ) v e + ( 1 x f ) v e v f + x e q f v f + 1 x f v f ( g K + δ f ) < k < B ( 1 q e + q f v e v f ) ( δ f δ e ) v e + ( 1 x f ) v e v f + x e B q f v f + 1 x f v f ( g K + δ f )
Or equivalently:
k ̲ ( ϰ ) = 1 η w e ( ϰ ) 1 q e ( ϰ ) [ 1 q e ( ϰ ) + q f v e ( ϰ ) v f ] + [ 1 x f v f [ δ f δ e ] ] v e ( ϰ ) + x e ( ϰ ) 1 η w e ( ϰ ) 1 q e ( ϰ ) q f v f + 1 x f v f [ g K + δ f ] k ¯ ( ϰ ) = 1 η 1 w f q f [ 1 q e ( ϰ ) + q f v e ( ϰ ) v f ] + [ 1 x f v f [ δ f δ e ] ] v e ( ϰ ) + x e ( ϰ ) 1 η 1 w f v f + 1 x f v f [ g K + δ f ]

Appendix D.5. Exercise 2—Savings Rate

With Equation (A7) we write:
GDP = C + I = D f I 1 η + I = D f η I 1 η
With the definition of s and the previous equation we have:
s = I GDP = [ 1 η ] I D f η I
Yet I = [ g K + δ ( k ) ] K , so:
s = [ 1 η ] [ g K + δ ( k ) ] D f K η [ g K + δ ( k ) ]
In addition, with Equation (A4) we have:
D f K = [ 1 x f ] Y f x e Y e K = [ 1 x f ] K v e Y e v f x e Y e K = [ 1 x f ] k v e v f x e k = 1 x f v f [ 1 v e k ] x e k = 1 x f v f [ 1 x f v f v e + x e ] 1 k
As a result we find an expression for s as a function of the technical parameters, g K , η , k (and ϰ ):
s = [ 1 η ] [ g K + δ f [ δ f δ e ] v e k ] 1 x f v f η [ g K + δ f [ δ f δ e ] v e k ] [ 1 x f v f v e + x e ] 1 k

Appendix E. Historical Evolution of v, xf and vf

Starting from the estimation of the total capital stock for the reference year ( K ( 2018 ) ), and with historical data regarding the total GDP and savings rate s, we find capital stock estimates from the time period 1990–2018.
K ( t + 1 ) = K ( t ) δ K ( t ) + I ( t ) ) = K ( t ) ( 1 δ ) + s ( t ) G D P ( t ) K ( t ) = K ( t + 1 ) s ( t ) G D P ( t ) 1 δ
With theses estimates for K ( t ) and the observed data for the GDP and savings rate the historical evolution of v is estimated. Then thanks the calibration procedure developed by Dupont et al. [24] and the yearly estimations of K and v, the historical evolution of the technical parameters of the FS can be retraced. These are shown in Figure A1. While q f has decreased by an average of 0.8% per year over the last 30 years, v, v f and x f have not seen any significant change.
Figure A1. Historical evolution of: (a): v ( t ) the capital intensity of the economy (i.e., the capital/GDP ratio), (b): q f ( t ) , the energy intensity of the FS, (c): x f ( t ) , the NEIC intensity of the FS and (d): v f ( t ) , the capital intensity of the FS.
Figure A1. Historical evolution of: (a): v ( t ) the capital intensity of the economy (i.e., the capital/GDP ratio), (b): q f ( t ) , the energy intensity of the FS, (c): x f ( t ) , the NEIC intensity of the FS and (d): v f ( t ) , the capital intensity of the FS.
Sustainability 13 05345 g0a1

Appendix F. Future Evolution of Ye and qf

We consider 2 scenarios for the future evolution of the energy demand and the energy intensity of the FS. The first scenario is the business as usual scenario (BAU), which is an extrapolation of past trends of Y e and q f (see Table A1). The second scenario is based on the sustainable development (SD) scenario of the IEA “World Energy Outlook” report ([25]). It represents an accelerated energy transition in order to achieve sustainable development goals in terms of energy access (universal access to modern energy by 2030), climate change (limiting the average temperature increase to below 2 °C) and air quality (significant reduction in air pollution-related deaths).
For BAU, extrapolating the annual growth rates observed over the period 1990–2018 (see Table A1), we obtain for 2040: Y e ( 2040 ) = 1.45 Y e ( 2017 ) , and q f ( 2040 ) = 0.85 q f ( 2017 ) .
Table A1. Historical annual growth rates of Y e and q f .
Table A1. Historical annual growth rates of Y e and q f .
1990–20181990–19992000–20092010–2018
Y e [%]1.631.022.051.16
q f [%]−0.58−1.850.56−0.72
The IEA SD scenario assumes a energy consumption for 2040 similar to today, and a higher annual rate of decline in the energy intensity of the economy than has been observed in the past. The IEA considers an annual decline rate for the primary energy intensity of GDP of 3.5%. Since only final energy is considered here, the rate of decrease of q f will be lower. If we consider an annual decrease rate of 2.25% (a rate never observed over the 1990–2018 period), we find q f ( 2040 ) = 0.6 q f ( 2017 ) .
In the long term, asymptotic values are given for Y e and q f , which are Y e l and q f l . For BAU, Y e l = 3 Y e ( 0 ) , while for SD, Y e l = Y e ( 0 ) . For both scenarios, q f l = 0.5 q f ( 0 ) .
The evolution of q f is estimated as follows:
q f ( t ) = q f l + ( q f ( 0 ) q f l ) ( 1 r q f ) t
which assumes that the difference between the initial value of q f and its limit value ( q f l ) decreases exponentially. Knowing q f in 2040, we can estimate r q f , the rate of change of the gap from q f ( t ) to its stationary value q f l . In the BAU framework, the evolution of Y e is modelled according to a logistic curve:
Y e ( t ) = Y e ( t 1 ) ( 1 + r Y e ( 1 Y e ( t 1 ) Y e l ) )
The rate r Y e is calibrated iteratively, starting from the value of Y e in 2040.
The values of r q f and r Y e for the two scenarios are shown in Table A2. The corresponding temporal evolutions of Y e and q f are shown in Figure A2.
Table A2. Calculated values of r q f and r Y e for both scenarios.
Table A2. Calculated values of r q f and r Y e for both scenarios.
r Y e [%] r q f [%]
BAU2.721.6
SD07
Figure A2. Future evolution of Ye and qf for the 2 scenarios: (a): Ye(t), (b): qf(t), (c): qf(t)/qf0.
Figure A2. Future evolution of Ye and qf for the 2 scenarios: (a): Ye(t), (b): qf(t), (c): qf(t)/qf0.
Sustainability 13 05345 g0a2aSustainability 13 05345 g0a2b

Appendix G. Link between v, g and gk

Using Equations (3), (17) and the definiton of the capital intensity of the economy ( v = K GDP ) we write:
Δ K = s GDP δ K = ( s v δ ) K
Hence the capital growth rate g K is given by:
g K = Δ K K = s v δ
Be g the GDP growth rate: g = Δ GDP GDP . One can write at first order g v = Δ v v = Δ K K Δ GDP GDP , which implies:
g = g K g v = s v δ g v
It shows that GDP growth is decomposed into the growth of the capital stock and the decline in the capital intensity of the economy.

References

  1. Budischak, C.; Sewell, D.; Thomson, H.; Mach, L.; Veron, D.E. Cost-minimized combinations of wind power, solar power and electrochemical storage, powering the grid up to 99.9% of the time. J. Power Sources 2013, 225, 60–74. [Google Scholar] [CrossRef] [Green Version]
  2. Ellabban, O.; Abu-rub, H.; Blaabjerg, F. Renewable energy resources: Current status, future prospects and their enabling technology. Renew. Sustain. Energy Rev. 2014, 39, 748–764. [Google Scholar] [CrossRef]
  3. Jacobson, M.Z.; Delucchi, M.A.; Bazouin, G.; Bauer, Z.A.F.; Heavey, C.C.; Fisher, E.; Morris, S.B.; Piekutowski, D.J.Y.; Vencill, T.A.; Yeskoo, T.W. Environmental Science the 50 United States. Energy Environ. Sci. 2015, 8, 2093–2117. [Google Scholar] [CrossRef]
  4. Jacobson, M. 100% Clean and Renewable Wind, Water, and Sunlight All-Sector Energy Roadmaps for 139 Countries of the World. Joule 2017, 1, 108–121. [Google Scholar] [CrossRef] [Green Version]
  5. Jacobson, M.Z. Roadmaps to Transition Countries to 100% Clean, Renewable Energy for All Purposes to Curtail Global Warming, Air Pollution, and Energy Risk. Earth’s Futur. 2017, 5, 948–952. [Google Scholar] [CrossRef]
  6. Davis, S.J.; Lewis, N.S.; Shaner, M.; Aggarwal, S.; Arent, D.; Azevedo, I.L.; Benson, S.M.; Bradley, T.; Brouwer, J.; Chiang, Y.M.; et al. Net-zero emissions energy systems. Science (80-.) 2018, 360. [Google Scholar] [CrossRef] [Green Version]
  7. Ram, M.; Bogdanov, D.; Aghahosseini, A.; Gulagi, A.; Oyewo, A.; Child, M.; Caldera, U.; Sadovskaia, K.; Farfan, J.; Barbosa, L.; et al. Global Energy System Based on 100% Renewable Energy-Power, Heat, Transport and Desalination Sectors; Lappeenranta University of Technology and Energy Watch Group: Lappeenranta, Finland; Berlin, Germany, 2019. [Google Scholar]
  8. Bogdanov, D.; Gulagi, A.; Fasihi, M.; Breyer, C. Full energy sector transition towards 100% renewable energy supply: Integrating power, heat, transport and industry sectors including desalination. Appl. Energy 2021, 283, 116273. [Google Scholar] [CrossRef]
  9. Moriarty, P.; Honnery, D. Feasability of a 100% Global Renewable Energy System. Energies 2020, 13, 5543. [Google Scholar] [CrossRef]
  10. Capellán-Pérez, I.; de Blas, I.; Nieto, J.; de Castro, C.; Miguel, L.J.; Carpintero, O.; Mediavilla, M.; Lobejon, L.F.; Ferreras-alonso, N.; Rodrigo, P.; et al. MEDEAS: A new modeling framework integrating global biophysical and socioeconomic constraints. Energy Environ. Sci. 2020. [Google Scholar] [CrossRef]
  11. WWF. The Energy Report 100% Renewable by 2050; Technical Report; WWF: Gland, Switzerland, 2011. [Google Scholar]
  12. Nieto, J.; Miguel, L.J.; Blas, I.D. Macroeconomic modelling under energy constraints: Global low carbon transition scenarios. Energy Policy 2020, 137. [Google Scholar] [CrossRef] [Green Version]
  13. Sers, M.R.; Victor, P.A. The Energy-emissions Trap. Ecol. Econ. 2018, 151, 10–21. [Google Scholar] [CrossRef]
  14. Dupont, E.; Jeanmart, H.; Possoz, L.; Fagnart, J.F.; Germain, M. Transition Énergétique et (dé) Croissance Économique; Université Catholique de Louvain, Institut de Recherches Economiques et Sociales (IRES): Ottignies-Louvain-la-Neuve, Belgium, 2017. [Google Scholar]
  15. Dale, M.; Krumdieck, S.; Bodger, P. Global energy modelling—A biophysical approach (GEMBA) Part 1: An overview of biophysical economics. Ecol. Econ. 2012, 73, 152–157. [Google Scholar] [CrossRef]
  16. Dale, M.; Krumdieck, S.; Bodger, P. Global energy modelling—A biophysical approach (GEMBA) Part 2: Methodology. Ecol. Econ. 2012, 73, 158–167. [Google Scholar] [CrossRef]
  17. Fagnart, J.F.; Germain, M.; Peeters, B. Can the Energy Transition Be Smooth? A General Equilibrium Approach to the EROEI. Sustainability 2020, 12, 1176. [Google Scholar] [CrossRef] [Green Version]
  18. Court, V.; Jouvet, P.A.; Lantz, F. Long-term endogenous economic growth and energy transitions. Energy J. Int. Assoc. Energy Econ. 2018, 39. [Google Scholar] [CrossRef]
  19. Režný, L.; Bureš, V. Energy Transition Scenarios and Their Economic Impacts in the Extended Neoclassical Model of Economic Growth. Sustainability 2019, 11, 3644. [Google Scholar] [CrossRef] [Green Version]
  20. Rye, C.D.; Jackson, T. A review of EROEI-dynamics energy-transition models. Energy Policy 2018, 122, 260–272. [Google Scholar] [CrossRef]
  21. Dupont, E.; Koppelaar, R.; Jeanmart, H. Global available wind energy with physical and energy return on investment constraints. Appl. Energy 2018, 209, 322–338. [Google Scholar] [CrossRef]
  22. Dupont, E.; Koppelaar, R.; Jeanmart, H. Global available solar energy under physical and energy return on investment constraints. Appl. Energy 2020, 257, 113968. [Google Scholar] [CrossRef]
  23. Fagnart, J.F.; Germain, M. Net energy ratio, EROEI and the macroeconomy. Struct. Chang. Econ. Dyn. 2016, 37, 121–126. [Google Scholar] [CrossRef]
  24. Dupont, E.; Germain, M.; Jeanmart, H. Estimate of the Societal Energy Return on Investment (EROI). Biophys. Econ. Sustain. 2021, 6. [Google Scholar] [CrossRef]
  25. IEA. World Energy Outlook 2020; Technical Report; IEA: Paris, France, 2020. [Google Scholar]
  26. EIA. Annual Energy Outlook 2020; EIA: Washington, DC, USA, 2020. [Google Scholar]
  27. Capellán-Pérez, I.; de Blas, I.; Nieto, J.; de Castro, C.; Miguel, L.J.; Mediavilla, M.; Carpintero, O.; Rodrigo, P.; Frechoso, F.; Cáceres, S. Guiding European Policy toward a Low-Carbon Economy. Modelling Sustainable Energy System Development under Environmental and Socioeconomic Constraints. D4.1 (D13) Global Model: MEDEAS-World Model and IOA Implementation ot Global Geographical Level; Technical Report; Pg. Marítim de la Barceloneta: Barcelona, Spain, 2017. [Google Scholar]
  28. OECD. Input-Output Tables; OECD: Paris, France, 2015. [Google Scholar]
  29. Kümmel, R.; Lindenberger, D.; Weiser, F. The economic power of energy and the need to integrate it with energy policy. Energy Policy 2015, 86, 833–843. [Google Scholar] [CrossRef]
  30. Ayres, R.U.; van den Bergh, J.C.J.M.; Lindenberger, D.; Warr, B. The underestimated contribution of energy to economic growth. Struct. Chang. Econ. Dyn. 2013, 27, 79–88. [Google Scholar] [CrossRef]
  31. Court, V.; Fizaine, F. Energy expenditure, economic growth, and the minimum EROI of society. Energy Policy 2016, 95, 172–186. [Google Scholar] [CrossRef]
  32. IPCC. Global Warming of 1.5 °C; Technical Report; Intergovernmental Panel on Climate Change: Geneva, Switzerland, 2018. [Google Scholar]
  33. IEA. Key World Energy Statistics 2020; Technical Report August; IEA: Paris, France, 2020. [Google Scholar]
Figure 1. Diagram of the economy, closed and with two sectors: the ES producing energy, and the FS producing general-purpose goods. Y i represent the output flows, E i represent the intermediate consumptions of energy (EIC), X i the non-energy intermediate consumptions (NEIC), and K i the stocks of capital used by each sector.
Figure 1. Diagram of the economy, closed and with two sectors: the ES producing energy, and the FS producing general-purpose goods. Y i represent the output flows, E i represent the intermediate consumptions of energy (EIC), X i the non-energy intermediate consumptions (NEIC), and K i the stocks of capital used by each sector.
Sustainability 13 05345 g001
Figure 2. The evolution of the technical parameters of the renewable energy subsector as a function of Y r e , the amount of renewable energy produced (in exajoules produced per year): (a): x r e ( Y r e ) , (b): v r e ( Y r e ) .
Figure 2. The evolution of the technical parameters of the renewable energy subsector as a function of Y r e , the amount of renewable energy produced (in exajoules produced per year): (a): x r e ( Y r e ) , (b): v r e ( Y r e ) .
Sustainability 13 05345 g002
Figure 3. Representation of the AS and of a TC, for a fictitious scenario with Y e and q f being fixed. The solid curves represent the bounds of the interval for k (Equation (28)), and the AS is the space between these two curves. The dashed curve shows a particular TC and the dots show the trajectory of the ET, thus the discrete pairs ( ϰ ( t ) , k ( t ) ) for a few successive years.
Figure 3. Representation of the AS and of a TC, for a fictitious scenario with Y e and q f being fixed. The solid curves represent the bounds of the interval for k (Equation (28)), and the AS is the space between these two curves. The dashed curve shows a particular TC and the dots show the trajectory of the ET, thus the discrete pairs ( ϰ ( t ) , k ( t ) ) for a few successive years.
Sustainability 13 05345 g003
Figure 4. Evolution of the technical parameters of the ES as a function of the penetration rate ϰ ( ϰ = 1 corresponds in this case to Y r e = Y e ( 2018 ) ): (a): v e ( ϰ ) , (b): x e ( ϰ ) , (c): q e ( ϰ ) .
Figure 4. Evolution of the technical parameters of the ES as a function of the penetration rate ϰ ( ϰ = 1 corresponds in this case to Y r e = Y e ( 2018 ) ): (a): v e ( ϰ ) , (b): x e ( ϰ ) , (c): q e ( ϰ ) .
Sustainability 13 05345 g004
Figure 5. Exercise 1, evolution as a function of θ of: (a): ϰ ( 2100 ) , (b): the asymptotical value ϰ l , and, only for the SD scenario and θ 0.86 , (c): the year of the end of the transition T f , and (d): the residual capital growth rate.
Figure 5. Exercise 1, evolution as a function of θ of: (a): ϰ ( 2100 ) , (b): the asymptotical value ϰ l , and, only for the SD scenario and θ 0.86 , (c): the year of the end of the transition T f , and (d): the residual capital growth rate.
Sustainability 13 05345 g005
Figure 6. Temporal evolution of economic variables for Exercise 2, scenario SD, with g K = 3% and θ = 0.418. (a): the renewable penetration rate ϰ ( t ) , (b): the savings rate s ( t ) , (c): the share of capital stock devoted to the ES μ ( t ) , (d): the energy price compared to its price in 2018 p ( t ) / p ( 2018 ) , (e): the share of the ES gross added value in GDP α ( t ) (blue curve), and the share of the ES net added value in NDP α n ( t ) (red curve), and (f): the FS output Y f ( t ) and its distribution: NEIC X ( t ) = X f ( t ) + X e ( t ) , final consumption C f ( t ) and investment I ( t )
Figure 6. Temporal evolution of economic variables for Exercise 2, scenario SD, with g K = 3% and θ = 0.418. (a): the renewable penetration rate ϰ ( t ) , (b): the savings rate s ( t ) , (c): the share of capital stock devoted to the ES μ ( t ) , (d): the energy price compared to its price in 2018 p ( t ) / p ( 2018 ) , (e): the share of the ES gross added value in GDP α ( t ) (blue curve), and the share of the ES net added value in NDP α n ( t ) (red curve), and (f): the FS output Y f ( t ) and its distribution: NEIC X ( t ) = X f ( t ) + X e ( t ) , final consumption C f ( t ) and investment I ( t )
Sustainability 13 05345 g006
Figure 7. Temporal evolution of economic variables for Exercise 2, scenario SD, with g K = 3% and θ = 0.418 (continuation and end). (a): the gross domestic product GDP(t) and its distribution between consumption C ( t ) and investment I ( t ) , (b): the distribution of the ES output Y e ( t ) , in relative terms, between self-consumption ( E e ( t ) / Y e ( t ) ), deliveries to the FS ( E f ( t ) / Y e ( t ) ), and final energy consumption ( C e ( t ) / Y e ( t ) ), (c): the capital intensity of the economy (i.e., the capital/GDP ratio) v ( t ) , (d): the capital growth rate g K ( t ) (constant) and its distribution between the GDP growth rate g ( t ) (decreasing) and the growth rate of the capital intensity of the economic g v ( t ) (increasing), (e): the ratios of the GDP to the total output of the economy (defined as Y T ( t ) = Y f ( t ) + p Y e ( t ) ) (blue curve), and the ratio between the NDP and GDP (red curve), and (f): the point-of-use societal EROI, EROI ( t ) .
Figure 7. Temporal evolution of economic variables for Exercise 2, scenario SD, with g K = 3% and θ = 0.418 (continuation and end). (a): the gross domestic product GDP(t) and its distribution between consumption C ( t ) and investment I ( t ) , (b): the distribution of the ES output Y e ( t ) , in relative terms, between self-consumption ( E e ( t ) / Y e ( t ) ), deliveries to the FS ( E f ( t ) / Y e ( t ) ), and final energy consumption ( C e ( t ) / Y e ( t ) ), (c): the capital intensity of the economy (i.e., the capital/GDP ratio) v ( t ) , (d): the capital growth rate g K ( t ) (constant) and its distribution between the GDP growth rate g ( t ) (decreasing) and the growth rate of the capital intensity of the economic g v ( t ) (increasing), (e): the ratios of the GDP to the total output of the economy (defined as Y T ( t ) = Y f ( t ) + p Y e ( t ) ) (blue curve), and the ratio between the NDP and GDP (red curve), and (f): the point-of-use societal EROI, EROI ( t ) .
Sustainability 13 05345 g007
Table 1. The calibration of the technical parameters and of the EROI for the reference year (2018).
Table 1. The calibration of the technical parameters and of the EROI for the reference year (2018).
v f [-]1.72
q f [MJ/US$ 2010]1.7
x f [-]0.47
v e [US$2010/(GJ/y)]39.07
q e [%]9.8
x e [US$2010/GJ]3.32
EROI [-]8.5
Table 2. An estimation of the technical parameters of both subsectors of the ES in 2018, corresponding to a penetration rate ϰ of 2%.
Table 2. An estimation of the technical parameters of both subsectors of the ES in 2018, corresponding to a penetration rate ϰ of 2%.
Energy SectorNon-Renewable
Subsector
Renewable
Subsector
q e [%]9.769.788.68
x e [US$2010/GJ]3.323.323.01
v e [US$2010/(GJ/y)]39.0726.65647.63
Table 3. Exercise 2: values for θ = 0 and θ = 1 of ϰ , s (in 2050, 2100 and T f ) and T f , for the two scenarios.
Table 3. Exercise 2: values for θ = 0 and θ = 1 of ϰ , s (in 2050, 2100 and T f ) and T f , for the two scenarios.
BAUSD
θ = 0 θ = 1 θ = 0 θ = 1
20502100 T f
2129
20502100 T f
2133
2050 T f
2085
2050 T f
2090
ϰ [%]16.846.110001710031.21000100
s [%]29.834.84424.425.941.130.740.123.936.8
Table 4. Minimum capital growth rate for a transition completed by 2070 (Exercise 2, scenario SD), and the corresponding savings rate at the end of the ET.
Table 4. Minimum capital growth rate for a transition completed by 2070 (Exercise 2, scenario SD), and the corresponding savings rate at the end of the ET.
θ g K [%] s ( T f ) [%]
02.9143.5
0.22.9643.1
0.4342.5
0.63.0442
0.83.0841.5
1.03.1241.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dupont, E.; Germain, M.; Jeanmart, H. Feasibility and Economic Impacts of the Energy Transition. Sustainability 2021, 13, 5345. https://0-doi-org.brum.beds.ac.uk/10.3390/su13105345

AMA Style

Dupont E, Germain M, Jeanmart H. Feasibility and Economic Impacts of the Energy Transition. Sustainability. 2021; 13(10):5345. https://0-doi-org.brum.beds.ac.uk/10.3390/su13105345

Chicago/Turabian Style

Dupont, Elise, Marc Germain, and Hervé Jeanmart. 2021. "Feasibility and Economic Impacts of the Energy Transition" Sustainability 13, no. 10: 5345. https://0-doi-org.brum.beds.ac.uk/10.3390/su13105345

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop