Relationship between Light Use E ﬃ ciency and Photochemical Reﬂectance Index Corrected Using a BRDF Model at a Subtropical Mixed Forest

: Light use e ﬃ ciency (LUE) is a key indicator of vegetation photosynthesis, which provides important insights into how vegetation productivity responds to environmental conditions. The photochemical reﬂectance index (PRI) is based on reﬂectance at 531 and 570 nm, which reﬂects the xanthophyll cycle process of plants under di ﬀ erent radiation conditions, and makes LUE related to plant optical characteristics. In this study, tower-based PRI and eddy covariance (EC) based LUEs were used to explore the ability of PRI to track LUE variations in a subtropical, evergreen mixed forest in South China. The results indicate that there is a stronger relationship between PRI and LUE, corrected by the bidirectional reﬂectance distribution function (BRDF), where R 2 = 0.46 before correction and R 2 = 0.60 after correction. Generally, PRI is able to capture diurnal and seasonal changes in LUE. Simultaneously, this study highlights a signiﬁcant correlation between LUE and PRI, but there is also a large seasonal di ﬀ erence in its correlation. The correlation in winter was signiﬁcantly stronger than summer. The strongest correlation is found in November (R 2 = 0.91) and the weakest is in July (R 2 = 0.34). Photosynthetically active radiation (PAR) had a strong inﬂuence on the LUE-PRI relationship, while vapor pressure deﬁcit (VPD) and air temperature (Ta) had negative inﬂuences on the relationship between LUE and PRI. Terrestrial laser scanning is used to retrieve the vertical structure of forest crown. Our results show that the vegetation canopy structure (i.e., e ﬀ ective leaf area index, LAI e ), extracted from terrestrial laser scanning (TLS) point data in subtropical mixed forests, had a weak inﬂuence on LUE. Our research suggests that environmental factors and vegetation canopy structures should be considered when using PRI to accurately estimate LUE.


Introduction
Gross Primary Production (GPP) represents the amount of organic carbon fixed by terrestrial plants through photosynthesis and is the driving force for multiple ecosystem functions [1][2][3]. An accurate Therefore, investigating patterns of LUE changes across different canopy structures can provide insights into how differences in forest structure affect LUE [11,37].
In this study, we investigated the relationship between PRI and LUE, obtained from measurements of an AMSPEC spectrometer and eddy covariance (EC) system in a subtropical mixed forest in China. The objectives of this study are: (1) to use the BRDF model to correct the PRI observed by a multi-angle spectrometer; (2) to study the effects of environmental and climatic factors on the relationship between PRI and LUE, and analyze the changes in the relationship between the PRI and LUE; (3) to use the TLS data to extract the vegetation canopy structure parameters, and to determine the influence of the vegetation canopy structure on LUE.

Study Area
The study was carried out in Dinghushan National Nature Reserve, Zhaoqing City, Guangdong Province (112 • Figure 1). According to long-term records of an adjacent weather station, the annual average temperature and precipitation were 20.9 • C and 1956 mm, respectively [38,39]. The dry and wet seasons are obvious. Rainfall occurs primarily during the months of April to September, with 80% precipitation in the whole year. The dry season was from October to March. The type of vegetation was subtropical evergreen conifer and broad-leaved mixed forest. The canopy height was 17 m and the main tree species were schima root-bark (Schima superba Gardner and Champion), Chinese chestnut (Castanopsjs chinensis Hance) and Chinese red pine (Pinus massoniana Lamb.). The soil taxonomy is red and yellow soils [40,41].
Remote Sens. 2020, 12, x FOR PEER REVIEW 3 of 18 In this study, we investigated the relationship between PRI and LUE, obtained from measurements of an AMSPEC spectrometer and eddy covariance (EC) system in a subtropical mixed forest in China. The objectives of this study are: (1) to use the BRDF model to correct the PRI observed by a multi-angle spectrometer; (2) to study the effects of environmental and climatic factors on the relationship between PRI and LUE, and analyze the changes in the relationship between the PRI and LUE; (3) to use the TLS data to extract the vegetation canopy structure parameters, and to determine the influence of the vegetation canopy structure on LUE.

