*Article* **Before Becoming a World Heritage: Spatiotemporal Dynamics and Spatial Dependency of the Soundscapes in Kulangsu Scenic Area, China**

**Zhu Chen 1,2 , Tian-Yuan Zhu 1,3, Jiang Liu 1,\* and Xin-Chen Hong <sup>1</sup>**


**Abstract:** Kulangsu is a famous scenic area in China and a World Heritage Site. It is important to obtain knowledge with regard to the status of soundscape and landscape resources and their interrelationships in Kulangsu before it became a World Heritage. The objective of this study was to explore the spatial dependency of the soundscapes in Kulangsu, based on the spatiotemporal dynamics of soundscape and landscape perceptions, including perceived sound sources, soundscape quality, and landscape satisfaction degree, and the spatial landscape characteristics, including the distance to green spaces, normalized difference vegetation index, and landscape spatial patterns. The results showed that perception of soundscape and landscape were observed in significant spatiotemporal dynamics, and the dominance of biological sounds in all sampling periods and human sounds in the evening indicated that Kulangsu scenic area had a good natural environment and a developed night-time economy, respectively. The green spaces and commercial lands may contribute to both the soundscape pleasantness and eventfulness. Moreover, the soundscape quality was dependent on the sound dominant degree and landscape satisfaction degree but not on the landscape characteristics. The GWR model had better goodness of fit than the OLS model, and possible non-linear relationships were found between the soundscape pleasantness and the variables of perceived sound sources and landscape satisfaction degree. The GWR models with spatial stationarity were found to be more effective in understanding the spatial dependence of soundscapes. In particular, the data applied should ideally include a complete temporal dimension to obtain a relatively high fitting accuracy of the model. These findings can provide useful data support and references for future planning and design practices, and management strategies for the soundscape resources in scenic areas and World Heritage Sites.

**Keywords:** soundscape mapping; soundscape quality; spatiotemporal dynamics; landscape satisfaction; landscape pattern; scenic area

## **1. Introduction**

Rapid urbanization has led to an increasing proportion of the population settling in cities worldwide [1]. Living in high-density built-up areas limits urban residents' access to nature and may expose them to certain environmental hazards, such as noise pollution [2]. Scenic areas possess abundant natural and cultural landscapes as well as infrastructure and facilities for human activities [3], which is a critical vehicle for people to interact with natural environment. The scenic areas located in urban areas are an important and special urban green space (UGS), providing a variety of important ecosystem services to urban residents, such as air purification [4], biodiversity conservation [5], nature-based recreation [6], and noise and quiet reduction [7]. Furthermore, they can improve the quality of urban environments, promote sustainable lifestyles, contribute to human health and

**Citation:** Chen, Z.; Zhu, T.-Y.; Liu, J.; Hong, X.-C. Before Becoming a World Heritage: Spatiotemporal Dynamics and Spatial Dependency of the Soundscapes in Kulangsu Scenic Area, China. *Forests* **2022**, *13*, 1526. https://doi.org/10.3390/f13091526

Academic Editor: Isabella De Meo

Received: 24 August 2022 Accepted: 11 September 2022 Published: 19 September 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

well-being [8], and even reduce mortality [9]. However, with the development of tourist industry, the conflicting interests e.g., between noise of tourist activities and experiencing the sounds of nature have become one of the focal issues in the management of scenic areas [10,11], because auditory perception is as important as visual perception in people's visiting experience [12,13].

Currently, scholars no longer equate the good acoustic environment with simply reducing noise levels [14,15], but rather emphasize the importance of people's subjective perception of the acoustic environment following the soundscape concept of *ISO 12913*-*1* [16]. Against this backdrop, more and more studies have explored and focused on the human perceived acoustic environment [17,18]. However, exploring soundscapes needs a multifaceted perspective because soundscape perception is associated with not only acoustic components but also non-acoustic factors [19]. Sound sources are the most basic and important acoustic components in creating soundscapes in urban areas [20], because the sound dominance in a soundscape is able to affect the spatiotemporal dynamics and the outcome of people's perception of the soundscape [21,22]. Regarding non-acoustic factors, landscape compositions account for a high proportion of the influence on soundscape perception, with can be summarized as subjective and objective aspects [18]. The former is mainly due to people's perception varied from landscape elements, such as naturalness [23], visual quality [24], urban contexts [25], features of landscape and architecture [17], audiovisual coherence [26], and infrastructure services [27]. The latter is related to the landscape characteristics, for example, accessibility [28], vegetation coverage [29], landscape spatial pattern [7], and biodiversity [30].

These relationships indicate that a soundscape may possess spatial dependencies on such factors. The spatial dependence is seen as a normal extension of the first law in geography [31]—"everything is related to everything else, but near things are more related than distant things". It may occur due to the spatial dimensions with regard to social-cultural contexts and economic factors [32]. Previous studies found that the temporal dimension is also a critical aspect for exploring the spatial dependence, and neglecting the temporal characteristics could lead to a misunderstanding of the "real" measure of spatial dependence over time [33,34]. However, to date only a few studies have investigated spatial dependence of soundscape quality but with some deficiencies. For instance, Hong and Jeon [35] explored the spatial dependence of urban soundscapes, nonetheless, solely on the perceived sound sources. Rice et al. [36] explored the spatial dependence of noise abatement on the features of protected areas, but they did not consider people's perception of the acoustic environment. In general, current studies have neither explored the spatial dependency of soundscape quality with the principal components, i.e., pleasantness and eventfulness [37], nor included variables in terms of landscape perception and green space features. Besides, the temporal characteristics of the spatial dependence of soundscape have not been effectively explored either. Exploring the spatial dependence of soundscape can clarify the interrelationship between soundscape perception and impact factors, which may help planners and managers identify the main disturbances to the acoustic environment in scenic areas and therefore find solutions and protection measures [35].

In 2016, Kulangsu was suggested as a cultural heritage by United Nations Educational, Scientific, and Cultural Organization (UNESCO), and subsequently was successfully listed as a World Heritage at the 41st World Heritage Congress in Krakow, Poland, in July 2017 [38]. Therefore, 2016 was a "landmark year" that represents a turning point for Kulangsu from China to the world. Kulangsu is a unique and historic international settlement, an important UGS, and a famous tourist attraction [39], the natural and cultural resources of which constitutes the unique soundscapes of outstanding universal value [40]. Unfortunately, such soundscapes are suffering from excessive disturbance and destruction by human activities, whether present, past, or even future. Accordingly, exploring the spatiotemporal dynamics of soundscapes and landscapes as well as the spatial dependence of soundscape quality in Kulangsu Island has significant implications for soundscape planning and management, and soundscape resources conservation, not only in UGS but also in the scenic area and

