Theory without practice: a reply to the note from Heijungs on the average versus marginal debate in Life Cycle Impact Assessment

1 Environmental Science, Radboud Institute for Biological and Environmental Sciences (RIBES), Radboud University Nijmegen, Nijmegen, The Netherlands 2 Ecological Systems Design, Institute of Environmental Engineering, ETH Zurich, Zurich, Switzerland 3 Industrial Ecology, Department of Energy and Process Engineering, Faculty of Engineering, Norwegian University of Science and Technology (NTNU), Trondheim, Norway


Introduction
recently discussed in his note "The average versus marginal debate in Life Cycle Impact Assessment: paradigm regained" how marginal and average effect factors may be derived in a multistressor context. He derived the average multistressor effect factor as follows: and where msPDF is the multistressor potentially disappeared fraction due to exposure to concentration of pollutants N and P, g(C N , C P ) is the mathematical expression of the combined (3) EF average,p = f 1 C P C P effect of pollutant N and P, and f1(C N ) and f2(C P ) are the mathematical expressions of the effect caused by the background concentration of pollutant N and P, respectively. It is, however, easy to see that this way of calculating the average effect factor is not going to work which Heijungs also admits by stating that "Only by making ad hoc allocations (…) would we be able to separate the contributions, as a first step in an averaging procedure." A second argument of Heijungs against the average effect factor approach is the presumed requirement that "an average, per-item value, multiplied by the number of items, should give the total value." In more practical terms, this means that, according to Heijungs: Heijungs showed with a theoretical example that Eq. 4 does not hold. Based on these two arguments-"ad hoc allocation" and "violation of the impact balance"-Heijungs concluded that "The average approach, as propagated by Huijbregts et al. (2011) … runs into unsurmountable problems in the case of non-additive combination rules…".
Here, we will proof that, in contrast to these statements from Heijungs, (i) ad hoc allocations are not required in the calculation of average effect factors and that (ii) the approximation of Eq. 4 is likely to hold in practice. Regardless the conceptual opportunities for both the marginal and average approaches, we conclude our reply by shortly discussing the practical obstacles to tackle before the marginal or average effect factor approaches may replace the more common linear approach. and There is no need to subdivide the impact between the pollutants upfront, as suggested in Eqs. 2 and 3 by Heijungs, one just need to know the change in the multistressor impact if the stressor of interest is excluded. In full equation mode, this means that the multistressor-based average effect factor can be calculated as: and where msPDF n-N and PDF n-P are the multistressor potentially disappeared fractions caused by the total set of pressures n without the contribution of the stressor under consideration (N or P), and PDF N and PDF P are the potentially disappeared fractions caused by respectively pollutant N and P.
In words, it means that one can calculate the average effect factor as the change in the multistressor Potentially Diseappeared Fraction (msPDF) of species by putting the stress of the pollutant of interest to zero. This version of the effect factor represents the average change in the msPDF over the concentration range from ambient to zero concentration of the pollutant of interest. Note that solution I is an average approach with the underlying assumption that the pollution level of all other  substances than the one assessed remain constant. This may be questionable, as also for those substances, we then assume a target when they are assessed, so indeed this is not a perfect approach. This means that average effect factors according to this solution cannot always be used to approximate the msPDF in a location or region. This is, however, also not the intention of Life Cycle Assessment studies where emissions related to a functional unit are evaluated.

Multistressor average effect factor derivation: solution II
As an alternative solution, the multistressor average effect factor can be calculated in a way it obeys the rule that "an average, per-item value, multiplied by the number of items, should give the total value": and where msPDF is the multistressor potentially disappeared fraction of species, ∑ i msPDF n−i equals in this example msPDF n-N + msPDF n-P . This way of allocating impacts to individual stressors is conceptually equal to the way average characterisation factors are derived for different types of land occupation with non-linear species area relationships (see Chaudhary et al. 2015). It is also consistent with the way the relative importance of individual stressors to the overall impact on species abundance is derived in biodiversity modeling (Schipper et al., 2020).

Practical example: eutrophication in European lakes and rivers
To show how the two solutions work in practice, we derived average effect factors for nitrate (NO3-) and total phoshorous (TP) for lakes and streams in Europe. We used response relationships for invertebrates, determined for the two stressors (i.e., NO3-and TP) in two freshwater types (i.e., lake and stream) combined with average monitoring concentrations in European river basins for year 2011, as reported in Azevedo et al. (2015). We found that the differences between the two average effect factor solutions are always smaller than a factor of 1.5 (Fig. 1). Given the many sources of uncertainty and variability in the derivation of characterisation factors for eutrophication (see e.g. Azevedo et al. 2013), we do not consider this factor of 1.5 relevant from a practical point of view.

Conclusion
In summary, there is no mathematical proof that a multistressor-based average effect factor is ill-defined as we have shown here. On the contrary, there are even two possible solutions for the issues identified by Heijungs. We also showed with a practical example of eutrophication impacts in European river basins that both average approaches give approximately the same results, implying that a violation of the impact balance by the approximation of solution I is not likely to occur. This means that there are, conceptually speaking, no winners or losers in the average vs marginal debate in Life Cycle Impact Assessment. Whether average and marginal (multi-stressor) effect factors can be calculated in practice with the data available, is, however, another question that remains to be answered. Although we were able to derive average effect factors of invertebrate species for the nitrogen and phoshorus pollution in European lakes and rivers, the situation is different for the vast majority of other stressors. In practice, most effect factors, particularly for chemicals, are derived by linear extrapolation without marginal nor average considerations (see e.g. Huijbregts et al. 2017). The linear assumption is simply an admission of our own ignorance on the functional form and shape of the response relations and the environmental background situation for the vast majority of pressures. The priority of future efforts should thus lie in obtaining reliable fate, exposure, and effects information.