Iñigo
Capellán-Pérez
*ab,
Ignacio
de Blas
ab,
Jaime
Nieto
ac,
Carlos
de Castro
ad,
Luis Javier
Miguel
ab,
Óscar
Carpintero
ac,
Margarita
Mediavilla
ab,
Luis Fernando
Lobejón
ac,
Noelia
Ferreras-Alonso
e,
Paula
Rodrigo
a,
Fernando
Frechoso
a and
David
Álvarez-Antelo
a
aResearch Group on Energy, Economy and System Dynamics, Escuela de Ingenierías Industriales, University of Valladolid, Paseo del Cauce s/n, 47011 Valladolid, Spain. E-mail: inigo.capellan@uva.es
bDepartment of Systems Engineering and Automatic Control, Escuela de Ingenierías Industriales, University of Valladolid, Paseo del Cauce s/n, 47011 Valladolid, Spain
cDepartment of Applied Economics, University of Valladolid, Av. Valle Esgueva 6, Spain
dDepartment of Applied Physics, Escuela de Arquitectura, University of Valladolid, Av Salamanca, 18, 47014, Valladolid, Spain
eCARTIF Foundation, Parque Tecnológico de Boecillo, Boecillo, Valladolid 47151, Spain
First published on 25th February 2020
A diversity of integrated assessment models (IAMs) coexists due to the different approaches developed to deal with the complex interactions, high uncertainties and knowledge gaps within the environment and human societies. This paper describes the open-source MEDEAS modeling framework, which has been developed with the aim of informing decision-making to achieve the transition to sustainable energy systems with a focus on biophysical, economic, social and technological restrictions and tackling some of the limitations identified in the current IAMs. MEDEAS models include the following relevant characteristics: representation of biophysical constraints to energy availability; modeling of the mineral and energy investments for the energy transition, allowing a dynamic assessment of the potential mineral scarcities and computation of the net energy available to society; consistent representation of climate change damages with climate assessments by natural scientists; integration of detailed sectoral economic structure (input–output analysis) within a system dynamics approach; energy shifts driven by physical scarcity; and a rich set of socioeconomic and environmental impact indicators. The potentialities and novel insights that this framework brings are illustrated by the simulation of four variants of current trends with the MEDEAS-world model: the consideration of alternative plausible assumptions and methods, combined with the feedback-rich structure of the model, reveal dynamics and implications absent in classical models. Our results suggest that the continuation of current trends will drive significant biophysical scarcities and impacts which will most likely derive in regionalization (priority to security concerns and trade barriers), conflict, and ultimately, a severe global crisis which may lead to the collapse of our modern civilization. Despite depicting a much more worrying future than conventional projections of current trends, we however believe it is a more realistic counterfactual scenario that will allow the design of improved alternative sustainable pathways in future work.
Broader contextBy substantially degrading the natural life-support systems and processes that sustain our existence, we are jeopardizing the survival of our societies. Holistic frameworks are required to assess the urgent and radical changes needed to change track and achieve a sustainable future. This study describes the open-source MEDEAS integrated assessment modeling framework, which has been designed to facilitate the assessment of policies for the best possible energy transition path during the next decades, with a focus on biophysical, economic, social and technological restrictions. The potentialities and novel insights that this framework brings are illustrated by the simulation of four variants of the continuation of current trends with the MEDEAS-World model. The consideration of alternative plausible assumptions and methods, combined with the feedback-rich structure of the model, reveal dynamics and implications absent in classical models. Our results show that the continuation of current trends will derive in biophysical scarcities and impacts which will most likely derive in regionalization, conflict, and ultimately global crisis, leading to the collapse of our modern civilization. Despite depicting a much more worrying future than conventional projections of current trends, this scenario seems a more realistic counterfactual scenario that will allow the design of improved alternative sustainable pathways in future work. |
In this context, environmental integrated assessment modeling stands among the most practical scientific approaches. It refers to any type of analysis that seeks to integrate multiple disciplines and dimensions aiming to capture interactions between human and natural systems – which tend to be complex, dynamic and highly non-linear-, with the goal of providing useful information for policy making. Integrated Assessment Models (IAMs) (or energy–economy–environment (E3) models), are computer programs that link an array of component models based on mathematical representations of information from various disciplines. A great diversity of environmental IAMs exists due to the different approaches used by modelers striving to capture the complex interactions and high uncertainties involved in the environment/economy/society interface, with the dominance of IAMs focused on climate change research.8–12
In fact, dozens, if not hundreds, of IAMs have been developed in the last few decades since the pioneer World3 was developed in the early 1970s.13,14 Despite great advances achieved in the field over the years,8,10 most IAMs (and especially those more policy-influential), share a core set of common assumptions whose validity is being disputed in the scientific discussion. First, IAMs are generally characterized by a rather sequential structure with limited feedbacks among the represented subsystems. The interconnectivity of modules has likely being constrained by the historical development of most IAMs through linkage of existing modules which were not originally designed for being interlinked.15 For example, natural science models must respect the laws of thermodynamics, while economic models often do not. Also, the discrepancy between the natural scientists’ understanding of ecological feedbacks and the representations of environmental damage found in IAMs (if any) is especially relevant for the case of climate change impacts. Most IAMs fail to capture the “potentially irreversible threat to human societies and the planet” stated, for example, in the Paris Agreement.4,9,16–21 Second, a lack of plurality in the methods to represent the economic dimension has been detected in the literature, dominated by assumptions of conventional general or partial equilibrium through optimization methods, perfect factor substitutability, as well as the widespread use of prices as indicators of scarcity. These simplifications fail to capture the relevance of sector complementarities within the economic structure, the socioeconomic system dynamics and the role of macroeconomic policies for sustainability governance.22–28 Third, the abundance of both fossil fuels and renewable energy sources is a default assumption in most of the prominent IAMs used for climate policy analysis; hence, future energy transitions are thus largely modeled as demand-driven transformations only constrained by available monetary investments.8,29,30 However, this assumption is disputed by studies in the literature showing that fossil fuels’ extraction might face significant constraints in the next few decades related with increasing geological restrictions as the quality of the resource decreases.30–32 Furthermore, a branch of literature is also showing that the replacement of fossil fuels in the current socioeconomic system by the large scale deployment of RES faces serious challenges in relation to biophysical factors such as intermittency or mineral and land requirements.33–41 Fourth, most IAMs disregard the implications that the future energy investments required to achieve the transition to renewables may have for the system.42–46 In fact, a favorable energy return on energy invested (EROI) (energy surplus) is a critical aspect of the viability of societies and has been associated with such fields as biology or anthropology as a key driver of increasing complexity and evolution for plants, animals and humans.47–49 Finally, Fifth, (the lack of) transparency has been highlighted as being an issue in the field of IAMs critically affecting credibility and robustness of the results disseminated.11,50,51
This study is dedicated to the description of the MEDEAS modeling framework which has been designed with the aim of facilitating the assessment of policies to lead to the best possible energy transition path during the next few decades, with a focus on biophysical, economic, social and technological restrictions and tackling the aforementioned limitations in the current state-of-the art of IAMs. MEDEAS models have been developed within the homonymous project (Modeling Energy System Development under Environmental and Socioeconomic constraints),‡ whose main aim is to provide policy makers and stakeholders with new modeling tools to better assess the impacts and limitations of the EU energy production/consumption system transition to a low-carbon sustainable socio-economy. The MEDEAS models focus on the analysis of strategic, long-term outcomes of the human–nature interface, and have been designed applying system dynamics, which facilitates the integration of knowledge from different perspectives and disciplines, as well as feedbacks from the different subsystems. These models are open-source and can therefore be used by scientists, stakeholders and policy makers. Models at three different geographical aggregated scales have been developed: global, EU and country-level. Given that these geographical scales share a common modeling approach, for the sake of simplicity, in this paper we focus on the global model version. The simulation of a business-as-usual (BAU) scenario within the global MEDEAS model (MEDEAS-World) activating/deactivating some key functionalities shows the potentialities and novel insights that this framework brings when simulating future sustainability pathways.
The paper is structured as follows: Section 2 describes the material and methods, Section 3 shows the simulation results of the 4 selected scenario cases, Section 4 discusses the obtained results and the novel insights of the proposed novel framework, and 5 concludes.
MEDEAS models are structured in nine main modules: economy, energy demand, energy availability, energy infrastructures and EROI, minerals, land-use, water, climate/emissions, and social and environmental impact indicators. Fig. 2 shows the conceptual schematic overview of the global-aggregated scale MEDEAS-W, including the main relationships between the different modules. The main characteristics of each module are:
Fig. 2 MEDEAS-world model schematic overview. The main variables connecting the different modules (represented by colored rectangles) are represented in italics and by solid arrows. The dashed arrow represents the exogenous driver inputs. EROI: energy return on energy investment. RES: renewable energy sources. Source: adaptation from Capellán-Pérez et al.52 |
• Economy: the economy is modeled assuming non-clearing markets (i.e., not forcing general equilibrium), demand-led growth and sector complementarity. Hence, production is determined by final demand and economic structure, combined with supply-side constraints such as energy availability. The economic structure is captured by the adaptation and dynamic integration of global WIOD input-output tables, resulting in 35 sectors (see ESI†) and 4 components of final demand.60
• Energy demand: final energy demand by sector and households is estimated through the projection of sectoral economic production and sectoral final energy intensities, considering efficiency improvements and inter-final energy replacements driven by policies and physical scarcity.
• Energy availability: this module includes the potential and availability of RES and non-renewable energy resources, taking into account biophysical and temporal constraints. In particular, the availability of non-renewable energy resources depends on both stock and flow constraints.61–63 In total, 25 energy sources and technologies, and 5 final energies are considered (electricity, heat, solids, gases and liquids), with large technological disaggregation. The modeling of energy availability is mainly based on the previous model WoLiM.64 The intermittency of RES is considered in the framework, computing endogenous levels of overcapacities, storage and new electrical grids, depending on the penetration of variable RES technologies.
• Energy infrastructures & EROI: this module represents the capacities for generating electricity and heat, considering planning and construction delays. The energy investments for RES to produce electricity are endogenously and dynamically modeled, which allows the Energy Return on Energy Investment (EROI) of individual technologies and the EROI of the whole energy system to be computed. The energy demand is affected by the variation of the EROI of the system. Transportation is modeled in great detail, differentiating between different types of vehicles for households, as well as freight and passenger inland transport.
• Minerals: minerals are required by the economy, including those required for the construction, operation and maintenance of alternative energy technologies. Recycling policies are available.
• Land-use: this module mainly accounts for the additional land requirements of RES.
• Water: this module projects water use by type (blue, green and gray) by economic sector and for households.
• Climate/emissions: the global model computes the climate change levels due to the greenhouse gas (GHG) emissions generated by human societies (non-CO2 emissions are exogenously set, taking RCP scenarios as reference 65). The carbon and climate cycle is adapted from C-ROADS.66,67 This module includes a damage function which translates increasing climate change levels into damages to human systems. In regional models, domestic GHGs are computed and climate damages are dependent on the global average surface temperature's change from the associated global model scenario.
• Social and environmental impacts: this module translates the “biophysical” results of the simulations into metrics related to social and environmental impacts. The objective of this module is to contextualize the implications for human societies in terms of well-being and environmental impacts for each simulation.
The modules have different levels of development, being the economy, climate and those related to energy the ones most detailed. Minerals, land-use and water modules are more stylized representations focused to compute social and environmental Impacts which however do not feedback to the rest of the system (see Fig. 2).
MEDEAS models dynamically operate as follows. For each period, a sectoral economic demand is estimated from exogenous pathways of expected Gross Domestic Product per capita (GDPpc), population and income distribution. The final energy demand required to fulfill production is obtained using hybrid (energy–economy) input–output analysis, combining monetary output and energy intensities by type of final energy. The final energy supply computed within the energy availability and energy infrastructures and EROI modules may satisfy (or not) the required demand. In case of energy scarcity of any final energy type, two phenomena occur simultaneously: (1) final energy replacement to shift towards more abundant final energy types, and (2) adaptation of the economic production to the available energy. The materials required by the economy with emphasis on those required by alternative green technologies are estimated, which allows eventual future mineral bottlenecks to be assessed. However, for the sake of simplicity, mineral availability does not constrain economic output in the current model versions. The new energy infrastructures require energy investments, whose computation allows us to estimate the variation of the EROI of the system, which in turn affects the final energy demand. The climate module computes the GHG emissions associated to the resulting energy mix (complemented by exogenous pathways for non-energy emissions), whose accumulation leads to a certain level of climate change, which in turn feeds back to the economy affecting final demand. Additional land and water requirements are accounted for. Finally, a set of social and environmental impacts is computed. The simulation of regional models requires exogenous inputs provided by the parent models, as well as defining assumptions on multiregional monetary and energy trade.
The remainder of the section focuses on the description of the main modules of MEDEAS-W.
Fig. 3 shows the schematic overview of the Economy module of the MEDEAS framework. The production of goods and services is determined by the variation over time of the aggregate demand. At the beginning of any period, the variation of final demand based on population and GDPpc change, together with income distribution scenarios, using labor and capital share over GDP, allows the expected GDP to be set, as well as the expected profits and the expected total wage bill. The final demand function enables the estimation of the monetary demand by component of final demand for each sector (see eqn (1) and ESI†). Household consumption (hh) and Gross Fixed Capital Formation (gfcf) are estimated with behavioral equations based on panel data regressions, while Government expenditures (ge) and changes in inventories (invent) remain as a constant share of the sectoral final demand. Total final demand by sector (fdi) is thus the sum of each institutional sector's demand (eqn (1)).
fdi = hhi + gfcfi + gei + inventi; i = ∈1,…,35 sectors | (1) |
Fig. 3 Schematic overview of the Economy module of the MEDEAS framework. The dashed arrow represents the exogenous driver inputs. IOA: input–output analysis. Source: own work. |
Households demand (hh) depends (eqn (2)) on its main source of income, wages, i.e. labour compensation (Lab). Investment, or gross fixed capital formation (gfcf), is determined (eqn (3)) by the expected profits, i.e. capital compensation (Cap). β2 coefficients in both equations were estimated through panel data regressions which proved consistent and well-fitted. β1i coefficients allow for considering the individual (each sector) effects. Interest rates were considered also as an explanatory variable of investment but were found not significant:
ln(hhi) = β0 + β1iSeci + β2ln(Lab) | (2) |
ln(gfcfi) = β0 + β1iSeci + β2ln(Cap) | (3) |
Thereafter, IOA is applied to compute the sectoral production xi required to fulfill the expected demand. IOA is based on the standard accountability identities whereby production equals intermediate and final demand. In this method, A is the squared technical coefficients matrix representing the combination of inputs required to produce the output (x), and thus A·x amounts to the intermediate monetary final demand (fd). Therefore, x = A·x + fd and thus, solving for x = (I − A)−1·fd, and renaming the first component as the so-called Leontief Matrix L, x = L·fd.¶73 The Leontief matrix L allows the model to estimate both direct and indirect production carriers of final monetary demand, measuring the production sensitivity to final monetary demand. See Nieto et al.72 for the analysis and discussion of different economic structure scenarios by setting a target A matrix to which the current one gradually evolves.
The estimation of the final energy demand (fed) in MEDEAS is performed through a general method based on projecting energy intensities by the final energy k (solids, liquids, gases, electricity and heat) (see Section 2.2.2, although bottom-up modelling is also possible, see Section 2.2.6). Thus, by multiplying the 35 sectors’ energy intensities, broken down by final energy type (eik, see eqn (4)) with the sectoral production xi, and adding the component of households obtained as their energy intensity (e_hhk, see eqn (5)) multiplied by its final consumption hh, we obtain the total final energy demand by final energy tfedk required to produce the economic output.
(4) |
(5) |
(6) |
(7) |
As mentioned above (see Section 2.2.1), the MEDEAS framework considers 5 types of final energies (solids, liquids, gases, electricity and heat) and 35 economic sectors, which allows us to calculate 35 × 5 sectoral energy intensities eik and 5 energy intensities associated to households consumption e_hhk. Hence, a total of 180 final energy intensities are handled in the MEDEAS framework. Except for inland transport sectors, which follow a bottom-up modeling (see Section 2.2.6) which could be eventually expanded to other sectors, the variation over time of each of the sectoral and households’ final energy intensities is modeled in a top-down approach.
The starting point for modeling the dynamic behavior of final energy intensities is the available historical data, taking as reference the data from the WIOD database environmental accounts (1995–2009)78 and complemented with the IEA balances.79 In fact, energy intensities have generally declined over the last few decades.80–82 The baseline trends assume variation rates of final energy efficiencies linked to historical technological change and economic growth,77 and are based on the model proposed by Schenk and Moll83 to explain the historical trend in energy intensity, considering biophysical and thermo-dynamical limits in the substitution of inputs.
The final energy intensities are assumed to change with reference to the baseline trends due to two factors:
(1) The variation in the energy intensity due to the improvement in energy efficiency, associated with the technology used in the consumption of each sector and type of final energy (Δeeffik).
(2) The variation in energy intensity due to the substitution of one type of final energy by another (Δesubik). In this case, the type of energy that is replaced decreases its final energy intensity and the type of energy that replaces the previous one increases its final energy intensity.
The variations of the sectoral final energy intensities are assumed to be driven by two main factors (see eqn (8) and (9)):
(a) Market factors related to the scarcity of each type of energy. The scarcity of a type of energy can lead to greater efforts to increase energy savings and improve efficiency, as well as the substitution of that type of energy by more abundant ones. Indicators of physical scarcity have been specifically developed to represent physical supply-demand imbalances. Market factors are modeled through the perception of final energy scarcity (PSk).
(b) Sociopolitical factors. Some policy measures, such as climate change mitigation policies, are promoting a greater effort in energy efficiency and fostering the substitution of fossil fuels.
Δeeffik = (PSk + effects of policies) × Maxeffik | (8) |
Δesubik = (PSk + effects of policies) × Maxsubik | (9) |
An estimate of the maximum variations (Maxeffik and Maxsubik) of the energy intensities has been obtained through a statistical analysis based on the historic data from WIOD.84 An analogous approach has been applied to modeling the final energy demand of households’ consumption.
The perception of scarcity of each type of final energy is dynamically estimated considering the eventual gap between supply and demand for each type of final energy εk. How to measure the scarcity or abundance of natural resources has been a controversial issue in economics for a long time.85 The mainstream approach considering prices as a reliable indicator of scarcity of natural resources has been criticized, given its theoretical and empirical weaknesses. Energy and mineral prices are subject to multiple influences (institutional framework, oligopolistic market structure, etc.), which prevent perfect competition from happening in both the short and long-term.26,86 Moreover, given the inertia and rigidities in the productive processes highly dependent on natural resources, important adjustments in the economic system are produced with quantity changes (instead of prices), as post-Keynesian approaches have highlighted.25 Thus, MEDEAS applies an alternative “biophysical” perspective to model inter-final energy substitution, which takes into account dynamically the extraction of natural resources and their physical scarcity (eqn (10)).61,87,88
(10) |
PSk(t) = scarcityk × SS + PSk(t − 1)/FF | (11) |
MEDEAS final energy category | NRE/RES | Primary sources of energy and technologies modeled in MEDEAS |
---|---|---|
a MEDEAS does not differentiate between different temperature levels for heat. | ||
Solids | NRE | Coal |
Peat | ||
Waste | ||
RES | Charcoal | |
Biomass (modern) | ||
Biomass (traditional) | ||
Liquids | NRE | Conventional oil |
Unconventional oil | ||
CTL (Coal to Liquids) | ||
GTL (Gas to Liquids) | ||
RES | Biomass (biofuels) | |
Gases | NRE | Conventional natural gas |
Unconventional natural gas | ||
RES | Biogas | |
Electricity | NRE | Natural gas |
Oil | ||
Coal | ||
Uranium | ||
Waste | ||
RES | Hydroelectricity | |
Geothermal power | ||
Biomass | ||
Oceanic (Wave, tidal, OTEC) | ||
Onshore wind | ||
Offshore wind | ||
Solar PV | ||
Solar CSP | ||
Biogas | ||
Heata | NRE | Coal |
Natural gas | ||
Oil | ||
Nuclear (cogeneration) | ||
Waste | ||
RES | Geothermal | |
Solar thermal | ||
Biomass | ||
Biogas |
Oil, natural gas, coal and uranium are the main non-renewable energy (NRE) resources considered in MEDEAS. Oil and natural gas are disaggregated between conventional and unconventional resources, distinction especially relevant to account for their different GHG emission factors. The availability of non-renewable energy resources in MEDEAS depends upon two constraints: stock (available resource in the ground) and flow (extraction rate of this resource). In fact, the constrained flow rates from deposits to consumers have been shown to be a more relevant limiting factor to follow demand than the remaining resource in situ. Technology can help regulate the extraction rate, but the latter is ultimately bound by geological-physical constraints.61–63,89
In order to consider the future availability of fossil fuels in the model, we performed a literature review of studies providing depletion curves over time, taking into account both stock and flow limit.63,90–104 These curves represent extraction levels compatible with geological constraints as a function of time. Depletion curves, rather than predictions, represent maximum extraction profiles for a fraction of the resource base estimated to be economically extractable in the future considering geological constraints: the actual rate of resource consumption might be affected by such variables as geopolitics, economic conditions and technology development. Given that demand is endogenously modeled for each resource, these depletion curves are transformed so as to be incorporated in the dynamic model, being converted into maximum production curves as a function of remaining resources (see Fig. 4 and Appendix B in Capellán-Pérez et al.,52 for details). The model user is free to select any depletion curve, introduce a new one, or even assume that no relevant constraints will affect the supply of any of the NRE resources in the simulation period.
Fig. 4 Example of the transformation of a depletion curve (left) in maximum extraction curve (right) suppressing time dependence as a function of remaining ultimately recoverable resources (URR). The depletion curve as example is total natural gas from Laherrère.105 Source: own work. |
Nine types of RES for electricity generation are modeled: hydro, solar PV, solar CSP, onshore wind, offshore wind, geothermal, biomass, oceanic and biogas; while for heat 4 sources are considered: solar, biomass, geothermal and biogas.
Bioenergy is modeled in four main categories: (1) traditional biomass, (2) modern solid biomass for heat and electricity, (3) residues, and (4) dedicated crops in current croplands and in marginal lands (i.e., land that was previously used for agriculture or pasture but that has been abandoned and not converted to forest or urban areas106) for biofuels. It is possible to activate the availability of third generation ligno-cellulosic fuels for a given year (see Capellán-Pérez et al.,52 for details). Special attention is devoted to the land requirements of RES technologies (see Section 2.2.7), given that the transition to RES will intensify the competition for land globally.38,107,108 In this sense, the potential of some RES is defined by the user as a function of the land assessed to be available in the future (solar on land, cropland for liquid biofuels and marginal lands), while for the rest of RES the potential is directly set in power terms by the user.
The selection of new alternative technologies still not competitive modeled in MEDEAS takes as reference the precautionary principle, which is the most robust approach in uncertainty contexts such as the one characterizing climate change and the sustainability crisis.68 Two criteria have been applied:52
(1) Given the need for urgent action to stabilize climate and reverse current unsustainable trends, focus on those technologies that are currently available and proven to be commercial. In fact, it has been shown that new technologies and energy systems take several decades to diffuse through the economy.109,110
(2) Ensure that the net energy balance of the considered technologies is positive, i.e., that the technology will be a “reasonable” net energy contributor to society. Although an energy system with EROI at the point of use <1–2 could still be used for some specific purposes, it would rather be an anecdotal technology; given the burden it would impose on the whole energy system (it would be an energy drain rather than a source).
In the light of these criteria, carbon capture and storage,111,112 Negative Emission Technologies (NETs),112–117 hydrogen, nuclear fast breeders and nuclear fusion are not included in the current MEDEAS framework.††118
This split is required to deal with the fact that data concerning mineral consumption globally are generally of low quality (lack of robust data and sectoral detail, etc.).40,119 Given the lack of data on material demand associated to the WIOD sectors, a stylized approach was applied to estimate the future demand of minerals by the rest of the economy, assuming a linear dependence with GDP.120 Historical data for the global GDP121 and the extraction of minerals122 for the period 1994–2015 were applied to estimate the parameters of the regressions, obtaining acceptable correlations in most cases.46
For each of the key “green” technologies, a representative technology was selected, taking into account the present and foreseen performance while avoiding those more likely to be affected by scarce minerals in the future.‡‡ Subsequently, a literature review was performed to identify the material intensity (kg MW−1) of each representative technology, including those related to additional grid requirements, for a total of 58 materials, of which 19 are minerals (aluminum, cadmium, chromium, copper, gallium, indium, iron, lead, lithium, magnesium, manganese, molybdenum, nickel, silver, tellurium, tin, titanium, vanadium and zinc). Selection criteria were based on considering all relevant materials so as to accurately estimate the embodied energy for the EROI estimation (see Section 2.2.5), with potential critical minerals having already been identified in the literature,40,127–130 as well as specific assessments.52 A comprehensive literature review, complemented by our own estimations, was performed to collate the most robust and accurate data concerning material requirements for each technology. In the case of uncertainty about potential double accounting, material requirements were not included. Hence, these estimations can be considered conservative/optimistic.
Given the existing uncertainties in future mineral availability,131–133 MEDEAS compares the cumulated primary demand of minerals (after accounting for recycling rates) with the current estimated level of their geological availability (reserves and resources) for a qualitative detection of risks to the future mineral supply. Hence, potential mineral scarcity in the current model version is not feed-backed and does not affect the rest of the model, as is the case with energy (see Section 2.2.1) (i.e., mineral supply is assumed to always fulfils demand).§§
The net energy available, and not the primary one, is the relevant factor to society.46–49 Acknowledging the difficulties of computing broad boundaries for the EROI,124,134–137 a first, conservative step is taken that focuses on the dynamic computation of the EROI of the system from a standard (EROIstsystem) approach. The standard approach includes the energy requirements to get energy (e.g., build, operate and dismantle power plants) but omits those energy requirements to deliver energy to the point of use of society (EROIfinal) as well as those energy requirements to build the machines and infrastructure required to construct the machines and infrastructure which allows to make the energy investments (EROIext). The dynamic EROIstsystem¶¶ is defined here as the ratio over time between the final energy delivered to society and the energy investments associated to the up-front energy costs of variable RES power plants before they start to deliver energy, the energy investments associated to operation and maintenance as well for variability management during operation, and the dismantling costs at the end of the lifetime.|||| Given the data intensiveness of this method,46 the focus is on variable RES given their much higher techno-sustainable potential comparing with dispatchable RES.143,144 (see eqn (12) and Fig. 5). The resulting net energy is given by eqn (13):
(12) |
(13) |
Fig. 5 Representation of the energetic metabolism of our society. Grey arrows refer to energy flows that are usable by human societies. The black arrow on the left-hand side (0) is a flux of materials with potential energy which can be transformed into usable energy. Dashed vertical arrows represent energy losses at each phase of the chain (transformation, storage and distribution losses). An exosomatic intermediary (arrows 2, 3, 4 and 5) is always required to transform the potential energy into useful exosomatic energy usable by society (1) (excluding non-energy uses) (size of arrows is not to scale). Source: Capellán-Pérez et al.46 |
The following assumptions are taken to compute the EROIstsystem:
• Dynamic and endogenous calculation of the EROIst of each RES variable technology for electricity generation computing the required energetic costs taking as a starting point the materials required in the construction, operation, maintenance and dismantling phases (see Section 2.2.3) and combining this data with the energy consumption per unit of material consumption from Life Cycle Analysis (LCA).145 Estimates are derived from physical inputs and no indirect estimates based on associated economic costs are considered. For RES dispatchables for electricity generation, EROIst values over the lifetime from the literature are taken as reference for the sake of simplicity.
• The EROIst of non-renewable energy sources, as well as for renewables other than electricity generation is conservatively146,147 assumed to be constant over time for the sake of simplicity.
• The energy requirements associated to the overcapacities, storage and construction of new electric grids related to the increasing penetration of variable RES technologies in the system are endogenously computed by the model.
In principle, we do expect technology improvements in the future, however we also foresee that there will be other relevant factors which in the medium term will tend to offset them as the RES progressively scale-up at large levels and gain substantial shares in the energy mix: (1) decreasing returns in the potential of renewables, i.e., after best places are occupied it is necessary to move to more uneconomical sites,148–150 phenomena which may be worsened by eventual land availability constraints in some cases; (2) the increase in energy requirements due to ore grade decrease of minerals with cumulated extraction;151–153 (3) thermodynamical limits to the continuous reduction of energy investments (e.g. related with limits to substitution); (4) the scarcity of some minerals in the future may drive the shift to more abundant minerals which in turn are generally characterized by a lower performance (e.g., replacement of Ag in solar PV34); and (5) thermodynamical limits from the side of generation, as for example the Benz law or the fact that there are absolute limits to the height of rotors for wind or the limits in the conversion from sunlight to electricity such as the Schokley–Queisser limit. Hence, the current performance parameters for RES technologies are taken as reference in the model as a first approximation and further work will be directed to model these opposed effects.
Summarizing, MEDEAS dynamically accounts for the EROIst of the system as shown in eqn (14):
(14) |
TFEC: total final energy consumption (excluding energy materials for non-energy uses).
TFEIRESelec: total final energy investments for renewable technologies of electricity generation.
TFEIstorageelec: total final energy investments for storage of electricity.
OEU: energy industry own-energy use (excepting for electric renewables to avoid double accounting).
g syst = final to primary energy sources ratio (1)/(0) (see Fig. 5), which make possible to deal with the challenge of computing the EROI of the full energy system in a dynamic framework with a set of final energies.46gtechn = 1 is defined for the computation of the EROI of each of the renewable technologies for the generation of electricity.
The variation of the EROI of the system would imply a variation in the energy intensity of the economic sectors linked to the generation, transformation and transport of energy, with implications for investments and available income to households. However, given the limitations of the WIOD sector categorization (e.g., sector 17 “Electricity, Gas and Water supply”, see ESI†), the effect of the variation of the EROI of the system is modeled instead as a variation of the total final energy required by the system in relation to a reference year (2015). Thus, a decrease (increase) of the EROI in relation to the reference year induces an increase (decrease) of the demand of total final energy. Hence, this approach does not capture the metabolic implications of the drop of the EROI of the system to very low levels.***154
Two types of household vehicles are considered: four-wheelers and two-wheelers; and four types of inland transport vehicles: heavy, light cargo, buses and trains. Four types of energies are considered: liquids, hybrid, gas (only natural gas, since GLP is considered a liquid fuel) and electric (which includes battery electric cars and plug-in hybrids). The focus is on those vehicle types which seem more realistic to represent a significant share of the future transportation mix based on current knowledge about technical performance and foreseen limits. Hence, some of the combinations such as: motorbikes based on gas or hybrid, purely electric heavy trucks,158–160 or hybrid and gas based trains have not been considered in this model version.
In this model version, we assume that the current patterns of mobility are maintained since cultural changes are far from conventional scenarios.161–164 Hence, the total number of vehicles is determined by economic demand (see Section 2.2.1).
The user can set policy targets in terms of target shares of a type of vehicle in a given year. The model translates these into changes to the corresponding final energy intensities of the Households and the Inland Transportation WIOD sector.
The corresponding derivative of the energy intensity by type of final energy k (liquids, electricity, gases) can be expressed as a function of either total households’ consumption (hh) or inland transport sector production (xinlandtransport) and their respective energy demand fed_hhk or fedinlandtransport,k of transport uses, considering the change in the share of vehicles and the technical efficiencies. Eqn (15) shows the expression applied for households’ liquids vehicles, describing the variation of eqn (5):
(15) |
Assuming the continuation of the current mobility patterns, the number and use of vehicles divided by households’ demand can be assumed to be constant (however, the modeling framework can be extended to account for modal shifts and demand-management options, see de Blas et al.165). Therefore, we might define the constants A1 and A2 (computed for the year 2015) (eqn (16)):
(16) |
(17) |
Fig. 6 represents the main elements of the Climate module in MEDEAS-W, which functions as follows: on the one hand, the anthropogenic CO2 emissions, endogenously generated within the model due to energy consumption and land-use changes, enter the carbon cycle, estimating the level of CO2 concentration in the atmosphere. The carbon cycle represents the dynamics between the carbon in the atmosphere, the biosphere (humus and biomass) and the ocean, including temperature feedbacks.
Fig. 6 Simplified structure of the Climate module in MEDEAS-W. The dashed arrow represents the exogenous driver inputs. Source: own elaboration. |
The emissions of the rest of the GHGs are mostly exogenous, except for CH4, for which the share of its emissions associated to the extraction, distribution and combustion of natural gas is endogenously calculated, considering its whole lifetime.176 The user can select the level of future emissions of the rest of the GHGs through the selection of their respective RCP scenario.65 The other GHGs cycles (CH4, N2O, PFCs, SF6 and HFCs) are also explicitly modeled in MEDEAS-W, including the mutual interactions, such as between CH4 and N2O.
The instantaneous warming of all GHGs (total radiative forcing) is aggregated through their respective radiative forcing coefficients, which allows the global-average temperature increase with relation to pre-industrial levels to be computed. The user can modify the equilibrium climate sensitivity, which is set by default to the standard value of 2.9 °C, in response to a doubling equivalent change in radiative forcing. Other climate impacts, such as ocean acidification and sea level change, are also calculated.
The modeling of climate change damages to human societies is a very challenging topic, given the complex uncertainties that pervade the coupled human–Earth system, the long time horizon of the problem, as well as the heterogeneous nature of climate impacts across regions, sectors and generations; additionally, there is no historical analogue for their assessment.17,177 Global surface temperature has already increased around +1 °C over pre-industrial levels and scientific evidence is increasingly indicating that the political targets of +1.5/+2 °C set in Climate Summits, such as that reached in Paris in 2015,178 do not represent a safe aspiration for the long-term sustainability of human societies.4,179,180 On the other hand, the review of the literature shows that there is a large gap between natural scientists’ understanding of climate impacts and their representation in IAMs, where temperature change is typically assumed to increase in business-as-usual scenarios between 4–5 °C by the end of the century, with negligible representation of the impacts and damages on humans and the biosphere.¶¶¶8,17–21
In this context, a novel methodology to introduce climate damages in MEDEAS-W, based on the concept of damage function, has been specifically developed, incorporating novel insights and methods proposed in the discussion/debate around the integration of climate change impacts in standard cost-benefit IAMs, including:
• Consistent representation of climate change damages with climate assessments by natural scientists assured through the application of the top-down inductive approach18,21,182 for the calibration of the damage function through informed-interpretation of “dangerous climate change”.
• Climate impacts affect the drivers of growth (or growth rates) and not just the level of production.17,18,20,21,177,183–185
Two options are available in MEDEAS-W for the user: (1) an energy losses function which reduces the final energy consumption available for production after climate change impacts (dependent on the level of CO2 concentrations), and (2) a monetary damage function which relates average global surface temperature change with the loss of GDPpc due to climate change damages. In both cases, the damage function is calibrated by applying an inductive method to assure consistency with climate change impact assessments. At current “low” climate change levels, the damages can be considered negligible, while at high concentrations/temperature increase, the damages are catastrophic. Non-linear functions are considered (logistic and parabolic, respectively) in order to represent the non-linearity between climate change and the related impacts. The user can also easily perform runs deactivating these damages and running the model without considering climate change impacts.||||||
Acknowledging that welfare is a multidimensional feature which cannot be reduced to a single variable,187 MEDEAS does not report a “unique” variable to measure welfare. Instead, the assessment of each scenario in social terms requires the concurrent analysis of a set of diverse variables from the different dimensions analyzed within the framework (see Table 2). With this complex perspective, trade-offs often arise between different dimensions. In particular, energy is a key variable of the MEDEAS framework which contains valuable information; in fact, an adequate energy supply has been identified as a key prerequisite for economic, cultural and social development in complex societies, especially at lower consumption levels,188–192 and as such it has been included in the SDG7 to “ensure access to affordable, reliable, sustainable and modern energy for all”.††††193
Dimension | Indicator |
---|---|
a Annual GDP per capita represents the per capita monetary measure of the market value of all final goods and services produced in a year. Hence, GDP represents a “purely economic” approach to measure social welfare. We stress that GDP per capita is not, and was never, designed as a measure of social or economic welfare, despite being the most common indicator of progress for policy-makers and Governments. In fact, above a certain level, reductions in GDP may be welfare enhancing.194–197 | |
Economy | GDP per capitaa |
Annual growth rate of real GDP per capita (8.1.1 indicator from the SDGs indicators193) | |
Energy intensity measured in terms of both primary and final energy and gross domestic product (GDP) (7.3.1 indicator from the SDGs indicators193) | |
Jobs associated to RES technologies | |
Energy | Final consumption per capita |
Primary consumption per capita | |
Electricity consumption per capita | |
Consumption of RES per capita | |
Share of RES in total final consumption | |
Annual penetration of RES in the total final and primary energy sources consumption | |
EROIst of individual RES electric technologies and of the whole system | |
Land-use | Land requirements for RES production |
Minerals | Mineral use per capita |
Cumulated primary mineral extraction vs. current reserves | |
Cumulated primary mineral extraction vs. current resources | |
Water | Total water and by type consumption per capita |
Share water use vs. renewable water resources | |
Climate change | CO2, CO2e and GHG emissions per capita |
CO2 emission per unit of value added (9.4.1. indicator from the SDGs Indicators193) | |
Atmospheric GHG concentration levels | |
Temperature increase over pre-industrial levels | |
Synthetic | Potential HDI level given energy use |
In this work, four scenario cases of the BAU narrative are simulated: the reference scenario (Ref) and three extreme sensibility cases deactivating two key features operating in the MEDEAS framework and which are subject to scientific uncertainty and discussion: (1) NRE and RES availability restrictions (_noER), and (2) climate change damages (_noCC): Ref_noER, Ref_noCC and Ref_noER_noCC (see Table 3). This scenario case architecture allows us to identify the implications of considering each hypothesis as well as their interaction, which will be discussed in Section 3.2.4.
Scenario case name | Restrictions to the availability of energy (NRE and RES)? | Climate change damages? |
---|---|---|
Ref | YES | YES |
Ref_noER | NO | YES |
Ref_noCC | YES | NO |
Ref_noER_noCC | NO | NO |
Fossil fuel resource abundance, understood as the vast geological availability accessible at an affordable price, is a default assumption in most IAMs and energy models.29,32,64,210 However, many uncertainties exist in relation to non-renewable energy endowments and their future availability related to technical, geological, economic and political factors. In fact, the IPCC-AR58 highlighted that future non-renewable energy availability has historically been under-studied under the assumption of future abundance. This mainstream view is confronted by the fact that biophysical constraints impose certain absolute limits on the use of non-renewable energies, which cannot be extracted at will.31,32,61,62 While the different estimations for conventional oil and gas tend to converge, respectively, to similar levels, the highest uncertainty is in the future development of unconventional fuels.30,31 The main issue is to what extent technological improvements will be able to compensate for the fact that pumping, due to the physical properties and geological characteristics of unconventional fuels, becomes substantially more energy consuming and slower.63,210–212 In this work, the updated forecasts produced by J. Laherrère,105 a senior geologist which has been analyzing the depletion of oil and gas for decades and whose estimates have been pretty consistent over time52,62 are applied. For coal, usually seen as a vast abundant resource, recent studies are pointing to potentially large overestimates in coal resource assessments as geologists uncover restrictions on the coal that is extractable.30,213 For the sake of simplicity, in this work, the highest estimate found in the peer-review literature is considered (high case from Mohr et al.,63). With relation to uranium, there is only one group of researchers estimating depletion curves of uranium using data regularly published by the Nuclear Energy Agency,92,95,104 so the most recent estimate from EWG92 is taken. These depletion curves are converted to maximum extraction curves following the method described in Section 2.2.3.
The techno-ecological potential of RES is also a controversial subject in the literature. The dominant view, which is reflected in most IAMs, considers renewables as a huge, abundant source of energy globally, the technological limits being generally assumed to be unreachable for decades, while the concern focuses on the economic, political or ecological constraints.35,37,143,214,215 This view is, however, challenged by an increasing body of literature which highlights the fact that the large-scale deployment of RES needed to replace fossil fuels in the current socioeconomic system faces serious challenges in relation to a diversity of factors: their integration in the energy mix due to their intermittency, seasonality and uneven spatial distribution;110,216,217 their lower energy density – which also exacerbates the global competition for land;36,38,107,108,110,218,219 their dependence on minerals and materials for the construction of power plants and related infrastructures that pose similar problems than non-renewable energy resources depletion;34,39,40,129 and finally their associated environmental impacts.220–225 All together, these factors tend to significantly reduce the techno-sustainable potential of RES.33–35,37,52,64,110,149,150,217,226,227 In this work, the techno-sustainable potentials of RES for the generation of electricity, heat and biofuels, considering the aforementioned factors estimated in Capellán-Pérez et al.,52 are applied.
Hence, there is a debate in the literature if energy resources could power the increasing demand of Humanity during the 21st Century. The above literature review shows that most models as well as institutional organizations lie in the first group, while the second view is in clear minority. Both visions are in fact related with different paradigms with relation to the assessment of scarcity of natural resources, i.e. relative228,229versus absolute87,230 scarcity.
With relation to climate change damages, models focusing on cost-effective policies generally omit the modeling of climate damages for the sake of simplicity,8,181 while those models including them as highly aggregated damages in cost-benefit IAMs have been shown to severely underestimate them.8,17–21,182,185,231,232 Here, we apply a damage function which has been calibrated assuming that, when the global average surface temperature change reaches +1.75 °C over pre-industrial levels, the climate damages on GDPpc would offset the expected GDPpc growth in that year.186 The calibration is justified by the fact that increasing scientific evidence is showing that a temperature increase of 1.5/2 °C may be dangerous for human societies,4,179,180,233 given that the current socio-economic paradigm requires a growing economic system to be functional.234,235
Hence, the 3 variants of the BAU reference scenario simulated in this paper explore the implications for the future of taking 2 radically different assumptions with relation to energy availability and climate change:
• Ref_noER is a world of unlimited energy resources to 2100 but subject to climate change damages as GHG cumulate in the atmosphere.
• Ref_noCC is a world of limited energy resources in which climate change does not impact human societies. This case is conceptually similar to other BAU scenarios performed with biophysical energy models such as STER, ECCO or GEMBA42,43,236 which do not consider constraints such as climate change damages or mineral availability.
• Ref_noER_noCC refers to a world of unlimited energy resources to 2100 and in which climate change does not impact human societies. Generally, these are the most common assumptions followed by the IAMs in the literature.
Table 4 depicts the description of the most relevant inputs and assumptions which characterize the BAU narrative in this work (Ref scenario) (see Appendix A for details). Three alternative scenario cases are implemented, relaxing the respective constraints. The simulation timeframe is set to 2100 as is customary when assessing climate change.
Scenario inputs & assumptions | BAU scenario | |
---|---|---|
a Limited to a maximum of +10% per year, given that very high exponential growth of early technologies cannot be maintained over time as the technology enters the mature phase. b Higher than historic trends (1–2% per year205,207), given that we assume that, even in a BAU context, to cope with a higher share of variable RES will require the promotion of large scale storage options such as PHS. | ||
Population (2015–2100) | SSP2 (+0.3% per year average over the period)201,202 | |
Expected GDPpc growth (2015–2100) | +2% per year (own estimation, see text) | |
Target labor share | 52% in 2050, constant thereafter (own estimation based on WIOD60,203) | |
A matrix (2100) | Global A matrix in last year available (2009) (own computation from WIOD60,203) | |
Efficiency improvements (final energy intensity) | Historical efficiency improvement trends by sector/households and final energy84 | |
Inland and households transport | ||
Oil dependence target per aggregated category in 2100: (continuation of oil dependence although at a decreasing share. Own estimations from 52,160,161). | ||
4-Wheel vehicles | 34% | |
2-Wheel vehicles | 0% | |
Heavy vehicles | 74% | |
Bus | 54% | |
Recycling rates of minerals (19 minerals) | Current recycling rates237 | |
Nuclear capacity | Constant current capacity208,209 | |
*Non-renewable energies depletion curves | ||
Oil | Laherrère (2018)105 | |
Natural gas | Laherrère (2018)105 | |
Coal | Mohr et al., (2015) High case63 | |
Uranium | EWG (2013)92 | |
GHG emissions from other gases than CO2 and CH4 from fossil fuel combustion | RCP8.5 pathway238 | |
*Climate change impacts | Damage function calibrated to +1.75 °C as “dangerous climate change level”186 | |
Renewables for electricity generation | ||
Technology | Annual capacity growth (2015–2100) (annual historic short-term averaged growth (2012–2015) 79,205,207)a | *Techno-sustainable potential 52 |
Hydroelectric | Annual historic short-term averaged growth (+3.8% per year) | 1 TWe |
Geothermal | Annual historic short-term averaged growth (+4.2% per year) | 0.3 TWth |
Bioenergy | Annual historic short-term averaged growth (+7.8% per year) | Shared potential for heat, liquids and electricity (30 EJ per year) |
Oceanic (wave, tidal and OTEC) | +10% per year | 0.05 TWe |
Onshore wind | +10% per year | 1 TWe |
Offshore wind | +10% per year | 0.25 TWe |
Solar PV | +10% per year | 100 Mha |
Solar CSP | +10% per year | |
Pumped Hydro Storage (PHS) | +5% per yearb | 0.25 TWe |
Renewables for heat generation | ||
Technology | Annual capacity growth (2015–2100) (commercial//non-commercial) (annual historic short-term averaged growth (2011–2014) 79,206,239)a | *Techno-sustainable potential 52 |
Solar thermal | +10% per year//+10% per year | Endogenous depending on urban land |
Geothermal | +5.2% per year//+7.7% per year | 4.4 TWth |
Modern solid bioenergy | +5.8% per year//+11.8% per year | Shared potential for heat and electricity (30 EJ per year) |
Bioenergy | ||
Type | Annual historic short-term averaged growth (2012–2015) 240 and IEA ETP (2017)204 | *Techno-sustainable potential 52 |
Conventional bioenergy | — | 30 EJ |
2nd generation cropland | +3.5% per year | 100 Mha |
3rd generation cropland (starting 2025) | 6% per year | |
Residues (starting 2025) | 6% per year | 25 EJ per year |
Marginal lands (starting 2025) | 6% per year | 386 Mha (∼4 EJ per year)106 |
Fig. 8 Final energy consumption per capita, gross domestic product per capita and system efficiency 1995–2100 for the four scenario cases. (A) Total final energy per capita (TFECpc, GJ per person per year). (B) Gross Domestic Product per capita (GDPpc, thousand $ per person, in USD chained linked volumes (1995)). (C) Total final energy consumption (TFEC) intensity (defined as TFEC/GDP, in J per $). (D) Total primary energy supply (TPES) intensity (defined as TPES/GDP, in J per $). (E) Primary to final energy ratio (defined as TPES/TFEC, dimensionless). (F) Dynamic energy return on energy invested (EROIst) of the system (dimensionless ratio). Representation of system EROI thresholds to sustain high levels of development in current industrial complex societies from Brandt154 (∼5:1) and Hall et al.,135 (∼10:1). Dollars correspond to 1995 US$. Comparison with the IPCC-AR5 range of baseline scenarios for 21008 in panels (B and D) median and 50%, 90% and 100% confidence intervals (values have been converted to MEDEAS equivalent units). Source: own elaboration. |
TFEC intensity levels decline in all scenario cases until mid-century. However, the high penetration of RES prevents further improvements of TFEC intensity beyond ∼2050: the material and energy investments required for the transition to RES imply a re-materialization of the economy that compensates for the efficiency improvements (Fig. 8C). On the other hand, the transition to RES allows for a continuous reduction of the TPES intensity, for all scenario cases, through the timeframe of the analysis. This is, however, not so fast as typically assumed in other models’ BAUs (see Fig. 8D).
Despite both the final (TFEC) and primary (TPES) energy intensities depicting a decreasing trend through the simulation period, the system efficiency, as measured by other indicators, does not improve. For example, the primary to final energy ratio (TPES/TFEC, see Fig. 8E) by 2100 is similar to current levels (1.5) for the scenario cases Ref_noER_noCC, Ref_noER and Ref_noCC, and only in Ref does the ratio improves to 1.25. Likewise, the EROI of the whole energy system also falls from current 12:1 levels to between 6 and 8:1 by 2100 (Fig. 8F). The latter range translates into an over-demand of 5–12% to supply the same level of net energy to society. An inverse relationship exists between the installation rate of new RES capacities and the EROI of the system, due to the up-front energy investments associated to the new RES capacities, which is more accentuated when the potential of those RES technologies with higher EROI, such as wind or hydro, is in full operation. Note that during a fast transition to RES, the EROI of the full system can be temporarily well below the weighted average of the static EROI of the technologies and their supporting systems. This happens in the Ref_noER case, where the EROI reaches a minimum by 2055 (see Fig. 8F) coinciding with the fastest RES growth in the system (see Fig. 7B), and recovering thereafter.46
It is noteworthy that in Fig. 7, 8A and B, the results for the Ref case differ substantially from the other three cases (Ref_noCC, Ref_noER, and Ref), while the results for all four cases are roughly similar in Fig. 8C–F. This is due to the fact that TPES, TFECpc and GDPpc are extensive variables which depend on the size of the human system, while the rest of variables (share of RES, intensities and EROI) are intensive variables characterizing the system which do not depend on its size. As a result, given that in the Ref_noER_noCC case the human sphere grows relentless within the biosphere, the extensive variables increase very fast. However, this variation is much smaller for the intensive variables since we are assuming a same common scenario BAU as a benchmark for the 4 simulated cases.
Fig. 9 Environmental impacts 1995–2100. (A) Total RES land requirements (million ha, Mha), as comparison, the current urban and arable land is represented;241 (B) CO2eq intensity of primary energy (tCO2e per J); (C) GHG emissions (GtCO2e per year); (D) temperature change over pre-industrial levels (°C). Comparison with the IPCC-AR5 range of baseline scenarios for 21008 in panels B, C and D: median and 50%, 90% and 100% confidence intervals (values have been converted to MEDEAS equivalent units). Source: own elaboration. |
The energy transition will require additional mineral requirements, besides the ones demanded by the rest of the economy. Table 5 shows that, despite the different levels of deployment of RES in the 4 simulations studied, some minerals are expected to undergo supply constraints in the future: copper, indium, manganese, molybdenum, nickel, silver, tin or tellurim.
When activating climate change damages, those increase progressively as GHGs build up in the atmosphere, ultimately destabilizing global climate, as seen in Section 3.2. Consequently, climate damages cancel out the expected GDPpc growth by the mid-century, causing a global recession afterwards in both scenario cases, reaching −3 and −3.5% GDPpc per year by 2100 for the Ref and Ref_noRE, respectively.
How do both limits to growth interact in the Ref scenario case? Energy scarcity is the most relevant constraint in this case from the 2020s to the early 2040s. During these years, as mentioned above, liquids scarcity hampers the expected GDPpc; renewables and efficiency improvements do not to fill the gap. Still, as seen in previous sections, the global economy continues to function and burn fossil fuels. The warming in the pipeline from GHG emitted in previous decades and the continued accumulation of GHGs in the atmosphere driving temperature increase, even in a likely context of crisis, makes that in the 2040s, climate change damages outpace energy scarcity as the main limit to growth, which ultimately stops having an effect in the 2050s. Climate damages exacerbate the tensions that harm the economy to such a level that, by then, energy availability is no longer an issue. In the late 2050s, climate change damages themselves hamper the expected GDPpc growth, and by 2100, climate damages cause GDPpc losses of ∼5% per year with relation to the expected GDPpc growth, i.e., GDPpc is decreasing at −3% GDPpc per year. Interestingly, although in the Ref scenario the GDPpc is not able to grow during most of the century, the rate of GDPpc reduction in the second half of the 21st Century is smaller than for the Ref_noER scenario case. This is due to the fact that there is a trade-off between fossil energy availability and climate change: if more fossil energy is available and burned by the economy, on the one hand it makes the economy work, but on the other it increases the level of GHGs in the atmosphere, worsening climate change, whose impacts affect the system with a delay due to the inertia of the climate.
GDPpc and TFECpc decrease from mid-century onwards, reaching levels below current requirements to cover basic needs in the current dominant socio-economic system.188,190 In the current context of globalization and growth-oriented economies,234,235 such a long period of crisis or recession would destabilize human societies, most likely driving them towards different socio-political regimes and thus altering the global geopolitical order. A BAU future where most countries maintain the growth-imperative in a context of likely biophysical scarcities (energy, minerals, land, etc.) would most likely boost conflict over the remaining available resources, deriving in a ‘regionalization scenario’ as identified by van Vuuren et al.:199
“Scenarios in this family assume that regions will focus more on their more self-reliance [sic], national sovereignty and regional identity, leading to diversity but also to tensions among regions and/or cultures. Countries are concerned with security and protection, emphasizing primarily regional markets and paying little attention to common goods. […] Among the more extreme scenarios are the conflict scenarios identified by the GSG [Global Scenario Group] (barbarization). A key issue in these scenarios is how much self-reliance is possible without becoming harmfully ineffective with respect to supranational issues of resource depletion and environmental degradation. The category, like others, includes different variants – but these are rarely optimistic with respect to global poverty reduction and environmental protection.”199
Hence, a widespread systemic global socioeconomic and environmental crisis is foreseen in the next few decades in the absence of fast and drastic global sustainability policies. This is in line with other assessments (e.g., ref. 4, 179 and 242–244). However, the novelty is the consistency of the outputs of the IAM with these assessments. It should also be acknowledged that the pressures on the system would likely be too strong for it to remain the same. Hence, at some point after the curves start to fall, the projections and the model lose their prediction capacity.
Still, two main differences with the BAUs from the literature are relevant: a higher TPES (which corresponds with higher TPES intensity levels) (Fig. 7A and 8D) and a lower CO2e intensity of primary energy (Fig. 9B). Despite a decline in TPES intensity in the 4 scenario cases over the simulated period, the reduction is found to be smaller than for other BAUs (22–27% vs. 65% for the median of the IPCC-AR5 baseline scenarios8). This can be explained by three features explicitly considered in the MEDEAS framework: the energy investments related to the deployment of electric renewables,44,46 biophysical and thermodynamic limits in the substitution of inputs in production in the medium and long-term,83,84 and the strong link between economic activity and the use of natural resources, given the difficulties to dematerialize economic production when accounting for sectors’ complementarity.72,245,246 A lower CO2e intensity of primary energy sources levels is obtained due to a higher penetration of RES in MEDEAS BAU scenario cases. This has more to do with scenario specification than modeling assumptions, given that the continuation of current high growth RES is assumed in the next few decades, unlike in other IAMs’ BAUs (e.g., ref. 8 and 247).
The obtained results for the Ref_noCC case for the energy outputs are broadly similar to the baselines from biophysical energy models such as STER, CORECCO or GEMBA,42,43,236 which project a decline in energy availability by the mid-21st Century. These models do not have a distinct economy module, and have a very simplified and stylized way to represent energy demand (typically through aggregated, sometimes biophysical, production functions), while the strength is on the representation of supply, and more specifically on the potential technical and biophysical constraints of supply to follow expected demand. Mineral constraints and climate change impacts are also not represented. The results from the Ref scenario could be conceptually compared with the baseline scenario of World3 (“standard run”),14,248 given that this model accounts for limited (stock) non-renewable energy availability and impacts through persistent pollution, which delivered an “overshoot and collapse” scenario in the first half of the 21st Century. It is noteworthy that this scenario has performed very well when comparing with 40 years of historical data.249,250
The transition to renewable energy sources will require large amounts of interrelated energetic, material and economic investments in order to build, operate and interconnect the new energy infrastructures. This finding is consistent with previous results from biophysical energy models such as STER, ECCO or GEMBA.42,43,236 However, most policy-relevant energy models and IAMs focus solely on monetary investments (e.g., IEA, IPCC, national governments, etc.). The modeling of the material and energy investments for the energy transition in MEDEAS allows us to dynamically assess potential mineral scarcities affecting most “green” technologies (with findings in the line of other works in the literature39,40,46,129) and to compute the net energy available to society. This approach reveals that there is a trade-off between urgent climate mitigation and the viability of the system, i.e., faster deployment rates of alternative energy systems tend to reduce the EROI of the system, highlighting the need to complement classical monetary costs (e.g., ref. 8 and 251–256) with biophysical quality indicators to analyze energy transitions.43,257,258 Capellán-Pérez et al.,46 showed that a rapid transition to RES would imply a re-materialization of the economy with the potential to counteract the historic trends of efficiency improvement.
The developed method also allows us to detect potential harmful situations of increasing gross energy output, while decreasing the net energy delivered to society, i.e., the so-called “energy trap”.45,138,259 In extreme cases, a too low EROI, even if the gross energy consumption is increasing, may even trigger a collapse of the full system, although this behavior is not captured in the current framework. In this sense, the EROI levels obtained by 2100 are between the thresholds identified in the literature in order to maintain a complex industrial society (5–10:1135,154).
Mainstream macroeconomic models – also most economic modules in IAMs – are typically based on computable general equilibrium forcing markets clearing, disregard of the sectoral structure of the economy (including the general assumption that energy can be readily substituted by manufactured capital), and optimization.28 However, the general assumption of perfect substitutability between factors has been widely criticized from an ecological economics perspective, which considers that complementarity better fits with reality.23,71,260–262 In fact, aggregated production functions fail to capture the relevance of economic structure in energy-environment interactions.22,263,264 Scholars have also made critical remarks concerning optimization as an unrealistic approach for modeling complex, dynamic systems where feedbacks and time matter.28,64,265 In MEDEAS, productive capacity is determined by demand and energy availability. Therefore, neither equilibrium nor optimization are imposed. Moreover, a more realistic approach to energy substitution is enabled by a rich modeling of the technological evolution rates and potentials of energy sources. By using an input–output framework, complementarity between productive factors is assumed, which allows assessing the direct and indirect effects in sectoral production, given an economic structure and the evolution of demand.73,266 In addition, IOA allows environmentally-extended hybrid approaches,73,267 which facilitates a much more accurate energy demand estimation, since it is conducted based on all sectors’ direct and indirect requirements by a vector of as many final energy sources as the modeler desires – in our case 5 types (see Section 2.2.3). Hence, the spillovers between sectors in IOA change the sectoral composition of energy demand. Because the energy requirements of each sector are different, both qualitatively – distinct energy sources – and quantitatively, economic structure can be of the utmost importance.
There has historically been a large discrepancy between natural scientists’ understanding of climatic damage and their representation in climate-focused IAMs. Despite the disruptive potential of future climate change, which threatens human societies as we know them nowadays through severely damaging our natural life-support systems, most IAMs either omit or underestimate climate damages. The most influential cost-benefit IAMs including climate damages, such as DICE, PAGE or FUND, calibrate their damage functions in their standard configuration to monetized estimates of damages which suggest that a century of climate change is about as bad for welfare as a year of standard economic growth (1-4% of Gross Domestic Product (GDP) loss).17,268 These discrepancies inevitably produce confusion in policy-makers and the citizens, as they receive contradictory messages from the scientific community in relation to the potential damages caused by climate change. Given the existing uncertainties and limitations to accurately represent the human–nature systems over the long-term, as well as the need for urgent action, we judge that the only honest and achievable objective of modeling is to strategically provide advice about the nature of the policy decisions to be made in the next few decades. In this context, MEDEAS proposes a novel, simple and transparent methodology to consistently integrate climate damages into the modeling framework with climate assessments by natural scientists.186 When surpassing 2 °C of global mean temperature increase, a collapse of the system is more likely than reaching 4 or 5 °C without affecting the socioeconomic dimension. However, it should also be borne in mind that reducing all climate impacts to a monetary value presents significant limitations, given that many impacts also affect non-market dimensions.8,17
Future work will be focused on the implications derived from the MEDEAS modeling approach, applying regional models both for BAU and policy-intervention scenarios.
As with any other model, MEDEAS models are not intended to precisely predict the future, but rather to understand system behavior and provide policy-guidance to the best options for the energy transition towards a sustainable economy. They are tools to explore strategies, informing qualitative insights through quantitative modeling. Likewise, the MEDEAS framework also presents a number of limitations and uncertainties such as one-way model integration, low quality data (e.g., IOT sectoral granularity and length of time series, IOT associated socioeconomic and environmental accounts, data on minerals’ availability and consumption, etc.), low development of some modules (e.g., land-use, water), etc. Note that the fix of some of the most important limitations would tend to worsen the obtained results, such as integrating the feed-back of mineral scarcity, the computation of the EROIext of the system instead of EROIst, capturing the metabolic implications of the drop of the EROI of the system to very low levels, etc. Further work is going directed¶¶¶¶ towards more dynamization in the Economy module, the creation of new modules to enable the endogenization of more dimensions (e.g., Population), deeper integration between modules (e.g., coherence between energy and monetary investments, endogenous change of the matrix A driven by energy mix and technological change, integration of land-use and water, etc.), expanding the representation of potential limits to growth (mineral availability, investments for the energy transition, labor supply, etc.), refining the modeling of climate change damages or the modeling of behavioral policies.
This paper describes the open-source MEDEAS integrated assessment modeling framework, which has been developed with the aim of informing decision-making to achieve the transition to sustainable energy systems with a focus on biophysical, economic, social and technological restrictions and applying alternative modeling assumptions and methods to a core set of common assumptions in IAMs whose validity is being disputed in the scientific discussion: lack of feedback, deficient representation of economic structure, dominance of equilibrium-optimization approaches, the widespread use of prices as indicators of scarcity, energy abundance presumption and neglect of the material and energy investments related to the transition to renewables.
The models are developed in system dynamics, which facilitates the integration of knowledge from different perspectives and disciplines, as well as feedbacks from the different subsystems. MEDEAS models are designed to run until 2060–2100, and are conceptually structured into nine modules: economy, energy demand, energy availability, energy infrastructures and energy return on energy invested (EROI), minerals, land-use, water, climate/emissions, and social and environmental impact indicators. The modules of Economy and Energy are the most detailed. The models include six main features which all together represent a step forward in the state-of-the-art of the field:
• Representation of biophysical constraints to energy availability
• Modeling of the material and energy investments for the energy transition, allowing to dynamically assess potential mineral scarcities and the net energy available to society
• Consistent representation of climate change damages with climate assessments by natural scientists
• Integrated representation of economic processes and biophysical limits to growth, including the integration of detailed sectoral economic structure (input-output analysis) within a system dynamics approach
• Energy shifts driven by physical scarcity
• Rich set of socioeconomic and environmental impact indicators
Hence, the MEDEAS framework incorporates two limits to growth that are rarely considered (even separately) in the literature: energy availability and climate change damages. This is a first step, given that human societies operate within many more biophysical constraints.1,169
The simulation of 4 business-as-usual (BAU) scenario cases, switching on/off these two limits to growth in the MEDEAS-World model illustrates the potentiality of the framework. The consideration of alternative plausible assumptions and methods, combined with the framework's feedback-rich structure, display dynamics absent from the classical models. In the absence of energy availability constraints and climate change damages, the results obtained throughout the 21st Century are broadly similar to those obtained by other IAMs’ BAUs in the literature.8 However, when activating any of these constraints, the results get completely modified. Unlike the current state-of-the-art, BAU ceases to be a scenario with low renewable penetration in the energy mix, Gross Domestic Product per capita (GDPpc) 4 to 8 times greater than today's and a 3.5–4.5 °C temperature increase by 2100. The MEDEAS-W results show in the same timeframe, instead, a BAU with a large penetration of renewables in the energy mix (60–80%) driving large requirements of minerals, energy investments and land, a persistent recession over the next decades with global-average GDPpc levels similar or below current levels caused by a combination of energy scarcity and climate damages, which prevent us from emitting GHGs at increasing rates. This contributes to keeping temperature increase levels <2.5 °C by 2100 (omitting the possibility of the occurrence of a fast tipping point in the climate system which would worsen the obtained results). Hence, a widespread systemic global socioeconomic and environmental crisis is foreseen in the next few decades in the absence of fast and drastic global sustainability policies. From the 6 ‘scenario families’ identified by van Vuuren et al.,199 in global environmental assessments, our results suggest that current trends increasingly correspond to the ‘regional competition’ storyline, where a focus on national sovereignty would give priority to security concerns and trade barriers, which may derive in conflict and even collapse scenarios. Despite depicting a much more worrying future than conventional BAUs, we do however believe it is a more realistic counterfactual scenario that will enable the design of improved alternative sustainable pathways in future work.
BAU | Business-as-usual. Scenario projecting future changes as the continuation of current trends |
C-ROADS | Climate rapid overview and decision support. https://www.climateinteractive.org/tools/c-roads/ |
CSP | Concentrated solar power |
Dmnl | Dimensionless |
EROI | Energy return on energy investment. Ratio between the energy delivered from a process divided by the energy required to get it over its lifetime |
ESOI | Energy stored on energy invested. Ratio between the energy stored in a storage device divided by the energy required to get it over its lifetime |
EU | European Union. EU-28 |
EV | Electric vehicle |
GCAM | Global change assessment model. http://www.globalchange.umd.edu/gcam/ |
GDP | Gross domestic product. Monetary value of final goods and services – that are bought by the final users – produced in a region in a given period of time |
GFCF | Gross fixed capital formation. Investments made by each component of the final demand in durable products produced by each sector |
GHG | Greenhouse gas. Gas that absorbs and emits radiant energy within the thermal infrared range |
HH | Households |
IAM | Integrated assessment model |
IEA | International energy agency |
IMAGE | Integrated model to assess the global environment. https://models.pbl.nl/image/ |
IOA | Input–output analysis. Quantitative economic model that represents the interdependencies between different sectors of a national economy or different regional economies |
IOT | Input–output tables |
IPCC | Intergovernmental panel on climate change. Intergovernmental body of the United Nations that is dedicated to providing the world with an objective, scientific view of climate change, its natural, political, and economic impacts and risks, and possible response options |
LCA | Life-cycle analysis. Technique to assess the required resource inputs (e.g. energy, materials, water, etc.) or related environmental impacts associated (e.g., GHG emissions) with all the stages of a product's life from raw material extraction through materials processing, manufacture, distribution, use, repair and maintenance, and disposal or recycling |
MEDEAS | Modeling energy system development under environmental and socioeconomic constraints. Name of the H2020 project (https://www.medeas.eu/) and by extension of the models developed in that framework |
NRE | Non-renewable energy |
pc | per capita |
PHS | Pumped hydro storage |
PV | Photovoltaic |
RCP | Radiative concentration pathway. Set of scenarios developed as a standard basis for near and long-term climate modeling experiments (http://tntcat.iiasa.ac.at:8787/RcpDb/) |
Ref | Scenario case used in this paper. With restrictions to energy availability; with climate change damages |
Ref_noCC | Scenario case used in this paper. With restrictions to energy availability; no climate change damages |
Ref_noER | Scenario case used in this paper. No restrictions to energy availability; with climate change damages |
Ref_noER_noCC | Scenario case used in this paper. No restrictions to energy availability; no climate change damages |
RES | Renewable energy source |
SDG | Sustainable development goal. Collection 17 global goals set in 2015 by the United Nations General Assembly and intended to be achieved by the year 2030 |
SSP | Shared socioeconomic pathway. Set of scenarios developed as a standard basis for projected socioeconomic global changes up to 2100 (https://tntcat.iiasa.ac.at/SspDb/) |
TFEC | Total final energy consumption. Final energy use from the side of demand |
TFES | Total final energy supply. Final energy use from the side of supply |
TPES | Total primary energy supply. Total final + generation + transport, distribution and transmission losses |
WIOD | World input–output database. http://www.wiod.org/ |
For population, the “middle of the road” of the Shared Socioeconomic Pathways (SSP2) was considered, given that this scenario is “intended to represent a future in which development trends are not extreme in either of the dimensions but rather follow middle-of-the-road pathways relative to the span of plausible outcomes for each element”.198 Combining medium fertility with medium mortality, migration and education scenarios for all countries, the SSP2 assumes a world of increasing population, but at a rate that was slowing down, reaching 9.2 billion people in 2050 and 9.0 in 2100 from the current 7.5.201,202
The expected gross domestic product per capita (GDPpc) is a key variable of the model, given the connection between economic activities and resource consumption. We select an intermediate value between historic global trends of +1.42% per year (1979–2014) – a period characterized by lower growth trends than those in the post-war “The Glorious Thirty” (1945 to 1975)72 – and the +2.5% average growth stated by the SSP2, i.e., +2% per year. This means that GDPpc would be expected to reach over 30000 US 1995 $ by 2100 from a current level of ∼7000 $.
In terms of income distribution, the historical decline of labor share, i.e., the part of national income allocated to wages, is extrapolated into the future, but at a reduced pace due to the increasing importance of emergent countries in the global GDP. As these countries, in their respective modernization processes, are improving their labor shares, it is to be expected that they partially compensate for the global decline trend in labor share in the future. Although China is following the opposite trend, this is due to its relatively higher current labor share. Hence, a value of 52% for the world average is assumed for 2050 from the current 56%, which is then kept constant until the end of the century. It is worth saying that scenarios with a capital share higher than the labor share are found in countries with a low developed welfare state and high rates of inequality indexes,269,270 while current EU and well developed welfare state countries depict values close to 60%.78
Given the complexities of consistently computing the dynamic behavior of the matrix A, the current global economic structure is assumed to remain constant over the next few decades in the context of a BAU scenario. This is supported by the relatively minor variations to the technical coefficients in the available historical data.60,203
Technological improvements are modeled top-down at sectoral level (except for transportation where it is done bottom-up), extrapolating the historical efficiency improvement trends by sector/households and final energies into the future.52,84
In inland and households transport, the continuation of the dependence on oil is assumed, although its share is reduced over time via the shift to electric, hybrid and gas vehicles. More optimistic prospects are expected for light vehicles than for heavy ones, based on current data and analyses.52,160,161 For the case of air and water transport, we consider that, within the BAU storyline, there would not be a substantial replacement of conventional fuels by low carbon alternatives, given the technological challenges involved which remain unresolved.155
Nuclear installed capacity and electricity generation are constant for over two decades in the absence of major new-build programs apart from China. In this context of increased ageing of existing reactors, increasing costs and general decline in interest in nuclear new-build, the continuation of current trends may drive global nuclear capacity downwards over the next few decades.208 This trend may be compensated for by eventual lifetime extension programs.209 For this reason, in the BAU, we consider that the nuclear capacity in operation globally will be constant. This assumption corresponds well with the projection of the international news agency for investors Bloomberg, which estimates, in its “New Energy Outlook 2018”, that by the mid-century nuclear will produce a level of electricity similar to current levels.209
The installation of RES technology capacity is growing globally, with large variations depending on the technology and the period of time analyzed. Particular conditions, such as the economic recession from the year 2008, have influenced their development pace. However, the fact that, when a technology starts to deploy, each additional power plant represents a high share in relation to the cumulative installed capacity should also be taken into account. This translates into explosive initial exponential deployment rates, which are softened over time as cumulated capacity increases. Annual historic short-term averaged growth (2012–2015) has generally been taken as the reference for electric,79,205,207 heat79,206,239 and bioliquid technologies.204,240 These rates are limited to a maximum of +10% per year, given that very high exponential growth of early technologies cannot be maintained over time as the technology enters the mature phase. Given the environmental impacts of conventional biofuels,33,271–274 it is likely that, in the future (BAU), current growth rates of this technology will not be maintained, and for these reasons, their growth is halved with relation to the historic short-term averaged growth.204
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9ee02627d |
‡ MEDEAS project webpage: https://www.medeas.eu/. |
§ MEDEAS models have been developed in system dynamics software Vensim DSS for Windows Version 6.4E (×32) and thereafter translated to Python. The latest versions of both Vensim and Python codes are freely available at: https://www.medeas.eu/model/medeas-model. Future updates of the models will be available at: http://geeds.eu/. |
¶ Boldface has been used to note matrices according to matrix notation. Nevertheless, subsequent equations are referenced to their index components. |
|| For a full description and technical parameters of the Economy module in MEDEAS-W see Nieto et al.72 |
** For a full description and technical parameters of the energy intensities modeling in MEDEAS see de Blas et al.84 |
†† For a full description and technical parameters and technology choice of the energy availability modeling in MEDEAS see Capellán-Pérez et al.52 |
‡‡ The selected representative technologies are the following: CSP with molten-salt storage without back-up, most efficient and used technology (back-up option is usually powered by non-renewable fuels such as natural gas);123 fixed-tilt silicon PV (best EROI)124 and subject to fewer mineral availability constraints (a weighted average is computed for some minerals taking into account the current share of PV sub-technologies such as thin-films);34 2 MW onshore wind turbines (currently the global average onshore wind turbine capacity is lower than 2 MW125); 3.6 MW offshore wind turbines taking as reference the current average size in Europe;125 and LiMn2O4 electric vehicle batteries given that, although they are less efficient than other alternatives (e.g. LiCoO2), the embodied energy for their fabrication is substantially lower.126 |
§§ For more details on the methodology and assumptions for the estimation of materials in MEDEAS, see Capellán-Pérez et al.46 |
¶¶ Note that dynamically accounting for energy magnitudes corresponds with power rather than to energy. Despite our dynamic approach significantly shortens the time step of the calculations (in the order of months), it represents still an average power over a certain time (which in conventional EROI studies corresponds with the lifetime of the technology). Hence, we decided to avoid the creation of a new term given that EROI is a concept widely used, and follow the terminology of “dynamic EROI” more commonly used in the literature75,138,139 (although other more recent works have coined new terms such as “power return on energy invested” (PROI)136,140 or “net external power ratio” (NEPR)141,142). Still, it should be kept clear that the different nomenclature refers ultimately to the same concept. |
|||| Note that there is a lack of standardization in the literature, and other works estimate the EROIst at the power plant instead rather than at the system level, use the energy delivered at the mouth-gate of the power plants rather than the one delivered at the final consumer, and do not incorporate the additional energy investments related with intermittency management. |
*** For more details on the methodology and assumptions for the estimation of the EROI of the system in MEDEAS, see Capellán-Pérez et al.46 |
††† For a full description and technical parameters of the modeling of transportation in MEDEAS, see Capellán-Pérez et al.52 |
‡‡‡ For more details on the land-use module in MEDEAS, see Capellán-Pérez et al.,52 and de Blas Sanz et al.53 |
§§§ Blue water refers to the consumption of surface and ground water; green water is the volume of rainwater consumed, mainly in crop production; and gray water is the volume of freshwater required to assimilate the load of pollutants based on existing ambient water quality standards. |
¶¶¶ The authors are aware that IAMs which do not include a damage function usually operate under the “cost-effectiveness” approach, i.e., they focus on finding a system configuration which minimizes the costs of the transition.8,181 However, by ignoring the costs of the non-action (already partially unavoidable), the results obtained under this method are, therefore, flawed. As an example, the investments required for the transition obtained in these IAMs are computed and communicated as mitigation costs. |
|||||| For a full description and technical parameters of the climate modeling in MEDEAS, see Capellán-Pérez et al.,52 and for the modeling of climate change damages in MEDEAS-W, see Capellán-Pérez and de Castro.186 |
**** Even though MEDEAS only estimates 1 out of 3 of the components of the HDI (neither life expectancy at birth, nor adult literacy or school enrolment are modeled), an alternative method to estimate the potential HDI that can be reached by a society has been developed. Assuming that the quality of life has a material dimension (minimum energy requirements), data for the final energy footprint of 40 countries for the timeframe 1995–2009 have been used to estimate a regression between the potential HDI levels and the energy use per person. |
†††† For more details on the social and environmental impact indicators, see Capellán-Pérez et al.,52 and de Blas Sanz et al.53 |
‡‡‡‡ This is different from projecting static current values into the future. For example, the dependence of transportation on oil is currently very high. However, given the current trends of expansion of alternative vehicles, it seems reasonable to assume that this dependence will decrease over the next few decades. |
§§§§ The variation of EROI at system level, as currently modeled, does not act as a separate limit to growth, given that in the absence of the other two constraints (as in the Ref_noER_noCC scenario case), the expected GDPpc is achieved. However, this is achieved at the cost of more energy expenditure and environmental impacts. |
¶¶¶¶ These limitations are being addressed in the ongoing H2020 project ‘LOCOMOTION’ (H2020-LC-CLA-2018-2, Project Number 821105, https://www.locomotion-h2020.eu), approved by the European Commission with the objective of improving the MEDEAS framework. |
This journal is © The Royal Society of Chemistry 2020 |