even other similar World Heritage Sites. Given this importance, the objective of this study is (1) to visualize and analyze the spatiotemporal dynamics with regard to soundscape and landscape perception, as well as the spatial landscape features in Kulangsu Island; and (2) to examine the spatial dependence of soundscape quality on these compositions. To this end, both global and local spatial regression methods were employed. importance, the objective of this study is (1) to visualize and analyze the spatiotemporal dynamics with regard to soundscape and landscape perception, as well as the spatial landscape features in Kulangsu Island; and (2) to examine the spatial dependence of soundscape quality on these compositions. To this end, both global and local spatial regression methods were employed.

destruction by human activities, whether present, past, or even future. Accordingly, exploring the spatiotemporal dynamics of soundscapes and landscapes as well as the spatial dependence of soundscape quality in Kulangsu Island has significant implications for soundscape planning and management, and soundscape resources conservation, not only in UGS but also in the scenic area and even other similar World Heritage Sites. Given this

#### **2. Materials and Methods 2. Materials and Methods**

#### *2.1. Study Area 2.1. Study Area*

The study area is located on Kulangsu Island (Figure 1a), a World Heritage and one of the National 5A level tourist attractions in China. The area of it is about 1.92 km<sup>2</sup> , with a length of 2.3 km from north to south and 1.6 km from east to west. The green spaces occupied more than one-third of the island (Figure 1b), which was the highest percentage (32.36%) of land use type (according to 2017 land use vector data of Kulangsu Island obtained from Xiamen Municipal Natural Resources and Planning Bureau). The types of green space include parks, squares, and woodlands [41]. The study area is located on Kulangsu Island (Figure 1a), a World Heritage and one of the National 5A level tourist attractions in China. The area of it is about 1.92 km², with a length of 2.3 km from north to south and 1.6 km from east to west. The green spaces occupied more than one-third of the island (Figure 1b), which was the highest percentage (32.36%) of land use type (according to 2017 land use vector data of Kulangsu Island obtained from Xiamen Municipal Natural Resources and Planning Bureau). The types of green space include parks, squares, and woodlands [41].

*Forests* **2022**, *13*, x FOR PEER REVIEW 3 of 21

**Figure 1.** Case study area: (**a**) location of Kulangsu Island in Fujian Province, China; (**b**) land use type with sampling points. **Figure 1.** Case study area: (**a**) location of Kulangsu Island in Fujian Province, China; (**b**) land use type with sampling points.

The present study was conducted from July 17, 2016, to July 21, 2016, in the application process of World Heritage. In July 2017, Kulangsu Island was officially inscribed on the World Heritage List. Based on the pilot study and relevant literature [7,12,42], 4 maincategories and 19 sub-categories sound sources were identified in the area (Table 1). The present study was conducted from 17 July 2016, to 21 July 2016, in the application process of World Heritage. In July 2017, Kulangsu Island was officially inscribed on the World Heritage List. Based on the pilot study and relevant literature [7,12,42], 4 maincategories and 19 sub-categories sound sources were identified in the area (Table 1).

**Table 1.** Classification of sound sources in the study area. **Table 1.** Classification of sound sources in the study area.


cessibility, and interviewed random passersby by questionnaire surveys on site in the morning (8:00–11:00), afternoon (13:00–16:00), and evening (17:00–20:00) of the days. A We selected 52 observation points according to the land use type, function, and accessibility, and interviewed random passersby by questionnaire surveys on site in the morning (8:00–11:00), afternoon (13:00–16:00), and evening (17:00–20:00) of the days. A team of 12 college students from landscape architecture faculty of Fuzhou University got involved to help conduct the field survey and distribute the questionnaire to random people in the study area. They were professionally trained before executing the field survey in order to ensure the quality of study. The data includes personal information and

We selected 52 observation points according to the land use type, function, and ac-

subjective perception ratings for the soundscape and landscape of the interviewees. A total of 703 valid questionnaires were returned, with 10 to 15 questionnaires on each sample site, indicating that enough and accurate results can be achieved [43,44]. The interviewees' information is shown in Figure S1 (Supplementary Material). Analysis of the questionnaire data in SPSS 25.0 showed that the alpha coefficient was 0.86, indicating that the reliability of the questionnaire data was good and suitable for further analysis.

#### *2.2. Data Collection*

## 2.2.1. Soundscape and Landscape Perception

Participants were asked to evaluate the sound source, soundscape quality, and landscape satisfaction degree of the environment, based on perceived indicators using a Likert 5-point scale (Table 2).


**Table 2.** Detailed information for each perceived indicator.

Notes: j is the jth sample, i is the ith source, and n is the sample size.

#### 2.2.2. Analysis of Landscape Characteristics

This study objectively quantifies the landscape features within Kulangsu in two aspects: (1) features of green spaces, including Distance to green spaces (DtGS) and Normalized difference vegetation index (NDVI), and (2) landscape spatial patterns, including patterns of the green space class and overall landscape. DtGS was calculated based on the "Euclidean distance" spatial analysis tool in GIS, which is measured in the projection unit of the raster and calculated from one cell center to the other cell center, and can objectively measure the spatial distance of the global scale [45], with the Equation (1) shows:

$$D = \sqrt{\left(\mathbf{x}\_2 - \mathbf{x}\_1\right)^2 + \left(y\_2 - y\_1\right)^2} \tag{1}$$

where (*x*1, *y*1) are the coordinates of one point, (*x*2, *y*2) are the coordinates of the other point; *D* is the Euclidean distance between (*x*1, *y*1) and (*x*2, *y*2).

We imported the collected Sentinel-2 level 2A images as an ensemble into the Google Earth Engine (GEE) JavaScript-based code editor environment [46], to calculate the value of NDVI for each time series image. Only images with less than 10% cloudiness in the study area were extracted, to ensure data integrity, and the calculated images were between March 2017 (the first Sentinel-2 available image in GEE) and January 2019. All data were top-of-atmosphere (TOA) images and had been atmospherically corrected [47]. The FMask algorithm [48] was used in multispectral instrument (MSI) data processed to identify cloud, cloud shadow, cirrus, and snow/ice observations. The *NDVI* values were calculated for

