Next Article in Journal
An Object Model for Seafloor Observatory Sensor Control in the East China Sea
Previous Article in Journal
Three-Dimensional Fluid–Structure Interaction Case Study on Elastic Beam
Previous Article in Special Issue
Beach Nourishment as an Adaptation to Future Sandy Beach Loss Owing to Sea-Level Rise in Thailand
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cost-Benefit Analysis of Adaptation to Beach Loss Due to Climate Change in Japan

1
School of Human Science and Environment, University of Hyogo, Himeji 6700092, Japan
2
Faculty of Humanities and Social Sciences, Yamagata University, Yamagata 9908560, Japan
3
International Research Institute of Disaster Science, Tohoku University, Sendai 9808572, Japan
4
Faculty of Urban Science, Meijo University, Nagoya 4618534, Japan
5
Institute for Creative Cities and Regions, University of Hyogo, Kobe 6512197, Japan
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2020, 8(9), 715; https://doi.org/10.3390/jmse8090715
Submission received: 26 August 2020 / Revised: 11 September 2020 / Accepted: 12 September 2020 / Published: 15 September 2020
(This article belongs to the Special Issue Future Coasts and Estuaries)

Abstract

:
To measure economic effects of changes in environmental quality caused by climate change in Japan, we estimate beach loss damage costs in Japan and in each prefecture and evaluate the economic effectiveness of hypothetical adaptation measures to restore sandy beaches. For analyses, we use a computable general equilibrium model (CGE) that integrates a utility function with environmental quality factors as an independent variable derived from a recreation demand function in a travel cost method (TCM). We use future projections of beach loss rates in 2081–2100 based on ensemble-mean regional sea-level rise (SLR) for four Representative Concentration Pathway (RCPs) scenarios (RCP2.6, RCP4.5, RCP6.0, and RCP8.5). The main findings of our study are presented as follows. (1) In 2081–2100, beach loss damage costs were estimated respectively as 398.54 million USD per year for RCP2.6, 468.96 (m.USD/year) for RCP4.5, 494.09 (m.USD/year) for RCP6.0, and 654.63 (m.USD/year) for RCP8.5. (2) For all RCPs, six prefectures for which the cost–benefit ratio exceeds 1.0 were Kanagawa, Osaka, Hyogo, Hiroshima, Saga, and Kumamoto. Our hypothetical adaptation measure of an artificial beach enhancement is expected to be quite effective as a public works project in these prefectures.

1. Introduction

According to the fifth Assessment Report (AR5) published by the IPCC [1], the medium-term and long-term countermeasures are expected to accommodate the possible impacts of climate change. The Ministry of the Environment (MOE) in Japan has discussed planning for a climate change adaptation policy. Some research projects in Japan such as S-8 [2] have forecast climate change effects by region and provide support for adaptative countermeasures. Assessment of climate change effects and the effectiveness of adaptation measures based on nationwide and regional climate change projection are needed.
Numerous attempts have been undertaken to evaluate the economic effects of climate change. Their evaluation methods are classifiable into two approaches: A partial equilibrium approach and a general equilibrium approach. The former method includes a travel cost method (TCM) and a contingent valuation method (CVM). These methods have been applied in some studies to quantify the economic value of the natural environment and ecosystems and the value of statistical life. Since these methods are partial equilibrium approaches, however, economic effects of changes in natural environment by climate change and environmental conservation policies on the whole economy cannot be captured. On the other hand, the latter method has a computable general equilibrium (CGE) analysis. Since a computable general equilibrium model explicitly formulates an objective function in economic agent, direct effects of climate change on economic activities of agent can be captured. In addition, since a CGE model treats all markets in the economy, indirect effects of climate change on the entire economy through changes in the behavior of agents can be captured. Using a CGE model, however, to measure the economic effects of climate change on the natural environment and ecosystems, formulation of the effects on them and estimation of their parameters in a model are necessary. As described above, numerous studies of economic evaluation of climate change have separately been analyzed by two approaches. Therefore, comprehensive assessments in a general equilibrium framework must be made through explicit linkage between a partial equilibrium approach and a general equilibrium approach.
By applying a recreation demand function to a general equilibrium model, for water reallocation issues in Nevada in the United States, Seung et al. [3] analyzed the effects of water reallocation on some recreation sectors and the agriculture sector. However, since the recreation demand function used in their study does not account for the generalized transportation cost, it is not consistent with a utility function. Ciscar et al. [4] comprehensively evaluated the economic and physical effects of climate change on the natural environment, ecosystem, and human society in Europe by treating four sectors as physical effect terms: Agriculture, coastal zone, flood, and tourism. Although their study produced estimates of respective physical effects from projected climate data under conditions of socioeconomic scenarios, and evaluated the projected economic effects by the economic model, it has no theoretical consistency between estimates of physical effect terms and economic models. In a general equilibrium analysis of waste problems in Japan, Miyata [5] derived a utility function consistent with a pre-formulated demand function from solving the integrability problem, and integrated externalities such as waste into a CGE model. For the sandy beach loss caused by climate change, Sakamoto and Nakajima [6] and Nakajima and Sakamoto [7] developed a CGE model that has a utility function consistent with a recreation demand function in a travel cost method by solving the integrability problem. Then, we extend the mode of describing the framework by Nakajima and Sakamoto [7] to simulate more realistic climate change scenarios.
In Japan, although numerous studies such as those of Mimura et al. [8] and Udo et al. [9] have been made of the physical effects of sea level rise and the beach loss caused by climate change, little is known about the economic effects of beach loss and adaptation strategies against climate change. From Table 1, Ohno et al. [10] and Sao et al. [11] evaluated recreation values for a sandy beach in Japan using a travel cost method. The former used beach loss rates calculated by Mimura et al. [8] and estimated the prefectural damage cost of a sandy beach. The latter used beach loss rates by Udo et al. [9], and evaluated not only the prefectural damage cost but also the effect of adaptation measures to restore a sandy beach. Since both studies were based on a partial equilibrium approach, they were unable to treat effects on prices and income caused by changes in environmental quality attributable to climate change. On the other hand, Sakamoto and Nakajima [6], Nakajima and Sakamoto [7], and Sao et al. [11] estimated beach loss damage costs attributable to climate change using a CGE model incorporating a utility function consistent with a recreation demand function in the travel cost method. However, since these studies used future projections by Mimura et al. [8], their economic assessments of beach loss became outdated. Although Sao et al. [11] evaluated adaptation measures for beach loss, their CGE model in these studies had only three goods: Composite goods, gasoline, and an expressway used to visit a sandy beach for recreation. As one might expect, their model framework was quite unrealistic.
Therefore, we sought to measure economic effects of changes in environmental quality attributable to climate change in Japan. Results were obtained using a CGE model that integrates a utility function with environmental quality factors as an independent variable derived from a recreation demand function in a travel cost method (TCM), we aim to estimate the damage cost of beach loss in each prefecture and in Japan and to evaluate the economic effectiveness of hypothetical adaptation measures to restore sandy beaches.