Study Area
The study was carried out in Dinghushan National Nature Reserve, Zhaoqing City, Guangdong Province (112°30′39″-112°33′41″E, 23°09′21″-23°11′30″N), a site of the Chinese Ecosystem Research Network (CERN) and the ChinaFLUX network. The flux tower (23°10′12″N, 112°32′4″E, 240 m above sea level), located in the core area of Dinghushan Nature Reserve, was established in 2002 ( Figure 1). According to long-term records of an adjacent weather station, the annual average temperature and precipitation were 20.9 °C and 1956 mm, respectively [38,39]. The dry and wet seasons are obvious. Rainfall occurs primarily during the months of April to September, with 80% precipitation in the whole year. The dry season was from October to March. The type of vegetation was subtropical evergreen conifer and broad-leaved mixed forest. The canopy height was 17 m and the main tree species were schima root-bark (Schima superba Gardner and Champion), Chinese chestnut (Castanopsjs chinensis Hance) and Chinese red pine (Pinus massoniana Lamb.). The soil taxonomy is red and yellow soils [40,41].

AMSPEC II Spectral Data
The vegetation canopy spectral data was obtained from the AMSPEC II spectral observation system which was installed at a flux tower height of 36 m. This was a multi-angle spectral observation system consisting of a UniSpec-DC Spectral Analysis System (PP Systems, Amesbury, MA, USA), upstream fiber, downstream fiber, an automatic tilt rotation, a data transmission system (PTU-D46-17P70T), and a data storage system under the tower. The data storage system continuously archived the multi-angle observations of vegetation canopy spectral reflectivity.
The upstream fiber was used to measure the downward radiation of the sun while the downstream fiber was used to automatically observe the spectral reflectance of the vegetation canopy. The spectral range of the spectrometer was 300-1100 nm with a spectral resolution of 3.3 nm. The range of viewing zenith angle was set as 42° to 62° and the viewing azimuth angle was set from 10° to 160° and 190° to 340° with an angle step of 10°. Data was archived at 15-min intervals over 24

AMSPEC II Spectral Data
The vegetation canopy spectral data was obtained from the AMSPEC II spectral observation system which was installed at a flux tower height of 36 m. This was a multi-angle spectral observation system consisting of a UniSpec-DC Spectral Analysis System (PP Systems, Amesbury, MA, USA), upstream fiber, downstream fiber, an automatic tilt rotation, a data transmission system (PTU-D46-17P70T), and a data storage system under the tower. The data storage system continuously archived the multi-angle observations of vegetation canopy spectral reflectivity.
The upstream fiber was used to measure the downward radiation of the sun while the downstream fiber was used to automatically observe the spectral reflectance of the vegetation canopy. The spectral Remote Sens. 2020, 12, 550 4 of 17 range of the spectrometer was 300-1100 nm with a spectral resolution of 3.3 nm. The range of viewing zenith angle was set as 42 • to 62 • and the viewing azimuth angle was set from 10 • to 160 • and 190 • to 340 • with an angle step of 10 • . Data was archived at 15-min intervals over 24 h each day. In order to analyze the LUE calculated by the EC flux data, the spectral observation data per half-hour was selected from 9:00-17:00 every day during April 2014, to March 2015. First, the spectrometer data was retrieved. Then, the whiteboard correction parameters and the value of the central wavelength were used to control the quality of data. The whiteboard was a standard diffuse reflection reference for measurement. After the outliers were removed, the canopy reflectance at the wavelengths of 531 and 570 nm were extracted from the archived data. The PRI value was calculated by Equation (1): where R 531 and R 570 represent the reflectivity of the vegetation canopy at 531 and 570 nm, respectively.