each image element and a time series *NDVI* image collection with a spatial resolution of 10 m is generated. The calculation formula is shown in Equation (2):

$$NDVI = \frac{NIR - RED}{NIR + RED} \tag{2}$$

where *RED* is the TOA values of the red band (630–680 nm); *NIR* is the TOA values of the near infrared band (845–885 nm). The *NDVI* takes values in the range of −1 to 1 [49].

To examine the landscape patterns, we considered landscape features in terms of area, density, shape, diversity, and aggregation on class and landscape levels [7,50]. Details of the calculated landscape spatial indices are shown in Supplementary Material (Table S1).

#### *2.3. Mapping Process*

Based on the questionnaire data, the mean values of ratings for soundscape and landscape perception were calculated for each observation point at each sampling time period, and the data were subsequently visualized in ArcGIS 10.7. By comparing the different interpolation methods provided in spatial analyst tools, the Inverse Distance Weighted (IDW) method was chosen to produce the soundscape and landscape perception maps to analyze their spatiotemporal dynamic characteristics. The IDW method is based on the spatial distance of the data points for weighted interpolation, and the closer the point is to the value, the greater the effect is. Conversely, the smaller the effect is [7,20,21].

The landscape index visualization is based on the moving window technique in Fragstats 4.2. The size and shape of the window can be defined by the user. The window is moved over each cell with positive value in the raster data, and the selected landscape index within the window is calculated. These values are then returned to the focal point (midpoint) of the cell, while a new continuous type of grid data is generated for each selected landscape index, where the cell values represent the "local neighborhood structure" [51]. According to the previous study [7], we set 175 m radius as the window size. In addition, we considered that Fragstats software gives negative values for cells near the edges and cells that are not fully included in the input grid window in the calculating process, which may lead to incomplete spatial data. Therefore, we created a buffer of the same size as the moving window (175 m) before inputting the grid data, thus minimizing the effect of boundary effects [52].

#### *2.4. Statistical Analysis*

(1) Principal Component Analysis (PCA). Based on the semantic attributes of the soundscape, the PCA method was applied to extract the principal components so as to obtain the determinants of overall soundscape quality. The eigenvalues of the extracted principal components were all greater than 1. The analysis was performed in SPSS 25.0.

(2) Multicollinearity Analysis. Before constructing a spatial regression model, it should be ensured that there is no multicollinearity problem among the independent variables included in the model [53]. The tolerance (TOL) and variance inflation factor (VIF) were used to perform the diagnosis of multicollinearity problems. TOL > 0.1 or VIF < 5 indicates that there is no multicollinearity problem between the analyzed [54]. Equations (3) and (4) were used to calculate TOL and VIF, respectively:

$$TOL = \frac{1}{VIF} \tag{3}$$

$$VIF = \frac{1}{1 - R\_j^2} \tag{4}$$

where *R* 2 *j* is the coefficient of determination for the regression analysis on all other variables.

(3) Spatial Regression Model. Spatial regression models capturing spatial dependence were used to examine the effects of perceived sound sources, landscape satisfaction degree, and landscape characteristics on perceived soundscape quality within Kulangsu Island.

Both the global spatial regression model—ordinary least squares (OLS) estimation model, and the local spatial regression model—geographically weighted regression (GWR) model were used (see Supplementary Material for details) [55], and the formulae of them are shown in Equations (5) and (6), respectively [56,57]. The explanatory ability and goodness of fit of these models were examined by the coefficient of determination (*R* 2 ) and the Akaike Information Criterion (AIC) [58,59]. All operations were performed in ArcGIS 10.7.

$$\mathbf{Y}\_{i} = \beta\_{0} + \beta\_{1}\mathbf{x}\_{1} + \beta\_{2}\mathbf{x}\_{2} + \dots + \beta\_{n}\mathbf{x}\_{n} + \varepsilon\_{i}\frac{1}{2} \tag{5}$$

where *Y<sup>i</sup>* is the dependent variable being explained; *x*1, *x*2. . . *x<sup>n</sup>* are the independent variable; *β*<sup>0</sup> is constant; *β*1, *β*2. . . *β<sup>n</sup>* are variable coefficients; *ε<sup>i</sup>* is the bias in estimating the coefficients.

$$Y\_i = \beta\_0(\mu\_{i\prime} v\_i) + \sum\_{k=1}^p \beta\_k(\mu\_{i\prime} v\_i) x\_{ik} + \varepsilon\_i(i = 1, 2, \dots, n) \tag{6}$$

where *Y<sup>i</sup>* is the explained dependent variable, *xik* is the independent variable at *k*th sample of *I*, (*µ<sup>i</sup>* , *v<sup>i</sup>* ) denotes the coordinates of the *i*th sample, *β<sup>0</sup>* (*µ<sup>i</sup>* , *v<sup>i</sup>* ) denotes the intercept of the *i*th sample, *β<sup>k</sup>* (*µi* , *v<sup>i</sup>* ) denotes the regression parameters of the *i*th sample, and ε*<sup>i</sup>* denotes the residuals of the model at the *i*th sample in estimating the coefficients. All data were normalized to 0 to 1 prior to analysis.

#### **3. Results**

#### *3.1. Spatiotemporal Dynamics of Soundscape and Landscape Perceptions*

#### 3.1.1. Soundscape Mapping

Figure 2 shows significant spatial and temporal variability in SDD for all four major sound sources. The SDD-GS is higher in the eastern and western parts, and relatively lower in the central part of the study area. The spatial distribution of SDD for BS, GS, and HS all changed significantly over time. The SDD-HS changed most significantly, and its high values covered almost entire area in the evening. However, the whole area was mainly dominated by BS in the morning. Compared with the SDD of the BS, GS, and HS, SDD-MS only changed slightly over time, and the overall distribution pattern was relatively stable. For the temporal variation of the mean values across land use types, the maximum mean values of SDD-BS in all periods were in the logistics and warehouse land. The mean values of SDD-GS were the always highest in development land, and the lowest in municipal land and logistics and warehouse land. The highest mean values of SDD-HS were found in logistics and warehouse land in P1 and P3, respectively, but in municipal utility land in P2. For SDD-MS, the highest mean values occurred in different land types in each period. The general trend of it was more diverse than the other three sounds.

