Estimation of Ecosystem Services Provided by Street Trees in Kyoto, Japan

: Street trees are integral components of urban green infrastructure. The importance of beneﬁts provided by street trees has motivated the development of various tools to quantify the value of ecosystem services. The i-Tree Eco is a widely applied method for quantifying urban forest structure, ecosystem services, and values. Since its ﬁrst release in 2006, i-Tree Eco has been successfully utilized in over 100 countries around the world. This study described one of the ﬁrst applications of the i-Tree Eco international project in Kyoto, Japan, by customizing the models and parameters to enhance the accuracy of analysis results. Kyoto’s street trees are prominently dominated by Ginkgo ( Ginkgo biloba L.), Trident Maple ( Acer buergerianum Miq.), Japanese Zelkova ( Zelkova serrata (Thunb.) Makino.), Tuliptree ( Liriodendron tulipifera L.), Flowering dogwood ( Cornus ﬂorida L.), London Planetree ( Platanus × acerifolia ), Plum/cherry ( Prunus spp.), and Weeping willow ( Salix babylonica ), which account for 92% of the 1230 sample trees and deliver ecosystem service beneﬁts at US$71,434.21 annually or US$58.07/tree/year. The annual value of each function was estimated at US$41.34/tree for carbon storage and sequestration, US$3.26/tree for stormwater runoff reduction, US$11.80/tree for adverse health mitigation effects, and US$1.67/tree for energy savings. The street tree species of Kyoto city that produce the highest average annual beneﬁts are among the largest trees currently in the population, including P. × yedoensis (US$225.32/tree), Z. serrata (US$123.21/tree), S. babylonica (US$80.10/tree), and P. × acerifolia (US$65.88/tree). Our results demonstrated a comprehensive understanding of street trees beneﬁts for Kyoto city, providing baseline information for decision-makers and managers to make effective urban trees management decisions, developing policy, and setting priorities.


Introduction
Urbanization, one of the most transformative trends in the 21st century, is accompanied by increasing populations and socioeconomic activities concentrated in cities. In recent decades, urbanization adversely impacts urban ecosystems and environmental quality through phenomena such as the urban heat island effect, air pollution, and alterations to hydrological systems [1,2]. Since the publication of the Millennium Ecosystem Assessment (MEA) [3] and The Economics of Ecosystems and Biodiversity (TEEB) report [4], ecosystem services have gained broader attention in many parts of the world [5,6]. Ecosystem services (hereafter referred to as ES) refer to the life-support functions performed by natural ecosystems that underpin humanity's most fundamental sources of well-being [7]. The strong desire to develop a sustainable urban environment that delivers the ES has encouraged policymakers and scholars to direct their attention to evaluating the potential of urban trees to mitigate environmental degradation [1].
Street trees are recognized as integral components of urban ecosystems, which can improve environmental quality by providing significant ecological benefits [8]. There is a growing body of literature that highlights the ES provided by street trees, such as storing carbon [9,10], regulating the air quality [11], as well as improving the streetscape and amenity [8,12]. Street trees also appear to be a feasible option for ameliorating the urban heat-island effect [13][14][15]. Moreover, trees play an important role in urban catchment hydrology through canopy interception and soil infiltration of rainfall [16,17].
Despite the high level of scholarly recognition of street tree benefits, many local governments fail to recognize the importance of street trees due to unknown economic values, while the costs of damage by trees such as leaf litter and infrastructure damage are widely reported [8,18]. Evidence of the economic worth of street trees in monetary terms is essential for decision-makers as it offers baseline information for long-term street-tree management and maintenance practices [19]. To understand the potential ES more fully and maximize the benefits of urban trees, several urban forest models have been developed and implemented in different cities [5,20]. The most frequently used model is the i-Tree software developed by the US Forest Service (www. itreetools.org), which provides a methodology using field data to assess the biophysical state and economic value of urban and community forests. Many studies based on the i-Tree software have demonstrated that the monetary values of ecosystem benefits provided by street trees exceed the annual cost of tree management. A study conducted in New York and Indianapolis showed that every $1 spent per year on tree-care-related expenditure resulted in US$5.6 and US$6.09 worth of ecosystem benefits, respectively, which was a greater economic benefit than that reported from any other city to date [21,22].
Although i-Tree tools have been extensively used across the US and European countries over the past 10 years, limited research using this method has been conducted in other parts of the world [1]. There are still uncertainties in applying US-based ES and associated benefit quantification models to other countries without appropriate modifications. For example, the fact that the acquisition of site-specific parameters is unavailable or insufficient will indicate possible inaccuracies of the results [5,23]. To improve model functionality in other countries, Hirabayashi et al. [24,25] conducted pilot studies in Japan, in which the accuracy of analysis results was enhanced by largely customizing the models and their parameters.
This study describes the application of i-Tree Eco with customized models and parameters for monetarization of the street tree benefits in Kyoto, Japan. The objectives of this study were as follows: (1). To demonstrate the procedure for implementing the i-Tree Eco project in Japan, including collecting all relevant data sources and explaining customizations conducted for each model. (2). To present objective data on the value of ecosystem services provided by street trees in Kyoto city as baseline data for evaluating returns on tree management investment. (3). To appeal for a better understanding of street trees as green assets, which could help ensure adequate tree maintenance and lead to future improvements in management and plans.

