Can Docked Bike-Sharing Systems Reach Their Dual Sustainability in Terms of Environmental Beneﬁts and Financial Operations? A Comparative Study from Nanjing, 2017 and 2023

: In this paper, we investigate the sustainability of docked bike-sharing in Nanjing in terms of environmental beneﬁts and ﬁnancial operations by comparing the data of March 2017 and March 2023 in Nanjing. We modify a community detection method, give and prove dynamic boundary conditions for the objective function of the heuristic algorithm, and realize the estimation of the rebalancing coefﬁcients for this mega-system, thus obtaining more accurate emission factors. We ﬁnd that there are signiﬁcant differences in the results obtained from environmental beneﬁt assessments over time. Further, there are also signiﬁcant differences at the national level. This may signify that the assessment data of one country’s system cannot give a direct reference for another country’s system. Second, we considered the economic basis required for the environmental beneﬁts of docked bike-sharing systems. We have calculated the sustainability of the system’s ﬁnancial operations by considering its revenues over the next nine years, including the cost of facility inputs, facility upgrades, dispatching costs, labor costs, maintenance costs, and the time value of money. The results show a 4.6-fold difference in emission factors between 2017 and 2023; comparing 2017 to 2023 (when demand loss has been severe), the investment in 2017 will be recouped 2 years later than in 2023. Switching distribution vehicles from fuel vehicles to electric trikes would severely deteriorate the operator’s key ﬁnancial metrics while only reducing the emission factor value by 8.64 gCO 2 eq/km, leading to an unsustainable system. This signals the potential for the ﬁnancial unsustainability, or even bankruptcy, of operators if the requirements for sustained emissions reductions from the bike-sharing system are divorced from the form of the economy on which it is sustainably operated. Finally, we consider the geographical patterns between environmental beneﬁts and ﬁnancial operations. We ﬁnd that ﬁnancial sustainability varies across geographic locations. Under ﬁnancial sustainability, we gave emission factors under the mix distribution vehicle scenario.


Introduction
Bike-sharing systems have grown rapidly in recent years and are considered important in reducing CO 2 emissions [1].Bike-sharing is considered to play a role in encouraging people to participate in energy saving and emission reduction [2].It promotes green transportation and contributes significantly to carbon neutrality [3].One study showed that Shanghai reduced 25,240 tons of CO 2 emissions due to the use of bike-sharing in 2016 alone [4].Docked bike-sharing systems are station-based, and they place bicycles in a network of stations and docks built throughout the city for users to pick up and use [5].Nanjing's docked bike-sharing system was established in 2015 [6], and those using the docked bikesharing system must rent bicycles from designated stations and then return them to docks at other bicycle stations [7].In China, docked bike-sharing systems are built and operated by the government, so in most cities, docked bicycles can be used for free for 30 to 45 min [8].The variables used in this paper and their corresponding interpretations can be viewed in the Abbreviations section.

Research Background
Currently, there are more studies on environmental benefit assessments based on bike-sharing systems.The very first study was based on the USA docked and dockless bike-sharing systems [9][10][11][12].In contrast, the study by Chen et al. on a dockless bike-sharing system in Beijing demonstrated completely different emission ratios than the above studies, especially regarding the rebalancing component [2].This suggests that bike-sharing systems vary from country to country, which can affect the potential environmental benefits, and thus case studies in different contexts are necessary for reliable estimation [13].In addition, there may be significant differences for different periods for the same bike-sharing system.For example, Wang et al. obtained significantly different emission factors by modeling seven different development stages of dockless bike-sharing in Beijing [14].
The entire life cycle of a shared bicycle includes multiple stages [13].At present, there is no standardized methodology for assessing the overall impact of the bike-sharing system [15].However, it mainly includes the production, use, management, and end-of-life of the infrastructure of the bike-sharing system [13], as well as the rebalancing operation and maintenance of the bike-sharing system [12].Similarly, fewer studies include nonsubjective estimates of rebalancing operations in the assessment of the entire life cycle of a shared bicycle [12].In addition, the study that included real rebalancing operations focused on a bike-sharing system located in New York City, USA [12] The results cannot be used directly because the rebalancing method, user habits, and distance distribution of stations are quite different from those of bike-sharing in China.As a result, there are fewer results for evaluating large-scale bike-sharing systems in China that include real rebalancing operations.In this paper, the rebalancing coefficients we obtain are quite different from those given in other studies [2,11].Moreover, there are also large differences between the rebalancing coefficients in different periods.Further, few studies comparatively assess the environmental benefits of the operation of shared bicycles in different periods.Moreover, most of the existing studies focus on the environmental benefits of bike-sharing systems and give policy recommendations based on the assessment results, and very few studies consider the economic basis required for their environmental benefits.If the sustainable emission reduction requirements of the bike-sharing system are detached from the economic form on which its orderly operation is based, leading to the financial unsustainability of the operator, the ultimate result may be the disappearance of the whole system; thus, retracting the emission reduction effect it has already achieved.Therefore, we assess the environmental benefits along with the financial and operational benefits for comparison.

Research Questions
This paper focuses on the following issues.First of all, this paper wants to construct an environmental benefit assessment of a large-scale bike-sharing system containing real rebalancing operations suitable for Nanjing.Second, this paper seeks to investigate whether docked bike-sharing systems can achieve financial operational sustainability while pursuing their environmental benefits, i.e., to consider the economic basis for the environmental benefits of docked bike-sharing systems.Finally, we want to find out if studies on the emission factors of bike-sharing systems differed at the national or city level.

Research Gaps Filled
In the research of this paper, we found that the results obtained from the environmental benefit assessment in different periods are significantly different, which is ignored by previous studies.Further, prior studies have focused only on the environmental benefits of the bike-sharing system, and few studies have considered the economic basis for the environmental benefits.The requirements for sustained emissions reductions from a bike-sharing system that are divorced from the form of the economy on which it is sustainably operated may lead to financial unsustainability or even bankruptcy for the operator.The research in this paper examines the environmental benefits of docked bike-sharing systems against the economic basis on which they are required to operate, giving more comprehensive advice.Based on this, we wanted to understand if there is a geographical pattern to the cut-off between environmental benefits and financial operations.Thus, we discussed it again in different communities.The results show that in different geographical locations, the system has different profit or loss results.More emission reductions can be achieved with financial sustainability by using a mixed distribution vehicle program.Finally, we concluded by comparing existing studies and found significant differences at the national level.This may indicate that the evaluation data of one country's system cannot be directly referenced for another country's system.Based on the above study, we give recommendations for the operation of the docked bike-sharing system in Nanjing.

Contributions of This Paper
The main contributions of this paper are as follows.First, we demonstrate dynamic boundary conditions for the objective function of the heuristic algorithm to estimate the rebalancing coefficients of the mega-system.Second, we calculated the sustainability of the system's financial operations by considering its revenues over the next nine years, including facility input costs, facility renewal costs, dispatch costs, labor costs, maintenance costs, and the time value of money.Third, we considered the impact of switching to a distribution vehicle that is considered more environmentally friendly on its environmental performance as well as its financial position to analyze the feasibility of further emission reductions.We have also explored geographic patterns between environmental benefits and financial operations, proposing mixed distribute vehicle programs that can balance the financial sustainability of the operator with further emission reductions.Finally, we give recommendations based on the differences in the comparisons with other studies as well as the conclusions of this paper.