PCA was used to extract the principal components of the six semantic attributes based on the varimax rotation method, for representing the soundscape quality. The Kaiser-Mayer-Olkin (KMO) test was 0.839 (KMO > 0.60), and the Bartlett's test of sphericity was 0.000 indicating highly significant (*p* < 0.001), which means that the dataset was suitable for PCA. A total of two principal components were extracted: Component 1 (Pleasantness) and Component 2 (Eventfulness), which explained 41.8% and 38.8% of the variance in the semantic attribute dataset, respectively (Table 3). The results of this analysis are in good agreement with the two-dimensional model of perceived affective quality of soundscape (Pleasantness-Eventfulness) proposed in the previous study [37]. This model can provide comprehensive information on soundscape perception.

**Figure 2.** Spatiotemporal distribution and temporal variation of the mean values of sound source dominance of the four main category sound sources. **Figure 2.** Spatiotemporal distribution and temporal variation of the mean values of sound source dominance of the four main category sound sources.


PCA was used to extract the principal components of the six semantic attributes based on the varimax rotation method, for representing the soundscape quality. The Kai-**Table 3.** Rotated component matrices of the PCA based on semantic attributes (numbers in parentheses represent explained variance).

theses represent explained variance). **Semantic Attribute Component 1: Pleasantness (41.8%) Component 2: Eventfulness (38.8%)** Pleasant 0.880 0.250 Comfort 0.906 0.236 Harmony 0.739 0.441 Vivid 0.491 0.703 Richness 0.295 0.867 The values of pleasantness and eventfulness for each observing point were calculated in SPSS 25.0, and then visualized via ArcGIS 10.7 (Figure 3). Both components had significant spatiotemporal characteristics, and appeared to be similar in their distribution patterns. In general, the values of pleasantness and eventfulness in the study area showed low and high values in the north and south, respectively, in all periods. Regarding mean values of the components in different periods for each land use type, the trends of pleasantness and eventfulness were similar. The maximum and minimum values occurred in the same land use type in all periods. The maximum values were found in development land all the time. The minimum values were found in municipal utility land in P1, and in logistics