The City of Kyoto
Kyoto City, the capital of Kyoto Prefecture, is located in the central part of Honshu, Japan (Figure 1), which has a humid subtropical climate with hot, humid summers (June-August), and cold, dry winters (December-February). Due to its basin topography, Kyoto experiences the most severe summer heat in Japan, with a mean daily maximum temperature of around 37 • C in August [26]. In the central urban area of Kyoto City, given the limited large green spaces, street trees are recognized as an irreplaceable green infrastructure for mitigating urban heat island effects and providing critical ecological connectivity to promote faunal abundance and diversity [27].
in accordance with the i-Tree Eco field manual [29]. Tree species and adjacent land use were verified in the field. Health condition, crown light exposure, and percent crown missing were estimated by one person by visual inspection. Total tree height, height to crown base, diameter at breast height (DBH), crown width was measured using a Trupulse 360 laser rangefinder (Laser Technology, Inc., CO, USA) and a diameter tape.

Estimation of Ecosystem Services By i-Tree Eco
i-Tree Eco estimates urban and community forest features such as biomass and leaf area and then, based on these properties, quantifies the ES provided by the forests. The system comprises of three components: model codes written in computer languages, parameters for the models (e.g., coefficients in model equations), and input data (e.g., weather data).
Currently, for a total of about 40 countries officially supported, input data such as tree species, weather, upper air, air quality, and location-related data (such as coordinates and population) are stored in the i-Tree server computers. Using these data, the analyses can be conducted in these 40 countries, with only the tree data prepared by users. Because the species, weather, and upper-air data stored in the server are globally available, even users in other countries can use i-Tree Eco with their local location, air quality, and precipitation data uploaded to the server via i-Tree Database (ITDB) (www. database.itreetools.org). However, this path only enables users to use their local data; it is impossible for the model itself and parameters to be modified for their use in their own countries.
The ES provided by trees that i-Tree Eco can estimate include (1) carbon storage and sequestration, (2) air pollution removal, (3) human health effects associated with air pollution removal, (4) heating and cooling energy savings in houses, and (5) stormwater runoff avoided. The models for (1), (2), and (5) are readily applicable to countries other than the US using the local data that are made available with ITDB. On the other hand, since

Collecting Tree Data
With reference to the basic plan for greening Kyoto City, the green corridor zone was selected as the research area, which embraces the east to west extent from Shirakawa-Dori Street to Nishioji-Dori Street and stretches north to south from Kitayama-Dori Street to Jyujo-Dori Street [27]. The total area is 48.85 km 2 , covering seven wards with a total population of 451,462 [28]. Field sampling of 1230 street trees was conducted from June to October 2018. Approximately 10% of street trees on the 41 streets were sampled at equal distances to generate the tree inventory ( Figure 1). Field measurements were performed in accordance with the i-Tree Eco field manual [29]. Tree species and adjacent land use were verified in the field. Health condition, crown light exposure, and percent crown missing were estimated by one person by visual inspection. Total tree height, height to crown base, diameter at breast height (DBH), crown width was measured using a Trupulse 360 laser rangefinder (Laser Technology, Inc., CO, USA) and a diameter tape.

Estimation of Ecosystem Services by i-Tree Eco
i-Tree Eco estimates urban and community forest features such as biomass and leaf area and then, based on these properties, quantifies the ES provided by the forests. The system comprises of three components: model codes written in computer languages, parameters for the models (e.g., coefficients in model equations), and input data (e.g., weather data).
Currently, for a total of about 40 countries officially supported, input data such as tree species, weather, upper air, air quality, and location-related data (such as coordinates and population) are stored in the i-Tree server computers. Using these data, the analyses can be conducted in these 40 countries, with only the tree data prepared by users. Because the species, weather, and upper-air data stored in the server are globally available, even users in other countries can use i-Tree Eco with their local location, air quality, and precipitation data uploaded to the server via i-Tree Database (ITDB) (www.database.itreetools.org). However, this path only enables users to use their local data; it is impossible for the model itself and parameters to be modified for their use in their own countries.
The ES provided by trees that i-Tree Eco can estimate include (1) carbon storage and sequestration, (2) air pollution removal, (3) human health effects associated with air pollution removal, (4) heating and cooling energy savings in houses, and (5) stormwater runoff avoided. The models for (1), (2), and (5) are readily applicable to countries other than the US using the local data that are made available with ITDB. On the other hand, since the models for (3) and (4) were originally developed based on the methods and data applicable to the US only, the application of these models to countries outside the US is very limited.
This study estimated (1) to (5) based on street tree measurements in Kyoto. The i-Tree Eco's limitations existing outside the US were greatly reduced by modifying the model codes, as well as parameters with the cooperation of an i-Tree developer. Table 1 presents the input data/parameters for each model with a notation of customization conducted in this study, and the following sections explain each model and customizations conducted for models of (3) and (5).

Carbon Storage and Sequestration
It was necessary to identify the correspondence between tree species found in Kyoto city and those in i-Tree Eco's species database to calculate carbon storage and sequestration into trees. Other than that, no customization for the codes, parameters, and data was necessary. Based on the property of each tree species, dry biomass for woody parts as well as leaves was calculated using the methods described in Nowak et al. [51,52]. Half of the dry biomass of trees was estimated as the carbon stored in trees. The growth of a tree was estimated for each tree species based on the health condition and planted site characteristics (e.g., crown light exposure) of the tree. The gross amount of carbon sequestered annually into a tree was then calculated from the difference in estimates of carbon storage between the current and next year [51,52].

Air Pollution Removal
Based on tree structures such as tree cover and evergreen percent in the study area as well as the leaf area index (LAI) estimated with i-Tree Eco, the removal of carbon monoxide (CO), nitrogen dioxide (NO 2 ), ozone (O 3 ), particulate matter less than 2.5 µm (PM 2.5 ), and sulfur dioxide (SO 2 ) were estimated as described by Hirabayashi et al. [53,54] and Nowak et al. [55]. Input data for the model (i.e., surface weather, upper air, and air pollutant data) from the local monitoring stations were employed here (Table 1). Although the parameters for the model were not optimized, the model itself was optimized to use local measurements of solar and net radiation rather than calculating these based on the extraterrestrial solar constant, coordinate, and other atmospheric properties (e.g., ozone depth and albedo), which is the default model implementation in i-Tree Eco.

