Assessing Environmental Control of Sap Flux of Three Tree Species Plantations in Degraded Hilly Lands in South China

: Prerequisite for selection of appropriate tree species in afforestation programs is to understand their water use strategy. Acacia mangium Willd., Schima wallichii Choisy, and Cunninghamia lanceolata (Lamb.) Hook are the three main vegetation restoration pioneer species in southern China, but no comparative research on the water use strategy of these three tree species have been reported. Our objective was to gain a detailed understanding of how photosynthetically active radiation (PAR), vapor pressure deficit (VPD), and soil water content (SWC) at different soil depths control the sap flux density ( J s ) in the dry and wet seasons. We measured the J s of these three tree species by using the thermal dissipation method in low subtropical China. We found that both S. wallichii and C. lanceolata differed clearly in their stomatal behavior from one season to another, while A. mangium did not. The canopy conductance per sapwood area of S. wallichii and C. lanceolata was very sensitive to VPD in the dry season, but not in the wet season. The J s of A. mangium was negatively correlated to SWC in all soil layers and during both seasons, while the other two species were not sensitive to SWC in the deeper layers and only positively correlated to SWC in dry season. Our results demonstrate that the three species have distinct water use strategies and may therefore respond differently to changing climate. measure the xylem sap flux density ( J s ) of these species. Our aim is to investigate how photosynthetically active radiation (PAR), vapor pressure deficit (VPD), soil water content (SWC) at different soil depths, and temperature affect J s of these three tree species in the dry and wet seasons.


Introduction
China, as the second largest plantation country in the world, has implemented large-scale planned afforestation since the 1950s, primarily focused on planting fast-growing and high production timber species during the 1970s and 1980s [1,2]. These programs have altered the coverage of forest in China from 16.0% in the 1980s to 21.66% in 2018 and this increase is expected to continue as China intends to extend afforestation for carbon storage and environmental protection [3,4]. Since the 1980s, an extensive afforestation campaign has been conducted in Guangdong Province, and consequently the forest coverage in this region had risen from 26.2% in 1979 to 58.9% in 2016 [5,6]. Although the plantation of pioneer tree species promotes wood yield and increases commercial profits, it consumes a large amount of water [7][8][9]. However, the use of water by plantations and the selection of appropriate species in many afforestation programs have not previously been considered [10].
Since the 1980s, exotic pioneer species such as Acacia mangium Willd, an evergreen broad-leaved tree with nitrogen-fixing, drought-tolerant, and fast-growing characteristics, have been introduced in southern China [11,12]. They can be effectively used due to the nitrogen fixation of the nodule and rapid growth in the early stage to improve soil conditions at the plantations [13]. Observations show that the growth, photosynthetic rate and water use of A. mangium plantations decreases after approximately 20-22 years after planting [12]. It can be expected that the use of water is essential to support the growth rate of A. mangium [14]. Therefore, understanding the water use strategy of the mature A. mangium has been a concern. Schima wallichii Choisy is a pioneer, evergreen, broad-leaved and fire-resistant tree species used for timber production and fire protection [15]. It has a wide ecological niche, is particularly common in disturbed and secondary forests [15] and is an important native species for ecological restoration in subtropical southern China. At present, however, there is little research on the water usage of S. wallichii in degraded and non-degraded ecosystems. Cunninghamia lanceolata (Lamb.) Hook, as an evergreen conifer, is the main traditional pioneer species that has been used for cultivation and wood production due to its rapid growth and excellent wood quality [16]. Moreover, C. lanceolata plantations also have positive impacts on ecosystem services, including carbon sequestration [17] and decreasing runoff [18]. According to FAO (2010) [1], in the past four decades, C. lanceolata plantations have expanded rapidly and their total planted area has now exceeded 12 million ha in southern China. Despite the positive ecosystem services they provide, there is still a lack of studies on the water use of C. lanceolata [19]. Selecting appropriate species and managing forest effectively can optimize the water use in plantations [20,21] as well as potentially maximize economic benefits. However, no comparative research on the water use strategy of these three common vegetation restoration pioneer tree species in southern China has been published.
In South China, human disturbance seriously affected the native evergreen broad-leaved forest, causing soil erosion and leading to vegetation degradation [22]. Since the mid-1980s, a vegetation restoration campaign has been carried out to replace the barren, hilly grasslands with pioneer species including exotic and native species [23]. A. mangium, S. wallichii, and C. lanceolata are the main vegetation restoration tree species, planted in 1985-1990, and they have developed into secondary subtropical forests after more than 30 years under a monsoon climate. These secondary forests are the product of long-term influences of a monsoon climate. Long-term observations have shown that the total annual precipitation has changed little since 1950 in southern China [24]. In recent decades, however, rising air temperature in southern China has shifted the rainfall patterns toward intensified rainfall (i.e., more floods) but has not raised the soil water content (SWC) in the wet season, while droughts are becoming frequent and more severe leading to declining SWC in the extended dry season [24].
In the current study, we set out to find which of these three tree species would be most suitable for planting in the degraded hilly lands in terms of adaptive water-use strategy, especially under the changing climate. In order to answer this question, we should know how the water use of these species is influenced by the environmental factors. We studied the water use of A. mangium, S. wallichii, and C. lanceolata in the degraded hilly lands of southern China using the thermal dissipation technique [25] to measure the xylem sap flux density (Js) of these species. Our aim is to investigate how photosynthetically active radiation (PAR), vapor pressure deficit (VPD), soil water content (SWC) at different soil depths, and temperature affect Js of these three tree species in the dry and wet seasons.

Site description
The experimental site (112°54′E, 22°41′N) is located in the Heshan National Field Research Station of Forest Ecosystem, Chinese Academy of Sciences, Guangdong Province, China. The average elevation of the station is 80 m. This region is dominated by a subtropical monsoon climate, with an annual average precipitation of 1700 mm, and temperature of 21.7 °C. The annual accumulated temperature above 10 °C is 7600 °C, and the annual solar radiation is 4350 MJ m −2 . There are distinct wet and dry seasons in this region. The wet season lasts approximately from April to September, while the dry season is approximately from October to March [26].
Our experiment was conducted in the three mono-species plantations of A. mangium, S. wallichii, and C. lanceolata (Table 1) from January to December in 2017. These three plantations were planted between 1985 and1990, located close to each other on different hill slopes with similar elevation (80 m), inclination, slope (20%-30%) direction, and soil properties ( Table 2). We randomly selected 15 trees of each tree species for the sap flow experiment (Table 3).  PAR, temperature, relative humidity (above 2 m), and wind speed data (above 10m) were measured (VAISALA MAWS301, Vaisala, Espoo, Finland) and data was provided by the Heshan National Field Research Station of Forest Ecosystem. The leaf area index (LAI) was measured by LAI 2000 (LI-COR, Lincoln, NE, USA) monthly. The linear distances from the station to the A. mangium, S. wallichii, and C. lanceolata sites were 300, 50 and 100 m, respectively. Hourly averages were used for all meteorological data and we calculated VPD (kPa) according to Campbell and Norman, 1998 [27]: where a, b, and c are constant parameters i.e., 0.611 kPa, 17.502 (unitless), and 240.97 °C, respectively, and RH is relative humidity (Units), T is temperature (Units).

Soil Water Content
The soil type at all the sites is lateritic red soil. The soil water content (SWC) was measured at 5, 10, 20, 30, 40, 50, 60, 90, 120, and 150 depths at three locations (up-, mid-, and downslope) at each site using soil moisture sensors (A755 GPRS, Adcon, Klosterneuburg, Austria). An Adcon A755 telemetry unit (Adcon) transmitted the data signals to a data server at 30 min intervals.