The BRDF Model
The PRI obtained from the AMSPEC observation was affected by the observation angle and physiological status. Assuming that the physiological state and lighting conditions of the plants remained unchanged every 15 min, and followed the method by Hilker et al. (2008a). The kernel-driven BRDF model has been widely used in ground measurement and satellite remote sensing. The semi-empirical kernel-driven model of vegetation canopy BRDF are generally comprised of a linear combination of three kernels describing isotropic, geometric, and volume scattering effects [42]. Isotropic scattering assumes that the leaves in the canopy were randomly distributed and have the characteristics of a Lambertian surface. Geometric optical scattering was due to the shape of the vegetation canopy. Volume scattering described the effect of scattering caused by the inside of the canopy [42]. The BRDF model could be decomposed into the sum of three weights of isotropic scattering, volume scattering, and geometric optical scattering. Among them, the Ross-Thick kernel was used for volume scattering, and the Li-Sparse kernel was used for geometrical optical scattering. It can be expressed as: where θ S , θ V , and ∅ are the solar zenith, view zenith, and relative azimuth angles, respectively, F i are a priori kernels based on either physical or empirical considerations, and k i are coefficients based on the observed value [11,43]. We used least square estimation to calculate the corresponding k i for each F i .

Terrestrial Laser Scanning Data
In order to extract the canopy structure factors of coniferous and broad-leaved mixed forest around the flux tower of Dinghushan Nature Reserve, a Trimble TX8 terrestrial laser scanner was used to scan the vegetation and to obtain the point cloud data. The scanner utilized a vertical rotating mirror based on a horizontal rotation with a scanning speed of 1 million points/second. The maximum standard measurement range was 120 m, the minimum measurement range was 0.6 m, and the noise of the measurement range was <2 mm. The angle range of the scanning field was 360 • × 317 • .
The forest canopy structure was scanned and a total of six stations were arranged. The distance between every two stations was greater than 10 m ( Figure 2). More than three target balls were placed between every two stations, and the three balls were not in the same row. The target ball was used as the feature point; the point cloud data stitching and registration between the two stations were completed. The point cloud data collected by the TLS scanning was like a hemispherical sphere with angular information. With the spectrometer at the center, the points within a radius of 10 m were selected for extraction of LAI e . First, the TLS data was divided into different parts according to the azimuth angle of 45 • . Second, we visually resolved each part and manually measured the height of the highest point on the ground, and manually removed the highest point and the points below the highest point. Finally, the residual points of the 3D point cloud data were projected onto the XOY plane. Third, the porosity of each part was calculated by the volume element method, and the second part of LAI e was calculated. Among them, based on Beer's law, the extinction coefficient was set to 0.5. In order to explore the effects of vegetation canopy structure on PRI and LUE, the LAI e under different classifications were calculated by the TLS data; the PRI spectral observation data of each record was graded based on the observed azimuth, divided into one every 45 • , and sliding was performed at an average of 10 • ; there were 33 levels.
Remote Sens. 2020, 12, x FOR PEER REVIEW 5 of 18 part of LAIe was calculated. Among them, based on Beer's law, the extinction coefficient was set to 0.5. In order to explore the effects of vegetation canopy structure on PRI and LUE, the LAIe under different classifications were calculated by the TLS data; the PRI spectral observation data of each record was graded based on the observed azimuth, divided into one every 45°, and sliding was performed at an average of 10°; there were 33 levels.

