1 Introduction

There is an urgent need to limit global warming to 2 °C or lower relative to pre-industrial levels. To reach this target, greenhouse gas (GHG) emissions need to be substantially reduced which particularly come from the burning of fossil fuels (IPCC 2014). To add to this challenge, large emerging economies, such as Brazil, China and India, are experiencing great economic and population growth (Alam et al. 2016; Haseeb et al. 2016; Zaman et al. 2016). This leads to a significant increase in consumption of resources, for instance in housing, vehicles and electronics. In the past years, an accelerated shift toward renewable (low-carbon) energy technologies can be observed (International Energy Agency 2015). The large-scale implementation of renewable energy technologies, however, comes with increased use of mineral resources, such as copper for photovoltaic systems and iron for wind power plants (Hertwich et al. 2014; Kleijn et al. 2011; Pehlken et al. 2017).

When comparing different electricity production technologies, trade-off evaluation with life cycle assessment is crucial (Prado-Lopez et al. 2016). To evaluate the trade-offs in using fossil vs. mineral resources to produce energy, surplus cost can be used as an indicator. The surplus cost potential of a resource is a suggested measure of the short-term economic scarcity of that resource assuming that its production costs will increase with continued production. Steen (1999) was the first to propose surplus cost as an indicator to assess the life cycle impacts of products, services and technologies. Despite this development from the 90s and others that followed, currently there is still extensive debate about natural resources in life cycle impact assessment (LCIA) (Dewulf et al. 2015; Drielsma et al. 2016; Lieberei and Gheewala 2017; Sonderegger et al. 2017; Steen and Palander 2016). The European Commission-Joint Research Centre-Institute for Environment and Sustainability (2011) indicated surplus costs as a promising approach to quantify abiotic resource scarcity but has, nevertheless, considered it not mature for recommendation in the evaluation of product life cycle impacts. Drielsma et al. (2016) also point to the potential of using economic scarcity as indicator for short-term availability. The surplus cost potential indicator can be seen as one example of how economic scarcity can be used in LCIA-based resource scarcity assessments. Recently, Vieira et al. (2016) derived surplus cost potentials (SCP) for metals, quantifying the average cost increase expected from future resource extractions. They described the relation between costs and extraction with cost-cumulative tonnage relationships and considered two different reserve estimates, i.e. (economic demonstrated) reserves and ultimate recoverable resource, in their SCP calculations. Ponsioen et al. (2014) derived SCPs for fossil fuels with cost-cumulative tonnage curves as well but focussed on marginal cost changes in the current situation and considered economic discounting in their calculations. Although substantial progress has been made in the further development of the SCP method, for practical application, the number of mineral resources needs to be expanded and the assessment of fossil fuels and mineral resources needs to be aligned (Klinglmair et al. 2014).

The goal of the study was to derive average surplus cost potentials for a large number of mineral and fossil resources. The SCPs were subsequently applied to derive surplus costs for non-renewable (hard coal, natural gas, oil, nuclear) and renewable (wind, solar, hydro, geothermal) electricity production in 56 countries.

2 Methods and data

2.1 Fossil resources

The SCP of fossil resource x was calculated by

$$ {\mathrm{SCP}}_x\left(\mathrm{fossil}\right)=\frac{\int_{{\mathrm{CFE}}_x}^{{\mathrm{MFE}}_x}\ \left({\varDelta C}_x\right)\ d{\mathrm{FE}}_x}{{\mathrm{R}}_x} $$
(1)

where C x is the production cost of fossil resource x (in USD/kg or USD/Nm3 x) determined via a cost-cumulative production relationship (see Eq. (2)), FE x is the future extraction (in kg or Nm3), R x is the total future production (in kg or Nm3), MFE x is the maximum amount available for extraction (in kg or Nm3) and CFE x is the current cumulative tonnage (in kg or Nm3) of fossil resource x extracted. There is a very important assumption of this model, namely that all reserves up to MFE will be extracted. The cumulative fossil resource extracted (CFE), the maximum amount available (MFE) and the actual reserve (MFE − CFE) for crude oil, hard coal and natural gas were taken from the International Energy Agency (2013) (see the Electronic Supplementary Material and Fig. 1).

Fig. 1
figure 1

