Comprehensive Evaluation of Healthcare Beneﬁts of Different Forest Types: A Case Study in Shimen National Forest Park, China

: Studies have shown that contact with nature plays a crucial role in the amelioration of human health. Forest therapy has recently received widespread attention as a novel and subsidiary treatment approach for stress recovery and health promotion. However, there is a lack of ample research on the comprehensive evaluation of the forest healthcare beneﬁts. Moreover, it is not entirely clear what kind of forest types and seasons are suitable for forest therapy activities and how healthcare forests should be constructed and managed. From September 2019 to January 2020 and May to August 2020., ﬁve forest types of Phyllostachys edulis forest, subtropical evergreen broad-leaved forest, Liquidambar formosana forest, Cunninghamia lanceolata forest, coniferous and broad-leaved mixed forest and a forestless control group in Shimen National Forest Park, Guangzhou City, Guangdong Province, China were selected. Variations in the character of negative air ion concentration, air oxygen content, human comfort index and phytoncide relative content were analyzed. Principal component analysis and systematic clustering were used to construct forest comprehensive healthcare index and evaluation grade in order to assess the healthcare beneﬁts of different forest types. In terms of negative air ion concentration, the subtropical evergreen broad-leaved forest was far ahead of the other forest types throughout the year, while the forestless control group was the worst. All stands reached the annual maximum in summer, followed by spring, autumn and winter. From the perspective of air oxygen content, summer > spring > autumn > winter, among them, all forest stands clearly exceeded the normal atmospheric oxygen content (20.9%) in the ﬁrst three quarters. Moreover, the air oxygen content of coniferous and broad-leaved mixed forest was the highest in ﬁve forest types; the forestless control group was the lowest. Judging from the human comfort index, in the whole year, all forest types, including the forestless group, were at the comfortable level and above. However, the ﬁve forest types still differed greatly in diverse seasons, among which Phyllostachys edulis


Introduction
In recent years, due to environmental contamination combined with intense social competition and heavier life pressure, a growing number of people are affected by lifestylerelated diseases and are in a sub-healthy state.Moreover, medical costs are increasing at home and abroad [1,2].Increasingly, people have reached a consensus on self-care.At the same time, traveling for health purposes is a global tourism trend.Numerous countries have identified health, medical, and wellness tourism as a key development area for the near or distant future.China has significant resources, as well as plans for developing health tourism [3].Furthermore, according to surveys, most people realize that the forest environment has a favorable influence on human physical and mental health [4].There has been increasing interest in forest therapy, defined as making contact with nature and taking in the atmosphere of the forest, as an environmentally friendly therapeutic approach for improving an individual's mental and physical relaxation and relieving stress [5].Studies conducted over the past twenty years have systematically illustrated the effects of forest therapy on human health, as well [6][7][8][9].Therefore, with the development of health tourism and the enhancement of people's awareness of healthcare, forest healthcare resources and function evaluation issue have been receiving unprecedented attention.
Forest healthcare resources are defined as all biological and abiotic factors with healthcare benefits in the forest environment, including negative air ions, phytoncide, air oxygen content, forest microclimate advantages, acoustic environment advantages, surface water environment, natural radiation level, biological population advantage, forest food and other forest healthcare factors.Negative air ions-the general names of negatively charged single gas molecules and light ion clusters in the atmosphere-are a vital indicator of forest healthcare benefits and are significant for the evaluation of air quality of different forest types.Early research has shown the beneficial effects of negative air ions in improving neuropsychological performance and treating mood disorders [10], alleviating depression [11], promoting sleep quality [12].In addition, under different influences of tree species, structures, ages, densities, heights, and so on, negative air ion concentrations of diverse forest types are significantly various.Phytoncides, first coined in 1928 by Boris P. Tokin, a Russian biochemist, are defined as natural volatile compounds derived from trees and plants to defend against decay or attack by herbivores [13].Previous research findings strongly suggest that tree-derived phytoncides have beneficial effects on antifungal, anti-inflammatory, antimicrobial, analgesic, anti-stress, antioxidant [14][15][16][17][18][19], human immune functions [20], even sleep quality improvement [13,21].The forest air has high oxygen content.Oxygen is our primary life support.The air we breathe is so vital that hypoxia would cause fatigue, memory loss, inattention, insomnia and other symptoms [22].Regular forest bathing activities can improve the sub-healthy status and make physical functions, intelligence and blood oxygen concentration reach the best states.Forest microclimate advantages refer to the human body's perception of the forest environment's temperature, relative humidity, wind speed and other factors, which is often measured by a human comfort index.The potential of the forest microclimate environment in regulating temperature is widely known, which is conducive to human leisure and healthcare [23,24].At the same time, evaluation of the human comfort index is one of the bases of the con-Forests 2021, 12, 207 3 of 26 struction of the scientific human settlement environment.There have been many studies on the evaluation and analysis of single-forest healthcare factors of negative oxygen ions, phytoncide, human comfort and the natural soundscape.However, the forest healthcare benefits reflected by a single index are one-sided; there have been few studies that have focused on evaluating forest healthcare benefits by multiple comprehensive indices.The aim of this study, in which the healthcare benefits of multiple environmental factors are considered comprehensively, therefore, was mainly to investigate three scientific problems as follows: (1) How can we comprehensively evaluate the healthcare benefits of diverse forest types?(2) Which forest types are suitable for people to carry out forest therapy activities?(3) How do we construct a healthcare forest, improve its healthcare benefits and optimize forest structure in a later period?In the article Suggestions on Promoting the Development of Forest Therapy Industry in 2019, The Chinese government mentioned that forest therapy bases are necessary to create a forest healthcare environment with obvious effects, targeted construction and improvement of stands with healthcare benefits.Hence, it is of great theoretical and practical significance to explore these problems.Shimen National Forest Park is one of the first batches of forest therapy bases in Guangdong Province, China.In this study, five different forest types and a forestless control group in Shimen National Forest Park were taken as the research objects, and the changes of environmental healthcare indicators such as negative air ion concentration, air oxygen content, forest microclimate (temperature, relative humidity and wind speed) and phytoncide were monitored.By using principal component analysis and system clustering method, the forest comprehensive healthcare index (FCHI) and evaluation criteria were constructed.The results will not only provide a theoretical basis for the comprehensive evaluation of the healthcare benefits of different forest types but also offer powerful data support for the future construction and structure optimization of healthcare forest, affording an essential selection basis of time and destination for forest therapy.

