Growth of Abies sachalinensis Along an Urban Gradient Affected by Environmental Pollution in Sapporo, Japan

Urban tree growth is often affected by reduced water availability, higher temperatures, small and compacted planting pits, as well as high nutrient and pollution inputs. Despite these hindering growth conditions, recent studies found a surprisingly better growth of urban trees compared to trees at rural sites, and an enhanced growth of trees in recent times. We compared urban versus rural growing Sakhalin fir (Abies sachalinensis (F. Schmidt) Mast.) trees in Sapporo, northern Japan and analyzed the growth differences between growing sites and the effects of environmental pollution (NO2, NOX, SO2 and OX) on tree growth. Tree growth was assessed by a dendrochronological study across a gradient from urban to rural sites and related to high detailed environmental pollution data with mixed model approaches and regression analyses. A higher growth of urban trees compared to rural trees was found, along with an overall accelerated growth rate of A. sachalinensis trees over time. Moreover, environmental pollution seems to positively affect tree growth, though with the exception of oxides OX which had strong negative correlations with growth. In conclusion, higher temperatures, changed soil nutrient status, higher risks of water-logging, increased oxide concentrations, as well as higher age negatively affected the growth of rural trees. The future growth of urban A. sachalinensis will provide more insights as to whether the results were induced by environmental pollution and climate or biased on a higher age of rural trees. Nevertheless, the results clearly indicate that environmental pollution, especially in terms of NO2 and NOX poses no threat to urban tree growth in Sapporo.


Introduction
The effects of air pollution and climate change on tree growth have been discussed ambiguously over the past decades. While several studies clearly link tree and forest damage with sensitivity to air and environmental pollution [1][2][3][4][5][6], some studies also named unfavorable climatic conditions with limited soil water availability [4,7], aggravated soil compaction, nutrient imbalances [8,9], and pests and disease infestation, as well as management errors as reasons for tree growth decline in forests [1,3,10,11]. On the other hand, several studies [12][13][14][15][16] found opposing positive effects of climate change conditions and increased environmental pollution on tree growth in forests. Presumed causes for the improved growth are higher nitrogen depositions, higher temperatures and higher CO 2 -concentrations, together with a longer growing season and changes in forest management [12,17]. However, limited water of a coniferous urban tree species growing in the boreal climate zone. This study focuses on the growth of urban and rural A. sachalinensis MAST. trees in Sapporo, northern Japan, and the effects of urbanization, climate change and environmental pollution. This town rapidly increased its population from about one to two million in the past five decades. The growth of urban and rural trees was assessed by dendrochronology and related to climate and environmental pollution. The following research questions were stated: 1.
What was the growth of A. sachalinensis in the past decades in the urban and rural areas of Sapporo, northern Japan? 2.
Are the growth trends of A. sachalinensis similar to worldwide trends of urban tree growth? 3.
Can differences in the growth of A. sachalinensis be found regarding the sampling sites?

Climate of Sapporo
Sapporo has about 2.0 million inhabitants and a size of 1121.12 km 2 , the biggest city of Hokkaido island, northern Japan (43 • 4 N, 141 • 21 E~43 • 3 43 N, 141 • 21 • 15 E). Its climate is characterized as cold, without a dry season and hot summers [39,40] with an average annual temperature of 7.8 • C and a precipitation sum of 1130 mm (mean of 1980-2012). Over the year, the highest temperatures occur in July, August and September, with August as the hottest month (average of 20 • C). The coldest months are January and February with below 0 • C. The highest amount of precipitation also occurs in August and September (Figure 1, Japan Meteorological Agency). of urban and rural A. sachalinensis MAST. trees in Sapporo, northern Japan, and the effects of urbanization, climate change and environmental pollution. This town rapidly increased its population from about one to two million in the past five decades. The growth of urban and rural trees was assessed by dendrochronology and related to climate and environmental pollution. The following research questions were stated: 1. What was the growth of A. sachalinensis in the past decades in the urban and rural areas of Sapporo, northern Japan? 2. Are the growth trends of A. sachalinensis similar to worldwide trends of urban tree growth? 3. Can differences in the growth of A. sachalinensis be found regarding the sampling sites?