Visualization of the derivation of the surplus cost potential following an average approach and a log-linear cost-cumulative production curve from the current cumulative fossil extracted (CFE) up to the global maximum fossil available for extraction (MFE)

The cost of fossil resource x was based on the cumulative tonnage extracted of that fossil resource. The relationship between the two variables was derived by fitting a log-linear curve through the cost and cumulative production:

$$ {C}_x=\frac{1}{a_x+{b}_x\bullet \ln \left({\mathrm{CFE}}_x\right)} $$
(2)

where a x and b x are, respectively, the intercept and slope of the log-linear distribution of the cost-cumulative production relationship for fossil resource x.

The data used to derive the cost-cumulative production relationships for crude oil, natural gas and hard coal was retrieved from the International Energy Agency (2013). The cost data reported by the International Energy Agency (2013) was in 1998 US dollars. The CPI Inflation Calculator (Bureau of Labor Statistics 2014) was used to convert the CFs into US dollars for 2013 (USD2013). For crude oil and natural gas, the data of the IEA includes the costs and amounts of fossil resources extracted up to now and expected for the future per production technique, e.g. conventional oil or oil in ultra-deep water. For hard coal, a production-cost curve is provided with free-on-rail (FOR) production costs up to the 2011 estimate of 728 gigatonnes (Gt) of global hard coal reserves (International Energy Agency 2013). The FOR production cost includes the cost of mining and delivering the coal from the mine plus surface handling, coal preparation or beneficiation, storage and loading costs. The International Energy Agency (2013) estimates require assumptions related to the different factors influencing production costs, such as differences in the specific resource category to be exploited, the location of current and future reservoirs and the evolution of capital costs for upstream and downstream production stages. Future operation costs particularly depend on assumptions concerning the development and use of new technologies to facilitate access to more resources (reducing costs) and the increased depletion of reservoirs (increasing costs). We refer the reader to the International Energy Agency (2013) report for more detailed information concerning the assumptions and uncertainties behind the derivation of the global cost curves (pages 227–234 of the report).

A log-linear least-squares fit was applied to the cost-cumulative tonnage data to derive the intercept parameter a and the slope parameter b of the cost-cumulative tonnage relationship for each fossil resource (see SI). The standard error of a and b and the explained variance (R 2) were also reported for each fit. The software R project for statistical computing version 3.2.1 (2015-06-18) (R Development Core Team 2015) was used for this purpose. The step to derive data points from the fossil resource charts from the International Energy Agency (2013) to derive the cost-cumulative production curves is described in detail in the Electronic Supplementary Material.

To express SCPs in USD2013 per gigajoules, we use the higher heating values used in the ecoinvent database (Jungbluth and Frischknecht 2010), namely 45.8, 19.1 and 38.3 MJ/kg or m3 for crude oil, hard coal and natural gas, respectively.

2.2 Mineral resources

Vieira et al. (2016) derived SCP values for 12 metals and platinum-group metals as a group with two reserve estimates, reserves and ultimate recoverable resource. Here, only ‘reserves’, defined as the amount of a mineral resource which can be economically extracted or produced at current prices and state of technology (U.S. Geological Survey 2016), were used. There are, however, mineral resources for which there are no SCPs available because cost and extraction data per mine are lacking. Here, we extended the number of metals and mineral resources covered by deriving a log-linear relationship between surplus cost potentials and prices for metals and mineral resources:

$$ {\mathrm{SCP}}_x\left(\mathrm{mineral}\right)={10}^a\times {\mathrm{price}}^b $$
(3)

where a and b are, respectively, the intercept and slope of the log-linear regression. In commoditized markets, as is the case for mineral and metal resources, price tends to converge toward cost (Crowson 2011; Dixon 2009; Tilton and Lagos 2007).

For the regression analysis, the mineral SCPs of Vieira et al. (2016) were used, while price data for 2013 was retrieved from Kelly and Matos (2013) except for the platinum-group metals and uranium. For palladium, platinum and rhodium, average price data for 2013 was retrieved from Kitco Metals Inc. (2015). The ESA spot U3O8 data (a weighted average of triuranium octoxide prices paid by EU utilities for uranium delivered under spot contracts during the reference year) published by the Euratom Supply Agency (2015) was used to calculate the price for uranium. The SCP and price values of the 12 resources and platinum-group metals as a group that were used as input to construct the regression model are reported in Table 1 and in the Electronic Supplementary Material, respectively.