Study Area
Shimen National Forest Park is located in the northeast part of Conghua District, Guangzhou City, Guangdong Province, China (23 • 36 50 "N~23 • 39 20" N, 113 • 46 16 "E1 13 • 49 17" E), see Figure 1.It has a total area of 26.36 square kilometers, with a forest coverage rate of 98.91%, which is China's first international forest bathing beach and also one of the first forest therapy tourism bases in Guangdong Province.The landform types of Shimen National Forest Park are mainly middle mountain, low mountain and hill, and its altitude ranges from 270 m to 1210 m.The terrain is high in the southeast, low in the northwest, surrounded by mountains in the north, east and south.Shimen National Forest Park is located at the northern edge of the tropic of cancer, with a south subtropical monsoon climate, which has wet, scorching summers and mild, sunny winters.The average air temperature fluctuates between 19.5 • C to 21.4 • C, where the coldest days are in January with an average temperature of 12.4 • C, and the hottest days are in July with an average temperature of 28.5 • C. Sunshine is plentiful all year;the mean annual solar radiation is approximately 440,870 J/cm 2 .The subtropical monsoon brings abundant rainfall with a long rainy season from May to August, and the annual average rainfall is about 1800 mm.

Plot Settings
There were five forest types and a forestless control group considered in this study; the details are presented in Table 1 and Figure 2. In the forestless control group, we selected three sample points in open spaces with no tree layer, mainly herbs and very close to the water (about 30 m).In the Phyllostachys edulis forest plots, we selected four sample points inside the forest.Its vegetation was mainly Phyllostachys edulis and a small number of herbs.Four sample points were selected inside the subtropical evergreen broad-leaved forest, which was mainly mingled with Machilus breviflora, Machilus Chekiangensis, Machilus chinensis, Diospyros morrisian, Lithocarpus corneus, Michelia maudiae, Sapium sebiferum, Manglietia pachyphylla, Bretschneidera sinensis, etc., with undergrowth of Melastoma candidum, Rhodomyrtus tomentosa, Melastoma dodecandrum, Sarcandra glabra, Lobelia angulata, Senecio scandens, Mucuna birdwoodiana, etc.In the Liquidambar formosana forest plots, we selected three sample points inside the forest.Its tree layer was exclusively Liquidambar formosana, with undergrowth of Lophatherum gracile.In the Cunninghamia lanceolata forest plots, we selected three sample points inside the forest.Its tree layer was exclusively Cunninghamia lanceolata, with undergrowth of Lophatherum gracile, Blechnum orientale, Pteris semipinnata, etc.Three sample points were set up inside the coniferous and broad-leaved mixed forest, in which the tree layer was mainly mingled with Liquidambar formosana (taking up 40%) and Cunninghamia lanceolata (accounting for 60%), with undergrowth of Callicarpae formosanae, Lophatherum gracile, Blechnum orientale, Nephrolepis auriculata, Gleichenia linearis, etc.

Plot Settings
There were five forest types and a forestless control group considered in this study; the details are presented in Table 1 and Figure 2. In the forestless control group, we selected three sample points in open spaces with no tree layer, mainly herbs and very close to the water (about 30 m).In the Phyllostachys edulis forest plots, we selected four sample points inside the forest.Its vegetation was mainly Phyllostachys edulis and a small number of herbs.Four sample points were selected inside the subtropical evergreen broad-leaved forest, which was mainly mingled with Machilus breviflora, Machilus Chekiangensis, Machilus chinensis, Diospyros morrisian, Lithocarpus corneus, Michelia maudiae, Sapium sebiferum, Manglietia pachyphylla, Bretschneidera sinensis, etc., with undergrowth of Melastoma candidum, Rhodomyrtus tomentosa, Melastoma dodecandrum, Sarcandra glabra, Lobelia angulata, Senecio scandens, Mucuna birdwoodiana, etc.In the Liquidambar formosana forest plots, we selected three sample points inside the forest.Its tree layer was exclusively Liquidambar formosana, with undergrowth of Lophatherum gracile.In the Cunninghamia lanceolata forest plots, we selected three sample points inside the forest.Its tree layer was exclusively Cunninghamia lanceolata, with undergrowth of Lophatherum gracile, Blechnum orientale, Pteris semipinnata, etc.Three sample points were set up inside the coniferous and broad-leaved mixed forest, in which the tree layer was mainly mingled with Liquidambar formosana (taking up 40%) and Cunninghamia lanceolata (accounting for 60%), with undergrowth of Callicarpae formosanae, Lophatherum gracile, Blechnum orientale, Nephrolepis auriculata, Gleichenia linearis, etc.

