# Effects of Vehicle Restriction Policies on Urban Travel Demand Change from a Built Environment Perspective

Authors: Xiaoyun Cheng; Kun Huang; Lei Qu; Tianbao Zhang; Li Li
Journal: Journal of Advanced Transportation
 (2020)

---

## Abstract

License plate restriction (LPR) policy presents the most straightforward way to reduce road traffic and emissions worldwide. However, in practice, it has aroused great controversy. This policy broke the original structure of the urban transportation mode, which needed some matching strategies to adapt to this change. Investigating this travel demand change is a challenging task because it is greatly influenced by features of the local built environment. Fourteen variables from four dimensions, location, land-use diversity, distance to transit, and street design, are used to depict the built environment; moreover, the severe collinearity underlies these feature variables. To solve the multicollinearity among the variables and high-dimensional problem, this study utilizes two different penalization-based regression models, the LASSO (least absolute shrinkage and selection operator) and Elastic Net regression algorithms, to achieve the variable selection and explore the impacts of the built environment on the change of travel demand triggered by the LPR policy. Travel demand changes are assessed by the relative variation in taxi ridership in each traffic analysis zone based on the taxi GPS data. Built environment variables are measured using the transportation network data and the Baidu Map Service points of interest (POI) data. The results show that regions with a higher level of public transportation service and a higher degree of the land mix have a stronger resilience to the vehicle restriction policy. Besides, the contribution rate of public transportation is stable as a whole, while the contribution rate of richness depends on specific types of land use. The conclusions in this study can provide in-depth insights into the influence of the LPR policy and underpin traffic complementary policies to ensure the effectiveness of LPR.


---

## Body

## 1. Introduction

In the past several decades, the drastic growth of private car ownership and usage has caused some vexing problems, such as traffic congestion and environmental pollution [1, 2]. Cars are indisputably an important contribution to photochemical smog pollution [3]. To improve the urban living environment, many cities have attempted to implement various traffic demand management (TDM) policies [4] to mitigate these problems. These strategies can be divided by purpose into two categories: (1) the control of private car ownership from the source through the license plate lottery or license plate auction policies implemented in Beijing, Shanghai, and Guangzhou [5] and (2) the reduction of private car usage through congestion pricing [6] and license plate restriction (LPR) policies. For years, LPR policies have been implemented in several cities worldwide, like South America (Bogota) [7] and Mexico City, the most populous city in North America [8]. Some cities with serious “urban disease” caused by air pollution have adopted joint countermeasures. For instance, Beijing implemented a combination of the license plate lottery and LPR policies.The LPR policy, also called road space rationing is the most straightforward way to reduce road traffic and emissions by limiting vehicle access to a road network based on the last digits of the license number on preestablished days. The one-day-per-week (ODPW) vehicle restriction policy does exactly what its name implies by forcing private vehicles off the road “one day per week.” A more stringent policy, the odd-and-even (OAE) policy, only permits vehicles to get on the road alternately based on the odd or even status of the last digit on their license plates. As implemented in Delhi, these two LPR policies are effective in the daytime on weekdays, and the following vehicles are exempted: all taxis and buses as well as cars operating on compressed natural gas (CNG) and electric power [9]. In practice, the OAE policy is usually taken as a “first-aid treatment” or a contingency plan [10] for special events or, during an occurrence of serious air pollution, to quickly reduce traffic volume, which was implemented in Beijing during the Olympic Games in 2008 [11–14] and when Xi’an and Langfang suffered severe air pollution [10]. Up to now, more than 10 cities in China have employed the ODPW policy as a routine TDM strategy and the number will increase with rapid motorization in further.However, in practice, the LPR policy has aroused great controversy and public ire [15]. It was found by managers and researchers that this policy was only effective in mitigating congestion and air pollution over a short period and the long-term effect has fallen short of expectations [16]. The numerical results also confirm the unintended consequences of LPR in many practical cases [17]. This is mainly because the LPR policies may not effectively make commuters transfer from restricted private car usage to public transit [18]. Conversely, such policies could stimulate purchasing a second car [18] and result in more second-hand car usage [19]. These measures used to evade the LPR policy exacerbate traffic congestions and air pollution and consume more parking spaces to accommodate the extra vehicles. Except for increasing the commuters’ acceptance of the LPR policy [16] from a subjective dimension, the primary goal of traffic managers today is to strategize how to guide LPR restricted travel demands towards a higher usage rate of public transit via external policy, because the built environment would affect residents’ behavior and change their social and economic activities [20]. Similarly, when the residents’ private vehicles are restricted, the choice of alternative modes of travel is directly influenced by the built environment. The following considers the impacts of the built environment on travel mode changes under the LPR policy.Three important issues need to be addressed when implementing the LPR policy: (1) What are the spatial characteristics of demand change caused by the LPR policy? (2) How does the built environment influence this demand change in different regions? (3) What targeted policies and strategies should urban and transportation managers take to adapt to changes in travel demand?Despite the rich studies scrutinizing the social and traffic effects of the vehicle restriction policy, related research is rarely reported on the travel demand change under the context of the LPR policy. This is partly because existing research is more focused on the effect of LPR policies on traffic congestion [9, 13] and air pollution mitigation [21]. However, the main obstacles are a lack of panel data to reflect citywide demand changes and an appropriate model to portray the complex relationship between the built environment and demand change.To fill the gap, this paper explores the built environment effect on the traffic mode transfer within the LPR policy context, leveraging taxi GPS trajectory data, because a taxi can provide a travel experience similar to that of the door-to-door service with a private car, and the change in taxi ridership can effectively reflect the transfer inclination. To achieve this work, two penalization-based regression models are employed to implement the variable selection of all built environment factors, since it allows the collinearity between variables. Then, we utilize Wilcoxon’s signed-rank test to scrutinize the regions with significant variation in travel demand.The remainder of this paper is organized as follows. The next section presents a literature review on the LPR policy and the relationship between the built environment and taxi ridership. Section3 introduces the data and variables. Section 4 presents the modeling approach. Section 5 presents the results. The final section highlights the key findings and offers implications for traffic demand management.


## 2. Literature Review

