Reductions in CO 2 Emissions from Passenger Cars under Demography and Technology Scenarios in Japan by 2050

: Climate policy requires substantial reductions in long-term greenhouse gas (GHG) emissions, including in the transportation sector. As passenger cars are one of the dominant CO 2 emitters in the transport sector, governments and the automobile industry have implemented various countermeasures, including decarbonization of fuels, more energy e ﬃ cient vehicles, and transport demand management. However, the total impact of these measures in the long term remains unclear. This study aims to clarify the CO 2 emissions reductions from passenger cars by 2050 in 1727 municipalities in Japan under a declining population. To estimate CO 2 emissions, we model travel behavior and tra ﬃ c situations reﬂecting the regional conditions of the municipalities, including population density and accessibility to public transport for the base year 2010. Assuming plausible scenarios for future populations and automobile technologies, we estimate CO 2 emissions from passenger cars. We estimate that CO 2 emissions will decline by 64–70% between 2010 and 2050, with automobile technologies playing the largest role. We ﬁnd that the impact of urban compaction is marginal at the national level but varies by municipality. These results imply that, given regional variations, all countermeasures, including technology and demand management, must be used to achieve the long-term target of CO 2 emissions reductions.


Introduction
The transport sector accounted for 23% of global CO 2 emissions in 2010, and its share of emissions is expected to increase at a higher rate than that of other sectors towards 2050 [1]. Road transport is a dominant source of emissions, accounting for 70% of transport sector emissions. Furthermore, 80% of road transport emissions come from passenger transport [2]. Consequently, substantial reductions in CO 2 emissions from passenger cars are required, and various countermeasures are expected to reduce these emissions. The countermeasures include reducing carbon intensity [3,4], improving the energy efficiency of vehicles [5][6][7][8], electrification of vehicles [9,10], and reducing vehicle weights [11], as well as demand management measures, such as traffic control and modal shift policies [12][13][14][15] and urban compaction [16][17][18][19][20][21].
The technological measures are expected to be effective everywhere if the technologies are diffused effectively in the market, whereas the impact of demand management measures varies depending on the regional conditions. The International Energy Agency (IEA) estimated that technological measures would reduce greenhouse gas (GHG) emissions in the transport sector by 54% from 2015 to 2060 under the "2 degree scenario" and by 83% for the same period in the "beyond 2 degree scenario" [22], in which a very ambitious technology deployment rate is assumed. In regard to demand management, automobile and air traffic by 25% globally by 2050. Their analysis was based on national level statistics and did not consider subnational conditions. Using Japanese statistics for 1980-2005, Matsuhashi and Ariga [24] clarified that an area with a higher population density emits less CO2 from automobiles. On the basis of this statistical relationship, they estimated CO2 emissions from automobiles for whole municipalities in Japan in 2030. They estimated that urban compaction would reduce CO2 emissions by 5% compared with a sprawled urban pattern. However, they did not consider the improvement of automobile technologies.
Studies analyzing mitigation impacts on a global scale usually simplify the transport demand conditions, and most do not consider spatial policies such as urban compaction. Meanwhile, the studies analyzing the impact of spatial policies usually do not consider vehicle technologies. To estimate the future potential to reduce CO2 emissions, we need to incorporate both technology and demand management. In this study, we estimate the CO2 emissions from passenger cars in 1727 municipalities in Japan by 2050 and analyze the factors affecting the emissions, considering both automobile technologies and travel demand. On the basis of this analysis, we clarify the impact of the scenarios on automobile CO2 emissions quantitatively and discuss the possible measures for substantial emissions reductions. The impact of the factors influencing emissions under depopulation, demographic change, and spatial shift of activities is not clarified sufficiently in the existing literature. Our research advances the CO2 mitigation studies in the transport sector by considering automobile technologies, demography, and spatial feature of regions simultaneously.

Analytical Framework and Method
Our analytical framework to estimate CO2 emissions is based on breaking down emissions into activity (A), modal structure (S), modal energy intensity (I), and carbon content of fuels (F), known as the ASIF approach [25,26]. This approach analyzes the emissions reduction potential of each factor and the interaction among the factors. In this study, we consider the energy sources as gasoline, diesel, electricity, and hydrogen. The well-to-wheel carbon contents of these sources are given by scenarios based on the literature. The energy efficiencies of vehicles by energy sources are also given by the scenarios. In the real world, emissions are affected by road speed. We quantify the regional road speed distribution based on traffic volumes and road supply. Traffic volumes are quantified based on demography and urban policies. To estimate the travel demand, we formulate the go out rate and the net trip rate by gender and age classes. Furthermore, we estimate the average occupancy rate of a vehicle and the average car trip length. The explanatory variables of these models are the urban conditions, including population and transportation facilities. We assume that there are three transport modes: cars, public transport, and walking/cycling. The modal shares of each are formulated using the explanatory variables for the urban conditions. The analytical framework is depicted in Figure 1. We denote CO 2 emissions from passenger cars in municipality i in year t as CO2 ti , which is formulated as follows: where x ti is the attribute vector of the municipality, ε t is the average real-world emissions factor of in-use vehicles, L is the trip length of passenger cars, ρ is the occupancy ratio, and N ck is the number of trips of a user of attribute k. N ck (x ti ) is expressed in the form: where tpp i is the total population, S tki is the population share of people with attribute k, gor tki is the go-out rate, ntg tki is net trips per person, and Pr ctki is the modal share of passenger cars. These functions and emission factors are determined in Section 4 using the data explained in the next section. As formulated in Equation (2), we use individual projection for the car travel demand. Some studies projected car travel demand using household projection because the purchase and the usage of cars are usually decisions by households, not individuals. The probable effect of using the individual projection instead of household projection is addressed in Section 6.

Data and Scenarios
In this study, we use the "nationwide person trip survey" [27] for travel pattern data and "mesh statistics and projections" [28] for data on the urban population, its density, and the urban area. Public transport facilities are given by the geographic information database of the Ministry of Land, Infrastructure, and Transport (MLIT) of Japan [28]. Road traffic and speed data are provided by the National Road Census [29]. In the following subsections, we provide a summary of these data and describe future scenarios of urban features and automobile technologies.

Nationwide Person Trip Survey
The nationwide person trip survey is conducted every 5 years and compiles travel behavior data, including trip purpose, travel modes, and number of trips, for selected cities in Japan. We use the data for 2010. A total of 130 cities were surveyed in 2010, including 70 large cities, in each of which 500 households are sampled. The data include age, gender, and travel behavior of the household members. The remaining 60 cities are small cities, in which 50 households were surveyed to collect the same data as for the large cities.

Population Projections and Spatial Distribution
We use the mesh population data [28] for the municipality population projections and spatial distributions. This data source provides the population by attribute for each 500 m × 500 m mesh, covering the entire land area of Japan. The base year of these data is 2010, and the mesh population in the base year is given by the national census. The future mesh populations are projected consistent with the population projections for all municipalities by the National Institute of Population and Social Security Research (IPSS) [30]. The mesh data include population projections by gender and age classes for all 5-year periods up to 2050. For reference, the population trajectory by age class from 1965 to 2065 compiled by the IPSS is shown in Figure 2.
On the basis of these spatial population data, we define the urban area as meshes with a population density of over 40 people/ha and calculate the urbanized area of the municipality. The urban population is the sum of the populations in the urban meshes, and the urban population ratio is the ratio of the urban population to the total population in the municipality. We define the "urban geometric moment" Sustainability 2020, 12, 6919 4 of 20 as the sum of the products of the mesh populations and the mesh distance from the population centroid of a municipality divided by the total population of the municipality. It is used as a representative index of the urban spatial scale.
Sustainability 2020, 12, x FOR PEER REVIEW 4 of 21 geometric moment" as the sum of the products of the mesh populations and the mesh distance from the population centroid of a municipality divided by the total population of the municipality. It is used as a representative index of the urban spatial scale.
Source: prepared by the author using IPSS data [30]

Public Transport Catchment Area
The location coordinates of railway stations and bus stops are obtained from the national land information system [28]. We assume that the catchment area of a bus stop (railway station) is a circular area within a 300 m (800 m) radius from the bus stop (station) based on the definitions provided by the MLIT. We calculate the catchment population as the sum of the population of the meshes that intersect with the catchment area. The merged catchment area of a bus stop and railway station is defined as the public transport catchment area. Figure 3 shows the catchment area in the case of the town of Chippubetsu, Hokkaido. The black polygon shows the administrative area of the town, the gray edged meshes indicate the out-of-catchment area, the small red circles show the bus stop catchment areas, the large red circles shows the rail station catchment areas, and the blue edged meshes intersecting with the catchment circles are defined as the meshes of the public transport catchment areas.

Road Traffic Census
The National Road Census [29] provides publicly available road statistics. The data for 2015 include traffic volume, peak traffic ratio, speeds for congested and noncongested travel time, road link length, and number of lanes for approximately 97,000 road links. Owing to data limitations, it is

Public Transport Catchment Area
The location coordinates of railway stations and bus stops are obtained from the national land information system [28]. We assume that the catchment area of a bus stop (railway station) is a circular area within a 300 m (800 m) radius from the bus stop (station) based on the definitions provided by the MLIT. We calculate the catchment population as the sum of the population of the meshes that intersect with the catchment area. The merged catchment area of a bus stop and railway station is defined as the public transport catchment area. Figure 3 shows the catchment area in the case of the town of Chippubetsu, Hokkaido. The black polygon shows the administrative area of the town, the gray edged meshes indicate the out-of-catchment area, the small red circles show the bus stop catchment areas, the large red circles shows the rail station catchment areas, and the blue edged meshes intersecting with the catchment circles are defined as the meshes of the public transport catchment areas.
Sustainability 2020, 12, x FOR PEER REVIEW 4 of 21 geometric moment" as the sum of the products of the mesh populations and the mesh distance from the population centroid of a municipality divided by the total population of the municipality. It is used as a representative index of the urban spatial scale.
Source: prepared by the author using IPSS data [30]

Public Transport Catchment Area
The location coordinates of railway stations and bus stops are obtained from the national land information system [28]. We assume that the catchment area of a bus stop (railway station) is a circular area within a 300 m (800 m) radius from the bus stop (station) based on the definitions provided by the MLIT. We calculate the catchment population as the sum of the population of the meshes that intersect with the catchment area. The merged catchment area of a bus stop and railway station is defined as the public transport catchment area. Figure 3 shows the catchment area in the case of the town of Chippubetsu, Hokkaido. The black polygon shows the administrative area of the town, the gray edged meshes indicate the out-of-catchment area, the small red circles show the bus stop catchment areas, the large red circles shows the rail station catchment areas, and the blue edged meshes intersecting with the catchment circles are defined as the meshes of the public transport catchment areas.

Road Traffic Census
The National Road Census [29] provides publicly available road statistics. The data for 2015 include traffic volume, peak traffic ratio, speeds for congested and noncongested travel time, road link length, and number of lanes for approximately 97,000 road links. Owing to data limitations, it is

Road Traffic Census
The National Road Census [29] provides publicly available road statistics. The data for 2015 include traffic volume, peak traffic ratio, speeds for congested and noncongested travel time, road link length, and number of lanes for approximately 97,000 road links. Owing to data limitations, it is Sustainability 2020, 12, 6919 5 of 20 difficult to estimate the impact of traffic on road speed for each link. As there are many observed links in this data set, we expect to capture the regional features of road speed and traffic by using the regionally pooled data.

Future Scenarios
The major drivers affecting CO 2 emissions include population change, urban structure change, and the progress of automobile technologies. In this study, the urban structure is defined by the spatial distribution of the population and, therefore, population change induces urban structure change. The baseline scenario of the spatial distribution of the population is given by the MLIT. We generate a compact scenario for the future urban structure. In this compact scenario, the total population for each municipality is the same as in the baseline scenario, but we assume that the populations decline in the meshes of the municipalities that have the smallest populations. In other words, the compact scenario assumes that the population becomes zero for the meshes with the smallest populations to meet the given municipality population decline. This is an extreme case of urban compaction. Figure 4 shows the spatial distribution of the population in 2010 and baseline and compact scenarios in 2050. In the case of the baseline scenario, the meshes of the low population density areas are distributed over a wide area. By contrast, in the compact scenario, we observe that the population meshes are located within a very limited area.
Sustainability 2020, 12, x FOR PEER REVIEW 5 of 21 difficult to estimate the impact of traffic on road speed for each link. As there are many observed links in this data set, we expect to capture the regional features of road speed and traffic by using the regionally pooled data.

Future Scenarios
The major drivers affecting CO2 emissions include population change, urban structure change, and the progress of automobile technologies. In this study, the urban structure is defined by the spatial distribution of the population and, therefore, population change induces urban structure change. The baseline scenario of the spatial distribution of the population is given by the MLIT. We generate a compact scenario for the future urban structure. In this compact scenario, the total population for each municipality is the same as in the baseline scenario, but we assume that the populations decline in the meshes of the municipalities that have the smallest populations. In other words, the compact scenario assumes that the population becomes zero for the meshes with the smallest populations to meet the given municipality population decline. This is an extreme case of urban compaction. Figure 4 shows the spatial distribution of the population in 2010 and baseline and compact scenarios in 2050. In the case of the baseline scenario, the meshes of the low population density areas are distributed over a wide area. By contrast, in the compact scenario, we observe that the population meshes are located within a very limited area. In this study, we consider five automobile technologies that have different CO2 emissions factors: (1) gasoline, (2) diesel, (3) hybrid electric vehicles (HEV), (4) electric vehicles/plug-in hybrid electric vehicles (EV/PHEV), and (5) fuel cell vehicles (FCV). The share of vehicles in-use by these technologies until 2030 are determined based on the "roadmap for EV/PHEV" [31] and the "new strategic roadmap for hydrogen and fuel cells" [32]. We calculate the sales share of the vehicles by technologies to ensure consistency with the in-use share, assuming a survival curve for automobiles under which their maximum duration is 20 years. After 2030, the sales share by technologies is determined based on the 2 °C Scenario (2DS) improve scenario in IEA Energy Technology Perspectives (ETP) 2012 [33]. Applying the survival curve to the sales scenario, we estimate the share of in-use automobiles by technologies. Figure 5 shows the sales share and the in-use vehicle shares for the five technologies. In this study, we consider five automobile technologies that have different CO 2 emissions factors: (1) gasoline, (2) diesel, (3) hybrid electric vehicles (HEV), (4) electric vehicles/plug-in hybrid electric vehicles (EV/PHEV), and (5) fuel cell vehicles (FCV). The share of vehicles in-use by these technologies until 2030 are determined based on the "roadmap for EV/PHEV" [31] and the "new strategic roadmap for hydrogen and fuel cells" [32]. We calculate the sales share of the vehicles by technologies to ensure consistency with the in-use share, assuming a survival curve for automobiles under which their maximum duration is 20 years. After 2030, the sales share by technologies is determined based on the 2 • C Scenario (2DS) improve scenario in IEA Energy Technology Perspectives (ETP) 2012 [33]. Applying the survival curve to the sales scenario, we estimate the share of in-use automobiles by technologies. Figure 5 shows the sales share and the in-use vehicle shares for the five technologies.
"The Plan for Global Warming Countermeasures" [34] sets a target for average fuel economy by 2030. We determine the improvement in fuel economy for new gasoline and diesel vehicles by referring to this plan. We assume that the fuel economies of gasoline and diesel technologies reach their upper limit in 2030 and do not change after this. On the basis of this assumption, we calculate the average fuel economy of the in-use vehicles. We fix the energy efficiency of the other technologies at the 2010 level, but the emissions from electric power generation are assumed to improve, reflecting the voluntary target of 370 g CO 2 /kWh by 2030 set by the Federation of Electric Power Companies of Japan [35]. After 2030, we have no information on the emissions reduction target for electric power generation, thus we simply fix the emissions factor in our scenario. For real-world fuel efficiency, we assume that HEV and EV/PHEV consume 67% more energy than is the case for test cycle fuel efficiency and that the other technologies, gasoline, diesel, and FCV, consume 43% more. The assumed CO 2 emission factors by automobile technologies and average emission factors of in-use vehicles are summarized in Table 1. "The Plan for Global Warming Countermeasures" [34] sets a target for average fuel economy by 2030. We determine the improvement in fuel economy for new gasoline and diesel vehicles by referring to this plan. We assume that the fuel economies of gasoline and diesel technologies reach their upper limit in 2030 and do not change after this. On the basis of this assumption, we calculate the average fuel economy of the in-use vehicles. We fix the energy efficiency of the other technologies at the 2010 level, but the emissions from electric power generation are assumed to improve, reflecting the voluntary target of 370 g CO2/kWh by 2030 set by the Federation of Electric Power Companies of Japan [35]. After 2030, we have no information on the emissions reduction target for electric power generation, thus we simply fix the emissions factor in our scenario. For real-world fuel efficiency, we assume that HEV and EV/PHEV consume 67% more energy than is the case for test cycle fuel efficiency and that the other technologies, gasoline, diesel, and FCV, consume 43% more. The assumed CO2 emission factors by automobile technologies and average emission factors of in-use vehicles are summarized in Table 1.