Observation Methods
Observations of sunny and stable weather-related data were conducted from September 2019 to January 2020 and May to August 2020.The field data were measured for three consecutive days at the end of each month.Due to the influence of the COVID-19 pandemic, research data from February to April 2020 are missing.Synchronous observations were performed in the morning (9:00-11:00), noon (12:00-14:00) and afternoon (15:00-17:00) of each test day.Four indices of negative air ion concentration (NAIC), air oxygen content (AOC), forest microclimate and phytoncide relative content (PRC) were measured in the five forested stands and the forestless control group.The negative air ion concentration was investigated by a unified corrected COM-3200PRO II negative ion monitor produced in Japan with an observation range of 10−1.999× 10 6 ions/cm 3 .Its measurement accuracy is ±10%; its mobility is ±10%.The data were measured 1.5 m from the ground, which was basically consistent with the adult breathing height.When the instrument was stable at each observation point, we read the data in four directions (east, south, west, north) and repeated it three times after zeroing.With a total of twelve data points, we obtained the negative air ion concentration by taking the average value.The

Observation Methods
Observations of sunny and stable weather-related data were conducted from September 2019 to January 2020 and May to August 2020.The field data were measured for three consecutive days at the end of each month.Due to the influence of the COVID-19 pandemic, research data from February to April 2020 are missing.Synchronous observations were performed in the morning (9:00-11:00), noon (12:00-14:00) and afternoon (15:00-17:00) of each test day.Four indices of negative air ion concentration (NAIC), air oxygen content (AOC), forest microclimate and phytoncide relative content (PRC) were measured in the five forested stands and the forestless control group.The negative air ion concentration was investigated by a unified corrected COM-3200PRO II negative ion monitor produced in Japan with an observation range of 10−1.999× 10 6 ions/cm 3 .Its measurement accuracy is ±10%; its mobility is ±10%.The data were measured 1.5 m from the ground, which was basically consistent with the adult breathing height.When the instrument was stable at each observation point, we read the data in four directions (east, south, west, north) and repeated it three times after zeroing.With a total of twelve data points, we obtained the negative air ion concentration by taking the average value.The oxygen content was measured by a TD6000-SH-O 2 oxygen content tester produced in China, which was repeated 3 times.Each instrument was corrected before measurement.The forest microclimate was measured with Kestrel 5500 hand-held meteorological instrument produced in the USA; synchronous monitoring data included three environmental factors, namely temperature (T), relative humidity (RH) and wind speed (V), which were each repeated three times.The  oxygen content was measured by a TD6000-SH-O2 oxygen content tester produced in China, which was repeated 3 times.Each instrument was corrected before measurement.
The forest microclimate was measured with Kestrel 5500 hand-held meteorological instrument produced in the USA; synchronous monitoring data included three environmental factors, namely temperature (T), relative humidity (RH) and wind speed (V), which were each repeated three times.The sampling height of wind speed was 2 m.The relative content of phytoncide was tested by the open, dynamic headspace adsorption collection method.The phytoncide collection devices provided by South China Botanical Garden, Chinese Academy of Sciences, which were hung at 1.5 m in the forest at the average breathing height of the human body, needed to work continuously for at least 8 h from 9:00 to 17:00.The experimental instruments mentioned above are shown in Figure 3.

Specific calculation methods
(1) Air oxygen content In order to eliminate the influence of altitude on air oxygen content, the measured value was corrected according to the fact that the air oxygen content decreases by 0.16% every 100 m above sea level [25].The correction formula used was: correction value = measured value + altitude/100 × 0.16 (2) Human comfort index Lu Dinghuang [26] used the relevant data of environmental hygiene methods to comprehensively consider the influence of temperature, relative humidity and wind speed on the human comfort index, which was frequently used.Therefore, the human comfort index proposed by Lu Dinghuang was used in this study, and the formula for calculating the human comfort index was as follows:

Specific Calculation Methods
(1) Air oxygen content In order to eliminate the influence of altitude on air oxygen content, the measured value was corrected according to the fact that the air oxygen content decreases by 0.16% every 100 m above sea level [25].The correction formula used was: correction value = measured value + altitude/100 × 0.16 (2) Human comfort index Lu Dinghuang [26] used the relevant data of environmental hygiene methods to comprehensively consider the influence of temperature, relative humidity and wind speed on the human comfort index, which was frequently used.Therefore, the human comfort index proposed by Lu Dinghuang was used in this study, and the formula for calculating the human comfort index was as follows: In the formula, HCI is the human comfort index, T is the temperature ( • C), RH is relative humidity (%), and Vis wind speed (m/s).The lower the HCI value, the higher the comfort level.Among them, "S ≤ 4.55" denotes a very comfortable level, "4.55 < S ≤ 6.95" represents a comfortable level, "6.95 < S ≤ 9.00" indicates an uncomfortable level, "S > 9. 00" means extremely uncomfortable.
(3) Phytoncide relative content The adsorption tube was properly sealed and preserved after the phytoncide gas was collected, eluted with dichloromethane solution in time, and then analyzed those processed samples with a gas chromatography-mass spectrometry (GC-MS).The test instrument model was Agilent 7890B-5977B/GC-MS, which is manufactured by Agilent in the USA.GC working conditions were as follows: The chromatographic column was HP-5MS elastic quartz capillary column (30 m, 0.25 mm, 0.25 µm).The temperature program for Volatile Organic Chemicals (VOCs) were: initial temperature 50 • C and kept for 2 min, 50 min for 10 min.The carrier gas was helium, and the sample injection volume was 1 µL.The split ratio was 20:1, the inlet temperature was 250 • C, and the split flow was 20 mL/min.MS working conditions were as follows: electron energy was 70 eV, transmission line temperature was 250 • C, and ion source temperature was 230 • C. The scanning mode was full ion scanning, and the quality scanning range was m/z 25-500 u.The total ion flow diagram was obtained by TM Software4.0,and the chemical information represented by each peak was confirmed and screened by searching the NIST standard spectrum library by computer.The volatile organic compounds were determined by consulting a relevant chemical dictionary or referring to relevant data for final confirmation and screening, and the relative content of phytoncide components of each plant was calculated by combining the peak area normalization method.