2. Methods and Data

2.1. Structure of Economic Model

We use the 2005 Input–Output table for Japan by MIC [15] as the reference dataset. Table 2 shows 30 sectors that we aggregated in our model. Economic influences comprise household, a production sector, an investment sector, government, and exports and imports.

2.1.1. Household

Figure 1 shows the consumption structure of household in our computable general equilibrium (CGE) model, where index R is used as household consumption for visiting a sandy beach and index H is used as household consumption excluding that for visiting a sandy beach. The set of all affordable bundles that satisfy a consumer’s budget constraint is derived from solving the basic problem of utility maximization. Then, consumption of the petroleum and coal products and transportation for visiting a sandy beach depends on a recreation demand function that incorporates travel cost (the petroleum and coal price, the price of goods and services supplied by the transport sector, and the value of time) and the sandy beach area. For details of derivation of utility function consistent with recreation demand function and definition of goods for visiting a sandy beach, see Appendix A and Appendix B.

2.1.2. Production Sector

Figure 2 shows that all production functions in the domestic production sector are assumed to have a nested function style. For the first step, labor L j and capital K j are aggregated into composite production factor V A j using a Cobb–Douglas production function. As the second step, to produce the gross domestic output Y j for the j-th production sector, the composite production factor is combined with intermediate inputs using a Leontief production function. In addition, the Cobb–Douglas production function allows us to describe substitution between labor L j and capital K j , while the Leontief production function does not between intermediate inputs X i j and composite production factor V A j [16].

2.1.3. Government Sector and Investment Sector

The government sector and investment sector are assumed to have behaviors modeled by Hosoe et al. [16]. The government earns revenues from an income tax, production tax, and indirect tax. Then, the government spends them on purchases of goods proportionately with the constant expenditure share. The structure of investment sector is the same as that of the government sector. In accordance with Hosoe et al. [16], the investment agent collects funds from the household, the government, and the foreign sector. Then, this virtual agent purchases investment goods proportionately with a constant share.

2.1.4. Export and Import

In accordance with Hosoe et al. [16], Figure 3 portrays the structure of the substitution between imports and domestic goods and that of the transformation between exports and domestic goods. Regarding imperfect substitution between imports and domestic goods, we adopt Armington’s assumption [17]. The i-th Armington-composite-good-producing sector aggregates domestic goods D i and imports I M i into composite goods Q i using a constant elasticity of substitution (CES) function. However, gross domestic output Y i is transformed into domestic goods D i and exports E X i using a constant elasticity of transformation (CET) function. Both parameters of elasticity of transformation σ D E X and elasticity of substitution σ D I M are assumed to be 2.0 exogenously.

2.2. Setting of Scenarios

2.2.1. Scenario of Beach Loss Caused by Climate Change

We use the results of the future projections of beach loss in Japan by Udo and Takeda [14] as the Beach loss scenario. Udo and Takeda [14] calculated the beach-loss rates in 2081–2100 relative to the reference period 1986–2005 based on ensemble-mean regional sea-level rise (SLR) for four Representative Concentration Pathway (RCPs) scenarios (RCP2.6, RCP4.5, RCP6.0, and RCP8.5) using 21 models of the Coupled Model Intercomparison Project Phase 5 (CMIP5). We specifically examine four scenarios in 2081–2100 based on ensemble mean SLR, and use 21 scenarios for uncertainty assessment using 21 models of CMIP5 (ACCESS 1.0, BCC CSM 1.1, CanESM2, CNRM CM5, CSIRO3.6.0, NOAA GFDL-ESM2M, NOAA GFDL-ESM2G, GISS-E2-R, HadGEM2-CC, HadGEM2-ES, INM-CM4, IPSL-CM5A-LR, IPSL-CM5A-MR, MIROC-ESM-CHEM, MIROC-ESM, MIROC5, MPI-ESM-LR, MPI-ESM-MR, MRI_CGCM3, NorESM-ME, NorESM1-M). In addition, Table 3 presents the minimum and maximum beach loss rates in each prefecture in our beach loss scenarios. Table 4 and Figure 4 show 47 prefectures in Japan.

2.2.2. Scenario of Adaptation Measure for Restoring Sandy Beach

For hypothetical adaptation measures related to beach loss, we assume that after erosion of coastal areas caused by the sea-level rise, the coastal area can be restored to its earlier state by implementation of adaptation measures such as a public works project using artificial beach enhancement. From considerations of data availability and comparison with earlier studies, we chose to use an average adaptation cost per unit area assumed by Sao et al. [12] for the scenario of adaptation measures for restoring sandy beaches. Sao et al. [12] collected data including those of 92 public works in 33 prefectures related to artificial beach enhancement, and assumed the average adaptation cost per unit area as 215.96 USD/m2 from available data for sandy beaches. However, since few public works projects are limited to artificial beach enhancement and since these projects include costs of protecting land unrelated to sandy beaches, it is noteworthy that the average adaptation cost that we assumed might be overestimated.
Finally, to estimate beach loss damage costs and to evaluate the economic effectiveness of hypothetical adaptation measures to restore sandy beaches, we measure the benefit as equivalent variation. For details of the definition of benefit, see Appendix C.

3. Results and Discussion

3.1. Economic Effects of Beach Loss