In the past decade, an unforeseen increase in private car ownership and usage [21] has caused severe traffic and environmental problems. The most direct and impactful TDM strategy, the license plate restriction (LPR) policy, has been widely implemented in major Chinese cities. A large number of researches on the effectiveness of that vehicle restriction policy and others have led to positive and negative findings, including urban traffic condition improvements [22–24], air pollution mitigation [22–24], safety efficacies [25], and acceptance and reactions of travelers [15], as well as illegal travel [8, 9, 14, 26].Reducing urban air pollution and the traffic congestion was the original intention of the urban managers and policy-makers who initiated the LPR policy. However, previous works related to the effects of all the vehicle restriction policy show an insignificant impact on air quality [23, 27]. Moreover, the long-run effectiveness of this policy will be weakened by vehicular volume adjustments [28] and the possibility of atmospheric chemistry that can reduce NO but elevate NO2, NOx, and O3 [29]. The direct effect of this policy on traffic conditions is evident to travelers and has been verified by researchers using various quantitative indicators and methods. Several studies have affirmed the LPR policy’s reduction in traffic volume [11–13, 30], which naturally leads to a reduction in traffic congestion [28]. Some studies investigated that the LPR policy led to a rise in travel speed based on the citywide records of the travel speed level [27] and the microwave detectors [13], respectively. Moreover, different levels of effectiveness were also explored between two typical license-plate-number-based vehicle restriction strategies, that is, ODPW and OAE. Mohan et al. [9] found the odd-even day traffic restriction policy forced a traffic demand shift from private vehicles to public transport. Liu et al. [10] compared these two policies from travel demand and traffic condition perspectives based on the massive license plate recognition (LPR) data in Langfang, China. They found the policy shift from ODPW to OAE significantly reduced the traffic volume and improved the travel speed in most of the road segments, in both peak hours and off-peak hours. But the effect was not as significant as expected due to the increase of illegal travel [14, 26] and antipolicy measures such as purchasing a second car [18], nor did the OAE substantially improve the control of traffic demand [31]. Therefore, Jia et al. [16] indicated that promoting the public’s acceptance is a vital prerequisite for implementing the LPR policy. Only the traffic demand shifts from private vehicles to public transport ensure the long-term effectiveness of the vehicle restriction policy [9, 12]. Thus, quantifying the demand change and understanding its determinants are important to realize the mode shift triggered by the vehicle restriction policy.Unfortunately, studies referring to exploring the spatial characteristic of the mode shift and its relationship with the built environment under the LPR policy are still lacking. Enormous literature has confirmed the critical role of the built environment on transit ridership [32, 33]. However, the issue of how restricted travel demands transfer among the different modes under different built environments has been largely neglected in the past years. This is essential for the manager and policy-makers to underpin or support complementary policies used to reduce residents’ beneficial losses (giving up unmanaged car use for a healthier environment) to the utmost extent. Moreover, Çetin and Yasin Er [34] found the transferred demand relies only on the taxi, a complement for other public transport modes. However, few researchers have attempted to investigate the potential problems of public transportation from the taxi standpoint. Therefore, it is necessary to explicitly reexamine the influence of the current built environment on the mode shift under the vehicle restriction policy from the perspective of taxi ridership.


## 3. Study Area and Data

### 3.1. Study Area

Serious fog and haze events have frequently occurred during the autumn and winter months in Xi’an, China [35]. Residents have been suffering a high level of atmospheric particulate matter (PM2.5 and PM10) and various health problems [36, 37]. Therefore, the ODPW policy was first initiated for this designated region on November 4 in 2016 by the Xi’an Environmental Protection Bureau and Municipal Transportation Bureau to reduce vehicle emissions. The ODPW policy was in effect during weekdays from 7 : 00 a.m. to 8 : 00 p.m. except on weekends. Specifically, this strategy restricted vehicles based on the two last digits of the license plates each day, and all new energy vehicles and public transport vehicles, such as taxi and bus, were not affected by this policy. Since the original intention of this policy is alleviating air pollution, the implementation program has not been routinized and depended on the air pollution degree during November. The detailed restricted days and the corresponding last digits are listed in Table 1.Table 1
The ODPW policy implementation in November 2016 in Xi’an.
Mon.Tue.Wed.Thu.Fri.Sat.Sun.123456———(4, 9)——78910111213———————14151617181920—(0, 5)(1, 6)(2, 7)(3, 8)——21222324252627(1, 6)——————282930—(4, 9)(0, 5)Note. The symbol “—” represents no restricted day. The bold values denote the date of license plate restriction and its corresponding last digits.We chose the main urban area of Xi’an within the scope of the LPR policy implementation as the study area of this research, as shown in Figure1. We adopted the officially defined traffic analysis zones (TAZs) as the spatial units to examine the built-environment effects on the spatial change of the taxi demand. TAZ is a basic geographic unit adaptive to transportation analysis; thus, the study area included 698 TAZs, and the dependent and explanatory variables were calculated for each TAZ.Figure 1
Vehicle restriction areas and study areas.The primary purpose of this paper is to explain the spatial change of travel demand led by the vehicle restriction policy. We chose two consecutive restricted and unrestricted weeks and used data from 10 weekdays to eliminate the interference of other factors in the results, such as changes in the regional population and car ownership. Sociodemographic characteristics can be considered static in a short time. Notably, another important factor that can greatly affect the results is the trip purpose. Some elastic demands, such as shopping, entertainment, and tourism, can be possibly changed or canceled once the private vehicles were forced off the road. So this study focused on the travel demand during the morning peak, i.e., from 7 : 00 to 9 : 00 a.m., as the travel plan of the commuting demand cannot be adjusted arbitrarily and has to choose an alternative traffic mode.


### 3.2. Data

The GPS-equipped taxi that is floating throughout the city is a suitable proxy of the detector used to investigate the actual effects of the LPR policy on the demand change. The taxi GPS data recorded the trajectories of about 11,555 urban taxicabs running across the city of Xi’an, which created about 21 million log records. Here the fleet size of the taxi did not change during the observation period. Each GPS record contains information of taxi ID, epoch time, position, velocity, and status (0 for vacancy and 1 for occupation) with a 30-second sampling interval. From the status, one can easily know if the taxi is occupied by the passenger or not at each trajectory sampling point.The number of taxi trips was observed to be fairly consistent every week except for external intervention [38] and exhibited reproducible patterns within weeks [39]. To depict the change induced by the external factors, we compared the daily ridership of two unrestricted weeks (blue lines) and a restricted week (orange line), as shown in Figure 2. Here, for retrieving the ridership data from GPS records, one can refer to [38, 40, 41]. We can intuitively see the evident increase resulting from the vehicle restriction policy. This illustrates that, on the whole, taxi ridership during the morning rush hour was relatively stable within a certain length of the observation period, and the vehicle restriction policy forced some of the travel demand to transfer to the taxi.Figure 2
Daily taxi ridership in 10 consecutive weekdays.To investigate the variation in effects of LPR with time, we illustrated the taxi ridership at the morning peak for these consecutive three weeks, as shown in Figure3. The ridership did not show a significant change on the first day of the ban, because most residents knew nothing or adopt a wait-and-see attitude towards this new policy. Therefore, the data collected on this first day could not be used in the final analysis.Figure 3
Time variations in taxi ridership under LPR policy.


#### 3.2.1. Relative Variation in Daily Taxi Ridership: Dependent Variable

We further investigate the local variation of taxi ridership in a fine-grained spatial perspective. The relative variation (RV) in ridership is used to measure the change of travel demand after the vehicle restriction policy and can eliminate the magnitude effects. The relative variation in each TAZ is calculated as shown in equation (1). Based on the ArcGIS platform, we visualized the different relative variations in TAZs adopting different colors, as shown in Figure 4. Blue stands for the decreased area and red for the increased area. The darker the color, the greater the change:(1)RVi=v¯after,i−v¯before,iv¯before,i,where RVi denotes the relative variation in taxi ridership in the ith TAZ. v¯after,i denotes the average taxi ridership after the vehicles restriction policy and v¯before,i denotes the average taxi ridership before the policy implementation in the ith TAZ.Figure 4
The relative variation of taxi ridership in TAZs.Figure4 shows the vehicle restriction policy has different effects in different regions. Even in two adjacent TAZs, the directions of this influence, i.e., the signs of the relative variation, may be completely opposite. Another result is drawn from the spatial distribution perspective; that is, the number of TAZs with the increased ridership is more than that with the decreased ridership in the southwest. Most TAZs have a decreasing tendency in the north. This complex impact of the vehicle restriction policy on the spatial variation of ridership is mainly caused by the built environment. Consequently, we attempt to explain these different effects from the built environment perspective.