Human Health Effects Associated with Air Pollution Removal
i-Tree Eco estimates avoided adverse health incidences and costs associated with changes in NO 2 , O 3 , PM 2.5 , and SO 2 concentration due to the removal of pollution by trees [55,56] with BenMAP [57] incorporated into i-Tree Eco. BenMAP was developed by the US Environmental Protection Agency by consolidating the human medical records and air quality measurements across the US with the knowledge gained from statistical analyses of those data.
Based on seven concentration change metrics, 13 adverse health endpoints can be analyzed with BenMAP. The concentration change metrics include the annual mean of daily 1 h maximum (1Max), daily mean for 8-10 a.m. (3Mean), daily mean for 6-9 a.m. (4Mean), daily maximum for 8 h moving average (8Max), daily mean for 9 a.m-4 p.m. (8Mean), daily mean (24Mean), and the quarterly mean of the daily mean (24MeanQ). When integrating BenMAP's 85 health effect functions into i-Tree Eco, the incidence multiplier, IM i (case/ppb/person or case/µg/m 3 /person), which is the incidence per unit change in concentration per person for each county, was estimated from Equation1 below. Similarly, the value multiplier, VM i (US$/ppb/person or US$/µg/m 3 /person), which is medical expenses per unit change in concentration per person for each county, was estimated from Equation (2).
This integration allowed for the calculation of reduction in adverse health incidence (∆I) and reduction in medical expenses (∆V) based on the ∆C i provided by the trees, P i, and these multipliers. For example, reduction in HAR associated with NO 2 for 0-to 99-years-olds can be calculated as follows: where P 58 , P 59 , and P 64 are populations for age groups 0-14, 15-64, and 65-99, respectively. ∆C 58 and ∆C 59 are both changes in 1Max of NO 2 , and ∆C 64 represents the change in 24Mean of NO 2 . To calculate the reduction in the medical expenses, the corresponding VM i is used in Equation (3). For Kyoto, the medical records for each age group associated with air quality measurements that could replace the BenMAP analyses were not readily available. Hence, the parameters, IM i and VMi, were modified for Kyoto based on the assumption that the response of humans to air pollution is the same, whether it is in the US or in Japan. In this process, the years for the maximum (baseline year) and minimum (control year) annual mean concentrations from 2010 to 2016 were first identified for each of the four air pollu-tants in Kyoto. All the US counties were then searched to identify those that had the closest concentration for baseline and control years as well as the change between the two years for each of the seven metrics for each air pollutant. A reference US county for each metric for each air pollutant was assigned to Kyoto through these processes. Adverse health incidence for each endpoint for each age group was adjusted based on the population ratio between Kyoto and the reference counties. Monetary values per incidence were adjusted based on the ratio (46%) of the mean medical expenses between Japan and the US [36]. For WLD and SLD, the monetary value per incidence was adjusted based on the ratio (65%) of the mean household income between Japan and the US [37]. For the monetary value for M, the median value for a statistical life (VSL) derived from the literature in Japan [38] was employed.

Heating and Cooling Energy Savings in Houses
Based on McPherson and Simpson [9] and Nowak et al. [58], i-Tree Eco estimates changes in heating and cooling demands for houses with two or fewer floors based on shade, windbreak, and transpiration effects by trees with 6 m or taller and located within 18 m from the house and windbreak effects by other buildings.
In i-Tree Eco, for 11 climate regions in the US, the base values of CO 2 emission change because of changes in demands for cooling and heating due to tree shade and windbreak effects are stored in look-up tables for the combination of the three house vintages (pre-1950/1950-80/post-1980), leaf type (deciduous/evergreen), tree height (6-10/10-15/15 m or taller), the distance between a tree and a house (0-6/6-12/12-18 m), and eight directions from a house to a tree. In addition, the base CO 2 emission changes because of changes in the heating and cooling demands affected by the transpiration from trees and the windbreak by buildings are stored in look-up tables for 10%, 30%, and 60% of trees and building covers combined in each climate region.
It was ideal if this kind of base values of CO 2 emission change due to tree effects were readily available in Kyoto, but it was not the case. Thus, based on the assumption that human's demand for cooling and heating is the same in the US and Japan, these look-up tables were used in this study by referencing a US climate region that best fits the climate in Kyoto. The approach requires four steps: (1) from 16 cities representing each US climate zone defined in McPherson [59], six candidate reference cities (RCs) were selected by comparing heating degree days (HDDs), cooling degree days (CDDs), and annual precipitation with the subject city (SC), Kyoto; (2) the root mean square error (RMSE) of climatic variables between the RCs and SC was calculated using Equation (4), where a, b, and c are positive weighting coefficients that add to 1.0, expressing the relative importance of each variable, (3) one RC with the minimum RMSE was selected as the reference city for Kyoto, and (4) a climate region in which the selected RC falls into is selected from i-Tree Eco's 11 climate regions [9].
In the US, the thermal resistance (R-value) for each part (wall, ceiling, window, floor, and foundation) of the default housing for each vintage for each climate region are defined [9], whereas, in Japan, the heat loss coefficient (Q-value) for an entire house is defined as an energy conservation standard [60]. By integrating housing parts with R-value into the entire house, the Q-value for the default house was calculated, which, in turn, was compared to the standard in Japan to identify the match of the house vintages between the US and Japan.
Building data for Kyoto was obtained from the Geospatial Information Authority of Japan [39]. Residential houses were identified with a footprint of 100 m 2 or less, based on an average footprint of the area [61], and the direction and distance between the closest tree to each house were calculated using a geographical information system (GIS). Trees less than 6 m in height were excluded in this study because they are too low to affect the energy use in nearby houses. Houses located farther than 18 m from the nearest street tree were excluded as they were too far from the tree for their energy use to be affected.

Avoided Stormwater Runoff
i-Tree Eco estimates avoided stormwater runoff based on Hirabayashi [62], in which storm water runoff for two scenarios, (1) with the current tree cover and (2) with no tree cover in the study area, were calculated. The difference between the two scenarios was considered avoided stormwater runoff because of the existence of the trees. Within the model, hourly precipitation, rainfall intercepted by tree leaves determined based on LAI, evaporation from leaves, rainwater dropped to the ground, infiltration of the pervious cover, and runoff from the impervious cover were calculated based on Wang et al. [63]. One limitation here is that because the soil information is not available in i-Tree Eco, it was assumed that all the rainwater reaching the pervious cover infiltrates into the ground, while the rainwater reaching the impervious cover all runoffs. Average impervious and pervious covers for the study area were determined based on JAXA's ALOS-2 land cover data [49], and it was assumed that these values were uniform across the study area. The valuation for the avoided runoff was performed with the default value in i-Tree Eco, which is 2.36 US$/m 3 for stormwater control facilities in the US.

Size Distribution
The size distribution (in terms of diameter at breast height, DBH) is a key factor in managing a resilient tree population, influencing present and future costs as well as the flow of ecological benefits [64].
The DBH structure of street trees in Kyoto city was distributed unevenly comparing to the "ideal" size distribution proposed by [65], with a preponderance of maturing (15-30 cm DBH), mature (30-45 cm DBH) street trees, which account for 43.90% and 33.30%, respectively. The distribution indicated that many of Kyoto city's street trees were planted 20 to 50 years ago, and they provide maximum benefits because of their size and condition. There is a paucity of young, small-diameter classes (0-15 cm DBH), where the proportion is 17.8% lower than the ideal for offset establishing-related mortality. The species most heavily represented in the large, mature classes (>30 cm DBH) included Z. serrata (66.70%), A. buergerianum (49.50%), P. × acerifolia (51.00%), and S. babylonica (76.30%). Notably, P. × yedoensis (31.00%) was present in old tree classes (>60 cm DBH), offering extensive ecological services in Kyoto city (Figure 2). Diameter at breast height (DBH) size distribution of predominant street tree species compared to an "ideal" distribution [65]. Note: The "Kyoto" classification represents the total of the sampled street trees.  Diameter at breast height (DBH) size distribution of predominant street tree species compared to an "ideal" distribution [65]. Note: The "Kyoto" classification represents the total of the sampled street trees.  Kyoto's sample street trees intercepted approximately 1699.39 m 3 of rainfall annually, and this effect was associated with the benefits of stormwater runoff reduction at US$4011.76 (Table 3). On average, each street tree contributed to an intercept of 1.38 m 3 of rainfall annually and provided a value of US$3.26. Z. serrata (US$7.10/tree), P. × yedoensis (US$4.83/tree), and P. × acerifolia (US$4.42/tree) had the greatest effect on stormwater runoff reduction benefits. Table 4 shows the information about US reference counties surveyed for BenMAP calculations. Based on air quality improvement, the avoided incidence of adverse health effects was estimated to be 3.8 cases and the associated economic value of US$14,515.05 annually or US$11.80/tree. The greatest amount of removal occurred with the O 3 and NO 2 pollutants, while the greatest value associated with removal was for PM 2.5 and O 3 ( Table 5). Most of these monetary values were dominated by the effects of reducing human mortality because BenMAP assigns the greatest value per incidence for human mortality by averaging US$3.9 million per incidence [66].  The best match city for Kyoto was Charleston, SC (RMSE = 1.47). Thus, the "US Southeast climate region" was used to model the environmental conditions in Kyoto using i-Tree Eco.

Human Health Effects
House Vintage Adjustment Table 6 summarizes the R-value to Q-value conversion using Post-1980 vintage houses in the southeast climate region as an example. The inverse of the R-value is the heat transfer coefficient, U-value. The U-value is then multiplied by the area of interest including walls, ceilings, windows, floors, and foundations to obtain the heat loss value. The sum of the heat loss is divided by the floor area to obtain the Q-value, which is 2.45. Table 6. R-value and Q-value conversion results (post-1980 vintage for the southeast region).  -1950 and 1950-1980 were also calculated with the above-described procedure. Table 7 summarizes the energy conservation standards in the US and Japan. The Q-value for post-1996 construction in Japan was 2.7, which is equivalent to that of the post-1980 in the southeast climate region. Therefore, post-1996 houses in Japan were treated as post-1980 houses in the US. Higher Q-values indicate lower insulation capacity. The

Energy Saving
Among the 1230 sampled street trees, 614 trees affected the energy consumption of 1 or 2-story houses by shading buildings, providing evaporative cooling, and blocking winter winds. The annual energy-related costs were estimated at US$2054.36 using customized parameters (Table 8).

Annual Net Benefits and Costs Annual Benefits
The total annual value of sample street trees was evaluated by summing the four different estimated ES benefits, which were calculated at US$71,434.21 annually or US$58.07 per tree on average ( Table 9). The largest benefits were carbon storage and sequestration, which accounted for 71.19% of the total benefits. In contrast, energy savings contributed the least to ES, at only 2.88% of the total benefits (Table 9). On a per tree basis, P. × yedoensis (US$225.32/tree), Z. serrata (US$123.21/tree), S. babylonica (US$80.10/tree), and P. × acerifolia (US$65.88/tree) produced significant benefits, whereas C. florida (US$11.78/tree), L. tulipifera (US$33.64/tree), P. jamasakura (US$43.31/tree), and G. biloba (US$43.74/tree) produced the least benefits (Table 10).

Expenditures
Based on the information provided by the Green Policy Promotion Office of the Kyoto City Construction Bureau, the municipality of Kyoto spent exceed US$4,500,000 to maintain the population of 50,000 street trees annually including the costs of pruning, cleaning the fallen leaves, and pest and disease control. The average annual street tree expenditure is estimated at US$90 (90 = 4,500,000/50,000) (Kyoto City Construction Bureau, personal communication). Results are reported in US dollars, economic data collected in Japanese Yen were converted to US dollars using an exchange rate of 1US dollar to 110 Yen (Table 9).