Table 1 Surplus cost potential (SCP) of fossil and mineral resources, the latter determined for reserves

The regression was applied to 57 additional mineral resources and for the groups garnets, gemstones, rare earths and zirconium mineral concentrates to derive SCP values from price data. The price data can also be found in the Electronic Supplementary Material.

2.3 Electricity production technologies

Data for 10 electricity technologies has been retrieved from ecoinvent v3.2, system model Recycled content (Ecoinvent Centre 2015), as implemented in SimaPro v8.2 (PRé Consultants bv 2015). Because there was no data to derive specific SCP values for peat and lignite, no SCP value for these resources to calculate the surplus costs of the various electricity technologies was included.

Ecoinvent provides life cycle inventory data per energy technology for a total of 56 countries. When various datasets for a specific technology within a country were available, the average surplus costs were calculated. For Canada, China and the USA, datasets per electricity production technology for specific geographical regions within each country are provided. All regions were considered equally important for the country average calculations. As a final step, the average and variation of the surplus costs of a technology between the countries were calculated. More information on the datasets analysed in this study can be found in the Electronic Supplementary Material.

3 Results

3.1 Characterization factors

The regressions between cumulative extraction of a fossil resource and the increase in production costs per unit of fossil resource extracted resulted in explained variances of 78% (for hard coal), 91% (natural gas) and 94% (for crude oil). The regression statistics for each fossil resource are reported in the Electronic Supplementary Material.

Table 1 shows the SCPs for the three fossil resources. The SCP value of crude oil is approximately 5.6 and 1.3 times larger than the SCP values of hard coal and natural gas, respectively.

As shown in Fig. 2, the surplus cost potential has a strong relationship with the price of a metal or mineral resource (explained variance of the regression is 92%). The intercepts of the log-linear relationship between price and SCPs is −0.70. This means that the SCP is typically a factor of 5 (100.70)) smaller than the price of a metal or mineral resource. Table 1 provides a list of SCP values for all resources, with a distinction between those empirically derived here (fossil resources), the mineral resources empirically derived by Vieira et al. (2016) and the mineral resources derived on basis of the price.

Fig. 2
figure 2

Log-linear relationships between the average price in 2013 (USD2013/kg x) and surplus cost potential (USD2013/kg x) derived for reserves (R). The 12 metals and platinum-group metals as a group for which Vieira et al. (2016) derived SCP values were used. Note: both axes are presented in a logarithmic scale

3.2 Surplus costs of electricity generation

In Fig. 3, the surplus cost per megawatt hour of electricity produced is shown for every electricity production technology considered. Per electricity technology, the average over all countries is reported as well as the minimum and maximum surplus cost. Detailed results can be found in the Electronic Supplementary Material.

Fig. 3
figure 3

Surplus cost of different sources of electricity (USD2013/MWh), calculated as the average of all country scores for that technology. The bars in this figure reflect the variability across countries considered showing the minimum and maximum country values

Hydroelectricity from run-of-river and from reservoir has the lowest surplus costs (0.3 to 0.4 USD2013/MWh). With around 2 to 3 orders of magnitude difference, the three fossil-based energy production technologies result in the highest overall surplus costs (23 to 148 USD2013/MWh). Nuclear, geothermal, photovoltaic and wind electricity production technologies result in intermediate results, namely between 1 and 8 USD2013/MWh. Note in this respect that the SCP for uranium obtained in our study is 0.02–0.06 USD2013/GJ, converting from kilogram to gigajoules by its heat value of 500 GJ/kg (World Nuclear Association 2014). This implies that the surplus cost potential of uranium as energy resource is 1.5–2.7 orders of magnitude lower than that of fossil resources.

From Fig. 3, it can also be seen that electricity production from wind shows the highest variability in SCP values, with approximately 2 orders of magnitude difference between the minimum and maximum value. The high variability in SCPs observed for wind may be explained by regional differences in wind intensity. For wind-powered electricity, Russia, Ukraine and South Africa have the largest surplus costs per megawatt hour of electricity produced, whereas Indonesia, Mexico and Turkey obtain the lowest surplus costs.