Literature Review
Bike-sharing programs have positive externalities [16].If 75% of the distance traveled by shared bicycles in 2015 had been caused by the substitution of cars, the CO 2 emissions from road traffic in Beijing would have been reduced by nearly 616,040 tons [16].In an assessment of shared bicycle use in Shanghai, it was obtained that 25,240 tons of CO 2 emissions was reduced in 2016 alone due to shared bicycle use [4].An assessment of bikesharing use in New York City also showed that between 2014 and 2017, the bike-sharing system reduced emissions by a total of 30,070 tons of carbon dioxide [17].
There are two types of bike-sharing systems: dockless bike-sharing systems and docked bike-sharing systems.They both have their user groups, with more than 40% of users only using one or the other [18].As a member of public transportation, studies have pointed out that they can compete with subway travel in certain areas and periods [19].Cycling distance, the number of commercial points of interest, metro ridership, and the distance to downtown are important variables in its integration with public transportation such as the metro [20].And people use bike-sharing services regularly during COVID-19 outbreaks [21].
For a docked bike-sharing system, infrastructure such as stations need to be established to provide services.Station density is an important positive factor for docked bike-sharing users to choose their destinations [22].Chen et al. gave a three-phase framework to determine the location of new bike stations [23].There are also studies using smartphone GPS data for station location design [24,25].Kuo et al. stated that for docked bike-sharing systems, it should not simply be a matter of installing more stations, but the focus should be on the overall improvement of service quality [26].Docked bike-sharing users are more likely to ride farther than dockless bike-sharing users [27].The distance traveled in a docked bike-sharing system follows a strong power law distribution [28].There is a cyclical pattern in their travel demand [29].Among the user groups, the proportion of elderly people is positively correlated with the use of docked bike-sharing [7].
Rebalancing is a very important activity for a bike-sharing system.Typically, the rebalancing process of a bike-sharing system can be described using a mixed integer linear programming model [30].Rebalancing is more of a response to morning and afternoon usage demand exceeding station capacity [31], and that demand can be stochastic [32] and needs to be described using, for example, a Markov process [33,34].Typically, rebalancing algorithms are designed to operate on stations with high demand by splitting the station [35,36].Most of the rebalancing models use heuristic algorithms for their solution process [37].For example, algorithms such as large neighborhood search algorithms, tabu search algorithms, and simulated annealing methods are used [38][39][40].The solution time of these algorithms is positively correlated with their search capability [41].At the same time, to impose some additional constraints on the rebalancing process, it may be necessary to design multi-objective optimization algorithms to encapsulate the impact of these constraints [34,42].In recent years, with the continuous development of machine learning technology, data-driven machine learning methods have been gradually applied to the optimization and solution process of rebalancing models [43,44].However, the rebalancing process has many problems and challenges waiting to be solved [45].
Zheng et al. noted that bike-sharing is currently an environmentally friendly practice as it achieves environmental benefits on all metrics except metal consumption [61].The promotion of bike-sharing can significantly reduce life cycle greenhouse gas emissions in the transportation sector [60].And the dockless bike-sharing system in Nanjing has been shown to obtain an emission reduction of 63.726 gCO 2 /eq/pkm [62].For the calculation of the emission factor of the bike-sharing system, a life cycle assessment was given by Bonilla et al. through the simulation data of the docked bike-sharing system and the dockless bike-sharing system [10].Luo et al. considered the impacts caused by rebalancing in their life cycle assessment of docked bike-sharing systems and dockless bike-sharing systems in the United States [11].Their rebalancing data for the dockless bike-sharing system was estimated based on the docked bike-sharing system, and they found that the main contribution of carbon emissions from bike-sharing was in the rebalancing phase.Kou et al. found that the bike-sharing system in the United States could only reduce less than 0.1% of the greenhouse gas (GHG) emissions from the transportation sector through their assessment [69].Chen et al. analyzed the full life cycle emissions of the dockless bike-sharing system of China's ofo dockless bike-sharing company and found that the main contribution of carbon emissions from shared bicycles is in the production stage [2].Mao et al. similarly assessed the full life cycle of shared bicycles in China and found that the production stage brings the greatest environmental impact, and its average contribution to the environmental impact is as high as 81.18% [13].Specifically, in Beijing, Wang et al. found in their assessment of the environmental benefits of dockless shared bicycles that they can obtain the desired environmental benefits with effective management [14].Finally, Chen et al. give the environmental benefits of shared bicycles in New York City based on realistically calculated rebalancing data [12].Table 1 gives a comparison between these studies.Rough estimates using empirical data or rough estimates using data from others' studies

Article Structure
The remainder of the paper is structured as follows.Section 2 covers materials and methods.Section 3 is the results.Section 4 is the discussion.Section 5 is the conclusion of the paper.

Data Description
In recent years, bike-sharing systems have become increasingly popular in cities around the world, where they play an important role in reducing air pollution and increasing sustainable transportation options in cities [60,70].In cities with large populations, long average commute times, and frequent use of public transportation, bike-sharing can have a greater impact on reducing air pollution than in cities with smaller populations [71].It significantly reduces commuting time and increases accessibility to work on a personal and spatial level [72].In an area with a well-developed bicycle infrastructure, bicycles may be faster than other modes of transportation [73].The freedom of bicycle parking makes it easy for residents and factory workers to use bicycles, leading to a high demand for the use of bike-sharing in residential and industrial areas [3].Bike-sharing is considered a convenient way to solve the "last mile" problem and connect users to the public transportation network [74].Further, during the COVID-19 pandemic, bike-sharing was considered to be a mode of transportation with a lower risk of infection [75].
Nanjing is located in China's Jiangsu Province, on the lower reaches of the Yangtze River, and is an important central city in eastern China, a nationally important research and education base, and a comprehensive transportation hub, with a current resident population of nearly 9.5 million.Nanjing began offering docked bike-sharing services in 2016.The data used in this paper are records of docked bike-sharing usage in Nanjing in 2017 as well as in 2023.The data contain information such as the time of borrowing and returning the bicycles, the station, and the serial numbers of the bicycles.In this section, data from March 2017 as well as March 2023 are used as examples to show the basic condition of the docked bike-sharing system in Nanjing.
Figure 1 gives a comparison of the use of docked shared bicycles in Nanjing over time.Figure 1a shows that the total daily mileage of Nanjing's docked bike-sharing system can reach two to five times as high in March 2017 as it did in March 2023.There were 5178 more docked shared bicycles in Nanjing in March 2017 than in March 2023, but both were located in the range of 30,000 to 40,000 bicycles.In addition, the number of stations in the docked bike-sharing system in Nanjing has increased by 386 in 6 years.Figure 1b gives the same conclusion as Figure 1a.The increase in the number of stations has not resulted in more frequent user use of the Nanjing docked bike-sharing system; on the contrary, there is two to five times less user demand for Nanjing docked bike-sharing in March 2023, seven years after the system was in operation, than in March 2017, the initial year of the system's operation.There is also a slight difference in the average cycling distance given in Figure 1c.Interestingly, the average distance per ride for users of docked bike-sharing in Nanjing is about 200 m shorter in March 2023, six years later, than in March 2017, six years ago.However, this may be due to COVID-19 epidemics or weather, for example.The average daily turnover per bicycle (the average number of times each bicycle is used per day) has changed from the original 3 to 4 times in March 2017 to about 1 time in March 2023, as shown in Figure 1d. Figure 1e shows that users also spent about 2 min less on average per ride in March 2023 than they did in March 2017.All of the above facts indicate a serious loss of user demand for the docked bike-sharing system in Nanjing.Its user demand has had a cliff-like decline compared to the early stage of the system construction, even though the number of station placements has been growing during these six years.Figure 2 gives the distribution of the average daily mileage traveled and the number of rides along with the cycling distance for the docked bike-sharing system in Nanjing.From Figure 2a, we can find that the peak of the number of rides occurs at trips ranging from 500 m to 2000 m, regardless of the period.Moreover, the total distance cycled in different intervals in March 2023 has a similar shape to the March 2017 distribution.However, there is a difference in the comparison of the number of rides in different cycling distance intervals given in Figure 2b.In March 2017, the peak of the number of rides was concentrated in the interval between 500 m and 1000 m, followed by a rapid decrease.In contrast, in intervals greater than 5000 m, the number of rides was approximately the same as in the 3500 to 4000-m interval.However, in Figure 2a, this interval contributes about 770 million kilometers of total cycling distance per month, which is the same as the 2000 to 2500-m interval.Conversely, come March 2023, there is a large change in the shape of the peak and distribution of the ridership.As can be seen in Figure 2b, the peak in cycling trips is in the interval of less than 500 m of cycling distance.Subsequently, the number of rides declined rapidly.This indicates that people's cycling habits have changed during these 6 years.In March 2023, not only the demand loss of Nanjing's docked bike-sharing system is severe, but people's cycling distance is also shorter and shorter, which corroborates with Figure 1c.Over the past six years, the share of shared bicycles in the "last mile" of less than one kilometer rides in Nanjing has increased by 9% compared to the total range.Given this dramatic change, perceptions of Nanjing's docked bike-sharing system, such as its sustainability in terms of CO 2 emissions and system profitability, may need urgent change.