#### 3.2.2. Built Environment: Explanatory Variables

We employed ArcGIS to develop built environment measures for each TAZ, leveraging the Baidu Map Service points of interest (POI) data, road network, and public transport network data. The built environment in this study is fully portrayed by 14 variables from location, land mix (diversity), distance to public transportation, and street design dimensions [42]. Notably, the density variables of population and employment are also very important for explaining the demand change, but these variables are not introduced here due to the inaccessibility of the most recent demographics. This lack of more recent data is due to the lack of a current citywide census, which is conducted every 10 years in most cities in China. The description and descriptive statistics of the 14 attribute variables used in this study are presented in Table 2. The location characteristic is represented by the distance from the centroid of the TAZ to the most important urban landmarks. Distance to the city center (Bell Tower) and distance to the nearest local central business districts (CBDs) are always key factors in assessing commuting behavior [27, 43, 44]. These two indicators were estimated using the Euclidian distance. Since the estimation accuracy of land use using POIs is higher than that using parcel-level land-use maps [45], we employed the POI-based mixed-use indices based on the Hill numbers [46–49] framework to characterize mixed-land use, following Yue et al.’s POI-based mixed-use strategy [50]. Hill numbers are commonly used to quantify species diversity in ecology and biogeography. They can achieve the measurements of multifaceted diversity of land use by calculating(2)Dq=∑i=1spiq1/1−q,where D is the diversity of land use and q is the order of diversity. s is the number of POIs categories present in a TAZ. pi represents the proportion of the ith categories to the total number of POIs within a TAZ. When q=0,1,and2, the equation changes to three different ways of measure, namely, richness, Shannon entropy, and the Simpson index.Table 2
Variable description and descriptive statistics of built environment.
DimensionsVariablesVariable descriptionMeanMax.Min.SDVIFLocationDistance to the city center (D_lc)The Euclidian distance from the centroid of each TAZ to Bell Tower (km)6.50114.6330.2593.3252.835Distance to the nearest local center (D_nlc)The Euclidian distance from the centroid of each TAZ to the nearest regional CBD (km)2.3157.7340.1271.3661.499Land mixRichness (R)The number of different land-use categories present in a TAZ14.56719.0000.0003.8286.674Shannon entropy (Sh_en)“It reflects the amount of order in both POIs categories and the number of POIs” (Yue et al., 2016. pp.9)8.56912.7740.0002.22646.705Simpson index (Si)It presents the evenness or concentration of POIs in a TAZ6.54910.6560.0001.94428.560Public transportationBus accessibility (B_acc)The route distance from the centroid of the TAZ to the nearest bus stop (km)0.3471.8360.0100.2363.676Transit accessibility (Tr_acc)The route distance from the centroid of the TAZ to the nearest metro station (km)1.5246.5320.1291.1292.734Coverage rate of bus stop within a radius of 300 m (cr_b_300)The ratio of bus stop service area within a radius of 300 m to the total area of the TAZ0.6361.0000.0000.2943.811Coverage rate of bus stop within a radius of 500 m (cr_b_500)The ratio of bus stop service area within a radius of 500 m to the total area of the TAZ0.8701.0000.0000.2535.750Coverage rate of metro station within a radius of 800 m (cr_m_800)The ratio of metro station service area within a radius of 800 m to the total area of the TAZ0.3181.0000.0000.4061.988Number of bus routes (N_br)The number of bus routes that pass through the TAZ9.4791.000.0013.111.227Street designDistance to the nearest minor road (D_nmir)Euclidian distance from the centroid of each TAZ to the nearest minor road (km)0.4142.5300.0000.3831.240Distance to the nearest main road (D_nmar)Euclidian distance from the centroid of each TAZ to the nearest main road (km)0.3722.9790.0000.3361.512Distance to the nearest expressway (D_nexpw)Euclidian distance from the centroid of each TAZ to the nearest expressway (km)1.1414.1000.0080.9271.474We assessed the public transport environment by distance-based accessibility, service coverage rate, and the number of stations. The index of distance to transit is described by the route distance from the centroid of the TAZ to the nearest bus stop or metro station [42]. We set up buffer zones with different radii for the bus stop and metro station and obtained the service area where the overlap part between the TAZ and the buffer zones is. Then, the service coverage rate within different radii was computed by the ratio of the service area to the total area of the TAZ. Furthermore, distance to expressways, distance to main roads, and distance to minor roads are computed in this study to measure the degree of convenience to get a taxi from the street environment dimension. It should be noted that although the density of road [27, 38, 51, 52] is commonly used to measure the street environment, this index cannot accurately reflect the opportunity to get a taxi in the TAZ-based spatial unit. We found that some TAZs had a road density of zero in our trial computation, even though adjacent roads to these TAZs were available. This was possibly due to the roads always being processed as TAZ boundaries.We delineated the built environment features as comprehensively as possible. However, the overlap of some dimensions and multiple measurements of the same attribute of the built environment [42, 53] created inevitable multicollinearity problems among variables. Therefore, a variance inflation factor (VIF) procedure was tested to diagnose the potential multicollinearity among the regressors, as shown in Table 3. Pearson’s correlation coefficients were also used to reflect the linear association of two variables, as illustrated in Figure 5. We can see that the most serious collinearity problem exists between indicators related to land mix. Meanwhile, the correlation between these indicators is also high.Table 3
Regression results of LASSO and Elastic Net models.
DimensionsVariablesLASSOElastic NetLocationRMSE0.2402480.239898Constant term0.2417690.259710Number of nonzero parameters68Distance to the city center00Distance to the nearest local center−0.016721−0.043576Land mixRichness−0.091567−0.081240Shannon entropy00Simpson index00Public transportationBus accessibility00Transit accessibility00Coverage rate of bus stop within a radius of 300 m−0.035096−0.034070Coverage rate of bus stop within a radius of 500 m0−0.013624Coverage rate of metro station within a radius of 800 m−0.030505−0.031962Number of bus routes0−0.000423Road designDistance to the nearest minor road00Distance to the nearest main road0.0223410.047424Distance to the nearest expressway0.1381790.121271Figure 5
Correlation analysis of variables based on Pearson’s correlation coefficient.


## 3.1. Study Area

