Introduction
Tropical forests are characterized by a warm and humid climate (Corlett 2011); however, there is currently little consensus on whether climate change will affect tropical forests. Paleoecological studies show that neotropical vegetation largely persisted after a 3 to 5 °C warming during the Paleocene–Eocene Thermal Maximum (Jaramillo et al 2010). However, this historical warming was short-lived and considerably slower than current warming and future warming predicted for the next century. A survey of the temperatures of broad-leaved forest land cover suggests that climatic warming could have severe consequences for tropical floras (Wright et al 2009). Closed-canopy forests are found in areas with a mean annual temperature below 28 °C, whereas areas with mean temperatures above 28 °C support shrubs and grasses instead of broad-leaved evergreen trees. Given that excessively high temperatures are typically associated with a high evaporative demand and dry climate, the absence of closed-canopy forests in areas with temperatures above 28 °C could also be a consequence of water limitation. This past record and the distribution of tropical forests suggest a temperature limit, and therefore the ecosystem sensitivity to this threshold needs to be further studied.
Photosynthetic performance, the basis for carbon sequestration and ecosystem production, is temperature dependent. In general, the light-saturated photosynthetic rate increases with temperature to a peak, which is followed by a decline (Sage and Kubien 2007). It has been suggested that current temperatures in regions supporting tropical forests are very close to or even exceed their photosynthetic optimum temperature (Topt) (Doughty and Goulden 2008). This is a potentially ominous warning sign for our warming Earth. Tropical forests store large amounts of carbon in biomass (Dixon et al 1994). Consequently, a slight perturbation in tropical carbon fluxes could have wide-ranging effects on global atmospheric CO2 concentrations (Anderegg et al 2015). Temperatures in excess of Topt could result in a sharp decline in photosynthetic carbon sequestration in tropical forests (Doughty and Goulden 2008, Vårhammar et al 2015). A decline in CO2 uptake by the forests could in turn result in an increase in atmospheric CO2, which would further accelerate warming through positive feedback.
Model simulations indicate that tropical forests are currently not at their high-temperature threshold. With the aid of widely-used process models, Lloyd and Farquhar (2008) showed that the temperature of tropical forests was still well below Topt. They argue that the apparent decrease in photosynthetic rate with increasing temperature is predominantly an indirect effect of stomatal closure (30%) and not a direct effect of warming on mesophyll processes (2%). Temperature increases could reduce photosynthesis, either directly through inhibiting the activity of photosynthetic enzymes and electron transport, or indirectly through decreasing stomatal conductance (Farquhar et al 1980). A linear increase in temperature could lead to an exponential growth of water vapour pressure deficit (D) (Campbell and Norman 1998), whereas stomatal aperture (conductance) decreases with increased D (Damour et al 2010). Since CO2 enters the mesophyll through stomata, intercellular CO2 (ci) and photosynthesis decrease when stomatal conductance (g) declines (Farquhar et al 1980). In contrast to the direct effects of temperature on the photosynthetic apparatus, a reduction in photosynthesis caused by increased stomatal resistance could be offset, at least partly, by elevated CO2 (Lloyd and Farquhar 2008). Elevated CO2 increases Topt by reducing photorespiration and stomatal resistance, which has a positive effect on the acclimation potential of photosynthesis. Moreover, stomatal closure reduces transpiration and subsequently reduces its cooling effect (Doughty 2011). This could in turn lead to excessively high temperatures at the leaf level, which could cause irreversible damage to the photosynthetic machinery (Berry and Björkman 1980, Doughty 2011).
In this study, to examine the potential effects of climate change on forest photosynthesis, we first quantified the Topt of ecosystem photosynthesis (ToptE) for seven tropical forests across different continents. We then analyzed the relationship between ToptE and mean growing season air temperature (Ta) to confirm the widely held consensus that these parameters increase simultaneously. Ecosystem physiological parameters were then inverted using a big-leaf analogized process model driven by ecosystem photosynthesis measurements. Further, we tested the hypothesis proposed by Lloyd and Farquhar (2008), which suggests that stomatal processes play a prominent role in determining Topt, and that increasing ambient CO2 concentrations will increase tropical forest Topt, which would imply that these forest are not as vulnerable to climate change as may have been indicated by Doughty and Goulden (2008). Finally, we discuss the implications of different climate warming scenarios on ecosystem photosynthesis in tropical forests.
Material and methods
Studied sites
Tropical rainforests are primarily distributed in the Amazon, Southeast Asia, and Africa. In the present study, we investigated seven tropical forests, four of which are located in Southeast Asia and three are in the Amazon (table 1). All three Amazonian sites are located near the equator (latitude ∼3°S): from west to east K34, K67, and K83. The four Asian rainforests are in two different locations: two sites (PDF and PSO) are near the equator (∼2°S or N) and the other two Thailand forests (MKL and SKR) are located at ∼14°N. The selected Asian rainforests are dominated by trees in the Dipterocarpaceae; the exception being the peat swamp forest of the PDF site. Canopy height typically exceeds 30 m, although in the peat forest the maximum height is approximately 26 m. The forest at the K83 site has previously been selectively logged for experimental purposes (Goulden et al 2004). All the studied forests have a year-round growing season, with the exception of MKL, in which a proportion of the trees shed their leaves during the late dry season. For additional detailed information on these sites and instrumentation please refer to the previously published studies of Hirata et al (2008) and Restrepo-Coupe et al (2013).
Table 1.Site information.
Lat. | Long. | Age | Hc | LAI | Ta | ppt | System | Anemometer | IRGA | Period | Country | |
---|---|---|---|---|---|---|---|---|---|---|---|---|
Amazonian | ||||||||||||
K34 | 2° 36' S | 60° 12' W | Primary | 30∼35 | 4.7 | 26.7 | 2286 | CP | Wind master, Gill | Li-6262, Li-Cor | 1999∼2006 | Brazil |
K67 | 2° 51' S | 54° 58' W | Primary | 35∼40 | 6.0 | 24.8 | 1811 | CP | CSAT3, Campbell | Li-6262, Li-Cor | 2002∼2006 | Brazil |
K83 | 3° 3' S | 54° 56' W | Selective logged | 35∼40 | 4.9 | 24.8 | 1811 | CP | CSAT3, Campbell | Li-7000/6262, Li-Cor | 2000∼2004 | Brazil |
SE Asia | ||||||||||||
MKL | 14° 34' N | 98° 50' E | ∼30 (2008) | 30 | 2∼3 | 27.5 | 1650 | CP | Wind master, Gill | Li-6262, Li-Cor | 2003∼2004 | Thailand |
2° 20' N | 114° 2' E | — | ∼26 | 5 | 26.3 | 2231 | OP | CSAT3, Campbell | Li-7500, Li-Cor | 2002∼2005 | Indonesia | |
PSO | 2° 58' N | 102° 18' E | Primary | 35∼45 | 6.52 | 25.3 | 1804 | OP | SAT550, Kaijo | Li-7500, Li-Cor | 2003∼2009 | Malaysia |
SKR | 14° 29' N | 101° 55' E | Mature | 35 | 3.5∼4.0 | 26.2 | 1240 | CP | Wind master, Gill | Li-6262, Li-Cor | 2001∼2003 | Thailand |
'—', no data available; Lat., Latitude; Long., Longitude; Age, stand age (yr); Hc, canopy height (m); LAI, leaf area index; Ta, mean annual temperature (°C); ppt, precipitation (mm); System, the open (OP) or closed path (CP) eddy covariance system; IRGA, the infrared gas analyzer model.
Eddy flux observations and data processing
The CO2 movement in the lower atmospheric boundary layer is primarily driven by turbulence that can be measured using the eddy covariance technique (EC) (Baldocchi 2003). Photosynthetic rates were quantified by examining the ecosystem–atmosphere CO2 exchange. The daytime CO2 exchange measured using EC apparatus is conceptually equal to net ecosystem photosynthesis (table 2).
Table 2.Terms (and their abbreviations) used at the leaf-level and the corresponding abbreviations used at ecosystem-level. Pg: leaf gross photosynthetic rate; GPP: gross primary production of ecosystem; RL: leaf respiration; RE: ecosystem respiration, which is the sum of autotrophic and heterotrophic respiration; TL: leaf temperature; Ta: air temperature near the top of the forest canopy; ciL: leaf intercellular CO2 concentration; ciE: bulk canopy intercellular CO2 concentration; DL: leaf-to-air water vapour deficit; DE: atmospheric water deficit; gs: stomatal conductance; Gc: canopy conductance; Pn: net leaf photosynthesis rate, where Pn = Pg − RL; NEE: net ecosystem exchange, where NEE = GPP − RE.
General description | Leaf | Ecosystem |
---|---|---|
Gross photosynthesis | Pg | GPP |
Dark respiration: Rd | RL | RE |
Temperature | T | Ta |
Intercellular CO2 concentration: ci | ciL | ciE |
Water vapour deficit: D | DL | DE |
Stomatal or canopy conductance: g | gs | Gc |
Maximum carboxylation rate: Vcmax | VcmaxL | VcmaxE |
Maximum electronic transport rate: Jmax | JmaxL | JmaxE |
Conductance sensitivity: g1 | g1L | g1E |
Temperature curve factor: S | SL | SE |
Temperature curve factor: H | HL | HE |
Net photosynthesis rate | Pn | NEE |
Optimal temperature: Topt | ToptL | ToptE |
We collected EC flux data for the seven forests from flux networks. The fluxes have a temporal resolution of 30 min or 1 h, and span at least two years. The major flux data used in this study include the following: net CO2 exchange (NEE, after storage flux correction), latent heat flux (LE), sensible heat flux (Hs), net radiation flux (Rn), and soil heat flux (G). In addition to flux data, we also used meteorological measurements, including air temperature (Ta, °C), relative humidity (hs, %), water vapour pressure deficit (DE, kPa), and soil water content (SWC, m3 m−3). For reproducibility, the data are available at the following sites:
AsiaFLUX dataset: https://db.cger.nies.go.jp/asiafluxdb/
Determining the optimum temperature (Topt) of photosynthesis
Determining of Topt could be done based on gross photosynthesis (GPP) dataset. The advantage of this way is reducing uncertainties related to respiratory processes which are most significant in eddy flux cases (Yi et al 2000, Yi et al 2004). However, a reliable method to obtain GPP by portioning NEE is currently unavailable either because light inhibition of leaf respiration or inconsistence of temperature dependency of autotrophic respiration (cf. Yi et al 2004 and reference therein).
Topt could also be determined by fitting a peak function to the temperature response of light-saturated photosynthesis (Lange et al 1974). In leaf-level studies, temperature is specified to leaf temperature (TL), and light in the leaf chamber is set to a saturating level during CO2 exchange measurements. There are, however, some modifications required when these equations are applied at the ecosystem level. Instead of leaf temperature, we used air temperature near the canopy level (Niu et al 2012). Therefore, in the present study, the ecosystem photosynthesis Topt (ToptE) was determined in terms of optimum air temperature. To determine values for light-saturated photosynthesis, we omitted all data points below site-specific saturating light levels. The site-specific light saturation point was calculated by applying a non-rectangular hyperbola to the stand-level photosynthesis–light response curve (Lasslop et al 2010).
There are several peak functions that could be used to fit the temperature response curve to determine Topt. We adopted a modified function of the model proposed by June et al (2004):
where NEEsat is the measured net ecosystem photosynthesis rate under light saturation (μmol m−2 s−1) (note that a positive NEE indicates photosynthesis uptake, in order to make it comparable to that of leaf-level conventions), Tk is the ambient temperature in degrees Kelvin, R is the gas constant, and NEE25 (μmol m−2 s−1), b, c, and ToptE (Kelvin) are fitted parameters.
The Farquhar–von Caemmerer–Berry (FvCB) model
A process-based photosynthesis model was used in this study. The model is a combination of the Farquhar–von Caemmerer–Berry photosynthesis (FvCB) model (Farquhar et al 1980) and the Ball–Berry stomatal conductance model (Ball et al 1987), with some additional parameterization information provided by von Caemmerer et al (2009). The detailed model equations are listed in table A1 in the appendix. We used an iteration method to solve intercellular CO2 (ci). The model was coded in the C++ environment and is available upon request.
The big-leaf analogy
The factors used in determining Topt can be summarized as photosynthetic biochemical (biochemical hereafter), respiratory, and stomatal processes (Hikosaka et al 2006, Lin et al 2012). In order to separate the relative contributions of each process, we implemented the FvCB model with a big-leaf analogy.
Firstly, the ecosystem as a whole was abstracted into a 'big leaf'. This is consistent with the philosophy of the EC method and enabled us to directly use the leaf-level FvCB model at an ecosystem-level. The EC system measures the gas exchange at ecosystem level analog to leaf chamber measurements do a small scale (table 2). Thus, daytime net ecosystem exchange (NEE) was regarded as equivalent to net photosynthesis rate (Pn) at the leaf level. At the ecosystem level, the air temperature near the canopy (Ta), air water vapour deficit (DE), and canopy bulk intercellular CO2 concentration (ciE) corresponded to leaf temperature (TL), leaf-to-air water vapour deficit (DL) and intercellular CO2 concentration (ciL) at the leaf level, respectively. Critically, ecosystem respiration (RE) at the ecosystem level was considered analogous to leaf respiration (RL) at the leaf level. After analogizing, ecosystem terms were derived that corresponded to the leaf terms, and the FvCB model was applied at the ecosystem level and driven by ecosystem measurements.
This type of big-leaf analogy differs from that of the big-leaf model (de Pury and Farquhar 1997), in that here the ecosystem as a whole was treated as a big leaf. Therefore, the parameters derived here for the ecosystem are not directly comparable to those used in leaf studies. The overall motivation for us in conducting this analogy was to make the parameter inversion as simple as possible but with necessary physiological considerations. Parameter inversion is very sensitive to initial parameter values because of many non-linear processes, and consequently it is more practical and helpful to construct simple models with certain assumptions than to execute a complex multi-layer model (Wang et al 2007).
Inverting photosynthesis parameters by combining the FvCB model and ecosystem fluxes
Parameters of the FvCB model at the ecosystem level were inverted from eddy flux observations after abstracting. For inversion, we used the Levenberg–Marquardt optimization algorithm. We also examined whether the inversion method would have a strong impact on the inverted parameters. A Bayesian statistical method, the 'adaptive population Monte Carlo approximate Bayesian computation' method (Lenormand et al 2013), was included for comparison. The algorithm pseudo-code of the Bayesian method was presented in Zeng et al (2017).
Results
The ToptE of tropical forests
The temperature dependence of light-saturated ecosystem photosynthesis (NEEsat) is shown in figure 1 general, there was a clear unimodal pattern for most of the sites. The ToptE determined by fitting equation 1 to the observations varied from 23.7 to 28.1 °C across sites.
A close relationship was found between mean annual air temperature (Ta) and ToptE (figure 2(a)). Since most tropical forests maintain a year-round growing season, the mean annual Ta could roughly be treated as the growth temperature. Therefore, tropical forests growing under higher growth temperature tend to have a higher ToptE. The slope of the linear relationship is close to one (1.12). ToptE was also related to mean air temperature under light saturated condition (figure 2(b)). When the sites with seasonally climate were omitted, a very close relationship was found between ToptE and mean air temperature under light saturated condition.
Contribution of physiological parameters to the change in ToptE across sites
The goodness-of-fit was shown in figure 3 when implemented the FvCB model to these datasets. In general, the model fitted results have a good relationship with that of observations. It suggests high reliability of these inverted parameters (table 3). Principal component analysis of these parameters identified three components that could explain over 86% of the variance. However, none of the three components were significantly correlated with ToptE (data not shown). A further correlation analysis showed that the activation energy of ecosystem respiration (RE) was the only parameter significantly correlated with ToptE (figure 4(a)). We found that two sites (PDF and MKL indicated by open circles in figure 4(b)) differ from the remaining sites with respect to the relationship between ToptE and stomatal sensitivity (g1E): with the exception of PDF and MKL, the sites show a negative correlation between ToptE and g1E. These two excluded sites have special water conditions which discussed latter in the discussion section.
Table 3.Parameters of a photosynthesis model inverted using a nonlinear regression method.
SiteRate at 25 °C (μmol m−2 s−1)Activation energy (J)SEHEg1EVcmaxEJmaxERdEVcmaxEJmaxERdEK34 | 294 | 190 | 9.03 | 65 340 | 76 340 | 10 000 | 706 | 218 782 | 2.31 |
K67 | 297 | 212 | 10.00 | 60 503 | 35 255 | 66 673 | 702 | 217 451 | 3.27 |
K83 | 290 | 215 | 9.64 | 63 009 | 33 047 | 39 721 | 408 | 126 451 | 2.97 |
MKL | 295 | 235 | 9.15 | 63 724 | 56 013 | 25 547 | 630 | 195 119 | 3.83 |
230 | 212 | 10.00 | 61 913 | 38 906 | 44 375 | 697 | 215 911 | 3.40 | |
PSO | 185 | 209 | 9.44 | 59 110 | 70 404 | 15 779 | 798 | 247 355 | 2.60 |
SKR | 252 | 239 | 9.59 | 62 342 | 56 916 | 35 102 | 710 | 220 120 | 2.69 |
Vcmax, maximum Rubisco activity; Jmax, maximum electron transport rate; Rd, dark respiration rate; S, a term similar to an entropy factor, H, the rate of decrease in the function above the optimum, g1, stomatal sensitivity factor.
The contribution of stomatal processes in determining ToptE
The PSO site, which has eight years' continuous flux data, was taken as an example of per-humid site to illustrate the contribution of stomatal processes in determine ToptE. The overall ToptE for the PSO site was 26.5 °C for the entire DE range, as shown by the dashed line in figure 5(a) When the whole dataset was divided into different DE levels, we found that the light-saturated photosynthesis rate (NEEsat) increased with Ta even at values greater than 30 °C (see data points and regression line in different colours, figure 5(a). Only at the highest DE level was NEEsat found to decrease with Ta. In this regard, ToptE (when NEEsat starts to decrease with an increase in Ta) should at least be 30 °C when DE is controlled. This contrasts with the value of 26.5 °C obtained using the full DE range and implies strong stomatal control of ToptE. We also inferred parameters of FvCB model for all DE sub-ranges (table 4). Most of these parameters showed a unimodal pattern along with increase of DE. A similar case was found in a site with strongly seasonal climate (figure 5(b)).
Table 4.The inverted parameters under different water vapor pressure deficit (DE) levels. This was carried out in the per-humid Pasoh site. This table could correspond to figure 4(a).
DE levels (kPa) | <0.7 | 0.7∼0.9 | 0.9∼1.1 | 1.1∼1.2 | 1.2∼1.3 | 1.3∼1.4 | 1.4∼1.5 | 1.5∼1.6 | 1.6∼1.8 | 1.8∼2.1 | >2.1 | |
---|---|---|---|---|---|---|---|---|---|---|---|---|
Rate at 25 °C (μmol m−2 s−1) | VcmaxE | 91 | 100 | 114 | 254 | 288 | 289 | 273 | 268 | 165 | 279 | 101 |
JmaxE | 125 | 123 | 107 | 110 | 140 | 118 | 129 | 129 | 107 | 91 | 102 | |
RdE | 0.10 | 0.51 | 0.10 | 1.91 | 8.10 | 5.67 | 8.50 | 8.98 | 6.59 | 5.00 | 2.90 | |
Activation energy (J) | VcmaxE | 62 714 | 59 957 | 60 577 | 64 505 | 68 994 | 70 017 | 66 453 | 66 916 | 63 578 | 66 943 | 59 348 |
JmaxE | 44 315 | 41 412 | 47 631 | 51 075 | 66 799 | 65 574 | 79 995 | 79 573 | 79 994 | 79 998 | 41 830 | |
RdE | 63 481 | 63 637 | 63 254 | 63 565 | 55 985 | 60 646 | 64 301 | 61 926 | 62 723 | 58 468 | 64 409 | |
SE | 742 | 729 | 756 | 798 | 1052 | 1126 | 903 | 896 | 859 | 874 | 742 | |
HE | 22 9765 | 22 5983 | 23 4348 | 24 7302 | 32 5870 | 34 8879 | 27 9685 | 27 7651 | 26 6097 | 27 0846 | 22 9800 | |
g1E | 11.00 | 10.51 | 13.00 | 20.43 | 36.50 | 30.72 | 13.85 | 14.31 | 10.97 | 3.44 | 2.08 |
Vcmax, maximum Rubisco activity; Jmax, maximum electron transport rate; Rd, dark respiration rate; S, a term similar to an entropy factor, H, the rate of decrease in the function above the optimum, g1, stomatal sensitivity factor.
Discussion
The ecosystem ToptE we quantified differs from leaf ToptL in several aspects. Leaf ToptL is specified to leaf temperature, not air temperature. The leaf surface is a direct light interceptor, which leads to stronger temperature variations in leaves than in ambient air, i.e. the transitional leaf temperature can easily reach 40 °C under full light (Doughty and Goulden 2008). In addition, the dark respiration term (Rd) at the ecosystem level (RE) is the sum of respirations from different organisms, litter, woody debris, and soil organic matter, whereas at the leaf level, the respiration term (RL) is specified to leaf respiration. Despite these differences, however, the ecosystem ToptE obtained in our study is very close to that of leaf ToptL. For example, two Costa Rican tropical forest species grown under a daily temperature of 27 °C showed a leaf ToptL of 27 °C (Vargas and Cordero 2013), which is close to the ecosystem ToptE value we determined for tropical forests. Similarly, Slot and Winter (2017) reported mean ToptE values of 30.4 °C and 29.2 °C for the upper-canopy leaves of 42 species in two lowland forests in Panama, which were close to the mean afternoon air temperature. The higher ToptL values reported by Slot and Winter (2017), compared with the values reported in the present study and those reported by Vargas and Cordero (2013), could be explained by the fact that upper canopy leaves experience higher light intensity and leaf temperatures compared to the whole canopy mean values. In the following sections, we discuss the possible mechanism of ToptE changes across sites, the contribution of stomatal processes to ToptE, and the implications of our findings.
The mechanisms of ToptE changes across sites
Our cross-site analysis shows that tropical forests growing in a warmer climate tend to exhibit higher ToptE (figure 2(a)). We separated the contributions of biochemical, respiratory, and stomatal processes to ToptE by means of parameter inversions. Respiratory process play a role in ToptE, as suggested by the significant relationship between ToptE and the activation energy (Ea) of respiration (RE) (figure 4(a)). This finding differs from those of leaf level studies, which indicate that leaf respiration plays a negligible role in ToptL (Lin et al 2012). At the ecosystem level, RE is the sum of the autotrophic respiration of all organisms (above- and below-ground) and heterotrophic respiration of soil organic matter and litter. At the leaf-level, however, RL reflects only leaf respiration, which typically represents a small fraction of net photosynthesis. These differences emphasize the importance of respiratory processes in studying ToptE at the ecosystem level, even though it is negligible at the leaf level.
It is known that activation energy (Ea) can represent the temperature sensitivity (Q10) of RE, the relationship of which can be expressed as follows:
The relationship between decreasing ToptE and increasing Q10 suggests that the RE of tropical forests with higher ToptE and mean Ta is less sensitive to Ta variations than forests with lower ToptE. This is consistent with previous reports showing that Q10 decreases with rising temperature, as a consequence of thermal acclimation (Tjoelker et al 2009, Slot and Kitajima 2015).
When we excluded PDF and MKL (two secondary forests with unique hydrological conditions) from the analysis, we identified a strong correlation between g1E and ToptE(figure 4(b)). The forest in the PDF site is drained peat swamp forest (Hirano et al 2007), which is generally waterlogged. By contrast, the MKL site experiences seasonal water deficits (Gamo et al 2005). Since stomatal conductance or g1E is highly sensitive to water availability, it seems appropriate to treat these two secondary forests as outliers when investigating ToptE–g1E relationships.
Theoretically, g1E, which is a stomatal sensitivity factor, would be expected to increase with increasing growth temperature, as does ToptE (Leuning 1990, Medlyn et al 2011). Our numeric simulation also supports this theoretic expectation for a specific site by checking ToptE with varied g1E (data not shown). Nevertheless, we found that at the ecosystem level for tropical forests, ToptE tends to decrease with increasing g1E (figure 4(b)). Since g1E is strongly correlated with the Ea of RE (Pearson's r = 0.96), it is would be difficult to state that the strong correlation shown in figure 4(b) is solely attributable to g1E or whether it is merely an indirect reflection of the Ea–ToptE relationship shown in figure 4(a). In situ warming experiments at both leaf and ecosystem levels might be helpful in reconciling these contrasts (Cavaleri et al 2015).
Interestingly, we found that biochemical processes did not play a significant role in ToptE changes across sites. Traditionally, Topt acclimation studies have primarily focused on biochemical processes (Hikosaka et al 2006). However, in the present study these key processes were found to make a negligible contribution to ToptE changes across sites. Subsequent to further confirmation that thermal acclimation of respiration rather than biochemical processes is a more important determinant of ToptE, these findings should be implemented in global change models (Lombardozzi et al 2015).
The role of stomatal processes in determining ToptE
Previous studies have shown that stomatal processes are potentially important (Lin et al 2012, Duursma et al 2014, Slot and Winter 2017), or even the most important factors determining ToptL (Lloyd and Farquhar 2008, Rowland et al 2015). Nevertheless, how stomatal processes control Topt is still not well understood. This uncertainty is partly caused by the confounding effects of temperature and water factors on Topt.
Relative humidity (hs) and vapour pressure deficit (D) are strongly dependent on temperature (Campbell and Norman 1998). When temperature rises, the saturated water vapour pressure will increase exponentially. This, in turn, will alter both hs and D, and hence stomatal conductance. Therefore, there is an indirect effect of temperature on Topt through its effect on stomatal conductance. This effect is illustrated in figure 5. The dashed line in figure 5(a) shows the temperature response curve represented in figure 1(f) under the full DE range. However, when we divided the whole dataset into different DE levels (as shown by the different colours in figure 5(a)), it was apparent that NEEsat increases with temperature within these subsets until the temperature exceed 30 °C . This indicates that the ToptE should be at least 30 °C if there is no DE limitation on the stomatal response, which is considerably higher than the value estimated from the entire DE range (26.4 °C). This analysis lends support to the idea that stomatal processes play a significant role in determining ToptE, as temperature may indirectly influence photosynthesis through changing D.
The implications under future climate change
In the future, the Earth's surface air is predicted to become richer in CO2 and higher in mean temperature (Corlett 2011). At present, however, there seems little consensus on how tropical forest ecosystem will respond under such a climatic scenario (see the review by Lloyd and Farquhar (2008)). Our findings, however, provide certain insights into how tropical forest ecosystems might respond and could serve as complement to previous studies in our pursuit of a more complete understanding future changes in forest photosynthesis.
Firstly, we revealed the role of stomatal limitation in determining ToptE at the ecosystem level, which is largely consistent with leaf-level findings. The role of stomatal effects in shaping ToptL have been well demonstrated in leaf-level measurements (Koch et al 1994, Ishida et al 1996, Carswell et al 2000, Slot and Winter 2017) and have been verified by a leaf-level model (Lloyd and Farquhar 2008). Our ecosystem flux analysis showed that without DE limitation on stomatal conductance, tropical forests could have a higher photosynthetic performance (NEEsat) under high Ta, as indicated by their increased ToptE (figure 5). The direct implication of this finding is that factors affecting stomatal conductance will contribute substantially to the modification of ToptE. Among these factors, the most prominent is ambient air CO2 concentration. Given unchanged moisture conditions (e.g. soil water or D), ToptE is expected to increase with CO2 and will decrease stomatal limitation on photosynthesis (Lloyd and Farquhar 2008). Accordingly, this can be considered as a positive signal for tropical forests given the prospect of increasing CO2 levels.
Secondly, tropical forests in environments with higher Ta tend to have higher ToptE (figure 2), which is consistent with growth chamber cultivation experiments (Kositsup et al 2009) and cross-season observations (Lange et al 1974). This pattern suggests potential acclimation of tropical forests to Ta, which is a further adaptive strategy that will increase the resilience of tropical forest given the predicted climate warming scenarios.
Thirdly, the significant relationship between the activation energy of respiration and ToptE implies the possible thermal acclimation of RE (figure 4(a)). The temperature acclimation of ecosystem respiration, i.e. the decrease in the sensitivity of respiration to temperature changes as growth temperature increases, would have a positive effect on net photosynthesis and lead to increases in ToptE.
Collectively, our findings indicate an optimistic future for tropical forests under the predicted conditions of global climate change. Nevertheless, some uncertainties remain. Firstly, increasing CO2 will reduce stomatal conductance and water losses and hence the cooling effect of transpiration. This could potentially result in excessively high leaf temperatures and consequently heat damage and declines in photosynthesis and carbon sequestration. Furthermore, because ecosystem respiration in the tropics and subtropics is generally more sensitive to warming than that of photosynthesis (Yi et al 2010, Zhang et al 2016), it remains unclear to what extent the warming-induced increase in night-time ecosystem respiration would offset the positive effect of thermal acclimation in photosynthesis on net carbon sequestration.
Conclusions and implications
In conclusion, we quantified ecosystem ToptE for tropical forests, which ranges from 23.7 to 28.1 °C. Moreover, we found that tropical forests with higher growth temperatures tend to have higher ToptE, suggesting the acclimation potential for many tropical forests. In contrast to previous studies, however, our results show that biochemical processes make only a minor contribution to the ToptE changes across sites. Instead, respiratory processes, which are generally negligible at the leaf level, play an important role in explaining ToptE variation across sites. Consistent with leaf level studies, stomatal processes are also critical in determining ToptE at the ecosystem level. Strong D and stomatal limitation on ToptE suggests that increasing CO2 concentrations may increase the ToptE of tropical forests.
Acknowledgments
The study was supported by the C-project of Talent, Hainan University and National Natural Science Foundation of China (31660142). Brazil flux data was produced with funds from the NASA LBA-DMIP project (# NNX09AL52G) (NASA), NASA LBA investigation CD-32, and the National Science Foundation's Partnerships for International Research and Education (PIRE). We thank two reviewers and an editor board member for their insightful and constructive comments on this work. The AsiaFLUX was acknowledged for permission on accessing and using their dataset.
: Appendix
Table A1.Equations used for the leaf photosynthesis biochemical model (FvCB).
Number | Equation |
---|---|
1 | |
2 | |
3 | |
4 | |
5 | |
6 | |
7 | |
8 | |
9 | |
10 | |
11 | |
12 |
Pn: net photosynthetic rate, Pc: Rubisco-limited photosynthesis, Pj: electron transport-limited photosynthesis, Ps: export-limited photosynthesis, Rd: dark respiration rate, ci: CO2 partial pressure at the carboxylating site, Oi: O2 partial pressure at the carboxylating site, Kc: Michaelis–Menten constant of Rubisco for CO2, Ko: Michaelis–Menten constant of Rubisco for O2, Vcmax: maximum Rubisco activity, Γ⁎: CO2 compensation point in the absence of Rd, J: electron transport rate, Jmax: maximum electron transport rate, Iabs: absorbed light, fabs: leaf absorbance (∼0.85), I: light intensity, f: correction factor for the spectral quality of light (∼0.15), gs: stomatal conductance, ca: ambient CO2 concentration, cs: leaf surface CO2 concentration, gb: laminar boundary layer conductance, Ws: wind speed, Lw: leaf width, hs: relative humidity, g0 and g1 are two model parameters, Pv is used for smoothing the transition between Pc, Pj, and Ps, Ea: activation energy, TK: temperature in degrees Kelvin, R: gas constant, S: term similar to an entropy factor, H: describes the rate of decrease in the function above the optimum.