Life Cycle Emission Factors for Docked Bike-Sharing System
Docked bike-sharing systems are considered to have positive benefits in terms of reducing carbon dioxide produced by the transportation sector [76].Life Cycle Assessment (LCA) provides a means to ensure the avoidance of shifting the burden of greenhouse gas (GHG) emissions by evaluating the use of natural resources by the docked bikesharing system and their environmental, economic, and social impacts throughout their life cycle [77].The life cycle carbon emissions of each docked shared bicycle include the total carbon emissions resulting from its four phases: manufacturing, operation, maintenance, and end-of-life [14] E li f ecycle bike = E Manu f acturing bike In this case, each station has both station and dock instruments, i.e.,

E
Operation bike , on the other hand, is closely related to the rebalancing of the docked bikesharing system in Nanjing and is given by the following equation: where d ride is the total docked shared bicycle cycling distance corresponding to the life cycle assessment period and E vehicle corresponds to the emission factor of the distribution vehicle.And f rebalancing is the rebalancing coefficient, whose approximation is given as follows: where d t rebalancing is the exact total rebalancing distance in period t and d t ride is the total cycling distance of docked shared bicycles in period t.Equation (4) approximates the rebalancing coefficient using the average rebalancing distance corresponding to each unit of cycling distance during period t to approximate the rebalancing mileage during the life cycle assessment.
Specifically, the life cycle emission factor for the docked bike-sharing system is where n bike , n dock , and n station represent the number of bicycles, docks, and stations, respectively, and B l represents the life expectancy of bicycles.

Rebalancing Coefficient
It has been shown that user demand for the docked bike-sharing system is uneven [78,79]; thus, rebalancing is necessary.For large-scale systems such as the docked bike-sharing in Nanjing, most of the existing studies have divided them into zones and then dispatched them [79].

Zoning for Rebalancing
Community detection algorithms are widely used for zoning in bike-sharing systems [79][80][81].Community detection is a large research area in complex networks, and many methods have been proposed in recent years.The well-known modularity methods are often used.However, even the directed modularity method has its weaknesses, i.e., it cannot detect the direction of the edges correctly.Moreover, it cannot detect patterns in the network [82].In addition, we need to find an algorithm that can specify the number of communities in advance since the bike-sharing company may have pre-determined the number of dispatch zones.Therefore, an EM-based community detection algorithm seems to fit the bill [83].However, the algorithm is only suitable for use in some unweighted directed networks, whereas if the usage data of Nanjing's docked bike-sharing system is constructed as a complex network, it is a directed weighted network.In this section, we will modify this community detection method to make it applicable to our goal of zoning and rebalancing.
Let ω ij denote the weighted adjacency matrix of the graph G with n nodes and a ij denote the unweighted adjacency matrix of the graph G with n nodes.We change a ij to ω ij in the iterative formula in the original model to obtain the following iterative formula: where c denotes the number of communities.θ rj is the probability that a node in community r has a link pointing to node i, and π r is the probability that a node is in community r when it is randomly selected in the network.q ir denotes the probability that a node i belongs to community r, and k i is the out-degree of node i.To obtain the final result of community detection, we need to iterate Equations ( 6)-(8).
To apply this method to the complex network of Nanjing's docked bike-sharing system, we define ω ij as the weighted adjacency matrix of the network, where the weights of the edges denote the number of borrowed bicycles.Specifically, the weight of an edge from node i to node j plus 1 indicates that a person borrowed a bicycle from station i and returned it at station j.There is a possibility that a person borrowed a bicycle from station i and returned it at station i.
However, there are two problems with the above-proposed method.Let us assume that in a given network G, there is a node j that has only out-degree but no in-degree.This is a common situation that exists in many networks, in which case θ rj will yield 0, which means that none of the communities will have connected edges pointing to node j from their nodes.Then in the next iteration step, q ir will yield 0 and the iteration process will collapse.This is because some nodes give connected edges with probability 0, which contradicts the assumptions of the method.The method assumes that the probability of each node belonging to each community is independent of other nodes.However, if the probability is 0, this means that it is no longer independent.Essentially, this method uses the most basic idea in probability theory, which is to express probability in terms of frequency.However, it is a fact that a frequency equal to 0 does not mean that the probability is equal to 0. Therefore, we rewrite Equation (7) as where ε is a very small positive random number close to 0. The second issue concerns weights.Although the method proposed above does not have theoretical problems, it may have computational problems.In a complex network constructed from the usage data of docked bike-sharing systems, a few nodes will have very large weights while others will have much lower weights, as shown in Figure 3, due to the uncertainty of their user requirements.This will make ∏ j θ ω ij rj very small; sometimes, the software may treat it as 0.Moreover, large weights dilute the importance of small weights because there is an order of magnitude difference between the weights, which can be found in the double logarithmic coordinate system in Figure 3.However, there are more nodes with small weights that also play an important role in the network.Therefore, we consider adjusting the weight of each node to make it more balanced.f (x) denotes the down-weighting function.f (x) must have the good property of being able to reduce the weight of each node to a finite interval [a, b].Let i be the node with the largest weight and j be the node with the smallest weight (except 0).a, b must satisfy Figure 3. Visualization of the weights of each connected edge in a network constructed using March 2017 data on the docked bike-sharing system in Nanjing.
Also, b − a must not be very large and all the weights must be evenly distributed in it.Apart from this, the definition of f (x) needs to be given according to the properties of the network and the problem faced.For example, if the weight gap between each node is large and almost all the nodes play an important role in the network and hence cannot be ignored, f (x) needs to be defined to narrow the gap.On the contrary, if the weight gap between each node is very small, or the weights are all close to 1, then f (x) needs to be defined to widen the gap and thus differentiate it from the unweighted network.Note that defining f (x) as a normalization process is inappropriate here.This is because it will not solve the computational problems described above.Moreover, large weights would make small weights less important since the order of magnitude difference remains large after the normalization process.For our bicycle network, we define where This function converts the weights to [1, 2].Ultimately, the improved method is given below To iterate Equations ( 13)-( 15), we define to denote the difference between two neighboring iterations.Here, t denotes the iteration step.If ∆ q < 10 −4 , we will stop the iteration.

Rebalancing Model for the Docked Bike-Sharing System
The Bike-sharing Rebalancing Problem (BRP) refers to the use of a series of distribution vehicles with maximum capacity constraints to pick up bicycles from stations with too high an occupancy level and deliver them to stations with too low an occupancy level and determine the vehicle routes to minimize the total cost of the entire dispatch process [30].Considering the shared bicycle is a special single type of good, the delivery demand at the bicycle redundant stations is considered zero and the pickup demand at the bicycle scarce stations is considered zero, so the shared bicycle rebalancing problem can be handled with the help of the PDVRP (Pickup-and-Delivery Vehicle Routing Problem) model.Here, we give a formulation of the BRP problem, which is built under the three-indicator formulation of the M-TSP (Multiple Traveling Salesman Problem) problem.K represents the set of distribution vehicles with index k; V is the set of warehouses and stations; and V 0 denotes the set of stations, N = |V 0 |; A = {(i, j)|i ∈ V, j ∈ V, i = j}.l ij represents the distance from station i to station j; x k ij is a 0 − 1 variable, with x k ij = 0 meaning that the kth vehicle does not pass through station i to station j, and x k ij = 1 vice versa; f k ij denotes the flow rate of the kth vehicle from station i to station j; u j denotes the imbalance degree of station j; Q represents the capacity of the vehicle; and v i is a dummy variable introduced to eliminate sub-loops.The specific formula is as follows: ∑ i∈V,i =h The formulation requires the minimization of the total distance traveled by distribution vehicles, where constraints (17) and (18) ensure that at least one distribution vehicle passes through and at most K distribution vehicles pass through each station; constraint (19) ensures that distribution vehicles entering a station will necessarily exit from that station; constraint (20) defines that all vehicles used depart from the warehouse; constraint (21) eliminates the sub-cyclic constraints; constraints (22) defines that the flow on each arc used does not exceed the vehicle capacity; constraint (23) strictly satisfies the delivery or pickup demand at each station; and constraint (24), constraint (25), and constraint (26) define the variable ranges.

