The Variation in Water Consumption by Transpiration of Qinghai Spruce among Canopy Layers in the Qilian Mountains, Northwestern China

: It is important for integrated forest-water management to develop a better understanding of the variation of tree transpiration among di ﬀ erent canopy layers in the forests and its response to soil moisture and weather conditions. The results will provide insights into water consumption by trees occupying di ﬀ erent social positions of the forests. In the present study, an experiment was conducted in the Qilian Mountains, northwest China, and 13 trees, i.e., 4–5 trees from each one of dominant (the relative tree height (H R ) > 1.65), subdominant (1.25 < H R ≤ 1.65) and intermediate-suppressed (H R ≤ 1.25) layers) were chosen as sample trees in a pure Qinghai spruce ( Picea crassifolia Kom.) forest stand. The sap ﬂux density of sample trees, soil moisture of main root zone (0 to 60 cm) and meteorological conditions in open ﬁeld were observed simultaneously from July to October of 2015 and 2016. The results showed that (1) The mean daily stand transpiration for the study period in 2015 and 2016 (July–October), was 0.408 and 0.313 mm · day − 1 , and the cumulative stand transpiration was 54.84 and 40.97 mm, accounting for 24.14% (227.2 mm) and 16.39% (249.9 mm) of the total precipitation over the same periods, respectively. (2) The transpiration varied greatly among canopy layers, and the transpiration of the dominant and codominant layers was the main contributors to the stand transpiration, contributing 86.05% and 81.28% of the stand transpiration, respectively, in 2015 and 2016. (3) The stand transpiration was strongly a ﬀ ected by potential evapotranspiration (PET) and volumetric soil moisture (VSM). However, the transpiration of trees from the dominant and codominant layers was more sensitive to PET changes and that from the intermediate-suppressed layer was more susceptible to soil drought. This implied that in dry period, such as in drought events, the dominant and codominant trees would transpire more water, while the intermediate-suppressed trees almost stopped transpiration. These remind us that the canopy structure was the essential factor a ﬀ ecting single-tree and forest transpiration in the dryland areas.


Introduction
Transpiration is an important component of water consumption from vegetated ecosystems worldwide and usually accounts for a large proportion of water consumption [1][2][3][4][5], especially for forest ecosystems in the water-limited arid regions such as northwest China [5,6]. Under the decade-long

Research Plot
In 2015, a 25 × 50 m 2 plot was set up in a pure spruce plot (38°33′12.6″ N, 100°17′7.9″ E; 2762 m a.s.l.), located on a northwest-facing slope with a slope gradient of 27° (Figure 1c). The forest in the plot consisted of Qinghai spruce ranging from 45 to 115 years in age [27]. The canopy density, i.e., the ratio of the total projected area of all trees to the total stand area, was 0.55. The tree density was 1120 trees·ha −1 . The tree height (H) ranged from 5.0 to 24.3 m, with an average of 10.6 ± 5.6 m. The tree diameter at breast height (DBH) ranged from 5 to 48 cm, with an average of 15.3 ± 9.1 cm. The canopy width ranged from 1.4 to 6.7 m, with an average of 3.6 ± 1.1 m. The leaf area index (LAI) ranged from 2.53 to 3.96 during the growing season (May-October), using the plant canopy analyzer (LAI-2200C, Li-Cor Inc., Lincoln, NE, USA).
The vegetation overstory was composed of only Qinghai spruce trees in the plot, but it had different canopy layers in vertical direction, i.e., different social positions for trees in the forest. Based In the Pailugou watershed, along the elevation gradient, the vegetation types are mountain forest grasslands, subalpine shrub meadows, and alpine meadows, respectively. The corresponding soil types are chestnut soil, grey cinnamon soil, subalpine shrub meadow soil, and alpine meadow soil, respectively [26]. The dominant tree species of mountain forest is Qinghai spruce (Picea crassifolia Kom.), which is distributed only on shaded (north-facing) or semi-shaded (northwest-and northeast-facing) slopes at elevations from 2500 to 3300 m a.s.l. Shrubs community consist mainly of Caragana jubata (Pall.) Poir., Potentilla fruticosa L., Berberis diaphana Maxim. Herbaceous vegetation consists mainly of Carex lanceolata Boott, Pedicularis spp., Stipa capillata.