Figure 5 shows damage costs presented by the projected beach loss for RCP2.6, RCP4.5, RCP6.0, and RCP8.5 in 2081–2100. These figures indicate that the higher the future temperature becomes, the more the damage cost of sandy beach increases. In 2081–2100, damage costs are estimated respectively as 398.54 million USD per year for RCP2.6, 468.96 (m.USD/year) for RCP4.5, 494.09 (m.USD/year) for RCP6.0, and 654.63 (m.USD/year) for RCP8.5.
Figure 5 and Table 5 present prefectural damage costs because of the projected beach loss in four RCPs in 2081–2100. For any RCPs, damage costs of four prefectures (Okinawa, Kanagawa, Niigata and Hyogo) accounts for about 40% to about 45% of the total damage cost to Japan. As shown in Table 5, for RCP2.6, damage costs of these four prefectures were estimated respectively as 22.85 million USD per year to 87.42 (m.USD/year). For RCP8.5, four prefectural damage costs are estimated as 33.87 (m.USD/year) to 87.42 (m.USD/year).
Figure 6 and Table 6 show prefectural damage costs per unit area attributable to beach loss in four RCPs in 2081–2100. Prefectures for which the damage cost per unit area is high were Kanagawa, Niigata, Toyama, Fukui, Kyoto, Osaka, Hyogo, Wakayama, Okayama, Hiroshima, Saga, Kumamoto, and Okinawa. Especially, damage costs per unit area in prefectures in western Japan or along the Inland Sea tend to be higher. For RCP2.6, damage costs per unit area to the tenth highest prefecture were Saga, Kumamoto, Kanagawa, Osaka, Hiroshima, Hyogo, Okayama, Okinawa, Fukui, and Toyama in order from the highest, estimated as 138.47 USD per unit area to 723.26 (USD/m2). For RCP8.5, damage costs per unit area to the tenth prefecture were Kanagawa, Saga, Kumamoto, Osaka, Okayama, Hyogo, Toyama, Hiroshima, Kyoto, and Wakayama, estimated as 178.34 to 765.76 (USD/m2).
For uncertainty assessment, we calculated 21 beach loss scenarios in 2081–2100 using 21 CMIP5 models. For RCP4.5, damage costs were estimated respectively as an average of 491.06 million USD/year, a minimum of 385.72 (m.USD/year), and a maximum of 739.50 (m.USD/year). In addition, Figure 7 shows prefectural damage costs per unit area using 21 CMIP5 models. As shown in Figure 7, although results of Kanagawa and Toyama have a large variance, those of many other prefectures have a small variance.

3.2. Cost–Benefit Analysis of Adaptation Policy

Figure 8 and Table 6 show cost–benefit ratios of adaptation policies for RCP2.6, for RCP4.5, for RCP6.0, and for RCP8.5 in 2081–2100. Prefectures in red in Figure 8 and shaded values in Table 6 have cost–benefit ratios larger than 1.0, i.e., the benefit from adaptation measures exceeds the cost because of beach loss. As described above, we assumed 215.96 USD/m2 of the average cost per unit area as the adaptation cost to restore a sandy beach. The higher the future temperature becomes, the more numerous the prefectures for which adaptation measures are cost-effective become. Especially in four prefectures along the Inland Sea, which are Osaka, Hyogo, Okayama, and Hiroshima, our hypothetical adaptation measure as a public works project of artificial beach enhancement is quite effective.
For RCP2.6, six prefectures for which the cost–benefit ratio exceeds 1.0 were Kanagawa, Osaka, Hyogo, Hiroshima, Saga, and Kumamoto. In contrast, for RCP8.5, eight prefectures with a cost–benefit ratio over 1.0 were Kanagawa, Toyama, Osaka, Hyogo, Okayama, Hiroshima, Saga, and Kumamoto.

3.3. Discussion

We compare the results of beach loss damage costs attributable to climate change with those of earlier studies. As shown in Table 1, Ohno et al. [10], Sakamoto and Nakajima [6], Nakajima and Sakamoto [7], and Sao et al. [11] used the future projection of beach loss calculated by Mimura et al. [8] and estimated the damage costs of sandy beach because of the sea level rise from 30 to 100 cm. These earlier studies estimated damage costs of the sea level rise as 247 to 832 (m.USD/year). It is apparent that the results of our study are slightly lower than those of earlier studies. Especially, although differences between results of Ohno et al. [10] and our study are larger, it is likely that these results became overestimated since Ohno et al. [10] formulated damage costs of beach loss by a proportional relation between the frequency of visiting the sandy beach for recreation and the sandy beach area.
Sao et al. [12] and Nakajima et al. [13] used the future projection of beach loss by Udo et al. [9] and respectively estimated damage costs for RCP2.6 and RCP8.5. Sao et al. [12] estimated them as 254–284 (m.USD/year) in 2031–2050 and 426–494 (m.USD/year) in 2081–2100. One reason for the difference between results reported by Sao et al. [12] and those of our study is that our general equilibrium approach reflects price changes and income changes that are not considered in the definition of consumer surplus derived from the partial equilibrium approach. Consequently, it is apparent that beach loss damage costs in our study are slightly lower than those found in earlier studies.
Table 7 shows the number of prefectures for which the cost–benefit ratio exceeds 1.0 in the adaptation scenarios using 215.96 and 182.76 USD/m2 as the average cost per unit area. Although the number of cost-effective prefectures between Sao et al. [11] and Sao et al. [12] is significantly different, Sao et al. [12] described that the difference between these studies resulted from the average adaptation cost per unit area. As described above, for the possibility that the average adaptation cost per unit area (215.96 USD/m2) that Sao et al. [12] assumed could be overestimated, we compared the effects of two adaptation costs. In both adaptation scenarios, the number of prefectures for which adaptation measures were cost–effective in our results was almost identical for all RCPs. Consequently, it is apparent that the results of our study are more robust than those of earlier studies. In addition, six prefectures for which the cost-benefit ratio exceeds 1.0 for all RCPs have a large damage cost despite the small area of their beaches. In other words, we consider the higher damage cost per unit area to be the reason why hypothetical adaptation measures are economically efficient.
Figure 9 portrays effects of two adaptation measures of RCP8.5 in 2081–2100. As described above, we demonstrated that the higher the future temperature becomes, the greater the number of prefectures for which adaptation measures are cost–effective. From Figure 9, the lower adaptation cost makes adaptation measures in Kyoto more effective. Additionally, one assumes that the adaptation cost becomes much lower, then we can say that adaptation measures in Niigata, Wakayama, Fukui, and Okinawa are potentially cost-effective. Therefore, it is apparent that the lower the average adaptation cost per unit becomes, the more numerous prefectures for which the adaptation measures are cost–effective become. Especially, in some prefectures along the Inland Sea such as Osaka, Hyogo, Okayama, and Hiroshima, our hypothetical adaptation measure of a public works project of an artificial beach enhancement is quite effective.

4. Conclusions

To assess the economic effects of changes in environmental quality caused by climate change in Japan, we used a computable general equilibrium model that integrates a utility function with environmental quality factors as independent variables derived from a recreation demand function in a travel cost method. Results show the estimated damage costs of beach loss in Japan and in the respective prefectures. We evaluated the economic effectiveness of hypothetical adaptation measures to restore sandy beaches. The findings obtained from this study are presented below.
  • Higher future temperatures will cause higher damage costs of sandy beaches. In 2081–2100, we estimated damage costs as 398.54 million USD per year for RCP2.6, 468.96 (m.USD/year) for RCP4.5, 494.09 (m.USD/year) for RCP6.0, and 654.63 (m.USD/year) for RCP8.5, respectively.
  • For all RCPs, six prefectures for which the cost–benefit ratio exceeds 1.0 were Kanagawa, Osaka, Hyogo, Hiroshima, Saga, and Kumamoto.
  • Higher future temperatures will bring high numbers of prefectures for which adaptation measures are cost–effective. Especially for four prefectures along the Inland Sea, which are Osaka, Hyogo, Okayama, and Hiroshima, our hypothetical adaptation measure of an artificial beach enhancement is expected to be quite effective as a public works project.
Further examinations can be expected to support further discussion. First, since we were unable to treat the adaptation cost endogenously, we will incorporate endogenous adaptation costs into our CGE model and evaluate the effectiveness of some adaptation measures. Secondly, since we evaluated only the recreation value (use value) estimated using TCM, we expect to develop a CGE model incorporating evaluation methods of non-use values.

Author Contributions

Project administration, K.N.; supervision, K.N. and N.S.; methodology, K.N., N.S., E.O. and M.M.; formal analysis, K.N.; resources, K.U. and Y.T.; data curation, K.N. and R.M.; writing—original draft preparation, K.N.; writing—review and editing, K.N. and N.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the Social Implementation Program on Climate Change Adaptation Technology (SI-CAT), the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan, and by the JSPS KAKENHI Grant-in-Aid for Scientific Research (C), grant number JP19K12448.

Acknowledgments

The authors would like to thank the SI-CAT and JSPS KAKENHI for their generous financial support.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Derivation of Utility Function Consistent with Recreation Demand Function

In accordance with Nakajima and Sakamoto [7], we present the expression for a utility function from a recreation demand function estimated using the zone travel cost method, and apply it to a computable general equilibrium model.
First, we assume that some regions are divided into N zones. Each zone has a natural environment, which is a sandy beach in our study. In addition, a recreation demand function from zone a to zone b is estimated using regression analysis as shown below. For elimination of negative estimation of visiting rate and treatment of heteroskedasticity attributable to different populations among zones, we employ a semi-logarithmic function. For details, see Cooper and Loomis [18].
ln ( z a b n a ) = γ ^ 0 + γ ^ 1 p a b + γ ^ 2 q b + ε ^ a b ,   a , b
p a b p 2 Y g a b + p 3 Y c a b + w t a b p N H ,   a , b
In Equations (A1) and (A2), z a b denotes the frequency of visits from zone a to zone b ; n a represents the population in zone a . ( z a b / n a ) signifies the visitation rate from zone a to zone b . Therein, p a b represents the travel cost per visit necessary to make a round trip from zone a to zone b . It is defined as a relative price of composite goods price p N H in Equation (A2). For it, p a b = p b a . Then, p 2 Y , p 3 Y , and w respectively denote the price of the petroleum and coal products sectors, the price of goods and services supplied by the transport sector, and the value of time as described below. Furthermore, g a b , c a b , and t a b , respectively represent the amount of gasoline per visit, expressway use per visit, and time per visit. Furthermore, q b signifies the natural environmental quality in zone b , which is a sandy beach in this study. By formulating the natural environmental quality as an explanatory variable explicitly, we can make a computable general equilibrium analyses by changes in environmental quality. In Equation (A2), γ ^ 0 , γ ^ 1 , and γ ^ 2 are estimated parameters. They are γ ^ 1 < 0 and γ ^ 1 > 0 ; ε ^ a b is a residual.
Since a demand function is derived from solving a utility maximization problem of household, a utility function exists corresponding to its demand function. This is known as an integrability problem. For details related to the integrability problem, one can consult works by Varian [19], Mas-Colell et al. [20], and Jehle and Reny [21]. By solving the integrability problem, a utility function with Equation (A1) as a recreation demand function and a budget constraint are derived as shown below.
u a = x a + 1 γ ^ 1 b z a b ( ln z a b ln Γ a b 1 ) ,   a
p N H x a + a ( p 2 Y g a b + p 3 Y c a b ) z a b = M a ,   a
Γ a b n a exp { γ ^ 0 + γ ^ 1 ( w t a b p N H ) + γ ^ 2 q b + ε ^ a b } ,   a , b
In those equations, u a denotes household utility in zone a , x a represents consumption of composite goods in zone a , and M a is the household income in zone a . Additionally, Γ a b signifies the number of visitors from zone a to zone b when the price of gasoline and toll fees between zone a and zone b equal zero ( p 2 Y = p 3 Y = 0 ). If γ ^ 1 < 0 , then it is z a b Γ a b in subjective equilibrium of household. The second term of the right-hand side in Equation (A3) is non-negative. For simplification, since our model relies on the assumption that there exists only a single household, a utility function and budget constraint are formulated by the summation of each variable with respect to the zone in Equations (A3) and (A4), as:
U = N H + 1 γ ^ 1 a b z a b ( ln z a b ln Γ a b 1 )
p N H N H + a b ( p 2 Y g a b + p 3 Y c a b ) z a b = M
where U represents utility ( = a u a ) , N H denotes consumption of composite goods ( = a x a ) , and M signifies income ( = a M a ) in household.
Parameters of a recreation demand function are estimated using the function form shown in Equation (A1). Ohno et al. [10] created a dataset with travel cost and the amount of traffic for visiting a sandy beach from MAFF [22] database and estimated parameters of a recreation demand function. Furthermore, environmental quality data are created as a logarithmic value of the sandy beach area by prefecture estimated by Mimura et al. [8]. Then, we employ the same data set as Ohno et al. [10] and the environmental quality data and estimate the parameters of our recreation demand function using the least-square method. Table A1 presents the estimated parameters. Since regression coefficients of travel cost are estimated as negative and since the coefficient of environmental quality is positive, the sign condition is satisfied. In addition, all variables were found to be significant. It is apparent that the coefficient of environmental quality is less than 1, and that the frequency of visits decreases gradually as the sandy beach area decreases. However, Ohno et al. [10] assumed that changing the sandy beach area changes the frequency of visits proportionally.
Table A1. Estimated parameters in the utility function.
Table A1. Estimated parameters in the utility function.
CoefficientEstimated Valuet-Value
γ 0 4.604 −7.575
γ 1 4.110 × 10 4 −14.029
γ 2 0.329 3.178
R ¯ 2 0.477
No. of observations227

Appendix B. Definition of Goods for Visiting a Sandy Beach

Let x 2 , a b represent gasoline consumption and x 3 , a b be the use of expressway needed to travel between zone a and zone b . Furthermore, the frequency of visiting them is denoted by z a b in independent variables of the utility function. The relations among them can be expressed as shown below.
x 2 , a b = g a b · z a b ,   a , b
x 3 , a b = c a b · z a b ,   a , b
Equations (A8) and (A9) can be shown as optimal solutions in a cost minimization problem based on a production function with the Leontief technology as presented below.
z a b = min { x 2 , a b g a b , x 3 , a b c a b } ,   a , b
From Equation (A10), it is apparent that a household produces a visit for a recreation site, and that Equation (A10) is a part of the utility function. It is presented as U Z a b in Figure 1.
Secondly, gasoline consumption and expressway use by a household are aggregated with respect to all zones, respectively, as Equations (A11) and (A12).
X 2 H = a b x 2 , a b
X 3 H = a b x 3 , a b
According to MIC [19], gasoline is produced in the “petroleum refinery” sector. The use of an expressway is produced in the “travel agency and other services related to the transport” sector, in the input–output table for Japan, which comprises 190 sectors. We assume that hypothetical sectors produce gasoline and expressways to visit a sandy beach.
Since we estimate parameters in a recreation demand function using travel cost data for Japan and annual traffic data among all prefectures, we can measure the annual gasoline consumption and a use of expressway to visit a sandy beach in a money metric by using these data. The annual gasoline consumption for visiting a sandy beach is divided from household consumption of the petroleum refinery products in the input–output table. Additionally, intermediate inputs and factor inputs in the petroleum refinery sector are divided at the same rate as household consumption. Similarly, we treat the annual use of an expressway for visiting a sandy beach. Other production sectors are aggregated as the composite goods sector. In summary, our model has three goods and sectors that include gasoline consumption for visiting a sandy beach, use of an expressway for visiting a sandy beach, and composite goods.

Appendix C. Definition of Benefit

Travel cost methods measure consumer surplus C S from Equation (A1) to evaluate natural environments economically. Actually, C S is proportional to the sum of the frequency of visits as shown below. In addition, p Y = ( p 1 Y , p 2 Y , p 3 Y ) is a price vector; q = ( q 1 , , q N ) is an environmental quality vector.
C S ( p Y , q ) a b p 2 Y g a b + p 3 Y c a c + w t a b n a exp { γ ^ 0 + γ ^ 1 ( l p N H ) + γ ^ 2 q b + ε ^ a b } d l = p N H γ ^ 1 a b z a b
In general, when environmental quality changes in q 0 q 1 in a travel cost method, from Δ C S = C S ( p Y 0 , q 1 ) C S ( p Y 0 , q 0 ) , the change in an environmental quality is evaluated by price p Y 0 fixed in the initial period. Although this evaluation method is a simple process, without fixing a price vector, the benefit by change in environmental quality should be measured by considering changes in the price vector. We aim at measuring benefits by considering changes in price and income.
We assume that we define the benefit as equivalent variation (EV). Since the indirect utility function converts the sum of the factor income and the consumer surplus into a composite goods term, it can be expressed as shown below.
V ( p Y , M , q ) = M + C S ( p Y , q ) p N H
Since the expenditure function is an inverse function of the indirect utility function with respect to income, it can be derived from solving for M , where V represents utility.
E ( p Y , V , q ) = p N H V C S ( p Y , q )
We assume that the beach loss causes changes in the sandy beach area of q 0 to q 1 . From the perspective of general equilibrium analysis, such a change in the exogenous variable affects prices and factor incomes. With this beach loss, it is assumed that the price system changes from p Y 0 to p Y 1 , and that the factor income changes from M 0 to M 1 . Then, the equivalent variation is shown as presented below.
E V = E ( p Y 0 , V ( p Y 1 , M 1 , q 1 ) , q 0 ) E ( p Y 0 , V ( p Y 0 , M 0 , q 0 ) , q 0 ) E V = p N H 0 [ M 1 + C S ( p Y 1 , q 1 ) p N H 1 M 0 + C S ( p Y 0 , q 0 ) p N H 0 ]
Next, to define a prefectural equivalent variation, we rewrite Equation (A13) by the summation with respect to zone a in Equation (A13) as follows.
C S b ( p Y , q ) a p 2 Y g a b + p 3 Y c a c + w t a b n a exp { γ ^ 0 + γ ^ 1 ( l p N H ) + γ ^ 2 q b + ε ^ a b } d l = p N H γ ^ 1 a z a b
We rewrite Equation (A16) using Equation (A17). We define the prefectural equivalent variation in zone b as shown below.
E V b = E b ( p Y 0 , V ( p Y 1 , M 1 , q 1 ) , q 0 ) E b ( p Y 0 , V ( p Y 0 , M 0 , q 0 ) , q 0 ) E V b = p N H 0 [ ( M 1 p N H 1 M 0 p N H 0 ) 1 γ ^ 1 { a ( z a b 1 z a b 0 ) } ]

