Toward more comprehensive environmental impact assessments: interlinked global models of LCIA and IAM applicable to this century

Despite the long-standing demand for research on dynamic lifecycle assessment (LCA) for policymaking, only a few studies have addressed this subject in conjunction with other systems thinking disciplines, such as energy systems models (ESMs) and integrated assessment models (IAMs), which have achieved tremendous success in assessing climate policies in future scenarios. This study explains our methodological advances in the global application of LCA incorporated in IAMs, specifically dose-response functions, models, and future scenarios. We investigate the application of the lifecycle impact assessment method based on endpoint modeling (LIME), developed under the current environmental situation in Japan and globally, to be consistent and compatible with our IAM, which comprises three resource balance models and a simplified climate model. The IAM endogenously generates most inventories consistent with energy scenarios and climate policies linked with the applied LIME. The IAM and LIME are formulated to minimize the discounted sum of supplying the cost of resources over their lifecycles (i.e., from development to end-of-life) to generate time evolutions for the endpoint impacts over this century on a global scale with/without the 2-degree Celsius (2DC) target in a 100% renewable energy scenario. Unlike existing LCA+ESM/IAM studies, which focus on power generation technologies and related (in)direct embedded energy consumption on a lifecycle basis, our model’s expansion to mineral and biomass resources, in addition to energy, has the following novel results: (1) The following inventories in the 2DC target are generally lower than those in business as usual (BAU): temperature and sea level rise, natural resource, and waste discharge; further, SOx emissions are significantly reduced by reducing coal production while increasing forestry. (2) The environmental impacts on the four endpoints of minerals, land use, and land-use change, with the exclusion of energy-related impacts, are significantly larger than those related to energy. (3) Finally, by ensuring inventory reduction, the 2DC target scenario can reduce overall endpoint impacts (by maximum around 20%), except the impacts on biodiversity resulting from forestry expansion to meet predetermined targets. Unlike mainstream IAM analyses, we incorporate LIME, instead of energy- and biomass-related resource and climate change impacts alone; our model thus provides a comprehensive perspective on various natural resources and their impacts on a lifecycle basis. The exclusion of the weighting process and retention of the four endpoints enable us to easily interpret the results. Further, this application of LCA to IAM enables us to further understand and assess natural resources and environmental impacts.


Background: explorative lifecycle assessment
Today, humankind faces new challenges in environmental externalities, as well as growth and scarcity (Simpson et al. 2005). Decade-long efforts have finally succeeded in globally mainstreaming climate change policy, with the initiation of the Intergovernmental Panel on Climate Change (IPCC), the International Resource Panel (IRP) for natural resources, and The Economics of Ecosystems and Biodiversity (TEEB) for biodiversity. However, newer issues are emerging, such as the disposal of plastic bags into the world's oceans. More comprehensive and holistic approaches are in demand since the environmental externalities have some common underlying economic activities, such as resource development, resource use, and consequential impacts.
Specific subjects have evolved new tools to tackle these issues. Systems thinking tools, energy systems models (ESMs) and/or integrated assessment models (IAMs), and lifecycle assessment (LCA) are commonly used to assess technology and the economy. They are complementary as well, in that the ESM/IAM performs dynamic (time evolution) in scenario analysis; however, it disregards indirect greenhouse gas (GHG) emissions on the bases of lifecycle. Therefore, LCA is better suited to assess the environment and geography. In addition, GHG accounting incorporates environmental impact assessment, whereas time evolutions or scenarios receive only basic static analysis and limited research effort.
In "Digesting the alphabet soup of LCA," Guinée et al. (2018) label the LCA studies that explore scenarios of potential futures as "eXplorative LCA," or XLCA. Our study falls under this recently named category. Our work contains multiple "X" scenarios. "Backcasting LCA" explores how normatively defined sustainability levels (e.g., planetary boundaries) can be met. We also explore 2-degree Celsius (2DC) targets under the 100% renewable energy scenario (Ren100), called "scenario-based LCA." We do this by combining our lifecycle impact assessment (LCIA) model-the lifecycle impact assessment method based on endpoint modeling (LIME)-with our three resource balance models (one of which is a bottomup type energy systems model (ESM)) and a simplified climate model, categorized as "integrated" LCA. Since our IAM is a forward-looking, nonpredictive, and perfect foresight model, our work involves an aspect of "anticipatory LCA." Finally, we illustrate the damages at the four endpoints defined in our LIME over the time horizon to the year 2100 through our IAM, which corresponds to "prospective LCA." With this definition of eXplorative LCA, this study undertakes a modeling exercise to interlink the LCIA and ESM/IAM. Here, the words "interlinked" and "integrated" are used interchangeably and distinguished imprecisely. In general, interlinked seems to satisfy situations before the completion of integration or to separate one model from the other, whereas integrated is used after the integration of the models.
Thanks to Guinée et al. (2018), we refrain from dwelling on the varieties of LCA modes and, thereby, avoid muddled definition-related discussions. However, our further discussions focus on the integrated LCA mode, integrated LCIA models and systems (e.g., energy scenarios), and the scenario-based LCA mode.
1.2 Literature review: integrating lifecycle assessment models and the Energy System Models and Integrated Assessment Models Vanderpaer and Gibon (2018) summarize prospective LCA from five attempts presented at a conference on the integration of energy scenarios into the LCA. They examine the challenges encountered in efforts to combine the two model domains (i.e., LCA and ESM/IAM) in terms of organizational, structural, technological, and methodological inconsistencies. Further, the authors present the perspectives of LCA practitioners on enabling data sharing, open sourcing, and best practices; however, due to their complexity and resource-intensive nature, they omit some perspectives. The current study includes these omitted perspectives as well.
There are three patterns of integration of LCA models and ESM/IAM: First, a process-based "cradle to grave" approach using LCA data (e.g., Tokimatsu and Kaya 2004;Hertwich et al. 2015;García-Gusano et al. 2017a, b;Volkart et al. 2018), which is the soft-linkage of these models in some instances. Second, multiplying the damage or economic factors (also called coefficient in other studies) in the LCIA database with inventories that are endogenously generated from the models through fully integrated, sometimes hard-linked, models (e.g., Pehl et al. 2018;Arvesen et al. 2018;Boubault et al. 2018;Tokimatsu et al. 2019). The third is a hybrid of the aforementioned and other categories.
The first approach, requiring process-based cradle-to-grave work, applies processes based on detailed inventory data for "input to and output from" the processes. These are either retrieved from LCA databases (Hertwich et al. 2015;García-Gusano et al. 2017a, b;Volkart et al. 2018) or are background data (Tokimatsu and Kaya 2004), for all the processes and stages following the blending of the database and model, along with the lifecycle of electricity generation (resource extraction, transport, distribution, process, construction, operation, decommission, and disposal). García-Gusano et al. (2017a) examine the Spanish power sector through the Integrated MARKAL-EFOM System (TIMES), Spain's inventories and mid-and end-points, using the ReCiPe 2018 LCA database. They further explore modeling exercises integrating the LCA and the ESMs to show indicators over time, ranging from technology-based, environmental externalities, and energy security (García-Gusano et al., 2016. Unlike these studies that run ESMs themselves, Hertwich et al. (2015) borrow power-mix scenarios from other studies to soft-link with several types of databases and models, including the LCA, to quantify energy and material requirements and environmental impacts simultaneously. Uniquely, they apply a hybrid LCA model (THEMIS) combined inputoutput framework with physical process lifecycle inventories.
Our prior studies focused only on nuclear power and the fuel cycle inventory at the global level (Tokimatsu and Kaya 2004) using an ESM, simply applying the Japanese nuclear fuel cycle database in Tokimatsu et al. (2006) (which generates Japanese nuclear power scenario manually, without applying ESMs), with minor modifications.
The primary merit of this type of analysis is that it can account for the time evolution (or dynamics) of inventories and the consequential impacts on the mid-and/or end-points of technology types, especially those that require massive energy in processing and operation, such as nuclear fuel cycles, and unlike those in which the most energy consumption is in the manufacturing stage (e.g., photovoltaic, wind).
The second approach, simply multiplying the factors, fully integrates LCA models with others by endogenously generating inventories from the ESMs and/or IAMs and then multiplying the factors from LCA models. This derives the metal requirements of power generation technologies on the lifecycle basis (Boubault et al. 2018), indirect energy consumption separated from direct energy in an LCA model  to apply global power generation mix , and the environmental external costs (Tokimatsu et al. 2019). Arvesen et al. (2018) conduct a methodological investigation in the precise separation of LCA factors for both direct and indirect allocation of energy consumption, compatible with the IAM framework. They applied the separated factors to two research directions: energy ) and endpoint impacts. These studies are similar to Boubault et al. (2018) in examining global power generation technologies; however, they chose the prospective and consequential LCA approach while Boubault et al. (2018) select the "allocation at the point of substitution" method and did not differentiate between the direct and indirect use of material flows and stocks.
The abovementioned LCA modes uncovered the third approach. Tokimatsu et al. (2016b) choose a soft-linkage approach, in which the model generates inventories, emission factors, and power production by multiplying the economic factors of LIME to obtain the time evolution of the costs and benefits of hypothetical power plant construction projects. Patrizio et al. (2017) apply background data from LCA (GEMIS, IMPACT), fully integrated as environmental external costs into a static ESM (BeWhere) without time evolution.
Some studies applying the cradle-to-grave lifecycle basis do not use system models; instead, they generate the scenarios manually Tokimatsu et al. 2006). Stamford and Azapagic (2014) assess power generation scenarios for the UK in 2070 in terms of the three pillars of sustainability: techno-economic, environmental, and social impacts. They use GEMIS, Ecoinvent, and Gabi as data sources. Other studies (Santoyo-Castelazo, Stamford, and Azapagic) are similar. Tokimatsu et al. (2006) develop their original Japanese nuclear power scenarios from 1960 to 2100 to show the time evolution of CO2 emissions and its intensity in the sector.

Objectives
Our study aims to show the future paths of endpoint damages caused by various impact categories at the global level in the twenty-first century and to offer a more comprehensive assessment of environmental impacts and resource use. The research question is as follows: what is the physical reduction in environmental impact (e.g., biodiversity) from a direct CO2focused climate policy (i.e., 2DC target)? The impacts covered are global warming (GW), mineral resources (MR), land use and land-use change (LU), waste (WA), local air pollution (AP), acidification (AC), and ozone depletion (OD). These impacts are due to resource management in all lifecycle stages (production, transformation, use, and end-of-use) of fossil fuels, uranium, metal (ferrous and non-ferrous base metals), non-metal (limestone), and bio-related products (agricultural products-rice, corn, and wheat-and wood products-logs, wood pulp, and timber).
Our study is different from previous studies in that it adopts a fully integrated LCA model in its IAM. The term "fully integrated" does not imply the inclusion of the weighting process in LIME; it excludes the weighting process in the LCIA framework and, instead, emphasizes Dose-Response Function (DRFs). One reason for this exclusion is that the damages are not monetized and presented in their own units (e.g., Disability Adjusted Life Years (DALY) on human health (HH)), and the monetary process is unfamiliar for some readers. Another reason for this exclusion is that it disregards the internalization of the environmental external cost already published by us (e.g., Tokimatsu et al. 2019) and focuses on the DRFs to reveal future trajectories of endpoint damages that did not appear in our modeling exercises.
Since we exclude the weighting process and environmental external cost, this study uses a conventional (internal) cost minimization model that endogenously generates these paths. Our IAM is distinct from others in that it uses techno-economic ESMs due to its comprehensiveness. This study utilizes the latest version of impact assessment based on LIME3 and the inclusion of the mineral resource balance and the simplified land-use models. Moreover, we apply the newest scenarios in climate policy (2DC target) originally investigated under Ren100, based on our prior modeling exercises .
This paper is structured as follows: Sect. 2 explains our methodology, focusing on DRF scenarios that are modified and blended with our IAM. The subsequent sections discuss the study's results, discussions, and conclusions. Appendix 4 provides a list of abbreviations, whereby readers can better understand the manuscript. This study interchangeably uses words with similar meanings but different nuances without providing precise definitions (e.g., interlinked, in the title of Sect. 4.4).

Methodology
We focus on the global application of the LIME framework combined with our systems models: the overall model in Sect. 2.1, the resource balance models (i.e., systems models) in Sect. 2.2, adjustment of the DRFs in LIME in Sect. 2.3, connection of adjusted LIME and the other models in Sect. 2.4, and scenarios in Sect. 2.5. Finally, we discuss the strengths and weaknesses of our strategy in terms of model choice in Sect. 4.3.
For brevity and the ease of understanding, Appendices 1 to 4 contain the following details: the original modelthe case before the global application of the LIME framework; details of DRFs; three resource balance models (Tables 2, 3, and 4; Fig. 1); and a list of abbreviations, parameters, subscripts, and variables. LIME is an LCIA method that includes the three main steps of LCIA (i.e., characterization, damage assessment, and weighting; see Appendix Figure 4). It determines the weights of fewer endpoints (damage-oriented approach), rather than impact categories requiring 10 items or more, using the economic evaluation methods.

Overall model illustration
In the prospective LCA, various ESMs/IAMs attempt to use such techniques for optimization, input-output, simulation, and general equilibrium, with different time spans (e.g., 2030, 2050, and 2100) and geographical coverage (i.e., country, region, and global). We developed a dynamic optimization model (i.e., minimize the sum of the discounted supply costs) with a computation time horizon from 2010 to 2150 with 10-year intervals and for 10 world regions. Figure 1 illustrates our overall modeling frameworks. The DRFs, as illustrated by arrows, connect the inventories in blue to the four endpoints in yellow (see also Table 1). With the weighting process excluded in the yellow box (LIME2 and LIME3), this study focuses on the direct, consequential damages at the endpoints of the inventory inputs in their physical units rather than in monetary units. The inventories are endogenously generated (only Ozone Depletion Substances (ODS) is given exogenously) from the four models in the blue box; namely, the simplified climate and the three resource balance models (energy systems, mineral resources, and bio-resources with land-use).
The objective function of this model is to minimize the sum of the discounted supply costs (including investment and fixed and variable operation and maintenance; see Appendix 3) of the resources over time from the three models (excluding the simplified climate). This setting implies that the world economy pursues economic efficiency. As environmental outcomes become significant social concerns, these are constrained or controlled; for example, greenhouse gas emissions or temperature increases caused by climate change. However, economic activities (i.e., resources production) are still adjusted according to constraint or control via their (private) cost minimization; thus, they do not internalize external environmental costs monetized from the outcomes.
We leave the investigation into weighting and internalization to other studies Tokimatsu et al. 2016a;Murakami et al. 2018;Tokimatsu et al. 2019). We excluded the internalizing external environmental costs-our original approach using LIME-to follow other studies to present the similarities and disparities of our results.