Multiple Indicators Comprehensive Evaluation Methods
(1) Data standardization According to the effects of each index on environmental quality, the above four indexes were divided into positive and negative indices.Positive indicators meant that the greater the value, the better the healthcare benefits.In this study, positive indicators were negative air ion concentration, air oxygen content and phytoncide relative content.The negative index meant that the greater the value, the worse the healthcare benefits.The negative index in this study was the human comfort index.In order to bring it into a unified evaluation system for comparison, the range normalization method was adopted for standardization, and the calculation formula was: negative indicator: indicator score = (maximum − current value)/(maximum-minimum) (2) Principal component analysis and weight determination Through standardization, the original value was unified to the dimensionless 0-1, which eliminated the difference of each index's influence on healthcare benefits.Using SPSS20.0, the original data were standardized and analyzed by principal component analysis, with a Kaiser-Meyer-Olkin (KMO) value of 0.504 and a significance of 0.000.According to the requirements of the KMO value range suitable for principal component analysis, the degree of principal component analysis in this study was "average", and the weight can be basically calculated by principal component analysis.As shown in Table 2, the original four environmental indicators could be summarized into two principal components.The cumulative contribution rate of these two principal components reached 85.237%, more than 80%, which can better reflect the information of all indicators.The two principal components were linearly combined as follows: (3) Establishing forest comprehensive healthcare index A comprehensive evaluation index of the forest environmental healthcare benefits was named the "forest comprehensive healthcare index" (FCHI), with the specific calculation formula of FCHI was as follows: (In the formula, NAIC i is the normalized value of negative air ion concentration at the i-th observation site.OC i is the normalized value of air oxygen content at the i-th observation site.HCI i is the normalized value of the human comfort index at the i-th observation site.PRC i w is as the normalized value of the phytoncide relative content at the i-th observation site) (4) Evaluation criteria of the forest comprehensive healthcare index In SPSS20.0, Ward and Euclidean distance methods were used to analyze FCHI values systematically, and a tree diagram of Ward connections was obtained (Figure 4).Combined with literature research and practical needs, the red line was taken as a reference line and divided into five grades according to the value from high to low.The strength of the comprehensive healthcare benefits represented by different numerical intervals was explained; thus, a grading standard of comprehensive evaluation index of the forest environmental healthcare benefits was gained (Table 3).SPSS20.0 was used for variance analysis, multiple comparisons, principal component analysis and systematic cluster analysis; Origin 2021 was used to draw the chart.Due to the impacts of the COVID-19 epidemic, the data of February-April 2020 are missing.Therefore, the spring data are only represented by May 2020.The summer data came from June to August 2020; the autumn data are derived from September to November 2019, and the winter data are derived from December 2019 and January 2020.