Discussion
This study in Kyoto is the first attempt to use the empirical data in quantifying the ES of street trees, and the results could provide the municipality with baseline values for future management.
Our results indicate that large-growing species with more leaf area, such as P. × yedoensis (174.17 m 2 /tree, US$225.32/tree), Z. serrata (255.66 m 2 /tree, US$123.21/tree), and P. × acerifolia (159.36 m 2 /tree, US$65.88/tree) appeared to be the most valuable species in Kyoto city, which reflects the importance of the leaf area is the driving force behind the ability of trees to offer ecological benefits for the community. Many previous studies have pointed out that tree canopy cover and leaf area play a key role in determining the delivery of ES [5,19,66]. The greater the tree canopy size, the greater the pollution removal and precipitation intercept, and the greater the value provided [5,21,22,66]. Maintaining the health and longevity of these large trees is critical to achieving high ES levels in Kyoto city.

Benefits-Cost Comparison
Regarding the benefits-cost ratio, it was demonstrated that the value of the annual benefits generated by street trees in Kyoto city did not outweigh tree-related expenditure. This finding could be explained by three factors: first, the most dominant species, G. biloba (US$43.74/tree), which accounted for 47.8% of the entire population, also had a low tree canopy cover and therefore provided a low level of benefit. In Kyoto city, the average leaf area of G. biloba (74.34 m 2 /tree) (DBH=26.1 cm, height = 8.55 m) is considerably lower than that of the same species in other cities. Peper et al. [67] predicted the size of 12 common street trees growing in Modesto, CA, USA, and showed that the leaf area of G. biloba was 235.29 m 2 /tree (DBH = 38.7 cm, height = 11.74 m). The difference may be primarily due to the street tree management practices implemented in Kyoto city. The conflict between street tree expansion and constrained planting space is prominent in Japan, resulting heavy pruning, which suppresses the growth of tree canopy [68].
Additionally, on a per-tree basis, expenditure for street trees in Kyoto city (US$90) is the highest compared with cities in the US and Europe, such as New York City, US (US$37), Santa Monica, US (US$53), Lisbon, Portugal (US$46) [21,69,70]. Considering the need for frequent pruning, the likelihood of increases in the burden of municipal expenditure is very high in Japan [68].
Furthermore, the estimation of the benefits presented in this study represents only a fraction of the comprehensive value of Kyoto's street trees. In light of the experience from other case studies conducted with the i-Tree Street model, property value accounted for the most important benefits in all cities [70]. Trees contribute many "other" intangible benefits that are difficult to translate into economic terms, such as biodiversity, beautification, increased human comfort, and sense of place, which cannot yet be factored into the i-Tree Eco model. Kyoto City is a world-famous tourist destination, and street trees in the City are considered a significant attraction for tourists as well as a benefit for urban residents. In recent years, various levels of government in Kyoto have become increasingly aware of the importance of street trees and have made a renewed investment in the median strip planting project. Additional research is needed in this area to provide reliable information on factors such as property attributes, market responses, and willingness to pay.

Advantages
i-Tree Eco is composed of three elements: model codes, model parameters, and input data. One advantage of i-Tree Eco is that it uses model parameters and input data globally applicable or flexibly substitutable to local site-specific values when running it in the region outside the US [51]. Those parameters/data globally applicable are denoted with "a" in Table 1, including the social cost of carbon and worldwide weather and upper air data stored in i-Tree server computers. Model parameters/input data that an international user can replace via ITDB are location-related data, precipitation, air quality, and energy prices (denoted with "b" in Table 1).
In addition, with the cooperation of an i-Tree developer, the model codes and parameters (denoted with "c" in Table 1) were modified in this study. As a result, it has become possible to estimate the air pollution removal and the avoided stormwater runoff by employing solar and net radiations measured in the vicinity of Kyoto rather than estimated values employed in i-Tree Eco.
Regarding the health effects, the analysis available through ITDB for typical international users is very limited. A regression equation for each of the four air pollutants (CO, NO 2 , O 3 , and PM 2.5 ) relating population density and monetary values avoided due to air quality improvement was constructed based on the county-based i-Tree Eco runs across the continental US [66]. What typical international i-Tree Eco users can do is to plug their population density into these regression equations to estimate the monetary values. It is impossible for them to calculate the avoided incidence and monetary values for each of the adverse health endpoints. By contrast, our study enabled us to quantify these detailed amounts by referencing counties in the US and adjusting their parameters for Kyoto.
Look-up tables utilized in i-Tree Eco to estimate the tree's effects on heating and cooling energy savings at houses across the US were developed based on lots of data collections, modeling, and analyses. These efforts were made possible thanks to funding and collaborations from national laboratories, federal agencies, non-profit organizations, local governments, and utility companies [9]. They also provided instruction on how to determine a reference US climate region for international applications of the method [59], and we followed it. Additionally, we developed a method to replace the US house vintages with those in Japan based on the energy conservation standards in the two countries. This is beyond what common international users can do using ITDB. The method we used in the study may be the best for now when using the look-up tables contained in i-Tree Eco.

