Leaf Phenology Variation within the Canopy and Its Relationship with the Transpiration of Populus tomentosa under Plantation Conditions

To provide a theoretical basis for developing intensive cultivation practices for Populus tomentosa plantations, the leaf phenology variation within the canopy and its relationship with transpiration of well-watered P. tomentosa tree in a dense plantation were investigated. The variation in canopy light interception, indicated by the ratio between net radiation under (Rn-u) and above (Rn-a) canopy, with leaf development was also studied. During the growing season, the achievement of maximum leaf number tended to be later in higher parts of the canopy. In the lower and middle canopy, the leaf number maximized earlier in the east-facing side than in the west-facing side, but this difference disappeared in the upper canopy. The Rn-u/Rn-a was stable in May, but declined and then varied steadily until late August. Generally, in May, the crop coefficient (Kcb) of the tree reached its highest level and was not correlated with leaf area (LA) in all layers (p > 0.05). However, it increased linearly (p < 0.001) with LA in the layers above a canopy height of 3 m from June to late August, and most of its variation was explained by LA in the 5–7 m layer. After late August, Kcb decreased linearly with decreasing LA in all layers (p < 0.001). Consequently, a temporal ecological strategy seems to be adopted by P. tomentosa leaves in different layers and azimuthal sides for efficient light acquisition. The contribution of the different canopy layers to tree transpiration can vary, with the leaves in the upper and all layers mainly controlling transpiration in summer and in spring and autumn, respectively.


Introduction
Poplar has been identified as one of the fastest growing trees in temperate zones, and is also one of the main tree species for developing fast-growing and high-yielding plantations in North China.Plantations of poplar can not only produce a large amount of wood and fiber resources, but also provide ecosystem and environmental services, such as watershed services, waste management, carbon sink, biodiversity and wildlife habitat, etc.At present, the poplar timber plantation area has reached 4.82 million ha, accounting for 7.66% of the total timber plantation area in China [1] and Forests 2018, 9, 603 2 of 17 ranking first in the world.However, the productivities and qualities of these poplar plantations are still very low, and they need to be improved urgently by environmentally friendly intensive cultivation practices, such as water-saving irrigation, fertilization, pruning, and weed control [2][3][4].
Understanding of the growth rhythm of the trees is the basis to improve plantation productivities by intensive cultivation practices that can be refined according to the characteristics of tree growth [5,6].The leaves are the primary photosynthetic organ for producing carbohydrate in most cases, so its phenology is closely related to tree's growth rhythm.Additionally, leaf phenology can alter the interception and allocation of light, microclimate, and the potential carbon-gain strategy within canopy [7][8][9][10][11].Thus, information on leaf phenology within the canopy will help to gain more profound insights into tree growth characteristics and to predict forest productivity [12].
However, since it is very difficult, time-consuming and laborious to closely observe the leaf phenology of trees in the field, most researchers have investigated the leaf development of forests by indirectly measuring leaf area index [13][14][15].There are also some researchers who did pioneering work on observing leaf dynamics in forest with the help of binoculars [16], mobile walk-up towers [17], meteorological towers [18] and a canopy walkway system [12].Although these studies investigated leaf dynamics in different canopy layers, their observations were mainly limited to variation in the areas of leaves from some representative branches, or branch units [9,16,17].This will definitely result in some bias of the actual leaf phenology within the canopy, because of the high heterogeneity in light availability along the canopy profile [19] which can result in tremendous variability in leaf development within the canopy [16].Therefore, observing the dynamics of all leaves is the most precise way to understand leaf phenology within the canopy.However, to our knowledge, no attempt has been conducted for such detailed work in forest plantations.In addition, previous works have provided deep insights into the vertical variation of leaf phenology within canopy [9,16,17,20], but the changes in the leaf phenology with azimuth at different canopy levels remain unclear and have not been investigated.
Knowledge of tree water use is also pivotal for developing efficient cultivation practices for promoting forest productivity and quality [5,21] and it can help to predict tree growth [22,23].Thus, a huge number of researchers have investigated the water use of plantations of many tree species [13][14][15]22].Whole tree water use is the summation of transpiration rates of all the leaves within the canopy, so transpiration of a tree is closely related to, and can be determined by, the leaf phenology and physiology in its canopy.At present, the control effects of stand leaf dynamics on tree water use have been well studied and understood by correlating transpiration of trees to stand leaf area index [14,15,24,25].However, the response of tree water use to the leaf areas in different canopy layers has been little investigated.Thus, it is still not clear how leaves in different canopy layers influence the seasonal transpiration variation of the trees in plantations.
The Chinese white poplar, Populus tomentosa Carr.(Salicaceae), plays a critical role in poplar timber plantations in the North China Plain [26], but the average actual plantation productivity of P. tomentosa in this region is still far lower than its potential productivity [5,27].Consequently, in order to provide a theoretical basis for developing high-efficiency intensive cultivation measures for P. tomentosa plantations, based on leaf phenology and tree water use, this study was conducted with three objectives: (1) to determine the leaf phenological variation with canopy height and azimuth in different canopy layers of P. tomentosa, (2) to investigate the variation of canopy light interception with leaf development in the canopy profile, and (3) to reveal the relationships between transpiration and the leaf areas in different canopy layers.To achieve these objectives, only one sample tree in a mature P. tomentosa plantation was selected as the research focus, due to the great difficulty in investigating the leaf phenology of a big tree that requires counting all the leaves within its canopy.However, although using only one sample tree for study may introduce a high risk of random effect, it provides the opportunity to thoroughly understand the leaf phenology variation within the canopy by counting all of the leaves.Besides, this method can also help to directly link the whole tree transpiration to the leaf development in different canopy parts.

Experimental Site Description and Experimental Plantation
The experimental site was located in Gaotang County, Liaocheng City, Shandong Province, China (36 • 58 N, 116 • 14 E, and elevation 27 m).In this area, the climate is a warm temperate monsoon.The average annual precipitation, air temperature, and free water-surface evaporation were 545 mm, 13.2 • C and 1880 mm, respectively.The deep soil, developed from quaternary alluvium, had a texture of silt and silt loam over the 0-120 and >120 cm soil depths, respectively.In the upper 40 cm layer, the soil has an average bulk density of 1.63 g cm −3 , pH of 8.1, 41.5 mg kg −1 available N, 76.8 mg kg −1 available K, 7.11 mg kg −1 available P, and 0.94% organic matter.
The experimental plantation was established in late March 2005 with the clone B301 of P. tomentosa.The trees were planted with an alternate narrow-(2 m) and wide-row (6 m) spacing scheme, with within-row spacing of 1 m (Figure 1).This led to a planting density of 2500 trees ha −1 .From 2010 to 2011, a subsurface drip irrigation experiment was installed in this plantation, and the experiment had three irrigation treatments, and one control non-irrigation treatment with each, including three replications [5].The present study was only conducted in the second replicated plot of the T25 treatment in 2010.In this treatment, the trees were irrigated when the average tensiometer-measured soil water potential at the 20 cm depth and 10 cm distant from a drip emitter reached −25 kPa.Due to the onset of the rainy season, and the high (shallow) groundwater level afterward [5], irrigation was only implemented from early April to early August, with total irrigation amount of about 336 mm.

Experimental Site Description and Experimental Plantation
The experimental site was located in Gaotang County, Liaocheng City, Shandong Province, China (36°58′ N, 116°14′ E, and elevation 27 m).In this area, the climate is a warm temperate monsoon.The average annual precipitation, air temperature, and free water-surface evaporation were 545 mm, 13.2 °C and 1880 mm, respectively.The deep soil, developed from quaternary alluvium, had a texture of silt and silt loam over the 0-120 and >120 cm soil depths, respectively.In the upper 40 cm layer, the soil has an average bulk density of 1.63 g cm −3 , pH of 8.1, 41.5 mg kg −1 available N, 76.8 mg kg −1 available K, 7.11 mg kg −1 available P, and 0.94% organic matter.
The experimental plantation was established in late March 2005 with the clone B301 of P. tomentosa.The trees were planted with an alternate narrow-(2 m) and wide-row (6 m) spacing scheme, with within-row spacing of 1 m (Figure 1).This led to a planting density of 2500 trees ha −1 .From 2010 to 2011, a subsurface drip irrigation experiment was installed in this plantation, and the experiment had three irrigation treatments, and one control non-irrigation treatment with each, including three replications [5].The present study was only conducted in the second replicated plot of the T25 treatment in 2010.In this treatment, the trees were irrigated when the average tensiometermeasured soil water potential at the 20 cm depth and 10 cm distant from a drip emitter reached −25 kPa.Due to the onset of the rainy season, and the high (shallow) groundwater level afterward [5], irrigation was only implemented from early April to early August, with total irrigation amount of about 336 mm.From the tree row in the west side of a tree belt in the middle of the experimental plot (Figure 1), one single representative tree was selected with a diameter at breast height = 9.05 cm, height = 11.6 m, under-branch height = 3.60 m, and average canopy width = 2.90 m.This tree of average size was selected for leaf number and area, trunk sap flux, and leaf physiological parameters measurements.

Leaf Number and Leaf Area
A steel scaffold (Figure S1) of about 12 m was installed around the sample tree at the end of April 2010 (Figure 1).The canopy of the sample tree at the beginning of the experiment was about 8 m length and was divided into eight layers by 1 m interval (Figure 1).The first-order branches and their sub-branches were also divided into different segments as belonging to different canopy layers according to their spreading trajectories.The different branch segments were distinguished by marking the branches using marker pen.All the first-order branches attached to a specific canopy layer were numbered and marked with waterproof labels, and the heights that they branched from the stem were defined as their branch heights, which was set to 0 at the bottom of the canopy (Figure 1).Detail information on the basal diameters and growth directions of these first-order branches are listed in Table 1.As illustrated in Figure S2, most first-order branches of the sample tree were located in four azimuthal directions of E, NE, N, and NW, which accounted for 76% of the total amount.From early May onward, measurements of leaf number and area were conducted about every 15-20 days.Each time, we climbed on the steel scaffold and counted every leaf within the canopy manually.The leaf number of the different branches and canopy layers were recorded.Then leaf samples, including ca. 50-60 leaves with different sizes, were randomly collected from the lower (canopy layer of 0-3 m), middle (3-6 m), and upper parts (>6 m) of the canopies of the adjacent trees, respectively.These leaf samples were taken to the laboratory for scanning using an Epson Twain Pro high-quality scanner at 400 dpi, and the images were analyzed using automatic image analysis software WinRHIZO (Regent Instruments Inc., Quebec City, QC, Canada) to obtain the projected area of leaves.The average area per leaf (cm 2 leaf −1 ) of the leaf samples from different canopy parts and the leaf number in various canopy layers were used to estimate the total leaf area (LA) of each first-order branch, and within a specific canopy layer of the sample tree.Then, the daily trends in the LAs of the different canopy layers were calculated by linear interpolation using the data of two successive measurements.In addition, in order to investigate the leaf phenology variation with azimuth in different canopy layers, we compared the dynamics of leaf number on first-order branches in east-(narrow row zone, including first-order branches in east, northeast and southeast directions) and west-facing (wide row zone, including first-order branches in west, northwest and southwest directions) sides (Figure 1; Table 1).Leaf number variation rates (number d −1 ) of different first-order branches and canopy layers between two successive measurement times were calculated by dividing the variation of leaf number by the number of days within this period.The leaf number decrease and increase rates of each first-order branch were summed respectively to obtain that of the whole tree, and their summation was the net leaf number variation rate.

Trunk Sap Flux
Trunk sap flux density (cm s −1 ) of the sample tree was measured using thermal dissipation probes (TDP30, Dynamax Inc., Houston, TX, USA).Two sets of probes were inserted into the sapwood about 1.3 m above the ground on the south and north aspects of the sample tree trunk, respectively, in early April 2010.Then, waterproofing sealant was placed around these probes to keep out any water, and reflective bubble insulation was used to wrap the stem and probes to reduce thermal gradients.Sap flux density was measured every 30 s, and the 10-min mean data were stored in an automatic data logger (Model DL2e, Delta-T Devices Inc., Cambridge, UK).During the experiment, sap flux was continuously monitored until the growing season ended.The measured sap flux and the sapwood area (calculated using our established allometric equation describing the variation of sapwood area with the trunk diameter of P. tomentosa [15]) were used to estimate the daily tree transpiration (mm d −1 ) according to Granier [28].Then the crop coefficient (K cb ) of FAO-56 for the sample tree was calculated by dividing the transpiration by the reference crop potential evapotranspiration (mm d −1 ), which was estimated using the meteorological data according to the Penman-Monteith equation [29].

Soil Water Content
At the end of April in 2010, one rectangular experimental plot of 4 m 2 was established within the average growing space around the sample tree, and the plot boundaries were along the middle lines between the sample tree and its adjacent trees in four directions (Figure 1).Twelve access tubes 2 m in length were installed within the experimental plot, and these were located in different distances from the trunk (detail of the installation scheme can be found in Xi et al. [23].A TDR tube probe (TRIME-IPH, IMKO Inc., Ettlingen, Germany) was then utilized to measure the soil volumetric water content (θ, cm 3 cm −3 ) dynamics within the plot from early April to early August with measurement intervals of 1-11 days.The θ profiles were measured every 10 or 20 cm intervals from the soil surface to 170 cm depth.Then the measured soil water data were used to calculate the soil water availability (r θ ) using Equation (1), where θ wp and θ fc are the θ at the soil water potential of wilting point (−1500 kPa) and field capacity (−10 kPa [30]), respectively.

Net Solar Radiation Above and under the Sample Tree Canopy
Two dome net radiometers were used to measure the net solar radiation above (R n-a ) and under (R n-u ) the canopy of the sample tree from the end of April to October in 2010.The radiometers were attached to the steel scaffold around the sample tree, with one about 1 m above the canopy (adjusted the radiometer installation height regularly with the growth of canopy) and the other one 1 m underneath the canopy (Figure 1).The net solar radiation was monitored every 10 min, and the data were stored in an automatic data logger.The ratio between net solar radiation under and above the canopy (R n-u /R n-a ) was calculated to reflect the amount of light intercepted by the canopy, with a higher R n-u /R n-a indicating lower light interception.

Leaf Physiological Parameters and Meteorological Factors
A Li-1600 steady-state porometer (Li-Cor Inc., Lincoln, NE, USA) was used to measure the leaf physiological parameters of transpiration rate (T r ) and stomatal conductance (G s ), plus photosynthetic photon flux density (Q n ) within the canopy.On 21 July (a sunny day) 2010, to measure leaf physiological parameters some three, four and three branches growing in different azimuthal directions were selected from the lower, middle and upper parts of the sample tree canopy, respectively.The measurements were taken every two hours from 08:00 to 18:00, and every time we measured the T r and G s of three fixed leaves from each sample branches and Q n .
Meteorological factors at the experimental site were measured using a weather station (Delta-T Devices Ltd., Cambridge, UK) located about 250 m from the experimental plot.The measured factors included solar radiation (R a ), photosynthetically active radiation, air temperature and relative humidity, wind speed and direction, and precipitation.

Statistical Analysis
Linear regression was applied to analyze the relationship between K cb and the LA in different canopy layers in different periods, with K cb and LA as the dependent and independent variables, respectively.All statistical analyses were undertaken using the Origin 9.0 (OriginLab, Northampton, MA, USA).

Variation of Soil Water Availability
In general, there was no apparent variation trend of r θ with depth throughout the 0-170 cm soil profile within the average growing space around the sample tree.In addition, r θ values in all layers were relatively high, with the seasonal average values varying from 0.73 to 0.92 (Figure 2).During the growing season, r θ varied intensively in 0-30 cm soil and kept relatively stable below 30 cm depth before the end of June.Afterward, with the arrival of the rainy season, r θ in all soil layers increased markedly and then kept stable near 1.0.

Leaf Physiological Parameters and Meteorological Factors
A Li-1600 steady-state porometer (Li-Cor Inc., Lincoln, NE, USA) was used to measure the leaf physiological parameters of transpiration rate (Tr) and stomatal conductance (Gs), plus photosynthetic photon flux density (Qn) within the canopy.On 21 July (a sunny day) 2010, to measure leaf physiological parameters some three, four and three branches growing in different azimuthal directions were selected from the lower, middle and upper parts of the sample tree canopy, respectively.The measurements were taken every two hours from 08:00 to 18:00, and every time we measured the Tr and Gs of three fixed leaves from each sample branches and Qn.
Meteorological factors at the experimental site were measured using a weather station (Delta-T Devices Ltd., Cambridge, UK) located about 250 m from the experimental plot.The measured factors included solar radiation (Ra), photosynthetically active radiation, air temperature and relative humidity, wind speed and direction, and precipitation.

Statistical Analysis
Linear regression was applied to analyze the relationship between Kcb and the LA in different canopy layers in different periods, with Kcb and LA as the dependent and independent variables, respectively.All statistical analyses were undertaken using the Origin 9.0 (OriginLab, Northampton, MA, USA).

Variation of Soil Water Availability
In general, there was no apparent variation trend of rθ with depth throughout the 0-170 cm soil profile within the average growing space around the sample tree.In addition, rθ values in all layers were relatively high, with the seasonal average values varying from 0.73 to 0.92 (Figure 2).During the growing season, rθ varied intensively in 0-30 cm soil and kept relatively stable below 30 cm depth before the end of June.Afterward, with the arrival of the rainy season, rθ in all soil layers increased markedly and then kept stable near 1.0.

Dynamics of Light Interception by the Canopy
As shown in Figure 3, R n-a was lower than R a throughout the growing season, and both of them reached their seasonally highest levels in May, and then declined and stayed relatively stable until mid-September, after which time they declined gradually.Relative to R a and R n-a , R n-u was lower and showed a similar variation trend from May to mid-September, but then it increased distinctly at the end of the growing season.

Dynamics of Light Interception by the Canopy
As shown in Figure 3, Rn-a was lower than Ra throughout the growing season, and both of them reached their seasonally highest levels in May, and then declined and stayed relatively stable until mid-September, after which time they declined gradually.Relative to Ra and Rn-a, Rn-u was lower and showed a similar variation trend from May to mid-September, but then it increased distinctly at the end of the growing season.From early to late May, the Rn-u/Rn-a was relatively stable.Afterward, it declined gradually in the following month and then remained steadily around 0.2 until late August.Thereafter, Rn-u/Rn-a increased sharply to about 0.55 within one month and then stayed stable again until the growing season end.

Leaf Number Variation in Different Branch Heights
There were considerable differences in the seasonality of leaf number on the branches from different branch heights (Figure 4).The leaf number variation of all branches exhibited a single peak, except for that of the branches at the bottom (0-1 m branch height) of the canopy, which had three peaks during the growing season.In the lower part (0-3 m branch height) of the canopy, the maximum leaf number appeared between early May and early June.With increasing branch height, the time when the leaf number increased to its maximum tended to be delayed, i.e., from early June to the middle of July.However, as for the branches at the top of the canopy (>7 m branch height), their leaves also developed relatively early, and maximum leaf number appeared in early June.
Obvious leaf drop on branches in different branch heights was initiated during the period from mid-July to early August, and by late August all leaves at the canopy bottom (0-1 m branch height) had fallen off (Figure 4).This was followed by branches in the 1-3 m branch height range, where all leaves fell off in mid-September.By averaging leaf number through the growing season, branches in 4-5 m branch height of the canopy were shown to possess the most leaves (1365 pieces), accounting for half (51%) of all leaves of the sample tree.Relative to branches in the lower canopy section (1-4 m From early to late May, the R n-u /R n-a was relatively stable.Afterward, it declined gradually in the following month and then remained steadily around 0.2 until late August.Thereafter, R n-u /R n-a increased sharply to about 0.55 within one month and then stayed stable again until the growing season end.

Leaf Number Variation in Different Branch Heights
There were considerable differences in the seasonality of leaf number on the branches from different branch heights (Figure 4).The leaf number variation of all branches exhibited a single peak, except for that of the branches at the bottom (0-1 m branch height) of the canopy, which had three peaks during the growing season.In the lower part (0-3 m branch height) of the canopy, the maximum leaf number appeared between early May and early June.With increasing branch height, the time when the leaf number increased to its maximum tended to be delayed, i.e., from early June to the middle of July.However, as for the branches at the top of the canopy (>7 m branch height), their leaves also developed relatively early, and maximum leaf number appeared in early June.
Obvious leaf drop on branches in different branch heights was initiated during the period from mid-July to early August, and by late August all leaves at the canopy bottom (0-1 m branch height) had fallen off (Figure 4).This was followed by branches in the 1-3 m branch height range, where all leaves fell off in mid-September.By averaging leaf number through the growing season, branches in 4-5 m branch height of the canopy were shown to possess the most leaves (1365 pieces), accounting for half (51%) of all leaves of the sample tree.Relative to branches in the lower canopy section (1-4 m branch height), those in the upper section (>4 m branch height) had much more leaves, and accounted for 77% of the total.

Leaf Number Variation in Different Canopy Layers
Seasonal leaf number variation showed great differences across the various canopy layers (Figure 5).Overall, the date when the leaf number increased to its maximum occurred later in the higher canopy layers.Throughout the canopy, the leaf number in the 6-7 m layer was always the highest during the growing season, and was followed in sequence by that in the 5-6, 4-5 and >7 m (but only after about mid-June) layers.The total season-average leaf numbers in these four layers accounted for 84% of that of the whole canopy.From early May to early June, leaf production rate of 35-54 leaves d −1 increased with time, and was about 3-54 times as great as the leaf falling rate of 1-11 leaves d −1 (Figure 6).Therefore, within

Leaf Number Variation in Different Canopy Layers
Seasonal leaf number variation showed great differences across the various canopy layers (Figure 5).Overall, the date when the leaf number increased to its maximum occurred later in the higher canopy layers.Throughout the canopy, the leaf number in the 6-7 m layer was always the highest during the growing season, and was followed in sequence by that in the 5-6, 4-5 and >7 m (but only after about mid-June) layers.The total season-average leaf numbers in these four layers accounted for 84% of that of the whole canopy.
branch height), those in the upper section (>4 m branch height) had much more leaves, and accounted for 77% of the total.

Leaf Number Variation in Different Canopy Layers
Seasonal leaf number variation showed great differences across the various canopy layers (Figure 5).Overall, the date when the leaf number increased to its maximum occurred later in the higher canopy layers.Throughout the canopy, the leaf number in the 6-7 m layer was always the highest during the growing season, and was followed in sequence by that in the 5-6, 4-5 and >7 m (but only after about mid-June) layers.The total season-average leaf numbers in these four layers accounted for 84% of that of the whole canopy.From early May to early June, leaf production rate of 35-54 leaves d −1 increased with time, and was about 3-54 times as great as the leaf falling rate of 1-11 leaves d −1 (Figure 6).Therefore, within From early May to early June, leaf production rate of 35-54 leaves d −1 increased with time, and was about 3-54 times as great as the leaf falling rate of 1-11 leaves d −1 (Figure 6).Therefore, within this period, the leaf number of the whole canopy increased greatly with time and reached its annual maximum in early June (Figure 6).Afterward, the rates of leaf production (4-9 leaves d −1 ) and leaf drop (5-13 leaves d −1 ) were relatively small and were almost identical until mid-July.So there was only little variation in the leaf number in this period.After mid-July, the sample tree appeared to stop producing leaves, while its rate of leaf drop increased abruptly and accelerated with time.Consequently, the number of leaves decreased rapidly during this period, with about 59 leaves falling off every day (Figure 6).
this period, the leaf number of the whole canopy increased greatly with time and reached its annual maximum in early June (Figure 6).Afterward, the rates of leaf production (4-9 leaves d −1 ) and leaf drop (5-13 leaves d −1 ) were relatively small and were almost identical until mid-July.So there was only little variation in the leaf number in this period.After mid-July, the sample tree appeared to stop producing leaves, while its rate of leaf drop increased abruptly and accelerated with time.Consequently, the number of leaves decreased rapidly during this period, with about 59 leaves falling off every day (Figure 6).

Leaf Number Variation in Different Azimuthal Sides
As illustrated in Figure 7, in both the lower (0-3 m branch height) and middle (3-6 m) parts of the canopy, the maximum leaf number of branches in the east-facing sides appeared about half a month earlier than that in the west-facing sides.In contrast, in the upper canopy (>6 m branch height), the leaves on branches on both sides increased to their maximum at the same time (early July).Overall, as to the whole canopy, the time when leaves on the branches increased to the maximum number was about twenty days earlier on the east-facing side.

Leaf Number Variation in Different Azimuthal Sides
As illustrated in Figure 7, in both the lower (0-3 m branch height) and middle (3-6 m) parts of the canopy, the maximum leaf number of branches in the east-facing sides appeared about half a month earlier than that in the west-facing sides.In contrast, in the upper canopy (>6 m branch height), the leaves on branches on both sides increased to their maximum at the same time (early July).Overall, as to the whole canopy, the time when leaves on the branches increased to the maximum number was about twenty days earlier on the east-facing side.
this period, the leaf number of the whole canopy increased greatly with time and reached its annual maximum in early June (Figure 6).Afterward, the rates of leaf production (4-9 leaves d −1 ) and leaf drop (5-13 leaves d −1 ) were relatively small and were almost identical until mid-July.So there was only little variation in the leaf number in this period.After mid-July, the sample tree appeared to stop producing leaves, while its rate of leaf drop increased abruptly and accelerated with time.Consequently, the number of leaves decreased rapidly during this period, with about 59 leaves falling off every day (Figure 6).

Leaf Number Variation in Different Azimuthal Sides
As illustrated in Figure 7, in both the lower (0-3 m branch height) and middle (3-6 m) parts of the canopy, the maximum leaf number of branches in the east-facing sides appeared about half a month earlier than that in the west-facing sides.In contrast, in the upper canopy (>6 m branch height), the leaves on branches on both sides increased to their maximum at the same time (early July).Overall, as to the whole canopy, the time when leaves on the branches increased to the maximum number was about twenty days earlier on the east-facing side.

Leaf Area Dynamics and Vertical Distribution
During the growing season, the LA variation amplitude in the 0-3 m canopy height was very small, but increased markedly with height above 3 m until reaching the canopy top, where it decreased greatly (Figure 8).In canopy above the 3 m height, the LA varied relatively little between 9-24 May, but from 24 May to 10 June it increased sharply by 64% to 265% in the different canopy layers.Then during the following one month (10 June-15 July), the LA in the different layers varied little and kept their highest levels for the growing season.After about 15 July, the LA in almost all layers started to decrease rapidly until the growing season ended.

Figure 7.
Variation in leaf number on first-order branches in the east-(narrow row zone) and westfacing (wide row zone) sides of the (a) upper (0-3 m branch height (Figure 1)), (b) middle (3-6 m) and (c) lower canopy parts (>6 m), and of (d) the whole canopy.

Leaf Area Dynamics and Vertical Distribution
During the growing season, the LA variation amplitude in the 0-3 m canopy height was very small, but increased markedly with height above 3 m until reaching the canopy top, where it decreased greatly (Figure 8).In canopy above the 3 m height, the LA varied relatively little between 9-24 May, but from 24 May to 10 June it increased sharply by 64% to 265% in the different canopy layers.Then during the following one month (10 June-15 July), the LA in the different layers varied little and kept their highest levels for the growing season.After about 15 July, the LA in almost all layers started to decrease rapidly until the growing season ended.In the lower part of the canopy (0-3 m canopy height), average LAs were similar in the different layers and were much lower than that in the other canopy layers.Above 3 m canopy height, the LA increased greatly and almost linearly with height.However, at the canopy top above 7 m height, the LA decreased sharply to a level similar to that in the middle canopy layer (Figure 8).
The seasonally averaged total LA of the sample tree was 15.86 m 2 (1.77-24.80m 2 ).Among different canopy layers, the LA within 6-7 m canopy height was the highest, accounting for one third (35%) of the total tree's LA.In addition, on average, 88% (75%-95%) of the total LA of the sample tree was allocated to its upper canopy sections (>4 m canopy height) (Figure 8).

Leaf Physiological Parameter Variation within the Canopy
On a typically sunny day on 21 July, when total leaf number and area were still at relatively high levels (Figures 5 and 8), Gs, Tr and Qn in both the lower and middle parts of the canopy exhibited a unimodal diurnal variation, with their peak values appeared between 12:00 and 14:00 (Figure 9).As for the upper canopy, both Gs and Qn showed bimodal variations, with their peak values occurring at both 10:00 and 14:00, but Tr only increased to its maximum at 14:00.Overall, with increasing canopy height, Gs, Tr and Qn increased markedly, especially between 10:00-15:30.Daily average values of Gs In the lower part of the canopy (0-3 m canopy height), average LAs were similar in the different layers and were much lower than that in the other canopy layers.Above 3 m canopy height, the LA increased greatly and almost linearly with height.However, at the canopy top above 7 m height, the LA decreased sharply to a level similar to that in the middle canopy layer (Figure 8).
The seasonally averaged total LA of the sample tree was 15.86 m 2 (1.77-24.80m 2 ).Among different canopy layers, the LA within 6-7 m canopy height was the highest, accounting for one third (35%) of the total tree's LA.In addition, on average, 88% (75%-95%) of the total LA of the sample tree was allocated to its upper canopy sections (>4 m canopy height) (Figure 8).

Leaf Physiological Parameter Variation within the Canopy
On a typically sunny day on 21 July, when total leaf number and area were still at relatively high levels (Figures 5 and 8), G s , T r and Q n in both the lower and middle parts of the canopy exhibited a unimodal diurnal variation, with their peak values appeared between 12:00 and 14:00 (Figure 9).As for the upper canopy, both G s and Q n showed bimodal variations, with their peak values occurring at both 10:00 and 14:00, but T r only increased to its maximum at 14:00.Overall, with increasing canopy height, G s , T r and Q n increased markedly, especially between 10:00-15:30.Daily average values of G s and T r in the upper canopy were 0.52 cm s −1 and 7.63 µg cm −2 s −1 , respectively, which were 17% and 87% (G s ), and 16% and 83% (T r ) higher than that in the middle and lower canopies, respectively.
and Tr in the upper canopy were 0.52 cm s −1 and 7.63 μg cm −2 s −1 , respectively, which were 17% and 87% (Gs), and 16% and 83% (Tr) higher than that in the middle and lower canopies, respectively.

Response of Transpiration to Leaf Areas in Different Canopy Layers
As shown in Figure 10, between 9-26 May, the Kcb of the sample tree reached its highest seasonal level with a mean of 0.71.Within this period, the LA in the 0-3 m canopy layer also reached its maximum, but in the layers above 3 m canopy height the LAs were still very low relative to their corresponding maximum values.Linear regression analysis showed non-significant correlated relationships existed between Kcb and the LAs in all canopy layers (p > 0.05) in this period.
From 27 May to 23 August, Kcb increased linearly (p < 0.001) with LA in canopy layers above 3 m height, although the increment was slight.Whereas, variation of LA in the 0-3 m canopy layer showed a non-significant correlation with that of Kcb (p > 0.05).In addition, the proportion of the Kcb variation during this period explained by LA was the highest in the 5-7 m canopy layers (R 2 = 0.300), which was followed in sequence by the 4-5 (R 2 = 0.257), 3-4 (R 2 = 0.235) and > 7 m (R 2 = 0.139) canopy layers.After late August, the Kcb and LA in all canopy layers (except the 0-1 m layer) decreased rapidly, and significant linear relationships were detected between them (p < 0.001).

Response of Transpiration to Leaf Areas in Different Canopy Layers
As shown in Figure 10, between 9-26 May, the K cb of the sample tree reached its highest seasonal level with a mean of 0.71.Within this period, the LA in the 0-3 m canopy layer also reached its maximum, but in the layers above 3 m canopy height the LAs were still very low relative to their corresponding maximum values.Linear regression analysis showed non-significant correlated relationships existed between K cb and the LAs in all canopy layers (p > 0.05) in this period.
From 27 May to 23 August, K cb increased linearly (p < 0.001) with LA in canopy layers above 3 m height, although the increment was slight.Whereas, variation of LA in the 0-3 m canopy layer showed a non-significant correlation with that of K cb (p > 0.05).In addition, the proportion of the K cb variation during this period explained by LA was the highest in the 5-7 m canopy layers (R 2 = 0.300), which was followed in sequence by the 4-5 (R 2 = 0.257), 3-4 (R 2 = 0.235) and > 7 m (R 2 = 0.139) canopy layers.After late August, the K cb and LA in all canopy layers (except the 0-1 m layer) decreased rapidly, and significant linear relationships were detected between them (p < 0.001).
and Tr in the upper canopy were 0.52 cm s −1 and 7.63 μg cm −2 s −1 , respectively, which were 17% and 87% (Gs), and 16% and 83% (Tr) higher than that in the middle and lower canopies, respectively.

Response of Transpiration to Leaf Areas in Different Canopy Layers
As shown in Figure 10, between 9-26 May, the Kcb of the sample tree reached its highest seasonal level with a mean of 0.71.Within this period, the LA in the 0-3 m canopy layer also reached its maximum, but in the layers above 3 m canopy height the LAs were still very low relative to their corresponding maximum values.Linear regression analysis showed non-significant correlated relationships existed between Kcb and the LAs in all canopy layers (p > 0.05) in this period.
From 27 May to 23 August, Kcb increased linearly (p < 0.001) with LA in canopy layers above 3 m height, although the increment was slight.Whereas, variation of LA in the 0-3 m canopy layer showed a non-significant correlation with that of Kcb (p > 0.05).In addition, the proportion of the Kcb variation during this period explained by LA was the highest in the 5-7 m canopy layers (R 2 = 0.300), which was followed in sequence by the 4-5 (R 2 = 0.257), 3-4 (R 2 = 0.235) and > 7 m (R 2 = 0.139) canopy layers.After late August, the Kcb and LA in all canopy layers (except the 0-1 m layer) decreased rapidly, and significant linear relationships were detected between them (p < 0.001).

Leaf Phenology Variation within the Canopy
To fully utilize the light resources so as to maximize growth within the growing season, plants usually develop different leaf phenological patterns within the canopy [9,11,16,17,20].This phenomenon was well demonstrated in our study by laboriously counting manually all the leaves of a mature P. tomentosa tree through the growing season.It was found that, except branches at the canopy top (>7 m branch height), the maximum leaf number tended to be reached later on branches attached to higher stems (Figure 4).A similar trend was more apparently reflected in the variation of leaf number at different canopy layers with canopy height, and the maximum leaf number in the top layer of canopy occurred last (Figure 5).This phenological asynchrony within the canopy indicates that P. tomentosa has adopted a "temporal ecological strategy" for efficient light acquisition [7].That is to say, earlier development and maturation of lower leaves allowed the lower canopy of P. tomentosa to fully intercept and utilize the higher irradiance in spring (Figure 3).However, it seems that P. tomentosa did not rely on its lower leaves to utilize the markedly increased irradiance in the lower canopy in autumn (Figure 3), as the time when all leaves had fallen off was earlier in the lower canopy.Similar leaf phenology variation within the canopy was also observed on many other tree species, such as Lyonia ovalifolia [11] and Acer saccharum Marshall [17] in temperate forests.
Although many researchers have investigated the leaf phenology variation with canopy height [9,12,[16][17][18]20], little attention has been paid to the variability in leaf phenology at different azimuths within different canopy layers.Previous studies showed that leaf phenology can be influenced by many abiotic and biotic factors, such as light, air temperature and humidity, wind speed, leaf water status, and ontogeny [8, 16,[31][32][33].However, as to the leaf phenology variation within the canopy, the light environment and leaf water status have been recognized as the two most critical determinants [16,31].In our plantation, due to the non-uniform characteristics of tree planting (Figure 1), there were great differences in the growing spaces of canopies at different azimuths, which could subsequently result in azimuthal variation in the light environment within the canopy.This might contribute to the earlier achievement of maximum leaf number in the east-facing sides (narrow row zone) of both the lower and middle canopies, relative to that in the west-facing sides (wide row zone) (Figure 7).As leaf water potential might have varied relatively little within the same canopy layer in a dense plantation, leaf water status difference might not be the reason for this phenomenon.On the contrary, maximum leaf number in both the east-and west-facing sides occurred simultaneously in the upper canopy (>6 m branch height) (Figure 7).This might be due to the relatively small difference in the irradiation among different azimuthal sides due to all sides were directly exposed to sunlight.Our observed phenomenon implies that the "temporal ecological strategy" also seems to be adopted by leaves in different azimuths of the same canopy layer for efficient light acquisition.In addition, the coexistence of leaf phenology synchronization, and asynchronization in different azimuthal sides within the canopy indicates that it will be biased to reveal the leaf phenology of a tree species growing under non-uniform planting scheme by just observing leaf development in a specific part of the canopy.
The present research is the first study to investigate thoroughly and precisely the leaf phenological variation of P. tomentosa within the canopy.Due to the great difficulty in observing leaf dynamics of big trees in a dense plantation, only one sample tree with average size was selected in this study.Thus, the leaf phenology observed on this single sample tree may not be completely applicable to other P. tomentosa trees.Nonetheless, we think our results are very important and valuable since, to our knowledge, no such attempts have been conducted to observe the dynamics of all leaves of a single big tree in the field.Although the data before early May were absent, our results indicated that the net leaf number variation rate of mature P. tomentosa exhibited a unimodal seasonal pattern with the peak occurring between late May and early June.This was different from the bimodal seasonal variation pattern, with peaks occurred in May and July, in tree diameter in the same plantation [5].
The main leaf production period of mature P. tomentosa under plantation conditions was before mid-June and its main leaf drop was initiated at about mid-July, which was roughly consistent with the leaf area index dynamics in the present plantation [15].This information will be useful for the cultivation of P. tomentosa plantations across similar sites to ours and for the understanding of its ecophysiological activities.

Leaf Area Vertical Distribution and Its Relationship with Canopy Light Interception
The LA in different canopy layers of the sample tree, and its variation, were strongly influenced by height.In the lower canopy parts (0-3 m canopy height), the LA in different layers was similar and varied relatively little.Whereas, above the 3 m canopy height, both the LA and its variation amplitude increased markedly with height, but decreased abruptly at the top of the canopy.This kind of vertical profile in LA might be determined by the decreasing light availability from canopy top to bottom (Figure 9), as the average widths of most canopy layers were similar (2-3 m) and showed no distinct variation trend with height (data not shown).This vertical LA distribution pattern can also help to explain the phenomenon of decreasing leaf water status with increasing canopy height, which had been observed in the present plantation by [34], and in other forests [10,17,35].This is because the amount of water transferred by sapwood would be shared by more leaves in the higher canopy layers.
There was an evident seasonality in the canopy light interception of the sample tree, reflected by R n-u /R n-a (Figure 3).This might be mainly controlled by the leaf development in the upper canopy because of its greater variability in the growing season (Figures 5 and 8).Thus, from about mid-June to mid-August, when LA in the upper canopy had reached its seasonal highest level, the amount of light intercepted by canopy also increased to its seasonal maximum.This subsequently resulted in R n-u decreasing to its annual lowest level.Overall, the average R a and R n-a were higher in spring than in autumn, but R n-u in these two seasons kept within similar ranges, which could be mainly due to the higher canopy light interception ratio in spring (Figure 3).This information on the seasonal R n-u and R n-a variation will help to understand and predict the seasonal dynamics of evapotranspiration and its components in dense plantations like ours [13,15].
Overall, relative to the situation in summer, both the LA of the sample tree and its vertical gradient within the canopy were lower in spring and autumn (Figure 8).According to the Lambert-Beer law, this might result in a smaller gradient in light availability from the canopy top to the bottom and higher light availability in lower canopy in spring and autumn.This was also reflected by the higher R n-u in these two seasons (Figure 3).As leaf transpiration rate and stomatal conductance were usually positively correlated with the surrounding light availability [10,11,36], the difference in leaf transpiration rates between the lower and upper canopies of trees in the present plantation might be smaller in spring and autumn.This speculated phenomenon has also been observed, or reflected in the data in other studies [11,36].Consequently, due to the seasonal variations of both the leaf transpiration rate and the LA along the canopy profile, the contribution of leaves in the lower canopy to the whole-tree transpiration might be higher in spring and autumn but lower in summer.

The Relationship between Transpiration and Leaf Area
To detect precisely the relationship between tree transpiration and leaf area, it is essential to eliminate the effects of soil water availability and other environmental factors [13,15].These requirements were achieved in this study by guaranteeing the high r θ in the rootzone around the sample tree (Figure 2) and by dividing its transpiration rate by the reference evapotranspiration to obtain K cb of the sample tree [29].
The K cb is a critical parameter for predicting plant water use on different scales and can provide insight into plant water relations in different environments [29].Thus, in order to obtain precise information on the seasonal K cb variation, many studies have related K cb to some easily measured plant development variables, such as leaf area, leaf area index, or canopy cover [13,25,37].Some researchers found that the relationship between K cb and plant variables was linear [13,37], while others, and our previous study in the same plantation, found it was nonlinear.The K cb tended to plateau as plant variables increased above a critical value [13,15,25,38,39].These contrasting conclusions were mainly obtained at the stand scale and based on the development variables of the whole canopies of plant individuals or populations.Whereas, in the present study, it was directly demonstrated that, on a single tree scale in a dense plantation, the K cb increased nonlinearly with increasing leaf areas of not only the whole tree (data not shown), but also different canopy layers.
In our previous study [15], the phenomenon that K cb of P. tomentosa plantation tended to plateau as the stand LAI increased above a threshold value was attributed to the fact that the stand sun-exposed LA in the upper canopy, which determines the plant water consumption [37], reached its maximum when the stand LAI was still small.However, based on the results of this study, it is more plausible to attribute that phenomenon to the seasonal variation in the vertical distribution of both the LA, and the sun-exposed LA within the canopy.In spring, as the total stand LA [15] and canopy light interception were small (Figure 3), a large portion of leaves from upper canopy to the bottom could be exposed directly to sunlight.That is, sun-exposed leaves were distributed throughout the whole canopy profile.This, combined with the high irradiation, high wind speed and low relative humidity [15] might result in the K cb increasing to its annual maximum level in late spring (May) (Figure 10), even though average tree LA was still relatively small at this time (Figure 8).After late May, both the stand LA (especially the LA of the upper canopy) and canopy light interception increased markedly while light in the lower canopy decreased greatly, which could subsequently lead to that sun-exposed leaves were mainly distributed on canopy top.As the LA of trees and its vertical distribution pattern changed little from June to late August, the sun-exposed LA, and dependent K cb , exhibited only a small variation during this period (Figure 10).Consequently, as for the whole growing season, the stand [15] and single tree (Figure 10) K cb varied relatively little with the LA after they both achieved their seasonal highest levels in late May or early June.

Contribution of Leaves in Different Canopy Layers to Tree Transpiration
Across different periods during the growing season, the relationships between K cb of the sample tree and LA of its different canopy layers changed greatly (Figure 10).This implies that the relative contributions of leaves in different canopy layers to the whole tree transpiration varied through the growing season.In May, the seasonal maximum levels of K cb and LA in the 0-3 m canopy layer occurred concurrently, while the LA in the upper layers stayed at relatively low levels (Figures 8 and 10), which strongly supports our speculation that the leaves in the lower canopy contributed greatly to the whole tree transpiration in spring.From June to late August, the non-significant correlation between LA and K cb indicated that the leaves in the bottom 3 m canopy layer only made a small contribution to tree water use.On the contrary, the leaves in layers above the 3 m canopy height exhibited an obvious contribution to the transpiration of the tree during this period, which was reflected by their significantly linear correlations (p < 0.001).In addition, according to the amount of variation in the K cb explained by the LA in different layers above 3 m, the leaves distributed within 5-7 m canopy height exerted the greatest control on the tree water use, which could be due to the coupling effects of both larger LA and higher leaf transpiration rates and stomatal conductances (Figure 9).After late August, due to the distinctly decreasing canopy light interception and the upper canopy LA (Figures 3 and 8), a large number of leaves in the lower canopy were directly exposed to sunlight, resulting in sun-exposed leaves that were distributed throughout the whole canopy profile again.Therefore, in this period, the leaves in all the canopy layers (except the bottom 1 m layer where all leaves had fallen off) showed a significant linear correlation with K cb (p < 0.001).That is, all of them contributed to transpiration.
Consequently, it can be deduced that the stand transpiration in a dense plantation may be mainly controlled by the total area and spatial distribution of the sun-exposed leaves within the canopy, and all the different canopy layers will exhibit different controls on the stand water relations at different periods.This information will provide a theoretical basis for further optimizing stand water-use

Figure 1 .
Figure 1.Schematic diagrams of the alternate wide-and narrow-row spacing planting scheme, the relative position of the sample tree, the installation patterns for dome net radiometers above and under the sample tree canopy, and the definitions of canopy layer and branch height used in this study.

Figure 1 .
Figure 1.Schematic diagrams of the alternate wide-and narrow-row spacing planting scheme, the relative position of the sample tree, the installation patterns for dome net radiometers above and under the sample tree canopy, and the definitions of canopy layer and branch height used in this study.

Figure 2 .
Figure 2. Seasonal variation of soil water availability (r θ ) at different soil depths within the average growing space of the sample tree.

Forests 17 Figure 2 .
Figure 2. Seasonal variation of soil water availability (rθ) at different soil depths within the average growing space of the sample tree.

Figure 3 .
Figure 3. Seasonal variation of solar radiation (Ra) in the experimental site, net solar radiation above (Rn-a) and under (Rn-u) the sample tree canopy, and the ratio between Rn-u and Rn-a.

Figure 3 .
Figure 3. Seasonal variation of solar radiation (R a ) in the experimental site, net solar radiation above (R n-a ) and under (R n-u ) the sample tree canopy, and the ratio between R n-u and R n-a.

ForestsForests
2017, 8, x FOR PEER REVIEW 8 of 17 branch height), those in the upper section (>4 m branch height) had much more leaves, and accounted for 77% of the total.

Figure 4 .
Figure 4. Variation of the leaf number on branches from different branch heights (Figure 1) of the sample tree.

Figure 5 .
Figure 5. Variation of the leaf number in different canopy layers and of the total canopy.

Figure 4 .
Figure 4. Variation of the leaf number on branches from different branch heights (Figure 1) of the sample tree.

Figure 4 .
Figure 4. Variation of the leaf number on branches from different branch heights (Figure 1) of the sample tree.

Figure 5 .
Figure 5. Variation of the leaf number in different canopy layers and of the total canopy.

Figure 5 .
Figure 5. Variation of the leaf number in different canopy layers and of the total canopy.

Figure 6 .
Figure 6.Leaf number decrease, increase and net variation rates of the sample tree in different periods during the growing season in 2010.The date format in different periods is "month.day".

Figure 6 .
Figure 6.Leaf number decrease, increase and net variation rates of the sample tree in different periods during the growing season in 2010.The date format in different periods is "month.day".

Figure 6 .
Figure 6.Leaf number decrease, increase and net variation rates of the sample tree in different periods during the growing season in 2010.The date format in different periods is "month.day".

Figure 7 .
Figure 7. Variation in leaf number on first-order branches in the east-(narrow row zone) and west-facing (wide row zone) sides of the (a) upper (0-3 m branch height (Figure 1)), (b) middle (3-6 m) and (c) lower canopy parts (>6 m), and of (d) the whole canopy.

Figure 8 .
Figure 8. Variation of the leaf areas with canopy height of the sample tree at different dates.The inset shows the cumulative leaf areas distribution plotted against the canopy height of the sample tree at different dates.The date format at different times is "year/month/day".

Figure 8 .
Figure 8. Variation of the leaf areas with canopy height of the sample tree at different dates.The inset shows the cumulative leaf areas distribution plotted against the canopy height of the sample tree at different dates.The date format at different times is "year/month/day".

Figure 9 .
Figure 9.Diurnal variation of (a) the leaf stomatal conductance (Gs), (b) leaf transpiration rate (Tr) and (c) photosynthetic photon flux density (Qn) in the lower, middle and upper canopy parts of the sample tree on 21 July 2010.

Figure 9 .
Figure 9.Diurnal variation of (a) the leaf stomatal conductance (G s ), (b) leaf transpiration rate (T r ) and (c) photosynthetic photon flux density (Q n ) in the lower, middle and upper canopy parts of the sample tree on 21 July 2010.

Figure 9 .
Figure 9.Diurnal variation of (a) the leaf stomatal conductance (Gs), (b) leaf transpiration rate (Tr) and (c) photosynthetic photon flux density (Qn) in the lower, middle and upper canopy parts of the sample tree on 21 July 2010.

Figure 10 .
Figure 10.Variation of the crop coefficient (K cb ) of the sample tree with leaf areas (LA) in different canopy layers.The date format in different periods is "month.day".

Table 1 .
Basal diameters and growth directions of first-order branches belonging to different branch height ranges.