Design Model Solution Algorithm and Solution Results
First, we try to use the solver for the exact solution.In this paper, we set the capacity Q of the distribution vehicles in the solution to be 50.At the very beginning, since we need to know the number of distribution vehicles and the number of selected vehicles is unknown until the optimal solution is solved, we determine the number of distribution vehicles (number of paths) according to the following equation: where N = |V 0 |.The imbalance of each station represents the difference between the bicycles parked at the end of the day's time at that station and the bicycles parked before the start of the day, with a positive value indicating |u i | more bicycles and a negative value indicating |u i | fewer bicycles.Figure 4 gives a visualization of the imbalance degree for each station on 1 March 2017, for the docked bike-sharing system in Nanjing.In addition to the above method, we give the number of vehicles K obtained by the heuristic algorithm based on the results of the heuristic algorithm presented below and reuse the solver for the rebalancing of that community.In addition, we define the warehouse of each community as one of the stations, which is calculated by the following equation: Due to the large size, the solver may not be able to obtain a satisfactory solution.Thus, we also design heuristic algorithms to solve satisfactory local optimal solutions for the shortest rebalancing path.First, we need to redesign the objective function to make it suitable for the solution of the heuristic algorithm.As shown in Figure 4, when setting the capacity of the distribution vehicle Q = 50, there are a considerable number of stations whose imbalance is larger than the capacity of the distribution vehicle.Thus, we split each station [35,36].For each station i except warehouse station u 0 , there are thus having At this point, each station is split into j logical stations with the same location, and they are re-notated as Here, there is always |µ| ≥ N − 1 that holds.
In the optimization process used in this paper, the distribution vehicle is allowed to carry a certain number of bicycles when it initially leaves the warehouse station.The number of bicycles to be initially carried needs to be derived from a calculation based on each distribution path.Under the condition of initially carrying bicycles, the sum of the cumulative imbalances calculated from each distribution path needs to satisfy the dynamic boundary conditions in Theorem 1.
Theorem 1.The dynamic boundary conditions under which a distribution vehicle can reach station i are given by the following inequality for the cumulative imbalance degree: where S i = ∑ i n=1 µ i represents the inventory of the distribution vehicle after rebalancing logical station i (which can take a negative value, excluding the number of bicycles carried when leaving the warehouse station), and let C start represent the number of bicycles carried by the distribution vehicle when it initially leaves the warehouse station.If S i does not satisfy Equation (32), the distribution vehicle needs to end that dispatch when it reaches station i − 1 and then return to warehouse station u 0 .
Proof of Theorem 1.It is easy to see that in every dispatching process, the following inequality always holds 0 thus, there are From the above two equations, it is easy to know By deforming the above two equations, we obtain as a result max In the proof, it follows from Equation (37) that It also holds that C start ≥ 0 is always true, and so the number of bicycles carried by the distribution vehicle on its initial departure from the warehouse station can be given by equation Consequently, we set the objective function as a set of sequences of logical stations µ i .Determine the interruption point θ i in each sequence M j by using the dynamic boundary conditions in Theorem 1.Let µ k 1 be the first interruption point θ 1 .When the first µ k i that does not satisfy the dynamic boundary conditions occurs in this sequence from the previous interruption point θ k−1 , it is noted as the next interruption point θ k , and M θ k is used to represent the segment of the sequence.Thus, the following can be obtained: which and d µ k i , µ k i+1 represents the distance between logical stations µ k i and µ k i+1 .
In this paper, we decided to use the Simulated Annealing (SA) method to calculate the shortest rebalancing path.The limiting behavior of the simulated annealing method has been well studied through Markov chains [84,85].In our present algorithm, the probability of acceptance of a non-optimal solution in the loop is the Boltzmann distribution The initial annealing temperature is taken as T 1 = 100.The coefficient of cooling for the annealing temperature is as follows: This cooling coefficient function will allow the algorithm to fully operate at high temperatures and slowly decline at low temperatures.In addition, after our attempts, the operation at temperatures greater than 60 is very important for the improvement of the quality of the final solution.The operators in the algorithm for generating new candidate solutions consist of the following three.The first is the swap operator, which randomly exchanges two logical stations in the sequence M j ; the second is the reverse operator; and the third is the transpose operator, which randomly divides the sequence M j into four segments of arbitrary length, i.e., M j = M 1 j M 2 j M 3 j M 4 j , and exchanges the positions of the two middle subsequences to obtain M new j = M 1 j M 3 j M 2 j M 4 j .One operator is randomly selected for the computation each time.The annealing is stopped when the optimal solution no longer changes after 150 consecutive annealing.

Modeling the Business of Docked Bike-Sharing Systems
Considering a project planning period of N l years, we discuss what the profit and loss of the business will look like over the next N l years if the docked bike-sharing system continues to operate at its current size for N l years.It is assumed that all infrastructure (stations) will be in service for S l years and that the bicycles will have a life expectancy of B l years.Based on the average size of the system, we examined the profit and loss scenarios of the firms over the next N l years.The basic information on the infrastructure size (total number of stations S, total number of bicycles B) is derived from the raw data on historical bicycle trips.

System Costs
The cost of a docked bike-sharing system is mainly composed of two aspects: the system construction cost and the operation and management cost.It is worth noting that this paper does not consider the cost of company administration and management, etc., but only the cost of hardware and dispatching.
The cost of system construction consists of two main aspects: the construction cost of the station and the procurement cost of the bicycles.The system construction stage is the main stage of consuming financial and material resources, characterized by large capital investment and short time.We assume that the construction cost of the station is S C and the procurement cost of bicycles is B C , so that the total construction cost of the system at the initial stage is The annual Operation and Management (O&M) Cost O total of a docked bike-sharing system is mainly the cost incurred by various tasks involved in the O&M process of the system, which mainly consists of rebalancing cost R C , maintenance cost M C , etc. Namely, Because rebalancing is often done by rebalancing employees and distribution vehicles, we divide the cost of rebalancing into two components: the cost of using distribution vehicles, R v c , and the cost of labor, R m c , which is First, consider the cost of using a distribution vehicle, which in this case can be trucks or e-trikes.It is assumed that in 2017, rebalancing operations were performed by trucks.For the year 2023, the distribution vehicle could be either trucks or e-trikes.Therefore, the cost of using distribution vehicles (trucks, e-trikes) can be obtained by the following equations, respectively where R P = 1 f rebalancing is the rebalancing coefficient, d y is the total monthly mileage traveled in year y, f p is the unit fuel price, f c is the fuel consumption per unit distance of the truck, t p is the truck procurement cost, and l t is the mileage of the truck at the end of its life.e p is the unit price of electricity, e d is the power consumption per unit distance, e c is the procurement cost of the e-trike, and e l is the end-of-life mileage of the e-trike.The power consumption per unit distance e d can be calculated based on the energy parameters of the e-trike used in the system.Assuming that the nominal battery voltage is e v , the battery capacity is e a , the motor power is e m , and the charging efficiency is η, first, we can calculate the range of the e-trike in a fully charged state In turn, we can obtain the power consumption per unit distance In addition, for the labor cost, we can obtain it by the following equation: where H c is the labor cost factor, v is the speed of the distribution vehicle (v t is the speed of the truck and v e is the speed of the e-trike), and H is the actual travel time of the distribution vehicle (the total rebalancing time minus the time of loading and unloading bicycles).
Next, we consider the maintenance cost of the docked bike-sharing system, M C .The maintenance cost consists of two main components, the cost of maintenance equipment, M v c , and the labor cost of maintenance personnel, M m c , i.e.,

Present Value Modeling of System Operational Benefits
Docked bike-sharing systems derive their revenues primarily from fees generated from riding.Next, we look at the profit and loss of the business for a certain period in the future, focusing mainly on when the business will be profitable.When analyzing the operational benefits of a docked bike-sharing system, we need to discount the operational benefits of the docked bike-sharing system in different years by taking into account the change in the time factor because of the long life cycle of the project and the varying monetary values at different stages.In this paper, we use four times the one-year LPR, Loan Prime Rate, as the interest rate r (the highest private lending rate allowed by Chinese law).The cumulative net benefit R n of the system in year y can be expressed as It is worth noting that in Equation ( 60), we also include the cost of bicycles, which is because bicycles are scrapped at the end of their useful life and the system needs to reintroduce bicycles in that year.

Results
In this section, we compare the emissions of Nanjing's 2017 docked bike-sharing system with those of the 2023 docked bike-sharing system to reveal whether the environmental benefits of Nanjing's docked bike-sharing system are sustainable over time.