Serious fog and haze events have frequently occurred during the autumn and winter months in Xi’an, China [35]. Residents have been suffering a high level of atmospheric particulate matter (PM2.5 and PM10) and various health problems [36, 37]. Therefore, the ODPW policy was first initiated for this designated region on November 4 in 2016 by the Xi’an Environmental Protection Bureau and Municipal Transportation Bureau to reduce vehicle emissions. The ODPW policy was in effect during weekdays from 7 : 00 a.m. to 8 : 00 p.m. except on weekends. Specifically, this strategy restricted vehicles based on the two last digits of the license plates each day, and all new energy vehicles and public transport vehicles, such as taxi and bus, were not affected by this policy. Since the original intention of this policy is alleviating air pollution, the implementation program has not been routinized and depended on the air pollution degree during November. The detailed restricted days and the corresponding last digits are listed in Table 1.Table 1
The ODPW policy implementation in November 2016 in Xi’an.
Mon.Tue.Wed.Thu.Fri.Sat.Sun.123456———(4, 9)——78910111213———————14151617181920—(0, 5)(1, 6)(2, 7)(3, 8)——21222324252627(1, 6)——————282930—(4, 9)(0, 5)Note. The symbol “—” represents no restricted day. The bold values denote the date of license plate restriction and its corresponding last digits.We chose the main urban area of Xi’an within the scope of the LPR policy implementation as the study area of this research, as shown in Figure1. We adopted the officially defined traffic analysis zones (TAZs) as the spatial units to examine the built-environment effects on the spatial change of the taxi demand. TAZ is a basic geographic unit adaptive to transportation analysis; thus, the study area included 698 TAZs, and the dependent and explanatory variables were calculated for each TAZ.Figure 1
Vehicle restriction areas and study areas.The primary purpose of this paper is to explain the spatial change of travel demand led by the vehicle restriction policy. We chose two consecutive restricted and unrestricted weeks and used data from 10 weekdays to eliminate the interference of other factors in the results, such as changes in the regional population and car ownership. Sociodemographic characteristics can be considered static in a short time. Notably, another important factor that can greatly affect the results is the trip purpose. Some elastic demands, such as shopping, entertainment, and tourism, can be possibly changed or canceled once the private vehicles were forced off the road. So this study focused on the travel demand during the morning peak, i.e., from 7 : 00 to 9 : 00 a.m., as the travel plan of the commuting demand cannot be adjusted arbitrarily and has to choose an alternative traffic mode.


## 3.2. Data

The GPS-equipped taxi that is floating throughout the city is a suitable proxy of the detector used to investigate the actual effects of the LPR policy on the demand change. The taxi GPS data recorded the trajectories of about 11,555 urban taxicabs running across the city of Xi’an, which created about 21 million log records. Here the fleet size of the taxi did not change during the observation period. Each GPS record contains information of taxi ID, epoch time, position, velocity, and status (0 for vacancy and 1 for occupation) with a 30-second sampling interval. From the status, one can easily know if the taxi is occupied by the passenger or not at each trajectory sampling point.The number of taxi trips was observed to be fairly consistent every week except for external intervention [38] and exhibited reproducible patterns within weeks [39]. To depict the change induced by the external factors, we compared the daily ridership of two unrestricted weeks (blue lines) and a restricted week (orange line), as shown in Figure 2. Here, for retrieving the ridership data from GPS records, one can refer to [38, 40, 41]. We can intuitively see the evident increase resulting from the vehicle restriction policy. This illustrates that, on the whole, taxi ridership during the morning rush hour was relatively stable within a certain length of the observation period, and the vehicle restriction policy forced some of the travel demand to transfer to the taxi.Figure 2
Daily taxi ridership in 10 consecutive weekdays.To investigate the variation in effects of LPR with time, we illustrated the taxi ridership at the morning peak for these consecutive three weeks, as shown in Figure3. The ridership did not show a significant change on the first day of the ban, because most residents knew nothing or adopt a wait-and-see attitude towards this new policy. Therefore, the data collected on this first day could not be used in the final analysis.Figure 3
Time variations in taxi ridership under LPR policy.


### 3.2.1. Relative Variation in Daily Taxi Ridership: Dependent Variable

We further investigate the local variation of taxi ridership in a fine-grained spatial perspective. The relative variation (RV) in ridership is used to measure the change of travel demand after the vehicle restriction policy and can eliminate the magnitude effects. The relative variation in each TAZ is calculated as shown in equation (1). Based on the ArcGIS platform, we visualized the different relative variations in TAZs adopting different colors, as shown in Figure 4. Blue stands for the decreased area and red for the increased area. The darker the color, the greater the change:(1)RVi=v¯after,i−v¯before,iv¯before,i,where RVi denotes the relative variation in taxi ridership in the ith TAZ. v¯after,i denotes the average taxi ridership after the vehicles restriction policy and v¯before,i denotes the average taxi ridership before the policy implementation in the ith TAZ.Figure 4
The relative variation of taxi ridership in TAZs.Figure4 shows the vehicle restriction policy has different effects in different regions. Even in two adjacent TAZs, the directions of this influence, i.e., the signs of the relative variation, may be completely opposite. Another result is drawn from the spatial distribution perspective; that is, the number of TAZs with the increased ridership is more than that with the decreased ridership in the southwest. Most TAZs have a decreasing tendency in the north. This complex impact of the vehicle restriction policy on the spatial variation of ridership is mainly caused by the built environment. Consequently, we attempt to explain these different effects from the built environment perspective.


### 3.2.2. Built Environment: Explanatory Variables

We employed ArcGIS to develop built environment measures for each TAZ, leveraging the Baidu Map Service points of interest (POI) data, road network, and public transport network data. The built environment in this study is fully portrayed by 14 variables from location, land mix (diversity), distance to public transportation, and street design dimensions [42]. Notably, the density variables of population and employment are also very important for explaining the demand change, but these variables are not introduced here due to the inaccessibility of the most recent demographics. This lack of more recent data is due to the lack of a current citywide census, which is conducted every 10 years in most cities in China. The description and descriptive statistics of the 14 attribute variables used in this study are presented in Table 2. The location characteristic is represented by the distance from the centroid of the TAZ to the most important urban landmarks. Distance to the city center (Bell Tower) and distance to the nearest local central business districts (CBDs) are always key factors in assessing commuting behavior [27, 43, 44]. These two indicators were estimated using the Euclidian distance. Since the estimation accuracy of land use using POIs is higher than that using parcel-level land-use maps [45], we employed the POI-based mixed-use indices based on the Hill numbers [46–49] framework to characterize mixed-land use, following Yue et al.’s POI-based mixed-use strategy [50]. Hill numbers are commonly used to quantify species diversity in ecology and biogeography. They can achieve the measurements of multifaceted diversity of land use by calculating(2)Dq=∑i=1spiq1/1−q,where D is the diversity of land use and q is the order of diversity. s is the number of POIs categories present in a TAZ. pi represents the proportion of the ith categories to the total number of POIs within a TAZ. When q=0,1,and2, the equation changes to three different ways of measure, namely, richness, Shannon entropy, and the Simpson index.Table 2
Variable description and descriptive statistics of built environment.
DimensionsVariablesVariable descriptionMeanMax.Min.SDVIFLocationDistance to the city center (D_lc)The Euclidian distance from the centroid of each TAZ to Bell Tower (km)6.50114.6330.2593.3252.835Distance to the nearest local center (D_nlc)The Euclidian distance from the centroid of each TAZ to the nearest regional CBD (km)2.3157.7340.1271.3661.499Land mixRichness (R)The number of different land-use categories present in a TAZ14.56719.0000.0003.8286.674Shannon entropy (Sh_en)“It reflects the amount of order in both POIs categories and the number of POIs” (Yue et al., 2016. pp.9)8.56912.7740.0002.22646.705Simpson index (Si)It presents the evenness or concentration of POIs in a TAZ6.54910.6560.0001.94428.560Public transportationBus accessibility (B_acc)The route distance from the centroid of the TAZ to the nearest bus stop (km)0.3471.8360.0100.2363.676Transit accessibility (Tr_acc)The route distance from the centroid of the TAZ to the nearest metro station (km)1.5246.5320.1291.1292.734Coverage rate of bus stop within a radius of 300 m (cr_b_300)The ratio of bus stop service area within a radius of 300 m to the total area of the TAZ0.6361.0000.0000.2943.811Coverage rate of bus stop within a radius of 500 m (cr_b_500)The ratio of bus stop service area within a radius of 500 m to the total area of the TAZ0.8701.0000.0000.2535.750Coverage rate of metro station within a radius of 800 m (cr_m_800)The ratio of metro station service area within a radius of 800 m to the total area of the TAZ0.3181.0000.0000.4061.988Number of bus routes (N_br)The number of bus routes that pass through the TAZ9.4791.000.0013.111.227Street designDistance to the nearest minor road (D_nmir)Euclidian distance from the centroid of each TAZ to the nearest minor road (km)0.4142.5300.0000.3831.240Distance to the nearest main road (D_nmar)Euclidian distance from the centroid of each TAZ to the nearest main road (km)0.3722.9790.0000.3361.512Distance to the nearest expressway (D_nexpw)Euclidian distance from the centroid of each TAZ to the nearest expressway (km)1.1414.1000.0080.9271.474We assessed the public transport environment by distance-based accessibility, service coverage rate, and the number of stations. The index of distance to transit is described by the route distance from the centroid of the TAZ to the nearest bus stop or metro station [42]. We set up buffer zones with different radii for the bus stop and metro station and obtained the service area where the overlap part between the TAZ and the buffer zones is. Then, the service coverage rate within different radii was computed by the ratio of the service area to the total area of the TAZ. Furthermore, distance to expressways, distance to main roads, and distance to minor roads are computed in this study to measure the degree of convenience to get a taxi from the street environment dimension. It should be noted that although the density of road [27, 38, 51, 52] is commonly used to measure the street environment, this index cannot accurately reflect the opportunity to get a taxi in the TAZ-based spatial unit. We found that some TAZs had a road density of zero in our trial computation, even though adjacent roads to these TAZs were available. This was possibly due to the roads always being processed as TAZ boundaries.We delineated the built environment features as comprehensively as possible. However, the overlap of some dimensions and multiple measurements of the same attribute of the built environment [42, 53] created inevitable multicollinearity problems among variables. Therefore, a variance inflation factor (VIF) procedure was tested to diagnose the potential multicollinearity among the regressors, as shown in Table 3. Pearson’s correlation coefficients were also used to reflect the linear association of two variables, as illustrated in Figure 5. We can see that the most serious collinearity problem exists between indicators related to land mix. Meanwhile, the correlation between these indicators is also high.Table 3
Regression results of LASSO and Elastic Net models.
DimensionsVariablesLASSOElastic NetLocationRMSE0.2402480.239898Constant term0.2417690.259710Number of nonzero parameters68Distance to the city center00Distance to the nearest local center−0.016721−0.043576Land mixRichness−0.091567−0.081240Shannon entropy00Simpson index00Public transportationBus accessibility00Transit accessibility00Coverage rate of bus stop within a radius of 300 m−0.035096−0.034070Coverage rate of bus stop within a radius of 500 m0−0.013624Coverage rate of metro station within a radius of 800 m−0.030505−0.031962Number of bus routes0−0.000423Road designDistance to the nearest minor road00Distance to the nearest main road0.0223410.047424Distance to the nearest expressway0.1381790.121271Figure 5
Correlation analysis of variables based on Pearson’s correlation coefficient.


## 3.2.1. Relative Variation in Daily Taxi Ridership: Dependent Variable

We further investigate the local variation of taxi ridership in a fine-grained spatial perspective. The relative variation (RV) in ridership is used to measure the change of travel demand after the vehicle restriction policy and can eliminate the magnitude effects. The relative variation in each TAZ is calculated as shown in equation (1). Based on the ArcGIS platform, we visualized the different relative variations in TAZs adopting different colors, as shown in Figure 4. Blue stands for the decreased area and red for the increased area. The darker the color, the greater the change:(1)RVi=v¯after,i−v¯before,iv¯before,i,where RVi denotes the relative variation in taxi ridership in the ith TAZ. v¯after,i denotes the average taxi ridership after the vehicles restriction policy and v¯before,i denotes the average taxi ridership before the policy implementation in the ith TAZ.Figure 4
The relative variation of taxi ridership in TAZs.Figure4 shows the vehicle restriction policy has different effects in different regions. Even in two adjacent TAZs, the directions of this influence, i.e., the signs of the relative variation, may be completely opposite. Another result is drawn from the spatial distribution perspective; that is, the number of TAZs with the increased ridership is more than that with the decreased ridership in the southwest. Most TAZs have a decreasing tendency in the north. This complex impact of the vehicle restriction policy on the spatial variation of ridership is mainly caused by the built environment. Consequently, we attempt to explain these different effects from the built environment perspective.


## 3.2.2. Built Environment: Explanatory Variables

