Next Article in Journal
Are Plant–Soil Feedbacks Caused by Many Weak Microbial Interactions?
Next Article in Special Issue
Drought Stress Might Induce Sexual Spatial Segregation in Dioecious Populus euphratica—Insights from Long-Term Water Use Efficiency and Growth Rates
Previous Article in Journal
Proton Microbeam Targeted Irradiation of the Gonad Primordium Region Induces Developmental Alterations Associated with Heat Shock Responses and Cuticle Defense in Caenorhabditis elegans
Previous Article in Special Issue
Climatic Warming-Induced Drought Stress Has Resulted in the Transition of Tree Growth Sensitivity from Temperature to Precipitation in the Loess Plateau of China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Impact of Urbanization on Tree Growth and Xylem Anatomical Characteristics

1
Aulin College, Northeast Forestry University, Harbin 150040, China
2
Key Laboratory of Sustainable Forest Ecosystem Management-Ministry of Education, School of Forestry, Northeast Forestry University, Harbin 150040, China
3
State Key Laboratory of Urban and Regional Ecology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China
*
Author to whom correspondence should be addressed.
Biology 2023, 12(11), 1373; https://doi.org/10.3390/biology12111373
Submission received: 21 September 2023 / Revised: 21 October 2023 / Accepted: 23 October 2023 / Published: 27 October 2023
(This article belongs to the Special Issue Dendrochronology in Arid and Semiarid Regions)

Abstract

:

Simple Summary

This study explored how urbanization (urban heat island effect and air pollution) affects urban tree growth. The radial growth and xylem anatomical characteristics of trees in urban and rural areas of Harbin were compared by means of tree-ring anatomy. Results showed that there were significant differences in the growth of both broadleaf trees and conifers between urban and rural areas. Urbanization may have the effect of slowing down tree growth, and conifers may be more sensitive to urbanization. Warming and drying climate in Harbin may be an important factor affecting tree growth. Our findings provide a theoretical basis for the future evaluation of the effect of urbanization on tree growth and ecological benefits, as well as the selection of tree species in urban forests.

Abstract

In the context of the intensification of global urbanization, how urbanization (urban heat island effect and air pollution) affects urban tree growth is not fully understood. In this paper, the radial growth and xylem anatomical characteristics of three different tree species (Quercus mongolica, Fraxinus mandshurica, and Pinus sylvestris var. mongolica) in urban and rural areas of Harbin were compared by means of tree-ring anatomy. The results showed that there were significant differences in the growth of both broadleaf trees and conifers between urban and rural areas. The vessel number, cumulative area of vessels, and theoretical hydraulic conductivity of all tree species in rural areas were higher than those in urban areas, indicating that urbanization may have the effect of slowing down growth. However, broadleaf trees in urban areas had higher vessel density and a greater percentage of a conductive area within xylem and theoretical xylem-specific hydraulic conductivity. The thickness of cell walls and cell wall reinforcement index of P. sylvestris var. mongolica were strongly reduced by air pollution, implying that it may be more sensitive to urbanization. Compared to Q. mongolica, F. mandshurica showed less sensitivity to urbanization. Warming and drying climate in Harbin may be an important factor affecting tree growth.

1. Introduction

With the continuous increase of the world population, the pace of world urbanization will continue to accelerate. It is estimated that by 2050, the world urbanization rate will increase from 56% in 2021 to 68% [1]. Human activities in cities inevitably lead to heat and pollution gas emissions, resulting in two important consequences of urbanization: urban heat island effect and air pollution [2]. Urban heat island effect and air pollution have potentially large impacts on people and vegetation in cities and surrounding areas [3,4].
Urban trees are the most important part of urban vegetation. On the one hand, they reduce air pollution, absorb carbon dioxide, regulate air temperature and humidity, reduce storm-water runoff, etc. [5]. On the other hand, they are heavily affected by the process of urbanization, especially the urban heat island effect and air pollution [6]. Many scholars have conducted a lot of research on the improvement of the urban environment via urban trees [7,8,9,10]. However, there are relatively few studies on the effects of the urban heat island effect and urban pollution on urban trees [3], especially on the anatomy of their trunk xylem. Studying the impact of urbanization on tree anatomy can help us observe changes in the impact of urbanization on tree growth from a more microscopic perspective. Due to the impact of urbanization, there are many differences in the environment between urban areas and surrounding natural areas, such as temperature, precipitation, and air pollution, which lead to different living conditions of trees [6]. High temperature and drought from the urban heat island effect may cause trees to create smaller vessels, while air pollution may prevent trees from forming vessels [11,12]. Therefore, the environmental differences between urban areas and surrounding natural areas and the impact of urbanization on urban vegetation growth have attracted more and more scholars’ attention [4].
Urban heat island (UHI) effect is the result of the physical properties of buildings and other structures, as well as the heat emitted by human activities [13]. The average UHI trend for global cities can reach 2.36 ± 0.69 °C at night and 1.60 ± 1.50 °C during the day [14]. UHI results in temperature differences between urban and rural environments, with consequences for urban vegetation, including altered plant phenology, extended growing seasons, and increased susceptibility to drought stress and insect infestation [2,15,16,17,18,19,20]. At the same time, many studies have proved that there is a correlation between the heat island effect and drought in urban areas [21,22,23,24,25,26]. Trees in cities respond to the urban heat island effect [27]. For example, Moser-Reischl et al. found urban tree growth is usually negatively affected by urban growing conditions [28]. Urban areas can create stressful environments for native trees, leading to increased mortality and ecosystem changes [29]. Rising temperatures and more frequent heat and drought events will exacerbate stress levels in urban trees [30]. The long-term impact of rising temperatures will reduce productivity, as warmer weather and slower tree growth exacerbate air pollution and seasonal climate change [31]. Compared to trees in rural environments, trees from urban areas are more sensitive and vulnerable to the climate [32,33].
Many studies have suggested that urbanization has a negative impact on tree growth. However, some studies also found that urbanization has a positive indirect impact on tree growth, although the land use change brought about by urbanization led to a direct reduction in the number of vegetation [4]. Some urban trees experienced accelerated growth in the past decades [6,34,35,36]. Urban trees take advantage of warm conditions despite the urban heat island effect and moisture limitations [5]. Other studies found no differences in radial growth dynamics between urban and rural trees, except for extreme climatic events in individual years, although the climate conditions in urban areas are warmer [37,38]. However, these studies basically focus on the study of radial growth and isotopes of trees and do not explain and discuss from the microscopic perspective of xylem anatomy.
Air pollutants refer to substances that, at certain concentrations in the atmosphere, may harm the health of organisms, affecting regional air quality and global climate [2,39]. The most frequently monitored air pollutants in urban areas include particulate matter (PM10, PM2.5, UFP), NO2, CO, SO2, and O3. Many studies suggest that air pollution negatively affects tree growth. Thus, tree rings can be used as archives of past air pollution [27,40,41,42,43,44,45]. The mechanism of air pollution affecting tree growth has been extensively studied [3,42,46,47,48]. The emission of SO2, NOx, CO2, dust, and other pollutants seriously interferes with the metabolism of trees [27]. Particulate matter can be adsorbed on leaf surface or cuticle of trees, or it can be metabolized by trees through the stomata or rhizosphere [48,49]. These changes will ultimately affect the photosynthetic rate and transpiration of trees [50,51]. Compared with urban heat island effect, air pollution directly affects tree growth, which is a major growth-limiting factor [5]. Air pollution increases the sensitivity of tree growth to climate factors [44]. Therefore, urban trees exposed to pollution sources are more vulnerable to extreme climate change [52]. Likewise, these studies only quantified tree growth through the radial width but did not explore the impact of pollution effects on tree growth from a more microscopic perspective.
Harbin is the provincial capital of Heilongjiang Province, the largest city in Northern China and an important manufacturing base [53]. In recent years, the urbanization process of Harbin has been very fast, and the environmental problems of the city have also become prominent, e.g., the urban heat island effect and air pollution caused by industrial emissions, road traffic, straw burning, and centralized heating in winter [54,55,56,57]. Quercus mongolica, Fraxinus mandshurica, and Pinus sylvestris var. mongolica are the main afforestation tree species in Harbin, and the amount of P. sylvestris var. mongolica and Q. mongolica tends to increase continuously. At the same time, as a high-latitude city in the northern hemisphere, Harbin is also one of the cities with the most significant temperature increases. However, whether urban heat island effect and air pollution in Harbin will aggravate the impact on tree growth in the context of climate warming is unclear. Therefore, related research is urgently needed.
This study aimed to investigate the impact of urbanization on radial growth and xylem anatomical characteristics of trees so as to provide a theoretical basis for evaluating the impact of urbanization on tree growth and ecological benefits and the selection of urban forest tree species in the future. To this end, we propose the following hypothesis: (1) There are differences in growth between urban and rural trees in Harbin, and the growth of trees in rural areas is better than that in urban areas (Harbin); (2) Air pollution will reduce tree growth, but urban heat island effect may promote tree growth; (3) Urban trees will have more and smaller vessels than rural trees. We temporarily ignore differences in soil, elevation, and natural and plantation forests in this study.

2. Materials and Methods

2.1. Study Areas