Sap Flow
We used home-made Granier's sensors (thermal dissipation probes, TDPs) to measure sap flow of 15 sample trees in each of the three plantations [25]. The TDP sensors consisted of a pair of 20-mmlong and 2 mm in diameter stainless steel probes. Each pair of probes was inserted on the northern side of the stem (1.3 m above ground) approximately 10 cm apart along the axis of the hydro-active xylem. The upper probe was heated by a DC supply of 120 mA, yielding a constant power of 0.2 W, while the lower probe remained unheated. We converted the instantaneous temperature difference between the two probes into a voltage value recorded by a data logger (Delta-T logger DL2e, Delta-T Devices Ltd, Cambridge, UK). The data were measured every 30 s and stored as 10 min averages [14]. Finally, we calculated Js (g m −2 s −1 ) following Granier (1987) [25]: where ΔTm is the temperature difference between the two probes obtained under zero flux conditions and ΔT is the instantaneous temperature difference. The "zero baseline" was determined on nights when VPD was zero or nearly zero for several hours, so that there was no driving force for sap flux.
ΔTm was determined separately for each tree over 7 days to avoid the underestimation of nocturnal sap flow [28]. The voltage value data was converted to Js by applying the Baseliner 4.0 program [29] (https://github.com/Coweeta/Baseliner4_exe).

Canopy Conductance per Sapwood Area
Because of the difficulty in obtaining the leaf area data of measured trees, canopy conductance per sapwood area (g, mmol m −2 s −1 ) was determined by using the following formula: which presumes that Js is equal to canopy transpiration, and Js is not influenced by hydraulic capacitance [30,31].

Data Analysis and Modeling
We analyzed both the daily and hourly relationships between the environmental conditions and Js since the daily analysis fails to catch fast responses to changing environmental conditions, while the hourly analysis in obscured by time lags. We concentrated more on the analysis of the daily data by looking at daily average values (which need not to occur at the same moment in time) to minimize the problems caused by the time lag due to the internal water stores, i.e., capacitance [32], and due to the slowness of the thermal dissipation probe method to respond to changes in flow rate [33,34]. We first used the daily average values of Js, PAR, VPD, and SWC, then after correcting the time lag between VPD, PAR and Js, we used the hourly average values minus the daily average values of Js, PAR, and VPD to analyze the hourly dynamics of Js within a day. We analyzed the relationship between the dependent variable Js and the explanatory variables PAR, VPD, and SWC in the daily and hourly (only PAR and VPD) tree-level data: The parameters β1, β2, and β3 describe the sensitivity of Js to the natural logarithm (better fit than other shapes) of PAR, VPD, and SWC, respectively (note that this is equivalent to fitting a power function Js = α (PAR) β1 (VPD) β2 (SWC) β3 ). All the statistical analyses were produced by R software (R Core Team; R version 3.5.3; RStudio version 1.1.463). The analysis for hourly data was made separately (1) for each species during the wet and dry seasons, (2) for small and large trees of each species, and (3) for each species with SWC measured at different soil depths. The size class in (2) was determined based on DBH separately for each species (note that in (1) and (2) the data were divided into two approximately similar size groups). The random term tree identity was added into the models to consider the dependence of observations for the same tree, and all values were transformed for the statistical analysis with natural logarithm. Possible differences in the relationship between Js and the explanatory variables VPD, PAR, and SWC between species, between wet and dry seasons, between small and large trees, and between models with SWC measured at different soil depths were compared with 95% confidence intervals (Tables 4-6); for example, if the confidence interval of the estimate of the effect of VPD on Js in the wet season overlapped with the estimate of the effect of VPD on Js in the dry season in A. mangium, no significant difference was reported. Temperature was left out of the model because temperature is highly correlated with another variable in the model, VPD. However, we compared the residuals of the model with temperature to test whether temperature had a direct (i.e. not through VPD) effect on sap flow. The effect of wind speed on Js was also analyzed, but it was not significant and was thus dropped from the final analyses.

Micrometeorology
The total amount of precipitation recorded was 1673 mm from January 2017 to December 2017, which was in the range of the long-term annual mean in the Heshan area. Rainfall in the wet season (1289 mm) accounted for 77% of the total annual precipitation (Figure 1a). The mean daily PAR was 223 (±9) µmol m −2 s −1 and 326 (±13) µmol m −2 s −1 during the dry and wet periods, respectively ( Figure  1c). The mean daily VPD was 0.55 (±0.03) kPa and 0.59 (±0.03) kPa in the dry and wet season, respectively (Figure 1d). The average daily temperature ranged from 8.7 to 31.8 °C during the whole measurement period (Figure 1f). The mean daily PAR, VPD, and temperature were all higher in the wet season than in the dry season. The difference in mean SWC of all layers between the wet and dry seasons was greater at the A. mangium and S. wallichii sites than at the C. lanceolata site (Figure 2). At the A. mangium site, average SWC of all layers was lower in the dry season (19%-28%) compared with the wet season (28%-33%) (Figure 2). At the S. wallichii site, average SWC of all layers ranged from 20% to 27% in the dry season and from 26% to 32% in the wet season, which can be compared with the average SWC of all layers of 20% to 26% in the dry period and 25% to 30% in the wet period at the C. lanceolata site (Figure 2). The leaf area index (LAI) of all the three sites decreased in the dry season compared with the wet season, and the decrease at the A. mangium site was larger (25%) than at the S. wallichii (18%) and C. lanceolata sites (10%) (Figure 1b).

Sap flow density and environmental effects
The recorded mean daily maximum Js was in the range of 8.  Table 3). The mean Js of A. mangium and S. wallichii was two times larger than that of C. lanceolata in magnitude (Figure 3). At these three sites, the highest Js occurred between April and October during the wet season, coinciding with the highest T and PAR measured during the year. The lowest values were recorded between December and March during the dry season. The Js was similar between the dry and wet seasons for A. mangium, while Js was different between the dry and wet seasons for S. wallichii and C. lanceolata (Figure 3). We compared the analysis results regarding sensitivity of Js to environmental variables separately during the wet and dry seasons and in small and large trees. No significant difference in the correlation coefficient of Js with PAR or VPD between small and large trees of any species was found, and therefore, the results are not shown. Figure 4 demonstrates raw data (i.e. without using the statistical model) of Js and single environmental drivers during the wet and dry seasons. Tables 4, 5 and 6 show results from the mixed models that consider the effect of all the environmental drivers simultaneously, thus showing the actual dependencies between the studied variables when the effects of all environmental drivers are considered in the model. Table 4 shows the results for the daily average values (abbreviation daily data) for all the variables and Table 5 shows the hourly average values minus daily average values (abbreviation hourly data) for all the variables. Note that to obtain the response of canopy conductance per sapwood area to the environmental factors, 1 should be subtracted from coefficient β2 (compare Equation 3 and Table 4 and 5), while coefficients β1 and β3 for the canopy conductance per sapwood area are the same as for Js. The Js of A. mangium was affected on the daily scale by PAR, VPD, and SWC in both seasons, and it was more sensitive to PAR in the dry season than in the wet season (Table 4). SWC affected Js negatively during both seasons and in all soil layers (Tables 4 and 6). On the daily scale, Js of S. wallichii was sensitive to VPD but the parameter β2 was nearly zero in dry season ( Table 4). Note that this indicates strict canopy stomatal control as stomata closure occurs concurrently with increases in VPD, while transpiration (Js used as a proxy) remains unchanged (Equation 3). Moreover, Js of S. wallichii was significantly correlated with SWC in the two topmost soil layers, SWC1 and SWC2 (Table 6), and the sensitivity of Js to SWC was significant only during the dry season (Table 4). On the daily scale, Js of C. lanceolata was strongly related to PAR and SWC but not affected by VPD in the dry season; in the wet season, it was positively affected by PAR and VPD ( Table 4). The Js of C. lanceolata showed no sensitivity to VPD and the parameter β2 was nearly zero, also indicating that its canopy conductance per sapwood area was very sensitive to VPD. Only SWC in the topmost soil layer (SWC1) significantly affected the Js ( Table 6). Js of S. wallichii and C. lanceolata was more sensitive to PAR but less sensitive to VPD in dry season than in wet season (Table 4). There were distinct differences among these three species in their sensitivity to the environmental factors. The Js of A. mangium was less sensitive to PAR and SWC but more sensitive to VPD than that of the other two species in the dry season, but less sensitive to VPD in the wet season (Table 4).
On the hourly scale, the results were similar to the daily scale, except for S. wallichii and C. lanceolata in dry season (see Table 4 vs. Table 5). The Js of S. wallichii and C. lanceolata was more sensitive to VPD but less sensitive to PAR on hourly than daily scale.
For all species the correlation of the model residuals with temperature was statistically significant in the dry season and for C. lanceolata in the wet season as well. The strongest correlation was observed in A. mangium ( Table 7). The effect of wind speed was also tested and no effect on Js was found. We compared the effects in one species between different seasons and in one season between different species based on the 95% confidence intervals (supplementary material Table S1). If they did not overlap, the difference was significant. The significant difference between seasons within one species is marked in bold. The significant difference between the species in one season is marked with a letter: Same letter indicates that the species behave in the same way and are not significantly different from each other. All estimates marked with an asterisk are significantly different from zero at p < 0.05. We compared the effects in one species between different seasons and in one season between different species based on the 95% confidence intervals (supplementary material Table S2). If they did not overlap, the difference was significant. The significant difference between seasons within a species is marked in bold. The significant difference between the species within a season is marked with a letter: The same letter indicates that the species behave the same way and are not significantly different from each other. All estimates marked with an asterisk are significantly different from zero at p < 0.05.