Eddy Covariance Measurement and the LUE Calculation
The flux tower was installed at 27 m, consisting of a three-dimensional ultrasonic anemometer (CSAT3, cAMpbell Scientific Ltd., USA) and fast response infrared CO2/H2O analyzer (Li-7500, Li-CorInc., USA) composition. The CO2 and H2O flux data were automatically archived with a CR500 data logger at a frequency of 10 Hz, and the average flux data was calculated at 30 minutes intervals. In addition, meteorological measurements above the vegetation canopy included canopy incident radiation, air temperature, saturated water vapor pressure, precipitation, soil moisture, etc. The net CO2 exchange (NEE) was calculated from the CO2 storage flux and turbulent flux obtained from flux observations, and the ecosystem respiration (Re) was obtained using an empirical relationship between nighttime requirements at a 50 mm depth and at soil temperature. Re is the ecosystem respiration during the day, obtained by deducting the exponential relationship between NEE at night and soil temperature at a 5 cm depth. GPP and LUE were estimated using the following equations: Photosynthetically active radiation (PAR) refers to the incident radiation at the top of the vegetation canopy, which can usually be obtained from meteorological observations. A fraction of the absorbed photosynthetically active radiation (FPAR) refers to the ratio of PAR absorbed by the green part of the vegetation canopy to total the PAR. It can be obtained from remote sensing data [5,44]. FPAR can be calculated as follows: where and are NDVI values for vegetation coverage of 5% and 98%, respectively. According to the characteristics of the vegetation canopy structure in the observation area, = 0.01, and = 0.5. and represent the minimum and

Eddy Covariance Measurement and the LUE Calculation
The flux tower was installed at 27 m, consisting of a three-dimensional ultrasonic anemometer (CSAT3, Campbell Scientific Ltd., Logan, UT, USA) and fast response infrared CO 2 /H 2 O analyzer (Li-7500, Li-CorInc., Lincoln, NE, USA) composition. The CO 2 and H 2 O flux data were automatically archived with a CR500 data logger at a frequency of 10 Hz, and the average flux data was calculated at 30 min intervals. In addition, meteorological measurements above the vegetation canopy included canopy incident radiation, air temperature, saturated water vapor pressure, precipitation, soil moisture, etc. The net CO 2 exchange (NEE) was calculated from the CO 2 storage flux and turbulent flux obtained from flux observations, and the ecosystem respiration (Re) was obtained using an empirical relationship between nighttime requirements at a 50 mm depth and at soil temperature. Re is the ecosystem respiration during the day, obtained by deducting the exponential relationship between NEE at night and soil temperature at a 5 cm depth. GPP and LUE were estimated using the following equations: Photosynthetically active radiation (PAR) refers to the incident radiation at the top of the vegetation canopy, which can usually be obtained from meteorological observations. A fraction of the absorbed photosynthetically active radiation (FPAR) refers to the ratio of PAR absorbed by the green part of the Remote Sens. 2020, 12, 550 6 of 17 vegetation canopy to total the PAR. It can be obtained from remote sensing data [5,44]. FPAR can be calculated as follows: where NDVI min and NDVI max are NDVI values for vegetation coverage of 5% and 98%, respectively. According to the characteristics of the vegetation canopy structure in the observation area, NDVI min = 0.01, and NDVI max = 0.5. FPAR min and FPAR max represent the minimum and maximum FPAR values, respectively, which are generally assumed to be 0.001 and 0.95, respectively. The FPAR was obtained by the normalized difference vegetation index (NDVI) estimation method [44], and the vegetation reflectivity in each band was obtained by an automatic multi-angle spectroscopy system. The NDVI can be calculated as follows: where R 850 and R 680 represent the reflectance of plant leaves at 850 and 680 nm, respectively, observed by an automated multi-angle spectroscopy system.  Figure 3a shows that the PRI did not change significantly at different observation angles, which was almost close to each other. Figure 3b shows the variation of PRI from different observation angles. The PRI values in the backward scattering direction were generally higher than the forward scattering direction. When the VAA is about 135 • , it is close to the solar incident azimuth angle (called backward scattering direction), and the value of PRI is greater than the VAA of 315 • (forward scattering direction at this time).

Variation of PRI at Different Observation Angles and the Result of BRDF Model
Remote Sens. 2020, 12, x FOR PEER REVIEW 6 of 18 The FPAR was obtained by the normalized difference vegetation index (NDVI) estimation method [44], and the vegetation reflectivity in each band was obtained by an automatic multi-angle spectroscopy system. The NDVI can be calculated as follows: where and represent the reflectance of plant leaves at 850 and 680 nm, respectively, observed by an automated multi-angle spectroscopy system.