In order to study the impact of urban heat island effect and air pollution on tree growth, we selected Harbin as the urban sampling areas (HRB), and Maoershan Forest Ecosystem National Field Scientific Observation and Research Station (LM) 90 km away from Harbin as the control sampling site (rural areas). It was divided into two gradients of air pollution in Harbin. The lightly polluted site is located in the Northeast Forestry University Urban Forestry Demonstration Base, Harbin, China (45°43′ N, 126°37′ E), and another site is on the side of Shuini Road, Daowai District, Harbin, China (45°50′ N, 126°43′ E), which is more polluted. Tree-ring sampling of Q. mongolica, F. mandshurica, and P. sylvestris var. mongolica was conducted at the above sites, respectively (Figure 1, Table A1).
The climate belongs to the mid-temperate continental monsoon climate, with long winters and short summers. The mean annual total precipitation is 569.1 mm. The mean annual temperature is 3.6 °C. The frost-free period is 140 d. The average elevation of Harbin is 151 m. The soil type of Harbin Experimental Forest Farm is black soil [58]. The pollutants discharged are mainly SO2, NOx, and smoke (powder) dust [59].
The main forest type of the Urban Forestry Demonstration Base of Northeast Forestry University (NF) is the artificial pure forest planted from the late 1950s to the early 1960s. There are about 30 pure artificial forests in the sampling areas, and the main tree species include Q. mongolica, P. sylvestris var. mongolica, Betula platyphylla, F. mandshurica, Juglans mandshurica, etc. [60]. Except for Q. mongolica oak forest, a large number of shrubs and herbs have grown under the plantations. The main shrubs include Sambucus williamsii, Ligustrum obtusifolium subsp. suave, Caragana arborescens, Syringa reticulata subsp. Amurensis, and Syringa oblata, and the herbs include Chelidonium majus, Hosta plantaginea, and Trigonotis peduncularis.
Shuini Road (SR), with an altitude of ~132–140 m, is near the Harbin Cement Company of Yatai Group. The sampling site was set in a park about 100 m away from the factory building. The park covered an area of about 30,000 square meters. The soil type is black-brown. The main tree species planted in the park were Salix matsudana and P. sylvestris var. mongolica. Sambucus williamsii, Lespedeza bicolor, Pilea pumila, Amphicarpaea edgeworthii, and Ixeris polycephala distribute under the forest.
The sampling site of the forest areas in contrast to Harbin is located at the National Field Scientific Observation and Research Station of Maoershan Forest Ecosystem in Heilongjiang Province (LM), Shangzhi, China (45°25′ N, 127°38′ E). It is 90 km southeast of Harbin. This area belongs to the temperate continental monsoon climate. The annual mean temperature is 2.8 °C. The frost-free period is between 120–140 d. The mean annual total precipitation is 723 mm, and the mean annual evaporation is 1093 mm. The zonal soil is dark brown forest soil.
The National Field Scientific Observation and Research Station of Maoershan Forest Ecosystem, Maoershan belongs to the branch of Zhangguangcai ridge, with an average altitude of 380 m. The main vegetation is the natural secondary forest formed by the broadleaf Korean pine forest after different degrees of disturbance (harvesting, management, burning, reclamation, etc.). The main tree species are Pinus koraiensis, F. mandshurica, Juglans mandshurica, Q. mongolica, Tilia amurensis, Betula platyphylla, etc.
The climate data were aggregated for both research areas. In general, the average annual temperature of HRB was 1.3 °C higher than that of LM, while the precipitation was 124.2 mm less than that of LM (Figure 2). From a yearly perspective, HRB had a higher temperature, maximum temperature, minimum temperature, and vapor pressure deficit (VPD), as well as lower precipitation (Figure A1). From 2014 to 2020, the Air Quality Index (AQI) of HRB was decreasing, indicating that air pollution is gradually improving. In general, HRB is warmer and drier than LM.

2.2. Sampled Tree Species

The sampled tree species of this research are Q. mongolica, F. mandshurica, and P. sylvestris var. mongolica. Q. mongolica (qm) is a widely distributed tree species in the broadleaf forests of Northeast China and is susceptible to aggravated drought [61]. It has great value as a highly exploited forest resource and is also included in the list of 100 forest plant species vulnerable to climate change [62]. F. mandshurica (fm) is one of the most important temperate hardwood species in the mixed broadleaf—Korean pine forest in Northeast China. It is highly sensitive to climate change and is extensively distributed but discontinuous in Northeast China, Eastern Russia, Northern Japan, and North Korea [63]. P. sylvestris var. mongolica (ps) is the most widely distributed evergreen conifer in the northeastern part of boreal Asian forests [64]. It is also one of the main test objects in dendrochronology research due its role in habitat formation, strong adaptability, distinct annual rings, and high sensitivity [42].

2.3. Sampling and Chronology Development

In 2021, we collected tree ring cores at three sites. Species sampled included two broadleaf trees (Q. mongolica, F. mandshurica) and one conifer (P. sylvestris var. mongolica). Smaller trees, such as those with DBH < 15 cm, were not collected because they might suffer from more competition or other ecological factors. Two or three cores were collected at breast height (~1.3 m) from healthy living individuals with an increment borer (inner diameter: 5.15 mm; Haglöf, Dalarna, Sweden). A total of 321 cores from 140 trees were collected in three sites. Table 1 shows the cores used in developing the chronology; an additional five cores from each tree species at each site were included in the anatomical analysis. The cores were fixed, air-dried, and sanded gradually until each ring was clearly visible. Each ring was then cross-dated to identify possible deletions or false rings, and an accurate calendar year was assigned using the skeleton-plot method under the microscope [65]. Subsequently, the width of the rings was measured with a Velmex measuring system (UniSlide A40 Series, New York, NY, USA) with an accuracy of 0.001 mm. The accuracy of the results was checked using the COFECHA program [66].
The chronologies were developed and the non-climatic effects on tree growth were removed using a modified negative exponential curve (ModNegExp) detrending method by dplR package [67]. The length, standard deviation (SD), mean sensitivity (MS), first-order autocorrelation (AC1), mean inter-series correlation (Rbar), signal-to-noise ratio (SNR), expressed population signal statistic (EPS) [68], and average growth rate (AGR) of the tree-ring chronologies were calculated using dplR package in the R statistical environment [67]. ANOVA was performed to observe the differences in growth between urban and rural trees. Results were tested using the least-significant difference method by agricolae package [69]. Correlation histograms were plotted using the Spearman correlation analysis method by linkET package [70] to analyze the correlation coefficients between monthly values, climatic factors, and radial growth. All analyses and plots were made with R statistical software (R version 4.2.2) [71].

2.4. Anatomical Characteristics Analysis

The tree ring cores were softened in an 89-degree water bath and preprocessed with a GSL1–microtome [72]. The sections were stained with safranin reagent, rested for 5 min, washed with water and 75% ethanol, and finally mounted with 33% glycerin [73]. A semiautomatic digital camera (E3ISPM05000KPA, Zhengzhou, China) connected to a microscope (XSP–GX6A, Shanghai, China) was used to convert the clear slice into a digitized image at 10× magnification with a spatial resolution of 1.52 pixels/µm. Subsequently, professional image analysis software such as Image–Pro Plus (version 6.1) and ROXAS (version 3.0) were used to automatically identify tree-ring image data and measure xylem anatomical characteristics [74].
Tree-ring anatomical parameters, including the mean ring width (MRW, μm), vessel number (VN), mean vessel density (VD, No.∙mm−2), cumulative area of all counted vessels (CTA, mm2), mean percentage of conductive area within xylem (RCTA, %), mean vessel area per rings (MVA, μm2), theoretical hydraulic conductivity (Kh, m4∙MPa−1∙s−1), theoretical xylem-specific hydraulic conductivity (Ks, m2∙MPa−1∙s−1), and overall mean hydraulic diameter (Dh, μm) of broadleaf trees were measured and calculated. For conifers, we calculated the mean ring width (MRW, μm), cell (tracheid) number (CN), mean cell density (CD, No.∙mm−2), cumulative area of all counted cells (CTA, mm2), mean percentage of conductive area within xylem (RCTA, %), mean cell lumen area per rings (MLA, μm2), theoretical hydraulic conductivity (Kh, m4∙MPa−1∙s−1), theoretical xylem-specific hydraulic conductivity (Ks, m2∙MPa−1∙s−1), overall mean hydraulic diameter (Dh, μm), overall mean thickness of all cell walls (CWT, μm), and overall mean cell wall reinforcement index (TB2) [75].
The one-way ANOVA and two-way ANOVA analysis of different tree species and sites were performed for each anatomical parameter to observe the differences in growth between urban and rural trees. Results were tested using the least significant difference method by agricolae package [69]. To evaluate which type of anatomical parameters were most influenced by climate change in urban and rural areas, we used principal component analysis to reduce dimensionality and find the dominant factor characteristics of the first and second principal components and analyze their correlation with climatic factors, respectively. Package factoextra and FactoMineR were used in this step [76,77]. Correlation histograms and heat maps were plotted using the Spearman correlation analysis method by linkET package [70] to analyze the correlation coefficients between monthly values, climatic factors, and anatomical parameters, as well as which factors might mainly affect tree ring parameters. All analyses and plots were made with R statistical software [71].

2.5. Climate Data