We employed ArcGIS to develop built environment measures for each TAZ, leveraging the Baidu Map Service points of interest (POI) data, road network, and public transport network data. The built environment in this study is fully portrayed by 14 variables from location, land mix (diversity), distance to public transportation, and street design dimensions [42]. Notably, the density variables of population and employment are also very important for explaining the demand change, but these variables are not introduced here due to the inaccessibility of the most recent demographics. This lack of more recent data is due to the lack of a current citywide census, which is conducted every 10 years in most cities in China. The description and descriptive statistics of the 14 attribute variables used in this study are presented in Table 2. The location characteristic is represented by the distance from the centroid of the TAZ to the most important urban landmarks. Distance to the city center (Bell Tower) and distance to the nearest local central business districts (CBDs) are always key factors in assessing commuting behavior [27, 43, 44]. These two indicators were estimated using the Euclidian distance. Since the estimation accuracy of land use using POIs is higher than that using parcel-level land-use maps [45], we employed the POI-based mixed-use indices based on the Hill numbers [46–49] framework to characterize mixed-land use, following Yue et al.’s POI-based mixed-use strategy [50]. Hill numbers are commonly used to quantify species diversity in ecology and biogeography. They can achieve the measurements of multifaceted diversity of land use by calculating(2)Dq=∑i=1spiq1/1−q,where D is the diversity of land use and q is the order of diversity. s is the number of POIs categories present in a TAZ. pi represents the proportion of the ith categories to the total number of POIs within a TAZ. When q=0,1,and2, the equation changes to three different ways of measure, namely, richness, Shannon entropy, and the Simpson index.Table 2
Variable description and descriptive statistics of built environment.
DimensionsVariablesVariable descriptionMeanMax.Min.SDVIFLocationDistance to the city center (D_lc)The Euclidian distance from the centroid of each TAZ to Bell Tower (km)6.50114.6330.2593.3252.835Distance to the nearest local center (D_nlc)The Euclidian distance from the centroid of each TAZ to the nearest regional CBD (km)2.3157.7340.1271.3661.499Land mixRichness (R)The number of different land-use categories present in a TAZ14.56719.0000.0003.8286.674Shannon entropy (Sh_en)“It reflects the amount of order in both POIs categories and the number of POIs” (Yue et al., 2016. pp.9)8.56912.7740.0002.22646.705Simpson index (Si)It presents the evenness or concentration of POIs in a TAZ6.54910.6560.0001.94428.560Public transportationBus accessibility (B_acc)The route distance from the centroid of the TAZ to the nearest bus stop (km)0.3471.8360.0100.2363.676Transit accessibility (Tr_acc)The route distance from the centroid of the TAZ to the nearest metro station (km)1.5246.5320.1291.1292.734Coverage rate of bus stop within a radius of 300 m (cr_b_300)The ratio of bus stop service area within a radius of 300 m to the total area of the TAZ0.6361.0000.0000.2943.811Coverage rate of bus stop within a radius of 500 m (cr_b_500)The ratio of bus stop service area within a radius of 500 m to the total area of the TAZ0.8701.0000.0000.2535.750Coverage rate of metro station within a radius of 800 m (cr_m_800)The ratio of metro station service area within a radius of 800 m to the total area of the TAZ0.3181.0000.0000.4061.988Number of bus routes (N_br)The number of bus routes that pass through the TAZ9.4791.000.0013.111.227Street designDistance to the nearest minor road (D_nmir)Euclidian distance from the centroid of each TAZ to the nearest minor road (km)0.4142.5300.0000.3831.240Distance to the nearest main road (D_nmar)Euclidian distance from the centroid of each TAZ to the nearest main road (km)0.3722.9790.0000.3361.512Distance to the nearest expressway (D_nexpw)Euclidian distance from the centroid of each TAZ to the nearest expressway (km)1.1414.1000.0080.9271.474We assessed the public transport environment by distance-based accessibility, service coverage rate, and the number of stations. The index of distance to transit is described by the route distance from the centroid of the TAZ to the nearest bus stop or metro station [42]. We set up buffer zones with different radii for the bus stop and metro station and obtained the service area where the overlap part between the TAZ and the buffer zones is. Then, the service coverage rate within different radii was computed by the ratio of the service area to the total area of the TAZ. Furthermore, distance to expressways, distance to main roads, and distance to minor roads are computed in this study to measure the degree of convenience to get a taxi from the street environment dimension. It should be noted that although the density of road [27, 38, 51, 52] is commonly used to measure the street environment, this index cannot accurately reflect the opportunity to get a taxi in the TAZ-based spatial unit. We found that some TAZs had a road density of zero in our trial computation, even though adjacent roads to these TAZs were available. This was possibly due to the roads always being processed as TAZ boundaries.We delineated the built environment features as comprehensively as possible. However, the overlap of some dimensions and multiple measurements of the same attribute of the built environment [42, 53] created inevitable multicollinearity problems among variables. Therefore, a variance inflation factor (VIF) procedure was tested to diagnose the potential multicollinearity among the regressors, as shown in Table 3. Pearson’s correlation coefficients were also used to reflect the linear association of two variables, as illustrated in Figure 5. We can see that the most serious collinearity problem exists between indicators related to land mix. Meanwhile, the correlation between these indicators is also high.Table 3
Regression results of LASSO and Elastic Net models.
DimensionsVariablesLASSOElastic NetLocationRMSE0.2402480.239898Constant term0.2417690.259710Number of nonzero parameters68Distance to the city center00Distance to the nearest local center−0.016721−0.043576Land mixRichness−0.091567−0.081240Shannon entropy00Simpson index00Public transportationBus accessibility00Transit accessibility00Coverage rate of bus stop within a radius of 300 m−0.035096−0.034070Coverage rate of bus stop within a radius of 500 m0−0.013624Coverage rate of metro station within a radius of 800 m−0.030505−0.031962Number of bus routes0−0.000423Road designDistance to the nearest minor road00Distance to the nearest main road0.0223410.047424Distance to the nearest expressway0.1381790.121271Figure 5
Correlation analysis of variables based on Pearson’s correlation coefficient.


## 4. Method

Given the multicollinearity and high-dimension problems, we utilized two penalization-based regression models [54], the LASSO and the Elastic Net (EN), to estimate the coefficients of the built environment. Compared to traditional linear regression models, the penalization-based approach has two outstanding advantages: the interpretability of the selected variables and the stability of a rigid regression estimation. Moreover, we also explored the difference in results between these two models, which is caused by different types of regularization in the penalty terms. At last, the root mean squared error (RMSE) of both models was reported to evaluate their performance.The LASSO (least absolute shrinkage and selection operator) approach presented by Tibshirani (1996) [55] can select the variable and infer the causal relationship simultaneously through shrinking some trivial estimated coefficients equal to zero. The basic framework of the variable selection minimizes the residual sum of squares subject to the sum of the coefficients’ absolute value being less than a constant. The LASSO estimator β^ [54] can be represented by the following Lagrangian form:(3)Lβ,λ=minβ∈ℜp1Ny−XβT22+λβ1,where equation (4) is the L1 norm penalty which is not differentiable at zero. The tuning parameter λ≥0, also known as the punishment coefficient, is used to balance the empirical error and the sparsity of the model parameter. When the parameter λ=0, the method turns to the general linear regression model:(4)β1=∑j=1pβj.Elastic Net [56] is an extension of LASSO by adding a quadratic part to the penalty, which uses a combination of LASSO and Ridge regression methods. The objective function of Elastic Net is defined by(5)Lβ,λ1,λ2=minβ∈ℜp1Ny−XβT22+λ1β1+λ2β22,where equation (6) is the L2 norm penalty. Elastic Net (EN) is a hybrid model and works between LASSO and Ridge regression by setting two tuning parameters λ1,λ2≥0. Let α=λ1+λ2 and L1_ratio=λ1/λ1+λ2. L1_ratio can be determined by minimizing the RMSE, as shown in Figure 6:(6)β22=∑j=1pβj2.Figure 6
Variation of RMSE withL1_ratio.All feature variable values were preprocessed by using theZ-score standardization and centered methods. The regularization parameters of these two regression models were estimated by performing the 10-fold cross validation (CV) [57] due to its outperformance of AIC, BIC, and 5-fold CV in our study. The nonzero regression coefficients and their corresponding feature variables in LASSO and Elastic Net are listed in Table 3.


## 5. Discussion

### 5.1. Regression Coefficients

