Biophysical Controls on Light Response of Net CO2 Exchange in a Winter Wheat Field in the North China Plain

To investigate the impacts of biophysical factors on light response of net ecosystem exchange (NEE), CO2 flux was measured using the eddy covariance technique in a winter wheat field in the North China Plain from 2003 to 2006. A rectangular hyperbolic function was used to describe NEE light response. Maximum photosynthetic capacity (P max) was 46.6±4.0 µmol CO2 m−2 s−1 and initial light use efficiency (α) 0.059±0.006 µmol µmol−1 in April−May, two or three times as high as those in March. Stepwise multiple linear regressions showed that P max increased with the increase in leaf area index (LAI), canopy conductance (g c) and air temperature (T a) but declined with increasing vapor pressure deficit (VPD) (P<0.001). The factors influencing P max were sorted as LAI, g c, T a and VPD. α was proportional to ln(LAI), g c, T a and VPD (P<0.001). The effects of LAI, g c and T a on α were larger than that of VPD. When T a>25°C or VPD>1.1−1.3 kPa, NEE residual increased with the increase in T a and VPD (P<0.001), indicating that temperature and water stress occurred. When g c was more than 14 mm s−1 in March and May and 26 mm s−1 in April, the NEE residuals decline disappeared, or even turned into an increase in g c (P<0.01), implying shifts from stomatal limitation to non-stomatal limitation on NEE. Although the differences between sunny and cloudy sky conditions were unremarkable for light response parameters, simulated net CO2 uptake under the same radiation intensity averaged 18% higher in cloudy days than in sunny days during the year 2003−2006. It is necessary to include these effects in relevant carbon cycle models to improve our estimation of carbon balance at regional and global scales.


Introduction
Vegetation productivity is the foundation of carbon sequestration and grain yield formation. For the limit on field observation techniques, early studies mainly focused on photosynthesis at the leaf level. The biochemistry and ecophysiology of leaf photosynthesis have been well understood and parameterized [1,2]. However, it is difficult to measure leaf photosynthesis over long periods and upscale photosynthetic rate from the leaf level to the canopy or ecosystem level because of the non-linear distribution of leaf area and radiation intensity within the vegetation canopy. With the development of micrometeorological techniques, especially the eddy covariance (EC) method, net photosynthetic rate of vegetation could be directly measured at the ecosystem level. The EC technique has been widely used in CO 2 flux measurements in forest, grassland and farmland ecosystems [3][4][5][6][7][8][9].
Plant photosynthesis is primarily driven by incident solar radiation. Light response models, including rectangular hyperbolic model [10,11] and non-rectangular hyperbolic model [12,13], were developed to describe the relationship between daytime net ecosystem CO 2 exchange (NEE) and photosynthetically active radiation (PAR). Besides solar radiation, factors influencing daytime NEE include environmental variables such as air temperature, vapor pressure deficit and soil water content and biological variables such as leaf area index and canopy conductance. These factors may be considered by (1) estimating the bias of simulated NEE as functions of influencing factors [14,15], or (2) revising light response parameters as functions of influencing factors [9,11,[16][17][18][19]. The effects of biophysical factors on light response parameters are often estimated using stepwise multiple linear regression models [9]. However, these models have not been used to assess the influence of variables on NEE residual so far.
Long-term observations and simulations have shown a worldwide decrease in surface solar radiation (global dimming) from 1950s to 1980s, with a partial recovery (brightening) in 1990s in some areas (e.g. high and middle latitude of the Northern Hemisphere) [20]. During the dimming period, direct radiation declined remarkably, whereas diffuse radiation enhanced. The increase in diffuse fraction of surface solar radiation may be attributed to increasing aerosol and/or cloudiness, as both factors tend to enhance scattering in the atmosphere [20]. Light responses of photosynthesis under various sky conditions were investigated in forests [8,[21][22][23] and croplands [24,25]. For the tall and dense vegetations, photosynthesis may be enhanced by diffuse radiation because diffuse light distributes more effectively within the canopy compared with direct light [26]. After considering the effects of direct and diffuse radiation on canopy photosynthesis in the model, Mercado et al. [27] found that the increases in diffuse fraction enhanced the global land carbon sink by 23.7% during the period from 1960 to 1999.
China is the largest wheat producer and consumer in the world [28] and continuously attempts to increase its production to ensure national food security. As one of large food production regions in China, the North China Plain (NCP) produces about half of the country's wheat [29]. In this study, CO 2 flux was measured continuously using the EC technique in a winter wheat field in the North China Plain for 4 years. The objectives are to (1) investigate NEE light response and the influencing factors, and (2) assess the effects of sky conditions on NEE light response. This study will improve our knowledge on the parameterization of carbon cycle models and the scenario analyses of carbon sink in the future under changing climate.

Study site
This study was conducted at Yucheng Comprehensive Experiment Station, Chinese Academy of Sciences (36u579N, 116u389E, 23.4 m). It is located at the North China Plain, with a temperate monsoon climate. Mean annual temperature is 13.1uC and annual solar radiation is 5242 MJ m 22 . Annual precipitation is about 528 mm. Soil organic content is 1.21% and pH value is about 7.9. The typical cropping system in this region is the biannual rotation with winter wheat and summer maize. In this study, winter wheat was planted in mid/late October and harvested in early/mid June. The detailed field management was described by Tong et al. [30].
Field observations CO 2 and latent heat fluxes were measured by the eddy covariance system with a three-dimensional sonic anemometer (model CSAT3, Campbell Sci. Inc., USA) and an infrared openpath CO 2 /H 2 O gas analyzer (model LI-7500, Li-Cor Inc., USA). The eddy covariance system, mounted at the height of 2.1 m, was used to measure 3-D wind speed, air temperature, humidity and CO 2 concentration above the canopy. Raw data were collected at 10 Hz and recorded by a CR5000 datalogger (model CR5000, Campbell Scientific Inc., USA). The CR5000 datalogger can store data and programs on a PC card. The card can be carried to a computer and the computer reads PC cards via its PCMCIA card slot directly.
Anemometers (model A100R, Vector, UK) and psychrometers (model HMP45C, Vaisala, Finland) were installed at heights of 2.05 and 3.25 m above the ground. Photosynthetically active radiation (PAR) was measured using a quantum sensor (model Li-190SB, Li-Cor Inc, USA). Solar radiation and net radiation (R n ) were measure by a pyranometer (model CM11, Kipp & Zonen, Delft, The Netherlands) and a net radiometer (model CNR-1, Kipp & Zonen, Delft, The Netherlands), respectively. Two soil heat flux plates were buried in the depth of 2 cm, one between rows and another between plants. Soil temperature was measured at 0, 5, 10, 30 and 50 cm depths. Soil water content (SWC) at the depths of 20 cm and 30 cm was measured with time domain reflectometers (TDR) (model CS616, Campbell Sci. Inc., USA). Rainfall was measured with a rain gauge (model 52203, Rm Young MI, USA). All meteorological data were recorded with a data logger (model CR23x, Campbell Sci. Inc., USA) and were stored at intervals of 30 min.
Biomass, leaf area index (LAI) and plant height were measured every 5 days during the growing season of winter wheat. There were three sampling plots (replications) for each measurement. Twenty plants were sampled continuously for each plot and the plots were selected randomly. Leaf area was measured with Leaf area instrument (model Li-3100, Campbell Sci. Inc., USA). In this study, the main growing season began when the cropland turned from the carbon source to the sink (5-day moving mean CO 2 flux from positive to negative), and it ended when shifted reversely.
Flux data quality control CO 2 and latent heat fluxes were calculated as follows: where r is air density, w' the vertical wind velocity, c' CO 2 concentration, l the latent heat of vaporization, E water vapor flux and q the specific humidity. Overbars indicate an averaging operation and primes denote deviations from the mean. Raw data were conducted by two dimension coordinate rotations [31] and Webb-Pearman-Leuning (WPL) correction [32] to obtain 30-min mean flux data. CO 2 and water vapor fluxes could be affected by rain and dew. The abnormal data were eliminated following the method used by Falge et al. [33]. In the daytime, data gaps were 29%, 8%, 10% and 16% in the growing seasons of winter wheat in 2003, 2004, 2005 and 2006, respectively. Gap filled data were not used in this paper.

Light response of NEE
In most ecosystems, the relationship between daytime NEE and PAR can be expressed by a rectangular hyperbolic function [33,34]: where a is the initial slope of the light response curve (initial light use efficiency), P max the maximum photosynthetic capacity, R d the daytime ecosystem respiration rate under dark conditions. The model (Eq. (3)) was fitted by the software ''Origin 7.0'' (Microcal Software Inc.). In this study, negative NEE means net CO 2 uptake by the cropland and positive NEE indicates net CO 2 emission from the cropland.
To study the environmental factors influencing daytime NEE besides PAR, NEE was computed using Eq. (3) and the residual (NEE r ) was obtained as measured NEE minus simulated NEE.
Stepwise multiple linear regression models were used to estimate the integrative influences of biophysical factors on NEE r . Statistical analyses were performed using the software ''SPSS 13.0'' (SPSS Inc.). According to the method used by Carrara et al. [17], Powell et al. [18] and Teklemariam et al. [7], air temperature (T a ), vapor pressure deficit (VPD), soil water content (SWC), leaf area index (LAI) and canopy conductance (g c ) were divided into many classes and mean NEE r was obtained for each class to show the trends in NEE r versus variables.
To investigate seasonal patterns of light response parameters, a rectangular hyperbola (Eq. (3)) was used to describe the light response of NEE every five days and a time series of a, P max and R d were obtained. Meanwhile, the biophysical variables in the daytime with the 5-day interval were averaged. Stepwise multiple linear regression analyses were applied to distinguish the key drivers for a, P max and R d . The relationships among light response parameters were also investigated at the 5-day scale.
Light response parameters varied among years. Standard errors (e 1 and e 2 ) were calculated for sunny and cloudy sky conditions and total standard error (e) was obtained for both sky conditions using the equation as follows: The difference between sunny and cloudy sky conditions would be pronounced if it was larger than the total standard error for two sky conditions.

Canopy conductance
According to Monteith and Unsworth [35], canopy conductance (g c ) was calculated as follows: where lE is the latent heat flux (W m 22 ), R n net radiation (W m 22 ), G soil heat flux (W m 22 ), g the slope of saturation vapor pressure to the air temperature curve (kPa K 21 ), r air density (mol m 23 ), C p the specific heat capacity of air (J mol 21 K 21 ), VPD vapor pressure deficit (kPa), c the psychrometric constant (kPa K 21 ), g c canopy conductance (m s 21 ), and g a aerodynamic conductance (m s 21 ) which can be given as follows [36]: where ra is aerodynamic resistance (s m21), z the height where wind velocity measured (m), d zero-plane displacement height (m), z0 roughness length (m), k the von Karman constant (0.41) and u wind velocity (m s21). z0 and d can be calculated as [37,38]: where h is the crop height (m).

Defining sunny and cloudy sky conditions
The clearness index (k t ) is used to describe sky conditions and the degree of impact of cloudiness on the solar radiation received at the Earth's surface [22]. It can be given by [39]: where S is global solar radiation (W m 22 ), S e the extraterrestrial irradiance at a plane parallel to the earth surface (W m 22 ), S sc the solar constant (1370 W m 22 ), t d the day of year, b the solar elevation angle, w the local latitude, d the declination of the sun, and v hour angle. The sky conditions were classified at a half-day scale because the days without clouds for the whole daytime were rare [39]. In this study, sunny skies were simply identified when the half-day mean k t was larger than a threshold in the morning and afternoon, respectively. 1.2 times of two-month running average k t was used as the threshold so that the number of obtained sunny half-days was close to the mean value observed in recent decades. Clearness index were plotted against solar elevation angles and fitted by cubic polynomials in the sunny mornings and afternoons, respectively ( Fig. 1). The data falling away from the major patterns were excluded [39]. Different from the studies in the forests [22,39], k t observed in the wheat field was lower in the afternoon than in the morning (Fig. 1a). Air pollution may reduce the atmospheric transparency, resulting in a small incident solar radiation at the land surface. The reduction was more evident in the afternoon than in the morning (Fig. 1b).

NEE Light response and influencing factors
The relationship between daytime NEE and PAR in a winter wheat field is shown in Figure 2. Light response curves in April and May were similar but both of them differed from those in March. From 2003 to 2006, P max was 46.664.0 mmol CO 2 m 22 s 21 , a 0.05960.006 mmol mmol 21 and R d 5.360.3 mmol CO 2 m 22 s 21 in April2May, two or three times as high as those in March (Table 1). However, the inter-annual coefficients of variation (CV) for a, P max and R d were great in March due to large variations of LAI, daytime mean T a , VPD and SWC among years. The variations in VPD and SWC resulted from significant changes in precipitation in March among years (Tables 1 and 2). In March 2005, a was so small and P max was so large that the light response curves were actually close to a line (Table 1 and Fig. 2). Table 3 indicates that NEE r declined with the increases in LAI, g c , T a and SWC (P,0.001) but increased with the increase in VPD (P,0.001). Positive/negative NEE r means measured NEE was higher/lower than simulated NEE; or measured CO 2 uptake was lower/higher than simulated CO 2 uptake. The impacts of biophysical factors on NEE r varied in different months. Factors influencing NEE r were sorted as LAI, SWC, g c , T a and VPD in March; g c , LAI, VPD and T a in April; and LAI, g c , VPD and T a in May. The effects of SWC on NEE r were significant in March but insignificant in April and May (Table 3).
Shifted trends in NEE r versus variables illustrated the limit in linear regression analysis (Fig. 3). With the increase in T a and VPD, NEE r turned to increase significantly (P,0.001) when T a was more than 25uC or VPD more than 1.121.3 kPa. With an increase in g c , the NEE r decline disappeared, or even turned into an increase (P,0.01) when g c exceed 26 mm s 21 in April or 14 mm s 21 in March and May ( Fig. 3 and Table 4). No reverse trends were found for NEE r versus LAI and SWC.

Light response parameters: seasonal variation and influencing factors
Light response parameters varied seasonally, with peaks in April or May ranged from 64.2 to 86.1 mmol CO 2 m 22 s 21 for P max , from 0.070 to 0.087 mmol mmol 21 for a and from 6.5 to 8.5 mmol CO 2 m 22 s 21 for R d during the year 200322006 (Fig. 4). The seasonal patterns of P max were similar with LAI but different from T a and VPD which have an increasing tendency from March to May. Compared with other years, P max peaked earlier in 2004 when LAI reached the maximum in mid-April due to fast warming in spring (Fig. 4). Figure 5 shows the scattered points of biophysical factors (LAI, g c , T a , VPD and SWC) versus NEE light response parameters (a, P max and R d ) at 5-day scale during the main growing season of winter wheat. Stepwise multiple linear regression models were used to assess their relationships ( Table 5). P max increased with the increase in LAI, g c and T a (P,0.001) but reduced with the increase in VPD (P,0.001). The factors influencing P max were sorted as LAI, g c , T a and VPD. a was proportional to ln(LAI), g c , T a and VPD (P,0.001). The impacts of LAI, g c and T a on a were larger than that of VPD. R d was proportional to T a (P,0.001). It was better to express the relationship between R d and T a using exponential equation instead of linear equation. However, the influences of SWC on all light response parameters were insignificant (Table 5).
During the growing season of winter wheat, a increased linearly with an increase in R d (P,0.001) (Fig. 6). a and R d enhanced firstly and then declined with the increase in P max . The maximum a and R d appeared when P max was around 50 mmol CO 2 m 22 s 21 . Their relationships could be expressed by quadratic polynomials (P, 0.001) (Fig. 6).   (Table 6). On average, P max , a and R d under cloudy sky conditions were 9%, 11% and 2% higher than those under sunny sky conditions, respectively. However, owing to large variation in a, P max and R d among years, their differences between two sky conditions were smaller than their total standard errors calculated by Eq. (4) ( Table 7), indicating unremarkable differences in light response parameters between sunny and cloudy sky conditions. NEE Light response in cloudy days differed from that in sunny days. The differences between two light response curves were significant in 2003 and 2005 but insignificant in 2004 and 2006. Cloudy and sunny light response curves were so close in 2004 and 2006 that their differences were within the confidence intervals due to scattered points of NEE versus PAR (Fig. 7). At the same PAR, simulated net CO 2 uptake was 33%, 13%, 23% and 8% (averaged 18%) higher under cloudy sky conditions than under sunny sky conditions in 2003,2004,2005 and 2006, respectively. The difference between two sky conditions was more than their standard errors ( Table 7), suggesting that net CO 2 uptake under cloudy sky conditions was significantly higher than that under sunny sky conditions.
As mentioned above, LAI and g c were key factors affecting NEE. After neglecting the effects of LAI on NEE during the period with large LAI, the NEE2g c relationships were investigated in several radiation classes (Fig. 8). The correlations between NEE and g c were described by logarithmic equations under strong, moderate and weak radiation, respectively (P,0.001). With the increase in ln(g c ), net CO 2 uptake enhanced more quickly under strong radiation than under low radiation. At the same g c , net CO 2 uptake was higher in cloudy days than in sunny days. Nevertheless, the differences between two regression lines were insignificant for all PAR classes because they were within the wide confidence intervals owing to scattered points of NEE versus ln(g c ) (Fig. 8).

Factors influencing NEE light response
In boreal and temperate forests, NEE r generally declines with increasing T a and VPD. These trends turn to rise when T a exceeds 20225uC [7,17] and VPD exceeds 1.121.3 kPa [7,18]. In this study, similar phenomena were observed beyond the same thresholds in a temperate cropland (Fig. 3 and Table 4). However, in the tropical forest, NEE r always increases with the increase in T a and VPD without reverse trends [40] (Loescher at al., 2003). It may be ascribed to high air temperature and humidity with a narrow range in the tropic region.
Net CO 2 uptake decreased with the increase in T a and VPD (Fig. 3), indicating that temperature and water stress occurred. Under high temperature (T a .25uC), photosynthesis was prohibited and soil and plant respirations were great, resulting in a lower net CO 2 uptake than the one simulated by Eq. (3). On the other hand, VPD controls photosynthetic rate through influencing stomatal closure. Under higher VPD, stomatal opening is smaller,  leading to the decrease of photosynthetic rate [41]. Moreover, the decline of photosynthetic rate was likely due to low leaf water potential caused by high transpiration rates [42]. Ecosystem respiration is large under higher VPD because of high air and soil temperatures [41]. Hence, net carbon uptake by the cropland dropped. In this study, lower VPD threshold in April than in May (Table 4) implied that wheat plants were more sensitive to the drought at the fast growing stage than at the senescence stage.
Around the critical value, the reverse trends in NEE r versus g c (Fig. 3 and Table 4) suggested a shift from stomatal limitation to non-stomatal limitation on daytime net CO 2 exchange. At dawn, dusk or in cloudy days, photosynthesis was inhibited by weak radiation even though g c was large. Higher g c threshold in April than in March and May (Table 4) indicated that net CO 2 uptake was more sensitive to stomatal behavior at the rapid growth stage than at the slow growth senescence stage. Therefore, NEE r was  Table 5. Effects of biophysical drivers (LAI, g c , T a , VPD and SWC) on light response parameters (P max , a and R d ) at the 5-day scale estimated by the stepwise multiple linear regression models during the growing seasons of 2003-2006. The meanings and units of biophysical factors and light response parameters were the same as Tables 1, 2     mainly controlled by g c instead of LAI in April. However, when g c was limited in colder March or drier May, LAI became the most important factor affecting NEE r in these months. The significant impacts of SWC on NEE r in March (Table 3) may be ascribed to large variation in soil moisture before and after irrigation at the turning green stage of winter wheat (Fig. 4). Nevertheless, the sufficient irrigation concealed these effects in the following two months (Table 3).

Factors influencing light response parameters
Photosynthetic light response parameters were not constant but varied seasonally [9,11,19]. For wheat canopy, mean P max in April2May observed in this study (49 mmol CO 2 m 22 s 21 ) was lower than the results (62267 mmol CO 2 m 22 s 21 ) obtained in a winter wheat field in central Germany [5] and a spring wheat field in Manitoba, Canada [43]. Mean a in April2May obtained in this study (0.056 mmol mmol 21 ) was lower than the value (0.063 mmol mmol 21 ) observed in a winter wheat field in central Germany [5], but higher than that (0.036 mmol mmol 21 ) obtained in a spring wheat field in Manitoba, Canada [43]. The magnitudes of light response parameters varied among different studies might be attributed to the discrepancy in temperature, humidity, canopy structure or/and other factors. P max had positive correlations with LAI [9,16,18,19], air temperature [17,18] and SWC [9,18] but negative correlation with VPD [9,44]. A similar phenomenon was found in this study ( Fig. 5 and Table 5). Among all factors, LAI was regarded as a key factor controlling P max [9,16,18]. P max increased linearly [18] or nonlinearly [16,19] with increasing LAI. The nonlinear relationship was expected as the leaves shade each other in the ecosystem with higher vegetation density [16]. In the cropland, the linear correlation between P max and LAI ( Fig. 5 and Table 5) illustrated a suitable density range for canopy light intercept.
LAI determines photosynthetic area while g c controls the photosynthetic intensity. Owing to the strong link between g c and photosynthetic rate, the influences of environmental factors (T a , VPD and SWC) on photosynthesis may be ascribed to their effects on g c [45]. Hence, factors influencing P max was sorted by LAI, g c , T a and VPD (Table 5). Our results in an irrigated wheat field differed from those obtained by Zhang et al. [9] who found that P max was only affected by LAI and VPD in a dry wheat field.
In this study, a was affected by LAI, T a , g c and VPD in an irrigated wheat field (Table 5), differing from the result reported by Zhang et al. [9] who found that a was only influenced by LAI in a dry wheat field. a represents weak light use efficiency by the plants. At dawn and dusk, the weak light is mainly composed of diffuse radiation. Under low LAI, less diffuse radiation was obtained by  For each year, NEE r was calculated at the same PAR using Eq. (3) and the light response parameters in Table 6. The mean values were obtained for two sky conditions and total standard error was computed using Eq. (4). The meanings of P max , a and R d were the same as Tables 1. Significance of the difference was '' * '' for P,0.05 if the absolute difference between two sky conditions was greater than the total standard error. doi:10.1371/journal.pone.0089469.t007 the short canopy. When LAI enlarged, the canopy structure was beneficial to receiving scatter light from each direction. Less diffuse radiation was intercepted by the sparse canopy when LAI decreased with leaf senescence. Moreover, at dawn and dusk, low radiation led to small g c . Under cold and moist conditions, a was more sensitive to T a and g c than to VPD. Generally, the key factor influencing respiration was temperature (Table 5). It likely changed to moisture under dry conditions. For instance, Zhang et al. [9] obtained that R d had positive correlations with LAI and soil moisture in a dry wheat field. In this study, R d was correlated with T a significantly but not related to VPD and SWC (Table 5 and Fig. 5) due to sufficient irrigation in the winter wheat field. Furthermore, the effect of LAI on R d was unremarkable (Table 5 and Fig. 5) because R d was composed of soil and plant respiration and only the later part was related to LAI [16].
After investigating the flux data observed in the grasslands and croplands over the world, Gilmanov et al. [19] pointed out that light response parameters were correlated with each other. Because P max was proportional to LAI, a and R d correlated to P max was the same as correlated to LAI ( NEE Light response and NEE2g c relations under various sky conditions P max in a wheat field was greater under cloudy skies than under sunny skies in most of years (Table 6), in agreement with the studies in the temperate forests [8,21,22]. Different from the temperate area, Zhang et al. [22] observed a reverse phenomenon in a subtropical forest: net CO 2 uptake under cloudy skies was less than that under sunny skies. It may be owing to much lower radiation intensity in cloudy days compared with sunny days in the subtropical forest. In the winter wheat field, differences of P max and a between two sky conditions were insignificant (Table 6) due to large variations in P max and a among years. The values of a in this study were consistent with those reported by Dengel and Grace [8] and Zhang et al. [22] for forests, but different from the Figure 8. Linear relationship between NEE and logarithmic canopy conductance (g c ) under sunny and cloudy sky conditions (solid lines) in a winter wheat field during the period of large LAI (0.8LAI max ,LAI,LAI max ). Between the upper and lower dotted lines were 95% confidence intervals. doi:10.1371/journal.pone.0089469.g008 results observed by Hollinger et al. [10], Rocha et al. [21] and Suyker et al. [24] who obtained a larger a under cloudy conditions in the forests and an irrigated maize cropland.
In spite of no significant differences in light response parameters between sky conditions, simulated net CO 2 uptake under cloudy sky conditions was greater than that under sunny sky conditions at the same PAR (Fig. 7). A similar phenomenon was observed in the temperate forests [8,21] and the irrigated cropland [24]. More CO 2 uptake by the canopy under cloudy skies than that under sunny skies was owing to many reasons. Firstly, the proportion of diffuse radiation increases under cloudy sky conditions and more light can reach below leaves of the canopies [46]. Photosynthetic rates of shaded leaves are promoted by the delivery of diffuse radiation [47]. In addition, compared with shaded leaves, the phenomenon of saturating photosynthesis easily happens for sunlit leaves because shaded leaves often illuminate brightly [48]. Secondly, canopy conductance was usually higher under cloudy sky conditions than that under sunny sky conditions. Diffuse radiation enhances canopy stomatal conductance mainly due to the reduction in VPD and blue light enrichment within the canopy during cloudy and overcast weather [49,[50][51][52]. Both low VPD [53] and blue light [54,55] can stimulate stomatal opening and then enhance photosynthetic rate.
In this study, NEE2g c relationships were expressed by logarithmic equations (Fig. 8), differing from the result obtained by Dengel and Grace [8] who pointed out that net CO 2 uptake enhanced linearly with the increase in g c . Under strong radiation, net CO 2 uptake enhanced quickly with increasing g c (Fig. 8), indicating that g c was the main factor affecting NEE. Compared with cloudy days, net CO 2 uptake was more sensitive to the variation in g c for higher temperature and VPD in sunny days.