Variation of PRI at Different Observation Angles and the Result of BRDF Model
PRI changed with observation angles; it was correlated with vegetation canopy structures and different environmental factors, produced under different observation angles. Figure 3a shows spectral observations at 10:30-10:45 on 13 July 2014, when the solar azimuth angle was 111° and the solar zenith angle was about 21°. Figure 3b shows spectral observations at 10:30-10:45 on 13 November 2014, when the solar azimuth angle was about 160° and the solar zenith angle was about 51°. Figure 3a shows that the PRI did not change significantly at different observation angles, which was almost close to each other. Figure 3b shows the variation of PRI from different observation angles. The PRI values in the backward scattering direction were generally higher than the forward scattering direction. When the VAA is about 135°, it is close to the solar incident azimuth angle (called backward scattering direction), and the value of PRI is greater than the VAA of 315° (forward scattering direction at this time).

Temporal Variation of the Relationship between PRI and LUE
Generally, PRI and LUE are positively correlated, and there is a negative correlation on rare days (R < 0) ( Figure 6). During 288 days, PRI and LUE were significantly associated for 84 days (R > 0.75). Overall, from 6 April 2014, to 1st August, 2014, the correlations of PRI and LUE start to decrease and then increase, then the correlations slightly decrease around 1 August 2014, and increase and stabilize