**Table 3.** Rotated component matrices of the PCA based on semantic attributes (numbers in paren-

Eventful 0.195 0.877

and warehouse land in P2 and P3. The mean value of pleasantness first reduced and then improved in logistics and warehouse land, nonetheless, it had opposite trends in each of the rest land use types. The trend of eventfulness was different from that of pleasantness only in commercial and business facility land, which continued to decline over time. logistics and warehouse land in P2 and P3. The mean value of pleasantness first reduced and then improved in logistics and warehouse land, nonetheless, it had opposite trends in each of the rest land use types. The trend of eventfulness was different from that of pleasantness only in commercial and business facility land, which continued to decline over time.

The values of pleasantness and eventfulness for each observing point were calculated in SPSS 25.0, and then visualized via ArcGIS 10.7 (Figure 3). Both components had significant spatiotemporal characteristics, and appeared to be similar in their distribution patterns. In general, the values of pleasantness and eventfulness in the study area showed low and high values in the north and south, respectively, in all periods. Regarding mean values of the components in different periods for each land use type, the trends of pleasantness and eventfulness were similar. The maximum and minimum values occurred in the same land use type in all periods. The maximum values were found in development land all the time. The minimum values were found in municipal utility land in P1, and in

*Forests* **2022**, *13*, x FOR PEER REVIEW 8 of 21

**Figure 3.** Spatiotemporal distribution and temporal variation of the mean values of pleasantness and eventfulness. **Figure 3.** Spatiotemporal distribution and temporal variation of the mean values of pleasantness and eventfulness.

#### 3.1.2. Landscape Satisfaction Degree 3.1.2. Landscape Satisfaction Degree

Figure 4 indicates the spatiotemporal distribution of each landscape satisfaction degree indicator. Values of SNL, SLD, and SHB were relatively high, while that of SSF was significantly lower. The mapping results show that the spatial differences were pronounced for SLD and SHB, respectively, but that for SSF was relatively small. However, the spatial distribution patterns of them changed significantly with time. In P2, high values of SNL, SLD, and SHB covered almost the whole area. Nonetheless in P3, the high Figure 4 indicates the spatiotemporal distribution of each landscape satisfaction degree indicator. Values of SNL, SLD, and SHB were relatively high, while that of SSF was significantly lower. The mapping results show that the spatial differences were pronounced for SLD and SHB, respectively, but that for SSF was relatively small. However, the spatial distribution patterns of them changed significantly with time. In P2, high values of SNL, SLD, and SHB covered almost the whole area. Nonetheless in P3, the high values of the first two weakened in the central and northern areas, but the high values of SHB still covered almost the whole area. The landscape satisfaction degree indicators showed presented the same trend of mean values in some land use types. The mean values of SNL and SLD in all types of land use, except for development land, showed the same trends. The trends of mean values of SHB and SSF were the same in all land use types. Both of their maximum mean values were found in logistics and warehouse land and development land in P2 and P3, respectively. The mean value of SVA increased and then decrease in all land use types, except for commercial and business facility land.

**Figure 4.** Spatiotemporal variation and temporal variation of the mean values of landscape satisfaction degree. **Figure 4.** Spatiotemporal variation and temporal variation of the mean values of landscape satisfaction degree.

values of the first two weakened in the central and northern areas, but the high values of SHB still covered almost the whole area. The landscape satisfaction degree indicators showed presented the same trend of mean values in some land use types. The mean values of SNL and SLD in all types of land use, except for development land, showed the same trends. The trends of mean values of SHB and SSF were the same in all land use types. Both of their maximum mean values were found in logistics and warehouse land and development land in P2 and P3, respectively. The mean value of SVA increased and then

decrease in all land use types, except for commercial and business facility land.

#### *3.2. Spatial Landscape Characteristics 3.2. Spatial Landscape Characteristics*

The DtGS and NDVI of the study area are shown in Figure 5. The result shows that the vast majority of the study area was close to green spaces (Figure 5a). The NDVI values showed low values in the central area (Figure 5b). The low NDVI values were also found in a few marginal areas of the study area. However, the NDVI values were generally at a high level, with more than 50% of the areas having the values over 0.5 [60], indicating a high vegetation coverage in Kulangsu Island. The DtGS and NDVI of the study area are shown in Figure 5. The result shows that the vast majority of the study area was close to green spaces (Figure 5a). The NDVI values showed low values in the central area (Figure 5b). The low NDVI values were also found in a few marginal areas of the study area. However, the NDVI values were generally at a high level, with more than 50% of the areas having the values over 0.5 [60], indicating a high vegetation coverage in Kulangsu Island. *Forests* **2022**, *13*, x FOR PEER REVIEW 10 of 21

**Figure 5.** The spatial features of green spaces: (**a**) distance to green spaces (DtGS) and (**b**) normalized difference vegetation index (NDVI). **Figure 5.** The spatial features of green spaces: (**a**) distance to green spaces (DtGS) and (**b**) normalized difference vegetation index (NDVI).

COHESION were found mainly in the central and southwestern areas.

Figure 6 shows the spatial distribution of the landscape indices at class level and landscape level. Except for SPLIT at the class level, all other landscape indices had signif-

at similar locations for CA and ED, but the latter covered a wider area. The high values of PD occurred only in a few areas near the northern and southern edges. The high values of IJI covered almost the entire area, while MESH and SPLIT showed lower levels within the area. For the landscape pattern of whole area, the high values of PD concentrated in the central and northern areas, and the high values of LSI were mainly found in the centraleastern area. The high values of SHDI were distributed similarly to those of IJI, which mainly gathered in the periphery of the study area, while the middle area showed a band of low values. SPLIT had a large area of low values in the area, and the high values of

Figure 6 shows the spatial distribution of the landscape indices at class level and landscape level. Except for SPLIT at the class level, all other landscape indices had significant spatial differences. For the landscape pattern of greenspaces, high values occurred at similar locations for CA and ED, but the latter covered a wider area. The high values of PD occurred only in a few areas near the northern and southern edges. The high values of IJI covered almost the entire area, while MESH and SPLIT showed lower levels within the area. For the landscape pattern of whole area, the high values of PD concentrated in the central and northern areas, and the high values of LSI were mainly found in the central-eastern area. The high values of SHDI were distributed similarly to those of IJI, which mainly gathered in the periphery of the study area, while the middle area showed a band of low values. SPLIT had a large area of low values in the area, and the high values of COHESION were found mainly in the central and southwestern areas. *Forests* **2022**, *13*, x FOR PEER REVIEW 11 of 21

**Figure 6.** Landscape spatial patterns in (**a**) class metrics and (**b**) landscape metrics. **Figure 6.** Landscape spatial patterns in (**a**) class metrics and (**b**) landscape metrics.

#### *3.3. Spatial Regression Models of Soundscape Quality 3.3. Spatial Regression Models of Soundscape Quality*

3.3.1. Multicollinearity Diagnostic Results

3.3.2. Global Spatial Regression Model

Spatial regression models, including OLS and GWR models, were built to examine the spatial dependence of soundscape quality (dependent variables) in Kulangsu Island. The explanatory (independent) variables include SDD of 4 main sound sources, landscape satisfaction degree, DtGS, NDVI, and landscape pattern indices. Spatial regression models, including OLS and GWR models, were built to examine the spatial dependence of soundscape quality (dependent variables) in Kulangsu Island. The explanatory (independent) variables include SDD of 4 main sound sources, landscape satisfaction degree, DtGS, NDVI, and landscape pattern indices.

To avoid possible covariance problems among the independent variables, we performed covariance diagnosis in the independent variables at the three periods and the whole day (Total), respectively, and selected the variables without multicollinearity in

The OLS model used the data on each grid to build relationships between the independent and dependent variables (Table 4). None of the landscape indices significantly influenced the soundscape quality at all periods. Only the variables related to SDD and landscape satisfaction degree were significant for the soundscape quality. Pleasantness was only influenced by SLD in P3, and eventfulness was only significantly influenced by SLD in both P2 and P3. The included independent variables had better performances in

in Table S2, Supplementary Material. Most variables passed the multicollinearity test, indicating rationality of the indicator selection. This analysis was performed in SPSS 25.0.

#### 3.3.1. Multicollinearity Diagnostic Results

To avoid possible covariance problems among the independent variables, we performed covariance diagnosis in the independent variables at the three periods and the whole day (Total), respectively, and selected the variables without multicollinearity in each period for the OLS regression. The results of multicollinearity diagnosis are shown in Table S2, Supplementary Material. Most variables passed the multicollinearity test, indicating rationality of the indicator selection. This analysis was performed in SPSS 25.0.

## 3.3.2. Global Spatial Regression Model

The OLS model used the data on each grid to build relationships between the independent and dependent variables (Table 4). None of the landscape indices significantly influenced the soundscape quality at all periods. Only the variables related to SDD and landscape satisfaction degree were significant for the soundscape quality. Pleasantness was only influenced by SLD in P3, and eventfulness was only significantly influenced by SLD in both P2 and P3. The included independent variables had better performances in explaining the pleasantness (*R* <sup>2</sup> was 0.653 and 0.649, respectively) and eventfulness (*R* <sup>2</sup> was 0.591 for both) in P1 and Total, compared with those in P2 and P3. Similarly, the results of the AIC demonstrate that the OLS models in P1 and Total had better fitting accuracy than those in P2 and P3.


**Table 4.** Results of the OLS model (\* *p* < 0.05, \*\* *p* < 0.01).

#### 3.3.3. Local Spatial Regression Model

The GWR model was to explore the spatial relationships between soundscape quality and the significant variables derived from the OLS model (Table 4). Before executing the GWR model, a kernel function needs to be selected for performing the geo-weighting algorithm, which is to estimate the local coefficients and their bandwidth sizes. There are two commonly used functions in GIS, namely the Gaussian fixed and adaptive kernel type. Based on the finding of previous study [35], we selected the Gaussian fixed kernel type to perform the GWR model.

The values of AIC of the GWR models in all period were significantly smaller than those of the OLS models (Table 5), indicating that the GWR model had better fitting results than the OLS model. The mean values of *R* <sup>2</sup> of the GWR models were higher than those of the OLS models for the all periods, nonetheless, except for the pleasantness models in P1 and Total. As shown in Figure 7, the *R* 2 spatial distributions of pleasantness model

in P1, and both pleasantness and eventfulness models in Total were relatively stationary compared with those in P2 and P3. Also, Table 5 indicates that the results of *R* <sup>2</sup> and AIC in P1 and Total were better than those in P2 and P3. These findings indicate that such spatial stationarity may achieve a better performance of the statistical model in terms of explanation of variance and goodness-of-fit. *Forests* **2022**, *13*, x FOR PEER REVIEW 13 of 21 **Table 5.** *R*<sup>2</sup> and AIC of the GWR model.



Note: To account for differences, the maximum and minimum values of the *R* <sup>2</sup> of the GWR models for the total dataset were retained to five decimal places. Note: To account for differences, the maximum and minimum values of the *R*<sup>2</sup> of the GWR models for the total dataset were retained to five decimal places.

**Figure 7.** Spatial distribution of local *R*<sup>2</sup> for pleasantness and eventfulness in the three and total sampling periods. **Figure 7.** Spatial distribution of local *R* 2 for pleasantness and eventfulness in the three and total sampling periods.

The spatial distribution of the coefficients of each independent variable are shown in Figure 8, the spatial stationarity also presented in pleasantness model in P1, and both pleasantness and eventfulness models in Total. In P1, SDD-MS and SVA exhibited the strongest negative and positive effects, respectively, on both pleasantness and eventfulness. In P2, SLD was positively related to both pleasantness and eventfulness in most areas, but their negative relationships appeared in few areas in the northwest. In addition, pleasantness was also positively influenced by SDD-GS, and the spatial differences were significant with the coefficients decreasing from northwest to southeast (from 0.808 to 0.139). In P3, both pleasantness and eventfulness were positively related to SLD only, with The spatial distribution of the coefficients of each independent variable are shown in Figure 8, the spatial stationarity also presented in pleasantness model in P1, and both pleasantness and eventfulness models in Total. In P1, SDD-MS and SVA exhibited the strongest negative and positive effects, respectively, on both pleasantness and eventfulness. In P2, SLD was positively related to both pleasantness and eventfulness in most areas, but their negative relationships appeared in few areas in the northwest. In addition, pleasantness was also positively influenced by SDD-GS, and the spatial differences were significant with the coefficients decreasing from northwest to southeast (from 0.808 to

the coefficients increasing from southwest to northeast. In the model of Total, all included

0.139). In P3, both pleasantness and eventfulness were positively related to SLD only, with the coefficients increasing from southwest to northeast. In the model of Total, all included variables had positive relationships with pleasantness and eventfulness, respectively. SLD was the variable with strongest effects on pleasantness and eventfulness, respectively. *Forests* **2022**, *13*, x FOR PEER REVIEW 14 of 21

**Figure 8.** Spatial distribution of local coefficients of significant variables in the three and total sampling period. **Figure 8.** Spatial distribution of local coefficients of significant variables in the three and total sampling period.

#### **4. Discussion 4. Discussion**

#### *4.1. Spatiotemporal Dynamics of Soundscape and Landscape in Kulangsu 4.1. Spatiotemporal Dynamics of Soundscape and Landscape in Kulangsu*

The results of the present study showed that both soundscape and landscape perception in Kulangsu Island were characterized by significant spatiotemporal dynamics. In terms of perceived sound sources, we calculated and visualized the SDD for each main sound source, to present the dominance of them in different temporal and spatial dimensions in Kulangsu (Figure 2). The results demonstrated that BS and HS were the dominant sounds during the three sampling periods. The dominance of BS may imply the good quality of ecological environment within Kulangsu, which can attract various vocal organisms such as birds and insects to congregate and communicate here [61]. This was also reflected the high level of biodiversity in the island. Such abundance of natural resources also provides a good opportunity for tourism development, with many urban residents and tourists choosing to visit the area for e.g., nature-based recreation [62]. We speculate that this maybe one of the reasons why HS was another dominant sound. The results of the present study showed that both soundscape and landscape perception in Kulangsu Island were characterized by significant spatiotemporal dynamics. In terms of perceived sound sources, we calculated and visualized the SDD for each main sound source, to present the dominance of them in different temporal and spatial dimensions in Kulangsu (Figure 2). The results demonstrated that BS and HS were the dominant sounds during the three sampling periods. The dominance of BS may imply the good quality of ecological environment within Kulangsu, which can attract various vocal organisms such as birds and insects to congregate and communicate here [61]. This was also reflected the high level of biodiversity in the island. Such abundance of natural resources also provides a good opportunity for tourism development, with many urban residents and tourists choosing to visit the area for e.g., nature-based recreation [62]. We speculate that this maybe one of the reasons why HS was another dominant sound.

We found that the HS-dominated area gradually expanded over time, and covered almost the entire area in the evening (Figure 2), indicating that this period was the most intensive period for people's activities. This survey was conducted during the summer of July, which is generally the hottest period of the year in Xiamen [63]. Accordingly, we speculate that one reason for this phenomenon may be due to the higher thermal comfort in the evening, compared to the built-up area. This is probably due to the good vegetation condition in Kulangsu, which has a significant effect on reducing temperature thus creating a comfortable microclimate. Generally, such cooling effects were most pronounced in the evening [64]. Another reason may be the development of the night-time economy in China in recent years, which has led to a new urban living pattern dominated by nighttime entertainment activities [65]. We found that the HS-dominated area gradually expanded over time, and covered almost the entire area in the evening (Figure 2), indicating that this period was the most intensive period for people's activities. This survey was conducted during the summer of July, which is generally the hottest period of the year in Xiamen [63]. Accordingly, we speculate that one reason for this phenomenon may be due to the higher thermal comfort in the evening, compared to the built-up area. This is probably due to the good vegetation condition in Kulangsu, which has a significant effect on reducing temperature thus creating a comfortable microclimate. Generally, such cooling effects were most pronounced in the evening [64]. Another reason may be the development of the night-time economy in China in recent years, which has led to a new urban living pattern dominated by night-time entertainment activities [65].

In the analysis of SDD mean values of the sound sources for different land use types

small in size and surrounded by large areas of green spaces. This also potentially reflects

In the analysis of SDD mean values of the sound sources for different land use types (Figure 2), high SDD of BS and HS unexpectedly appeared several times in logistics and warehouse land. We suspect that this may be because this type of land use was generally small in size and surrounded by large areas of green spaces. This also potentially reflects that both animals and visitors were more willing to get close to the natural landscapes. The high GS-SDD mean values were mainly found in development land. This may be due to the harbors and waterfront spaces that were close to the sea in such land use type and therefore the GS especially water and wave sounds were more dominant. Such results are similar to the findings of previous studies located in coastal areas [66].

A similarity was found in the spatiotemporal distribution patterns of pleasantness and eventfulness of the soundscape (Figure 3), which implies they may be potentially interrelated. This finding is the same as previous soundscape studies located in urban built-up areas and urban forests [20,22]. For mean values of pleasantness and eventfulness across land use types (Figure 3), the maximum values of them were always found in development land, and mainly concentrated in the south according to the mapping result. This is possibly because the development land to the south were close to large areas of green spaces and commercial lands, and therefore subject to the "radiation effect" from such land types. The development land, on the other hand, were relatively smaller and thus may receive higher aggregated values of pleasantness and eventfulness. Another possible reason was that the development land encompassed many waterfront spaces where people prefer to stay, to experience nature and response higher values of pleasantness. Moreover, the consequent increase in foot traffic may also lead to high values of eventfulness [67].

As for landscape satisfaction degree, we found that the spatiotemporal distribution patterns of SNL, SLD, SHB, and SVA were similar to those of pleasantness and eventfulness in P1, but they were relatively different in P2 and P3 (Figures 3 and 4). This finding suggests that the harmony or congruency between the soundscapes and landscapes in Kulangsu was still deficient especially in the afternoon and evening, and therefore it has room to improve in future landscape planning and design. We speculate that this audio-visual inconsistency may be due to the increasing spatial dominance of HS in these periods (Figure 2), which negatively affected the audio-visual perception [68]. Regarding the spatiotemporal distribution of SSF, it was significantly lower compared to the other landscape perception indicators. However, the distribution pattern of SSF was also similar to those of pleasantness and eventfulness in the morning and afternoon (Figures 3 and 4). This similarity indicates that people may prefer to use these service facilities especially in these two periods. We suppose that this is because the temperatures in daytime are generally high in Kulangsu. This finding is also in line with the results of the previous study, which found infrastructure and facilities were an important factor contributing to the soundscape quality in scenic areas [27].

The analysis of "objective" landscape characteristics also illustrates the superior natural environmental characteristics of Kulangsu (Figures 5 and 6), which was consistent with the finding derived from the spatiotemporal distribution of SDD-BS (Figure 2). These results indicated that Kulangsu had good green space accessibility and a high proportion of vegetation cover, which allows visitors to easily approach and experience nature and contributes to the comfortable microclimate here [60]. We found that although the high values of SDD-BS appeared in different areas over time, they almost always fell in the areas with low DtGS or high NDVI, indicating the interrelationship between biological sounds and green space features. Regarding the landscape pattern, the results of class metrics indicate that the green space patches were large and morphologically intact, with a high degree of connectivity, and did not have significant fragmentation. In addition to SPLIT, the spatial distributions of CA, PD, ED, IJI, and MESH were more or less similar to those of pleasantness and eventfulness in the three periods, suggesting that diverse and coherent structure of green space patches may contribute to soundscape quality [69]. Based on the results of the landscape metrics, we found that the structure of the patches in the center of Kulangsu was more complex and homogeneous, and the connectivity between patches was

good. This may be due to the high density of historical architectures in the central area [70], but such buildings still maintain a good continuity. The edge of the study area, on the other hand, had a richer patch composition and different types of patches may be interspersed between similar patches. We found that this complexity tends to accompany areas where commercial land occurred. Based on such, we suggest that the remediation planning of scenic areas could focus on preventing the fragmentation effects on the landscape especially from dispersed commercial activities [71].

#### *4.2. Spatial Dependencies of Soundscape Quality*

Most indicators with regard to perceived sound sources, landscape satisfaction degree, the features of green spaces, and landscape patterns passed the test of multicollinearity on explaining the soundscape quality (Table S2), and then included in the OLS and GWR models. Interestingly, the results of OLS models (Table 4) show that the soundscape quality had no spatial dependence on all objective landscape characteristics, namely DtGS, NDVI, and landscape metrics on class and landscape levels. Both pleasantness and eventfulness of the soundscape were significantly dependent on the perceptual indicators of soundscape and landscape, i.e., SDD-BS, SDD-GS, SDD-MS, SLD, SHB, and SSF, in different periods. This finding was not exactly the same as previous studies on the relationship between soundscape elements and landscape patterns [7,72]. We speculate that this discrepancy may arise from two aspects. On the one hand, the research objectives and methods were different. The present study aimed to explore the spatial dependence of soundscape quality on many types of variables rather than the correlations between two variables. The spatial dependence can indicate the propensity for nearby locations to interact and share similar properties, which is also an essential part of modeling soundscape [35]. Therefore, the regression model accounting for the spatial dependence can explain more about the relationships between the dependent and independent variables, and predict the dependent variable based on such, while the correlation does not necessarily imply causation [73]. On the other hand, the study area was different. Our study was devoted to exploring the soundscape and landscape resources of the Kulangsu scenic area, whereas the previous studies were in a multifunctional urban area [7,20,72]. Kulangsu Scenic Area is rich in natural and cultural resources and locally distinctive sound sources, which is basically different from the urban area.

The result of soundscape spatial dependency on SDD-BS, SDD-GS, and SDD-MS suggests that introducing more birdsongs and vegetation sounds, through enriching the vegetation types and density, may contribute to the biodiversity and reduction of mechanical noise [74], therefore improving the pleasure and structure diversity of the soundscape [37]. Besides, according to the dependence on SLD and SHD, the soundscape quality can also be enhanced by optimizing the spatiotemporal features of landscape design, with regard to the lightscape (especially in the evening) and cultural innovation [24,75]. For instance, adding interesting lightscapes designed according to the natural and cultural environment characteristics of Kulangsu. The type and application of the light sources should be able to integrate into but not affect the natural environment [76]. Furthermore, the historical culture of Kulangsu can be combined with modern elements to revitalize old buildings by injecting new functions, such as cultural and creative business or exhibitions [77]. Also, the "Piano Island", one of the famous titles of Kulangsu, can be used to create e.g., a theme-specific "musical environment" in different public open spaces by combining classical and contemporary featured piano music. The significant dependency of soundscape quality on the landscape perception indicators implies that only considering the effects of sound features on soundscape quality was insufficient in the past study [35]. This result also proves that the landscape compositions are indispensable factors in creating a soundscape [17,43]. The temporal differences of spatial dependence of the soundscape in Kulangsu can contribute specific information to the soundscape planning and design strategies for scenic areas and Heritage Site.

The GWR model was found to have a better goodness-of-fit than the OLS model for all periods according to the results of AIC, which was consistent with the previous study [35]. However, the *R* 2 results of GWR model performing pleasantness in P1 and Total were lower than those of OLS model (Tables 4 and 5). Interestingly, the spatially stationary relationships were found in the results of *R* <sup>2</sup> and coefficients of the variables in pleasantness model in the morning and both pleasantness and eventfulness models in total dataset, respectively (Figures 7 and 8). This finding suggests that the pleasantness models had better fitting accuracy if they were spatially stationary, but the explanatory power of the regression line for the dependent variable became weaker, which implies that the nonlinear relationship may be able to better explain the spatial dependence of soundscape pleasantness. This is also a good proof that the spatial stationarity helps to verify and select the suitable statistical methods in predicting variables [78,79]. As shown in Figures 7 and 8, the total dataset models of pleasantness and eventfulness exhibited the best spatial stationarity, indicating the data used should ideally be of a time-varying nature that may contain as complete a temporal dimension as possible. For instance, the collected data should ideally cover morning, afternoon, evening, and even in the night or early morning.

#### *4.3. Limitations and Future Studies*

The limitations of the present study may mainly stem from the data collected and the model employed. Although vegetation coverage was considered in this study, the classification of green spaces was not in sufficient detail in terms of the vegetation type. The outcomes may therefore vary due to more detailed classification of vegetation cover. However, a high spatial resolution vegetation type data in China is difficult to obtain in general. Future studies are suggested to utilize data regarding vegetation cover types, such as coniferous forests, deciduous forests, shrubs, and grasslands, to explore the spatial relationship between their features and soundscape quality.

Moreover, the questionnaire data collected in 2016 was not up-to-date, nonetheless such data was still of research interest and significance, because this investigation can provide valuable historical research data for further exploring the environmental changes and soundscape planning and management. Given this importance, we suggest that future research could explore the spatiotemporal evolution of soundscape quality in Kulangsu in conjunction with the results of the present study. Based on e.g., the concept of DPSIR model [80], it is possible to predict the future state of soundscape resources in landscape planning.

The spatial model used in this study was based on a linear regression algorithm. Therefore, although the results indicated that landscape characteristics did not have significant effects on soundscape quality, this does not account for the possible non-linear relationships between these variables [81]. It has been proposed that the non-linear models may be able to explain soundscape perception better than the linear models [82]. However, the nonlinear prediction models related to soundscape are still in their infancy, and most of them are only able to predict the values of output variables without spatiotemporal dynamic features [22,83,84]. Based on such, we recommend future studies could use the non-linear models to explore the soundscape spatiotemporal properties and to predict soundscape perception based on such. The models, such as the cellular automata, artificial neural networks, and Fuzzy-logic models [85,86]. Although the linear model is not perfect, the spatial regression model used in this study still presented a relatively better performance in spatially explaining the soundscape quality than the conventional regression model [20].

#### **5. Conclusions**

This study provided indispensable data with regard to the spatiotemporal dynamics of soundscape and landscape perception as well as the spatial landscape characteristics in Kulangsu. Furthermore, the spatial dependence of soundscape quality was examined. We found that sound dominant degree, soundscape quality, and landscape satisfaction degree

all showed high spatiotemporal dynamics, and the distance to greenspaces, NDVI, and landscape patterns presented obviously spatial variations. Moreover, soundscape quality had spatial dependence on the sound dominant degree and landscape satisfaction degree, but not on the objective landscape characteristics. Such spatial dependence of soundscape quality provides useful suggestions for the design of soundscape and landscape. The GWR model had better goodness-of-fit than the OLS model, and the results with spatial stationarity of the GWR models suggest that applied data should consider as complete a time dimension as possible in exploring the spatial dependence of soundscape quality, which also guidance for soundscape modeling. The findings of this study will allow landscape planners to determine what factors the quality of a soundscape depends on, which gives a relatively comprehensive consideration of selected indicators for modeling soundscapes. They are also meaningful for the bringing soundscape evaluation into planning and design practices, especially for the development of specific planning goals and design strategies, which aim to protect, restore, and optimize the soundscape resources in the scenic areas and World Heritage Sites.

**Supplementary Materials:** The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/f13091526/s1, Figure S1: Social/demographical/behavioral information of the interviewees; Table S1: Calculated landscape spatial indices at class and landscape levels; Table S2: Results of multicollinearity diagnosis for pleasantness and eventfulness at each sampling period; detailed concepts and definition in the statistical analysis. Reference [87] is cited in the supplementary materials.

**Author Contributions:** Conceptualization, Z.C., T.-Y.Z. and J.L.; methodology, Z.C. and J.L.; Software, Z.C. and T.-Y.Z.; validation, Z.C., T.-Y.Z., J.L. and X.-C.H.; formal analysis, Z.C. and T.-Y.Z.; investigation, Z.C. and T.-Y.Z.; resources, Z.C., T.-Y.Z. and J.L.; data curation, Z.C. and T.-Y.Z.; writing—original draft preparation, Z.C.; writing—review and editing, Z.C. and J.L.; visualization, Z.C. and T.-Y.Z.; supervision, J.L.; project administration, Z.C., J.L. and X.-C.H.; funding acquisition, J.L. and X.-C.H. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the National Natural Science Foundation of China (51508101 and 52208052), Fujian Provincial Department of Science & Technology (2017J01694), and the Program of Humanities and Social Science Research Program of Ministry of Education of China (21YJCZH038).

**Data Availability Statement:** Data available on request due to restrictions, e.g., privacy or ethical.

**Acknowledgments:** The first author, Zhu Chen, would like to thank the China Scholarship Council for the support of a doctoral scholarship (grant number: 202108080105). All the authors would like to thank the anonymous reviewers for their valuable comments, which have greatly improved the quality of this paper. We thank Christina von Haaren and Johannes Hermes for their very helpful suggestions in this paper.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


MDPI St. Alban-Anlage 66 4052 Basel Switzerland Tel. +41 61 683 77 34 Fax +41 61 302 89 18 www.mdpi.com

*Forests* Editorial Office E-mail: forests@mdpi.com www.mdpi.com/journal/forests

MDPI St. Alban-Anlage 66 4052 Basel Switzerland

Tel: +41 61 683 77 34 Fax: +41 61 302 89 18

www.mdpi.com ISBN 978-3-0365-6273-5