In Fig. 4, the average relative contribution to the surplus costs per resource for each electricity source is shown. Mineral resources that contribute with at least 5% of the total SCP score per technology are individually shown and the rest is grouped in the series ‘Remaining mineral resources’. The contribution for all three fossil resources is always displayed. We found that the surplus costs of fossil resource use are always higher compared to mineral resources. For the fossil-based electricity generation technologies, the relative contribution to the total surplus costs by fossil resources is always larger than 99%. For the renewable electricity production technologies, mineral resources contribute 8 to 29% to the overall surplus costs. Nuclear electricity production technologies show the largest contribution from mineral resource use (33–34%), mostly resulting from uranium. The mineral resources contributing with at least 5% of the total SCP are aluminium, nickel and zinc.

Fig. 4
figure 4

Relative contribution of fossil and mineral resources to the surplus costs of various sources of electricity production. Mineral resources that contribute with at least 5% of the total SCP score per technology are individually shown and the rest is grouped in the series ‘Remaining’. The contribution for all three fossil resources is always displayed

4 Discussion

4.1 Comparison with other surplus cost methods

Existing LCIA methods focused on quantifying the surplus costs resulting from production of abiotic resources include EPS (Steen 1999), ReCiPe (Goedkoop and De Schryver 2009), Ponsioen et al. (2014) and the economic resource scarcity potential (ESP) (Schneider et al. 2014). Although ESP includes the economic dimension, this method encompasses various elements aimed at identifying altogether supply risks associated with resources, aligning thus with criticality approaches. Because the ESPs per resource were not explicitly provided in Schneider et al. 2014, the quantitative comparison with the approach here proposed was not possible to carry out. In EPS and ReCiPe, characterization factors (CFs) expressed as surplus costs were derived for mineral and fossil resources, whereas Ponsioen et al. (2014) derived SCPs specifically for fossil resources. The SCPs for EPS and ReCiPe were retrieved from Steen (1999), Goedkoop and De Schryver (2009) and Roerbech et al. (2014) to be compared with the SCPs here derived combined with those derived by Vieira et al. (2016) (see Fig. 5). The SCPs used for EPS and ReCiPe are available at the Electronic Supplementary Material.

Fig. 5
figure 5

Surplus cost potential (each it its original unit) from EPS (Steen 1999), ReCiPe (Goedkoop and De Schryver 2009) and those here derived combined with Vieira et al. (2016) for reserves, the latter in the x-axis. Both axes are displayed in logarithmic scale

In the EPS method (Steen 1999), future extraction costs in combined characterisation and weighting factors (ELU) per kilogram were calculated using future sustainable production technology costs. EPS covers a total of 69 resources, 65 mineral and 4 fossil (Roerbech et al. 2014; Steen 1999), for which 41 also have SCPs. The SCPs here derived are in general lower (4 to 7 out of 41) than the CFs derived in EPS. For a specific resource, the SCPs here derived can be up to 2 orders of magnitude larger (for cesium) and 5 orders of magnitude lower (e.g. cadmium).

The characterization factors for resources expressed as surplus costs, expressed in US dollars from base year 1998, in ReCiPe were compared to both sets of SCPs derived in this paper. ReCiPe covers 20 mineral and 5 fossil resources. There were 20 resources covered by both methods. When comparing the values from ReCiPe and the SCPs here derived, there is a maximum of 2 orders of magnitude difference. Although both ReCiPe and the surplus cost method here proposed both assess the extra costs to be paid for extracting future resources, there are key differences between both methods. Our SCP calculations apply mine-specific data for determining relationships between the extraction of resources and the increase in production costs, whereas ReCiPe adopted one constant value. Another important difference is the discounting of future costs applied in the ReCiPe method. In the method applied here, it was chosen not to apply discounting because, according to Hellweg et al. (2003), discounting across generations because of pure time preference should not be applied in LCA.

Ponsioen et al. (2014) derived three sets of characterization factors expressed as surplus costs for fossil fuels only using a method which combines a marginal modelling approach, future production scenarios and discounting. The surplus cost values derived in Ponsioen et al. (2014) with zero discounting (egalitarian perspective) for crude oil, hard coal and natural gas were, respectively, 15.4, 7.4 and 3.4 USD2013 per gigajoule (converted with Bureau of Labor Statistics (2014) from US dollars from base year 1998 to 2013). These SCP values are in the same order of magnitude as those obtained in this paper. The difference in the SCPs obtained can be explained by the difference in modelling approach (average vs marginal) and by the use of more recent data in this paper compared to Ponsioen et al. (2014).