Limitations and Future Directions
Despite the advantages brought by customizing i-Tree Eco's model codes, parameters, and input data, there still exists several limitations and uncertainties in the assumptions and processes conducted in this study, which in turn affect the validity of the results.
The modeled carbon values are estimates based on tree growth allometric equations from the US, and the carbon estimate error includes the uncertainty of using biomass equations and conversion factors [51]. It is well known that the growth characteristics of street trees are greatly affected by different management practices and climate conditions. An estimation uncertainty was present in the international case study using surrogate US species data. Currently, there are limited studies that have compiled growth equations for street trees in Japan [71]. Future research is needed to develop growth curves to understand the biomass equations for urban trees in Kyoto city with greater accuracy.
With regards to the health effects assessed, adverse health incidences and associated monetary values reduced because of the air quality improvement were estimated based on BanMAP, assuming that the response of humans to air pollution is the same in Japan and in the US. Although the process taken in this study to adjust the health effects in the reference US counties for Kyoto may be plausible when the established methods like BenMAP lack in Japan, there exist uncertainties in the assumption. For instance, susceptibility to air pollution may be affected by many factors such as genetic background, race, ethnicity, and culture [72] to name a few, which vary between Japan and the US. One future research direction is to explore epidemiologic data and develop methods in Japan that enable analyses like BenMAP.
In the process of the heating and cooling energy saving calculations, the matching of the climate region and the house vintages between Japan and the US were not perfect. Hence, uncertainties exist there. It is desirable to develop a similar means to quantify the trees' effects on household energy savings based on the climate and house characteristics in Japan.
Lin et al. [73] performed thorough sensitivity analyses on several components of i-Tree Eco (i.e., carbon storage/sequestration, bio-emissions, and dry deposition of air pollution) to identify important input variables for each analysis. It is suggested that increasing the accuracy of these important variables is an effective way to reduce uncertainty in the model output. Unfortunately, energy savings, health effects, and stormwater reductions that we assessed in this study were not included in these analyses. Therefore, it will be a great addition to the i-Tree research and user communities if the sensitivity analyses on these three components are performed in the future.

Conclusions
This study describes the first sample tree inventory data-based street tree ES assessment in Kyoto, Japan, by customizing the model and parameters of the i-Tree Eco model. The results presented in this study should be considered first-order estimations of the ES, since they were unable to validate against the ground truth due to a lack of such data. Despite that, treating our results as a reference value, they contribute new knowledge on the structure, function, and value of Kyoto's street trees.
Street tree survival, growth, and management in Kyoto city pose a unique set of problems because the majority of street trees are growing in a stressful urban environment that has been impacted and constrained by construction for many years. To maintain the flow of benefits the city currently enjoys, management recommendations derived from this analysis are as follows: Continue investing in intensive maintenance of large-stature mature trees to prolong the lifespans of tree species such as P. × yedoensis, Z. serrata, and P. × acerifolia.
To maximize the trees' potential for reducing energy consumption and ensure longterm net benefits from continuous levels of tree canopy cover, heavy pruning should be discontinued, and planting strips should be advocated in new street tree plans.
It is recommended that diversification be continued to reduce dependence on species such as G. biloba, while concentrating on selecting tree species that can tolerate restricted site conditions, avoiding unnecessary pruning and management costs.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.