It can be seen from Table3 that different numbers of features were selected by different models. The feature variables corresponding to the nonzero regression coefficients in the Elastic Net model are two more than those in the LASSO model, which confirmed the sparsity of LASSO results caused by the L1 norm penalty [57]. Figure 7 shows the different regression paths of coefficients under different penalty terms. The horizontal axis represents tuning parameters and the vertical axis represents the coefficient of variation.Figure 7
Coefficient paths of LASSO and Elastic Net regression models (letα=λ in LASSO and α=λ1+λ2 in Elastic Net).As shown in Table3, the RMSE of the Elastic Net is slightly lower than that of LASSO, which means the Elastic Net model performed better in this study. Therefore, further analysis will be based on the results of the Elastic Net. Compared to the Elastic Net, the regression coefficients of coverage rate of bus stop within a radius of 500 m and number of bus routes were shrunk to zero by LASSO. Meanwhile the influence degree and direction of other features estimated by Elastic Net were consistent with the results of LASSO. This is due to the same mechanism of parameter estimation between these two regression models. Overall, four dimensions of the built environment had an influence on the relative variation (RV) of taxi ridership in varying degrees. Eight out of 14 variables appear to have significant influences on relative variations (RVs). Some variables of location, land mix, and public transportation dimensions have negative effects on the RV; they are distance to the nearest local center, richness, coverage rate of bus stop within a radius of 300 m, coverage rate of bus stop within a radius of 500 m, coverage rate of the metro station within a radius of 800 m, and number of bus routes. Meanwhile features of distance to the nearest main road and distance to the nearest expressway in the road design dimension have positive influences on the RV.Specifically, for variables pertaining to public transportation, the RV is primarily influenced by thecoverage rate of the bus stop within a radius of 300 m, secondarily by the coverage rate of the metro station within a radius of 800 m, and thirdly by the coverage rate of the bus stop within a radius of 500 m. These variables are all inversely related to RV, which means that relative variations in taxi ridership are generally greater in locations that have inconvenient public transportation. This is reasonable because more people have to choose a taxi when their private cars are restricted on road in the environment of inconvenient public transportation. It should be noted that the rail transit network in Xi’an has not formed yet and only four subway lines are in operation. Therefore, the degree of urban rail transit influence will be far greater than the results of current analysis due to its overwhelming advantage in traffic capacity. Note that two highly correlated variables cr_b_300 and cr_b_500 are permitted into the Elastic Net regression model, but the latter was insignificant in the LASSO model. This makes sense because these two indicators reflect different levels of public transportation service. Regions with a higher level of public transportation service have a stronger ability of vehicle restriction policy response. For the land-use dimension, both two models kept the richness. From the location perspective, distance to the nearest local center has a significant influence effect and distance to city center does not. This suggests that a convenient living environment can mitigate the traffic contradiction caused by the LPR policy. More concretely, to reduce the negative externalities of the LPR policy, improving the service level of the public transport network is the primary mean of transferring the travel demand restricted by the LPR policy to public transit. In linear regression, a higher value of the absolute magnitude indicates a greater marginal effect. Therefore, for the policy-maker and transportation planner, the priority is to increase the basic requirements for the service coverage of public transportation. Then the walking distance to the bus stop is considered to be further shortened by increasing the bus stop density, that is, the coverage rate of the bus stop within the radius of 300 m. Besides, planning more urban subcenters to improve the convenience of the living environment is a long-term solution.The above conclusions are not only consistent with previous findings [9, 10] but improve and broaden the existing literature from the research perspective and data type. Although the vehicle plate recognition data [10] is more accurate and has a higher coverage rate than personal trip surveys [26, 58] in evaluating road network performance, the data collection spots are limited to the intersections equipped with detectors. Compared with this data, the floating car data in this study can be theoretically used to obtain traffic conditions at any time and any location on a road network. Therefore, we conducted a finer-grained spatial analysis in travel demand transfer based on taxi GPS data, which has not yet been reported to the best of our knowledge.


### 5.2. Contribution of Different Variables

To investigate the relative impact of different variables, we analyzed the different contributions of selected built environment variables to the relative variation of taxi ridership. For a given variable, its contributions differ from the TAZ. Hence, we compute the mean contribution and its associated standard deviation for each built environment variable. The results are illustrated with the error bar in Figure8. The circle represents the average contribution level of its corresponding variable and the length of the bar represents its volatility of contribution in different regions. The longer the bar is, the more dispersed contribution of a variable is. We can see that features pertaining to public transportation all have a stable contribution as a whole. This suggested that improving the level of public transportation service will be effective in any region. However, the contrition rate of richness with the longest bar has large differences among regions. Therefore, further research about specific types of land use will be necessary before the related strategies can be proposed.Figure 8
Contribution rate of different built environment variables.


### 5.3. Significantly Influenced Regions

Regions with great variations in taxi ridership are significantly influenced by the vehicle restriction policy and need more effective improvement measures to meet the transferred travel demand. Combined with Wilcoxon’s signed-rank test, a nonparametric counterpart of the pairedt-test, we further scrutinized the regions with statistically significant differences, as shown in Figure 9. We screened out the top ten increased TAZs and the top three (all) decreased TAZs with a statistically significant difference as the priority areas for improvement. The number of these selected 13 TAZs and their corresponding built environment parameters are listed in Table 4. The rank is derived based on the absolute value of the relative variation. Overall, taxi ridership throughout the main city presented a growth trend. This suggests that the vehicle restriction policy imposed tremendous pressure on the taxi industry, which is positioned as an auxiliary public transport mode. The only three decreased areas were those that generally have the high value of selected built environment variables, which means that public transportation handled the transferred travel demand from the individual travel mode.Figure 9
Spatial distribution of relative variation in taxi ridership on Friday.Table 4
Significantly influenced regions ordered by the variation in taxi ridership.
TendencyNo. of TAZsRelative variationImportant indexes of the built environmentValueRankcr_b_300cr_b_500cr_m_800Increase710.978310.51251.00000.00004620.875020.428910.99981110.785730.00000.00000.0000500.665440.03030.18380.00002380.631650.91541.00000.00002050.614560.73051.00000.00005050.613970.57701.00000.0000920.609880.41960.80920.00006640.602990.13770.57770.96544540.555010110Decrease175−0.439210.17040.31870.3119266−0.394620.61760.99910.0000305−0.159630.94421.00000.6630Note. cr_b_300 is the abbreviation for the coverage rate of bus stop within a radius of 300 meters; cr_m_800 is the abbreviation for the coverage rate of the metro station within a radius of 800 meters.


## 5.1. Regression Coefficients