Research Plot
In 2015, a 25 × 50 m 2 plot was set up in a pure spruce plot (38 • 33 12.6" N, 100 • 17 7.9" E; 2762 m a.s.l.), located on a northwest-facing slope with a slope gradient of 27 • (Figure 1c). The forest in the plot consisted of Qinghai spruce ranging from 45 to 115 years in age [27]. The canopy density, i.e., the ratio of the total projected area of all trees to the total stand area, was 0.55. The tree density was 1120 trees·ha −1 . The tree height (H) ranged from 5.0 to 24.3 m, with an average of 10.6 ± 5.6 m. The tree diameter at breast height (DBH) ranged from 5 to 48 cm, with an average of 15.3 ± 9.1 cm. The canopy width ranged from 1.4 to 6.7 m, with an average of 3.6 ± 1.1 m. The leaf area index (LAI) ranged from 2.53 to 3.96 during the growing season (May-October), using the plant canopy analyzer (LAI-2200C, Li-Cor Inc., Lincoln, NE, USA).
The vegetation overstory was composed of only Qinghai spruce trees in the plot, but it had different canopy layers in vertical direction, i.e., different social positions for trees in the forest. Based on the relative tree height (H R ), that was, the ratio of the height of each tree to the average tree height in the plot [20,[30][31][32], all 140 Qinghai spruce trees in this research plot were divided into three layers in the vertical direction, i.e., the dominant (H R > 1.65), codominant (1.25 < H R ≤ 1.65) and intermediate-suppressed layers (H R ≤ 1.25). The dominant layer, the top canopy layer, was composed of dominant trees, which were defined as the tallest trees in the plot. The codominant layer, the second canopy layer, was composed of codominant trees, which constituted a group of similar trees whose tree height were slightly lower than dominant trees. The intermediate-suppressed layer, the third canopy layer, was composed of intermediate and suppressed trees, which were located under larger trees, and their crowns were mostly or totally overlapped by neighboring crowns [30][31][32].
Trees from the dominant layer were the tallest trees in the plot, with an average tree height of 18.1 ± 2.9 m. There were 21 dominant trees, with only 15% of all trees in the plot, but their total sapwood area was 8735.32 cm 2 , accounting for 41.7% of the total plot sapwood area. The average tree height of the codominant and intermediate-suppressed layers was much smaller than that of the dominant layer, with decreases of 4.2 and 11.4 m, respectively. The number of the intermediate-suppressed layer trees was the largest, 84 trees, 4 times that of the dominant layer, but the total sapwood area was the smallest, with a value of 5455.45 cm 2 , accounting for 26.1% of the total plot sapwood area (Table 1). According to the relationship between the sapwood area and tree diameter at breast height (DBH) [6] in this Qinghai spruce stand, the sapwood area (A s , mm 2 ) per Qinghai spruce tree was accurately calculated using Equation (1): where DBH (mm) is the tree diameter at breast height. According to the relationship between the leaf area and DBH and tree height (H) [33] in this research area, the leaf biomass (W L , kg) per Qinghai spruce tree was calculated using Equation (2): The dominant shrub in the plot was mainly Potentilla glabra Lodd., and the average height was 0.65 ± 0.05 m at a coverage of approximately 4%. The herbaceous vegetation was mainly Carex lanceolata Boott, and the average height was 5.8 ± 2.2 cm at a coverage of approximately 17%. Moss covered the forest floor, which was Abietinella abietina, and the average thickness was 8.4 ± 3.1 cm, while the coverage was approximately 30%.
The soil in the plot was thick and highly permeable, with an average soil thickness of 70 cm. The soil type was mountainous grey cinnamon soil and its texture was loam, with sand content of 34%, a total soil porosity of 71.2%, and a bulk density of 0.83 g·cm −3 in the root zone (0 to 60 cm soil layer) [28,34]. Within the 0 to 60 cm root zone, the saturation moisture content was approximately 0.69 m 3 ·m −3 and the field capacity was approximately 0.58 m 3 ·m −3 [28,34].

Sap Flow Measurement and Transpiration Calculation
In order to describe the variation in transpiration among different canopy layers in a Qinghai spruce plot, four or five representative trees (i.e., healthy, average shape) from each canopy layer were chosen as sample trees according to the percent of total sapwood area of each canopy layer in the plot, which characters varied with the diameter at breast height (DBH) from 10 to 38 cm, the tree height from 6.1 to 20.5 m, the canopy width from 2.4 to 6.7 m, the canopy thickness of the difference between the tree height and the height of the first live branch from 4.0 to 15.3 m, the sapwood area from 6741 to 53,396 mm 2 , and the leaf biomass from 5.67 to 32.06 kg ( Table 2). The sap flow density of thirteen sample trees was measured at breast height (1.3 m above ground level) using thermal diffusion probes (SF-L, Ecomatik, Munich, Germany) consisting of four sensors that were 20 mm in length and 2 mm in diameter (S 0 , a heated sensor powered by a constant current at a voltage of 12 V; S 1 , an unheated sensor; and S 2 and S 3 , reference sensors). These sensors were inserted into the outer 20 mm of the xylem at breast height on the north-facing side of trunks (not the sun-exposed side). Before insertion, the outer bark was peeled off, and each probe was coated with silicone gel to ensure good thermal contact between the probe elements and sapwood. After insertion, the exposed cambium was covered with silicone gel to reduce evaporation from the wood surface and then covered with aluminum foil to avoid physical damage and thermal influences from solar radiation. Data were recorded by a data logger (DL2, Delta-T Devices, Cambridge, UK) every 10 min.
The sap flow density (J s , mL·cm −2 ·min −1 ) (flow per unit of sapwood area) of an individual sample tree was calculated by Equations (3) and (4) [4,13]: where d tact is the difference in instantaneous temperature ( • C) between the heated and reference sensors; d tmax is the value of d tact at which the sap flow density is nil or near zero; and T 1-0 , T 1-2 , and T 1-3 are the temperature differences between sensors S 1 and S 0 , S 1 and S 2 , and S 1 and S 3 , respectively. The azimuthal variation in sap flux density of each tree was related with tree crown architecture [35] and crown exposure to the sun [36,37], both of which were thought to be closely related to the canopy density [5]. When the canopy was highly closed, the azimuthal variation in sap flux density was thought to be negligible [5]. In this study, highly closed forests (with canopy densities of 0.55) was chosen as research plot, thus we assumed the azimuthal variation was small for each tree.
For Qinghai spruce trees, their sapwood depth normally was less than 40 mm and its sap flow density varied with sapwood depth by the Gaussian distribution function with a peak at the depth of 20 mm [29]. In this study, our probes were only 20 mm long, which could not measure the sap flow density within the sapwood deep exceeding 20 mm [13]. Thus, the J s in the 20 to 40 mm was calculated by Gaussian regression equation (Equation (5) where d (mm) is sapwood depth from the cambium, and J s (Kg·day −1 ) is daily sap flow density at d mm depth in the sapwood. To quantify the daily and monthly variation in transpiration among canopy layers, the daily transpiration (F i , mm·day −1 ) of each canopy layer per unit area was calculated by Equation (6): where J ci (mL·cm −2 ·min −1 ) is the average J s value of the sample tree of canopy layer i, with i = 1, 2 or 3, representing the dominant, codominant and intermediate-suppressed layers, respectively; A ci is the cumulative sapwood area of all trees in canopy layer i in the plot; and S (cm 2 ) is the projected area of the plot. The daily stand transpiration (T, mm·day −1 ) per unit area was calculated by Equation (7): Qinghai spruce trees began growing on approximately 25 May 2015 and 20 June 2016, and ended growing on approximately 21 August 2015 and 14 September 2016 [28]. The better observation on transpiration should be carried on the whole growing season (May-October). Unfortunately, we could not collect the data throughout the whole growing season, although we carried out the observation in two growing seasons, since the thermal diffusion probes did not work and record the data on the early period of growing season. This may be due to that the lower temperature leaded to the unstable electricity supply and worse data. Therefore, the whole experiment was conducted from July to October of 2015 and 2016.
The vapor pressure deficit (VPD, kPa) was calculated from air temperature (T, • C) and relative humidity (R h , %) using Equation (8) [38]: According to the Food and Agriculture Organization of the United Nations (FAO) Penman-Monteith equation, the daily potential evapotranspiration (PET, mm·day −1 ), as a comprehensive meteorological index of many weather parameters (air temperature, solar radiation, VPD, and wind speed), can reflect the evaporative demand determined by the air conditions and available energy and thus has an important effect on transpiration [13]. The daily PET was used to estimate by using the Equation (9) [39]: where ∆ (kPa· • C −1 ) is the slope of the vapor pressure versus temperature relationship; R n (MJ·m −2 ·day −1 ) is the net radiation; G (MJ·m −2 ·day −1 ) is the soil heat density; γ (0.053 kPa· • C −1 ) is the psychrometric constant; T ( • C) is the mean air temperature; U (m·s −1 ) is the wind speed; e s (kPa) is the saturation vapor pressure; and e a (kPa) is the actual vapor pressure.
In 2015 and 2016 (July-October), the volumetric soil moisture of main root zone (the 0 to 10 cm, 10 to 20 cm, 20 to 40 cm, and 40 to 60 cm soil layers) was continuously monitored at a single location by using soil moisture sensors (5-TE, Decagon, Pullman, WA, USA). Data were collected every 10 min by Due to the lack of facilities, nine points were evenly distributed within the plot at fixed positions to revise VSM 1 . Soil samples were collected at each point with a time interval of approximately 10 days, by using an aluminum box with a volume of 100 cm 3 , from the 0 to 10 cm, 10 to 20 cm, 20 to 40 cm, and 40 to 60 cm soil layers, to calculate the mean volumetric soil moisture of the 0 to 60 cm soil layer (VSM, m 3 ·m −3 ). Based on these data, VSM was estimated by using the following relationship (Equation (10)): .96, n = 12, and P < 0.01 (10)

Data Analysis
In this study, the stand transpiration was calculated based on the cumulative transpiration of the four canopy layers, and the transpiration of each canopy layer was calculated based on the sapwood area and the measured sap flow density. The transpiration response to PET and volumetric soil moisture (VSM) were analyzed on daily basis using the measured data. The transpiration model was determined by polylinear regression analysis using the Statistical Product and Service Solutions (SPSS) version 21.0 (IBM Inc., Chicago, IL, USA).
The whole experiment was conducted from July to October of 2015 and 2016. Unfortunately, we could not collect the whole data in this period since the equipment abruptly failed and unstable electricity supplied. Therefore, to calculate the daily and monthly cumulative transpiration, the missed data was reasonably interpolated by the linear relationship between daily transpiration of each canopy layer and daily potential evapotranspiration (PET) built in this study (Table 3). To determine the impact of soil moisture limitations on transpiration, we divided VSM into two groups (VSM ≥ 0.21 m 3 ·m −3 and VSM < 0.21 m 3 ·m −3 ). VSM = 0.21 m 3 ·m −3 was the threshold value, which was determined based on the deviation from linearity of the relationship between PET and stand transpiration. Well-watered soil conditions occurred when VSM was more than 0.21 m 3 ·m −3 , while water-limited soil conditions occurred when VSM was less than 0.21 m 3 ·m −3 .
The different transpiration responses to PET and VSM among the canopy layers were analyzed on a daily basis. The standardized coefficients of the regression equation [40,41] could reflect the important degree of transpiration response to environmental factors. In addition, the larger the standardized coefficients was, the stronger the relationship between transpiration and environmental factors was [42].
To further examine the impact of soil moisture on transpiration in the period under water-limited soil conditions, we defined the reduction in transpiration as follows (Equation (11)):

Environmental Characteristics during the Growing Season
In the Pailugou watershed, precipitation variability was high from July to October of 2015. There were 45 rainy days with a total precipitation of 227.2 mm, of which 82.2% of rainy days was less than 10 mm·day −1 . Rainy days occurred mainly from 1 to 8 July, from 1 to 13 August, and from 2 to 11 September, with a total precipitation of 77.5, 31.0 and 60.4 mm, respectively (Figure 2a). However, precipitation was distributed more evenly from July to October of 2016. There were 60 rainy days with total precipitation of 249.9 mm in 2016 (July-October). 86.7% of the rainy days occurred mainly from July to September (Figure 2b).

Stand Transpiration and Its Vertical Distribution among Canopy Layers
The stand transpiration fluctuated with the daily PET. The stand transpiration varied in the ranges of 0.047-1.098 mm·day −1 and 0.082-0.870 mm·day −1 , with averages of 0.408 and 0.313 mm·day −1 , respectively, from July to October of 2015 and 2016 (Figure 2k,l).     On cloudy (e.g., 19 August 2015) and rainy days (e.g., 7 August 2015), the daily stand transpiration was much lower than that on a sunny day (e.g., 4 August 2015), decreasing by 0.36 and 0.69 mm, respectively ( Figure 3a); the daily stand transpiration significantly decreased for the dominant and codominant layers, where the intermediate-suppressed layer only slightly declined. However, the transpiration ratio for the three canopy layers showed little variation with the changing weather conditions; the ratio was high for the dominant and codominant layers, accounting for 81.30% and 80.64%, respectively, of the stand transpiration on cloudy and rainy days, and the ratio of the intermediate-suppressed layer was low, accounting for 18.70% and 19.36%, respectively, of the stand transpiration on cloudy and rainy days (Figure 3b).
The stand transpiration varied with each month. From July to October in 2015 and 2016, the cumulative stand transpiration was 54.84 and 40.97 mm, accounting for 24.14% and 16.39% of the total precipitation over the same period, respectively. In 2015, the cumulative stand transpiration monotonously decreased from 18.75 mm in July to 9.15 mm in October, and it exhibited little variation in September and October, with a difference of only 0.27 mm (Figure 4a). However, in October 2016, the stand transpiration was lower than that in other months (Figure 4b).

Daily Stand Transpiration Responses to PET and VSM
Transpiration is strongly dependent on potential evapotranspiration (PET) and volumetric soil moisture (VSM) [4,13,16,43], especially in arid and semiarid areas where the soil moisture is rather low and the evaporative demand is high [4,6].
In the period of well-watered soil conditions (VSM ≥ 0.21 m 3 ·m −3 ) and rainy days, the daily stand transpiration responded to the daily PET following a linear relation ( Figure 5); that is, the daily stand transpiration increased with increasing daily PET according to Equation (12). 0.120 PET 0.008 p T    (R 2 = 0.87, n = 66, and P < 0.01) (12) where Tp (mm·day −1 ) is the simulated stand transpiration, and PET (mm·day −1 ) is the daily potential evapotranspiration.
In the period of water-limited soil conditions (VSM < 0.21 m 3 ·m −3 ), there are clear deviations from linearity in the relationship between PET and stand transpiration ( Figure 5). Figure 6 shows the difference between the measured transpiration (T) and simulated transpiration (Tp), and a substantial negative deviation occurred during this period, indicating that transpiration was limited by soil moisture. When VSM was less than 0.21 m 3 ·m −3 , the daily stand transpiration responded to VSM following a linear relation except on rainy days ( Figure 7); that is, the daily stand transpiration decreased with the decreasing VSM based on Equation (13).
where VSM (m 3 ·m −3 ) is the mean volumetric soil moisture of the 0 to 60 cm soil layer.
On rainy days, the measured and calculated transpiration values were very similar from 1 to 13 August and from 2 to 6 September in 2015, even when VSM remained low ( Figure 6). Precipitation creates a temporary water reservoir that can be absorbed by trees in an uninhibited fashion [16].

Daily Stand Transpiration Responses to PET and VSM
Transpiration is strongly dependent on potential evapotranspiration (PET) and volumetric soil moisture (VSM) [4,13,16,43], especially in arid and semiarid areas where the soil moisture is rather low and the evaporative demand is high [4,6].
In the period of well-watered soil conditions (VSM ≥ 0.21 m 3 ·m −3 ) and rainy days, the daily stand transpiration responded to the daily PET following a linear relation ( Figure 5); that is, the daily stand transpiration increased with increasing daily PET according to Equation (12).
T p = 0.120 × PET − 0.008 R 2 = 0.87, n = 66, and P < 0.01 (12) Forests 2020, 11, 845 11 of 17 where T p (mm·day −1 ) is the simulated stand transpiration, and PET (mm·day −1 ) is the daily potential evapotranspiration.   In the period of water-limited soil conditions (VSM < 0.21 m 3 ·m −3 ), there are clear deviations from linearity in the relationship between PET and stand transpiration ( Figure 5). Figure 6 shows the difference between the measured transpiration (T) and simulated transpiration (T p ), and a substantial negative deviation occurred during this period, indicating that transpiration was limited by soil moisture. When VSM was less than 0.21 m 3 ·m −3 , the daily stand transpiration responded to VSM following a linear relation except on rainy days ( Figure 7); that is, the daily stand transpiration decreased with the decreasing VSM based on Equation (13).
where VSM (m 3 ·m −3 ) is the mean volumetric soil moisture of the 0 to 60 cm soil layer.     On rainy days, the measured and calculated transpiration values were very similar from 1 to 13 August and from 2 to 6 September in 2015, even when VSM remained low ( Figure 6). Precipitation creates a temporary water reservoir that can be absorbed by trees in an uninhibited fashion [16].

Different Responses of the Daily Transpiration to PET and VSM among Canopy Layers
In the period of well-watered soil conditions and rainy days, the daily transpiration of the dominant and codominant layers increased linearly with increasing PET, and the coefficients of determination (R 2 ) of the regression relation were 0.880 and 0.797, respectively. Conversely, the daily transpiration of the intermediate-suppressed layer hardly changed with PET (R 2 = 0.613). The difference in transpiration among the three canopy layers increased rapidly with increasing PET (Figure 8). The relationship between the dominant-layer transpiration and PET was more sensitive than that of any other layer, and the standardized coefficient (K d ) of the regression equation was 0.938 (Table 3).

Different Responses of the Daily Transpiration to PET and VSM among Canopy Layers
In the period of well-watered soil conditions and rainy days, the daily transpiration of the dominant and codominant layers increased linearly with increasing PET, and the coefficients of determination (R 2 ) of the regression relation were 0.880 and 0.797, respectively. Conversely, the daily transpiration of the intermediate-suppressed layer hardly changed with PET (R 2 = 0.613). The difference in transpiration among the three canopy layers increased rapidly with increasing PET (Figure 8). The relationship between the dominant-layer transpiration and PET was more sensitive than that of any other layer, and the standardized coefficient (Kd) of the regression equation was 0.938 (Table 3).
In the period of water-limited water conditions, the daily transpiration decreased rapidly with decreasing VSM for the dominant and codominant layers, which declined slightly for the intermediate-suppressed layer (Figure 8). The R 2 of each regression relation exceeded 0.95, but the intermediate-suppressed layer transpiration responses to VSM was more notable than that of any other layer, and the standardized coefficient was 0.980.    In the period of water-limited water conditions, the daily transpiration decreased rapidly with decreasing VSM for the dominant and codominant layers, which declined slightly for the intermediate-suppressed layer (Figure 8). The R 2 of each regression relation exceeded 0.95, but the intermediate-suppressed layer transpiration responses to VSM was more notable than that of any other layer, and the standardized coefficient was 0.980.

Tree Dominance Determines Its Role in Forest Transpiration
In related studies, trees were distinguished by DBH, and larger trees were found to generally transpire more than smaller trees [24,26], but it was also found that the mean sap flow of trees was not strongly correlated with DBH [18,22,44]. In addition to DBH, tree dominance in forest was also the main influencing factor on transpiration [22], because it was related to the ability to acquire living resources, including water and energy resources. To quantify tree dominance in forests, the forest canopy could be divided into three layers, i.e., the dominant, codominant and intermediate-suppressed layers according to the difference of tree height between trees and their neighboring trees [31,32].
Trees from the dominant layer had a higher sap flow density than that of trees from the suppressed layer [20,21,45]. In this study, the dominant-layer trees were further proven to transpire more, and the dominant layer was the largest contributor to the stand transpiration, which accounted for 55.65% and 55.75% of the cumulative stand transpiration, respectively, in 2015 and 2016 (Figure 4), although only 15% of the trees in the plot occurred in the dominant layer (Table 1). This occurred due to a larger conductive sapwood area and leaf biomass of the dominant-layer trees of a forest located in the upper layer (Table 1), and to the intense solar radiation and large vapor pressure deficit in the upper layer [21,32,46]. This key contribution of dominant-layer trees was also found by in previous research, such as the study of Ceulemans et al. [47] in an Abies amabilis stand in western Washington, USA, and that of Sun et al. [48] in a Pinus koraiensis plantation, Northeast China.

Effects of PET and VSM on the Daily Transpiration among Canopy Layers
Transpiration in forest stands is mainly controlled by PET and VSM [4,13,16,38,43]. The dominant and codominant layers occupy favorable social positions in the forest and generally obtain more resources, especially light and soil moisture [32]. They are usually considered to be less vulnerable to drought [46], i.e., PET increase and soil water decrease.
However, in this study, it was found that the dominant and codominant layers were more sensitive to the changing PET and VSM. In the period of well-watered soil conditions and rainy days of 2015, the daily transpiration of the dominant and codominant layers increased linearly with increasing daily PET; conversely, the daily transpiration of the intermediate-suppressed layer hardly changed with daily PET (Figure 8). Similarly, in 2016, the transpiration of each canopy layer increased linearly with increasing daily PET, but the dominant and codominant layers were more sensitive to PET. Trees for the dominant and codominant layers lie in the upper and second layers of the forest, and take more favorable the spatial positions in the stand, greater exposure to light, and the stronger driving force that might increase sap flow than the intermediate-suppressed trees [16,20,21].
In the period of water-limited soil conditions, the daily transpiration of the dominant and codominant layers decreased rapidly with the soil moisture decline, while that of the intermediate-suppressed layer declined slightly (Figure 8). Trees from the intermediate-suppressed layer seem to be limited more by soil moisture and less sensitive to PET. Moreover, under the water-limited soil conditions, the average daily transpiration was reduced by 35.6%, 29.8% and 46.5% in the dominant, codominant and intermediate-suppressed layers, respectively, compared to the transpiration at the same PET in the period of well-watered soil conditions (estimated by Equation (11)). This occurred due to the shallow rooting system of the spruce trees and the soil thickness of no more than 60-100 cm in the forest in this study. For Qinghai spruce trees, with up to 95% of their active roots (fine roots for water and nutrient uptake) distributed within the upper 45 cm of the soil profile [49], the tree roots of the intermediate-suppressed layer were undeveloped, so that small trees could not absorb more soil water from the deeper soil layers and had to limit their water consumption, i.e., transpiration, when soil water was limited.
In this study, we only analyzed the variation of transpiration among canopy layers under the different weather and month conditions, and its responses to PET and VSM from July to October of 2015 and 2016. Therefore, it is difficult to build a transpiration model for all year due to limited data. Thus, it is necessary to follow the transpiration of Qinghai spruce with more years of observations and more sample trees with different canopy layers in future studies, to develop an integrated model describing the complex response of transpiration to the changing environment. In addition, there was the obviously vertical variations in weather conditions and forest structure (PAR, VPD, temperature, LAI) in a natural forest. Thus, in future studies, we will pay more attention to the conditions for different social positions, i.e., the vertical variations in weather condition and their effects on the transpiration of different canopy layers.

Implications for Tree Water Use
In this study, the transpiration of the intermediate-suppressed trees almost stopped when the soil water content below 0.15 m 3 ·m −3 , whereas the transpiration of the dominant and codominant trees became smaller but continued. This seems to suggest that, when severe drought occurs, the dominant and codominant trees can take water and keep photosynthesis activity with controlling transpiration to some extent, while the intermediate-suppressed trees are more likely to suffer from water stress. The intermediate-suppressed trees are generally shaded by the dominant and codominant trees from the more intense solar radiation and VPD, and the shade might enable intermediate-suppressed trees to survive, although severe soil drought has a strong impact on water use in the intermediate-suppressed trees. If extreme drought occurs after cutting the dominant and codominant trees, suppressed trees would be exposed to higher atmospheric evaporative demand and would also suffer from more severe soil drought because of an increase in forest-floor evaporation in a prolonged dry season [46]. This study reminds us that the upper trees need more water to be alive in any water conditions.

Conclusions
This study revealed the stand transpiration of Qinghai spruce and its variation among canopy layers and how PET and VSM affected transpiration at stand and canopy layer levels in the Qilian Mountains, Northwestern China. From July to October in 2015 and 2016, the cumulative stand transpiration was 54.84 and 40.97 mm. The largest contributor to the stand transpiration was the dominant layer, accounting for 55.65% and 55.75% of the cumulative stand transpiration, respectively, in 2015 and 2016. In another word, the larger trees were larger than their neighbors, which occupied a favorable social position in the forest, and transpired more.
The transpiration was sensitive to variations in PET and VSM in the Qinghai spruce forest. When the mean VSM of the 0 to 60 cm soil layer dropped below 0.21 m 3 ·m −3 , the daily stand transpiration was often limited by VSM, and it declined according to a significant positive linear correlation with VSM. However, the daily stand transpiration was limited by the daily PET in the period of well-watered soil conditions (VSM ≥ 0.21 m 3 ·m −3 ) and rainy days, and it increased following a significant positive linear correlation with the daily PET. Moreover, the daily transpiration was more sensitive to PET changes in the dominant and codominant layers than in the intermediate-suppressed layers, while trees from the intermediate-suppressed layer seem to be limited more by soil moisture. Overall, the change in stand transpiration was mainly influenced by the dominant and codominant layers. Thus, a sustainable forest structure, considering the canopy layer, is important for forest-water management in the dryland regions.