Discussion
We found that both native species (S. wallichii and C. lanceolata) differed clearly in their stomatal sensitivity (at the daily scale) from one season to another, while the exotic A. mangium did not (Table  4). In the dry season, the parameter β2 of both native species was close to zero (Table 4). This means tight stomatal control, i.e. high sensitivity of g to VPD in S. wallichii and C. lanceolate. Earlier studies had also shown a sensitive canopy stomatal regulation in some tropical trees, resulting in partial canopy stomatal closure at VPD >1 kPa [35]. This implies that tropical trees had a conservative water use strategy when exposed to short-term drought [36,37]. High sensitivity of canopy stomatal control to VPD has been linked to conditions of water stress [38,39]. This result could be explained by consequent drought stress in S. wallichii and C. lanceolata during the dry season. For these two species, the positive correlation between Js and SWC indicates that water availability was limiting their water use in the dry season (Table 4). In contrast, the difference of stomatal behavior between the wet and dry seasons was not clear for A. mangium (Table 4). This may be because increasing SWC had a negative effect on the Js of A. mangium during both seasons (Table 4). One possible reason for this could be that A. mangium was suffering from too high soil water content. During the wet season, the Js of two native species was more sensitive to VPD than that of A. mangium (Table 4), indicating that A. mangium was more sensitive in canopy stomatal control to VPD than the native species were.
On the hourly scale, the dynamics of actual transpiration rate, which responds quickly to changes in environmental conditions, might be considerably different from those of the sap flow rate due to the time lag between these two variables [40,41]. The response of Js to the environmental conditions were similar to the daily scale when analyzed on the hourly scale, except in the case of S. wallichii and C. lanceolata in dry season (see Table 4 vs. Table 5). The canopy conductance per sapwood area response of these species was much stronger in daily than in hourly scale during the dry season. One reason for this difference in behavior could be the role of internal water stores within the tree and the decoupling of transpiration rate from the sap flow rate [42]. S. wallichii and C. lanceolata seemed to be suffering from drought stress during dry season (Table 4), and diurnal usage and refilling of internal water stores maybe a strategy under drought condition [43]. It is also known that the role of internal water stores increases when the total water usage is low [44,45], and the lowest water usage occurred in S. wallichii and C. lanceolata in the dry season (see the intercept, parameter in (α), in Table 4).
The Js of S. wallichii and C. lanceolata was very sensitive to PAR during the dry season, even more sensitive than that of A. mangium (Table 4). Liu et al. (2014) [46] found the same result i.e., that the Js of C. lanceolata was more sensitive to PAR than VPD during the dry season in Eastern China. A possible explanation for this could be the lower light levels during the dry season in winter as the light-response curve of photosynthesis is steeper at low light levels [47]. The greater sensitivity of the Js of C. lanceolata to light during the dry season may also be linked to their suffering from drought stress, as drought stress has been found to increase the sensitivity of canopy stomatal control to light [48].
Increasing SWC had a negative effect on the Js of A. mangium during both seasons ( Table 4), suggesting that soil water content could have been too high for A. mangium. The Js was similar between the dry and wet seasons for A. mangium (Figure 3). One would expect the Js to be higher during the wet season because of the higher level of light. However, this was not the case and could be due to too much water in the soil in the wet season causing oxygen deficiency and preventing optimal root functioning [49]. It has been reported that A. mangium is intolerant of excessively wet soil and grows better on sites with good drainage [50]. We also found that Js of A. mangium was sensitive to temperature. This finding corresponds to the report of Booth and Hong (1991) [11] that low temperatures in winter reduce the growth of A. mangium in Guangzhou (Guangdong Province). The effect of SWC on Js was significantly positive for S. wallichii and C. lanceolata only during the dry season. For C. lanceolata in eastern China, Liu et al. (2014) [46] also found that SWC (at depths 5, 10, 15 and 40 cm) had positive effect on Js in winter (dry season).
We also analyzed the effect of SWC at different soil depths on Js (Table 6), because SWC at different layers makes a different contribution to plant water use due to a variety of root system depths for different tree species [51]. The Js of A. mangium was sensitive to SWC in all layers, but the native species (S. wallichii and C. lanceolata) were not sensitive to SWC in deeper soil layers ( Table 6). This indicates that A. mangium may have deeper roots than both native species and could extract water down to at least 0-150 cm depth. Although Pan et al. (1996) [52] described that the roots of A. mangium were shallow (most roots concentrated in the upper 28cm), several other authors' reports support our results. For example, the rooting depth of A. mangium has been observed to be deeper than 1 m [53] and the fine roots to reach down to 12 m in monospecific stands of A. mangium [54]. Fine roots of C. lanceolata trees have been reported to be mainly distributed in the 0-20 cm soil layer in subtropical China [55], which is consistent with that C. lanceolata may only extract water from the soil layer at 0-30 cm depth (Table 6), while S. wallichii seemed to use water from 30-60 cm depth (Table 6). Also, Zeng (1988) [56] reported that most of the roots of S. wallichii in south China were concentrated in the 30-50 cm soil layer. Such differences in rooting depth may also explain why Js was reduced in S. wallichii, and C. lanceolata but not in A. mangium (Figure 3) during dry season when soil water availability was low. A. mangium is likely the most deeply rooted among the three species studied, and therefore had the least response to dry season [57].
The mean daily maximum Js of A. mangium at 9-55 g m −2 s −1 was lower than the values reported for A. mangium at Sardinilla, Panama (81 g m −2 s −1 ) of larger mean DBH [58], but similar to the values measured in an A. mangium plantation in South China (11-80 g m −2 s −1 ) of similar DBH [12]. The maximum Js of S. wallichii was 17 g m −2 s −1 in DBH = 15.5 cm trees in the Lesser Himalayas of Central Nepal [59], which was similar to the daily maximum Js of S. wallichii (11-24 g m −2 s −1 ) in the same DBH at our site. The Js of C. lanceolata was much smaller than that of the two other species, and its peak value was less than half that of the other species (Table 3). The reason for the low Js of C. lanceolata might be that C. lanceolata is a coniferous species, which generally uses less water in relation to broadleaved trees [60]. Moreover, the canopy stomatal conductance has been found to be lower in C. lanceolata than in broad-leaved trees, which also leads to lower Js in C. lanceolata stems [61]. This could also be explained by the thick wax layer on the leaf surface of C. lanceolata [61]. Although Zhang et al. (2016) [62] found a similar Js of C. lanceolata with similar DBH to our results in Southwest China, the Js of C. lanceolata was much smaller than reported for other locations in southern China [61,63]. Moreover, Li and Ritchie (1999) [64] reported that C. lanceolata requires well-drained soil and that it will not grow well in wet conditions. Jøker (2000) [65] also found that the best growth of C. lanceolata is obtained on well-drained loamy soils. According to Li and Ritchie (1999) [64], there are only a few places in South China which offer suitable terrain for C. lanceolata to grow well. Our results of mean daily maximum Js in the range 9-55 g m −2 s −1 for the A. mangium site and 11-39 g m −2 s −1 for the S. wallichii site both fall within the lower range of 4-139 g m −2 s −1 reported for subtropical and tropical tree species [58,60,66]. Low Js was the result of continuously high air humidity at the study site (Figure 1e), with VPD rarely exceeding 1 kPa (Figure 1d), which was supported by earlier studies on different species in the same area [67] and on a perhumid tropical forest of Sulawesi, Indonesia [68].
The annual rainfall (1600+ mm/year) of the observation site is enough to keep SWC at over 20% even in the cool-dry season, which might not induce severe evapotranspiration depression for any of the species at present. However, long-term observations have shown that droughts are becoming frequent and more severe, leading to declining SWC in the extended dry season [24]. Hu et al. (2018) [26] also showed that the dry season extended from November-March during the 1964-1983 period to October-April during the 2001-2010 period in the Heshan area (our sites), and they found that the extended dry season reduced transpiration in Schima superba (the same family as S. wallichii). Based on our results, we speculate that the exotic species A. mangium, which prefers higher temperatures and drier soil, will become more suitable in Southern China, whereas the more drought-sensitive native species S. wallichii and C. lanceolata could suffer more severely from water stress under future meteorological and soil conditions.

Conclusions
We found distinct differences in how PAR, VPD and SWC controlled the water use of the three plantation species. S. wallichii and C. lanceolata differed clearly in their stomatal behavior and Js at the daily scale from one season to another, while A. mangium did not. S. wallichii and C. lanceolata may suffer from water stress in dry season, but for A. mangium there seemed to be even an excess amount of water in the soil in both seasons. Extended dry season in the research area will likely impact the ecosystem hydrologic cycle and thus the water use of the studied species in the future.
Supplementary Materials: The following are available online at www.mdpi.com/xxx/s1, Table S1: 95% confident interval of model relating daily sap flow density (Js) to environmental factors in different seasons and species.