Temporal Variation of the Relationship between PRI and LUE
Generally, PRI and LUE are positively correlated, and there is a negative correlation on rare days (R < 0) ( Figure 6). During 288 days, PRI and LUE were significantly associated for 84 days (R > 0.75). During each month, there are some differences in the relationship between PRI and LUE ( Table  1). The linear correlations between PRI and LUE are the smallest from June to August. The smallest correlation occurred in July 2014 (R 2 = 0.   Table 2 shows the response of LUE and PRI to PAR, VPD, and Ta in the wet season (from April 2014 to September 2014) and dry season (from October 2014 to March 2015). In the dry season, LUE was affected differently by PAR, VPD, and Ta, among which PAR had the greatest impact on LUE (R 2 = 0.85), followed by VPD (R 2 = 0.26), and the effect of Ta on LUE can be ignored. In the wet season, PAR had the greatest effect on LUE (R 2 = 0.58). VPD and Ta had no significant difference in its affects on LUE. The correlation of LUE and PAR was higher in the dry season than in the wet season. In the dry season, similar to the response of LUE, PRI had the greatest response to PAR (R 2 = 0.67), and VPD had the second greatest effect on PRI (R 2 = 0.24), while temperature had no effect on PRI (Table 2). In the wet season, the responses of PRI to PAR and Ta were more obvious (R 2 = 0.59 and R 2 = 0.55), while During each month, there are some differences in the relationship between PRI and LUE ( Table 1) Table 2 shows the response of LUE and PRI to PAR, VPD, and Ta in the wet season (from April 2014 to September 2014) and dry season (from October 2014 to March 2015). In the dry season, LUE was affected differently by PAR, VPD, and Ta, among which PAR had the greatest impact on LUE (R 2 = 0.85), followed by VPD (R 2 = 0.26), and the effect of Ta on LUE can be ignored. In the wet season, PAR had the greatest effect on LUE (R 2 = 0.58). VPD and Ta had no significant difference in its affects on LUE. The correlation of LUE and PAR was higher in the dry season than in the wet season. In the dry season, similar to the response of LUE, PRI had the greatest response to PAR (R 2 = 0.67), and VPD had the second greatest effect on PRI (R 2 = 0.24), while temperature had no effect on PRI (Table 2). In the wet season, the responses of PRI to PAR and Ta were more obvious (R 2 = 0.59 and R 2 = 0.55), while VPD had a weaker effect on PRI (R 2 = 0.29). Overall, LUE and PRI had similar responses to environmental factors.  Figure 7a and b illustrate the effects of air temperature on LUE and PRI, observed by an AMSPEC spectrometer. Figure 7c reveals the change in the correlation between the PRI and the LUE as a function of temperature. During the study, the air temperature ranged from 5 • C to 35 • C. When the air temperature is less than 20 • C, there is no correlation between temperature and LUE and PRI (Figure 7a,b). When the temperature is higher than 20 • C, the LUE decreases gradually with the increase of temperature. The correlation between PRI and LUE increases and then decreases with rising air temperature. When the temperature is between 10 • C to 15 • C, the correlation between PRI and LUE is at its largest (R 2 = 0.78). The changes in LUE and PRI, as affected by VPD and the changes in the LUE and PRI as a function of VPD, are shown in Figure 7d-f. LUE and PRI decreased with VPD, and the correlation between PRI and LUE decreased with increasing VPD. When VPD < 0.5, the best correlation between PRI and LUE (R 2 = 0.22) occurred. When VPD ranges from 1.5 and 2, there is no correlation between LUE and PRI.

The Effects of Environmental Factors on the Correlation between PRI and LUE
Remote Sens. 2020, 12, x FOR PEER REVIEW 10 of 18 VPD had a weaker effect on PRI (R 2 = 0.29). Overall, LUE and PRI had similar responses to environmental factors.   Figure 7a and b illustrate the effects of air temperature on LUE and PRI, observed by an AMSPEC spectrometer. Figure 7c reveals the change in the correlation between the PRI and the LUE as a function of temperature. During the study, the air temperature ranged from 5 °C to 35 °C. When the air temperature is less than 20 °C, there is no correlation between temperature and LUE and PRI (Figure 7a and b). When the temperature is higher than 20 °C, the LUE decreases gradually with the increase of temperature. The correlation between PRI and LUE increases and then decreases with rising air temperature. When the temperature is between 10 °C to 15 °C, the correlation between PRI and LUE is at its largest (R 2 = 0.78). The changes in LUE and PRI, as affected by VPD and the changes in the LUE and PRI as a function of VPD, are shown in Figure 7d, e, and f. LUE and PRI decreased with VPD, and the correlation between PRI and LUE decreased with increasing VPD. When VPD < 0.5, the best correlation between PRI and LUE (R 2 = 0.22) occurred. When VPD ranges from 1.5 and 2, there is no correlation between LUE and PRI.  Figure 8 shows the overall pattern of LUE simultaneously responding to air temperature, VPD, and PAR. Different temperature conditions may correspond to different radiation conditions, and different air humidity. Additionally, it was also established, in turn, that LUE showed a decreasing trend when radiation increases. The LUE also decreased slightly with the increase of VPD. However, LUE was not affected by air temperature.

Effects of Canopy Structural Complexity
The results of forest canopy structure parameters are shown in Table 3. In the azimuth angle of 0-360°, the forest canopy structure showed an obvious difference in spatial space. Table 3 lists the spatial distribution of LAIe and LUE. The space usage data of LUE is predicted by AMSPEC-PRI. The  Figure 8 shows the overall pattern of LUE simultaneously responding to air temperature, VPD, and PAR. Different temperature conditions may correspond to different radiation conditions, and different air humidity. Additionally, it was also established, in turn, that LUE showed a decreasing trend when radiation increases. The LUE also decreased slightly with the increase of VPD. However, LUE was not affected by air temperature.  Figure 8 shows the overall pattern of LUE simultaneously responding to air temperature, VPD, and PAR. Different temperature conditions may correspond to different radiation conditions, and different air humidity. Additionally, it was also established, in turn, that LUE showed a decreasing trend when radiation increases. The LUE also decreased slightly with the increase of VPD. However, LUE was not affected by air temperature.

Effects of Canopy Structural Complexity
The results of forest canopy structure parameters are shown in Table 3. In the azimuth angle of 0-360°, the forest canopy structure showed an obvious difference in spatial space. Table 3 lists the spatial distribution of LAIe and LUE. The space usage data of LUE is predicted by AMSPEC-PRI. The

Effects of Canopy Structural Complexity
The results of forest canopy structure parameters are shown in Table 3. In the azimuth angle of 0-360 • , the forest canopy structure showed an obvious difference in spatial space. Table 3 lists the spatial distribution of LAI e and LUE. The space usage data of LUE is predicted by AMSPEC-PRI. The spatial distribution is slightly different. Figure 9 illustrates the relationship between LAI e and LUE. Obviously, as the LAI e increases, the LUE decreases slightly and then increases (R 2 = 0.13, p = 0.07).

BRDF Model Correction
The heterogeneity of the forest canopy structure resulted in directional spectroscopic observations of the vegetation canopy reflectivity. Compared to the traditional vertical observation method, using multi-angle observations makes it easier to capture the spatial structure characteristics of the vegetation canopy. Using the BRDF model correction, the spatial structure characteristics of vegetation can be more accurately observed through spectral observations [45], can eliminate the R² = 0.13 p=0.07

BRDF Model Correction
The heterogeneity of the forest canopy structure resulted in directional spectroscopic observations of the vegetation canopy reflectivity. Compared to the traditional vertical observation method, using multi-angle observations makes it easier to capture the spatial structure characteristics of the vegetation canopy. Using the BRDF model correction, the spatial structure characteristics of vegetation can be more accurately observed through spectral observations [45], can eliminate the uncertainty brought by the sun and observation angles, and increase the comparability of different observations. Both the observation angle and the sun angle were changed before the model was corrected. The changing angle will cause the variance of the observation data to increase. After the model is corrected, the angle is fixed, the uncertainty caused by the angle factor will no longer exist, and the estimation accuracy of the LUE and PRI will be improved [11]. This result is also consistent with the results obtained by Jia et al. [11]. Compared with the results obtained by Jia et al., the R 2 of exponential adjustment in our study is lower. This may be due to the different types of forests. However, spectroscopic observations are also affected by environmental factors such as radiation and temperature, which are also factors that affect the accuracy of PRI in estimating LUE [10,24,[45][46][47].

Response of PRI and LUE to Environmental Factors
The growth environment (light intensity, temperature, and moisture) in which the plant is located affects the physiological state of the plant [38]. For subtropical evergreen forests, PRI and LUE exhibit very complex changes due to the changing weather environment. During the whole study period, PRI showed a significant correlation with LUE ( Figure 4). Significant positive correlations between PRI and LUE were found on some days ( Figure 5), months (Table 1), as well as throughout the study period ( Figure 5). These results are consistent with previous studies conducted on canopy scales and different temporal scales [10,11,25,43,47]. PRI can be used to track LUE in subtropical evergreen coniferous and broad-leaved mixed forest. PRI can be used to track LUE in subtropical evergreen coniferous and broad-leaved mixed forests. A better correlation between PRI and LUE requires seasonal and monthly scale analysis ( Figure 6 and Table 1). Obviously, the correlation between PRI and LUE in the dry season is higher than that in the wet season. From April to March of the next year, the correlation between PRI and LUE shows a trend of decreasing first and then increasing. The relationship between the PRI and LUE is mainly affected by environmental factors. Among them, radiation is the main driver of plant photosynthesis, and vegetation has higher photosynthesis under high light and optimal temperature conditions ( Table 2). This is consistent with previous results [11]. However, under low radiation conditions, high LUE in plants was observed, which represented the best physiological growth state of vegetation light energy utilization. In contrast, LUE was lower under higher radiation conditions. At the same time, when the vegetation is in the optimal temperature condition, the radiation intensity increases, while the LUE decreases. Therefore, during the study period, the half-hour scale observation of radiation intensity plays a dominant role, in comparison to the temperature. This conclusion is consistent with the results reported by Jia et al. (2008). While illumination is the determinant of the correlation between PRI and LUE, in the dry season, PRI and LUE are also sensitive to VPD. For example, as VPD increases, LUE also increases, indicating that LUE is limited to a certain extent to a high VPD environment. This is consistent with the results of Zhang Qian et al. (2015) [10]. The difference is that PRI and LUE are not affected by temperature during the dry season, indicating that the temperature does not interfere with the characterization of LUE by PRI. Therefore, in the subtropical evergreen coniferous and broad-leaved mixed forest, both radiation and VPD will become important factors affecting the relationship between PRI and LUE. The different meteorological conditions corresponding to different months alter the relationship between PRI and LUE. For subtropical forests, finding an empirical equation for the relationship between PRI and LUE for a full growing season is still a problem.
In our study, compared to temperate deciduous forests, the LAI e changes in the subtropical mixed forests in the Dinghushan area are small [40], the phenological changes are not significant, and the changes in pigment content in the leaves of subtropical coniferous and broad-leaved mixed forests are small. This is the reason why there is a weak dynamic change of LUE.

Effects of Vegetation Canopy Structure on LUE
For subtropical mixed forests, canopy structure is not a main factor affecting LUE change. The estimation of PRI is complicated by the light conditions in the tree crown, especially for the closed crown with complex light conditions, which affects the estimation of LUE [48]. Therefore, we attempted to combine optical data with lidar data in order to infer the impact of forest structure on plant photosynthesis. The spatial distribution difference of LAI e does not show a great impact on the change of LUE. In our study, we believe that the uneven terrain factors in the study area have an impact on the accuracy of the LAI e extracted by the TLS. For the topographic relief, the data of TLS cannot obtain the complete point cloud data of the canopy, especially the upper part of the canopy, which affect the LAI e extraction results. In addition, the study area is located in the subtropical region with a humid and hot climate. The cloud layer is thick and has a strong scattering and refraction effect on the solar light [41]. Most of the light reaching the crown is scattered light, which leads to more multiple scatterings and more extensive light conditions [11,47,48]. This result is inconsistent with other results [11,15,49]. The higher the structural complexity of the canopy, the more it will pass through the canopy, resulting in more light interception, resulting in more multiple scatterings and wider radiation conditions [50][51][52]. The reason for the inconsistency with the results found by Jia et al (2018) may be that the direct radiation received by the boreal forest is higher than that of tropical forests, and PAR and illumination conditions are one of the most important factors affecting LUE [11].

Conclusions
Multi-angle tower-based optical spectra data were observed with EC flux data during the period from April 2014 to March 2015 for the purpose of studying the ability of PRI in estimating the LUE of a sub-tropical planted mixed forest stand of Southern China. The following conclusions could be drawn from this study: 1. With the increase of VPD, Ta and PAR, the daily scale PRI, and LUE showed a downward trend, and the response of LUE to changes in these bioclimatic factors were more sensitive than PRI. The correlation of PRI and LUE at half-hour scale was the most significantly affected by PAR, followed by temperature in the wet season. The correlation of PRI and LUE at a half-hour scale was mostly affected by PAR, followed by VPD in the dry season.
2. During the study period, there was a significant linear relationship between LUE and PRI on both the half-hour and daily scale. There was a significant difference in correlation between PRI and LUE throughout the growing season. The correlation was the strongest in November, 2014 (R 2 = 0.91), and the weakest in July, 2014 (R 2 = 0.24).
3. For subtropical forests, the spatial difference of the canopy is small, and has weak effects on LUE.
In addition to these findings, some shortcomings in our research deserve some attention. For example, the measurements of some pigments (e.g., xanthophyll and chlorophyll) can provide more useful information to further analyze seasonal changes in subtropical mixed forests. In our subsequent work, it will be necessary to pay attention to the influence of the vegetation canopy structure on LUE.