Climate of Sapporo
Sapporo has about 2.0 million inhabitants and a size of 1121.12 km 2 , the biggest city of Hokkaido island, northern Japan (43°4' N, 141°21' E ~ 43°3'43' N, 141°21°15'E). Its climate is characterized as cold, without a dry season and hot summers [39,40] with an average annual temperature of 7.8 °C and a precipitation sum of 1130 mm (mean of 1980-2012). Over the year, the highest temperatures occur in July, August and September, with August as the hottest month (average of 20 °C). The coldest months are January and February with below 0 °C. The highest amount of precipitation also occurs in August and September (Figure 1, Japan Meteorological Agency). To further analyze the effect of climate on tree growth, we calculated the de Martonne-index [43] on the basis of monthly precipitation and monthly temperature from 1983 to 2012. The snow-free period is from mid-April to early November.

Environmental Pollution Data
Environmental pollution data (nitrogen dioxide NO2, oxides OX, sulfur dioxide SO2 and nitric oxide NOX) in Sapporo have been measured at several stations across the city ( Table 1). Elevation of each monitoring site is about 50-200 m a.s.l. The highest concentrations of NO2 were found at Kita-1jyo and Tsukisamu-South, with the heavily trafficked region having a of mean concentration of 20.8 ppm. The highest OX concentrations were found in Yamahana and Atsubetsu-the SE located suburb or close to reserved forests. The average concentrations of all stations was 26.9 ppm. For SO2 we found mean concentrations of 4.4 ppm (highest concentrations at stations Sapporo Middle Part and West To further analyze the effect of climate on tree growth, we calculated the de Martonne-index [43] on the basis of monthly precipitation and monthly temperature from 1983 to 2012. The snow-free period is from mid-April to early November.

Environmental Pollution Data
Environmental pollution data (nitrogen dioxide NO 2 , oxides O X , sulfur dioxide SO 2 and nitric oxide NO X ) in Sapporo have been measured at several stations across the city ( Table 1). Elevation of each monitoring site is about 50-200 m a.s.l. The highest concentrations of NO 2 were found at Kita-1-jyo and Tsukisamu-South, with the heavily trafficked region having a of mean concentration of 20.8 ppm. The highest O X concentrations were found in Yamahana and Atsubetsu-the SE located suburb or close to reserved forests. The average concentrations of all stations was 26.9 ppm. For SO 2 we found mean concentrations of 4.4 ppm (highest concentrations at stations Sapporo Middle Part and West District with heavy traffic), while for NO X the overall mean concentration was about 39.9 ppm (highest concentrations at stations Tsukisamu-South, near the reserve forest, and Kita-1-jyo, with heavy traffic). The average measured temperature of all stations across Sapporo was 8.6 • C. The warmest temperatures occurred at southeast Atsubetsu, in the east of Sapporo and at Yamahana, in the central south of Sapporo. Further, carbon monoxide CO was measured at the station Kita-1-jyo from 1986 to 2012, the mean value was 0.9 ppm. The values of NO 2 and O X vary evenly around the mean, while SO 2 , NO X and temperature data are dominated by outliers, possibly induced by measurement errors or extreme conditions; e.g., at main roads. Table 1. Environmental pollution (nitrogen dioxide: NO 2 ; oxides: O X ; sulfur dioxide: SO 2 ; and nitrous gases: NO X ) and climate data of several weather stations across Sapporo, Japan from 1983 to 2012, and variance from the mean value (+ higher than mean, -lower than mean). Added * to station name means suburb or closest to a green area (e.g., near an agriculture field, forested park or reserved forests).  The change of the environmental pollution values in Sapporo over six time periods ranging from 1983 to 2012 is displayed in Figure 2. Although NO X , NO 2 and SO 2 concentrations were decreasing, the O X concentrations were increasing. The regression line showed a high coefficient of determination for NO X , SO 2 and O X (R 2 > 0.65) over time; however, for NO 2 it was lower (R 2 < 0.25). For NO X and NO 2 , there was a peak concentration recorded in the period 1998 to 2002. The overall highest concentrations of environmental pollution were found for NO X , with values up to 50 ppm. The lowest concentrations were observed for SO 2 .  (1983-1987, 1988-1992, 1993-1997, 1998-2002, 2003-2007 and 2008-2012) of several weather stations across the city of Sapporo. Given are the best fitting regression lines with regression coefficients and coefficients of determination R 2 .

Sample Tree Species
The study focused on Sakhalin fir, A. sachalinensis, an evergreen conifer species originating of the Sakhalin islands and southern Kurils, Russia. The species also occurs in northern Hokkaido, Japan. It prefers moist climates with cool summers and mild winters, though it faces problems if exposed to waterlogged soils [44,45]. The shade tolerance of A. sachalinensis is very high and the growing rate low [46]. A screening experiment for 18 species native to Japan by Yamaguchi, et al. [47] showed, that most Abies species native to Japan are classified into the intermediate ozone sensitivity type (responses to AOT40; 16-30 ppm h), A. sachalinensis prefers slightly acidic soils with pHs around 5 [48].

Data Collection
Across the city of Sapporo, 109 A. sachalinensis individual trees were chosen for data sampling at six sites along a gradient from the city center to the suburbs of Sapporo and a forest area outside of the city (Figure 3).  (1983-1987, 1988-1992, 1993-1997, 1998-2002, 2003-2007 and 2008-2012) of several weather stations across the city of Sapporo. Given are the best fitting regression lines with regression coefficients and coefficients of determination R 2 .

Sample Tree Species
The study focused on Sakhalin fir, A. sachalinensis, an evergreen conifer species originating of the Sakhalin islands and southern Kurils, Russia. The species also occurs in northern Hokkaido, Japan. It prefers moist climates with cool summers and mild winters, though it faces problems if exposed to waterlogged soils [44,45]. The shade tolerance of A. sachalinensis is very high and the growing rate low [46]. A screening experiment for 18 species native to Japan by Yamaguchi, et al. [47] showed, that most Abies species native to Japan are classified into the intermediate ozone sensitivity type (responses to AOT40; 16-30 ppm h), A. sachalinensis prefers slightly acidic soils with pHs around 5 [48].

Data Collection
Across the city of Sapporo, 109 A. sachalinensis individual trees were chosen for data sampling at six sites along a gradient from the city center to the suburbs of Sapporo and a forest area outside of the city ( Figure 3).
The trees in sample plot 1 and plot 2 (Shirahata-yama 1 and Shirahata-yama 2) together with plot 6 (Misumai of Hokkaido University Forests) were all classified as rural ( Table 2). The trees of plot 3 (Hokkaido University Nursery in Sapporo) were growing along street canyons and were therefore classified as urban. The trees of sample plots 4 and 5 (both Hitsujigaoka-7 site) were classified as suburban trees; however, for model development, were merged with the urban trees. The sample trees at the urban and suburban plots were typically trees growing in cities, standing along streets, in front of buildings and at squares. The trees at the forest sites were party planted, mostly due to esthetic reasons and not for timber production, since the wood of A. sachalinensis gains low prices in wood markets. The stem density at the rural plots spans from 2000 (plot 2), over 2200 (plot 6), to 2500 (plot 1) stems per hectare. To the best of our knowledge, the trees at plots 3 (urban), and 2 to 4 (rural to suburban) were not further managed. Trees of plots 1 (rural), as well as of plots 5 to 6 (suburban to rural) faced typhoon events or thinning. The soil nutrient status was slightly different across the sample plots. While urban and suburban plots had higher nitrogen (N) content and bulk densities than rural plots, the calcium (Ca) and magnesium (Mg) content was higher at rural plots.
The phosphorus (P) and potassium (K) content was not consistent over plot classification; however, highest concentrations were found for P at Shirahata-yama 2 and for K at Hitsujigaoka-7 2nd plots. The trees in sample plot 1 and plot 2 (Shirahata-yama 1 and Shirahata-yama 2) together with plot 6 (Misumai of Hokkaido University Forests) were all classified as rural ( Table 2). The trees of plot 3 (Hokkaido University Nursery in Sapporo) were growing along street canyons and were therefore classified as urban. The trees of sample plots 4 and 5 (both Hitsujigaoka-7 site) were classified as suburban trees; however, for model development, were merged with the urban trees. The sample trees at the urban and suburban plots were typically trees growing in cities, standing along streets, in front of buildings and at squares. The trees at the forest sites were party planted, mostly due to esthetic reasons and not for timber production, since the wood of A. sachalinensis gains low prices in wood markets. The stem density at the rural plots spans from 2000 (plot 2), over 2200 (plot 6), to 2500 (plot 1) stems per hectare. To the best of our knowledge, the trees at plots 3 (urban), and 2 to 4 (rural to suburban) were not further managed. Trees of plots 1 (rural), as well as of plots 5 to 6 (suburban to rural) faced typhoon events or thinning. The soil nutrient status was slightly different across the sample plots. While urban and suburban plots had higher nitrogen (N) content and bulk densities than rural plots, the calcium (Ca) and magnesium (Mg) content was higher at rural plots. The phosphorus (P) and potassium (K) content was not consistent over plot classification; however, highest concentrations were found for P at Shirahata-yama 2 and for K at Hitsujigaoka-7 2 nd plots.  Before increment core collection, data on tree structure and the site conditions were recorded, including diameter at breast height, 1.3 m (dbh), tree height (h), height to the crown base (cb), crown radius in four directions (N, E, S, W), tree position (coordinates and altitude), site condition, tree vitality, and open surface area of the unpaved area around the tree in four directions. Based on the recorded data, the average crown radius cr (Equation (1)) and crown diameter cd (cd = cr × 2), the crown projection area (CPA) (Equation (2)) and the crown volume (cv) (Equation (3)) were calculated.
The increment core collection was conducted at each tree. Two cores opposite to each other were extracted at a height of 1.3 m in northern and western directions with a 5 mm diameter increment corer (Haglöf Sweden AB; Långsele, Västernorrland, Sweden).

Core and Data Processing
The cores were processed by mounting on wooden boards with regard to the grain direction. Thereafter, the cores were sanded until the highest visibility of the cross-sectional area, and then polished with progressively finer sandpaper from grit size 180 up to 800. Annual tree-ring widths of the cores were measured with a digital positioning table with a resolution of 1/100 mm (Rinntech e.K., Heidelberg, Germany). Crossdating and synchronization of ring-width data were accomplished by the software TsapWin (Rinntech e.K., Heidelberg, Germany) using standard dendrochronological methods [49][50][51].
All following analyses were conducted with the package dplR of R [52]. The biological age trend (higher growth of younger trees) in the ring-width data was removed by a double detrending procedure applied to all series (modified negative exponential curves and cubic smoothing splines with 20 year rigidity, 50% wavelength cutoff). The resulting index series contained only year-to-year variability associated with fluctuations in climate [49,50] producing dimensionless ring-width indices (RWIs). In a final step, autocorrelation was removed by autoregressive (AR) models (maximum order of three) and the series were averaged using Tukey's biweight robust mean. This reduces bias caused by extreme values. Mean sensitivity was calculated as assessment of chronology quality; it depends on the year-to-year variability and was employed as a measure of variability. All further analyses of climate-growth correlations were conducted with the resulting chronologies. From the chronologies, the ages of the analyzed trees were derived. When the exact age of the tree was not clear (missing tree pith, among others), the age was back calculated based on the un-detrended average growth rate of the last ten years and the dbh of the tree.

Statistics
Data on previously measured tree structures (dbh, height and crown values) were tested for significant differences between groups of urbanity (urban versus rural versus suburban and urban versus rural). Since the assumptions of normality and homogeneity were not met, the Kruskal-Wallis test with pairwise testing for significance (pairwise Wilcox-test with Bonferroni-Holm p-value correction) was applied for three groups and the Kruskal-Wallis test for two groups. Statistics were done in R, version 3.3.3 [53].

Trend Analysis (Long-Term Trends)
Using the R package lme4, two linear mixed models of the following forms Equations (4) and (5) were developed to assess the influence of the time of age, growth (before 1960 and since 1960) and urbanity (urban-rural) on the annual basal area (response variable) derived from increment cores. To differentiate between the two growths-trend relevant periods (before 1960 and since 1960), we introduced the dummy variable recent, with 1 indicating each observation later than 1959 and 0 otherwise, which is in accordance with growth trends, since approximately the 1960s, that have been identified for forest trees.
In Equations (4) and (5) the basal area is the response variable for the jth of n i observations in the ith of M groups or clusters, and a 1 , . . . , a n and b 1 , . . . , b n are the fixed effects with the a parameter's components of the intercept and the b parameter's components of the slope. When a 1 in Equation (4) differed significantly from 0, the age-basal area relationship before 1960 had a different intercept than that since 1960. In Equation (5), any differences in the intercepts would indicate that the intercept of urban trees was not the same as that for rural trees. The parameter b 1 in both equations have an analogous meaning to that of the slope. The c parameters are random effects, which are assumed to have a normal distribution. These random effects cover statistical dependencies, which are due to the nested data structure. The errors ε ij are assumed to have an independent, identical distribution. With both models, we try to answer the research questions 1, 2 and 3 to derive the growth trend of A. sachalinensis in Sapporo in view of time and location.

Results
For the 109 measured A. sachalinensis trees, an average dbh of 34.4 cm with a mean age of 59 years was found. The mean tree height was 17.8 m and the average cd was 5.9 m. Moreover, an average cv of 1621 m 3 was found. Table 3 displays the measured and calculated tree structures of A. sachalinensis in Sapporo classified by the sampling sites. The trees at rural plot 2 (Sharahata-yama 2) and rural plot 6 (Misumai) were the oldest with an average age of around 100 years. The according tree characteristics dbh, cd, cpa, crown length and cv were greatest at both plots as well. Youngest trees were found at suburban plots 4 and 5, Hitsujigaoka 1st and Hitsujigaoka 2nd, both in the Hitsukigaoka-7 site, with a mean age of around 35 years. Tree on those plots were, therefore, the smallest. Table 3. Number of sampled trees per plot n with urbanity classification (urban, suburban, rural) and tree characteristics (age, dbh, tree height, crown start, crown diameter, crown projection area CPA, crown length and crown volume) of A. sachalinensis in Sapporo, Japan. A further classification of A. sachalinensis to urban, suburban and rural growing trees illustrated that greatest tree structures were mostly found for rural trees (with exception of crown start). Urban and rural tree structures were significantly different, with suburban trees showing intermediate size (crown start), or were similar to the urban trees (dbh, tree height, cd, cpa and cv). Only the crown lengths of suburban trees were more similar to those of rural trees (Table 4).

Short-Term Growth Trends
When analyzing the overall growth trends of A. sachalinensis in Sapporo, a strong age-trend could not be found (Figure 4a). At younger ages, A. sachalinensis did not show a better radius growth than older trees displayed. However, with the exception of 2011, a decrease in growth can be observed during the last ten years. After double-detrending, these trends cannot be found in Figure 4b   When looking in more detail at the growth of A. sachalinensis in Sapporo, analyzing the growth chronologies of trees at urban, suburban and rural settings, distinct trends can be found ( Figure 5). The raw radial growth curves of rural trees displayed a far higher age of rural trees than of urban and suburban trees (Figure 5a). Urban and suburban chronologies started the earliest, in 1985. When regarding their growth, rural trees showed lowest stem growth rates, while urban and suburban trees indicated clearly higher values. For all chronologies, a decline was found for the last ten years, which was most pronounced for suburban and urban trees. After the detrending process and after the agetrend removal (Figure 5b), the growth chronologies appear more similar. However, common positive  When looking in more detail at the growth of A. sachalinensis in Sapporo, analyzing the growth chronologies of trees at urban, suburban and rural settings, distinct trends can be found ( Figure 5). The raw radial growth curves of rural trees displayed a far higher age of rural trees than of urban and suburban trees (Figure 5a). Urban and suburban chronologies started the earliest, in 1985. When regarding their growth, rural trees showed lowest stem growth rates, while urban and suburban trees indicated clearly higher values. For all chronologies, a decline was found for the last ten years, which was most pronounced for suburban and urban trees. After the detrending process and after the age-trend removal (Figure 5b), the growth chronologies appear more similar. However, common positive or negative pointer years for all three chronologies could not be found. The decline in the past years was observed as well, albeit a steep increase in growth was found in 2011. The basal area increment (Figure 5c) indicates a higher level of BAI of rural trees due to their higher age and longer growing period. Both urban and suburban trees had a marked increase in BAI since their planting, exceeding rural trees between 2005 and 2008. The statistical characteristics of the displayed chronologies can be found in Table 5. While the highest growth rate was found for urban and suburban trees, rural trees were growing on average the least. Mean sensitivity was 0.23, with highest for the urban chronology and lowest for the suburban chronology (0.19). Autocorrelation was found highest for rural trees, at 0.68 and least for urban trees, at 0.33.  The statistical characteristics of the displayed chronologies can be found in Table 5. While the highest growth rate was found for urban and suburban trees, rural trees were growing on average the least. Mean sensitivity was 0.23, with highest for the urban chronology and lowest for the suburban chronology (0.19). Autocorrelation was found highest for rural trees, at 0.68 and least for urban trees, at 0.33.

Long-Term Growth Trends
Tree growth in Sapporo has shown an enhanced growth for the past decades compared with the growth before 1960. This is in accordance with the overall global trend of urban trees found by [17] ( Figure 6a). Young A. sachalinensis grew similarly, regardless of their planting time; however, after reaching an age of around 50 years, trees planted in more recent decades showed an enhanced growth, also passing the overall worldwide average. Urban A. sachalinensis trees illustrated a far better growth compared to rural trees, in particular after an age of 60 years, urban tree growth exceeded rural trees' growth by far (Figure 6b). In comparison to the worldwide results on urban and rural tree growth by Pretzsch, Biber, Uhl, Dahlhausen, Schütze, Perkins, Rötzer, Caldentey, Koike, van Con, du Toit, Foster and Lefer [17], A. sachalinensis rural trees showed a lower average growth rate. Younger trees were growing less than the worldwide average, though, at the old age of more than 100 years, A. sachalinensis exceeded the growth rate of other studies trees. Urban trees had a similar growth rate to those trees measured worldwide until an age of 50 years, and then an extremely high growing rate was found for A. sachalinensis in Sapporo. reaching an age of around 50 years, trees planted in more recent decades showed an enhanced growth, also passing the overall worldwide average. Urban A. sachalinensis trees illustrated a far better growth compared to rural trees, in particular after an age of 60 years, urban tree growth exceeded rural trees' growth by far (Figure 6b). In comparison to the worldwide results on urban and rural tree growth by Pretzsch, Biber, Uhl, Dahlhausen, Schütze, Perkins, Rötzer, Caldentey, Koike, van Con, du Toit, Foster and Lefer [17], A. sachalinensis rural trees showed a lower average growth rate. Younger trees were growing less than the worldwide average, though, at the old age of more than 100 years, A. sachalinensis exceeded the growth rate of other studies trees. Urban trees had a similar growth rate to those trees measured worldwide until an age of 50 years, and then an extremely high growing rate was found for A. sachalinensis in Sapporo. The described trends of Figure 6a,b can also be found by using a mixed model (Table 6a,b). While age had an overall significantly positive effect on growth of A. sachalinensis, the time of growth (before or since 1960) acted as negative driver, though it was not significant. In combination with age, a significant, positive interaction was found again. Therefore, the older trees were, the better was the tree growth since 1960. A similarly positive effect for age on tree growth as a single term was found when looking at the influence of urbanity. Urbanity had a highly positive influence on growth, indicating a better growth of urban A. sachalinensis compared to rural trees. In combination, however, urban trees showed a more enhanced growth with higher age than rural trees. Table 6a. Results of the linear mixed model on the annual basal area increment (mm 2 yr−1) of all analyzed trees (response variable) with the individual tree coded as random effect, and the fixed effects were log (age), time of growth (period before 1960 and period since 1960) and their interaction. The described trends of Figure 6a,b can also be found by using a mixed model (Tables 6 and 7). While age had an overall significantly positive effect on growth of A. sachalinensis, the time of growth (before or since 1960) acted as negative driver, though it was not significant. In combination with age, a significant, positive interaction was found again. Therefore, the older trees were, the better was the tree growth since 1960. A similarly positive effect for age on tree growth as a single term was found when looking at the influence of urbanity. Urbanity had a highly positive influence on growth, indicating a better growth of urban A. sachalinensis compared to rural trees. In combination, however, urban trees showed a more enhanced growth with higher age than rural trees. Table 6. Results of the linear mixed model on the annual basal area increment (mm 2 yr −1 ) of all analyzed trees (response variable) with the individual tree coded as random effect, and the fixed effects were log (age), time of growth (period before 1960 and period since 1960) and their interaction.

Parameter
Fixed

Relationships of Growth with Environmental Pollution
The growth was also related to the measured environmental pollution and climate data of Sapporo (Figure 7).

Relationships of Growth with Environmental Pollution
The growth was also related to the measured environmental pollution and climate data of Sapporo (Figure 7). Radius increment, concentrations of environmental pollutions and climate data were averaged in six periods of five years. Over the analyzed periods, strong correlations of growth with NOX, SO2, OX and CO were found (R 2 > 0.8). However, weak relationships were found between growth and NO2, temperature, precipitation and the de Martonne-index (R 2 < 0.4). When looking at the strengths of the relationships, a strong positive effect of higher environmental pollution on tree growth in Sapporo Radius increment, concentrations of environmental pollutions and climate data were averaged in six periods of five years. Over the analyzed periods, strong correlations of growth with NO X , SO 2 , O X and CO were found (R 2 > 0.8). However, weak relationships were found between growth and NO 2 , temperature, precipitation and the de Martonne-index (R 2 < 0.4). When looking at the strengths of the relationships, a strong positive effect of higher environmental pollution on tree growth in Sapporo was derived for NO X , NO 2 , SO 2 and CO. An increase in concentrations of these gases was strongly related to an increase in growth. On the other hand, for the O X concentration and the climate data temperature, precipitation and the de Martonne-index, a negative relationship with growth was observed. A higher concentration of O X as well as higher temperatures, more precipitation and a higher de Martonne-index was correlated with decreased growth.

Discussion
The studied A. sachalinensis trees in Sapporo, Japan were, across all sampling sites, an average age of 59 years old. This age is comparatively high for an urban tree species, since studies such as Nowak, et al. [54] and Skiera and Moll [55], stated a very low life span of 15 years on average and a high mortality rate of 6.6% per year. However, the most frequently planted trees in cities are deciduous tree species, illustrating a very high annual mortality rate of 12.6% for Acer negundo L., 7.4% for Robinia pseudoacacia L. and 6.8% for Fraxinus pennsylvanica Marsh., compared to results for coniferous species (3.2% for Picea abies (L.) Karst. and 0% for Pinus strobus L.) [54]. In a Japanese forest stand, a competition driven mortality rate of 1.7% per year was found for A. sachalinensis, with an overall age range of 33 to 187 years and a mean age of 104 years in the Nakagawa Experiment Forest of Hokkaido University, located in the northern part of the island [56]. Moreover, a maximum age of 249 years for canopy trees of A. sachalinensis was derived [56]. This is in line with the maximum age we found in our study, as a maximum age of 250 years was found for trees sampled at the rural Misumai plot. The oldest trees were found at the rural plots Shirahata-yama 2 and Misumai. The trees at suburban and urban plots were markedly younger. These plots were in the rural outskirts of the city and have been established earlier than the other urban and suburban sites. The average age of all sampled trees across Sapporo was younger compared to the study of Hiura, Sano and Konno [56]; however, when excluding the fairly young urban and suburban trees, an average age of 82 years for the trees at the forested sites can be derived. When only looking at plots Shirahata-yama 2 and Misumai, the average age of the trees was 100 years, which is very similar to the study about A. sachalinensis growing in a Japanese forest [56]. All in all, the found growth rates and ages seem to be in line with the results of other studies on coniferous species, urban trees and A. sachalinensis.
Tree age is correlated with tree structures like diameter at breast height (dbh), tree height and crown parameters, such as diameter, height, projection area and volume [57], and the older trees at rural plots Shirahata-yama 2 and Misumai had significantly greater dbh, height and other crown parameters than the trees growing at suburban and urban plots. Due to their young age, though, urban and suburban trees show a high growth rate. While the trees sampled at the forest-sites are more than double as old as the trees at suburban and urban plots, the difference in dbh is only 14 cm and 12 cm respectively. Overall, the height growth of A. sachalinensis in Sapporo is faster compared to trees measured in a northern Japanese forest stand [56]; however, urban trees often can reach greater tree sizes than forest trees due to more unobstructed growth and less competition [58].
The average raw radius growth rate of A. sachalinensis in Sapporo, northern Japan is with 4.2 mm/year, in line with other studies analyzing the growth of this species, such as Hiura, Sano and Konno [56], who found a growth rate of 1.0-7.0 mm/year, and Umeki [59] describing 3.1 mm per year, both studying A. sachalinensis at forest-sites in Japan. In urban areas, the growth rates of coniferous species have seldom been measured. Peper, et al. [60] provided data on the dbh increments between the ages 15 and 30 years for three coniferous species in Santa Monica, California, US. The transformed data indicate radius increments of 5.7 mm per year for Pinus canariensis C. Sm., 4.7 mm for Ficus microcarpa L. f. and 8.4 mm for Cedrus deodara (D. Don.) G. Don., which all are higher than those found for A. sachalinensis in Japan. When accounting for the similar age, the growth rate F. microcarpa is comparable to the radius increments of urban and suburban trees of A. sachalinensis. The highest growth was found for urban and suburban A. sachalinensis trees, while trees at rural plots grew least. This might be induced by the younger age of urban and suburban trees.
For all further analysis, the age difference was accounted for in model development by including the age of the trees or by an age detrending of the raw tree ring widths. A similar year-to-year variability (mean sensitivity) was found for the urban, suburban and rural chronologies, albeit trees at the urban sites showed higher mean sensitivity. This seems logical, since trees growing in extreme environments often display stronger sensitivity values than trees of sites with less stress and better water supply [61]. The higher autocorrelation of rural trees indicated, that radial growth was strongly influenced by the conditions of the preceding year, due to preformed shoot primordia in this species [62]. The lower autocorrelation of urban trees might be caused by the different conditions at the urban sites compared to forest sites. The urban environment seems to affect trees positively, while rural trees are more often faced with stressful conditions, such as possibly less nutrient inputs and higher ozone concentrations over 60~80 ppb (nmol mo −1 ). Information on the soil nutrient status showed, that the soil at rural plots has higher phosphorus and magnesium concentrations; however, the nitrogen content is lower. Moreover, the soil compaction is slightly less at rural sites.
The radius and BAI chronology of urban trees also indicated a growth decline in the past decades, though the steep increase in growth during the past year before sample taking. In total, the chronology of the urban trees showed a decreased variability and a higher stability over the past decades. The chronologies of urban and suburban trees illustrated high radius growth rates and strong increases in BAI. Though urban and suburban chronologies highlighted an opposing pattern to the rural trees in past decades (decrease in radius increment and BAI), this is in particular visible for the suburban chronology.
Trees of A. sachalinensis in Sapporo showed an enhanced growth rate during the most recent decades compared to former times (before 1960), as well as a better growth of urban trees compared to rural trees. These findings are in line with the trends found for R. pseudoacacia in Santiago de Chile [63] and Quercus nigra Mill. in Houston, Texas [42]. However, results of A. sachalinensis might be biased by the strong age and growth differences of urban and rural trees, since no sampled urban tree was older than 47 years (planted around [1970][1971][1972][1973][1974][1975]. Trees of the category "before 1960" were only rural trees, which illustrated a lower growth rate than urban trees. The different age of the trees was included in the model performance, though the great difference might still bias the results slightly. In total, the results of tree growth in Sapporo follow the overall worldwide trend [17] found for trees in metropolises, which describes an enhanced growth of urban trees, possibly induced by a longer growing period, higher CO 2 -concentrations and warmer temperatures compared to rural sites [15,41,64]. The accelerated growth of urban trees was most pronounced for the boreal zone. Urban trees of this zone, possibly gain most by higher temperatures, by higher nutrient inputs and by prolonged growing periods, especially when further water availability is sufficient [17]. An extended growing period by climate change, though, will not affect A. sachalinensis, as much as it might affect other deciduous species. Since A. sachalinensis is a coniferous species, it can use photosynthesis throughout the whole year, and much longer than deciduous species anyway. Guardans [2], on the other hand, described an increased climate change and environmental pollution sensitivity of trees in the boreal area due to temperature and water stress. This was also found by Rötzer, et al. [65] for urban trees in temperate cities. Due to the relative high precipitation amounts in Sapporo with more than 1000 mm per year, though, this pollution sensitivity to temperature and water stress might not be true for A. sachalinensis. Moreover, correlating radial growth of this species with precipitation and a climate index, the de Martonne-index revealed a negative effect of increased precipitation on growth. Similarly, higher temperatures were negatively correlated with growth as well. Together with a negative influence of a higher de Martonne-index, it seems that possibly, a high water-availability and warmer temperatures have surprisingly negative effects on growth, albeit correlation coefficients were low (R 2 < 0.4). As this species is adapted to cold boreal climate zones, the negative relationship with temperature seems likely. Moreover, Temperate Plants Database [66] and Chittendon [44] described A. sachalinensis as a species growing well on moist soils, but it is, however, sensitive to water logging. If soils are not well-drained, growth and vitality declines might follow. Since high run-off rates are typical for urban soils, a stronger waterlogging effect during the snow melting season in rural areas could occur.
In summary, neither a prolonged growing season, nor higher temperatures induced the better growth of urban trees. A further correlation with environmental pollution revealed that NO X , SO 2 , CO and NO 2 had strong positive effects on urban tree growth. Only with O X , a negative correlation was found. Related to Figure 7, SO 2 is originated from diesel cars with the use of relatively low-quality gas, and its trend is similar to the emission of NO gas. These (SO X and NO) decreases were accompanied by an increase in O 3 . This O 3 is partly provided from transboundary O 3 from windward regions [47]. Firs growing in Japan, such as Abies sp., were classified as mid-tolerant tree species in terms of elevated O 3 , among 25 tree seedlings, through a big size screening text with Open Top Chambers [46,47]. The fir especially, is sensitive at the needle flushing time of mid-April to early May, when O 3 concentration usually reaches 50-70 ppb. Several studies reported a negative effect of ozone concentrations to tree growth [47,[67][68][69][70]. Due to photosynthesis inhibition found for street trees [71], later bud break times and reduced leafs per bud were found for F. crenata Blume [47,72]. Moreover, these findings overlap with the results of Gregg, et al. [73], who described a strong negative effect of ozone concentration on the growth of clonal poplar. A reduced ozone exposure at urban sites was stated as reason for better growth of urban plants compared to higher ozone concentration and reduced growth of rural plants [73].
However, a decrease in environmental pollution in Sapporo was found during past years, though O X was increasing. Therefore, a monitoring of A. sachalinensis in future will provide further insights. The bias by the age differences of the urban trees in urban and rural surroundings might then be excluded, and the influence of atmospheric concentrations and risks of waterlogging could be analyzed in more detail.

Conclusions
A. sachalinensis trees in Sapporo, northern Japan showed distinct growth patterns related to their growing site. Very old rural trees had lower growth rates than younger urban trees. Moreover, an enhanced growth of trees during past years was found. Growth relationships with environmental pollutions were mainly positive. Only O X seemed to be a major cause of reduced rural tree growth, possibly together with higher tendencies of waterlogged soils, though this was not analyzed in our study. In total, environmental pollution and warmer temperatures, as well as further growing conditions at urban sites (e.g., small planting pits and compacted soils) did not affect tree growth as negatively as could be expected. Funding: This research was funded by the AUDI Environmental Foundation for funding this study (project 5101954: "Reaktionskinetik von Bäumen unter Klimaveränderungen"-"Reaction kinetics of trees under climate change").