4.2 Limitations

Our analysis is not without uncertainties which are further discussed below. First, we assumed that cost increases are expected from future extractions following the cumulative cost curves derived from current practice. Although there are several effects of resource extraction leading to cost increase, such as the need to use lower ore grade, more remote and more difficult to process deposits, new technology and innovations and new discoveries can lead to cost decreases (Humphreys 2013; Svedberg and Tilton 2006; Tilton and Lagos 2007; Yaksic and Tilton 2009). For fossil resources, future technology development is accounted for in the estimates made by the International Energy Agency (2013). For mineral resources, technological innovation, economies of scale and new discoveries may partly offset higher costs of extracting resources (Crowson 2012; Curry et al. 2014). Reserves beyond 20–30 years of consumption are rarely identified (Yaksic and Tilton 2009), meaning that new mining projects with lower costs than other running mines in the future may occur. For instance, over the past 30 years, copper resources have more than doubled (Tilton and Lagos 2007). However, over more than a century, copper extraction costs show neither great rise nor fall (Svedberg and Tilton 2006). Despite these observations, potential cost reductions due to new technology development were not included in the SCP calculations of minerals because this is unpredictable (Tilton and Lagos 2007; Yaksic and Tilton 2009). Neglecting these potential future cost reductions in the mineral SCP calculations implies an even lower contribution of mineral resources to the electricity resource footprint.

For the calculation of SCPs for fossil resources, future cost and production estimates from the International Energy Agency (2013) were used. Of course, these are forecasts determined on basis of the depletion of current basin sites and the development and use of new technologies for each specific fossil fuel, both of which are uncertain (International Energy Agency 2013).

The choice to integrate until MFE supposes that all the resource is consumed in the end, although there will be exponentially rising costs. This may not be realistic for resources that can be substituted by more abundant alternative resources. For mineral resources, this assumption can be tested by considering two types of reserve estimates: (1) reserves, which include only identified reserves that are presently economically viable and are likely to be consumed in the short-term future, and (2) ultimate recoverable resource (URR), which is by definition the ultimate quantity of economically viable reserves over all of human time. As an attempt to estimate URR, also called ultimately extractable reserves (UER) and extractable global resource (EGR), as explored by UNEP (2011a), a rough estimate of URR defined as 0.01% of the total amount in the crust to 3 km depth will be used. However, it should be noted that this estimate is highly uncertain because of unknown future circumstances that will continue to influence and modify economic viability and so the final URR estimate. SCPs derived with URR as reserve estimate are on average a factor 3 higher with a maximum of a factor of 5.7 larger than those derived with reserves (see Electronic Supplementary Material). Using the SCPs based on URR, we found a larger contribution of mineral resources compared to SCPs derived with reserves. For the renewable electricity production technologies, mineral resources contribute 8–29% (reserves scenario) and 17–46% (ultimate recoverable resource scenario) to the overall surplus costs. For nuclear technologies, the difference is up to 20% more contribution for URR from mineral resources. For fossil resources, only one reserve estimate was used, based on International Energy Agency (2013) estimates. Depending on climate policies and the future development of renewable energy resources, it is unclear whether all the fossil reserves will be actually consumed by society.