Rebalancing Coefficient
For the docked bike-sharing system in Nanjing, there were 1086 stations built in 2017 and as many as 1476 stations in 2023.We begin with zoning for the 2017 and 2023 Nanjing docked bike-sharing systems.
Figure 5 gives the results of community detection based on the above-modified community detection algorithm.Figure 5a gives the results of the community detection of the docked bike-sharing network in Nanjing in March 2017.We chose to divide the 1086 stations into 10 communities.As can be seen in Figure 5a, there is a partial overlap in the edge location of each community.Therefore, as shown in Figure 5b, we manually adjusted the stations at the edge positions of each community so that they no longer overlap.Figure 5c gives the results of the community detection of the docked bike-sharing network in Nanjing in March 2023.We chose to divide the 1476 stations into 12 communities.However, the detection of communities resulted in 13 communities after geo-labeling.One of the communities contains nodes that are geographically distant from each other and are split by another community.Therefore, we define that community as two different communities.In addition, as in Figure 5b, we also manually adjusted the stations at the edge locations of each community so that they no longer overlap, and its final result is shown in Figure 5d.Our rebalancing operation will be based on the results of this community detection and will be performed separately within each community.Next, we solve the rebalancing coefficients based on the rebalancing model.Table 2 gives the results of the exact solution using real data from 1 March 2017, for the docked bikesharing stations in Nanjing.We perform rebalancing operations within each community.Gurobi 10.0.1 successfully solves the shortest rebalancing path for only one community.For all the other 9 communities, we encountered incidents where the Python 3.9.7 kernel died and restarted the kernel on its own after the program had been running for some time.One possible reason given by Gurobi is that the RAM size of the workstation is not large enough to solve the objective function.Meanwhile, none of the other 9 communities were even able to find a set of feasible solutions by Gurobi 10.0.1 before the Python kernel died.The results of SA's performance with this data are given in Table 3.It can be found that SA can give more satisfactory local suboptimal solutions, where the total absolute imbalance u total is calculated by the following equation: which represents the sum of the absolute values of the imbalances within each community.Total time (parallel) 1  1:07:41.4 1 The SA algorithm uses parallel computing techniques.For Gurobi, we did not solve the 10 communities in parallel since parallel computing techniques would have been utilized in the solution process within a single solution problem.

