Response to Commentary: Further clarifications on a more basic modeling framework for life cycle methods

COPYRIGHT © 2023 Schaubroeck. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms. Response to Commentary: Further clarifications on a more basic modeling framework for life cycle methods

In response to the commentary of Heijungs and Suh (2022), some further clarifications are provided concerning the more basic modeling framework brought forward in a previous article published in this journal (Schaubroeck, 2022a). Since Heijungs and Suh (2022) do not highlight clearcut errors, as will be made apparent, the main messages and content of the particular work (Schaubroeck, 2022a) still stand and hold, but its content could benefit from further explanation. Detailed replies to the comments of Heijungs and Suh (2022) are presented in the Supplementary Table 1. Replies to comments related to the book by Heijungs and Suh (2002) are restricted to the latter. Clarifications concerning the modeling framework and its relevance are presented in the following three points. Note that a preprint version of this article has been uploaded online (Schaubroeck, 2022b).
First, it is good that Heijungs and Suh (2022) again bring forward that the matrix inversion approach (Linear inverse modeling) is not the only computational structure/framework for LCA. Hence, other or more basic modeling frameworks, such as the one presented by Schaubroeck (2022a), can be brought forward and can be of relevance.
Second, when it comes to the notation style of the formula of Equation (3) in the work of Schaubroeck (2022a), Heijungs and Suh (2022) question having a function with a set of output variables and a set of input variables, but some references point out this possibility, for example. In fact, one needs to rather consider the particular function as p p : V n → W n (considering V n and W n are co-domains with n and m type of elements), for which a unique set of input variable amounts will provide a unique set of outcome variable amounts. Yet, indeed, the formulation of the equation in question is somewhat of an uncommon way of formulating things, especially in the field of LCA. To further clarify, the following two alternative versions of the equation are presented (using parameters as specified in the original work): In this equation, each parameter F i has a corresponding function p i that only has this parameter as a single output variable, which is more conventional. One can also just write it as a collection: When it comes to the type brackets and their order, namely, f (. . . {. . . }), this is also somewhat unconventional. The second set of brackets "{. . . }" is strictly speaking redundant, but it was used to distinguish the set of flows separately from time and other conditions. Moreover, an example is presented for further clarification, as also advised by Heijungs and Suh (2022). Consider the case related to Figure 1 on heat and electricity in the original work (Schaubroeck, 2022a). That figure depicts industrial processes, but there is, in fact, also the implicit process of energy supply selection, which represents dynamics and is more interesting to study. Using Equation (3) from the original work (Schaubroeck, 2022a), one can derive the following formula: electricityfromgrid, heatfromboiler, electricitycfromCHP, heatfrom CHP} = energysupplyselectionprocess p electricity, heat ,

timepoint
( If we consider a simplistic version where one can readily switch between the energy sources and there is no restriction due to fixed ratios of the Combined Heat and Power generator (CHP), the following oversimplified mathematical expression can be obtained for each timepoint: To put it simply, only if both heat and electricity are simultaneously asked, the CHP is used, otherwise, electricity is obtained from the grid and heat from a boiler.
Third, the main advantage of the presented framework in the original work (Schaubroeck, 2022a) is its openness to further developments, including non-linear developments. The value of the associated abstractness seems questioned by Heijungs and Suh (2022), but it is the abstractness that makes it general and open. Also note that, in the original work (Schaubroeck, 2022a), it is not claimed to have introduced new specific non-linear frameworks, and alternatively, reference is made to many. Moreover, many other advantages of the original work (Schaubroeck, 2022a) have also been overlooked. To further emphasize and clarify this, some benefits/novelties of the original work (Schaubroeck, 2022a) are listed here: • An expression in mathematical terms of an open and general framework in the context of LCA, allowing to guide further users and developers to be more free in modeling development. • Specification of how this computational framework is aligned with a more basic and open framework for LCA, as recently specified by Schaubroeck et al. (2022), and can be specifically applied for attributional and consequential LCA. • An elaboration on dynamic effects beyond a specification over time, and how this is limited in the matrix inversion approach. In fact, this aspect is the most elaborated and novel. • A specification of a case on co-product dynamics is illustrated in Figure 1 of the original work (Schaubroeck, 2022a). • Pointing out further possibilities to improve the frameworks of DyPLCA (Pigné et al., 2020) and Temporalis (Cardellini et al., 2018): "For non-linear modeling, process equations could be readily changed, possibly grouped by types of processes to facilitate execution. In the case of dynamics, specifically for graph-based search tools, dynamic graphs could be considered that can alter over time (Holme, 2015;Vernet et al., 2022)". In the original work, it has been forgotten to mention that these frameworks also include a complete temporal differentiation, with a complete temporal database for a version of ecoinvent 3.2 in the case of DyPLCA. This benefit should be considered when evaluating the extent of limited practical flow coverage with the given computational power.
Overall, with these further clarifications, I hope to have addressed the points raised by Heijungs and Suh (2022) concerning the framework.

Author contributions
The author confirms being the sole contributor of this work and has approved it for publication.