In Vieira et al. (2016), SCP values were derived for 12 metals and for the group platinum-group metals using mine-specific tonnage and production costs. To allow the applicability of the SCP method within LCA practice, SCP values were derived for over 60 metals, minerals and mineral resource groups extrapolating from price data from 2013. Prices were used because a good correlation was found between the 13 SCP values explicitly derived in Vieira et al. (2016) and their prices. However, we would like to emphasize that this does not imply that there is a causative relationship between current price and SCP. Also, despite the high correlation found, some of the metals deviate from the regression line. For instance, the SCP derived empirically for zinc for reserves is 9.0 USD2013/kg whereas the value determined with the price correlation is a factor of 22 smaller (0.4 USD2013/kg). This implies that metals which are co-mined predominantly in zinc deposits, as is the case for indium, may also be underestimated when their SCP values are derived on basis of the price correlation. For photovoltaic electricity production, the technology is often associated with indium (Öhrlund (2011). Indium has the largest contribution to the total surplus cost per megawatt hour for photovoltaic but still smaller than 1%. If the SCP for indium is directly derived from the zinc SCP via (PRICEindium / PRICEzinc) × SCPzinc, the SCP of indium reserves increases from 1.0 × 102 to 2.6 × 103 USD2013/kg. Although this is a significant increase in SCP, the contribution of indium to the total surplus costs of the PV example above becomes 2% with fossil resources remaining the dominant factor.

Also, there are various limitations with using prices. Prices fluctuate due to factors other than mining and milling costs, which are the ones considered in the SCP approach. For instance, if demand rises much faster than the production capacity of mines, prices increase even with constant production costs because of lack of equilibrium between supply and demand. The price data used for extrapolation is average data for global trade of that metal or mineral. This means that price of both primary as well as secondary material, when applicable, is considered. This is particularly relevant for metals where supply from recycling and the difference between the prices of primary and secondary sources are large (Reck and Graedel 2012). For instance, ‘base metals’ such as iron, copper and zinc have end-of-life recycling rates of above 50% (UNEP 2011b) implying that the price of these scrap metals in also lower than from primary sources. Finally, using the average prices of 2013 instead of long-run trend prices as a predictor for mineral SCPs is also a point for discussion. Prices reflect market and geopolitical conditions and policy, thus other than resource scarcity, at a certain moment in time. There can be significant differences between real prices and long-run trend prices (Cuddington and Nülle 2014; Rossen 2015; Tilton and Lagos 2007). Yaksic and Tilton (2009) argue that long-run trend prices offer the most useful insights regarding mineral depletion, and Tilton and Lagos (2007) state that the long-run future equilibrium price can be assumed to converge toward production cost all other things being equal. For this reason, the correlation between SCPs and long-run trend prices was analysed as well (see Electronic Supplementary Material). According to Rossen (2015), the long-run price trend lasts longer than 70 years so average prices for a period larger than 70 years were used, whenever available. The log-linear relationship between surplus cost potential and the long-run trend price of a metal or mineral resource has an explained variance of 92% and an intercept of −0.50. This means that the SCP is typically a factor of 3 (100.50) smaller than the price of a metal or mineral resource, whereas when using the real price of 2013, it was a factor of 5. Also, comparing prices from 2013 and long-term (historical) prices, we found a linear relationship with a high explained variance (R 2 = 0.93) (see Fig. S2 in the Electronic Supplementary Material).

The datasets used for the electricity production technologies are from ecoinvent v3.2. In life cycle inventory libraries, data for capital goods is only roughly modelled and thus incomplete and of relatively poor quality so the surplus costs here calculated may be underestimated (Arvesen and Hertwich 2012). Also, in the calculation of the surplus costs, SCP values for lignite and peat were set at zero because there was no data to derive specific SCP values for these fossil fuels. If we adopt the SCP value of hard coal for lignite and peat as a first proxy, the surplus costs of all electricity technologies considered change less than 8% (see Electronic Supplementary Material for details).

5 Conclusions

The Surplus Cost Potentials calculated for fossil and mineral resources in our study can be considered as a useful step toward a coherent comparison of resource scarcity caused by renewable and non-renewable electricity production technologies. Of course, SCP is an extra life cycle impact assessment indicator next to others that cover other environmental effects, when comparing different electricity production technologies. Fossil and mineral resources have been assessed in the same way, both based on cost-cumulative production relationships as basis to derive SCP values. There are, however, two main differences between them: (1) the SCPs for fossil resources have been derived using past as well as future costs whereas the mineral SCPs have been derived by using current cost data of individual mines only, and (2) for minerals, two reserve estimates were used whereas for fossils only one estimate for future production was used. Fossil fuels always dominate surplus costs of electricity production compared to mineral resources, even for renewable technologies that do not require burning of fossil fuels to produce electricity. The electricity production technologies fuelled by fossil resources result in the largest surplus costs and hydropower from reservoir and run-of-river have the lowest surplus costs per megawatt hour produced. This case study shows that the surplus cost method facilitates the evaluation of trade-offs between mineral and fossil resource use in life cycle assessment.