Calculation of Life Cycle Emission Factors
For the docked bike-sharing system in Nanjing, where [2] E Manu f acturing bike = 31.3711kgCO 2 eq/bike, ( E Scrap bike = 3.7339 kgCO 2 eq/bike.
Bicycles have a lifespan of 3 years [13].For the carbon emissions from the stations, we distribute them equally to each bicycle according to the lifetime of the station, which is taken to be 10 years [11].While [12] Êstation = 2008.1 kgCO 2 eq/station, ( From the above calculations, it can be found that the emission factor for the docked bike-sharing system in Nanjing in 2023 is 4.6 times higher than that in 2017.This is caused by the fact that the utilization of the docked bike-sharing system in Nanjing was significantly higher in 2017 than in 2023.As a snapshot, the average daily turnover rate per bicycle in 2017 was 4.04 times, while that value dropped to 1.17 times in 2023.

Comparison of Operations of Docked Bike-Sharing Systems in Nanjing, 2017 vs. 2023
Table 4 gives the specific values of the parameters of the system size.Table 5 gives the specific values of the parameters of the system cost.According to the formula given in Section 2.4.1 and the values of the parameters in Table 5, we can calculate the use cost R v c and labor cost R m c of the distribution vehicles of the docked bike-sharing system in Nanjing in 2017 and 2023, as shown in Table 6.According to the statistics of bicycle maintenance in each month, the average maintenance cost of a bicycle is 5.8969 yuan, where the data of the number and category of maintenance equipment are from Chen et al. [2] and the price data are from 1688.com.The number of bicycles that are used 1-31 times in a month of Nanjing's docked bike-sharing system in 2017 and 2023 is shown in Figure 6a.Based on the correspondence between the number of days riding shared bicycles and the number of repairs given by Chen et al. [2], we obtain the relationship between the number of repairs and the number of bicycles in one month of the docked bike-sharing system in Nanjing in 2017 and 2023 (as shown in Figure 6b).Thus, in year y, the number of annual maintenance bicycles of Nanjing's bike-sharing system, m y , can be calculated, and further, the cost of maintenance equipment of Nanjing's bike-sharing system in year y, M v c , can be calculated.Assuming that one person can only repair 20 bicycles a day and that the labor cost factor is the same as the labor cost factor in the rebalancing operation, the labor cost of the repairer can be calculated as M m c .The specific calculations are given in Table 7.At this point, we can calculate the initial construction cost and the average annual operating cost of the system based on the 2017 and 2023 sizes, as shown in Table 8.According to the reality, in 2017: trips under 30 min are free, and for trips over 30 min, each trip is priced according to the cycling time, and the price increases by RMB 1 for every additional 15 min of cycling time.In 2023: all trips under 30 min are RMB 1, and for trips over 30 min, the price increases by RMB 1 for every additional 15 min.From this, we can calculate the annual revenue P of the system in 2017 and 2023, respectively, as shown in Table 9.It is worth noting that in Equation ( 60) for = 4, 7, we also include the cost of bicycles, which is because bicycles are scrapped after 3 years of use and the system needs to reintroduce bicycles; and in Equation ( 59) for y = 2, 3, 5, 6, 8, 9. Taking LPR = 3.45% [90], we can calculate the profit and loss of the system for the next 10 years based on the system size and usage in 2017 and 2023, respectively, as shown in Table 10, where we consider two scenarios of truck rebalancing and e-trike rebalancing, respectively, in 2023.The results in Table 10 show that based on the size and use of the docked bike-sharing system in Nanjing in 2017, the business can be profitable in the ninth year.Based on the size and usage of the docked bike-sharing system in Nanjing in 2023, in the case of trucks, the enterprise has been losing money in the first five years until the beginning of the sixth year.As for the case of e-trikes, the enterprise has been losing money.

Comparison of Emission Factors
This paper shows that in 2023, compared with 2017, the demand for the use of the docked bike-sharing system in Nanjing is lost severely, and people's cycling distance is getting shorter and shorter.Seven years after the completion of Nanjing's docked bikesharing system, people seem to have faded their enthusiasm for shared bicycles.With roughly the same number of bicycles placed, the turnover rate per bicycle has dropped within these six years to one-fourth of what it was originally.
As a result, the sustainability of Nanjing's docked bike-sharing system changed significantly.The emission factor did not continue the lower value of 22.02 gCO 2 eq/km from 2016.In 2023, the emission factor for Nanjing's docked bike-sharing system increased by a factor of 4.6 to 101.42 gCO 2 eq/km.
In contrast, if e-trikes are used for the rebalancing operation, the emission factor decreases to 92.79 gCO 2 eq/km.The roles of the different components in the emission factor of the docked bike-sharing system in Nanjing are given in Figure 7.In Figure 7, it can be found that rebalancing only accounts for a very small part of its emission factor.In 2017, it only accounted for 10.79%.This compares to 8.86% in 2023.After changing the distribution vehicles from fuel trucks to e-trikes, the percentage decreases to 0.39% in 2023.However, in terms of emission factor, the action only reduces 8.62 gCO 2 eq/km, which is a small result.It can be seen that in China, the main emissions of the docked bike-sharing system are concentrated in the facilities of the system.Among them, the production, maintenance, and scrapping of bicycles contribute the most to the emissions, followed by the contribution of the docks and, finally, the stations.They make up the bulk of the amortized CO 2 from cycling.The specific proportionality between them is given in Table 11.

Environmental Benefits and Operational Sustainability
Second, we discussed the operational sustainability of the docked bike-sharing system in Nanjing.The results show that without considering the costs of company administration and management, and only considering the costs of hardware and dispatching, the payback cycle of Nanjing's docked bike-sharing system is very long.Users were very motivated to ride in 2017, and the first half hour of each ride was free for users.Perhaps this is one of the key reasons why user motivation was significantly greater in 2017 than in 2023.If user motivation were to be significantly sustained, it would take 9 years to recoup the cost of the investment and make a profit under the 2017 program without expanding the system.User motivation for the Nanjing docked bike-sharing system was clearly not sustained, and there was a significant loss of user demand for the system in 2023.In addition, the charging policy was changed and the free-cycling hours were removed, but the system was still being expanded all the time.Surprisingly, however, the number of years to recoup investment costs was reduced to six years in that unfavorable scenario.This may have been due to the coverage of costs by the price increase and the fact that dispatch costs have dropped significantly with the loss of user demand.The loss of user demand, while causing a significant increase in the emission factor, receives some cost-effective compensation.
However, is it the loss of user demand that led to the elimination of free cycling hours to cover the increased costs, or is it the elimination of free cycling hours to cover the costs that led to the loss of user demand?Construction of Nanjing's docked bike-sharing system began on 31 March 2015, and representatives of China's first dockless bike-sharing companies, ofo, and Mobike, entered Nanjing on 9 January, 2017 and 13 January, 2017, respectively [91].In the early stage of their entry into the market, there were numerous activities such as "free cycling for a week" and "recharge 100 RMB and get 210 RMB extra", etc. [91] A 9 March, 2017 news release showed that the top five dockless bike-sharing companies in terms of deployment had placed a total of 104,000 volumes of dockless shared bicycles in Nanjing [91].This paper uses data from the docked bike-sharing system in Nanjing in March 2017, with 39,302 bicycles placed in that month.The initial period of the most intense competition among bike-sharing companies also became the period of the highest usage of shared bicycles.Subsequently, in 2018, the number of shared bicycles in Nanjing rose to 449,000 [92]; that year, Nanjing limited the total number of shared bicycles to 300,000-380,000 [92]; at the same time, due to the disorderly expansion of bikesharing systems and excessive placement, a large number of "shared bicycle graveyards" appeared [93]; also at the end of the same year, ofo was on the verge of bankruptcy, with a large amount of money owed to its users and suppliers that remains outstanding (including the authors of this article) [94].Subsequently, in 2019, bike-sharing operators generally began to increase prices [95].In 2021, only three major players remained in the dockless bike-sharing market, namely Hello Bike, Meituan Bike, and Qingju Bike [96].In 2022, the bike-sharing system ushered in the second wave of price increases [97].In addition, as shown in Figure 8, although the number of demands in 2023 is going to drop significantly compared to 2017, the difference in the distribution of cycling time across periods is not as big.Even after the removal of the first half-hour free time, the percentage of first half-hour rides in March 2023 is comparable to March 2017, as shown in Figure 8b.To sum up, the bike-sharing system has gone through the early stage of subsidizing at a loss as well as capturing the market and is now entering the stage of seeking profitability, and non-sticky users may be lost in large numbers due to the price hike.Therefore, we believe that firstly the recession of the "bike-sharing fever" led to a decline in demand, and then the bikesharing operators generally began to increase prices, falling into a spiral of falling demand and rising prices.Further, after replacing the distribution vehicles with e-trikes, the Nanjing docked bikesharing system is not operationally sustainable in any way.Not only has the system failed to recover its costs within nine years, but the losses have continued to mount each year.This is likely due to a one-fifth drop in the capacity of dispatched e-trikes, resulting in a significant increase in the number of distribution vehicles, distances, and dispatch workers needed.However, the operation only reduces the emission factor by 8.62 gCO 2 eq/km, which is only 8.50% of dispatching with fuel trucks.It severely worsens the operator's key financial metrics without reducing CO 2 emissions by much.

Impact of COVID-19 on the Docked Bike-Sharing System
The COVID-19 epidemic affected all market sectors.The decline in demand in 2023 compared to 2017 may be partly attributed to the COVID-19 epidemic.During the COVID-19 epidemic, many studies observed a decline in demand for bike-sharing systems [98,99].Therefore, in addition to the decline in demand due to the recession of the "bike-sharing fever", a portion of the demand may have been lost due to the COVID-19 epidemic.The loss of users due to the popularity of COVID-19 will reduce the revenue of the system and at the same time reduce the total number of bicycle miles cycled, thus affecting the rebalancing of the system.The relationship is complex, but we can see from the results that lower demand for use in 2023 than in 2017 is rewarded with shorter payback years.In other words, the loss of demand accelerates the recovery of investment costs.
Specifically, COVID-19 primarily impacts the annual O&M costs O total of public bicycles.As shown in Equation (50), O total includes the rebalancing cost R C and maintenance cost M C .For the rebalancing cost R C , the outbreak of COVID-19 will have an impact on both the distribution vehicle usage cost R v c and the labor cost R m c , as shown in Equation (51).The distribution vehicle usage cost R v c is mainly determined by the rebalancing mileage, which is related to the overall usage of the docked bike-sharing system.The outbreak of COVID-19 has had a significant impact on the users' use of docked bike-sharing systems, which is mainly manifested in the changes in users' demand, average cycling distance, and cycling time [98][99][100][101].Similarly, the labor cost R m c is related to the overall usage of the docked bikesharing system.The rebalancing workload of the docked bike-sharing system is directly linked to labor costs, as shown in Equation (56).In addition, COVID-19 also has a direct impact on the labor cost factor H c .Therefore, if the effect of COVID-19 on the quantification of the model is considered, the coefficients α rebalancing COV ID−19 and α H c COV ID−19 need to be added to adjust the rebalancing cost and labor cost.Specifically, Equations ( 52), (53), and (56) need to be rewritten as However, in addition to COVID-19's impact on the demand for docked bike-sharing systems, government policies and initiatives, as well as operators' pricing strategies, all have an impact on demand, which interacts with each other and may need to be differentiated by further quantitative research to be quantified separately.Labor costs, on the other hand, are more complex.It is related to the labor contracts offered by the operators.For example, in Nanjing, where the docked bike-sharing system is government-led, the labor contracts with workers may be of the infinite duration type, meaning that H c may not be affected by the COVID-19 epidemic, whereas, if labor dispatch type of employment is used (which is more common in privately operated bike-sharing companies), H c could be significantly affected by the COVID-19 epidemic.In addition, there may also be correlations with inflation, unemployment, etc., which interact with COVID-19 outbreaks, and H c may need to be differentiated by further quantitative studies.

Comparison of Rebalancing Solutions
In Section 3.1, we found that the exact algorithm is not up to the task of solving the rebalancing calculation, while SA can solve the rebalancing model.In this section, we will compare different rebalancing algorithms.The rebalancing algorithms in this section will all be compared on the same dataset, using the same criteria.The dataset is the data from the Nanjing docked bike-sharing system on 1 March 2017, and all the solutions are programmed using Python 3.9.7 and implemented using Anaconda on a Win10 system with anIntel(R) Core(TM) i9-10900K CPU (Intel Corporation, Santa Clara, CA, USA) and 64 GB RAM.This includes the results of solving with the solver in Table 2 (Python 3.9.7 programming call to Gurobi 10.0.1) and the results of the SA algorithm in Table 3.The stopping criterion is the same as the SA algorithm except for the exact algorithm.In this section, we additionally implement the following four rebalancing algorithms for computing the shortest rebalancing path.
The first one is the Genetic Algorithm (GA).The number of chromosomes per generation P is an important metric for genetic algorithms.The suggested value of P is C ≤ P ≤ 2C [102], where C is the chromosome length.In this study, we take For the mutation rate, although some studies suggest taking [103] However, in our attempts we found that for this study, the best results were obtained when P m = 1.At this point, we sort the parent and child at the same time, select the new offspring, and terminate when the iteration is 500 generations.The crossover operator and fitness evaluation function here are the same as Goldberg et al. [104].The mutation operator is the reverse operator, also known as "2-Opt", which operates by deleting the dispatch paths between two random bicycle stations and reconnecting them so that they cross [105].The results are shown in Table 12.∆ represents the difference between the current solution and the optimal result among all methods (including the exact algorithm, if any).The ∆ values for the SA algorithm are given in Table 13.The second is the Immune Algorithm (IA), which has an antibody-to-antibody affinity threshold T ρ = 0.7 and a diversity evaluation index α = 0.95.The rest of the settings are identical to the GA algorithm.The results of the IA algorithm are given in Table 14.
Table 14.Results of the solution using IA.min f units are in meters in the table.

min f
Time ∆

Community
The third is the Tabu Search algorithm (TS).The tabu search algorithm is a metaheuristic stochastic search algorithm, which starts from an initial feasible solution, selects a series of specific search directions (moves) as a trial, and chooses the move that achieves the most changes in the value of a specific objective function.To avoid falling into local optima, a flexible "memorization" technique is used in TS search by creating a tabu table.
The tabu table records and selects the optimization process that has been performed and guides the direction of the next search step.The results of the TS algorithm are given in Table 15.The fourth type is the large neighborhood search algorithm (LNS).The large neighborhood search algorithm improves the quality of the solution step by step by alternating between the destroy operator and the repair operator after giving the initial solution.where the destroy operator randomly selects an α destory proportion of damage points, whereas the repair operator uses a greedy algorithm to repair to the position that minimizes the cost.In this study, α destory = 20%.The results of the LNS algorithm are given in Table 16.
Summarizing the above results, it can be found that although the exact algorithm solves better than the heuristic algorithm, the exact algorithm is unable to solve all of the rebalancing problems within the acceptable time and the acceptable equipment cost.A comparison of the rebalancing algorithms that appear in this paper is shown in Summary Table 17.And among all the five rebalancing algorithms, LNS has a very fast convergence rate but the quality of the solution is poorer.The SA algorithm obtains the best quality of the solution, but its convergence rate is the only algorithm that is more than one hour.GA and IA have a convergence rate comparable to TS, but the quality of the solution may be slightly inferior to the TS algorithm.Therefore, the SA algorithm can be chosen if a higher-quality solution is desired.For an acceptable and feasible solution in a shorter time, one can try the TS algorithm, which takes half the time of the SA algorithm.In this paper, we choose to use the SA algorithm to solve for the rebalancing coefficients.

Comparison with Existing Studies
Most of the emission factors for bike-sharing systems given by existing studies are from the USA as well as Beijing, China, and their results show great regional differences.Luo et al. [11] give emission factors for docked bike-sharing systems in the USA ranging from 26 gCO 2 eq/km to 147 gCO 2 eq/km, with a median value of 65 gCO 2 eq/km, and the emission factors for dockless bike-sharing systems ranging from 78 gCO 2 eq/km to 160 gCO 2 eq/km, with a median value of 118 gCO 2 eq/km.The emission factor for a docked bike-sharing system in New York, USA in 2020 given by Chen et al. [12] is 98.17 gCO 2 eq/km.The data-integrating multiple docked bike-sharing systems in the USA for 2016 given by Kou et al. [69] show an emission factor of 128 gCO 2 eq/km, whereas Wang et al. [14] give an emission factor of 65.16 gCO 2 eq/km for the docked bike-sharing system in Beijing, China, and for the dockless bike-sharing system, the emission factor is 315.06 gCO 2 eq/km.This implies that the emission factors of the bike-sharing systems vary greatly geographically.In this paper, the emission factor of the docked bike-sharing system in Nanjing in 2017 is given as 22.02 gCO 2 eq/km, and for 2023, the emission factor is 92.84 gCO 2 eq/km or 102.70 gCO 2 eq/km.
In addition, the proportions of the components in the emission factors of the docked bike-sharing systems in China and the USA also differ significantly.In the emission factors for the 2016 docked bike-sharing system in New York, USA, given by Luo et al. [11], the contribution of the bike manufacturing phase is only 5%, while the contribution of the rebalancing phase is 32.6%; for the average value in the USA, the contribution of the bike manufacturing phase is only 5.5%, the contribution of the bicycle manufacturing phase is only 5.5%, and the contribution of the rebalancing phase is 35.9%; for the USA dockless bike sharing system, the contribution of the bicycle manufacturing phase is 17.3% and the contribution of the rebalancing phase is 72.6%.In another study, Chen et al. [12] give an emission factor for the docked bike-sharing system in New York, USA, with a contribution of 11.3% for the bike manufacturing stage and only 11.8% for the rebalancing stage.On the contrary, for China, in the emission factor given by Chen et al. [2] for Beijing's 2017 dockless bike-sharing system, the contribution of the bicycle manufacturing stage was 83.8% and the contribution of the rebalancing stage was 4%.Wang et al. [14] gave a contribution of 83% from the bicycle manufacturing stage and 5% from the rebalancing stage in the emission factor for the dockless bike-sharing system in Beijing, while the contributions of the bicycle manufacturing stage and the rebalancing stage in the emission factor for the docked bike-sharing system located in Beijing were 40.9% and 17.8%, respectively.In our study, the contribution of the bicycle manufacturing stage of Nanjing's docked bike-sharing system in 2017 was 59.18%, the contribution of the rebalancing stage was 10.78%, the contribution of the bicycle manufacturing stage of Nanjing's docked bike-sharing system in 2023 was 51.78%, and the contribution of the rebalancing stage was 8.88%.It can be seen that the compositional share of the emission factors of bicycle sharing in China and the United States is very different.However, in terms of individual countries, different studies for China and the United States have yielded similar emission factor shares, respectively.
As for the emission factor share of the docked bike-sharing station facilities, the 2016 data share of the New York docked bike-sharing system in the United States given by Luo et al. [11] is 52.9%, and the average share across the United States is 51.1%.The data share of the New York docked bike-sharing system in the United States given by Chen et al. [12] is 66.2%.The data percentage of Beijing's docked bike-sharing system given by Wang et al. [14] is 35%.In this paper, the percentage of this data for Nanjing, China's docked bike-sharing system is 30.3% in 2017; the percentage of this data is 39.33% in 2023.Again, the existing studies produce a huge gap between China and the United States; and for China or the United States only, the data between different studies are similar and can be corroborated with each other.

Geographic Patterns between Environmental Benefits and Financial Operations
In the above section, we divided all stations in 2017 into 10 communities; similarly, all stations in 2023 were divided into 13 communities.In this section, we explore the operational revenues and costs within each community to analyze the geographic patterns of financial operations.Specifically, we calculate the revenues and costs within each community separately, including system construction costs and O&M costs.The benefits of each cross-community trip are divided equally between the two communities connecting the trip; similarly, each cross-community trip is divided equally between each community.Based on this, the total system construction and maintenance costs can be allocated to each community based on the ratio of the total number of trips supported by each community to the total number of trips.
Table 18 gives the operational benefits of the system by community according to the 2017 benchmark.We portray the geographic distribution of docked bike-sharing systems in Nanjing by community.It can be found that there is a huge difference in the geographical distribution of system operation.Among them, communities 1, 2, 3, 7, 8, and 10 all have negative benefits in the 9-year operation period.In particular, the losses of communities 1, 3, 7, and 10 increase year by year.Communities 4, 5, and 6 have positive returns in the third year.In contrast, community 9 has a longer payback time.From Figure 5b, it can be found that those communities that can realize a profit during the operation period are all located in the most prosperous main city of Nanjing.
Table 19 gives the sub-community operational benefits of using trucks for rebalancing based on the 2023 baseline.Communities 3, 4, 5, 7, 10, and 13 are unable to achieve a positive return within the nine-year operating period, and their losses all widen each year.Communities 1, 6, 9, and 12 can achieve positive returns in the third year; communities 2 and 8 take five years to recover their investment costs; and community 11 takes six years to become profitable.From Figure 5d, it can also be found that those communities that can realize a profit during the operation period are all located in and around the most prosperous main city of Nanjing.Compared to 2017, however, the area covered by non-loss-making communities in 2023 has expanded.

Community
Table 20 gives the sub-community operational benefits of rebalancing using e-trikes based on the 2023 baseline.Although the results in Table 10 have shown that the operational program is not sustainable, communities 1, 8, 9, and 12 can make a profit over the 9-year operational period, while the rest of the communities continue to lose money.This suggests that customizing the rebalancing scheme by sub-communities may allow for better environmental benefits.In other words, using e-trikes for rebalancing in communities 1, 8, 9, and 12 and trucks for rebalancing in the rest of the communities to further reduce the carbon emissions of the system are viable options.Under this scenario, the emission factors of the Nanjing docked bike-sharing system are E li f ecycle bike,2023,mix = 99.55 gCO 2 eq/km. ( Table 20.Operational benefits of sub-communities using e-trikes for rebalancing in 2023.

Community
The operational benefits of the system under this scenario are given in Table 21.In comparison with Equations ( 73) and (74), this emission factor is smaller than the emission factor for rebalancing with all trucks and larger than the emission factor for rebalancing with all e-trikes.However, with this emission factor, the operator can achieve financial sustainability.

Policy Proposals
Through the above research, we can find that the reduction in demand makes the emission factor larger.Thus, in terms of environmental protection and carbon dioxide emission reduction, the encouragement of green travel behaviors using docked bike-sharing systems can protect the environment and reduce carbon dioxide emissions.Consequently, the government should increase its initiatives for green travel behaviors such as bike-sharing.
However, the results in Table 10 point to the fact that greater demand for docked bike-sharing systems may result in longer payback periods for operators.Consequently, increasing the use of docked bike-sharing systems on top of the existing ones, such as through government initiatives, would hurt the operators.Therefore, the government should provide certain incentives or compensation to the operators in the context of promoting low-carbon mobility.
In addition, the benefits of some emission-reduction actions, such as mandating operators to use more environmentally friendly e-trikes for rebalancing operations, cannot be matched with the costs incurred by the operators (Figure 9).As a result, the government may not need to give more guidance on rebalancing.However, the operational benefits of the system may have a geographical pattern.Therefore, if trucks and e-trikes are scheduled for dispatch according to geographic distribution (e.g., community detection in this paper), it may be possible to achieve both environmental benefits and financial operational sustainability.In Section 4.5, we point out that rebalancing operations account for a small percentage of the emission factor of bike-sharing in China.Instead, the government should impose certain restrictions on bike manufacturing and station setup manufacturing, such as the use of more environmentally friendly materials, to further reduce emissions from docked bike-sharing systems.Finally, in the discussion in Section 4.5 we found that research on the environmental benefits of bike-sharing systems varies greatly at the national level.This means that for different countries, the results of existing studies cannot directly guide decision-making in other countries.In contrast, for different cities within a country, the results are similar.Therefore, cross-references and policy recommendations between cities within countries are possible.In the case of this study, it is instructive for docked bike-sharing systems in other cities in China.However, this study does not apply to cities outside of China, such as those in the United States.Cities in the USA need to refer to their country-specific studies when developing policies for docked bike-sharing systems.

Conclusions
In this paper, the sustainability of Nanjing's docked bike-sharing system is investigated by comparing the data from March 2017 and March 2023 for its system.The sustainability of the system is divided into two parts: one is its environmental sustainability, which is mainly measured using CO 2 emissions, and the other is its financial sustainability, which is calculated by modeling its revenue.
In calculating the emission factor of the docked bike-sharing system in Nanjing, we simulated the whole rebalancing process.Since the system is very large, we first modified a community detection method to realize the region zoning.Second, since the solver could not work properly under such a large-scale system, we used three heuristic algorithms for comparative calculations.To be able to apply these heuristics successfully, we rewrite the objective function into the form of a suitable heuristic algorithm and and prove dynamic boundary conditions for the new objective function.Ultimately, the simulated annealing method is chosen by comparison for the subsequent rebalancing coefficient calculation.We define the rebalancing coefficients and approximate the rebalancing coefficients by solving the one-week dispatch journey for estimating the approximate daily rebalancing journeys.Eventually, the emission coefficients of the docked bike-sharing system in Nanjing in 2017 and 2023 are obtained, and there is a 4.6-fold difference between them.In addition, whether one uses e-trikes or fuel truck dispatching has less impact on the emission factor of the docked bike-sharing system in Nanjing.
In estimating the financial sustainability of Nanjing's docked bike-sharing system, we considered the system's revenues over the next nine years, including facility input costs, facility replacement costs, dispatch costs, labor costs, maintenance costs, and the time value of money.The results show that in 2017, when user demand for use of Nanjing's docked bike-sharing system is high, it will take eight years for the investment to be recouped; conversely, in 2023, when there is a significant loss of demand, the investment will be recouped in only six years.Additionally, switching distribution vehicles from fuel vehicles to e-trikes would result in an unsustainable system.After the conversion, it would severely deteriorate the operator's key financial indicators with only a reduction in the emission factor value of 8.62 gCO 2 eq/km, and not only would the investment cost be unrecoverable, but the losses would widen year after year.We also explore the geographic patterns between environmental benefits and financial operations and find that arranging a mixed form of dispatch based on the detection of communities can balance environmental benefits and financial operations to reach dual sustainability.In addition, this paper compares the performance of some rebalancing algorithms and selects the best-performing one for rebalancing.Further, we find that the study on the emission factor of bike-sharing can obtain similar results within countries, while it varies greatly between countries.
This paper has the following shortcomings.First, the model constructed in this paper can only qualitatively describe the impact of COVID-19 on docked bike-sharing systems.Since the impact of COVID-19 on the docked bike-sharing system is at the overall level, quantifying that impact involves multiple aspects of research and requires more sensitive data that are harder to obtain (specific labor contract information, operators' pricing strategies, etc.), and thus we are currently unable to quantitatively calculate that impact.Second, due to limitations in time, computing power, and data access, this paper only uses nearly two months (8 weeks) of data per year.The number of data in this paper may not be large enough.More data may be able to obtain more accurate results and exclude more hidden random factors.

Abbreviations
Variables and abbreviations used in this paper and their corresponding explanations:

Figure 1 .
Figure 1.(a) The total number of miles traveled daily by docked shared bicycles in Nanjing.(b) The total number of times docked shared bicycles are used daily in Nanjing.(c) The average number of miles traveled per trip daily by docked shared bicycles in Nanjing.(d) The daily turnover rate of docked shared bicycles in Nanjing.(e) Average time per trip daily for docked shared bicycles in Nanjing.

Figure 2 .
Figure 2. (a) Distribution of the average daily mileage traveled by docked shared bicycles with cycling distance in Nanjing.(b) Distribution of the total number of daily uses of docked shared bicycles with cycling distance in Nanjing.

Figure 4 .
Figure 4. Imbalance at each station of the docked bike-sharing system on 1 March 2017, in Nanjing.

Figure 5 .
Figure 5. (a) Community detection results for the March 2017 complex network.(b) Community detection results for the complex network in March 2017 after manual adjustment.(c) Community detection results for the March 2023 complex network.(d) Community detection results for complex networks in March 2023 after manual adjustment.

Figure 6 .
Figure 6.(a) Relationship between the cumulative number of bicycle cycling days and the number of repairs.(b) Relationship between the number of repairs in a month and the number of bicycles in the docked bike-sharing system in Nanjing in 2017 and 2023.

Figure 7 .
Figure 7. CO 2 emissions per kilometer of cycling from sub-sources.

Figure 8 .
Figure 8.(a) Distribution of the number of docked shared bicycle rides with cycling time in Nanjing.(b) Percentage distribution of the number of docked shared bicycle rides in Nanjing with cycling time.

Figure 9 .
Figure 9. Use e-trike to rebalance shared bicycles.(a) Use e-trike to rebalance shared bicycles at night; (b) Use e-trike to rebalance shared bicycles during the day.
Life cycle emission factor for the docked bike-sharing system E Manu f acturing bike Life cycle emission factors for the manufacturing stage of docked shared bicycles E Maintenance bike Life cycle emission factors for the operational phase of the docked bike-sharing system E Operation bike Life cycle emission factors for the maintenance phase of the docked bike-sharing system E Scrap bike Life cycle emission factors for the scrapping stage of docked shared bicycles

Table 1 .
Comparison of emission factor studies of bike-sharing systems.

Table 2 .
Results of exact solution using Gurobi 10.0.1.minf units are in meters in the table.And when Q = 50, the distribution vehicles are fuel trucks.Due to the requirement of CO 2 emission reduction, we also calculated the rebalancing coefficients for dispatching by e-trikes in March 2023.At this time, Q = 10 and

Table 3 .
Results of the solution using SA.min f units are in meters in the table.

Table 4 .
Values of system scaling-related parameters in 2017 and 2023.

Table 5 .
The values of the relevant parameters.

Table 6 .
Initial construction costs and average annual operating costs in 2017 and 2023.

Table 7 .
Number of maintenance bicycles and cost of maintenance equipment in 2017 and 2023.

Table 8 .
Initial construction costs and average annual operating costs for 2017 and 2023.

Table 9 .
System annual revenues for 2017 and 2023.

Table 11 .
Contributions of rebalancing, bicycles, docks, and stations to the emission factors of a docked bike-sharing system in Nanjing.

Table 12 .
Results of the solution using GA.min f units are in meters in the table.

Table 13 .
Results of the solution using SA.min f units are in meters in the table.

Table 15 .
Results of the solution using TS.min f units are in meters in the table.

Table 16 .
Results of the solution using LNS.min f units are in meters in the table.

Table 17 .
Summary of the rebalancing algorithms that appear in this paper.

Table 18 .
Operational benefits of sub-communities in 2017.

Table 19 .
Operational benefits of sub-communities using trucks for rebalancing in 2023.

Institutional Review Board Statement:
Not applicable.Not applicable.Restrictions apply to the availability of these data.Data were obtained from Nanjing Public Bicycle Co. and are available from the authors with the permission of Nanjing Public Bicycle Co.The authors declare no conflict of interest.
When a node is randomly selected in the network, the probability that it is in the community r q irProbability that a node i belongs to community r k iThe out-degree of node i − 1 variable, representing whether or not the kth vehicle passes through the connecting edge from station i to station j f k ij Flow of kth vehicle from station i to station j u j