Characteristics of Negative Air Ion Concentration in Different Forest Types in Different Seasons
According to the analysis of variance, the negative air ion concentration in different forest types showed a strongly significant difference (p = 0.000 < 0.01) in different seasons (Figure 5).In spring, the order according to the negative air ion concentration mean value (±SD) was: subtropical evergreen broad-leaved forest (7394 ± 1811 ions/cm 3 ) > Phyllostachys edulis forest (2644 ± 1044 ions/cm 3 ) > Cunninghamia lanceolata forest (2461 ± 709

Characteristics of Air Oxygen Content in Different Forest Types in Different Seasons
According to the analysis of variance, the air oxygen content of different stands demonstrated strongly significant differences (p = 0.000 < 0.01) in different seasons (Figure 6).In spring, the order was: coniferous and broad-leaved mixed forest

Characteristics of Human Comfort Index in Different Forest Types in Different Seasons
It is clear that the temperature in the five forest types and forestless control group was highest in summer, followed by spring and autumn, and lowest in winter (Figure 7).In addition, the relative humidity was highest in spring, followed by summer and winter, and lowest in autumn from the vast majority of the forest types, excluding Phyllostachys edulis forest, in which relative humidity in autumn was higher than in winter, as can be seen from Figure 8.The wind speed of the forestless control group was significantly higher than the other five stands throughout the year.Simultaneously, the wind speed of most stands in spring and summer was higher than that in autumn and winter (Figure 9).

Characteristics of Human Comfort Index in Different Forest Types in Different Seasons
It is clear that the temperature in the five forest types and forestless control group was highest in summer, followed by spring and autumn, and lowest in winter (Figure 7).In addition, the relative humidity was highest in spring, followed by summer and winter, and lowest in autumn from the vast majority of the forest types, excluding Phyllostachys edulis forest, in which relative humidity in autumn was higher than in winter, as can be seen from Figure 8.The wind speed of the forestless control group was significantly higher than the other five stands throughout the year.Simultaneously, the wind speed of most stands in spring and summer was higher than that in autumn and winter (Figure 9).
The human comfort index of different forest types and forestless control group were calculated by formula 2 (HCI = 0.6 × |T − 24| + 0.07 × |RH − 70| + 0.5 × |V − 2|) mentioned above.Through variance analysis, the human comfort index of five forest types and forestless control group indicated a highly significant difference (p = 0.000 < 0.01) in different seasons (Figure 10).In spring, the order was: Phyllostachys edulis forest  (6.20).Throughout the whole year, all forest types and the forestless group were at a comfortable level and above.

Evaluation of Comprehensive Healthcare Benefits of Different Forest Types
The values of four environmental healthcare indicators in different seasons were normalized in accordance with the abovementioned formulas (3) and (4).At the same time, in terms of the forest comprehensive healthcare index (FCHI) mentioned in the research methods, the comprehensive healthcare benefits of each stand were evaluated (Table 8).From the table below, according to the FCHI values, the five forest types were all at level I in spring and summer, while the forestless control group at level II.This meant that in spring and summer, the comprehensive healthcare benefits of five forest types were all very strong, which was extremely beneficial to human health.Meanwhile, the comprehensive healthcare benefits of the forest control group were strong as well.In autumn, the subtropical evergreen broad-leaved forest and coniferous and broad-leaved mixed forest were at level I, which had extremely strong healthcare benefits.In addition, Phyllostachys edulis forest, Liquidambar formosana forest and Cunninghamia lanceolata forest were at level II, which had strong healthcare benefits.Moreover, the forestless control group was at the III level; its' healthcare benefits were weak.In winter, Liquidambar formosana forest and Cunninghamia lanceolata forest and coniferous and broad-leaved mixed forest were at level II, which had a strong healthcare benefit.The subtropical evergreen broad-leaved forest and Phyllostachys edulis forest were at level III, and their healthcare benefits were weak.The forestless control group was at level V, having null healthcare benefits.

Evaluation of Comprehensive Healthcare Benefits of Different Forest Types
The values of four environmental healthcare indicators in different seasons were normalized in accordance with the abovementioned formulas (3) and (4).At the same time, in terms of the forest comprehensive healthcare index (FCHI) mentioned in the research methods, the comprehensive healthcare benefits of each stand were evaluated (Table 8).From the table below, according to the FCHI values, the five forest types were all at level I in spring and summer, while the forestless control group at level II.This meant that in spring and summer, the comprehensive healthcare benefits of five forest types were all very strong, which was extremely beneficial to human health.Meanwhile, the comprehensive healthcare benefits of the forest control group were strong as well.In autumn, the subtropical evergreen broad-leaved forest and coniferous and broad-leaved mixed forest were at level I, which had extremely strong healthcare benefits.In addition, Phyllostachys edulis forest, Liquidambar formosana forest and Cunninghamia lanceolata forest were at level II, which had strong healthcare benefits.Moreover, the forestless control group was at the III level; its' healthcare benefits were weak.In winter, Liquidambar formosana forest and Cunninghamia lanceolata forest and coniferous and broad-leaved mixed forest were at level II, which had a strong healthcare benefit.The subtropical evergreen broadleaved forest and Phyllostachys edulis forest were at level III, and their healthcare benefits were weak.The forestless control group was at level V, having null healthcare benefits.In general, the comprehensive healthcare benefits of five forest types and a forestless control group was very strong in spring and summer, followed by autumn and the worst in winter.From the comprehensive data of the whole year, the coniferous and broad-leaved mixed forest is most suitable for the forest recreational and healthcare activities, second, by the subtropical evergreen broad-leaved forest.Moreover, the comprehensive healthcare index of the coniferous and broad-leaved mixed forest was at the "I" very strong level in spring, summer and autumn, except in winter at the "II" strong level.While subtropical evergreen broad-leaved forest at the "I" very strong level in the first three quarters, "III" weak level in winter.Our results demonstrated large seasonal differences in NAIC of five forest types and forestless control group in the Shimen National Forest Park.The seasonal variation of NAIC in Phyllostachys edulis forest, subtropical evergreen broad-leaved forest and forestless control group were summer > spring > autumn > winter, while Liquidambar formosana forest, Cunninghamia lanceolata forest and coniferous and broad-leaved mixed forest were summer > spring > winter > autumn.Generally speaking, NAIC in summer and spring were higher than that in autumn and winter.This finding contrasts with some previous studies in temperate regions that reported the higher NAIC in the summer and autumn, while the lower in the spring and winter [27].This difference can be explained by the source of negative air ions, which naturally occurred from ultraviolet radiation, photosynthesis, etc. [28].The stronger solar radiation and photosynthesis, the higher NAIC [29].In addition, According to the generation mechanism of negative air ions, the main forms of negative air ions were O 2 ) n, etc. and negative air ions can be formed only with certain water content in the environment [30,31].The relationship between NAIC and humidity has long been a focus of research; different scholars [32][33][34][35][36] have studied this dynamic problem and have reached a more general conclusion, namely, that the relationship between NAIC and humidity was positive, more important, the most significant meteorological factor affecting NAIC was the air relative humidity.When the humidity is high, the amount of OH − •(H 2 O)n increases, and with humidity increased, NAIC rises continuously [29].In summer, strong solar radiation, high relative humidity, vigorous plant physiological activities and strong photosynthesis promote the production of negative air ions.Meanwhile, Shimen National Forest Park, located in the south subtropical zone, warmed up earlier in spring when trees entered the growing season, with lush foliage and enhanced photosynthesis.Therefore, its observation results were higher than those in autumn and winter when trees were close to withering, with heaps of fallen leaves and reduced photosynthesis.As can be seen from Figure 6, the relative humidity in spring was highest, followed by summer, which was 15-30% higher than that in autumn and winter.Thus, the NAIC in summer and spring was higher than that in autumn and winter.

Stands Difference
The different forest types were widely separated in terms of NAIC.It was clear that the NAIC of five forest types was remarkably higher than that in the forestless control group.This finding was consistent with most academics [29,37,38].The NAIC in the subtropical evergreen broad-leaved forest was approximately 2-4 times that of the other four stands, which may be related to natural forest community structure, forest area, canopy density and forest community leaf area index, etc.Previous studies have found that the higher canopy density [39], the larger community leaf area index [40], and the more complicated forest community structure [41], the higher NAIC.This was also consistent with the research conclusions of some researchers [41][42][43]: the NAIC of the natural forest was higher than that of the artificial forest.The NAIC in the mature forest was the highest, followed by middle-aged forest and the lowest in young forest.Therefore, as the subtropical evergreen broad-leaved forest is a natural and mature forest, there was no doubt that NAIC produced by it was much higher than that in the other four artificial stands.
In the artificial forest, the NAIC of Phyllostachys edulis forest was significantly higher than that of the other three stands, which was consistent with the research conclusion of Qin Jun et al. [44].This perhaps accounts for the larger leaf area index and higher biomass and canopy density of Phyllostachys edulis forest.Plants can produce higher NAIC through photosynthesis and gas exchange [45].Phyllostachys edulis forest with lush branches and straight trunks were beneficial to photosynthesis and gas exchange.However, there was very little distinction between the NAIC of Cunninghamia lanceolata forest and coniferous and broad-leaved mixed forest (p > 0.05).In spring and autumn, the NAIC of Cunninghamia lanceolata forest and coniferous and broad-leaved mixed forest were higher than Liquidambar formosana forest (p < 0.05), while in summer and winter, the three forest types without differences.This conclusion is in contrast with some scholars [37,42,46] who have advocated the NAIC of the broad-leaved forest is higher than coniferous and broad-leaved mixed forest and coniferous forest.A possible cause might be Liquidambar formosana forest in the young forest, with lower canopy density, which was not conducive to producing negative air ions.

Characteristics of Air Oxygen Content Differences in Different Forest Types
Overall, the variation of air oxygen content in different forest types showed consistency throughout the different seasons, namely, summer > spring > autumn > winter.This finding is in accord with the views of Li Chunhua [47].The air oxygen content has a very significant positive correlation with air temperature and light intensity [48], a significant negative correlation with air humidity [49].The correlation coefficient between air oxygen content and meteorological factors was air temperature > air humidity > light intensity [48].Seasonal variation of air oxygen content in the forest community was not only related to meteorological factors but also related to its growth cycle.In spring and summer, when trees are in the growing season with high light intensity, sufficient sunlight, larger canopy density and vegetation coverage, vegetation has strong photosynthesis and releases more oxygen into the atmospheric environment.This is particularly the case in spring, with the highest humidity all around the year, in which air oxygen content is lower than that in summer.Although there was no deciduous period for vegetation in Shimen National Forest Park, and the forest community still carries out photosynthesis, though the physiological activity of plants in autumn and winter is correspondingly weakened, and the light time and temperature are reduced gradually, which retards the activity of photosynthetic enzymes.These factors all lead to the seasonal differences in oxygen content in the forest communities.Meanwhile, relevant studies had shown that the oxygen content in the air generally accounts for 20.9% of the air volume [22,50].The air oxygen content of different forest types in spring, summer and autumn were higher than 20.9%, except winter.This also reflected that Shimen National Forest Park was just like an enormous forest oxygen bar, whose forest environment was conducive to human health.
As can be seen from Figure 4, the air oxygen content in different forest types was as follows: coniferous and broad-leaved mixed forest > Cunninghamia lanceolata forest > Liquidambar formosana forest > Phyllostachys edulis forest > subtropical evergreen broadleaved forest > forestless control group, which coniferous and the broad-leaved mixed forest was the highest all year round.Consistent with the relationship of the air oxygen content with tree age, the oxygen content was positively correlated with the forest ages in the middle and young forest stage but negatively correlated after the trees mature.Therefore, the oxygen concentration in different forest types of the environment increased with the increase of the forest's age and gradually decreased with the age of trees after they matured [25].Hence, subtropical evergreen broad-leaved forest as the mature forest, which ability to release oxygen was worse than the other four stands.

Characteristics of Human Comfort Index Differences in Different Forest Types
In general, compared to the forestless control group, five forest types done well in the human comfort index.The forestless control group lacked the shelter of trees, the temperature and humidity changed greatly, the wind speed was higher; indisputably, its environment was less comfortable in contrast with the forest surroundings.People reached a consensus in the conclusion that the forest environment plays a better role in regulating microclimate than forestless [23,51,52].The human comfort index of the five stands varied significantly in different seasons.In spring and summer, the Phyllostachys edulis forest and subtropical evergreen broad-leaved forest were more comfortable than Liquidambar formosana forest, Cunninghamia lanceolata forest and coniferous and broadleaved mixed forest.While in autumn and winter, it was in reverse.Among them, the human comfort index of Phyllostachys edulis forest and subtropical evergreen broad-leaved forest were most comfortable in summer, while Liquidambar formosana forest, Cunninghamia lanceolata forest and coniferous and broad-leaved mixed forest are done best in autumn.Most researchers manifested that different microclimate elements had various influences on the human comfort index, among which air temperature had the greatest impact on the human comfort index [24,53,54].The attitude of Phyllostachys edulis forest and subtropical evergreen broad-leaved forest was approximately 800 m, which cold in winter and cool in summer.It is mainly affected by its temperature, which forest circumstances were more comfortable in warm seasons and less comfortable in chill seasons.Liquidambar formosana forest, Cunninghamia lanceolata forest and coniferous and broad-leaved mixed forest were located in about 290-450 m.In summer, three forest types with high-temperature, which was less comfortable than that in autumn.

Characteristics of Phytoncide Relative Content Differences in Different Forest Types
Biogenic volatile organic compounds vary widely between plant species and vegetation types [55].This conclusion was also found in this paper.There were great differences in the phytoncide relative content and species volatilized by five stands in different seasons, from the mean, summer > spring > autumn > winter.The temperature had a crucial effect on the release of phytoncide in terms of environmental factors, and there was a positive correlation between phytoncide volatilization and temperature [56][57][58][59][60].In general, phytoncide relative content gradually increased when the temperature rose.The main reason was that the activity of synthetase was affected by temperature.This directly led to the obvious seasonal variation in phytoncide emissions.In addition, Light had an effect on plant photosynthesis rate, transpiration rate and stomatal conductance [61].The light intensity had a sharp influence on the release of phytoncide, and its release rose with the increase of light intensity [60,62,63].Hence, summer with higher temperature and ample sunlight, in which five forest stands emitted more phytoncide, followed by spring.Furthermore, there were great differences in phytoncide substances emitted from forests in different seasons.This confirmed the views of Xu Jiehua [64]: different enzymes have different responses to temperature, and the components of phytoncide volatilized from forests are various under different temperature conditions.Among the five stands, the species and relative content of phytoncide in the subtropical evergreen broad-leaved forest were apparently higher than those in the other stands, which may be related to the mature and natural secondary forest community, with rich vegetation types, complete community structure, high forest coverage and canopy density, which can produce more species and relative content of phytoncide.The majority of scholars [65][66][67] commonly asserted that tree species play a decisive role in the phytoncide components and the emission rate.Tree age and plant development stage also had a certain impact on phytoncide components and their content.

Conclusions
The forest comprehensive healthcare index (FCHI) was established by using the principal component analysis method, including negative air ion concentration, air oxygen content, human comfort index and phytoncide relative content.Moreover, the evaluation criteria were divided by using the method of system cluster analysis.It provides a theoretical basis for the comprehensive evaluation of the forest healthcare benefits, which is conducive to in-depth research of the forest health environment in the future.To sum up, the FCHI of five stands, which ranged evidently in different seasons, preceded a distinctly forestless control group in the whole year.The forest environment has a favorable influence on the human body and mind, so it is suggested that citizens go to the forest environment persistently for healthcare and physical and mental washing.Moreover, five forest types were all suitable to conduct forest therapy in spring and summer, followed by autumn, and the FCHI in winter were poorer.From the comprehensive data of all years, the coniferous and broad-leaved mixed forest is most suitable for the forest recreational and healthcare activities, followed by the subtropical evergreen broad-leaved forest.Therefore, when planning and constructing the forest therapy base in Shimen National Forest Park in the future, the coniferous and broad-leaved mixed forest should be allocated more in the stand transformation to improve forest healthcare benefits.We should focus on protecting and developing the landscape resources of the subtropical evergreen broad-leaved forest and make rational use of their health activity space.
In this paper, the comprehensive healthcare benefits of different stands in Shimen National Forest Park were analyzed and discussed in depth.Future research work can be carried out from the following aspects: (1) To fully explore the healthcare benefits of the forest, besides the environmental indicators monitored in this paper, other environmental indicators should also be considered, such as air particulate matter and air bacteria content.
(2) The healthcare benefits of the forest ultimately affect the human body.Therefore, we should link forests with human health and study the physical and mental healthcare benefits of different forest communities.(3) In this study, all data were monitored in the forest, but the forest edge is the primary area where people carry out forest recreation activities, so it is necessary to pay attention to the experimental gradient inquiry of the forest edge healthcare benefits subsequently.

Figure 1 .
Figure 1.Location of the Shimen National Forest Park.

Figure 1 .
Figure 1.Location of the Shimen National Forest Park.

Forests
wind speed was 2 m.The relative content of phytoncide was tested by the open, dynamic headspace adsorption collection method.The phytoncide collection devices provided by South China Botanical Garden, Chinese Academy of Sciences, which were hung at 1.5 m in the forest at the average breathing height of the human body, needed to work continuously for at least 8 h from 9:00 to 17:00.The experimental instruments mentioned above are shown in Figure3.

2. 3 . 4 .
Data Processing and Analysis SPSS20.0 was used for variance analysis, multiple comparisons, principal component analysis and systematic cluster analysis; Origin 2021 was used to draw the chart.Due to the impacts of the COVID-19 epidemic, the data of February-April 2020 are missing.Therefore, the spring data are only represented by May 2020.The summer data came from June to August 2020; the autumn data are derived from September to November 2019, and the winter data are derived from December 2019 and January 2020.

Figure 4 .
Figure 4. Results of system cluster analysis-tree diagram.

Figure 4 .
Figure 4. Results of system cluster analysis-tree diagram.

Figure 5 .
Figure 5.Comparison of negative air ion concentration in different forest types.A-Phyllostachys edulis forest; B-Subtropical evergreen broad-leaved forest; C-Forestless control group; D-Liquidambar formosana forest; E-Cunninghamia lanceolata forest; F-Coniferous and broad-leaved mixed forest.The different lowercase letters indicate significant difference among the stand types in the same time (P＜0.05).

Figure 5 .
Figure 5.Comparison of negative air ion concentration in different forest types.A-Phyllostachys edulis forest; B-Subtropical evergreen broad-leaved forest; C-Forestless control group; D-Liquidambar formosana forest; E-Cunninghamia lanceolata forest; F-Coniferous and broad-leaved mixed forest.The different lowercase letters indicate significant difference among the stand types in the same time (p < 0.05).

Figure 6 .
Figure 6.Comparison of air oxygen content in different forest types.A-Phyllostachys edulis forest; B-Subtropical evergreen broad-leaved forest; C-Forestless control group; D-Liquidambar formosana forest; E-Cunninghamia lanceolata forest; F-Coniferous and broad-leaved mixed forest.

Figure 6 .
Figure 6.Comparison of air oxygen content in different forest types.A-Phyllostachys edulis forest; B-Subtropical evergreen broad-leaved forest; C-Forestless control group; D-Liquidambar formosana forest; E-Cunninghamia lanceolata forest; F-Coniferous and broad-leaved mixed forest.The different lowercase letters indicate significant difference among the stand types in the same time (p < 0.05).

Figure 7 .
Figure 7.The temperature changes of different forest types in different seasons.A-Phyllostachys edulis forest; B-Subtropical evergreen broad-leaved forest; C-Forestless control group; D-Liquidambar formosana forest; E-Cunninghamia lanceolata forest; F-Coniferous and broad-leaved mixed forest.

Figure 7 .
Figure 7.The temperature changes of different forest types in different seasons.A-Phyllostachys edulis forest; B-Subtropical evergreen broad-leaved forest; C-Forestless control group; D-Liquidambar formosana forest; E-Cunninghamia lanceolata forest; F-Coniferous and broad-leaved mixed forest.The different lowercase letters indicate significant difference among the stand types in the same time (p < 0.05).

Figure 8 .
Figure 8.The relative humidity changes of different forest types in different seasons.A-Phyllostachys edulis forest; B-Subtropical evergreen broad-leaved forest; C-Forestless control group; D-Liquidambar formosana forest; E-Cunninghamia lanceolata forest; F-Coniferous and broad-leaved mixed forest.The different lowercase letters indicate significant difference among the stand types in the same time (p < 0.05).

Figure 8 .
Figure 8.The relative humidity changes of different forest types in different seasons.A-Phyllostachys edulis forest; B-Subtropical evergreen broad-leaved forest; C-Forestless control group; D-Liquidambar formosana forest; E-Cunninghamia lanceolata forest; F-Coniferous and broad-leaved mixed forest.

Figure 9 .
Figure 9.The wind speed changes of different forest types in different seasons.A-Phyllostachys edulis forest; B-Subtropical evergreen broad-leaved forest; C-Forestless control group; D-Liquidambar formosana forest; E-Cunninghamia lanceolata forest; F-Coniferous and broad-leaved mixed forest.The different lowercase letters indicate significant difference among the stand types in the same time (p < 0.05).

Figure 9 .
Figure 9.The wind speed changes of different forest types in different seasons.A-Phyllostachys edulis forest; B-Subtropical evergreen broad-leaved forest; C-Forestless control group; D-Liquidambar formosana forest; E-Cunninghamia lanceolata forest; F-Coniferous and broad-leaved mixed forest.

Figure 10 .
Figure 10.Comparison of the human comfort index in different forest types.A-Phyllostachys edulis forest; B-Subtropical evergreen broad-leaved forest; C-Forestless control group; D-Liquidambar formosana forest; E-Cunninghamia lanceolata forest; F-Coniferous and broad-leaved mixed forest.

Figure 10 .
Figure 10.Comparison of the human comfort index in different forest types.A-Phyllostachys edulis forest; B-Subtropical evergreen broad-leaved forest; C-Forestless control group; D-Liquidambar formosana forest; E-Cunninghamia lanceolata forest; F-Coniferous and broad-leaved mixed forest.The different lowercase letters indicate significant difference among the stand types in the same time (p < 0.05).

Figure 11 .
Figure 11.Comparison of phytoncide species and relative content in different forest types.A-Phyllostachys edulis forest; B-Subtropical evergreen broad-leaved forest; C-Forestless control group; D-Liquidambar formosana forest; E-Cunninghamia lanceolata forest; F-Coniferous and broad-leaved mixed forest.

Figure 11 .
Figure 11.Comparison of phytoncide species and relative content in different forest types.A-Phyllostachys edulis forest; B-Subtropical evergreen broad-leaved forest; C-Forestless control group; D-Liquidambar formosana forest; E-Cunninghamia lanceolata forest; F-Coniferous and broad-leaved mixed forest.

Table 1 .
Overview of the sample plots.

Table 1 .
Overview of the sample plots.

Table 2 .
Principal component score coefficient matrix and index weight of all indicators.

Table 3 .
Criteria for the forest comprehensive healthcare index (FCHI) grades.

Table 3 .
Criteria for the forest comprehensive healthcare index (FCHI) grades.

Table 8 .
Normalized value of indicators and comprehensive healthcare benefits evaluation.