Three resource balance models
The resource balance models are formulated through linear programming by dynamic optimization according to perfect foresight (i.e., demand and unit cost of variables are given deterministically) to minimize the discounted supply costs covering the overall lifecycle of resources and fully integrated bottom-up technologies. Other studies contain the details of the models (Tokimatsu et al. 2016c(Tokimatsu et al. , 2017a(Tokimatsu et al. , 2017b. Appendix 3 describes the model linkages, objective functions, and variable lists. The three models are fully integrated, and they consider representative resources from the three sectors. The energy systems model fully covers primary and secondary energy resources via some 100 supply-side technologies excluding end-use, unlike the TIMES/TIAM model families. The mineral resource balance model uses ferrous and non-ferrous base metals (iron and steel, aluminum, copper, lead, and zinc) and limestone for cement production, expressing the entire all-material lifecycle. The bio-resource and simplified land use cover the representative agricultural (rice, corn, and wheat; pork, chicken, lamb, and beef) and wood-based products (paper, wood lumber, boards, and fuel) based on five land patterns and its changes. These models express "balance" by supplying the resources to satisfy the demand generated from exogenously provided gross domestic product (GDP) and population scenario (i.e., SSP-2). Therefore, these models are not of the "general equilibrium" type nor do they encompass all sectors; however, they present sufficient resolutions for their consequential impacts.  1 Overall structure of our IAM used in the current study. Doseresponse functions (DRFs), which are illustrated by arrows, connect the inventories in blue to the four endpoints in yellow. The inventories are endogenously generated (except Ozone Depletion Substances (ODS), which is given exogenously) from the four models in the blue box, that is, the simplified climate model and the three resource balance models (energy systems, mineral resources, and bio-resources with land use). The white boxes are the compartments used by existing studies; most of them use energy systems models (ESMs) and midpoint impacts in a specific country or region till the middle of this century, whereas all the others shown in blue and yellow (including mineral resources and bio-related resources with endpoint impacts) are extensions of our global scale model that is valid up to the end of this century.
The models on energy systems and mineral balance determined the flow and stocks of energy and material using the least-cost mechanism over their lifecycle. In the simplified land-use model, the allocations of resource flow are provided by parameter settings, whereas the amount of produced resources is determined by the area's costs for land use, landuse changes, and production yields. The three resource balance models endogenously generate inventories related to resources (production, disposal, land use, and land-use changes) and their consequential outcomes (emissions of GHGs, NOx, and SOx), where increases in temperature and the sea level as a result of global warming are computed ex-post by the simplified climate model. The generated inventories (shown in blue in Fig. 1) are given to the DRF, corresponding to the arrows shown in the yellow part of the figure, to calculate the four endpoint damages.

Outline of DRFs
The DRFs are generally expressed in the following aggregated forms (see details in Appendix 2) to modify the factors (DC j ) provided by LIME2&3, using the terms [A j,t , B t ] to apply our IAM: The first term (DF j,t ) represents the damage functions expressed by per capita GDP (y j,t ), global mean temperature rise (TR t ), and sea level rise (SR t ). This corresponds to the damage functions (DF) for global warming (GW) at all endpoint impacts (i.e., human health (HH), natural resources (NR), and biodiversity (BD)).
The second term [DC j ⋅A j,t ⋅B t ] denotes the damage of all impact categories other than GW, with damage factors (DC j ) and their adjustment terms (A j,t , B t ). The damage factors (i.e., constant level) are sourced from their lists provided in LIME2 and LIME3.
Since these lists include the environmental damage of either almost all the countries (LIME3) or Japan alone (LIME2), the first term A j,t reflects the order of magnitude of the endpoint damage in accordance with the Shared Socioeconomic Pathway (SSP) (van Vuuren et al. 2012). It follows that (i) population density (PD j,t ) for HH, (ii) per capita GDP (y j,t ) for NR, and (iii) primary productivity (PP) and potential NPP (PN j,t ) are determined by the Chikugo model (Uchijima and Seino 1985) using temperature and precipitation data from the Model for the Assessment of Greenhouse Gas Induced Climate Change (MAGICC) (Wigley 2010). These adjustments, which correspond to the first bracketed term in Eq. (3), are ratios of those levels in the target year and region divided by the values in the year 2000 (for Japan and in the LIME2 case). In addition to this ratio adjustment, two further adjustments are required (OD j,t in the second term) for MR (using LIME3) and (S(α j ) in the third term) BD in WA (LIME2), the former for ore degradation (OD j,t ) over time in copper, lead, and zinc and the latter (S(α j )) by a ratio of population density of endangered species in the target country divided by the values for Japan, corresponding to the second and third terms in Eq. (3), respectively.
The B t term expresses the increase in biodiversity impact in accordance with the reduction in area β(TR t ) of natural land (i.e., forestland and grassland) through the land development by TR j,t . This adjustment is required to determine the biodiversity impact in any impact category (i.e., mineral resources (MR), land use and land-use change (LU), and waste (WA)). The damage in NR can be broken down into the following three areas: agricultural products, land resources of submerged land area created via the SR t [m], and energy resources according to ambience (i.e., cooling, heating). The damage function is formulated from damage factors [ton/°C, km 2 /m, degree day/°C] per unit increase in TR t and SR t , multiplied by their increase. For energy resources, the relative increase in damage is also adjusted by degree days, heat demand GDP intensity, and per capita energy consumption expressed by functions of y j,t . Unlike other endpoint impacts, their physical units are converted into monetary terms by multiplying their values with the market prices in LIME3. To reflect this notion, combining our resource balance models, we use shadow prices obtained endogenously, except for land area.

Damage functions for global warming impacts
The damage factors in BD are sourced from Tang et al. (under review), which was revised based on the lists in LIME3 for the biodiversity impact of global warming (DC j,i = GW,e = BD ) [EINES/yr] on a grid (approximately 10 km 2 ). From this, we aggregated our regional damage factors.

Damage factors for impacts other than global warming
The impact categories included here are MR, LU, WA, AP, AC, and OD. The PP is a common endpoint for all impact categories, while BD is for MR, LU, and WA; HH and NR are for AP, AC, and OD (see Table 1). The damage factors are sourced from LIME3 for MR and LU, while the others are from LIME2. The ore degradations in copper, lead, and zinc are tied to cumulative mine production, endogenously computed in a mineral balance model and integrated into our overall model over time. This degradation leads to the discharge of huge overburdens of deep mining and land development, resulting in the damages on PP and BD.
The damage factors (DC j in Eq. 2) are divided into two terms. One is atmospheric concentration rise (e.g., NOx, SOx) by emission increase (denoted as dC/dE) and the other is DALY loss by atmospheric concentration rise (dDALY/dC). Hence, the DALY loss by increased emissions is expressed as follows: The first term, dC/dE (concentration increase by unit emission), is developed by employing both the plume and puff models, which assume a linear relationship between concentration and emission (note that "substance emission speed" is used instead of "emission" in these models), while a nonlinear assumption was made between concentration and local weather conditions (e.g., wind direction, wind velocity, and atmospheric stability). In the development of LIME, dC/dE is weighted average of results by the model calculations under the three weather conditions; however, the adjustment overlooks the condition changes in location and point in time, implying that adjustment of this term (dC/dE) simply follows the employed models.
The second term (dDALY/dC) is divided into the following three terms: Each term has the following units: (case person −1 yr −1 micro-g −1 m 3 ), (person m −2 ), and (DALY case −1 ); hence, dDALY/dC has the unit (DALY m −2 yr −1 micro-g −1 m 3 ). The first term, "contribution of concentration to increase of rate of death or disease," is a constant, implying that death or disease rate is linear and proportional to concentration and that it is common across background socioeconomic scenarios. The third term, "DALY per death or disease" is also independent from the scenarios, while only the second term, population density, is directly linked with the scenarios. This is why is the population density is adjusted (A j,t in the Eq. 3).

Adjustment for biodiversity
B(β(TR t )) is required for all BD-related impact categories, while S(α j ) is used only for WA based on LIME2. S(α j ) is an adjustment term between the Japanese values and the target region, regarding the increase in extinction risk of certain species, where α is the ratio of population density of the endangered species in the target country, divided by the Japanese values. We found correlations between the numbers of endangered species and their increased risk per 1 ha. To simplify the regression analysis, we used domestic mesh data (10 km 2 ) in numbers between those listed in the Red Data Book and the increase in the extinction risk of Japanese vascular plants identified by the environmental impact assessments. The function S(α j ) is calibrated using Japanese soil and gravel industry data (Itsubo and Inaba 2014).
B denotes the extinction risk increase on a cumulative basis from the original state to a particular point in time. Since the extinction risk is expressed by a log-linear relation with the number of vascular plant species in the Japanese Red Data Book, B can be derived by an integral up to the area change ratio β (between two states: the current and original) being proportional to the number of species. We adopt three relational patterns between extinction species and β from Thomas et al. (2004); the best data were chosen to tie with the corresponding temperature rise (TR t ).

Connecting Dose-Response Functions with inventories generated endogenously from other models
The overall endpoint damages (ED e,i,j,t ) are generally expressed as the following aggregated form: (1) for GW, DR e,j,t,i = GW,v is a direct route of TR t and SR t (interpreted as parts of inventories) from the simplified climate model; (2) with the exception of GW, DR e,j,t,i = GW,v is multiplied by the inventories (Inv j,t,i m GW,v ) unrelated to GW after summing all the inventories (v). We then generate consequential endpoint damages using DRF.
Appendix 2 contains detailed specific formulations. For GW, formulations are given to each endpoint (i.e., HH, NR, and BD). NR is further divided into three resources (i.e., agriculture, land, and energy) to which formulations are given. Except for GW, damage factors from LIME2 and LIME3, which are created to multiply inventories, are simply multiplied.

Scenarios investigated by this study
Our 100 renewable energy scenario (Ren100) was developed from an energy scenario dominated by gas and renewables after allowing additional settings and assumptions. These are as follows: (i) radical expansion in the allowable maximum share of renewables in the power sector, (ii) phase-out of nuclear and fossil fuel in the production of hydrogen and liquid fuel, in addition to the power and heating sectors by 2100 ). The Ren100 scenario was considered inherently incompatible with nuclear power; hence, the nuclear phase-out is in addition to the phase-out of fossil fuels in the energy sectors. To our knowledge, detailed 100% renewable energy scenarios have not been adequately addressed at the global level (Føyn et al. 2011;Hong et al. 2015), although there are many such studies at local or national levels (Connolly et al. 2011;Krajačic et al. 2011;Liu et al. 2011;Mathiesen et al. 2011;Østergaard et al. 2011a, b).
Based on the Ren100 scenario, the two climate policy scenarios are as follows: (1) business as usual (BAU) with no emission control on GHGs (including non-carbon GHGs), and (2) 2DC. Cumulative emissions in 2DC are zero over the computational time horizon (from 2010 to 2150) allowing the offset of peak-out and negative emissions by giving the 350-ppm scenario cumulative emission constraints from Wigley-Richels-Edmonds (WRE) (Wigley et al. 1996). We cap the cumulative emissions constraints in the 2DC and restart for computation from BAU solutions in each energy scenario. We computed these model runs for this study according to modifications in the model framework. Figure 2 depicts the inventories that were endogenously generated from our three resource balance models interlinked with the simplified climate model. Earlier LCA+ESM studies provide inventories of coal, NOx, and SOx in Fig. 2b. None of the other inventories in Fig. 2 has appeared in earlier studies, and they are included for the first time in the current study, which is one of its novel aspects. Figure 2a indicates the climate-related inventories (TR, SR) for global warming that was computed ex-post by GHG emissions from the three resource balance models. Figure 2b shows the inventories (m, w, p) endogenously computed from the models for mineral balance and energy systems plus exogenously given ODS scenarios (o). The left bar denotes the inventories for m (production of coal and limestone) and w (waste from concrete and metals) corresponding to the impact categories for mineral resources and waste, whereas the right bar indicates p (the pollution generated by NOx and SOx) for acidification and local air pollution, and o (ODS, negligible) for ozone depletion. Figure 2c represents inventory u, the area of land-use patterns with the total amount (left bar) and areal changes from the 20 patterns in land-use change (right bar).

Inventories
Our energy scenario is based on Ren100; hence, in Fig.2a, TR under BAU increases to 2.5 DC and is simultaneously stabilized at 0.55 m; further, Fig.2b depicts the lower rate of production of coal and limestone and shows a dramatically decreasing trend for the production rates of NOx and SOx. Figure 2c shows the expansion of forestland to absorb atmospheric CO 2 or of biomass energy with carbon capture and storage to meet the 2DC target. We input all these inventories into the equations listed in Sect. 2 and Appendix 2 to derive endpoint damages.  Figure 3a is dominated by global warming (yellow), followed by local air pollution (orange) due to the emission of SOx from coal production and negligible ozone layer depletion. We exclude particulate matter, such as PM 10 and PM 2.5 , since these cannot be endogenously estimated using our energy model.

Endpoint damages
Undernutrition caused by global warming has the largest impact since it has the highest baseline mortality rate. It is followed by malaria, coastal floods, and diarrhea. The overall damages seem almost stable over the timeline; regardless of the increase in temperature, the implication of damage per degree rise in temperature [DALY/deg C] shows a declining trend. One of the two reasons for the declining trend in the baseline mortality rate from onesixth to one-half is the increase in per capita income; thus, the values for malnutrition, diarrhea, and malaria are close to zero. The other reason is that the percentage in the total population is assumed to follow a dramatic declining trend in the SSP scenario for malnutrition and diarrhea, currently declining to 5% and 20% from 20% and 50%, respectively, in the largest region.
Impacts other than global warming are determined by inventory levels, rather than damage factors; the emission level of ODS is 1/20 of SOx and NOx. The impacts on local air pollution decline due to the phasing out of coal in the energy mix. A comparison of (a) global warming only (yellow parts), which is treatable according to existing LCA+ESM studies, and (b) the overall impacts shown in Fig. 3 revealed that the level of overall impacts declines (1) over time and (2) in 2DC Fig. 2 Inventories endogenously generated from our IAM. a Climaterelated inventories for increases in temperature (TR; blue, left) and sea level (SR; pink, right) for both climate policy scenarios (BAU; dashed, 2DC; solid) generated from the simplified climate model by inputting the endogenously calculated greenhouse gas emissions of both CO 2 and non-CO 2 using the energy systems model. These are not considered by existing LCA+ESMs studies. b Inventories of resource production and discharged waste (left axis); consequential emissions of ODS, NOx, and SOx (right axis) endogenously generated from the models for mineral balance and energy systems; and exogenously provided ODS scenario.
Unlike existing LCA+ESM studies, which consider energy-related inventories (i.e., coal, NOx, and SOx) alone, our model includes the mineral resource model and expands other inventories, some of which are significantly high levels of inventories. c Areas in land use (left solid bar and axis) and land-use changes (right hatched) generated using bio-resources and the simplified land-use model. Three land-use patterns (from bottom to top: forestry, grassland, and arable land) are shown, which are the main causes of land-related damages. All 20 patterns are illustrated for land-use change, among which forest is dominant. These land-related inventories have not been mentioned by existing LCA+ESM studies compared with BAU due to the significant reduction in local air pollution in the beginning and the slight increase in global warming later. This observation is a novel aspect of our study and is obtained by the addition of models, other than ESMs.
In Fig. 3b, the bar on the far left of the graph presents natural resources damages in monetary terms, including the resource rental by user cost (left axis), whereas the damages to resources in agriculture, land, and energy in terms of ambience (from left to right) are in physical units (along the right axis). The figure shows that the metrics are dominated by blue (damage to the land resource from global warming), followed by local air pollution and energy, whereas the others are invisible. The damage to land resource types (terrestrial or wetland) in monetary units is comparative because the larger proportion of submerged areas in terrestrial land has lower land cost. The shadow prices for copper and energy for heating are sufficiently high to reflect large-scale resource rental and energy resource damages, respectively, whereas the prices of On comparing (a) the energy related only (global warming and local air pollution, in blue, green, and yellow) by ESMs with (b) the others (ozone layer depletion and mining, in red) shown in this figure, we observe BAU and 2DC and clarify that (1) energy related shows little difference in the four columns, whereas (2) "other than local air pollution" (in red) in the "all" column (the extreme left) indicates moderate reduction. It is noted that independent studies on the impact assessment for each sector (i.e., agriculture, land, and energy) have been conducted on different and inconsistent backgrounds. In contrast, our results are consistent because they are generated by our own modeling, which is a strength of our study. Figure 3c depicts all the biodiversity damages, most of them dominated by land-use changes (in green) and some by global warming (in yellow). Damage from global warming increases in accordance with TR. The LU from grassland to forest, arable land, and other land types has the largest impact since it represents the largest areal transformation and highest damage factors. Comparing (a) global warming only (in yellow) with (b) other than global warming (LC in green) in this figure, we observe that the overall impacts strengthen over time and 2DC is slightly smaller than BAU in global warming only, whereas there is a slightly larger difference between 2DC and BAU after 2060 regardless of the changes in the breakdown of LC (green). The addition of the green items in this figure is a novel aspect of our IAM. Figure 3d displays all the damages from primary productivity in the left bar (including LU, left axis); they are dominated by land use, particularly in urban land areas (human settlements) and grasslands because of the large area in use. The value highlighted in the darkest green color at the top of the left bar is magnified in the right bar, excluding land use (excluding LU, right axis), where positive and negative aspects generally correspond to other (unmanaged land) to urban (human settlement) and loss of grassland due to larger level damage factors. It is apparent that the energy-related impacts of ESMs (i.e., impacts of SOx and NOx on local air pollution; denoted in orange) are negligible, whereas all the other impacts are from the newly added models in our IAM.

Interpretation of our results in relation to extant literature
To the best of our knowledge, no other study is directly comparable with ours. Therefore, it is important that we present the plausibility of our results. In this section, we interpret our results by comparing them with the general conclusions of other studies. We follow three endpoint results (Fig. 3): breakdown of human health, biodiversity, and comparison with other LCIA models.
The World Health Organization (WHO) conducts distinct studies on human health. Its Global Health Risks (World Health Organization 2009) considers environmental problems such as unsafe water, sanitation, hygiene, indoor smoke from solid fuels, lead exposure, urban outdoor air pollution, and global climate change. The World Health Organization (2009) predicts that the impacts of global warming and local air pollution on human health are comparable, and this prediction is similar to our findings. These impacts are highlighted by the World Health Organization (2014)  Next, we consider biodiversity. In The Millennium Ecosystem Service (MEA 2006) and The Economics of Ecosystems and Biodiversity (TEEB 2012) reports, the main impacts on biodiversity are from land use and global warming. Many studies on biodiversity impacts are an assessment of extinction speed, while relatively few examine the LCA. The global number of species extinction is reported to be in the tens of thousands (Pereira et al. 2010), while our results show 200 to 400 species. Our breakdown shows that this extinction is the result of land use, which has a greater impact than global warming. This confirms the plausibility of our result. In the field of LCA, the LC-IMPACT (2019) model treats acidification and water resource for biodiversity impacts, which are beyond the scope of our modeling framework.
Finally, we describe the dependency of our results on the LCIA model. The results for the endpoints of our study are based on LIME, whose impacts must be similar when applying other models (e.g., ReCiPe (2018) or LC-IMPACT), since the three models have similar damage coefficients. Human health impacts from global warming in LIME (Tang et al. 2019) are of a similar order to those found by De Schryver et al. (2009), who employed ReCiPe and LC-IMPACT. Regarding biodiversity impacts arising from global warming, LIME only treats plant species, while ReCiPe includes multiple taxon. However, impacts from plant species alone are of a comparable order in both models (Tang et al., under review). Although results on the four endpoint impacts seem to change to some extent if impact assessment methods (e.g., ReCiPe (2018) or LC-IMPACT (2019)) are applied instead of LIME, we expect that the scenario comparison of BAU and 2DC are change to a limited degree.

Novel aspects of this study
Here, we summarize the original aspects of this study in terms of the study's methodology, focus on LCA, and results sourced from the methodology.

Methodology
The key methodological originality of this study, compared with existing LCA+ESM studies, is that it extended the framework comprehensively (see Fig. 1) to include the (1) the simplified climate model, (2) the mineral resources model, (3) the simplified land-use model, (4) endpoint impacts by LIME, and (5) a hard-linked characteristic. Many of the existing studies employed ESMs such as TIMES/TIAM model families (see the review by Volkart 2017), which are not linked to the simplified climate model to compute rises in temperature and sea level; hence, GHG emissions from ESMs were combined with the LCI dataset (e.g., ecoinvent) in some studies (e.g., Volkart et al. 2018;Boubault et al. 2018) and/or with the LCIA model (e.g., ReCiPe) in others (e.g., García-Gusano et al. 2017a, b;Pehl et al. 2018).

Focus on lifecycle assessment
The studies discussed in Sect. 4.2.1 mostly focus on power generation technologies alone, whereas our study considers the overall GHG emissions from energy systems. A model for land use and land-use change has only been considered by Pehl et al. (2018), and it uses MAgPIE in a soft-linked manner; however, the authors consider GHG emissions alone in their assessment. Contrastingly, our study evaluates endpoint impacts on biodiversity and primary productivity in a hard-linked manner. Further, other studies did not use the mineral resource model, whereas our study is unique in that the mineral model (Tokimatsu et al. 2017b) is used in our framework.
The focus on LCA varies among studies. Some emphasize the cradle-to-grave nature of LCA, as well as indirect emissions in plant construction or embedded emissions during energy resource production and transport in LCI and both methodological investigation and model analysis (e.g., Pehl et al. 2018;Arvesen et al. 2018;Volkart et al. 2018), whereas other studies focus on the "outcomes" of using LCI or LCIA, such as metal requirements (e.g., Hertwich et al. 2015;Boubault et al. 2018) and endpoint impacts (e.g., García-Gusano et al. 2017a, b). Unlike the studies by Arvesen et al. (2018) and Weidema et al. (2018), this study does not focus on linking LCI and LCIA in mid-point damages.
By applying another impact assessment method, midpoint damages can be examined using more than 10 environmental categories (global warming, acidification, eutrophication, etc.; see García-Gusano et al. 2017a, b) in an uncomparable manner. Hence, it is impossible to compare scenarios by mid-point in a quantitative manner. On the contrary, it is easier to interpret results and compare various environmental issues by using an assessment of four endpoint damages, which was the method used in this study. One advantage to indicating the results using a single index by the weighting factors for four endpoints is that the cost-benefit can be expressed in monetary units. However, an inherent disadvantage is that this method has significant uncertainty caused by the (i) aggregation of multiple factors into a single index; (ii) omission of detailed and concrete environmental information; and (iii) use of methods such as interviews or surveys of the general public, which are more subjective approaches than natural science methods. Hence, showing results on four endpoints is a reasonable alternative.

Results
The originality of our results is based on the extension of our model framework. Most of the existing LCA+ ESM studies reveal results on the power generation technologies of their life stages (i.e., construction, operation, and decommission) related to their (in)direct LCI inventories (e.g., GHG emissions and metal requirement), whereas our study encompasses highly expanded and general inventories that are multiplied by damage factors modified from LIME3. Hence, the only overwrapping inventories are those from ESM (e.g., coal production and emissions of SOx and NOx); the remaining inventories, such as increases in temperature and sea level, material-related inventories (production of limestone and waste discharge), and land-use changes, are left untreated in existing studies. Only a few studies (e.g., García-Gusano et al. 2017a, b) examine the overwrapping endpoint impacts on human health and ecosystem services; contrarily, our study focuses more on natural resources (Fig. 3b) and primary productivity (Fig. 3d).
Moreover, apart from LCA+ESM studies, many publications on ESMs/IAMs investigate the co-benefit of climate change mitigation efforts on local air pollution (e.g., Bollen et al. 2010). They focused on reductions in GHG, SOx, and NOx emissions, and their cost-benefit provided no results in terms of endpoints for human health damage using common units (e.g., DALY; see Fig. 3a). Therefore, our results are not at all similar to the results of such co-benefit studies.
Furthermore, most studies on climate change mitigation using ESMs/IAMs focused on policy assessment through scenario comparisons of economic efficiency or cost-benefit. Our results are based on not only climate change mitigation efforts and their co-benefit but also other impact categories, in terms of endpoint damages, which clarify that the impacts of mining, land use, and land-use change are significant, as well. Specifically, the 2DC scenario in this study, which heavily depended on renewables, promoted forestry and focused on trade-off by expending these efforts. We recommend that climate change mitigation efforts be assessed according to not only the efforts' direct effects but also more comprehensive perspectives on environmental issues.

Question of balance regarding our Dose-Response Functions approach
To apply the cradle-to-grave philosophy, a question of balance should always be addressed while margining large models to investigate the dynamic mode of LCA. We developed the damage factors in LCIA by adding the consequential damages divided by the released inventories (causes). We estimated the damages using detailed, specific models to conduct relevant impact assessments. The damages retained short-to long-term horizons depending on the impacts. Consequently, the damage factors intentionally eliminated time evolution to provide "averaged," aggregated damages per unit of inventory; thus, LCA practitioners can easily calculate damages using basic arithmetic. This is an extremely simplified aspect of our model for practitioners, whereas studies such as Arvesen et al. (2018) and Volkart et al. (2018) present much more complicated model.
LCIA models consistently provide multiple impact pathways in the linear form to assess potential impacts aggregated over time and space, instead of the actual impacts of a specific source on a specific receptor (endpoint). This is very different from the so-called damage functions in climate-welfare economy IAMs. IAMs with no specific relationships between the actual impacts of a specific source on a specific endpoint simply link the economic damage in GDP loss with a rise in temperature. Intrinsically, the investigation in the dynamic mode of LCA proceeds in the opposite direction such that it disaggregates the impacts. This disaggregation with time expansion is plausible, suitable, and compatible for long-lived, stock-based phenomena, such as global warming, whereas we do not expand on other impacts in our model and they remain unchanged to apply the factors. This is because our model has a 10-year interval, which is sufficient to approximate the consequential impacts in practical LCA exercises.
Another dimension, besides time, is the spatial expansion that we explain here. The development of the damage factors covering 197 countries for the LIME3 also required huge resources (e.g., time, people, and budget). Our approach for the category impacts using LIME2 seems somewhat rough or adhoc, as we only adjust the levels of the endpoints by order of magnitude. The next step seems to be to adjust not only the endpoints but also make an additional amendment to the DRF or develop country-level damage functions, as in LIME3.
Skipping such strategies otherwise requires tremendous resources to operate all the models related to the impacts via soft-linkage, after blending the underlying socio-economic, energy, and climate policies. Hence, our adjustment of DRF in LCIA to be compatible with IAM is a well-balanced approach.

Strengths and weaknesses of our approach
Overall, our approach can provide more unique insights than those in "conventional" and "mainstreaming" approaches in both disciplines of LCA and IAM. In addition to the originality and comprehensiveness of this study, as described in Sect. 1.2, this section highlights the strengths and weaknesses of our approach, mainly sourced from the LCIA modeling framework.
First, the strength of many European studies is in their technology assessments, in particular, in power generation analysis based on the New Energy Externalities Development for Sustainability (NEEDS) project (Schenler et al. 2008). This approach demonstrates the three pillars of sustainability-social, techno-economic, and environmentin detail, for power generation technologies. On the other hand, our Japanese approach has advantages in terms of comprehensive coverage of environmental impacts and resources and the weighting process based on LIME.
Second, we compute environmental damages on the four endpoints within the framework of LIME, consistent with scenarios in the model, although sectoral impact assessments (e.g., human health, agriculture, biodiversity) are more comprehensive. For example, comprehensive human health damage includes various causes unrelated to the environment, such as tobacco or alcohol use, which we exclude. Our analysis shows the sum of the impacts of both local air pollution and global warming over time under the same energy and climate scenarios, while sector impact assessments are conducted separately under uncommon scenarios and assumptions. The same elements correspond to other impact categories; hence, an advantage of our analysis is in its consistent assessment within (limited) environmentrelated damages under common and flexible scenarios and assumptions.
Regarding strengths and weaknesses, the former (strength) is available and traceable to the open-public method we developed in this research for modification of DRFs to apply open-public damage factors in LIME. Hence, our modified damage factors can be reproduced, even rearranged, to apply ESMs/IAMs or inventory scenarios borrowed from other studies. The weakness of our approach is its complexity, resource-and time-intensiveness, and closed-source nature, particularly as Vanderpaer and Gibon (2018) argue that the importance of views by and usefulness for LCA practitioners, open-source software, and data sharing are necessary for such prospective LCA studies.
Another weakness of our analysis is that it uses no common assumptions with the other "comprehensive" and "sectoral detailed" analyses. Thus, our results will be difficult to validate or compare with similar or related studies. For example, some studies investigate 100% renewable energy with the 2DC target settings investigated here. Moreover, our methodology for impact assessment, DRF, seems to be similar to those employed in FUND; however, distinctive differences in the LIME framework separate the damage assessment in physical terms from the monetary valuation, by employing the weighting process (i.e., damage assessment in monetary terms). However, in contrast to this study, the FUND model incorporates highly aggregated damage functions in monetary terms, neither at the endpoints nor in physical units.

Co-benefit of climate policy: outcome of interlinking models of lifecycle assessment and integrated assessment
This article focuses on integrating models of integrated assessment and lifecycle assessment. This section describes another outcome in addition to that in Sect. 4.3 (specifically, dynamic impact assessment). The other outcome is to show the cobenefit from global warming mitigation. We prepared two sets of climate policies, BAU and 2DC, using the cumulative emission cap. This results in transformation of energy for reducing GHG emissions and land-use changes for absorbing and releasing carbon through forestry. Reduction in GHG emissions lowers the coal consumption rate, which, in turn, reduces the SOx and NOx emissions. The decrease in coal consumption rate then brings co-benefits in addition to global warming for greater reduction in human health impacts in 2DC than in BAU (see Fig. 3a). The increase in the forestry area implies that the transformation of other land-use patterns is suppressed, leading to relative reduction in biodiversity (Fig  3c). On the contrary, impacts on natural resource (Fig. 3b) all relate to global warming; hence, we can expect no co-benefit.

Conclusions
This study illustrates the future paths of the four endpoint damages on a global level over this century, caused by various impact categories endogenously tied to our hard-linked, fully Integrated Assessment Model. This study describes the DRFs while excluding the weighting process in the LCIA framework. Further insights, not discussed here, can be identified, especially in terms of sustainable development aspects. We leave these to the other studies mentioned in the introduction section. The studies in the XLCA modes (i.e., LCA and ESMs, IAMs) promise that they can provide unforeseeable and fruitful insights.
In recent decades, there have been suggestions to use the models for LCA and systems analysis (Jolliet et al. 2004); however, this has led to some controversy, especially in the weighting process in the LCIA framework (Hordijk 1995). Stimulated by such controversy, we have pioneered studies in this area. Currently, the situation seems to have dramatically changed with many studies published, triggered by the NEEDS project, which is expected to widen the XLCA modes. The strength of the European studies lies in technology assessment, especially power generation based on the NEEDS project, while our Japanese studies have advantages in their comprehensive coverage of environmental impacts and resources and the weighting process based on LIME.
LIME3 is a "modeling methodology," not software, and it provides "lists of damage factors, weighting factors, and economic factors" and is downloadable from . Adding another category requires further investigations based on the modeling framework; hence, in this regard, it is inflexible. However, this article described the precise details of the application to our IAM, which can also be useful to any ESM/IAM. The ESM side requires an open source tool (e.g., OSeMOSYS) for better researcher engagement. Such trends seem to have accelerated over the past several years after being handled by closed scientific communities for decades. The authors believe XLCA modes studies may significantly contribute to expanding our understanding of society. We strongly encourage research scientists in both disciplines to participate in studies on the XLCA modes.
on Japanese conditions and modeling; however, the latest version (LIME3; LCA Society of Japan 2019;  was developed on a global scale. Details of the model developments via the three steps can be found in the literature cited elsewhere in this study, as such descriptions cover several hundred pages and are outside the scope in this study, especially LIME2 (Itsubo and Inaba 2012, 2013. LIME3 enables the comparison and aggregation of environmental impacts among 193 countries on a global scale. These include global warming, air pollution, photochemical ozone creation, water use, fossil fuel, mineral resources, and forest resources tied to land use (those underlined are included in our model here). We calculated the damage factors, connecting potential damages of the four endpoints with a unit amount of inventory release on the country-to-country basis of 193 countries, using the dose response function (DRF).

LIME3 in our IAM
In this study, damage by global warming is based on LIME3 for human health and biodiversity, including further investigations and modifications by Tang et al. (2019), while natural resources are sourced from LIME2 (Itsubo and Inaba 2012). The damage caused from CO2 emissions was calculated for each disease (malnutrition, diarrhea, cardiovascular disease, malaria, and coastal and inland flooding) (Tang et al. 2018a;Tang et al. 2019). Biodiversity damage from climate change was assessed by statistical analysis of vegetation distribution changes of 250 species of plants in Japan (Tang et al. 2018c); however, Tang et al. (under review) have expanded the grid data used to over 280 thousand taxa (see Appendix 2).
Damages from global warming to crop production (rice, wheat, and corn) and energy resource consumption is categorized by ambience (cooling or heating) in terms of global mean temperature increase and for land resources, in terms of submerged areas (onshore, wetland, tundra, and desert) by increase in sea level. The damage to resources (crops, energy) is valued using shadow prices obtained in our three resource balance models, while those for areas are referenced from Fankhauser (1998).
Damages from land use and land-use change to biodiversity in LIME3 were assessed using the incremental risk of extinction of 2000 species of vascular plants in Japan, as per each type of land use and transformation (Yamaguchi et al. 2018). We adopted Yamaguchi et al.'s (2018) methodology for 193 countries, using International Union for Conservation of Nature and Natural Resources (IUCN) data. Land-use damage to primary productivity is measured as the differential between potential NPP and current NPP, using the landuse pattern at the country level, which is disaggregated into the country level from those at the global level (Haberl et al. 2007) after identifying overlapping global vegetation patterns and national borders. LIME provides the following definitions for this impact category. Inventories in land use and land-use change are areas in occupation (or appropriation) in certain patterns (e.g., forestry) and in changes (e.g., from forest to human settlement). The time for appropriation of land use is assigned to every land-use pattern, while the time for restoration to arrive at the original productivity level is given in land-use change. Damages by these focused on degradation of potential NPP and increasing risk of extinction for vascular plants (disregarding numerous impacts associated due to a lack of qualitative information for impact assessment).
Damages by mineral resources in LIME3 are assessed for three endpoints (natural resources, biodiversity, and primary productivity). The natural resources are valued using the user cost method. The mining-transformed land area consists of the mine site, the deposit site from accumulated overburden discharged from mining, and the waste disposal site for tailings generated from ore concentration and coal pretreatment. This transformed area is formulated in terms of weight of mined resources, ore grade, ore body thickness, and the socalled hidden flow. The number of endangered species and primary productivity are estimated based on the transformed area. The above modeling is based on country and regional data aggregated from site-specific data such as mining method (underground, surface), ore body thickness, and ore grade.
Our modeling does not address the following damages in LIME3. The trans-boundary effects were taken into account for PM 2.5 , both primary and secondary produced air pollutants (i.e., SOx, NOx, black carbon, and organic carbon) (Tang et al. 2018d) and photochemical oxidant substances (Tang et al. 2018b) using a global chemical transport model. Agricultural irrigation water constrained by freshwater consumption was combined with malnutrition damage by crop production loss through irrigation, via global changes in food production and trade (Motoshita et al. 2018).

LIME2 in our IAM
The following damages included in our model are based on LIME2: waste (Itsubo and Inaba 2014), local air pollution and acidification (Itsubo and Inaba 2013), and ozone layer depletion (Itsubo and Inaba 2012). Damage factors for waste are only included in biodiversity and primary productivity at landfilling sites in the final disposal of metallic and cement waste. These are considered relatively common in countries other than Japan, although the damage factors provided for numerous types of waste in LIME2 are more countryspecific damage estimates. Discharged wastes are converted to area changes for the landfilling site, which is transformed into endpoint damages in a similar manner as mineral resources.
The following damages are included: local air pollution (on human health) by both SOx and NOx, and acidifications (on natural resources and primary productivity) by NOx. For local air pollution, respiratory disease is caused by an increase in pollutant concentration in SOx (immobile, point source from chimney) and NOx (both immobile and mobile from automobiles), both as primary and secondary pollution, through emissions. The concentration increase in emissions is computed for both short (0.2-20 km) and long (200-2000 km) distances with a source-receptor relationship (SRR).
For the acidification via acidic deposition from the air to the natural environment (terrestrial surface ground, and sea), terrestrial vegetation for primary production, and wood and fishery production (catch of Salmonidae) for natural resources are included in LIME2 (i.e., hard to quantify impacts, such as cultural assets, are excluded). Examples of main acid substances emitting hydrogen ion (H + ) are sulfuric acid (H 2 SO 4 ) and nitric acid (HNO 3 ), which are changed via an alteration process from emissions of SOx and NOx, respectively. Those substances are deposited into the natural environment via dry gaseous, particulates or wet precipitation, such as rain. An increase in acidic deposition via computed SRR is used to evaluate increasing endpoint impacts. Decrease in the primary production rate is through an increase in the soil solution Al 3+ caused by a pH decrease in soil solution, whereas the Salmonidae decrease is through the lowering of the intrinsic growth rate via a pH decrease in land water.
The damages from ODS in LIME2 are human health (skin cancer, cataract), natural resources (production in agriculture and wood), and primary productivity (both terrestrial and aquatic ecosystems). ODS emissions improve results, with an increase in UV-B reaching the ground by an ozone decrease in the stratosphere from a concentration increase in both the troposphere and stratosphere. Wavelength corrections of UV-B (human health (erythema), terrestrial plants (PAS), and phytoplankton (PhAS)) are used to decrease their global levels from latitude band levels combined with geographical distribution in the population and agricultural production. This, in turn, estimates their global level damages corrected by lifetime ODS atmosphere.

Excluded weighting process
The impact of weighting factor on each endpoint provides an environmental external cost expressed as a single index; we use a conjoint analysis to obtain the weighting factors for the endpoints. Therefore, we can calculate the monetary value of the environmental impact of one unit of each endpoint based on the environmental valuation techniques. Information on this can be found elsewhere Murakami et al. 2018;Tokimatsu et al. 2016a). As a caveat, we have excluded this weighting process from this study.
Appendix 2. Detailed mathematical formulations in the LIME framework applied on a global scale over time 1. Outline Table 1 provides a causal chain, DRF, connecting inventory (v), impact category (i), and endpoint (e) in our modeling framework. Appendix 4 provides lists of the meanings of subscripts and abbreviations. The check marks (✓) indicate damage factors (,) corresponding to the inventory, impact categories, and endpoints. The numbers in subscript (e.g., 7 for m) represent the number of minerals (i.e., bauxite, coal, copper, iron ore, lead, limestone, and zinc). Since MR contains three endpoints (NR, PP, and BD), the number of damage coefficient (DC) categorized in MR equals 21, provided from LIME.
This appendix provides detailed formulations of DRF and the aggregated form. Details of the development of damage function (DF) and DC in LIME are available in related studies.  [°C]. Both DPD and BMR are expressed as functions of y j,t, while RR is tied to TR t . These DPD, BMR, and RR are derived from regression analysis based on WHO data after rearranging the original WHO 21 regions for our model (10 regions). The formulation of DF and the regression analysis are based on past studies (Tang et al. 2018a), which were updated based on 2014 estimates from the WHO (Tang et al. 2019). RR, excluding malnutrition, is blended with the SSP scenario based on scenarios for the years 2030 and 2050 in the WHO (2014) report. RR for malnutrition is formulated as a step function using those scenario levels and BMR in 2030 and 2050. BMR in coastal flooding areas is given a fixed number commonly used in scenarios. DPD is considered the percentage of the population with the target diseases. We express the impacts on agricultural products (i.e., rice, corn, and wheat) by multiplying damage factors with global temperature rise, as follows: where f food (rice, corn, and wheat) k the three sub categories in NR (f,l,r) ASP f,j,t [USD/Million ton]: shadow prices of the regional agricultural products at time t AGP f [Million ton/yr/°C]: parameters for the damages in global level per unit temperature rise RAP f,j [%] parameters to allocate the global damages to region j The ASP f,j,t is derived from marginal changes in land-use cost divided by those in production.
The AGP f is derived from changes in agricultural production between 1990 and when CO 2 doubled, divid- Fig. 4 Overall structure of LIME ed by the global mean temperature rise between the years.
The RAP f,j is derived from the regional products differential between 1990 and when CO 2 doubled, divided by the temperature rise.

Land resources by sea level rise
The damage to land resources is expressed by submerged area in terms of the unit rise of the sea level, multiplied by the height of the rise in sea level and the areal value, expressed as follows. Unlike the two other categories (i.e., agricultural products and energy), the constant value of the unit price of the land area is applied instead of the shadow value derived from our model.
where l land, tundra, desert, and wetland LAP l,j [USD/km 2 ]: unit price of land area SMA l [km 2 /m]: submerged land area in a unit rise of sea level, derived from gridded data for coastal lines above 0.5 m, from the sea level; it excludes those with high population density SR t [m]: height of sea level rise at time t

Energy resources
Impact on energy resources is measured by the economic damage of energy, expressed by the following equation. The first term (SPR) is the shadow price of energy consumed in heating and cooling. The second term is the adjustment by the relative increase in "degree days" (in both heating and cooling) from 1990 multiplied by the third term, per capita energy consumption (RCC) and the population number (N). The fourth term is adjustments for energy consumption, both in ratios of degree days and heat demand per GDP intensity per capita in regions and time, compared with Japan in 1990.
The global mean temperature rise is reflected in the degree day increase multiplied by its unit temperature rise, as in Eq. (11). A simplified regression analysis on energy consumption per capita for cooling and heating with regard to per capita GDP suggested that we apply the two patterns of functions, as in Eqs. (12,13). The damage factors in BD are sourced from Tang et al. (under review), which was revised from lists in LIME3 for the biodiversity impact of global warming (DC j,i = GW,e = BD ) [EINES/yr] in a grid basis (approximately 10 km 2 ), with the regional damage factors aggregated. DC j,i = GW,e = BD is developed from patterns from five global climate models (GCMs), two with/without plant species mobility and three plant taxa (gymnosperms, angiosperms, and pteridophytes with numbers 1052, 268000, and 12000, respectively).
The factor is estimated from the difference in the reciprocal of time to extinction averaged by the number of species divided by temperature difference, then multiplied and summed with the numbers of the three taxa (i.e., 281,052). The differences are between those levels in RCP 8.5 and in RCP 2.6 over 100 years . The time to extinction is determined so that a potential habitat area in 1970 is extinguished (assuming it reaches 1 km 2 ) by a speed with ratio between the areas in 2070 and 1970.
3. Other than global warming: Adjustment for BD by derivations of B(β) and S(α j ) The derivation of the damage factor increase (B(β)) is as follows. Through a simple regression analysis, we found a good correlation between E and V in vascular plants from the Japanese Red Data Book, where E indicator to express species extinction, where E = ∑ s (1/T a, s − 1/T b, s ), corresponding to the sum of inverse extinction time (T) differential before and after losing one species. V number of endangered vascular plant species Then e(β), the ratio of E 2 and E 1 in the current and original state, is expressed using habitat area R and its areal change ratio β, since V can be considered proportional to R.
where β = ΔR 1 /R 2 , R 2 = R 1 − ΔR 1 , V 2 = V 1 − ΔV 1 Using the above, B(β) can be derived from an integral from 0 to β, as follows. This function is applicable in β j,t (areal changes in forestry plus grassland with time evolution in region j) and β(TR j,t ).
Next, we defined S(α j ) as a magnifying term to adjust between Japan and the target region, regarding the increase in extinction risk of species in which α is defined as the ratio of population density of endangered species in the target country, divided by the Japan values. The function S(α j ) is calibrated using data for Japanese soil and its gravel industry, from which 1.5 and 4.39 are sourced, and a, b are obtained by the regression.
Appendix 3. The three resource balance models

Outline
The three resource balance models use dynamic linear programming at the global level to equate supply and demand for resources. The model includes production, transportation, conversion, final production, disposal, and recycling. The energy model (Tokimatsu et al. 2016c) and the non-fuel mineral resources balance model (Tokimatsu et al. 2017b), the landuse model for biomass and food (Tokimatsu et al. 2017a) calculate the discounted sum of supply costs to determine the structure of demand and supply balance of fuel minerals (i.e., energy) and non-fuel minerals.
The manufacturing sectors of the global economy are proxied using three representative manufacturing sectors: electricity and machinery, construction and building and automobile industry. Non-fuel materials exist as in-use stocks of product goods during the assumed products lifetime, after which they are out-of-use stocks, and finally disposed or recycled.
The land-use model (Tokimatsu et al. 2017a) considers five categories of land use (forestry, grassland, cropland, urban, others) and 20 possible land-use changes between uses. Landuse change is generated through evolutions in demand (food, biomass resources, and urban land use) and supply (land rent, land conversion, biomass and food production) of land.
The three models are integrated (Fig. 5) based on a detailed set of technology options, concerning fuel and non-fuel mineral resources, land use and land-use change. The three resource balance models use dynamic linear programming at the global level to equate supply and demand for resources. The model includes production, transportation, conversion, final production, disposal, and recycling. The left side of Fig. 5 indicates resource supply, while the right side shows the demand side and end-use products and waste disposal. The upper section illustrates mineral and material flows; the lower section shows biomass and food flows, and the middle section shows energy flows.

Direct cost formulation from the resource balance models
The supply costs are generated by mining, milling, dressing, smelting, and refining for non-fuel minerals, the conversion process of electrical or chemical energy, and transportation among the ten global regions to final demand. Supply costs TC is divided into energy systems cost EC, non-fuel energy resource cost NFC, and land and agriculture related cost LC. Therefore, the total cost TC is  Annual power generation of power plant PWh/yr subscripts rg, rg1 The 10 regional allocations yr, yr1 The time step gr The grade of the resources ie The primary energy and the secondary energy ip The energy conversion plants ih The heat use technologies id The sectoral heat demand iv The battery types (in transport sector) it The load duration pattern (base, middle, peak) is The final energy demand (power, heat, transport) ic The carbon capture and storage technologies degree day per unit temperature rise, derived from its rise in TR t = 1.5°C LAP l,j [USD/km 2 ] unit price of land area RAP f,j [%] parameters to allocate the global damages to region j in agricultural products SMA l [km 2 /m] submerged land area in a unit rise of sea level, derived from gridded data for coastal lines above 0.5 m height, from sea level, excluding those with high population density α ratio of population density of endangered species in the target country, divided by the Japan values 3. Subscripts a, b, c coefficients used in functions applied by regression analysis d disease (malaria, diarrhea, malnutrition, cardiovascular diseases, coastal flooding and inland flooding) e endpoints (human health (HH), natural resource (NR), primary productivity (PP), and biodiversity (BD)) f foods (rice, corn, and wheat) a magnifying term to adjust between Japan and the target region for the increase in extinction risk of species SPR r,j,t [USD/toe] shadow price of energy resources derived from marginal changes in the energy systems cost, divided by changes in energy resource consumption SR t [m] height of sea level rise at time t t 10-year time step (t = 2010, 2020, …, 2150) TR t [°C] global mean temperature rise at time t V number of endangered vascular plant species y j,t [USD/cap] per capita GDP at time t β ratio of land area change in forestry plus grassland Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.