The climate data were from the National Meteorological Science Data Center (China Meteorological Data Network, http://data.cma.cn/, access on 31 December 2020). Climate variables include air temperature (TEM), precipitation (Prec), ground temperature (GTEM), maximum air temperature (Tmax), minimum air temperature (Tmin), and relative humidity (RH). Standardized Precipitation Evapotranspiration Index (SPEI) was calculated using the SPEI package [78], and Vapor Pressure Deficit (VPD) was calculated by Equation (1).
V P D = 0.61078 × e 17.27 × T E M T E M + 237.3 × ( 1 R H )
Two sampling sites (NF and SR) shared the same climate data of Harbin Meteorological Station, and the climate data of LM site were from Shangzhi Meteorological Station. Air pollution data were obtained from China Air Quality Online Monitoring and Analysis Platform (https://www.aqistudy.cn/, access on 28 February 2022) and Harbin Ecological Environment Bureau (http://www.harbin.gov.cn/, access on 31 January 2022). The air quality index data of LM was from Shangzhi Monitoring Station. Since Shangzhi is a city, it will be slightly higher than that in the Maoershan forest areas. Variables include Air Quality Index (AQI), PM2.5, PM10, NO2, CO, SO2, and O3–8 h.

3. Results

3.1. Radial Growth of Trees in Urban and Rural Areas

Tree-ring chronology lengths ranged from 39 years to 77 years (Figure 3 and Table 2). F. mandshurica and Q. mongolica from rural areas are natural forests with an average age of more than 70 years, while those from urban areas are plantations with an average age of 50–60 years. The average age of P. sylvestris var. mongolica from heavily polluted areas is even shorter at just under 40 years. The results of the statistical characteristics of the chronologies showed that there was a good mean inter-series correlation (Rbar) between most of the tree samples, indicating a consistency in the variation of ring width between different cores. The first-order autocorrelation coefficient (AC1) of the chronologies fluctuated widely, indicating that the growth status of different tree species in the previous year had a different impact on the current year. The mean sensitivity (MS) of each chronology was 0.28, indicating that the sample had a strong signal of high-frequency fluctuations. Except for P. sylvestris var. mongolica in heavily polluted areas, the signal-to-noise ratio (SNR) of the chronologies was high, and the expressed population signal statistics (EPS) reached above 0.85, indicating that the signal contained in the sample size collected in this study was basically representative of the overall characteristics [79].
The results of the ANOVA showed that the ring width of both F. mandshurica and Q. mongolica in urban areas were significantly lower than those in rural areas, but there was no significant difference in the ring width of P. sylvestris var. mongolica (Figure 4). When the chronology from 1974 to 2020 was segmented, it was found that there was no significant difference between urban and rural F. mandshurica from 1974 to 1989, but there was a significant difference between them from 1990 to 2020 (Figure A2A,D,G). Similarly, the differences between urban and rural Q. mongolica were progressively significant from 1974 to 2020 (Figure A2B,E,H). Both broadleaf trees had significantly higher ring width in rural areas than in urban areas, but the situation was very different for conifers. From 1974 to 1989, there was a significant difference between urban and rural P. sylvestris var. mongolica, but this difference disappeared from 1990 to 2004 (Figure A2C,F). From 2005 to 2020, the ring width of P. sylvestris var. mongolica in rural areas was again significantly higher than in urban areas, but the ring width in heavily polluted areas was higher than that in lightly polluted areas (Figure A2I).
The chronology of F. mandshurica was overall positively correlated with monthly SPEI, with a higher but non-significant coefficient during the growing season (Figure 5A). The chronology of urban Q. mongolica was positively correlated with monthly SPEI and had a significant positive correlation with monthly SPEI during the growing season, while the chronology of rural Q. mongolica had a non-significant positive correlation with monthly SPEI only during the growing season (Figure 5B). The chronology of P. sylvestris var. mongolica in lightly polluted areas was significantly and positively correlated with monthly SPEI, while the chronology of P. sylvestris var. mongolica in rural areas was significantly and negatively correlated with monthly SPEI (Figure 5C). The chronology of P. sylvestris var. mongolica in heavily polluted areas, on the other hand, was positively correlated with monthly SPEI during the year before the end of the growing season, with a higher but non-significant coefficient during the growing season (Figure 5C).

3.2. Xylem Anatomical Characteristics of Trees in Urban and Rural Areas

For F. mandshurica, VD, RCTA, and Ks in urban areas were significantly higher than those in rural areas, while VN, CTA, and Kh in rural areas were significantly higher than those in urban areas (Figure 6). Other parameters had no significant differences between urban and rural sites. For Q. mongolica, VD, RCTA, MVA, and Ks in urban areas were significantly higher than those in rural areas, while VN, CTA, Kh, and Dh in rural areas were significantly higher than those in urban areas. The common points of parameters of F. mandshurica and Q. mongolica in comparison between urban and rural areas were VD, RCTA, Ks, VN, CTA, and Kh. The difference was that Q. mongolica had significant differences in MVA and Dh, while F. mandshurica did not. The results showed that there were significant differences in most xylem anatomical parameters between urban and rural trees (Figure 6).
The two-way ANOVA of the effect of site (urban and rural) and species (Q. mongolica and F. mandshurica) on the anatomical parameters of xylem of trees showed that the effects of sites on most parameters were significant (Table 3). Except Ks and RCTA, species also had significant effects on anatomical parameters. Similarly, the interaction between sites and species had important influence on anatomical parameters, as most of the parameters were significant except Ks (Table 3).
It can be seen from the comparison of areas with different urban pollution levels (heavy pollution—SR, and light pollution—LM) that air pollution has led to a significant decline in CN, CTA, and Kh of P. sylvestris var. mongolica, while Dh in urban areas was significantly higher than that in rural areas (Figure 7). In the two pollution concentration gradients in Harbin, the differences between CWT, TB2, and Kh were more significant than those between urban and rural areas (Figure 7).

3.3. PCA Dimensionality Reduction

The principal component analysis of the main xylem anatomical parameters of the three species in urban and rural areas showed that (Figure 8 and Figure 9) the first component (PC1) of F. mandshurica in urban areas (46.4%) was mainly composed of Kh, MVA, CTA, and Dh (Figure 8A), while the PC1 in rural areas (51.5%) was mainly composed of Ks, Kh, RCTA, and CTA (Figure 8B). The second component (PC2) of F. mandshurica in urban areas (30.4%) was mainly composed of RCTA, VD, and Ks (Figure 8A), while the PC2 in rural areas (26.1%) was mainly composed of VN (Figure 8B). Compared with urban and rural areas, the main anatomical parameters exchanged in PC1 and PC2, and their contributions, were completely opposite (Figure 8A,B).
The PC1 of Q. mongolica in urban areas (52.8%) was mainly composed of VN, CTA, Kh, and VD (Figure 8C), while the PC1 in rural areas (43.7%) was mainly composed of Ks, MVA, RCTA, and Kh (Figure 8D). The PC2 of Q. mongolica in urban areas (27.9%) was mainly composed of Ks and RCTA (Figure 8C), while the PC2 in rural areas (25.3%) was mainly composed of VN (Figure 8D). Like F. mandshurica, compared with the urban and rural areas, the main anatomical parameters exchanged in PC1 and PC2, and their contributions, were completely opposite (Figure 8C,D).
The contribution of each anatomical parameter to PC1 and PC2 of P. sylvestris var. mongolica was similar at both urban pollution levels (Figure 9A,B). The PC1 of P. sylvestris var. mongolica in two urban sites (62.1% and 63.1%) were mainly composed of MLA, Ks, RCTA, Dh, and CWT (Figure 9A,B), while the PC1 in the rural site (56.3%) was mainly composed of MLA, Kh, RCTA, CTA, CD, and TB2 (Figure 9C). The PC2 of P. sylvestris var. mongolica in two sites (30.9% and 20.5%) were mainly composed of CTA and CN (Figure 9A,B), while the PC2 in rural areas (28.9%) was mainly composed of Dh, CWT, and Ks (Figure 9C). Similar to F. mandshurica and Q. mongolica, the contribution of anatomical parameters of the xylem of P. sylvestris var. mongolica to the principal component axis is opposite in urban and rural areas (Figure 9A–C).

3.4. The Relationship between Xylem Anatomical Characteristics of Trees and Climate

The PC1 (Kh, MVA, CTA, Dh) of the main anatomical parameters of F. mandshurica xylem in urban mainly represented hydraulic efficiency and safety (Figure 8A). The PC1 of F. mandshurica xylem in urban areas was not significantly correlated with monthly SPEI (Figure 10A). The PC1 (Ks, Kh, RCTA, CTA) of the main anatomical parameters of F. mandshurica xylem in rural areas mainly represented hydraulic efficiency and safety (Figure 8B). The PC1 of F. mandshurica xylem in rural areas was positively correlated with monthly SPEI, and the non-growing season coefficient was higher but not significant (Figure 10A). The PC2 of the main anatomical parameters of F. mandshurica xylem in urban and rural areas mainly represented vessel growth (Figure 8A,B). The relationship between PC2 of the main anatomical parameters of F. mandshurica and monthly SPEI was similar to that of PC1. That is, the PC2 of F. mandshurica xylem in rural was positively correlated with monthly SPEI (Figure 10B). Similarly, the non-growing season coefficient was higher but not significant. The PC2 of F. mandshurica xylem in urban was significantly positively correlated with monthly SPEI at the end of the growing season (Figure 10B).
The relationship between the principal components of the main xylem parameters of Q. mongolica in rural and monthly SPEI was quite different from that of F. mandshurica (Figure 11A,B). The PC1 (VN, CTA, Kh, VD) of the main anatomical parameters of Q. mongolica xylem in urban areas mainly represented vessel growth (Figure 8C). The PC1 of Q. mongolica xylem in urban areas was positively correlated with monthly SPEI, and the autumn coefficient was higher but not significant (Figure 11A). The PC1 (Ks, MVA, RCTA, Kh) of the main anatomical parameters of Q. mongolica xylem in rural areas mainly represented hydraulic efficiency and safety (Figure 8D). The relationship between PC1 of the main anatomical parameters of Q. mongolica and monthly SPEI in rural areas was not significant (Figure 11A). The PC2 (Ks and RCTA) of the main anatomical parameters of Q. mongolica xylem in urban areas mainly represented hydraulic efficiency and safety (Figure 8C). The PC2 of Q. mongolica xylem in urban areas was not significantly correlated with monthly SPEI (Figure 11B). The PC2 (VN, Dh, CTA, and Kh) of the main anatomical parameters of Q. mongolica xylem in rural areas mainly represented vessel growth (Figure 8D). The PC2 of Q. mongolica xylem in rural areas was positively correlated with monthly SPEI and significant in non-growing season (Figure 11B).
The relationship between the principal components of main xylem parameters of P. sylvestris var. mongolica and monthly SPEI was similar between PC1 and PC2 (Figure 12A,B). The PC1 (MLA, Ks, RCTA, Dh, and CWT) of the main anatomical parameters of P. sylvestris var. mongolica xylem in urban areas mainly represented hydraulic efficiency and safety (Figure 9A,B). The PC1 of P. sylvestris var. mongolica xylem in urban was positively correlated with monthly SPEI, and the autumn coefficient was higher but not significant (Figure 12A). The PC1 (MLA, Kh, RCTA, CTA, and TB2) of the main anatomical parameters of Q. mongolica xylem in rural areas mainly represented hydraulic safety (Figure 9C). The PC1 of P. sylvestris var. mongolica xylem in rural areas was negatively correlated with monthly SPEI, and the spring coefficient was significant (Figure 12A). The PC2 (CTA and CN) of the main anatomical parameters of P. sylvestris var. mongolica xylem in urban areas mainly represented hydraulic safety (Figure 9A,B). The PC2 of P. sylvestris var. mongolica xylem in urban areas was not significantly correlated with the SPEI (Figure 12B). The PC2 (Dh, CWT, and Ks) of the main anatomical parameters of Q. mongolica xylem in rural areas mainly represented hydraulic efficiency (Figure 9C). The PC2 of P. sylvestris var. mongolica xylem in rural areas was negatively correlated with monthly SPEI but was not significant (Figure 12B).
In general, the PC2 of the main anatomical parameters of Q. mongolica xyelm and the main anatomical parameters of P. sylvestris var. mongolica xylem in rural areas were more sensitive to SPEI, but Q. mongolica was mostly positively correlated, while P. sylvestris var. mongolica was mostly negatively correlated (Figure 11B and Figure 12A).
Except MVA, the anatomical parameters of Q. mongolica in rural areas were mainly negatively correlated with monthly SPEI, while the anatomical parameters of Q. mongolica in urban areas were mainly positively correlated with monthly SPEI (Figure A3). However, the correlation between F. mandshurica and monthly SPEI was much weaker; only Dh in rural areas and VD in urban areas were significantly positively correlated in the last months of the year (Figure A3). The responses of P. sylvestris var. mongolica to SPEI can be divided into two categories: growth anatomical parameters (CN, CD, CWT, and TB2) were generally positively correlated, while hydraulic anatomical parameters (CTA, RCTA, MLA, Kh, Ks, and Dh) were generally negatively correlated (Figure A3). Their response months were also different. The responses of P. sylvestris var. mongolica in urban areas to SPEI were concentrated in the second half of the year, while the responses of P. sylvestris var. mongolica in rural areas were distributed from the previous growing season to the first half of the year (Figure 12 and Figure A3).

4. Discussion

The study found significant differences in radial growth and xylem anatomical characteristics between urban and rural populations in both broadleaf trees and conifers. The changes of VN, CTA, and Kh were basically stable in all tree species, i.e., higher in rural areas than in urban areas (Figure 6). The difference of Q. mongolica between urban and rural areas was quite large, while F. mandshurica did not respond strongly to climatic factors. P. sylvestris var. mongolica may be more strongly affected by urbanization than the two broadleaf tree species.

4.1. Effect of Urban Heat Island on Anatomical Characteristics of Trunk Xylem

The urban heat island effect seems to have a negative impact on the growth of trees in Harbin. However, it has been debated whether the heat island effect promotes or hinders tree growth. Some studies have shown that long-term temperature increases in urban forests ultimately reduce productivity and slow tree growth [31]. Urban trees are sensitive to high temperatures and drought, which is thought to be one of the main factors limiting radial tree growth [33]. However, quite a few researchers have found that the urban heat island effect promotes tree growth. It has been noted that urban trees have experienced accelerated growth since the 1960s and tend to grow faster than trees in rural areas because urban trees are able to take advantage of warm conditions without being limited by water availability [5,35,36]. Street trees benefit from high temperatures and low precipitation in built-up urban areas [6]. It was also found that despite warmer climatic conditions in urban areas, there was no difference in the growth of urban and rural trees [37,38]. Affected by the heat island effect, the urban environment of Harbin is warmer and drier than the surrounding rural environment (Figure A1). For urban trees, warming and drying in Harbin decreased xylem anatomical parameters such as VN, CTA, and Kh. However, VD, RCTA, and Ks of broadleaf trees and Dh of conifers seem to be promoted by the urban heat island effect, which may be able to compensate for the decrease in hydraulic conductivity caused by the reduction in the number and total areas of vessels. We found that the growth decline of P. sylvestris var. mongolica in Harbin appeared to be more serious at the end of the 20th century and the beginning of the 21st century (Figure 3), which was considered to be the period when the urban heat island effect was strongest [54]. However, broadleaf trees did not show a similar trend, which may be due to the compensating mechanism mentioned above. In response to the seasons, Gillner et al. (2014) found that temperature and moisture availability from April to July of the current year and the summer and autumn of the previous year were the main determinants of radial growth. Broadleaf tree species were more sensitive to summer heat and drought in urban areas (summer signal), while conifers were suitable for analyzing urban heat islands in late winter and early spring (winter signal) [33]. For radial growth, our study was broadly consistent with this conclusion but not for anatomical parameters, implying that radial growth and anatomical characteristics did not correspond exactly.
Our results showed that the urban heat island effect in Harbin led to a decrease in the vessel area of all tree species (Figure 6C and Figure 7E). Islam et al. found that diffuse porous Chukrasia tabularis decreased vessel area under the stress environment, particularly in hot and dry conditions [80]. E. grandis seedlings grown at higher temperatures had a 40% reduction in the lumen area of xylem vessels [81]. At the same time, the urban heat island effect in Harbin resulted in increased vessel density of broadleaf tree species (Figure 6B). Islam et al. found that, under hot and dry conditions, not only the vessel area (hydraulic efficiency) decreased, but also the vessel density (hydraulic safety) increased [80]. Gea-Izquierdo et al. showed that trees adapted to more xeric conditions by reducing radial growth and hydraulic diameters, as well as increasing vessel density, but xylem was not adapted to respond to long-term temperature increases [82]. Under heat treatment, Salix pulchra Cham. had higher vessel density than ambient-growing plants [83]. However, some researchers have published contrary conclusions. Subfossil oaks from the Younger Dryas cold period, up to 3.5 °C cooler than present temperatures, had more vessel frequencies than living oaks [84]. In Central Europe, the vessel density of pedunculate oak (Quercus robur L.) was negatively correlated with the temperature of the current year in floodplain sites [85]. In Harbin, in order to adapt to a drier, warmer urban environment, trees sought a trade-off between maximizing conductivity and creating safety margins against cavitation and embolism, resulting in smaller vessel area and greater vessel density [80].
The urban heat island effect in Harbin also significantly increased the theoretical xylem-specific hydraulic conductivity (Ks) of both broadleaf tree species, while that of conifer species was not significantly increased (Figure 6G and Figure 7I). Similarly, the specific xylem conductivity of Salix pulchra Cham. after heat treatment was 2.5 times higher than that of the ambient-growing plants [83]. For ponderosa pine, larger tracheid diameters at elevated temperatures caused an increase in xylem-specific hydraulic conductivity [86]. However, Voelker et al. (2012) found that Younger Dryas cold period subfossil oaks had a higher theoretical xylem conductivity due to their higher vessel density compared to the current oaks growing under much warmer conditions. In Harbin, the larger vessel frequencies of urban trees compensated for the decrease in vessel dimensions, so the theoretical xylem conductivity was always higher than that of rural trees [84].

4.2. Effects of Air Pollution on the Anatomical Characteristics of Trunk Xylem

Many studies have shown that environmental pollution will have a negative impact on tree growth [5,45]. PM and Al, Ba, and Zn in the air have significant effects on the interannual change of urban tree growth [5]. Moser-Reischl et al. also found a decline in Khaya senegalensis growth in downtown Hanoi compared to rural trees [28]. This may be due to low groundwater levels and high pollution. Among the two pollution concentration gradients in Harbin, CN, CTA, CWT, Kh, and TB2, all decreased significantly in the heavily polluted area (Figure 7). Among these anatomical parameters, the differences in CWT and TB2 were greater than those of urban and rural controls (Figure 7C,D). These differences were mainly caused by differences in pollution because the climate of the two sampling sites is relatively close. For P. sylvestris var. mongolica, pollution caused a strong decrease in CWT and TB2 parameters related to cell wall formation, but urban heat islands did not induce such changes. Pollution may affect the physiological activities of urban trees, leading to a decline in tree growth. Sensuła et al. (2015) found that the pollution emissions would cause serious disturbances to tree metabolism, resulting in a sharp decline in tree growth. Despite the apparent long-term decline in tree growth caused by pollution, tree growth will resume after the pollution slows [44]. Air pollution and drought in the 1950s destabilized pine populations and caused more severe damage to the trees [42]. Conversely, environmental pollution also appears to have a positive effect on tree growth. [34]. In Harbin, air pollution reduced the vessel number of trees (Figure 6A and Figure 7A). However, distinctly different results were published. Rao et al. found that the characteristics of xylem of trees growing under pollution stress had a relatively large number of vessels [87]. Coal-smoke pollutants led to an increase in the number of vessels with reduced sizes of Euphorbia hirta L. [88]. This may be because the large amount of dust pollution in Harbin hinders the photosynthesis and transpiration of trees, thus slowing down the formation of vessels and causing their numbers to decline.
The vessel density results in Harbin were more consistent with the conclusions of most studies. That is, air pollution increased the vessel density of trees in heavily polluted areas (factory areas) rather than lightly polluted areas (urban forests) (Figure 6B and Figure 7B). Sukumaran found that the vessel density of Cassia occidentalis increased to a very significant level under the influence of air pollutants [89]. Prosopis spicigera Linn., growing under the influence of combined air pollutants, had a significantly higher vessel density [90]. In Brazil, the vessel density of Cecropia glazioui Sneth. was significantly higher in polluted sites than in unpolluted and intermediate sites [91]. However, Mahmooduzzafar et al. showed that air pollution from coal burning at the Badarpur Thermal Power Plants, New Delhi, reduced the vessel density of Indian blackberry trunks (Syzygium cumini Skeel.) [92]. In Harbin, the increasing vessel density of wood under the stress of polluted environments resulted in a lower vulnerability of polluted trees [93].
Air pollution in Harbin also reduced the vessel area (Figure 6C and Figure 7E). Although the cumulative area of all vessels (CTA) decreased, the mean percentage of vessels (RCTA) in the xylem of broadleaf trees increased (Figure 6D), which corresponds to a reduction in the increment of tree growth (Figure 4). Mahmooduzzafar et al. found that the vessel area on transverse surfaces of Indian blackberries was lower than that in unpolluted areas [92]. However, Sukumaran showed that the vessel area of Cassia occidentalis significantly increased under the influence of air pollutants [89]. In addition, coal smoke caused an increase in the percentage of wood vessels of Prosopis cineraria in Delhi, India [12]. In Harbin, air pollution inhibited pigment concentration, NR activity, and sugar content, and it promoted stomatal index and nitrate and sulfur contents and resulted in low stomatal conductance, leading to a decrease in net photosynthetic rate and, thus, a reduction in the vessel area [12].

4.3. Uncertainty in Interpretation of Results of Urban and Rural Growth and Xylem Anatomical Differences

Our research did not take into account the differences between Harbin and its surrounding rural in soil, elevation, natural forest, and plantation forest, which would lead to differences in the growth and anatomical characteristics of trees between urban and rural areas. F. mandshurica and Q. mongolica in the rural areas (Maoershan) are natural forests with an average age of more than 60 years, while those in the urban area (Harbin) are planted forests with an average age of 50 to 60 years. Natural forests may grow slower than plantations due to competition. However, trees in natural forests are more resistant to pollution. In addition, there is a difference in altitude between Maoershan and Harbin. Trees grow higher in the rural area (Maoershan) than in the urban area (Harbin). As the altitude rises, the temperature will drop, which will reduce the growth of trees. However, the lower altitude of Harbin will cause the temperature to rise. At the same time, lower altitudes can worsen droughts and may also affect growth. Furthermore, the rural area (Maoershan) is a natural coniferous forest soil with rich nutrients, while the urban area (Harbin) is an unnatural soil after human disturbance, and the soil fertility and texture are worse than those in rural. Therefore, from the perspective of soil condition, the quality of urban soil will be much lower than that of rural soil, which may have negative effects. We ignored these differences for now, but the impact of these differences is still significant.
In addition, there are still many deficiencies in this study. The Harbin research areas and rural areas are 87 km apart, but the two sites in Harbin are only 14 km apart, and the detailed climate data of these two sites have not been obtained yet. Currently, the pollution situation of the two sampling sites can only be inferred from the surrounding environment. Pollution data for rural areas (Shangzhi City) are also too scant. The research method in this study is not deep enough, and many contents remain unexplored. Future research in this direction can be started from the differences of physiological processes and chemical composition of trees caused by urbanization.

5. Conclusions

There were significant differences in radial growth and xylem anatomical characteristics between urban and rural populations in both broadleaf trees and conifers. On parameters such as VN, CTA, and Kh, trees in rural areas were higher than those in urban areas. However, broadleaf trees in urban areas generally had higher VD, RCTA, and Ks, which was not observed in conifers. Meanwhile, the overall mean thickness of all cell walls (CWT) and the overall mean cell-wall reinforcement index (TB2) of P. sylvestris var. mongolica were strongly reduced in heavily polluted areas. This may mean that conifers are more sensitive to urbanization because they do not have the mechanisms to counteract the decrease in the vessel number, area, and hydraulic conductivity caused by the urban heat island effect and air pollution, and they are strongly affected in cell-wall formation. Among the two different broadleaf tree species, F. mandshurica had almost no significant correlation with climatic factors, and it may be more adaptable to the urban environment than Q. mongolica. Warming and drying in Harbin may be an important factor affecting tree growth.

Author Contributions

X.W. and X.G. conceived the study. X.G., B.Z., Z.C., W.S. and Z.L. conducted field samples and collected other data. X.G. and X.W. undertook all lab work and data analyses and wrote the first draft of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (42177421 and 41877426) and the Fund of Eco-meteorological Innovation Open Laboratory in Northeast China, China Meteorological Bureau (stqx2018zd02).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Datasets are available on request. The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Acknowledgments

We are very grateful to the staff of the forestry bureau workers, who helped us work in the field.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A

Table A1. Sampling site information.
Table A1. Sampling site information.
SiteSpeciesForestSlopeOrientationAltitude
UrbanNFF. mandshuricaArtificial
pure
forest
Level ground141 m
Q. mongolica141 m
P. sylvestris var. mongolica141 m
SRP. sylvestris var. mongolica132–140 m
RuralLMF. mandshuricaNatural
secondary
forest
GullyDark slope383 m
Q. mongolicaUphillSunny slope395–420 m
P. sylvestris var. mongolicaMiddle and
downhill
Sunny slope361 m
Figure A1. The change of annual mean (A), mean ground surface (B), maximum (C) and minimum (D) temperatures, total precipitation (E), mean SPEI (F), mean vapor pressure deficit (VPD) (G), and mean air quality index (AQI) (H) in urban (HRB) and rural (LM) areas from 2002 to 2020. The fitting equation and p-value were shown.
Figure A1. The change of annual mean (A), mean ground surface (B), maximum (C) and minimum (D) temperatures, total precipitation (E), mean SPEI (F), mean vapor pressure deficit (VPD) (G), and mean air quality index (AQI) (H) in urban (HRB) and rural (LM) areas from 2002 to 2020. The fitting equation and p-value were shown.
Biology 12 01373 g0a1
Figure A2. Comparison of the ring width of F. mandshurica (A,D,G), Q. mongolica (B,E,H), and P. sylvestris var. mongolica (C,F,I) among urban (SR—heavy pollution, NF—light pollution) and rural (LM) sites from 1974–1989 (AC), 1990–2004 (DF), and 2005–2020 (GI). The significance values of the ring width are indicated. NS. no significant difference; * p < 0.05; ** p < 0.01; *** p < 0.001.
Figure A2. Comparison of the ring width of F. mandshurica (A,D,G), Q. mongolica (B,E,H), and P. sylvestris var. mongolica (C,F,I) among urban (SR—heavy pollution, NF—light pollution) and rural (LM) sites from 1974–1989 (AC), 1990–2004 (DF), and 2005–2020 (GI). The significance values of the ring width are indicated. NS. no significant difference; * p < 0.05; ** p < 0.01; *** p < 0.001.
Biology 12 01373 g0a2
Figure A3. Correlation analysis between the xylem anatomical parameters and monthly SPEI of F. mandshurica (A,B), Q. mongolica (C,D), and P. sylvestris var. mongolica (E,F) in rural (A,C,E) and urban (B,D,F) areas. The significance values of the xylem anatomical characteristics are indicated. * p < 0.05; ** p < 0.01; *** p < 0.001.
Figure A3. Correlation analysis between the xylem anatomical parameters and monthly SPEI of F. mandshurica (A,B), Q. mongolica (C,D), and P. sylvestris var. mongolica (E,F) in rural (A,C,E) and urban (B,D,F) areas. The significance values of the xylem anatomical characteristics are indicated. * p < 0.05; ** p < 0.01; *** p < 0.001.
Biology 12 01373 g0a3

References

  1. Martin, W. United Nations Human Settlements Programme (UN-HABITAT). In Max Planck Encyclopedia of Public International Law; Oxford University Press: Oxford, UK, 2022; pp. 1–13. [Google Scholar]
  2. Oke, T.R.; Mills, G.; Christen, A.; Voogt, J.A. Urban Climates; Cambridge University Press: Cambridge, UK, 2017; ISBN 9781139016476. [Google Scholar]
  3. Gregg, J.W.; Jones, C.G.; Dawson, T.E. Urban Ozone Depletion: Why a Tree Grows Better in New York City. Nature 2003, 424, 183–187. [Google Scholar] [CrossRef] [PubMed]
  4. Zhang, L.; Yang, L.; Zohner, C.M.; Crowther, T.W.; Li, M.; Shen, F.; Guo, M.; Qin, J.; Yao, L.; Zhou, C. Direct and Indirect Impacts of Urbanization on Vegetation Growth across the World’s Cities. Sci. Adv. 2022, 8, eabo0095. [Google Scholar] [CrossRef] [PubMed]
  5. Locosselli, G.M.; de Camargo, E.P.; Moreira, T.C.L.; Todesco, E.; Andrade, M.D.F.; de André, C.D.S.; de André, P.A.; Singer, J.M.; Ferreira, L.S.; Saldiva, P.H.N.; et al. The Role of Air Pollution and Climate on the Growth of Urban Trees. Sci. Total Environ. 2019, 666, 652–661. [Google Scholar] [CrossRef]
  6. Dahlhausen, J.; Rötzer, T.; Biber, P.; Uhl, E.; Pretzsch, H. Urban Climate Modifies Tree Growth in Berlin. Int. J. Biometeorol. 2018, 62, 795–808. [Google Scholar] [CrossRef]
  7. Locosselli, G.M.; Chacón-Madrid, K.; Zezzi Arruda, M.A.; Pereira de Camargo, E.; Lopes Moreira, T.C.; Saldiva de André, C.D.; Afonso de André, P.; Singer, J.M.; Nascimento Saldiva, P.H.; Buckeridge, M.S. Tree Rings Reveal the Reduction of Cd, Cu, Ni and Pb Pollution in the Central Region of São Paulo, Brazil. Environ. Pollut. 2018, 242, 320–328. [Google Scholar] [CrossRef] [PubMed]
  8. Rahman, M.A.; Stratopoulos, L.M.F.; Moser-Reischl, A.; Zölch, T.; Häberle, K.H.; Rötzer, T.; Pretzsch, H.; Pauleit, S. Traits of Trees for Cooling Urban Heat Islands: A Meta-Analysis. Build. Environ. 2020, 170, 106606. [Google Scholar] [CrossRef]
  9. Rötzer, T.; Rahman, M.A.; Moser-Reischl, A.; Pauleit, S.; Pretzsch, H. Process Based Simulation of Tree Growth and Ecosystem Services of Urban Trees under Present and Future Climate Conditions. Sci. Total Environ. 2019, 676, 651–664. [Google Scholar] [CrossRef]
  10. Salmond, J.A.; Tadaki, M.; Vardoulakis, S.; Arbuthnott, K.; Coutts, A.; Demuzere, M.; Dirks, K.N.; Heaviside, C.; Lim, S.; MacIntyre, H.; et al. Health and Climate Related Ecosystem Services Provided by Street Trees in the Urban Environment. Environ. Health Glob. Access Sci. Source 2016, 15, 95–111. [Google Scholar] [CrossRef]
  11. Islam, M.; Rahman, M.; Bräuning, A. Xylem Anatomical Responses of Diffuse Porous Chukrasia Tabularis to Climate in a South Asian Moist Tropical Forest. For. Ecol. Manag. 2018, 412, 9–20. [Google Scholar] [CrossRef]
  12. Iqbal, M.; Mahmooduzzafar; Aref, I.M.; Khan, P.R. Behavioral Responses of Leaves and Vascular Cambium of Prosopis cineraria (L.) Druce to Different Regimes of Coal-Smoke Pollution. J. Plant Interact. 2010, 5, 117–133. [Google Scholar] [CrossRef]
  13. Parker, D.E. Urban Heat Island Effects on Estimates of Observed Climate Change. WIREs Clim. Chang. 2010, 1, 123–133. [Google Scholar] [CrossRef]
  14. Varquez, A.C.G.; Kanda, M. Global Urban Climatology: A Meta-Analysis of Air Temperature Trends (1960–2009). NPJ Clim. Atmos. Sci. 2018, 1, 32. [Google Scholar] [CrossRef]
  15. Mimet, A.; Pellissier, V.; Quénol, H.; Aguejdad, R.; Dubreuil, V.; Rozé, F. Urbanisation Induces Early Flowering: Evidence from Platanus acerifolia and Prunus cerasus. Int. J. Biometeorol. 2009, 53, 287–298. [Google Scholar] [CrossRef] [PubMed]
  16. Xing, X.; Zhang, M.; Li, K.; Hao, P.; Dong, L. Spatial Heterogeneity of First Flowering Date in Beijing’s Main Urban Area and Its Response to Urban Thermal Environment. Int. J. Biometeorol. 2022, 66, 1929–1954. [Google Scholar] [CrossRef]
  17. Kabano, P.; Lindley, S.; Harris, A. Evidence of Urban Heat Island Impacts on the Vegetation Growing Season Length in a Tropical City. Landsc. Urban Plan. 2021, 206, 103989. [Google Scholar] [CrossRef]
  18. Melaas, E.K.; Wang, J.A.; Miller, D.L.; Friedl, M.A. Interactions between Urban Vegetation and Surface Urban Heat Islands: A Case Study in the Boston Metropolitan Region. Environ. Res. Lett. 2016, 11, 054020. [Google Scholar] [CrossRef]
  19. Cregg, B.M.; Dix, M.E. Tree Moisture Stress and Insect Damage in Urban Areas in Relation to Heat Island Effects. J. Arboric. 2001, 27, 8–17. [Google Scholar] [CrossRef]
  20. Zipper, S.C.; Schatz, J.; Singh, A.; Kucharik, C.J.; Townsend, P.A.; Loheide, S.P. Urban Heat Island Impacts on Plant Phenology: Intra-Urban Variability and Response to Land Cover. Environ. Res. Lett. 2016, 11, 054023. [Google Scholar] [CrossRef]
  21. Kumar, R.; Mishra, V.; Kumar, R.; Mishra, V. Sustainability of Smart Cities under Climate Variability and Climate Change in India. AGUFM 2015, 2015, B33E-0783. [Google Scholar]
  22. Huang, S.; Zhang, X.; Yang, L.; Chen, N.; Nam, W.H.; Niyogi, D. Urbanization-Induced Drought Modification: Example over the Yangtze River Basin, China. Urban Clim. 2022, 44, 101231. [Google Scholar] [CrossRef]
  23. Dong, L.; Mitra, C.; Greer, S.; Burt, E. The Dynamical Linkage of Atmospheric Blocking to Drought, Heatwave and Urban Heat Island in Southeastern US: A Multi-Scale Case Study. Atmosphere 2018, 9, 33. [Google Scholar] [CrossRef]
  24. Winguth, A.M.E.; Kelp, B. The Urban Heat Island of the North-Central Texas Region and Its Relation to the 2011 Severe Texas Drought. J. Appl. Meteorol. Climatol. 2013, 52, 2418–2433. [Google Scholar] [CrossRef]
  25. Standard, L.C. The North Central Texas Urban Heat Island Magnitude and Its Relation to Severe and Extreme Drought. Master’s Thesis, The University of Texas at Arlington, Arlington, TX, USA, 2021. [Google Scholar]
  26. Nimac, I.; Herceg-Bulić, I.; Žuvela-Aloise, M.; Žgela, M. Impact of North Atlantic Oscillation and Drought Conditions on Summer Urban Heat Load—A Case Study for Zagreb. Int. J. Climatol. 2022, 42, 4850–4867. [Google Scholar] [CrossRef]
  27. Sensuła, B.; Wilczyński, S.; Opała, M. Tree Growth and Climate Relationship: Dynamics of Scots Pine (Pinus sylvestris L.) Growing in the Near-Source Region of the Combined Heat and Power Plant During the Development of the Pro-Ecological Strategy in Poland. Water Air Soil Pollut. 2015, 226, 220. [Google Scholar] [CrossRef]
  28. Moser-Reischl, A.; Uhl, E.; Rötzer, T.; Biber, P.; van Con, T.; Tan, N.T.; Pretzsch, H. Effects of the Urban Heat Island and Climate Change on the Growth of Khaya senegalensis in Hanoi, Vietnam. For. Ecosyst. 2018, 5, 37. [Google Scholar] [CrossRef]
  29. Wilde, E.M.; Maxwell, J.T. Comparing Climate-Growth Responses of Urban and Non-Urban Forests Using L. tulipifera Tree-Rings in Southern Indiana, USA. Urban For. Urban Green. 2018, 31, 103–108. [Google Scholar] [CrossRef]
  30. Gillner, S.; Bräuning, A.; Roloff, A. Dendrochronological Analysis of Urban Trees: Climatic Response and Impact of Drought on Frequently Used Tree Species. Trees-Struct. Funct. 2014, 28, 1079–1093. [Google Scholar] [CrossRef]
  31. Ordóñez, C.; Duinker, P.N. Assessing the Vulnerability of Urban Forests to Climate Change. Environ. Rev. 2014, 22, 311–321. [Google Scholar] [CrossRef]
  32. Hunt, A.; Watkiss, P. Climate Change Impacts and Adaptation in Cities: A Review of the Literature. Clim. Chang. 2011, 104, 13–49. [Google Scholar] [CrossRef]
  33. Schneider, C.; Neuwirth, B.; Schneider, S.; Balanzategui, D.; Elsholz, S.; Fenner, D.; Meier, F.; Heinrich, I. Using the Dendro-Climatological Signal of Urban Trees as a Measure of Urbanization and Urban Heat Island. Urban Ecosyst. 2022, 25, 849–865. [Google Scholar] [CrossRef]
  34. Moser-Reischl, A.; Rötzer, T.; Biber, P.; Ulbricht, M.; Uhl, E.; Qu, L.; Koike, T.; Pretzsch, H. Growth of Abies sachalinensis Along an Urban Gradient Affected by Environmental Pollution in Sapporo, Japan. Forests 2019, 10, 707. [Google Scholar] [CrossRef]
  35. Moser, A.; Uhl, E.; Rötzer, T.; Biber, P.; Dahlhausen, J.; Lefer, B.; Pretzsch, H. Effects of Climate and the Urban Heat Island Effect on Urban Tree Growth in Houston. Open J. For. 2017, 07, 428–445. [Google Scholar] [CrossRef]
  36. Pretzsch, H.; Biber, P.; Uhl, E.; Dahlhausen, J.; Schütze, G.; Perkins, D.; Rötzer, T.; Caldentey, J.; Koike, T.; van Con, T.; et al. Climate Change Accelerates Growth of Urban Trees in Metropolises Worldwide. Sci. Rep. 2017, 7, 15403. [Google Scholar] [CrossRef]
  37. Kukarskih, V.V.; Devi, N.M.; Bubnov, M.O.; Komarova, A.V.; Agafonov, L.I. Urban Heat Island of Ekaterinburg: Does It Affect Radial Growth of the Scots Pine? J. Sib. Fed. Univ.-Biol. 2022, 15, 264–278. [Google Scholar] [CrossRef]
  38. Shumilov, O.I.; Kasatkina, E.A.; Timonen, M.; Herva, H.; Kirtsideli, I.; Kanatjev, A.G.; Shumilov, O.I.; Kasatkina, E.A.; Timonen, M.; Herva, H.; et al. “Urban Heat Island” Effect on Tree Growth at Several Cities of Northern Europe. EGUGA 2010, 12, 1422. [Google Scholar]
  39. Yadav, I.C.; Devi, N.L. Biomass Burning, Regional Air Quality, and Climate Change. In Encyclopedia of Environmental Health; Elsevier: Amsterdam, The Netherlands, 2019; pp. 386–391. ISBN 9780444639523. [Google Scholar]
  40. Ballikaya, P.; Marshall, J.; Cherubini, P. Can Tree-Ring Chemistry Be Used to Monitor Atmospheric Nanoparticle Contamination over Time? Atmos. Environ. 2022, 268, 118781. [Google Scholar] [CrossRef]
  41. Binda, G.; Di Iorio, A.; Monticelli, D. The What, How, Why, and When of Dendrochemistry: (Paleo) Environmental Information from the Chemical Analysis of Tree Rings. Sci. Total Environ. 2021, 758, 143672. [Google Scholar] [CrossRef]
  42. Kukarskih, V.V.; Devi, N.M.; Bubnov, M.O.; Komarova, A.V.; Agafonov, L.I. Radial Growth of Scots Pine in Urban and Rural Populations of Ekaterinburg Megalopolis. Dendrochronologia 2022, 74, 125974. [Google Scholar] [CrossRef]
  43. Sensuła, B.; Wilczyński, S. Climatic Signals in Tree-Ring Width and Stable Isotopes Composition of Pinus sylvestris L. Growing in the Industrialized Area Nearby Kędzierzyn-Koźle. Geochronometria 2017, 44, 240–255. [Google Scholar] [CrossRef]
  44. Sensuła, B.; Wilczyński, S. Records of Anthropogenic Pollution in Silesia Captured in Scots Pine Tree Rings: Analysis by Radiocarbon, Stable Isotopes, and Basal Area Increment Analysis. Water Air Soil Pollut. 2022, 233, 143. [Google Scholar] [CrossRef]
  45. Sensuła, B.; Wilczyński, S.; Monin, L.; Allan, M.; Pazdur, A.; Fagel, N. Variations of Tree Ring Width and Chemical Composition of Wood of Pine Growing in the Area Nearby Chemical Factories. Geochronometria 2017, 44, 226–239. [Google Scholar] [CrossRef]
  46. Fowler, D.; Cape, J.N.; Unsworth, M.H.; Mayer, H.; Crowther, J.M.; Jarvis, P.G.; Gardiner, B.; Shuttleworth, W.J. Deposition of Atmospheric Pollutants on Forests. Philos. Trans. R. Soc. Lond. B Biol. Sci. 1989, 324, 247–265. [Google Scholar] [CrossRef]
  47. McDonald, A.G.; Bealey, W.J.; Fowler, D.; Dragosits, U.; Skiba, U.; Smith, R.I.; Donovan, R.G.; Brett, H.E.; Hewitt, C.N.; Nemitz, E. Quantifying the Effect of Urban Tree Planting on Concentrations and Depositions of PM10 in Two UK Conurbations. Atmos. Environ. 2007, 41, 8455–8467. [Google Scholar] [CrossRef]
  48. Prajapati, S.K. Ecological Effect of Airborne Particulate Matter on Plants. Environ. Skept. Crit. 2012, 1, 12–22. [Google Scholar]
  49. Grantz, D.A.; Garner, J.H.B.; Johnson, D.W. Ecological Effects of Particulate Matter. Environ. Int. 2003, 29, 213–239. [Google Scholar] [CrossRef]
  50. Wang, X.; Wang, C.; Wu, J.; Miao, G.; Chen, M.; Chen, S.; Wang, S.; Guo, Z.; Wang, Z.; Wang, B.; et al. Intermediate Aerosol Loading Enhances Photosynthetic Activity of Croplands. Geophys. Res. Lett. 2021, 48, e91893. [Google Scholar] [CrossRef]
  51. Wang, B.; Wang, Z.; Wang, C.; Wang, X.; Li, J.; Jia, Z.; Li, P.; Wu, J.; Chen, M.; Liu, L. Field Evidence Reveals Conservative Water Use of Poplar Saplings under High Aerosol Conditions. J. Ecol. 2021, 109, 2190–2202. [Google Scholar] [CrossRef]
  52. Sæbø, A.; Benedikz, T.; Randrup, T.B. Selection of Trees for Urban Forestry in the Nordic Countries. Urban For. Urban Green. 2003, 2, 101–114. [Google Scholar] [CrossRef]
  53. Song, S.; Xu, D.; Hu, S.; Shi, M. Ecological Network Optimization in Urban Central District Based on Complex Network Theory: A Case Study with the Urban Central District of Harbin. Int. J. Environ. Res. Public Health 2021, 18, 1427. [Google Scholar] [CrossRef]
  54. Li, W. Effect of Urbanization on Urban Heat Island: A Case Study in Harbin. Environ. Sci. Manag. 2020, 45, 58–62. [Google Scholar]
  55. Wu, H.; Liu, Y.; Liuli, L.; Liu, Y.; Liu, B.; Zheng, T.; Wang, P. Estimation method of air pollution load of straw burning in Harbin. Huanjing Kexue Xuebao/Acta Sci. Circumstantiae 2020, 40, 3803–3810. [Google Scholar] [CrossRef]
  56. Yu, M.; Cao, L.; Han, Z.; Han, F.; Liu, Y.; Wang, X. Study on the Influence of Urbanization Process on Air Temperature in Harbin. HEILONGJIANG Meteorol. 2021, 38, 27–29. [Google Scholar] [CrossRef]
  57. Luo, Y.; Liu, S.; Gao, Y.; Zhang, M.; Xuan, L. Analysis of PM2.5 Pollution Characteristics during Heating Period in Harbin City. J. GreenSci. Technol. 2021, 23, 42–44. [Google Scholar] [CrossRef]
  58. Chang, J. The Study on Investigation of Soil in Urban Green Space and Effects for the Improving in Harbin. Master’s Thesis, Northeast Forestry University, Harbin, China, 2015. [Google Scholar]
  59. Qiu, Y.; Li, L.; Lu, H.; Li, J.; Wang, X.; Jiang, L.; Su, J. Estimation of Air Pollution Emission from Cement Plant in Beijing Based on Thermal Anomaly Detection. Res. Environ. Sci. 2020, 33, 2265–2271. [Google Scholar] [CrossRef]
  60. Zhang, J.; Liu, B. Research on Population Dispersal of Quercus mongolica in the City Forestry Demonstration Base of Northeast Forestry University. Mod. Agric. Technol. 2013, 135–137. [Google Scholar] [CrossRef]
  61. Zhang, X.; Lv, P.; Xu, C.; Huang, X.; Rademacher, T. Dryness Decreases Average Growth Rate and Increases Drought Sensitivity of Mongolia Oak Trees in North China. Agric. For. Meteorol. 2021, 308–309, 108611. [Google Scholar] [CrossRef]
  62. Moon, N.H.; Moon, G.H.; Chun, J.H.; Shin, M.Y. Dendroclimatological Analysis and Tree-Ring Growth Prediction of Quercus mongolica. For. Sci. Technol. 2020, 16, 32–40. [Google Scholar] [CrossRef]
  63. Chen, L.; Zhu, L.; Liu, S.; Lei, P.; Yuan, D.; Li, Z.; Wang, X. Xylem Features Detrending Methods Matter: A Case Study on Earlywood Vessels of Fraxinus mandshurica. Ecol. Indic. 2021, 130, 108041. [Google Scholar] [CrossRef]
  64. Li, W.; Jiang, Y.; Dong, M.; Du, E.; Wu, F.; Zhao, S.; Xu, H. Species-Specific Growth-Climate Responses of Dahurian Larch (Larix gmelinii) and Mongolian pine (Pinus sylvestris var. mongolica) in the Greater Khingan Range, Northeast China. Dendrochronologia 2021, 65, 125803. [Google Scholar] [CrossRef]
  65. Stokes, M.A.; Smiley, T.L. An Introduction to Tree-Ring Dating; UAPress: Tucson, AZ, USA, 1996. [Google Scholar]
  66. Holmes, R.L. Computer-Assisted Quality Control in Tree-Ring Dating and Measurement. Tree-Ring Bull. 1983, 43, 69–78. [Google Scholar]
  67. Bunn, A.G. A Dendrochronology Program Library in R (DplR). Dendrochronologia 2008, 26, 115–124. [Google Scholar] [CrossRef]
  68. Wigley, T.M.L.; Briffa, K.R.; Jones, P.D. On the Average Value of Correlated Time Series, with Applications in Dendroclimatology and Hydrometeorology. J. Appl. Meteorol. Climatol. 1984, 23, 201–213. [Google Scholar] [CrossRef]
  69. de Mendiburu, F. Agricolae: Statistical Procedures for Agricultural Research. 2021. Available online: https://cran.r-project.org/web/packages/agricolae/index.html (accessed on 27 April 2023).
  70. Huang, H. LinkET: Everything Is Linkable. 2021. Available online: https://github.com/Hy4m/linkET (accessed on 27 April 2023).
  71. R Core Team. R A Lang. Environ. Stat. Comput. R Found. Stat. Comput. Vienna, Austria. 2021. Available online: http://www.R-project.org/ (accessed on 27 April 2023).
  72. Gärtner, H.; Lucchinetti, S.; Schweingruber, F.H. New Perspectives for Wood Anatomical Analysis in Dendrosciences: The GSL1-Microtome. Dendrochronologia 2014, 32, 47–51. [Google Scholar] [CrossRef]
  73. Gärtner, H.; Cherubini, P.; Fonti, P.; von Arx, G.; Schneider, L.; Nievergelt, D.; Verstege, A.; Bast, A.; Schweingruber, F.H.; Büntgen, U. A Technical Perspective in Modern Tree-Ring Research—How to Overcome Dendroecological and Wood Anatomical Challenges. JoVE (J. Vis. Exp.) 2015, 2015, e52337. [Google Scholar] [CrossRef]
  74. von Arx, G.; Carrer, M. ROXAS—A New Tool to Build Centuries-Long Tracheid-Lumen Chronologies in Conifers. Dendrochronologia 2014, 32, 290–293. [Google Scholar] [CrossRef]
  75. Peters, R.L.; Balanzategui, D.; Hurley, A.G.; von Arx, G.; Prendin, A.L.; Cuny, H.E.; Björklund, J.; Frank, D.C.; Fonti, P. RAPTOR: Row and Position Tracheid Organizer in R. Dendrochronologia 2018, 47, 10–16. [Google Scholar] [CrossRef]
  76. Lê, S.; Josse, J.; Husson, F. FactoMineR: An R Package for Multivariate Analysis. J. Stat. Softw. 2008, 25, 1–18. [Google Scholar] [CrossRef]
  77. Kassambara, A.; Mundt, F. Factoextra: Extract and Visualize the Results of Multivariate Data Analyses. 2020. Available online: https://cran.r-project.org/web/packages/factoextra/index.html (accessed on 27 April 2023).
  78. Beguería, S.; Vicente-Serrano, S.M. SPEI: Calculation of the Standardised Precipitation-Evapotranspiration Index. 2017. Available online: https://cran.r-project.org/web/packages/SPEI/index.html (accessed on 27 April 2023).
  79. Li, Z.; Chen, W.; Wei, J.; Maierdang, K.; Zhang, Y.; Zhang, S.; Wang, X. Tree-Ring Growth Responses of Liaodong Oak (Quercus wutaishanica) to Climate in the Beijing Dongling Mountain of China. Shengtai Xuebao 2021, 41, 27–37. [Google Scholar] [CrossRef]
  80. Islam, M.; Rahman, M.; Bräuning, A. Long-Term Hydraulic Adjustment of Three Tropical Moist Forest Tree Species to Changing Climate. Front. Plant Sci. 2018, 9, 1761. [Google Scholar] [CrossRef]
  81. Thomas, D.S.; Montagu, K.D.; Conroy, J.P. Temperature Effects on Wood Anatomy, Wood Density, Photosynthesis and Biomass Partitioning of Eucalyptus grandis Seedlings. Tree Physiol. 2007, 27, 251–260. [Google Scholar] [CrossRef]
  82. Gea-Izquierdo, G.; Fonti, P.; Cherubini, P.; Martín-Benito, D.; Chaar, H.; Cañellas, I. Xylem Hydraulic Adjustment and Growth Response of Quercus canariensis Willd. to Climatic Variability. Tree Physiol. 2012, 32, 401–413. [Google Scholar] [CrossRef] [PubMed]
  83. Gorsuch, D.M.; Oberbauer, S.F. Effects of Mid-Season Frost and Elevated Growing Season Temperature on Stomatal Conductance and Specific Xylem Conductivity of the Arctic Shrub, Salix Pulchra. Tree Physiol. 2002, 22, 1027–1034. [Google Scholar] [CrossRef]
  84. Voelker, S.L.; Noirot-Cosson, P.-E.; Stambaugh, M.C.; McMurry, E.R.; Meinzer, F.C.; Lachenbruch, B.; Guyette, R.P. Spring Temperature Responses of Oaks Are Synchronous with North Atlantic Conditions during the Last Deglaciation. Ecol. Monogr. 2012, 82, 169–187. [Google Scholar] [CrossRef]
  85. Tumajer, J.; Treml, V. Response of Floodplain Pedunculate Oak (Quercus robur L.) Tree-Ring Width and Vessel Anatomy to Climatic Trends and Extreme Hydroclimatic Events. For. Ecol. Manag. 2016, 379, 185–194. [Google Scholar] [CrossRef]
  86. Maherali, H.; DeLucia, E.H. Interactive Effects of Elevated CO2 and Temperature on Water Transport Inponderosa Pine. Am. J. Bot. 2000, 87, 243–249. [Google Scholar] [CrossRef] [PubMed]
  87. Rao, K.S.; Rajput, K.S.; Srinivas, T. Seasonal Cambial Anatomy and Development of Xylem in Dalbergia Sissoo Growing under the Influence of Combined Air Pollutants. J. Sustain. For. 2004, 18, 73–88. [Google Scholar] [CrossRef]
  88. Gupta, M.C.; Ghouse, A.K.M. Effects of Coal-Smoke Pollutants from Different Sources on the Growth, Chlorophyll Content, Stem Anatomy and Cuticular Traits of Euphorbia hirta L. Environ. Pollut. 1987, 47, 221–229. [Google Scholar] [CrossRef]
  89. Sukumaran, D. Effect of Air Pollution on the Anatomy Some Tropical Plants. Appl. Ecol. Environ. Sci. 2014, 2, 32–36. [Google Scholar] [CrossRef]
  90. Rajput, K.S.; Rao, K.S.; Kim, Y.S. Cambial Activity and Wood Anatomy in Prosopis spicigera (Mimosaceae) Affected by Combined Air Pollutants. IAWA J. 2008, 29, 209–219. [Google Scholar] [CrossRef]
  91. Alves, E.S. The Effects of the Pollution on Wood of Cecropia glazioui (Cecropiaceae). IAWA J. 1995, 16, 69–80. [Google Scholar] [CrossRef]
  92. Mahmooduzzafar; Hegazy, S.S.; Aref, I.M.; Iqbal, M. Anatomical Changes in the Wood of Syzygium Cumini Exposed to Coal-Smoke Pollution. J. Food Agric. Environ. 2010, 8, 959–964. [Google Scholar]
  93. Gupta, M.C.; Iqbal, M. Ontogenetic Histological Changes in the Wood of Mango (Mangifera indica L. cv Deshi) Exposed to Coal-Smoke Pollution. Environ. Exp. Bot. 2005, 54, 248–255. [Google Scholar] [CrossRef]
Figure 1. Three sample sites distributed in Heilongjiang province, Northeast China (A). NF (slightly polluted urban site)–Northeast Forestry University Urban Forestry Demonstration Base, Harbin, northeast China; SR (heavily polluted urban site)–Shuini Road, Harbin, Northeast China; LM (rural site)–Maoershan Forest Ecosystem National Field Scientific Observation and Research Station, Shangzhi, Northeast China. (B) Temperature differences among the three sampling sites (2015); (C) PM2.5 concentration differences among the three sampling sites (2001).
Figure 1. Three sample sites distributed in Heilongjiang province, Northeast China (A). NF (slightly polluted urban site)–Northeast Forestry University Urban Forestry Demonstration Base, Harbin, northeast China; SR (heavily polluted urban site)–Shuini Road, Harbin, Northeast China; LM (rural site)–Maoershan Forest Ecosystem National Field Scientific Observation and Research Station, Shangzhi, Northeast China. (B) Temperature differences among the three sampling sites (2015); (C) PM2.5 concentration differences among the three sampling sites (2001).
Biology 12 01373 g001
Figure 2. Comparisons of monthly mean temperature (A), annual mean temperature (B), monthly mean precipitation (C), and annual total precipitation (D) between urban (HRB) and rural (LM) (1951–2020). Standard deviation was shown with error bars. The significance values of the annual climate data are indicated. *** p < 0.001.
Figure 2. Comparisons of monthly mean temperature (A), annual mean temperature (B), monthly mean precipitation (C), and annual total precipitation (D) between urban (HRB) and rural (LM) (1951–2020). Standard deviation was shown with error bars. The significance values of the annual climate data are indicated. *** p < 0.001.
Biology 12 01373 g002
Figure 3. Tree-ring chronologies of F. mandshurica (fm), Q. mongolica (qm) and P. sylvestris var. mongolica (ps) xylem among urban (SR—heavy pollution, NF—light pollution) and rural (LM) sites.
Figure 3. Tree-ring chronologies of F. mandshurica (fm), Q. mongolica (qm) and P. sylvestris var. mongolica (ps) xylem among urban (SR—heavy pollution, NF—light pollution) and rural (LM) sites.
Biology 12 01373 g003
Figure 4. Comparison of the ring width of F. mandshurica (A), Q. mongolica (B), and P. sylvestris var. mongolica (C) among urban (SR—heavy pollution, NF—light pollution) and rural (LM) sites. The significance values of the ring width are indicated. NS. no significant difference; ** p < 0.01; *** p < 0.001.
Figure 4. Comparison of the ring width of F. mandshurica (A), Q. mongolica (B), and P. sylvestris var. mongolica (C) among urban (SR—heavy pollution, NF—light pollution) and rural (LM) sites. The significance values of the ring width are indicated. NS. no significant difference; ** p < 0.01; *** p < 0.001.
Biology 12 01373 g004
Figure 5. Correlation analysis between the tree-ring chronologies and monthly SPEI of F. mandshurica (A), Q. mongolica (B), and P. sylvestris var. mongolica (C) among urban (SR—heavy pollution, NF—light pollution) and rural (LM) sites. The dash line indicates the 0.05 significant level.
Figure 5. Correlation analysis between the tree-ring chronologies and monthly SPEI of F. mandshurica (A), Q. mongolica (B), and P. sylvestris var. mongolica (C) among urban (SR—heavy pollution, NF—light pollution) and rural (LM) sites. The dash line indicates the 0.05 significant level.
Biology 12 01373 g005
Figure 6. Comparison of mean vessel number (VN, A), vessel density (VD, B), cumulative area of all counted vessels (CTA, C), mean percentage of conductive area within xylem (RCTA, D), mean vessel area per rings (MVA, E), theoretical hydraulic conductivity (Kh, F), theoretical xylem-specific hydraulic conductivity (Ks, G), and overall mean hydraulic diameter (Dh, H) of the xylem of F. mandshurica (fm) and Q. mongolica (qm) between urban (NF) and rural (LM) areas. The significance values of the xylem anatomical characteristics are indicated. NS. no significant difference; * p < 0.05; ** p < 0.01; *** p < 0.001.
Figure 6. Comparison of mean vessel number (VN, A), vessel density (VD, B), cumulative area of all counted vessels (CTA, C), mean percentage of conductive area within xylem (RCTA, D), mean vessel area per rings (MVA, E), theoretical hydraulic conductivity (Kh, F), theoretical xylem-specific hydraulic conductivity (Ks, G), and overall mean hydraulic diameter (Dh, H) of the xylem of F. mandshurica (fm) and Q. mongolica (qm) between urban (NF) and rural (LM) areas. The significance values of the xylem anatomical characteristics are indicated. NS. no significant difference; * p < 0.05; ** p < 0.01; *** p < 0.001.
Biology 12 01373 g006
Figure 7. Comparison of cell number (CN, A), cell density (CD, B), cell-wall thickness (CWT, C), cell-wall reinforcement index (TB2, D), cumulative total cell area (CTA, E), mean percentage of conductive area within xylem (RCTA, F), mean lumen area (MLA, G), theoretical hydraulic conductivity (Kh, H), theoretical xylem-specific hydraulic conductivity (Ks, I), and mean hydraulic diameter (Dh, J) of P. sylvestris var. mongolica xylem among urban (SR—heavy pollution, NF—light pollution) and rural (LM) sites. The significance values of the xylem anatomical characteristics are indicated. NS. no significant difference; * p < 0.05; ** p < 0.01; *** p < 0.001.
Figure 7. Comparison of cell number (CN, A), cell density (CD, B), cell-wall thickness (CWT, C), cell-wall reinforcement index (TB2, D), cumulative total cell area (CTA, E), mean percentage of conductive area within xylem (RCTA, F), mean lumen area (MLA, G), theoretical hydraulic conductivity (Kh, H), theoretical xylem-specific hydraulic conductivity (Ks, I), and mean hydraulic diameter (Dh, J) of P. sylvestris var. mongolica xylem among urban (SR—heavy pollution, NF—light pollution) and rural (LM) sites. The significance values of the xylem anatomical characteristics are indicated. NS. no significant difference; * p < 0.05; ** p < 0.01; *** p < 0.001.
Biology 12 01373 g007
Figure 8. PCA analysis of xylem anatomical characteristics of F. mandshurica (A,B) and Q. mongolica (C,D) in urban (A,C) and rural (B,D) areas. Different colors represent the contribution of different factors.
Figure 8. PCA analysis of xylem anatomical characteristics of F. mandshurica (A,B) and Q. mongolica (C,D) in urban (A,C) and rural (B,D) areas. Different colors represent the contribution of different factors.
Biology 12 01373 g008
Figure 9. PCA analysis of xylem anatomical characteristics of P. sylvestris var. mongolica in urban (A—heavy pollution, B—light pollution) and rural (C) sites. Different colors represent the contribution of different factors.
Figure 9. PCA analysis of xylem anatomical characteristics of P. sylvestris var. mongolica in urban (A—heavy pollution, B—light pollution) and rural (C) sites. Different colors represent the contribution of different factors.
Biology 12 01373 g009
Figure 10. Correlation analysis between the first (A) and second (B) principal components of xylem anatomical parameters and monthly SPEI of F. mandshurica in urban (NF) and rural (LM) areas. The dash line indicates the 0.05 significant level.
Figure 10. Correlation analysis between the first (A) and second (B) principal components of xylem anatomical parameters and monthly SPEI of F. mandshurica in urban (NF) and rural (LM) areas. The dash line indicates the 0.05 significant level.
Biology 12 01373 g010
Figure 11. Correlation analysis between the first (A) and second (B) principal components of xylem anatomical parameters and monthly SPEI of Q. mongolica in urban (NF) and rural (LM) areas. The dash line indicates the 0.05 significant level.
Figure 11. Correlation analysis between the first (A) and second (B) principal components of xylem anatomical parameters and monthly SPEI of Q. mongolica in urban (NF) and rural (LM) areas. The dash line indicates the 0.05 significant level.
Biology 12 01373 g011
Figure 12. Correlation analysis between the first (A) and second (B) principal components of xylem anatomical parameters and monthly SPEI of P. sylvestris var. mongolica among urban (SR—heavy pollution, NF—light pollution) and rural (LM) sites. The dash line indicates the 0.05 significant level.
Figure 12. Correlation analysis between the first (A) and second (B) principal components of xylem anatomical parameters and monthly SPEI of P. sylvestris var. mongolica among urban (SR—heavy pollution, NF—light pollution) and rural (LM) sites. The dash line indicates the 0.05 significant level.
Biology 12 01373 g012
Table 1. Sampling sites, tree species, and sample information.
Table 1. Sampling sites, tree species, and sample information.
SiteSpeciesNo. of TreesNo. of CoresMean DBHMean Age
UrbanNFF. mandshurica153025.44059
Q. mongolica142923.52951
P. sylvestris var. mongolica153031.85351
SRP. sylvestris var. mongolica153022.59335
RuralLMF. mandshurica153035.64762
Q. mongolica132634.08769
P. sylvestris var. mongolica153035.04043
Table 2. The length, standard deviation (SD), mean sensitivity (MS), first-order autocorrelation (AC1), mean inter-series correlation (Rbar), signal-to-noise ratio (SNR), expressed population signal statistic (EPS), and average growth rate (AGR) of the tree-ring chronologies.
Table 2. The length, standard deviation (SD), mean sensitivity (MS), first-order autocorrelation (AC1), mean inter-series correlation (Rbar), signal-to-noise ratio (SNR), expressed population signal statistic (EPS), and average growth rate (AGR) of the tree-ring chronologies.
Length (Year)SDMSAC1RbarSNREPSAGR
NF_fm65 (1956–2020)0.2760.2970.5800.46713.1210.9291.690
LM_fm77 (1944–2020)0.1950.2050.5770.42411.0310.9172.252
NF_qm57 (1964–2020)0.2350.3570.2330.50514.2950.9351.443
LM_qm76 (1945–2020)0.1840.2620.1360.59719.2470.9511.869
SR_ps39 (1982–2020)0.2410.3520.3280.1893.2640.7652.579
NF_ps58 (1963–2020)0.2610.2260.7860.52916.8310.9442.542
LM_ps47 (1974–2020)0.1860.2610.2690.57420.210.9533.022
Table 3. Two-way ANOVA analysis of the effects of site (urban and rural) and species (F. mandshurica and Q. mongolica) on xylem anatomical parameters.
Table 3. Two-way ANOVA analysis of the effects of site (urban and rural) and species (F. mandshurica and Q. mongolica) on xylem anatomical parameters.
DfSum SqMean SqF-Valuep-Value
VDSite159259210.090.002
Species120,65220,652351.62<0.001
Site × species12912914.960.029
VNSite18044804444.61<0.001
Species139,13839,138217.06<0.001
Site × species14256425623.6<0.001
CTASite10.1970.197163.12<0.001
Species10.6690.669554.68<0.001
Site × species10.0050.0054.090.047
DhSite170057005100.14<0.001
Species16794679497.13<0.001
Site × species16819681997.48<0.001
KhSite16.88 × 10−136.88 × 10−13110.09<0.001
Species11.46 × 10−121.46 × 10−12233.83<0.001
Site × species18.41 × 10−148.41× 10−1413.46<0.001
KsSite10.0980.09838.01<0.001
Species10.0030.0031.050.309
Site × species10.0020.0020.890.349
MVASite11.03 × 1081.03 × 10826.34<0.001
Species12.30 × 1092.30 × 109588.80<0.001
Site × species13.77 × 1073.77 × 1079.650.003
RCTASite1610.7610.7103.96<0.001
Species11.51.50.260.614
Site × species1111.0111.018.89<0.001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gao, X.; Zhao, B.; Chen, Z.; Song, W.; Li, Z.; Wang, X. The Impact of Urbanization on Tree Growth and Xylem Anatomical Characteristics. Biology 2023, 12, 1373. https://doi.org/10.3390/biology12111373

AMA Style

Gao X, Zhao B, Chen Z, Song W, Li Z, Wang X. The Impact of Urbanization on Tree Growth and Xylem Anatomical Characteristics. Biology. 2023; 12(11):1373. https://doi.org/10.3390/biology12111373

Chicago/Turabian Style

Gao, Xiaohui, Binqing Zhao, Zecheng Chen, Wenqi Song, Zongshan Li, and Xiaochun Wang. 2023. "The Impact of Urbanization on Tree Growth and Xylem Anatomical Characteristics" Biology 12, no. 11: 1373. https://doi.org/10.3390/biology12111373

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop