The Intergovernmental Panel on Climate Change (IPCC) reviewed many integrated assessment models (IAMs) to evaluate the impacts of climate change policies (IPCC  ). Although many researchers have begun to measure the impacts of mitigating greenhouse gas (GHG) emissions for the 1.5˚C and/or 2˚C target in the Paris Agreement using IAMs, there are few studies focusing on adaptation under the Paris Agreement. One reason is that it is difficult to measure the cost of adaptation to climate change by IAMs. However, some studies analyze adaptation to climate change by global economic models in other situations. Therefore, it is important to clarify regional adaptation and economic effects of the long-term temperature targets of the Paris Agreement and its climate policies in the world.
Latest IAMs research on climate change applies two scenarios for calculating impacts of climate change: Shared Socio-Economic Pathways (SSPs) and Representative Concentration Pathways (RCPs). The SSPs consist of five socio-economic scenarios with mitigation and adaptation to climate change and include GDP, population and GHG emissions scenarios. Each SSP scenario represents a possible societal future. The RCPs consist of four climatic scenarios on radiative forcing, and each RCP corresponds to a GHG mitigation policy. Integrated Assessment Modeling Consortium (IAMC), an IAM community, suggests that IAM researchers apply combined SSP-RCP scenarios to make it easier to assess literatures within and across research communities as well as across studies at different scales and in different regions (van Vuuren et al.  ). In this study, we measure regional impacts of adaptation to climate change for the Paris Agreement under the SSPs and RCPs combination scenarios using a global economic model in which each region has multiple economic sectors. We first estimate the parameters of our Evaluation Model for Environmental Damage and Adaption (EMEDA) on adaptation under various scenarios to decompose the total costs and damages of climate change into climate change damages and adaptation costs. Second, we simulate EMEDA to obtain optimal adaptation levels and economic impacts of climate change under the SSP-RCP combination scenarios for the 2.0˚C target. Finally, we discuss the 1.5˚C target case which has lower climate impacts and aggressive adaptation.
By calculating the global economic benefits and costs of aggressive adaptations to global warming using the PAGE model, Hope et al.  showed that aggressive adaptation was optimal. Rosenzweig and Parry  dealt with the adaptation benefits and/or costs of climate change in the agricultural sector. Later, AD-DICE (de Bruin et al.  ) was developed as a one-region world economic model in which the total costs of climate change are decomposed into adaptation costs and residual costs. Based on AD-DICE, AD-RICE model (de Bruin et al.  ) was developed as a multi-regional model with one sector per region. In this research, we use the dynamic EMEDA since it includes multiple sectors per region (Sakaue et al.  ).
The dynamic EMEDA consists of eight regions, each of which has each eight sectors (see Table 1). OECD8 consists of EU countries, Canada, Australia and New Zealand, and FSU_EEurope consists of the Former Soviet Union (FSU) and other European countries. OAsiaOceania includes Asian and Oceanian countries while OAmerica represents the countries of Latin America. General economic data and GHG data in the year 2004 are from GTAP7. To calculate coefficients of change in CO2 emissions, results from the Asia-Pacific Integrated Model
Table 1. Regions and sectors considered in dynamic EMEDA based on GTAP7.
(AIM) are used (ICA-RUS  ). GDP growth and population growth until the year 2100 under each SSP scenario are from the SSP database (IIASA  ). CO2 emissions until the year 2100 under RCP scenarios are from the RCP database (IIASA  ).
In dynamic EMEDA, the value-added function with damage of global warming is as:
where is a sector, is a region, and is time. is a value-added production function without global warming damages, is capital, and is labor. is a mitigation cost as a percent of GDP and is a sea level rise (SLR) and non-SLR damage function modified from RICE 2010 (Nordhaus  ). A mitigation cost as a percent of GDP is as:
where is the rate of CO2 emissions decline, is a CO2 emissions coefficient and is an EMEDA parameter of abatement cost functions. Damage function is as:
where is the rise in atmospheric temperature (in ˚C compared to the year 1900), is the level of SLR caused by temperature rise (in meters compared to the year 2000), and , , and are parameters.
Following AD-DICE (de Bruin et al.  ) and AD-RICE (de Bruin et al.  ) models, modified damage function for adaptation is as:
where is adaptation level, is gross damage of global warming as a percent of GDP and is adaptation cost as a percent of GDP. Gross damage of global warming as a percent of GDP is as:
where , , , , and are parameters. These are natural extensions of the damage function in RICE 2010. Based on AD-DICE, adaptation cost as a percent of GDP is as:
where and are parameters. In this study, we set the SSP2, “middle-load” socio-economic scenario, as our standard scenario to simulate SSP-RCP combinations. We estimate all parameters, , , , , , , and , by minimizing the sum of discounted squared errors where the discount rate is three percent in the SSP2. Following de Bruin et al.  and Hsiang and Narita  , we adopt several constraints for and : 1) ; 2) ; 3) ; and 4) . Table 2 shows the estimated parameters in the SSP2 scenario. Since the challenge of adaptation is different under SSP scenarios, we assume that reflects the difference in adaptation technology of various SSP scenarios. For example, is reduced (increased) by 30% in SSP1 (SSP3) compared to SSP2.
Table 2. Estimated parameters on adaptation for the SSP2.
We analyze three SSP and two RCP combination scenarios with business as usual (BaU) scenario, a reference scenario in which no region reduces CO2 emissions: SSP1-RCP2.6, SSP1-RCP4.5, SSP1-BaU, SSP2-RCP2.6, SSP2-RCP4.5, SSP2-BaU, SSP3-RCP2.6, SSP3-RCP4.5 and SSP3-BaU. Then, we discuss the 1.5˚C target case under SSP2 scenario, which is SSP2-1.5˚C.
The path of temperature rise in SSP1-BaU, SSP2-BaU, SSP3-BaU, SSP2-RCP2.6 and SSP2-RCP4.5 are shown in Figure 1 left. The highest temperature rise scenario is SSP3-BaU, in which temperature rises by about 4˚C in the year 2100. In the case of SSP2-RCP2.6, temperature rises are moderate with an increase of around 2˚C. In SSP2-RCP4.5, temperature rise amounts to about 2.8˚C above 1900 temperatures by the year 2100. There is little difference in sea level rise across scenarios until about 2070, when changes of different magnitudes start to be seen (Figure 1, right).
3.1. Regional Optimal Adaptation Level
Each region executes adaptation for reducing the total damages from climate change. An optimal adaptation level of a region is defined as that which a region chooses for minimizing its total damages from climate change with respect to an adaptation level. Therefore, we obtain optimal adaptation levels by region at time by minimizing the modified damage function with respect to adaptation level (de Bruin et al.  ). SSP-RCP combination optimal adaptation levels in EMEDA are shown in Figure 2. Our model shows higher damages from climate change compared to AD-DICE. This is because we use climate module, damage function and cost function with the latest pathways.
In each scenario, adaptation level increases as temperature rises because of increment of gross damage of climate change. Focusing on RCP scenarios, the optimal adaptation level in RCP2.6 is lower than that of BaU, since global warming damages are lower in RCP2.6. In the year 2100, the optimal adaptation level ranges from 0.4 to 0.7 in the BaU, from 0.3 to 0.5 in the RCP2.6 and from 0.4 to 0.6 in the RCP4.5. This indicates that gross damages by global warming are
Figure 1. Left: temperature rise (˚C above 1900) and Right: sea level rise (meters above 2000).
Figure 2. Adaptation levels by SSPs-BaU, -RCP4.5 and -RCP2.6 scenarios.
reduced by approximately half though adaptation. Under BaU, the optimal adaptation level in SSP1 is lower than in SSP2 while higher in SSP3. However, RCP2.6 and RCP4.5 scenarios are less different among SSPs because reduction in CO2 emissions largely affects SSP scenarios.
Comparing regions, the optimal adaptation level in Africa is lower than other regions in the year 2100 under SSP2-BaU. Changes in adaptation levels are faster in Africa since they have rapid economic development and population growth.
OECD countries also have lower adaptation levels because the costs of adaptation are higher (Figure 3). Interestingly, the optimal adaptation level in China is higher in several scenarios since China suffers more damages from sea level rise.
3.2. Regional Adaptation Costs
In the results of adaptation costs, we mainly discuss the SSP2-combination scenarios: SSP2-BaU, SSP2-RCP4.5 and SSP2-RCP2.6. SSP2-RCP2.6 is the 2.0˚C target scenario achieving the CO2 emissions reductions proposed by the Paris Agreement, while SSP2-RCP4.5 is the most modest CO2 emissions reduction scenario among all SSP2-combination scenarios.
Regional adaptation costs in SSP-RCP4.5, SSP2-RCP2.6 and SSP2-BaU are shown in Figure 3. Large regional differences in adaptation costs can be found in SSP2-BaU. Adaptation costs in OAsiaOceania and Africa amount to more than one percent of GDP in the year 2100 while those in Japan, USA, OECD8, FSU_EEurope and Latin America are less than 0.5 percent of GDP. This indicates that adaptation costs in developing countries are higher than in other countries. This is consistent with adaptation funds for developing countries in the Paris Agreement (UNFCCC  ) and similar to AD-RICE results (de Bruin et al.  ). There are smaller regional differences in SSP2-RCP2.6. Chinese adaptation costs are higher than in other regions because of higher damages from sea level rise. In SSP2-RCP4.5, there are some differences among regions. Adaptation costs in Africa, China, and OAsiaOceania are higher than other regions while the costs in Africa, China, and OAsiaOceania are less than one percent of GDP in the year 2100.
Figure 4 shows global total damages and costs (TOTAL), mitigation cost (MTGC), adaptation cost (ADPC), gross damage without adaptation (GDMG), reduced damage with adaptation (RDMG) in SSP2-BaU, SSP2-RCP2.6 and SSP-RCP4.5, respectively. Reduced damage with adaptation is calculated by the product of gross damage and adaptation level. Total damages and costs are given by the sum of mitigation cost, adaptation cost and reduced damage with adaptation. Compared to adaptation cost, reduced damage of climate change in SSP2-BaU is more than double adaptation cost in SSP2-BaU. Total damages and costs in SSP2-RCP2.6 are more than in other scenarios because of higher mitigation cost. That is, CO2 emissions reduction with adaptation in SSP2-RCP2.6 creates a great burden to the world. Total damages and costs in SSP2-RCP4.5 are slightly more than that in SSP2-BaU because of higher mitigation cost.
3.3. Aggressive Adaptation for the 1.5˚C Target
Finally, we consider the SSP2-1.5˚C combinations scenario by calculating aggressive adaptation level using revised AD-DICE (ICA-RUS  ) temperature path of SSP2-1.5˚C. The aggressive adaptation level is defined as adaptation level under the CO2 emissions path of SSP2-RCP2.6 which keeps reduced damage of
Figure 3. Regional adaptation costs of the world (percent of GDP) of SSP2-BaU, -RCP4.5 and -RCP2.6 scenarios.
Figure 4. Mitigation and adaptation costs of the world (percent of GDP) of SSP2 scenarios in BaU, RCP4.5 and RCP2.6.
climate damage below the climate damage calculated by the temperature path of SSP2-1.5˚C and the damage function given by Equation (3).
Regional adaptation levels and costs in SSP2-1.5˚C are shown in Figure 5. We find that SSP2-1.5˚C requires more adaptation levels than other SSP2 combination scenarios (Figure 3) even if CO2 emissions in SSP2-RCP2.6 can be attained. The adaptation level of each region peaks around the year 2070 since the damages from global warming in SSP2-1.5˚C increase slower than in SSP2-RCP2.6. Compared to SSP2-BaU, aggressive adaptation level in SSP2-1.5˚C is higher than optimal adaptation level in SSP2-BaU. This means that under the 1.5˚C target each region pays more adaptation costs than in the cases of high temperature rise. Interestingly, regional adaptation level in SSP2-1.5˚C tends to be higher in Japan, USA and Latin America than OAsiaOceania, China and Africa. This tendency is different to other SSP2 scenarios. One reason is that damage from climate change in developed countries increase more rapidly as temperature rises from 1.5˚C to 2.0˚C.
Most regional adaptation costs in SSP2-1.5˚C are higher than in SSP2-BaU, SSP2-RCP4.5 and SSP2-RCP2.6. For example, it amounts to more than one percent of GDP in the 2070s in Japan and FSU. In addition, adaptation cost of OECD8, including EU, amounts to about one percent of its GDP in 2070s. One reason for these high costs is that their economic growth rates are lower than other regions. If a strong target is attained by adaptation, in a region with lower economic growth, adaptation costs are relatively higher.
Figure 5. Adaptation levels and regional adaptation costs (percent of GDP) in the SSP2-1.5˚C scenario.
Figure 6. Damages of climate change, mitigation and adaptation costs of the world (percent of GDP) in the SSP2-1.5˚C scenario.
Damages from climate change, mitigation costs and adaptation costs in SSP2-1.5˚C of the world are shown in Figure 6. The difference between adaptation costs and reduced damages in SSP2-1.5˚C is smaller than in other SSP2 combination scenarios. This indicates that under the 1.5˚C target adaptation costs are high in addition to high mitigation costs for the SSP2-RCP2.6. Therefore, when considering a strong target, it is important to calculate both adaptation costs as well as the costs of damage from global warming.
This research examines regional impacts of adaptation to climate change under the Paris Agreement using a multi-sector model with SSP-RCP combination scenarios.
We firstly find that adaptation costs of Asia and Africa in the 2.0˚C target are higher compared to the other regions in the 2.0˚C target. Thus, some developing countries would have to pay higher adaptation costs than other countries.
The adaptation costs of Japan, EU and Latin America, which are higher than in other countries, amount to one percent of GDP under the SSP2-1.5˚C scenario. This indicates that with the 1.5˚C target, these regions would have higher adaptation costs in comparison to the 2.0˚C target.
To achieve the goal of the Paris Agreement, we need to know not only mitigation costs, which many IAM researchers have simulated, but also adaptation costs. There are large differences among regional adaptation costs which vary by climate change target scenario. Therefore, we should consider both mitigation and adaptation costs for executing policies for challenging climate change.
*This research was supported with funding from the Environment Research and Technology Development Fund (S-10-4) of the Ministry of the Environment, Japan. We would like to thank Yosuke Arino for helpful comments.
 Intergovernmental Panel on Climate Change (IPCC) (2014) Climate Change 2014: Mitigation of Climate Change. Contribution of Working Group III to the Fifth Assessment. Report of the Intergovernmental Panel on Climate Change. Cambridge University Press, Cambridge.
 van Vuuren, D.P., Kriegler, E., O’Neill, B.C., Ebi, K.L., Riahi, K., Carter, T.R., Edmonds, J., Hallegatte, S., Kram, T., Mathur, R. and Winkler, H. (2014) A New Scenario Framework for Climate Change Research: Scenario Matrix Architecture. Climatic Change, 122, 373-386.
 de Bruin, K.C., Dellink, R.B. and Agrawala, S. (2009) Economic Aspects of Adaptation to Climate Change: Integrated Assessment Modelling of Adaptation Costs and Benefits. OECD Environment Working Papers 6.
 Sakaue, S., Yamaura, K. and Washida, T. (2015) Regional and Sectoral Impacts of Climate Change under International Climate Agreements. International Journal of Global Warming, 8, 463-500.
 (2015) The Human Choices to Climate Change Risk of Global Scale. Integrated Climate Assessment-Risks, Uncertainties and Society (ICA-RUS). (In Japanese)