References

  1. IPCC. Climate Change 2014 Mitigation of Climate Change: Working Group III Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar]
  2. S-8 Climate Change Impact and Adaptation Research Project Team. Climate Change “Impacts on Japan”: Comprehensive Impact Assessment and Adaptation Measures Based on New Scenarios; The Environment Research and Technology Development Fund (ERTDF), Ministry of the Environment, Japan’s Strategic Research and Development Domain S-8 Comprehensive Study on Impact Assessment and Adaptation for Climate Change 2014 Report; National Institute for Environmental Studies: Tokyo, Japan, 2014. [Google Scholar]
  3. Seung, C.K.; Harris, T.R.; Englin, J.E.; Netusil, N.R. Inpacts of water reallocation: A combined computable general equilibrium and recreation demand model approach. Ann. Reg. Sci. 2000, 34, 473–487. [Google Scholar] [CrossRef]
  4. Ciscar, J.-C.; Iglesias, A.; Feyen, L.; Szabo, L.; Regemorter, D.V.; Amelung, B.; Nicholls, R.; Watkiss, P.; Christensen, O.B.; Dankers, R.; et al. Physical and economic consequences of climate change in Europe. Proc. Natl. Acad. Sci. USA 2011, 108, 2678–2683. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Miyata, Y. A General equilibrium analysis of waste-economic system: A CGE modeling approach. J. JSCE D 1995, 12, 259–270. [Google Scholar] [CrossRef] [Green Version]
  6. Sakamoto, N.; Nakajima, K. Development of general equilibrium model consistent with travel cost method. J. JSCE G 2012, 68, II_217–II_228. (In Japanese) [Google Scholar] [CrossRef] [Green Version]
  7. Nakajima, K.; Sakamoto, N. General equilibrium approach consistent with travel cost method for economic evaluation of beach erosion by climate change. In Proceedings of the 53rd of European Congress of the Regional Science Association International, Palermo, Italy, 27–31 August 2013. [Google Scholar]
  8. Mimura, N.; Inoue, K.; Ikusebashi, S.; Izumiya, T.; Nobuoka, N. Title in Japanese only (Sunahama-ni-taisuru Kaimen-jousyou-no Eikyou-hyouka(2): Yosoku-model-no Datousei-no Kensyou-to Zenkoku-kibo-no Hyouka). J. Coast. Eng. 1994, 41, 1161–1165. (In Japanese) [Google Scholar] [CrossRef]
  9. Udo, K.; Takeda, Y.; Yoshida, J.; Mano, A. Future projections of beach erosion in Japan using sea level change data of the MIROC5 model. J. JSCE G 2013, 69, I_239–I_247. (In Japanese) [Google Scholar] [CrossRef] [Green Version]
  10. Ohno, E.; Hayashiyama, Y.; Morisugi, H.; Nohara, K. Global warming damage cost of sandy beach loss: Travel cost method approach. Glob. Environ. Res. 2009, 14, 291–297. (In Japanese) [Google Scholar]
  11. Sao, H.; Morisugi, M.; Ohno, E.; Sakamoto, N.; Nakajima, K.; Morisugi, M. Examination of regional impacts of adaptation policies and sand erosion damage due to climate change. J. JSCE G 2013, 69, I_249–I_257. (In Japanese) [Google Scholar] [CrossRef]
  12. Sao, H.; Tomoda, G.; Mori, R.; Morisugi, M.; Ohno, E.; Nakajima, K.; Sakamoto, N. Estimation of loss of recreation value due to sand beach erosion and effect of adaptation policy. J. JSCE G 2017, 73, I_191–I_199. (In Japanese) [Google Scholar] [CrossRef]
  13. Nakajima, K.; Sakamoto, N.; Ohno, E.; Morisugi, M.; Mori, R. Cost benefit analysis of adaptation against beach erosion due to climate change. J. JSCE G 2018, 74, I_425–I_436. [Google Scholar] [CrossRef]
  14. Udo, K.; Takeda, Y. Projections of future beach loss in japan due to sea-level rise and uncertainties in projected beach loss. Coast. Eng. J. 2017, 59, 1740006. [Google Scholar] [CrossRef] [Green Version]
  15. Ministry of Internal Affairs and Communications (MIC). 2005 Input–Output Tables for Japan. 2009. Available online: https://warp.da.ndl.go.jp/info:ndljp/pid/8943652/www.stat.go.jp/english/data/io/2005/zuhyou/ioe05103.xls (accessed on 10 August 2020).
  16. Hosoe, N.; Gasawa, K.; Hashimoto, H. Textbook of Computable General Equilibrium Modelling: Programming and Simulations; Palgrave Macmillan: Hampshire, UK, 2010. [Google Scholar]
  17. Armington, P.S. A theory of demand for products distinguished by place of production. Int. Monet. Fund Staff Papers 1969, 16, 159–176. [Google Scholar] [CrossRef]
  18. Cooper, J.; Loomis, J. Testing whether waterfowl hunting benefits increase with greater water deliveries to wetlands. Environ. Resour. Econ. 1993, 13, 545–561. [Google Scholar] [CrossRef]
  19. Varian, H.R. Microeconomic Analysis, 3rd ed.; W. W. Norton & Company, Ltd.: London, UK, 1992. [Google Scholar]
  20. Mas-Colell, A.; Whinston, M.D.; Green, J.R. Microeconomic Theory; Oxford University Press: New York, NY, USA, 1995. [Google Scholar]
  21. Jehle, G.A.; Reny, P.J. Advanced Microeconomic Theory, 2nd ed.; Addison Wesley: New York, NY, USA, 2000. [Google Scholar]
  22. Ministry of Agriculture, Forestry and Fisheries (MAFF). The 2003 (11th) Census of Fisheries of Japan. 2004. Available online: https://www.e-stat.go.jp/en/stat-search/file-download?statInfId=000023622426&fileKind=0 (accessed on 10 August 2020).