Model
In this section, we explain the models that we construct to estimate traffic volume and CO2 emissions, using the data for 130 municipalities described in Section 3. For each municipality, we estimate: (1) trip generation, (2) modal split, (3) vehicle occupancy ratio, and (4) trip length by car. Multiplying the estimated vehicle-kms and CO2 emissions factors, we estimate the total volume of CO2 emissions. Using this model, we estimate the future emissions for every 5 years from 2020 to 2050 and decompose the emissions by factors of technologies, road speed, population, demography, and urban structure under the scenarios explained above.

Model Estimation
First, we model the go-out rate and the net trip generation by gender and age classes using linear regression analysis. Second, the modal split is represented by a nested logit model. The lower level model is the choice between car and public transport, which is in the motorized transport nest. The

Model
In this section, we explain the models that we construct to estimate traffic volume and CO 2 emissions, using the data for 130 municipalities described in Section 3. For each municipality, we estimate: (1) trip generation, (2) modal split, (3) vehicle occupancy ratio, and (4) trip length by car. Multiplying the estimated vehicle-kms and CO 2 emissions factors, we estimate the total volume of CO 2 emissions. Using this model, we estimate the future emissions for every 5 years from 2020 to 2050 and decompose the emissions by factors of technologies, road speed, population, demography, and urban structure under the scenarios explained above.

Model Estimation
First, we model the go-out rate and the net trip generation by gender and age classes using linear regression analysis. Second, the modal split is represented by a nested logit model. The lower level model is the choice between car and public transport, which is in the motorized transport nest. The upper level model is the choice between motorized and nonmotorized (walking/cycling) transport. The third and the fourth models are linear regression models used to model vehicle occupancy ratio and car trip length, respectively. All the above regression models are estimated using stepwise regression and the backward elimination approach to maximize the Akaike information criterion.

Go-Out Rate Model
By stepwise regression, we obtain the following formula to estimate the go-out rate gor tki of people with attribute k in municipality i in year t: where tpp ti is the total population, urbA ti is the area of the urban area, popD ti is population density, and α gk0 , α gk1 , α gk2 , and α gk3 are parameters. The estimated parameters are shown in Table 2. NA indicates rejected variables.
The intercept is the go-out rate when all other variables are zero. This parameter is lower for elderly people and females. The parameters of the total population are significantly positive except for males aged 65-74 years and females over 75 years. This means that most classes of people tend to go out from their homes in larger cities.

Net Trip Generation Model
Net trip generation refers to the average number of trips of a person who goes out from home. The net trip number ntg tki is described as follows: where α nk0 , α nk1 , α nk2 , and α nk3 are parameters. The estimated parameters are shown in Table 3.
The intercept indicates the net trip number when all explanatory variables are zero. The intercept is highest for males aged 65-74 years. There are no clear patterns in the explanatory variables based on the attributes of people. For instance, the number of trips is larger when the municipality population is larger for groups of males aged 15-64 years, males over 75 years, females less than 65 years, and females over 75 years. However, the variable does not affect males and females aged 65-74 years, and it reduces the net trips of males less than 15 years. We estimate the aggregated modal share using the following nested logit model: where Pr ctki is the modal share of cars, V ctki and V wtki are the utilities of driving a car and of walking/cycling, respectively, pBs ti and pRs ti are the populations of the bus stop and the railway station catchment areas, respectively, and ups ti is the urban population share. IV kti is an inclusive value of the lower level model expressed as IV kti = log(1 + exp(V ctki )). Table 4 shows the estimated parameters of the lower level model, the choice set of which is cars and public transport. The intercept indicates the car dummy, and it is positive for all categories of people. This means that people tend to choose cars more than public transport. Among the explanatory variables, the urban population share is selected for attributes of males over 15 years and females between 15 and 74 years. This parameter is negative for these people. This means that the public transport share tends to be larger when the urban population share of the municipality is larger. For males and females aged 14 years or less, the population density parameter has a negative value. This means that public transport tends to be chosen more by these young people as population density becomes higher. These young people do not drive themselves but would be passengers accompanied by other people.  Table 5 shows the estimated parameters of the upper level model. The dummy parameter of walking/cycling is negative except for people aged 14 years and under. In other words, people over 15 years tend to choose motorized transport, and people under 14 tend to choose walking/cycling. The effects of the other explanatory variables differ by people's attributes. The inclusive value (IV) is only selected for males aged 15-64 years. The sign of the parameter is negative, which means that as the utility of the motorized mode becomes higher, the probability of choosing walking/cycling becomes lower. The signs of the other explanatory variables are positive, which means that as the population within the bus catchment area becomes higher and the population density increases, the share of those choosing walking/cycling becomes higher. Except for the intercept parameter, all variables are rejected Sustainability 2020, 12, 6919 9 of 20 for females aged over 75 years, indicating that none of the explanatory variables affect the choice of walking/cycling for this group of people.

Vehicle Occupancy Rate Model
The vehicle occupancy rate ρ ti is estimated using the following form: where ugm ti is an urban geometric moment, as defined above. The estimated parameters, shown in Table 6, indicate that the occupancy rate becomes higher as the total population becomes larger, the urban geographic moment becomes smaller, and the urban population share becomes higher. Average car trip length is estimated using the following form: where tap ti is the population in the public transport catchment area, which is described as tap ti = tpp ti × S tr,ti , and S tr,ti is the population share of the public transport catchment area. The estimated parameters of the model are shown in Table 7. From the table, the base car trip length is about 11.8 km. The trip length is shorter when the urban population share is higher and the population of the public transport catchment area is larger. Intuitively, the population of the public transport catchment area seems to have no relation with car trip length, but this index may be used as a proxy variable of urban size. The MLIT reported that the CO 2 emissions per vehicle-km ε are a function of the average travel speed v in the form: where the unit of ε is gCO 2 /km and v is km/h. The parameters for gasoline and diesel passenger cars in 2000 are given in Table 8 [36].  Figure 6 shows the emissions curve over the road speed for gasoline vehicles using Equation (10). The emissions factor takes a minimum value around 60-70 km/h and higher values at lower speed ranges. This nonlinearity of the emissions factor requires consideration of road congestion in the emissions estimation. If we use average speed only, then we would underestimate the CO 2 emissions.  If the road census has sufficient sample links in a municipality, then the speed distribution of the municipality can be obtained. Figure 7 shows the road speed distribution of Tokyo, Nagoya, and Osaka compiled using speed and traffic data from the national road traffic census. The horizontal axis indicates the road speed, and the vertical axis indicates the total traffic volume of the speed range. In these three cities, the peak of the distribution is around 20 km/h. The red dashed line is the lognormal probability density function fitted to the data using the maximum likelihood method. These figures indicate that the log-normal distribution can represent the road speed distribution. Hereafter, therefore, we apply the log-normal distribution to represent the speed distribution for all municipalities. As shown later, the effect of a road speed shift on the emissions factor is standardized by using the change rate of the emissions factor according to the road speed shift. Therefore, the effect has no dimension, and the level of vehicle emissions in Equation (10) does not affect future emissions.
If the road census has sufficient sample links in a municipality, then the speed distribution of the municipality can be obtained. Figure 7 shows the road speed distribution of Tokyo, Nagoya, and Osaka compiled using speed and traffic data from the national road traffic census. The horizontal axis indicates the road speed, and the vertical axis indicates the total traffic volume of the speed range. In these three cities, the peak of the distribution is around 20 km/h. The red dashed line is the log-normal probability density function fitted to the data using the maximum likelihood method. These figures indicate that the log-normal distribution can represent the road speed distribution. Hereafter, therefore, we apply the log-normal distribution to represent the speed distribution for all municipalities.
If the road census has sufficient sample links in a municipality, then the speed distribution of the municipality can be obtained. Figure 7 shows the road speed distribution of Tokyo, Nagoya, and Osaka compiled using speed and traffic data from the national road traffic census. The horizontal axis indicates the road speed, and the vertical axis indicates the total traffic volume of the speed range. In these three cities, the peak of the distribution is around 20 km/h. The red dashed line is the lognormal probability density function fitted to the data using the maximum likelihood method. These figures indicate that the log-normal distribution can represent the road speed distribution. Hereafter, therefore, we apply the log-normal distribution to represent the speed distribution for all municipalities. We estimate the parameters of the log-normal distribution for each municipality to fit with its speed distribution using the maximum likelihood method. If traffic volume and road supply affect the speed distribution, then we can observe the effect of them on the estimated parameters. Figure 8 shows the plot of two parameters of the log-normal distribution over the average lane traffic of municipalities. We estimate the parameters of the log-normal distribution for each municipality to fit with its speed distribution using the maximum likelihood method. If traffic volume and road supply affect the speed distribution, then we can observe the effect of them on the estimated parameters. Figure 8 shows the plot of two parameters of the log-normal distribution over the average lane traffic of municipalities. In this figure, each point indicates the lane traffic and the estimated parameter of a municipality, and the green line is the regression line. μ is the mean parameter and it tends to decrease as lane traffic increases. σ is the standard deviation parameter, which tends to increase as lane traffic increases. This means that as the average lane traffic becomes larger, the speed distribution of a municipality becomes slower and deviates more. The estimated parameters of the regression line are shown in Table 9. Using the relationship among traffic, road speed, and the emissions factor, we can estimate the impact of traffic volume on the average CO2 emissions factor ε. As the emissions factor reflects the progress of automobile technology, we need to distinguish these effects.
First, we consider the case where automobile technology is fixed. We denote p(V) as a probability density function of traffic speed V and f(V) as the emissions factor at speed V. Then, the average emission factor is expressed as follows: In this figure, each point indicates the lane traffic and the estimated parameter of a municipality, and the green line is the regression line. µ is the mean parameter and it tends to decrease as lane traffic increases. σ is the standard deviation parameter, which tends to increase as lane traffic increases. This means that as the average lane traffic becomes larger, the speed distribution of a municipality becomes slower and deviates more. The estimated parameters of the regression line are shown in Table 9. Using the relationship among traffic, road speed, and the emissions factor, we can estimate the impact of traffic volume on the average CO 2 emissions factor ε. As the emissions factor reflects the progress of automobile technology, we need to distinguish these effects.
First, we consider the case where automobile technology is fixed. We denote p(V) as a probability density function of traffic speed V and f (V) as the emissions factor at speed V. Then, the average emission factor is expressed as follows: (11) Traffic change affects the speed distribution p(V), which is characterized by µ and σ. Of course, the speed distribution is affected not only by the traffic volume but also by the other factors, including density of intersections, land use, and landscape. Therefore, the relation between traffic volume and the speed distribution parameters deviates, as shown in Figure 8. In this study, we assume that the factors other than the traffic volume are fixed for each municipality and that the sensitivities of the traffic volume to these parameters are the same for all municipalities. Then, the speed distribution parameters can be calculated by the following equations using the average lane traffic q i for municipality i: Here, the suffix 0 indicates the base year, 2010, and α µ and α σ are the slopes of the regression lines in Figure 8. Denoting the emissions factor by Equation (11) at time t as ε t ri and the estimated emissions factor with fixed technology at time t + 1 as ε t+1 ri , the change rate of the emissions factor by the traffic speed change α t+1 ri is given as: We denote the change rate of the emissions factor by automobile technology from t to t + 1 as α t+1 T . Then, the emissions factor at time t + 1, considering both traffic volume and technology, can be expressed as follows: Here, we assume that technological progress is disseminated to all municipalities at once. Substituting Equation (15) into Equation (1), we can estimate the CO 2 emissions reflecting technological progress and travel speed change.

Effect of Technological Factor
The primary approximation of the change in CO 2 emissions in response to the technological factor is expressed in the following form:

Effect of Road Speed Change
The change rate of the emission factor by road speed change α t+1 ri differs by municipality. Therefore, the total CO 2 emissions change by road speed change is expressed as follows: Sustainability 2020, 12, 6919 13 of 20

Effect of the Total Population
The variable of total population tpp ti appears in L, ρ, and N. Therefore, the primary approximation of the change in CO 2 emissions by the total population change is derived as follows: where: ∂N ctki ∂tpp ti = S tki ·gor tki ·ntg tki ·Pr ctki +tpp i ·S tki ·α gk1 ·ntg tki ·Pr ctki +tpp i ·S tki ·gor tki ·α nk1 ·Pr ctki (21)

Effect of Demographic Composition
The demographic composition can be represented by the share of population attributes S tki . The effect of a change of this share on CO 2 emissions can be expressed as follows:

Results
Using the above model, we estimate the total CO 2 emissions in 2010 to be 104 MTCO 2 . Comparing this with the observed volume of 106 MTCO 2 , our result underestimates emissions by around 1.6%. Although there are no official statistics on CO 2 emissions from passenger cars by municipality, the Ministry of Environment (MoE) of Japan provides an estimate for 2010 [37]. For the 1555 municipalities that appear both in our data and in the MoE estimate, we find that the correlation between our and the MoE's estimates is 0.976. Thus, despite the approaches differing, the estimates are highly correlated. Given the high accuracy of the total CO 2 emissions estimation, we can say that our approach provides an adequate estimation of CO 2 emissions from passenger cars. Figure 9 shows the estimated national total CO 2 emission from 2010 to 2050 and the contribution of various factors to emission changes. In the figure, the red line shows the trajectory of estimated total CO 2 emissions, which begin at 104 MTCO 2 in 2010 and are estimated to fall to 32.7 MTCO 2 in 2050, a 64% reduction. The technology factor contributes most to this decline, accounting for 43.6% of the reduction in emissions in 2050. The second largest contributor is the total population, which accounts for 19.7% of the emissions reduction, followed by demography, which contributes 3.4%. The urban structure factor, shown by the green line, is estimated to result in a 4.3% increase in total emissions. This means that, under the baseline scenario, suburbanization induces longer trips and greater car usage, both of which increase emissions. The dark gray area indicates the speed factor, which, by reducing congestion, decreases emissions by 1.9%.

CO2 Emissions Changes by Factors Consistent with Total CO2 Emissions Change
The above formulation of CO2 emissions change by factor X is a primary approximation (for convenience, we denote it as ΔCO2| ), and the sum of the formulas does not generally equal the total emissions reductions estimated by Equation (1). We denote the adjusted CO2 emissions change by factor X as ΔCO2 | and define it by the following equation:

Results
Using the above model, we estimate the total CO2 emissions in 2010 to be 104 MTCO2. Comparing this with the observed volume of 106 MTCO2, our result underestimates emissions by around 1.6%. Although there are no official statistics on CO2 emissions from passenger cars by municipality, the Ministry of Environment (MoE) of Japan provides an estimate for 2010 [37]. For the 1555 municipalities that appear both in our data and in the MoE estimate, we find that the correlation between our and the MoE's estimates is 0.976. Thus, despite the approaches differing, the estimates are highly correlated. Given the high accuracy of the total CO2 emissions estimation, we can say that our approach provides an adequate estimation of CO2 emissions from passenger cars. Figure 9 shows the estimated national total CO2 emission from 2010 to 2050 and the contribution of various factors to emission changes. In the figure, the red line shows the trajectory of estimated total CO2 emissions, which begin at 104 MTCO2 in 2010 and are estimated to fall to 32.7 MTCO2 in 2050, a 64% reduction. The technology factor contributes most to this decline, accounting for 43.6% of the reduction in emissions in 2050. The second largest contributor is the total population, which accounts for 19.7% of the emissions reduction, followed by demography, which contributes 3.4%. The urban structure factor, shown by the green line, is estimated to result in a 4.3% increase in total emissions. This means that, under the baseline scenario, suburbanization induces longer trips and greater car usage, both of which increase emissions. The dark gray area indicates the speed factor, which, by reducing congestion, decreases emissions by 1.9%.   Figure 10 shows the estimated CO 2 emissions and the contribution of the factors under the compact scenario. The total CO 2 emissions are estimated at 30.6 MTCO 2 in 2050, which is a 70% reduction from the 2010 level. The reduction rate under the compact scenario is 6% higher than that of the baseline scenario. In 2050, the technology factor contributes the largest share of this reduction at 40.6%, followed by the total population factor (17.8%), the urban structure (6.6%), the demography (3.1%), and the road speed factor (2.6%). Under the baseline scenario, the urban structure results in a 4.3% increase in emissions and, thus, the urban compaction scenario results in a substantial reduction in the emissions caused by the urban structure factor. Furthermore, the road speed factor reduces emissions more under the compact scenario than under the baseline scenario. In this study, urban compaction reduces the share of car usage and the vehicle travel length, while the road provision is assumed to be fixed. Therefore, there is less average lane traffic and a higher road speed under the compact scenario. It should be noted that the result would differ if urban compaction resulted in road traffic concentrating at certain links and inducing congestion.
Sustainability 2020, 12, x FOR PEER REVIEW 16 of 21 Figure 10 shows the estimated CO2 emissions and the contribution of the factors under the compact scenario. The total CO2 emissions are estimated at 30.6 MTCO2 in 2050, which is a 70% reduction from the 2010 level. The reduction rate under the compact scenario is 6% higher than that of the baseline scenario. In 2050, the technology factor contributes the largest share of this reduction at 40.6%, followed by the total population factor (17.8%), the urban structure (6.6%), the demography (3.1%), and the road speed factor (2.6%). Under the baseline scenario, the urban structure results in a 4.3% increase in emissions and, thus, the urban compaction scenario results in a substantial reduction in the emissions caused by the urban structure factor. Furthermore, the road speed factor reduces emissions more under the compact scenario than under the baseline scenario. In this study, urban compaction reduces the share of car usage and the vehicle travel length, while the road provision is assumed to be fixed. Therefore, there is less average lane traffic and a higher road speed under the compact scenario. It should be noted that the result would differ if urban compaction resulted in road traffic concentrating at certain links and inducing congestion.  The municipalities in mountain areas have smaller populations and emit fewer emissions compared with urban municipalities in 2010, and the total volumes of emissions reduction is smaller in the future. The per capita reductions in volume of CO2 emissions are smaller in large cities and larger in small cities. In small cities, the share of car usage is high, making the impact of emission reduction technologies high. Conversely, in large cities where the modal share of public transport is relatively high, the impact of automobile emission reduction technologies on per capita CO2 emissions is relatively small. The volume of CO2 emission reductions per vehicle-km is high in large cities. The road speed in large cities is relatively slow as a result of congestion and, therefore, technological improvements have a stronger effect in terms of reducing the absolute volume of emissions. At the same time, road speed will be faster and the emissions factor will be improved when traffic is reduced. Combined with these factors, the CO2 emission per vehicle-km is reduced more in large cities than in small cities. In summary, the impact of CO2 emission reduction differs by region, and different indices of emission changes correspond to different regional features. Details of these regional differences are provided in the Supplementary Materials.  The municipalities in mountain areas have smaller populations and emit fewer emissions compared with urban municipalities in 2010, and the total volumes of emissions reduction is smaller in the future. The per capita reductions in volume of CO 2 emissions are smaller in large cities and larger in small cities. In small cities, the share of car usage is high, making the impact of emission reduction technologies high. Conversely, in large cities where the modal share of public transport is relatively high, the impact of automobile emission reduction technologies on per capita CO 2 emissions is relatively small. The volume of CO 2 emission reductions per vehicle-km is high in large cities. The road speed in large cities is relatively slow as a result of congestion and, therefore, technological improvements have a stronger effect in terms of reducing the absolute volume of emissions. At the same time, road speed will be faster and the emissions factor will be improved when traffic is reduced. Combined with these factors, the CO 2 emission per vehicle-km is reduced more in large cities than in small cities. In summary, the impact of CO 2 emission reduction differs by region, and different indices of emission changes correspond to different regional features. Details of these regional differences are provided in the Supplementary Materials. Figure 12 shows the differences in the three indices between baseline and compact scenarios in 2050. The compact scenario derives lower CO 2 emissions than does the baseline scenario, and its impact is higher in large cities and lower in small cities. In other words, urban compaction is effective in reducing CO 2 emissions in urban areas but less effective in rural areas. The effect of urban compaction in reducing CO 2 emissions per capita is relatively weak at the center of large cities, but it is strong in local cities and in the suburbs of large cities. Population density in the center of large cities is already high, and the effect of urban compaction may be limited there. Local cities or suburbs of large cities have more capacity to reduce emissions per capita by urban compaction. The difference of emissions per vehicle-km indicates that urban compaction boosts reductions in emissions. In this study, urban compaction reduces the traffic volume by more than the baseline scenario and, as a result, traffic congestion is alleviated further. Therefore, the emissions per vehicle-km are reduced more in a compact scenario.
Sustainability 2020, 12, x FOR PEER REVIEW 17 of 21 Figure 12 shows the differences in the three indices between baseline and compact scenarios in 2050. The compact scenario derives lower CO2 emissions than does the baseline scenario, and its impact is higher in large cities and lower in small cities. In other words, urban compaction is effective in reducing CO2 emissions in urban areas but less effective in rural areas. The effect of urban compaction in reducing CO2 emissions per capita is relatively weak at the center of large cities, but it is strong in local cities and in the suburbs of large cities. Population density in the center of large cities is already high, and the effect of urban compaction may be limited there. Local cities or suburbs of large cities have more capacity to reduce emissions per capita by urban compaction. The difference of emissions per vehicle-km indicates that urban compaction boosts reductions in emissions. In this study, urban compaction reduces the traffic volume by more than the baseline scenario and, as a result, traffic congestion is alleviated further. Therefore, the emissions per vehicle-km are reduced more in a compact scenario.

Discussion and Conclusions
In this study, we estimated that CO2 emissions from passenger cars in 2050 are reduced by 64% and 70% under baseline and compact scenarios, respectively, compared with 2010 emissions. In 2016, the Japanese government set a long-term target of an 80% reduction in GHG emissions by 2050 in the  Figure 12 shows the differences in the three indices between baseline and compact scenarios in 2050. The compact scenario derives lower CO2 emissions than does the baseline scenario, and its impact is higher in large cities and lower in small cities. In other words, urban compaction is effective in reducing CO2 emissions in urban areas but less effective in rural areas. The effect of urban compaction in reducing CO2 emissions per capita is relatively weak at the center of large cities, but it is strong in local cities and in the suburbs of large cities. Population density in the center of large cities is already high, and the effect of urban compaction may be limited there. Local cities or suburbs of large cities have more capacity to reduce emissions per capita by urban compaction. The difference of emissions per vehicle-km indicates that urban compaction boosts reductions in emissions. In this study, urban compaction reduces the traffic volume by more than the baseline scenario and, as a result, traffic congestion is alleviated further. Therefore, the emissions per vehicle-km are reduced more in a compact scenario.

Discussion and Conclusions
In this study, we estimated that CO2 emissions from passenger cars in 2050 are reduced by 64% and 70% under baseline and compact scenarios, respectively, compared with 2010 emissions. In 2016, the Japanese government set a long-term target of an 80% reduction in GHG emissions by 2050 in the

Discussion and Conclusions
In this study, we estimated that CO 2 emissions from passenger cars in 2050 are reduced by 64% and 70% under baseline and compact scenarios, respectively, compared with 2010 emissions. In 2016, the Japanese government set a long-term target of an 80% reduction in GHG emissions by 2050 in the Plan for Global Warming Countermeasures [34]. In this plan, the sectoral reduction targets are not determined. However, if we apply this target to passenger cars, an additional 10-14% reduction in emissions is required. Depopulation and aging in Japan reduce transport sector activities and consequently reduce CO 2 emissions, with the impact of these changes estimated to be about half of the reduction achieved via technological progress. To achieve the 2050 target, further emissions reduction measures are required.
We demonstrated that the expected population reduction will reduce the population density, which in turn will induce higher car dependency and longer car trip lengths. As a result, the urban structure factor increases the CO 2 emissions in the baseline scenario. Therefore, we need a measure to increase population density under the depopulating situation. As shown in our analysis, even an extremely drastic urban compaction scenario would reduce emissions by only 6%. This result clarifies that control of the urban space solely has a quite limited effect on CO 2 emissions reductions.
Reducing traffic volumes alleviates road congestion in the urban areas, which contributes to reducing CO 2 emissions, but the estimated impact is only 2%. In this study, urban compaction is estimated to reduce vehicle-kms and CO 2 emissions. However, if compaction and the resulting concentration of traffic aggravate congestion, then part of the emission reduction expected from urban compaction would not be achieved. Thus, urban compaction measures need to be combined with other countermeasures, including traffic control, land use management, and public transport provision, to alleviate road traffic concentration.
Our spatial analysis indicated that the per capita CO 2 emission reduction rate is higher for the local cities where the modal share of passenger cars is high. In these cities, improvements in automobile technologies will have a higher impact. The emission reduction effect of urban compaction is high in local cities and in the suburban municipalities of large cities but smaller at the center of large cities and very minor in municipalities in mountainous regions. These results indicate that the effects of emission reduction measures differ by regions. Therefore, the countermeasures of transport demand management should reflect regional features rather than applying uniform national measures. For instance, the combination of road speed improvements and road pricing to curb car traffic demand can be effective in reducing CO 2 emission in large cities. However, they may not be effective in small towns and low-density cities, where people are more dependent on cars and have fewer alternatives and therefore cannot change their travel behavior in response to these measures. In this case, road pricing simply increases user costs but may not reduce CO 2 emissions. Urban compaction may have a small effect in the center of large cities, which are already high density, and in mountainous areas, urban compaction has only a marginal effect owing to the sparse development patterns and low population density. Conversely, urban compaction in local cities or suburbs of large cities that have urban sprawl patterns will be effective in increasing population density and in reducing automobile dependency and car trip length. Demand management measures, including urban compaction policies, impose burdens on people to induce them to alter their behavior. Therefore, it is important to determine an effective framework for emission reduction measures that consider differences in impacts by regions.
In this study, we did not analyze the effect of a carbon intensity reduction by energy source and other demand side management measures, such as teleworking. The emissions factor for electricity is assumed to be fixed after 2030. If the share of renewable energy is increased substantially, then the emissions from using electricity will be reduced. As assumed in this study, if electric vehicles are widely disseminated, then a reduction in the emissions factor for electricity will have a considerable effect on CO 2 emissions from automobiles. Furthermore, if synthetic fuels derived from bioenergy or by a Fischer-Tropsch process can be delivered at a reasonable cost, then the internal combustion engine cars achieve zero well-to-wheel emissions. However, although these technologies may have a substantial impact, their progress is highly uncertain. Transport demand management measures such as teleworking and modal shifts, which may have a large impact on CO 2 emissions reductions, are not fully considered in this study because the feasibility of these measures is also very uncertain. Teleworking is becoming popular among private companies owing to the current infectious disease epidemic, but it is unclear whether this working style will continue. We need to carefully understand the impact of people's behaviors, including their commuting patterns and location choices, that may affect travel behavior.
It should also be noted that our projections of car usage are based on the individual level, not the household level. Some studies use household-level projection because the purchase and the usage of cars are usually the decision by the household [38,39]. Household attributes, such as size, number of children, and income, affect the decision in purchasing cars and usage of cars, especially for accompanied travel. Households are composed of individuals, thus households have a wider variety of types than individuals. The projection of the future composition of household types for all municipalities is another big challenge. In Japan, depopulation and aging are expected to reduce the average household size, but it is unclear how the future shift of household composition affects car usage besides the aging of individual people. It should be addressed in future studies. Additionally, our demography scenario is based on IPSS and MLIT projections, which assumes a relatively conservative number of foreign immigrants (about 70 thousand net immigration annually in 2035). It is also possible to increase the foreign immigrants by the change of the Japanese immigration policy in the future, however, it is unpredictable. Our approach can be applied to the cases of the other demographic scenarios, and the increase of immigration will of course multiply the CO2 emission from passenger cars.
In summary, even with the assumption of substantial diffusion of CO 2 emissions reduction technologies, we find that the emissions reductions from passenger cars will not reach 80% in 2050. We need an additional 10-14% reduction to achieve the 80% target. In addition to the countermeasures considered in this study, we need to apply other measures to achieve the reduction target. There is no panacea of mitigation measures for climate change and all possible measures must be considered simultaneously. This maxim is also applicable in the CO 2 emissions reduction in the passenger car sector.