It can be seen from Table3 that different numbers of features were selected by different models. The feature variables corresponding to the nonzero regression coefficients in the Elastic Net model are two more than those in the LASSO model, which confirmed the sparsity of LASSO results caused by the L1 norm penalty [57]. Figure 7 shows the different regression paths of coefficients under different penalty terms. The horizontal axis represents tuning parameters and the vertical axis represents the coefficient of variation.Figure 7
Coefficient paths of LASSO and Elastic Net regression models (letα=λ in LASSO and α=λ1+λ2 in Elastic Net).As shown in Table3, the RMSE of the Elastic Net is slightly lower than that of LASSO, which means the Elastic Net model performed better in this study. Therefore, further analysis will be based on the results of the Elastic Net. Compared to the Elastic Net, the regression coefficients of coverage rate of bus stop within a radius of 500 m and number of bus routes were shrunk to zero by LASSO. Meanwhile the influence degree and direction of other features estimated by Elastic Net were consistent with the results of LASSO. This is due to the same mechanism of parameter estimation between these two regression models. Overall, four dimensions of the built environment had an influence on the relative variation (RV) of taxi ridership in varying degrees. Eight out of 14 variables appear to have significant influences on relative variations (RVs). Some variables of location, land mix, and public transportation dimensions have negative effects on the RV; they are distance to the nearest local center, richness, coverage rate of bus stop within a radius of 300 m, coverage rate of bus stop within a radius of 500 m, coverage rate of the metro station within a radius of 800 m, and number of bus routes. Meanwhile features of distance to the nearest main road and distance to the nearest expressway in the road design dimension have positive influences on the RV.Specifically, for variables pertaining to public transportation, the RV is primarily influenced by thecoverage rate of the bus stop within a radius of 300 m, secondarily by the coverage rate of the metro station within a radius of 800 m, and thirdly by the coverage rate of the bus stop within a radius of 500 m. These variables are all inversely related to RV, which means that relative variations in taxi ridership are generally greater in locations that have inconvenient public transportation. This is reasonable because more people have to choose a taxi when their private cars are restricted on road in the environment of inconvenient public transportation. It should be noted that the rail transit network in Xi’an has not formed yet and only four subway lines are in operation. Therefore, the degree of urban rail transit influence will be far greater than the results of current analysis due to its overwhelming advantage in traffic capacity. Note that two highly correlated variables cr_b_300 and cr_b_500 are permitted into the Elastic Net regression model, but the latter was insignificant in the LASSO model. This makes sense because these two indicators reflect different levels of public transportation service. Regions with a higher level of public transportation service have a stronger ability of vehicle restriction policy response. For the land-use dimension, both two models kept the richness. From the location perspective, distance to the nearest local center has a significant influence effect and distance to city center does not. This suggests that a convenient living environment can mitigate the traffic contradiction caused by the LPR policy. More concretely, to reduce the negative externalities of the LPR policy, improving the service level of the public transport network is the primary mean of transferring the travel demand restricted by the LPR policy to public transit. In linear regression, a higher value of the absolute magnitude indicates a greater marginal effect. Therefore, for the policy-maker and transportation planner, the priority is to increase the basic requirements for the service coverage of public transportation. Then the walking distance to the bus stop is considered to be further shortened by increasing the bus stop density, that is, the coverage rate of the bus stop within the radius of 300 m. Besides, planning more urban subcenters to improve the convenience of the living environment is a long-term solution.The above conclusions are not only consistent with previous findings [9, 10] but improve and broaden the existing literature from the research perspective and data type. Although the vehicle plate recognition data [10] is more accurate and has a higher coverage rate than personal trip surveys [26, 58] in evaluating road network performance, the data collection spots are limited to the intersections equipped with detectors. Compared with this data, the floating car data in this study can be theoretically used to obtain traffic conditions at any time and any location on a road network. Therefore, we conducted a finer-grained spatial analysis in travel demand transfer based on taxi GPS data, which has not yet been reported to the best of our knowledge.


## 5.2. Contribution of Different Variables

To investigate the relative impact of different variables, we analyzed the different contributions of selected built environment variables to the relative variation of taxi ridership. For a given variable, its contributions differ from the TAZ. Hence, we compute the mean contribution and its associated standard deviation for each built environment variable. The results are illustrated with the error bar in Figure8. The circle represents the average contribution level of its corresponding variable and the length of the bar represents its volatility of contribution in different regions. The longer the bar is, the more dispersed contribution of a variable is. We can see that features pertaining to public transportation all have a stable contribution as a whole. This suggested that improving the level of public transportation service will be effective in any region. However, the contrition rate of richness with the longest bar has large differences among regions. Therefore, further research about specific types of land use will be necessary before the related strategies can be proposed.Figure 8
Contribution rate of different built environment variables.


## 5.3. Significantly Influenced Regions

Regions with great variations in taxi ridership are significantly influenced by the vehicle restriction policy and need more effective improvement measures to meet the transferred travel demand. Combined with Wilcoxon’s signed-rank test, a nonparametric counterpart of the pairedt-test, we further scrutinized the regions with statistically significant differences, as shown in Figure 9. We screened out the top ten increased TAZs and the top three (all) decreased TAZs with a statistically significant difference as the priority areas for improvement. The number of these selected 13 TAZs and their corresponding built environment parameters are listed in Table 4. The rank is derived based on the absolute value of the relative variation. Overall, taxi ridership throughout the main city presented a growth trend. This suggests that the vehicle restriction policy imposed tremendous pressure on the taxi industry, which is positioned as an auxiliary public transport mode. The only three decreased areas were those that generally have the high value of selected built environment variables, which means that public transportation handled the transferred travel demand from the individual travel mode.Figure 9
Spatial distribution of relative variation in taxi ridership on Friday.Table 4
Significantly influenced regions ordered by the variation in taxi ridership.
TendencyNo. of TAZsRelative variationImportant indexes of the built environmentValueRankcr_b_300cr_b_500cr_m_800Increase710.978310.51251.00000.00004620.875020.428910.99981110.785730.00000.00000.0000500.665440.03030.18380.00002380.631650.91541.00000.00002050.614560.73051.00000.00005050.613970.57701.00000.0000920.609880.41960.80920.00006640.602990.13770.57770.96544540.555010110Decrease175−0.439210.17040.31870.3119266−0.394620.61760.99910.0000305−0.159630.94421.00000.6630Note. cr_b_300 is the abbreviation for the coverage rate of bus stop within a radius of 300 meters; cr_m_800 is the abbreviation for the coverage rate of the metro station within a radius of 800 meters.


## 6. Conclusions

Given the intrinsic relationship between taxis and public transport, this study stands by the position of the citywide taxi services to rethink the shortcomings of public transportation service under the vehicle restriction policy. The two types of penalization-based regression models, the LASSO and Elastic Net regression models, were employed to explore how the built environment influences the mode shift of commuting travel, using taxi GPS data covering the periods before and after the ODPW policy implementation in Xi’an.The contrastive analysis demonstrates that the vehicle restriction policy indeed triggered an overwhelming amount of restricted private car travel demand transferred to a taxi, which was supposed to be shared by public transportation. The regression analysis shows that the relative variation is most strongly related todistance to the nearest local center and secondarily to coverage rate of the bus stop within a radius of 300 m. These findings highlight the vital role of improving public transport services in implementing vehicle restriction policies, which is consistent with the results of previous studies [9, 12]. Furthermore, the marginal benefit and contribution rate of selected variables can provide the implications to the regional differences in the effectiveness of the strategies.Besides, it is confirmed that ineffective enforcement also contributed to the unsatisfactory results of LPR [18]. Therefore, a set of reasonable penalty mechanisms should be established to strengthen the effectiveness of the vehicle restriction policy. For instance, illegal travel was captured by surveillance and punished by driving offense points and traffic violation fines. However, for travelers with a relatively high-income level, the deterrent effect of this penalty measure will be offset by higher time value [10]. Thus, it is necessary to link to the individual credit system to restrain illegal travel behavior.In the future, a growing number of cities like Xi’an could adopt the LPR policy to rein in automobile use and will have to experience the adaptive phase of implementing this policy. Therefore, the findings of this study could help policy-makers predict the possible influence of the LPR policy on urban wide travel demand and refine the policy of travel demand management. However, there also exist a few limitations in this study. First, the other types of travel demand data, such as AFC data, bus IC card data, and ride-sourcing data [59], should be used for data fusion analysis to understand the mode transfer under the LPR policy more precisely. Second, as previously mentioned, population and employment in the demographics should be introduced into the regression model in further research.



---

*Source: 9848095-2020-07-17.xml*