Figure 1. Structure of household.
Figure 1. Structure of household.
Jmse 08 00715 g001
Figure 2. Structure of the production sector.
Figure 2. Structure of the production sector.
Jmse 08 00715 g002
Figure 3. Structure of exports and imports.
Figure 3. Structure of exports and imports.
Jmse 08 00715 g003
Figure 4. 47 prefectures in Japan.
Figure 4. 47 prefectures in Japan.
Jmse 08 00715 g004
Figure 5. Prefectural damage costs attributable to beach loss in four representative concentration pathway (RCPs) scenarios in 2081–2100 (ad).
Figure 5. Prefectural damage costs attributable to beach loss in four representative concentration pathway (RCPs) scenarios in 2081–2100 (ad).
Jmse 08 00715 g005
Figure 6. Prefectural damage costs per unit area for RCP2.6 (a) and RCP8.5 (b) in 2081–2100.
Figure 6. Prefectural damage costs per unit area for RCP2.6 (a) and RCP8.5 (b) in 2081–2100.
Jmse 08 00715 g006
Figure 7. Prefectural damage cost per unit area using the beach loss rates calculated by 21 CMIP5 models.
Figure 7. Prefectural damage cost per unit area using the beach loss rates calculated by 21 CMIP5 models.
Jmse 08 00715 g007
Figure 8. Cost–benefit ratios of adaptation policies in four RCPs in 2081–2100.
Figure 8. Cost–benefit ratios of adaptation policies in four RCPs in 2081–2100.
Jmse 08 00715 g008
Figure 9. Comparison of cost–benefit ratios of two adaptation scenarios for RCP8.5 in 2081–2100.
Figure 9. Comparison of cost–benefit ratios of two adaptation scenarios for RCP8.5 in 2081–2100.
Jmse 08 00715 g009
Table 1. Comparison with damage costs reported from earlier studies.
Table 1. Comparison with damage costs reported from earlier studies.
Study [10][6,7][11][12][13]Our study
Method TCMCGE + TCMCGE + TCMTCMCGE + TCMCGE + TCM
Future projection [8][8][8][9][9][14]
Climate model ------MIROC53 models21 models
Sea Level RiseSLR: 30 cm (A)522247290------
(B)56.6%56.6%26.5–97.0%------
SLR: 65 cm (A)753440530------
(B)81.7%81.7%53.9–100.0%------
SLR: 100 cm (A)832551--------
(B)90.3%90.3%--------
(A) Damage cost/(B) Rate of coastal erosionRCP2.62031–2050(A)------254116–147–184--
(B)------11.9–74.6%11.9–74.6%--
2081–2100(A)------426335–402–440399
(B)------25.6–100.0%25.6–100.0%22.6–100.0%
RCP4.52031–2050(A)--------142–150–197--
(B)--------11.8–69.2%--
2081–2100(A)--------410–462–471469
(B)--------28.0–100.0%26.8–100.0%
RCP6.02031–2050(A)------------
(B)------------
2081–2100(A)----------494
(B)----------27.2–100.0%
RCP8.52031–2050(A)------284174–179–252--
(B)------13.9–83.8%13.9–83.8%--
2081–2100(A)------494615–644–644654
(B)------36.1–100.0%36.1–100.0%37.0–100.0%
Unit: 1 million USD/year.
Table 2. Sector classification.
Table 2. Sector classification.
SectorCode SectorCode
1AgricultureAGR16Information and communication equipmentICE
2ForestryFRS17electronic componentsELC
3FisheryFSH18Transportation equipmentTRE
4MiningMIN19Precision instrumentsPRI
5FoodsFOD20Other manufacturing productsOMF
6Textile and paper productsTEX21ConstructionCNS
7Chemical productsCPR22ElectricityELY
8Petroleum refinery productsPET23CommerceCOM
9Coal productsCOL24Finance and insuranceFIN
10Ceramic, stone and clay productsCER25Real estateEST
11Iron and steelIRN26TransportTRP
12Non-ferrous metalsNFM27Facility service for road transportRTP
13Metal productsMTL28Public administrationPBA
14General machineryMCH29AccommodationsACM
15Electrical machineryELM30Other servicesSRV
Table 3. Future projections of beach loss rate in Japan by Udo and Takeda (2017).
Table 3. Future projections of beach loss rate in Japan by Udo and Takeda (2017).
RCP2.6RCP4.5RCP6.0RCP8.521 Models
Mean66.975.577.588.576.5
[min., max.][22.6, 100.0][26.8, 100.0][27.2, 100.0][37.0, 100.0][28.8, 100.0]
Table 4. 47 prefectures in Japan.
Table 4. 47 prefectures in Japan.
PrefectureCode PrefectureCode PrefectureCode
1HokkaidoHKD17IshikawaISK33OkayamaOKY
2AomoriAMR18FukuiFKI34HiroshimaHRS
3IwateIWT19YamanashiYMN35YamaguchiYGC
4MiyagiMYG20NaganoNGN36TokushimaTKS
5AkitaAKT21GifuGIF37KagawaKGW
6YamagataYGT22ShizuokaSZK38EhimeEHM
7FukushimaFKS23AichiACH39KochiKOC
8IbarakiIBR24MieMIE40FukuokaFKO
9TochigiTCG25ShigaSIG41SagaSAG
10GunmaGNM26KyotoKYT42NagasakiNGS
11SaitamaSTM27OsakaOSK43KumamotoKMT
12ChibaCHB28HyogoHYG44OitaOIT
13TokyoTKY29NaraNAR45MiyazakiMYZ
14KanagawaKNG30WakayamaWKY46KagoshimaKGS
15NiigataNGT31TottoriTTR47OkinawaOKW
16ToyamaTYM32ShimaneSMN
Table 5. Prefectural damage costs attributable to beach loss in four RCPs in 2081–2100.
Table 5. Prefectural damage costs attributable to beach loss in four RCPs in 2081–2100.
Pref.RCP2.6RCP4.5RCP6.0RCP8.5 Pref.RCP2.6RCP4.5RCP6.0RCP8.5
1HKD9.7412.0512.5317.1325SIG0.000.000.000.00
2AMR1.441.751.782.5226KYT5.507.037.7312.26
3IWT2.222.862.914.1327OSK2.883.523.675.59
4MYG2.663.033.125.0128HYG22.8525.9327.3333.87
5AKT2.363.013.135.1129NAR0.000.000.000.00
6YGT3.123.894.076.3130WKY12.9417.9818.6126.42
7FKS2.352.913.095.5631TTR0.490.650.650.95
8IBR6.839.2610.3213.9732SMN4.195.355.599.26
9TCG0.000.000.000.0033OKY4.245.655.7212.23
10GNM0.000.000.000.0034HRS14.5614.5614.5614.56
11STM0.000.000.000.0035YGC15.9915.9915.9915.99
12CHB11.2715.4117.3724.7036TKS1.181.501.612.44
13TKY0.761.141.181.4137KGW3.804.565.049.37
14KNG39.2953.5062.1187.3238EHM7.8011.3012.3715.24
15NGT26.8431.0131.8251.4539KOC1.682.452.705.55
16TYM6.649.5510.1218.9240FKO3.775.556.3110.63
17ISK1.962.602.804.1341SAG7.177.177.177.17
18FKI9.8710.9811.3713.5542NGS15.0116.2016.2016.20
19YMN0.000.000.000.0043KMT15.4715.4715.4715.47
20NGN0.000.000.000.0044OIT2.893.593.846.94
21GIF0.000.000.000.0045MYZ6.139.2111.2917.47
22SZK10.4612.3312.9217.8246KGS11.5115.3715.7718.42
23ACH4.615.726.048.1747OKW87.4287.4287.4287.42
24MIE8.6411.5212.3923.97 Total398.54468.96494.09654.63
Unit: 1 million USD/year.
Table 6. Prefectural damage costs per unit area in four RCPs in 2081–2100.
Table 6. Prefectural damage costs per unit area in four RCPs in 2081–2100.
Pref.RCP2.6RCP4.5RCP6.0RCP8.5 Pref.RCP2.6RCP4.5RCP6.0RCP8.5
1HKD6.557.107.228.5125SIG0.000.000.000.00
2AMR6.446.556.576.8626KYT136.12145.50150.02183.80
3IWT26.3929.8830.1438.2727OSK409.10427.33431.78493.81
4MYG26.4227.2827.4932.5328HYG232.02248.90256.99298.80
5AKT15.2216.2316.4320.1129NAR0.000.000.000.00
6YGT68.9172.4373.3385.1230WKY101.59126.07129.48178.34
7FKS30.8332.9433.6345.5031TTR6.126.316.316.68
8IBR62.2165.0166.2970.9332SMN41.0644.4745.2358.68
9TCG0.000.000.000.0033OKY174.81198.67200.03360.31
10GNM0.000.000.000.0034HRS239.17239.17239.17239.17
11STM0.000.000.000.0035YGC92.5592.5592.5592.55
12CHB77.4082.5085.0795.6136TKS12.7714.0814.5618.88
13TKY4.875.555.646.1137KGW50.8554.1156.3181.33
14KNG487.64556.09602.93765.7638EHM54.7271.3877.3294.63
15NGT102.64107.68108.68137.0139KOC8.419.8310.3618.45
16TYM138.47164.64170.42285.8840FKO31.1936.8039.5359.50
17ISK14.8015.8216.1518.5941SAG723.26723.26723.26723.26
18FKI139.76146.55149.00163.7142NGS78.7584.9684.9684.96
19YMN0.000.000.000.0043KMT560.73560.73560.73560.73
20NGN0.000.000.000.0044OIT39.6743.5045.0068.39
21GIF0.000.000.000.0045MYZ29.8736.4441.7761.63
22SZK41.9343.4043.8848.1246KGS15.4018.7419.1321.86
23ACH40.0042.5143.2648.7347OKW148.66148.66148.66148.66
24MIE45.1451.7653.9792.82 Total57.5959.6361.2470.62
Unit: USD/m2.
Table 7. Comparison with numbers of cost-effective prefectures of earlier studies.
Table 7. Comparison with numbers of cost-effective prefectures of earlier studies.
Study[11][12] [13] Our study
MethodCGE + TCMTCM CGE + TCM CGE + TCM
Future projection[8][9] [9] [14]
Climate model--MIROC5MIROC5MRI-CGCM3HadGEM2-ES21 models
Average adaptation cost (USD/m2)182.76215.96 215.96/182.76 215.96/182.76
SLR: 30 cm 17----------
SLR: 65 cm 20----------
RCP2.62031–2050--24/44/44/4
2081–2100--16/75/56/66/6
RCP4.52031–2050----4/44/44/4--
2081–2100----6/86/66/76/7
RCP6.02031–2050------------
2081–2100----------6/7
RCP8.52031–2050--24/44/44/4--
2081–2100--18/98/98/98/9
Unit: The number of cost-effective prefectures.

Share and Cite

MDPI and ACS Style

Nakajima, K.; Sakamoto, N.; Udo, K.; Takeda, Y.; Ohno, E.; Morisugi, M.; Mori, R. Cost-Benefit Analysis of Adaptation to Beach Loss Due to Climate Change in Japan. J. Mar. Sci. Eng. 2020, 8, 715. https://doi.org/10.3390/jmse8090715

AMA Style

Nakajima K, Sakamoto N, Udo K, Takeda Y, Ohno E, Morisugi M, Mori R. Cost-Benefit Analysis of Adaptation to Beach Loss Due to Climate Change in Japan. Journal of Marine Science and Engineering. 2020; 8(9):715. https://doi.org/10.3390/jmse8090715

Chicago/Turabian Style

Nakajima, Kazunori, Naoki Sakamoto, Keiko Udo, Yuriko Takeda, Eiji Ohno, Masafumi Morisugi, and Ryuta Mori. 2020. "Cost-Benefit Analysis of Adaptation to Beach Loss Due to Climate Change in Japan" Journal of Marine Science and Engineering 8, no. 9: 715. https://doi.org/10.3390/jmse8090715

APA Style

Nakajima, K., Sakamoto, N., Udo, K., Takeda, Y., Ohno, E., Morisugi, M., & Mori, R. (2020). Cost-Benefit Analysis of Adaptation to Beach Loss Due to Climate Change in Japan. Journal of Marine Science and Engineering, 8(9), 715. https://doi.org/10.3390/jmse8090715

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop