Variations in Amazon forest productivity correlated with foliar nutrients and modelled rates of photosynthetic carbon supply

Lina M. Mercado, Sandra Patiño, Tomas F. Domingues, Nikolaos M. Fyllas, Graham P. Weedon, Stephen Sitch, Carlos Alberto Quesada, Oliver L. Phillips, Luiz E. O. C. Aragão, Yadvinder Malhi, A. J. Dolman, Natalia Restrepo-Coupe, Scott R. Saleska, Timothy R. Baker, Samuel Almeida, Niro Higuchi, Jon Lloyd


The rate of above-ground woody biomass production, WP, in some western Amazon forests exceeds those in the east by a factor of 2 or more. Underlying causes may include climate, soil nutrient limitations and species composition. In this modelling paper, we explore the implications of allowing key nutrients such as N and P to constrain the photosynthesis of Amazon forests, and also we examine the relationship between modelled rates of photosynthesis and the observed gradients in WP. We use a model with current understanding of the underpinning biochemical processes as affected by nutrient availability to assess: (i) the degree to which observed spatial variations in foliar [N] and [P] across Amazonia affect stand-level photosynthesis; and (ii) how these variations in forest photosynthetic carbon acquisition relate to the observed geographical patterns of stem growth across the Amazon Basin. We find nutrient availability to exert a strong effect on photosynthetic carbon gain across the Basin and to be a likely important contributor to the observed gradient in WP. Phosphorus emerges as more important than nitrogen in accounting for the observed variations in productivity. Implications of these findings are discussed in the context of future tropical forests under a changing climate.

1. Introduction

Recent research has found large-scale variations in stem growth rates (the rate of new wood production into both boles and branches, WP) as well as tree turnover rates (mean rates of tree recruitment and mortality) across Amazonia. In general, the forests of western Amazonia are more dynamic, with younger and faster growing trees which have lower wood density than those of the central and eastern Amazon [13]. The scale of the variability is substantial, with WP varying by more than a factor 2 [1].

Different hypotheses have been proposed to explain this spatial variability in WP [1]. First, owing to the proximity of the Andes, forest soils in western Amazonia tend to be richer in nutrients than their counterparts in central and eastern Amazon [4,5], and therefore greater soil fertilities might explain the higher WP found in the west. Indeed, WP has been found to be related to soil P and soil N concentrations and soil C : N ratios for 59 sites across the Amazon Basin, and also partly related to the amount of rainfall and its spatial distribution [6]. Trees with higher foliar nutrient contents are generally associated with these more fertile soils of the western part of the Basin [7], and one possibility is that the influence of soil fertility on WP is exerted via effects of foliar [N] and/or foliar [P] on the canopy-level gross primary productivity (GP).

A second possibility is that GP actually varies very little across the Amazon Basin, but that different patterns of carbon allocation to respiration and/or other above- and below-ground organs could also explain the observed variations in stem growth [1]. However, faster growing trees on high-nutrient soils from western sites seem to allocate nearly the same proportion of productivity to above and below ground as their counterparts in the slower growing central and eastern Amazon forests [8]. This result does not, however, exclude the possibility that the slower growing eastern forests expend a greater proportion of their GP on respiration than on growth [9].

In this study, we test the first hypothesis by quantifying the extent to which variations in simulated GP across the Amazon Basin can explain observed variations in WP. To investigate the likely variability in GP, we undertake a basin-wide application of a tropical forest canopy gas-exchange model, which has already been calibrated and validated at various sites across the Amazon Basin [10].

Under the assumption of nitrogen (N) limitation, leaf photosynthesis is usually modelled based on the measured linearity between photosynthetic capacity and foliar N content. This reflects the large investment of foliar nitrogen in photosynthetic machinery [11]. Nevertheless, for tropical ecosystems, it has been suggested that phosphorous (P) rather than nitrogen (N) may constrain productivity in lowland tropical rain forests. A relative abundance of N in tropical rainforests has been suggested from foliar and soil δ15 N measurements [1214] as well as by high rates of nitrogen oxide emissions and considerable losses of nitrogen through leaching processes in many tropical forest systems [9]. This is consistent with the suggestion first made by Vitousek [15] that most tropical forests may be phosphorus- rather than nitrogen-limited and consistent with the few available studies showing a close correlation between photosynthesis and foliar [P] for tropical forest species [16]. In a cross-biome analysis of the influence of P on the linear relationship between the light-saturated rate of photosynthesis (Amax) and foliar [N], it was found that the slope of such linear relationship increases with leaf [P] [17]. This suggests that in P-limited ecosystems, the relationship between Amax and foliar [N] is constrained by low [P] availability.

Most recently, Domingues et al. [18] implemented P limitation into the main photosynthetic parameters of the Farquhar & von Caemmerer [19] photosynthesis model. The parametrizations for maximum RuBisCO activity (Vmax) and electron transport capacity (Jmax) used leaf-level nutrient and photosynthesis measurements, taken across a precipitation gradient incorporating different types of woody species in West Africa. Initial tests showed that the model (which allows for both Vmax and Jmax to be limited by either N or P) can also successfully predict leaf-level photosynthetic rates for tropical trees in the Cameroon, Bolivia and Australia [18].

Here, we use an ecosystem canopy-scale photosynthesis model, validated using flux tower data from five sites in the Brazilian Amazon [10], to simulate GP for 38 sites across the Amazon Basin incorporating possible nitrogen and/or phosphorus limitations using the parametrizations developed by Domingues et al. [18] and in situ measurements. Specifically, we examine to what extent the simulated GP with constraints of foliar [P] and/or [N] can explain the observed variability in stem growth rates.

2. Methods

(a) Data

(i) Sites

We simulate GP for 38 primary lowland rainforest sites located in Colombia, Venezuela, Ecuador, Peru and Brazil (figure 1). For all sites foliar nutrients, specific leaf area (S) and leaf area index (L) data had already been collected as part of the work of the RAINFOR Consortium. Site descriptions are given in table 1.

View this table:
Table 1.

Description of sites. Latitude (lat) and longitude (long) are given in decimal degrees, altitude (alt) in m, precipitation (precip) in mm and temperature (temp) in degrees Celsius. Soil types are taken from Quesada et al. [5] and remaining information is taken from Patiño et al. [20]. Sites in bold are nitrogen-limited according to equation (2.2).

Figure 1.

Rainforest site locations used in this study. Black and grey symbols denote N- and P-limited sites, respectively, according to Domingues et al. [18] parametrization (equation 2.2).

(ii) Leaf nutrient data (nitrogen and phosphorous) and specific leaf area

Nitrogen and phosphorus concentrations for upper canopy leaves used to derive canopy photosynthetic parameters (Vmax and Jmax) are as in Fyllas et al. [7] with around 20 trees sampled per plot, with average values for each (usually) 1 ha plot calculated here with a species abundance weighting viz:Embedded Image 2.1where Embedded Image represents the plot-level estimate for the average value of any parameter, Θ, Ns is the number of times that any species, s, occurs in the plot, ns is the number of times that species was actually sampled in the plot, m is the total number of species sampled in the plot and Θt,s represents individual measurement of parameter Θ on tree number t of species, s. Effectively, equation (2.1) gives a species abundance-weighted estimate for the plot-level average value of Θ, with the estimate of Embedded Image taking into account the fact that different species have different characteristic nutrient concentrations [7] as well as vastly varying relative abundances in different plots.

(iii) Soil phosphorous

Observed total soil phosphorus concentrations [21] available at 33 out of the 38 studied sites are taken from Quesada et al. [5].

(iv) Leaf area index

Leaf area index (L) was derived from around 20 hemispherical photographs per site [22], these usually having been taken at the end of the rainy season (S. Patiño 2001–2004, unpublished data). For the JARU site, L was taken from the estimate of Meir et al. [23]. We assume a constant L for each site, throughout the simulations.

(v) Atmospheric carbon dioxide concentration

Six-monthly atmospheric [CO2] from 1982 to 2001 (340 to 373 µmol mol−1) were used as input data [24]. Although significant diurnal variations in [CO2] may occur both above and within tropical forest canopies [25], neglecting this diurnal variability should not have a major impact on the overall estimate of GP.

(vi) Climatology

The WATCH Forcing Data were used to provide sub-daily meteorology, and are based on the ERA-40 reanalysis regridded to half-degree with adjustment of meteorological variables for changes in elevation plus monthly adjustments based on observations [26,27]. The model input variables are solar radiation (downward shortwave radiation flux), wind speed, air temperature and specific humidity at three-hourly time steps, for the period 1980–2001. The monthly adjustments include for the downward shortwave radiation flux allowance for cloud fraction and the effects of both seasonal- and decadal-variations in atmospheric aerosol-loading, as well as the use of observed monthly average temperature and diurnal temperature range. The 38 sites are covered by 25 half-degree grid squares in the WATCH Forcing Data as some of the sites are adjacent to each other (figure 1).

(vii) Data for validation

Simulated values of mean annual gross primary productivity in the absence of foliar respiration (Embedded Image, defined as GP plus canopy foliar respiration, RC [28]) are compared against bottom-up estimates based on individual measurements of the different components of the forest carbon-cycle at the five sites where data were available. The individual components included branch production, litterfall, tree growth, fine root production and autotrophic- and heterotrophic-respiration [8,2931] (Y. Malhi, L. E. O. C. Aragão & D. B. Metcalfe 2005–2006, unpublished data). The methodology was standardized across all sites, in eastern Amazonia at Caxiuana (CAX-06) and western Amazonia at Tambopata, Peru (TAM-04 and TAM-05), but for sites located in central Amazon at Manaus (MAN-K34, denoted here MAN-01) and at Tapajos (TAP-04), Embedded Image was derived from component measurements presented in the literature [31].

Eddy correlation estimates of Embedded Image have been calculated by subtracting daytime measured net ecosystem exchange from ecosystem respiration (measured during night). The method to derive ecosystem respiration from night-time measurements is explained elsewhere [3234]. Hourly night-time net ecosystem exchange measurements include a correction during periods of low turbulence using a friction velocity threshold [35]. A measure of the night-time friction velocity uncertainty of fluxes is presented using the upper and lower bounds of the friction-velocity threshold used for corrections as in Saleska et al. [34].

Mean annual values of Embedded Image from eddy correlation at MAN-01, TAP-04 and an additional tower at Jaru (southwestern Amazon) were calculated from monthly mean values [36] taking the mean of all available years with an uncertainty bound explained above. Mean annual Embedded Image at MAN-05 and CAX-06 are taken from table 12 in Malhi et al. [31].

(viii) Stem wood production (WP) and basal area growth ΔB data

Recent estimates for WP and basal area growth ΔB were available for 35 sites taken from Quesada et al. [4], and supplemented in a few cases with unpublished data from the RAINFOR database [37]. The methodology for these measurements is described in detail elsewhere [1,38] and is based on tree-by-tree records of long-term diameter growth over multiple intervals that total on average greater than 10 years per plot, and applies mean measured species-level wood densities, stand-level allometric models and census-interval corrections to estimate stand-level wood production.

(b) Model, modelling approach and parametrization

The canopy-scale photosynthesis model [10,39] uses the C3 leaf photosynthesis model from Farquhar & von Caemmerer [19] and the sun and shade approach from de Pury & Farquhar [40] to scale from leaf- to canopy-level, and assumes an optimization of stomatal behaviour [41]. The model has been calibrated and evaluated at five eddy covariance sites in the Brazilian Amazon.

Photosynthetic parameters (Vcmax and Jmax) were parametrized using two approaches. First, using eqn (2) from Domingues et al. [18] which includes both N and P constraints on photosynthesis and a leaf structure term S, specific leaf area, in (cm2 g−1). These ‘min{N : P}’ relationships are:Embedded Image 2.2where Embedded Image and Embedded Image are Vmax and Jmax expressed on a leaf dry-weight basis (µmol g−1 s−1), subsequently converted to area basis using the specific leaf area S, with [N]DW, [P]DW expressed in mg g−1 DW. The empirical coefficients from equations (2.2) areEmbedded Image andEmbedded Image with a coefficients in (µmol mg−1 s−1), b coefficients in (µmol cm−2 s−1) and c coefficients in (µmol g−1 s−1).

The remaining parameter values, from the photosynthesis model, are assumed invariant across all sites and are as follows: temperature sensitivities of Vmax and Jmax, SJ (693.124 J mol−1 K−1) and HJ (220 000 J mol−1), respectively; the curvature factor (0.7 unit less) of the potential rate of electron transport equation, quantum yield of photosynthesis (0.35 mol electrons mol−1 photons) and λ for the stomatal conductance model being 1200 mol mol−1.

(i) Scaling up from leaf- to canopy-level

Canopy-level values of Vmax, Jmax and leaf respiration Rd were estimated as the integral of the vertical profile of their leaf-level values over the entire canopy leaf area (L) following the method of de Pury & Farquhar [40], as implemented and applied in Mercado et al. [39] and Mercado et al. [10]. Canopy-level maximum carboxylation activity of RuBisCO Embedded Image (µmol m−2 s−1), is estimated asEmbedded Image 2.3where Embedded Image is the estimated RuBisCO activity at the top of the canopy obtained by the relevant substitutions of and Embedded Image, Embedded Image and/or Embedded Image into equation (2.2) and kP is the extinction coefficient for photosynthetic capacity that defines how photosynthetic capacity decreases with the cumulative L downwards from the top of the canopy. Low values of kP mean shallow profiles of photosynthetic capacity, and thus for any given values for Embedded Image and L then Embedded Image increases with decreasing kP.

The upper canopy photosynthetic capacity parameter estimates, Embedded Image and Embedded Image, were estimated using equation (2.2) with Rd being estimated as a constant fraction of Embedded Image (0.022) [10]. We calculate the extinction coefficient kP using the equation derived by Lloyd et al. ([28], fig. 10) from a compilation of data from different sources from broadleaf forest and trees, viz:Embedded Image 2.4where Embedded Image is in the units of µmol m−2 s−1.

From equation (2.4), ‘kP’ increases with increasing Embedded Image. This means that plants with high upper-leaf photosynthetic capacities have relatively steeper vertical canopy photosynthetic profiles than plants with low values of upper-leaf photosynthetic capacities which have shallower profiles.

We incorporate an inhibition of leaf respiration with light [42] where for any leaf in the canopy, Rd is reduced by 30 per cent when its incident irradiance, I, is higher than 10 µmol m−2 s−1.

(ii) Outline of simulations and analysis

  • — Canopy photosynthesis is simulated using the parametrizations given by the ‘min{N : P}’ equations and model evaluation is assessed with available observations from single components of the ecosystem carbon balance (the so-called ‘bottom up’) and also from eddy correlation.

  • — Additionally, the following sensitivity tests are performed. (i) A control simulation at which Embedded Image and Embedded Image are held constant (to the mean value for all sites given by the ‘min{N : P}’ equations), in order to assess the climate effect on simulated GP when compared with simulations with variable parameters, (ii) a set of simulations assuming N limitation, and (iii) P limitation at all sites using the ‘min{N : P}’ equations.

  • — Relationships between simulated GP (under all model configurations), observed leaf nutrients (N,P), total soil P and observed stem wood production (WP) and basal area growth (ΔB) are provided (at 33 sites where all data are available) and linear-adjusted Pearson correlation coefficient (R2) is used to assess how the above given variables are related.

3. Results

(a) Leaf nutrients

As shown in figure 2, Embedded Image varies by a factor of 3–4 when expressed on either an area or dry-weight basis and with Embedded Image varying less across the 38 sites. The highest Embedded Image content was found at some of the western sites in Ecuador and Peru while the lowest values are found in eastern Brazil. Lowest Embedded Image values were at some of the sites in northeast Venezuela and were highest in Ecuador, north and south Peru. As Embedded Image increases, Embedded Image : Embedded Image ratios decrease from ca 40 to 10.

Figure 2.

Variability of foliar Embedded Image and Embedded Image on (a) area basis, (b) dry-weight basis and (c) their N : P ratio. Data taken from Fyllas et al. [7].

(b) Relationship between foliar (N,P) and soil (P) nutrients to observed stem wood production (WP) and basal area growth (ΔB)

Figure 3 illustrates relationships between measured Embedded Image and Embedded Image in leaves on a dry-weight and area basis with WP (top row) and ΔB (bottom row). Additionally, the figure shows the relationships between measured total soil P with WP (top row) and ΔB (bottom row). Pearson's adjusted R2 shows significant correlation (p < 0.05) between foliar P, on both dry-weight and area basis with WP and ΔB. Total soil P is also significantly correlated with ΔB, but apparently not with WP (p = 0.062). Foliar N is not correlated with WP, but is correlated with ΔB, but only when expressed on an area basis. This suggests that P is more important than N at explaining some of the observed variability in stem-wood production and basal area growth. However, foliar P and foliar N are correlated with each other (R = 0.54, p < 0.001). Using partial correlation, foliar P is significantly correlated with WP (R =+0.49, p = 0.004) controlling for foliar N. On the other hand, partial correlation of foliar N with WP is not significant (R = −0.06, p = 0.737) after controlling for its correlation with foliar P. Using a canopy photosynthesis model, we now explore the extent to which P limitation to photosynthesis is more important than N limitation as implied by figure 3.

Figure 3.

Relationships of foliar N and P (on dry-weight and area basis) [7] and soil P content [5] against (a) observed stem growth (WP), (b) basal area growth (ΔB) [6] at 33 sites. Filled and open symbols correspond to N- and P-limited sites, respectively, according to equation (2.2). The triangle symbols correspond to sites with no available soil phosphorous data.

(c) Relationship between simulated net carbon uptake and observed stem wood production (WP) and basal area growth (ΔB)

Simulated mean annual Embedded Image, averaged over the period 1980–2001 for all sites, as shown in table 2, ranges from 30.9 Mg C ha−1 a−1 for CAX-01 in eastern Amazonia to 41.2 Mg C ha−1 a−1 for JAS-05 in Ecuador.

View this table:
Table 2.

Simulated mean annual Embedded Image in Mg C ha−1 a−1 for the period 1980–2001 using ‘min {N : P}’ (equation 2.2) and ‘PA-only’ (equation A 1) model parametrizations.

In figure 4, the relationships between simulated GP and observed WP (top row) and ΔB (bottom row) are illustrated for four model configurations with various combinations of Embedded Image and Embedded Image. For the first simulation, Embedded Image and Embedded Image are taken as invariant across all sites using the mean dataset values obtained from equation (2.2). This simulation thus shows the effect of climate alone on simulated GP. The remaining three simulations allow for between-site climate variability, but with different assumptions regarding the nature of nutrient limitations on tropical forest tree photosynthesis viz: (i) N limitation, (ii) P limitation, and (iii) both N and P limitation with Embedded Image and Embedded Image calculated using the ‘min{N : P}’ relationship of equation (2.2) (figure 4).

Figure 4.

Relationships of simulated Embedded Image (using the ‘min{N : P}’ relationship under four model configurations) against (a) observed stem growth (WP), and (b) basal area growth (ΔB) [6] at 33 sites. Each of the four model configurations correspond to a column: (i) with invariant Embedded Image and Embedded Image, assuming (ii) N limitation, (iii) P limitation and (iv) both N and P limitation across all studied sites. Filled and open symbols correspond to N and P limited sites, respectively, according to equation (2.2). The triangle symbols correspond to sites with no available soil phosphorous data.

Simulated GP with invariant photosynthetic parameters showed no significant correlation with either WP or ΔB. There was, however, a significant correlation between GP simulated under N-limitation, under P-limitation, and under both N- and P-limitation with both WP and ΔB. From the relationships with WP, simulated GP under P-limitation (adjusted R2 = 27.4%) explains nearly the same observed variability as when both N and P are limiting (adjusted R2 = 26.25%). From the relationships with ΔB, simulated GP under P-limitation (adjusted R2 = 48.2%) explains less observed variability than when both N and P are limiting (adjusted R2 = 60.5%). This means that, when explaining basal area growth, N becomes important as well. This is consistent with figure 3, based on observations; it is clear that P is more important than N in accounting for variations in canopy photosynthesis that are relatable to the observed variability in stem wood production. However, in some circumstances to explain the basal area growth rates, a relatively low N availability is important.

(d) Model evaluation

Comparison of observed and modelled Embedded Image, GP and RC at the few available sites (table 3 and figure 5) shows that simulated values are close to the observations, with simulated RC closer to the bottom-up estimates (−5%). Comparison between Embedded Image derived from bottom-up and eddy correlation shows much more variability in the bottom-up values. Simulated Embedded Image using the ‘min{N : P}’ relationships is on average 14 and 7 per cent higher than the mean from bottom-up and eddy covariance estimates, respectively.

View this table:
Table 3.

Model evaluation of simulated gross primary productivity in the absence of foliar respiration (Embedded Image defined as GP plus canopy foliar respiration, RC), GP and RC in Mg C ha−1 a−1 using available observations from ‘bottom-up’ and eddy correlation. Upper and lower bounds of uncertainty correction for eddy correlation values are related to a night-time friction velocity threshold used [35]. Data sources for each site are described in the methods.

Figure 5.

Model evaluation of simulated gross primary productivity (Embedded Image), and leaf respiration (RC) using available observations derived from (a) eddy correlation and (b,c) bottom-up. Error bars for eddy correlation measurements correspond to uncertainty [35] and correspond to standard error for bottom-up estimates. Data sources for each site are described in the methods.

4. Discussion

(a) Model parametrization and implications for simulated GP

Although the parametrization of Domingues et al. [18] was developed using a dataset based on savannah and forest trees from West Africa, it performed surprisingly well, in as much that the predicted GP correlated well with the observed values of WP and ΔB across the Amazon Basin. Although it might be argued that this was more-or-less inevitable given the already strong correlation between foliar P and these growth parameters (figure 3), what is important here is that this result was obtained through the model results showing that most sites should be represented as phosphorus- rather than nitrogen-limited in order to account for the maximum range in woody growth rates (figure 4). Only a few sites were predicted to be nitrogen-limited (black in figure 1, filled circles in figures 3 and 4). These were located on either alisols, fluvisols or cambisols—all of which are relatively high in available soil phosphorus [5] or, in one case, on a podzol, this being a soil type for which nitrogen availability may be unusually low [5,12]. Nevertheless, despite some suggestion of N limitations (at seven sites), there was little practical difference between the ‘P-limited’ sites and ‘min{N : P}’ cases (figure 4, third and fourth columns) in their ability to predict either of the growth metrics investigated. This is because even when N-limitation was predicted, the modelled photosynthetic parameters under P-limitation were only marginally greater. For example at TIP-03, BOG-02 and YAN-02, it appears that N and P are co-limiting and Embedded Image derived from Embedded Image and Embedded Image were 54.5 and 55.2 µmol m−2 s−1 at TIP03, 56.2 and 59.8 µmol m−2 s−1 at BOG-02 and 50.3 and 50.6 µmol m−2 s−1 at YAN-02, respectively. At the remaining four sites, more N-limitation was obtained (only three included in R2 calculations), corresponding values for Embedded Image derived from Embedded Image and Embedded Image were 59.7 and 69.1 µmol m−2 s−1 at JAS-05, 57.8 and 65.3 µmol m−2 s−1 at TAM-06, 50.3 and 56.4 µmol m−2 s−1 at CUZ-03 and 52.4 and 65.3 µmol m−2 s−1 at SCR-04. Finally, application of the Domingues et al. [18] parametrization to our dataset showed P-limitation at sites with foliar N : P ratios higher than 17. This is also the threshold value suggested for P-limitation in plants [43,44].

Indeed, it may be that because the small number of N-limited sites are used in the simple (but independent) parametrization of photosynthesis based on foliar [P] in area basis, the ‘[P]A only’ equation of Mercado et al. [10] (equation A 1) gives as good relationship between GP and both WP and ΔB (figure 6) as the more complex ‘min{N : P}’ model of Domingues et al. [18]. Additionally, the related model predictions of GP based on foliar [N] on area basis ‘[N]A only’ (equation A 1) performs particularly poorly (figure 6). Indeed, the simple ‘[P]A only’ equation actually gives rise to slightly higher R2 than the ‘min{N : P}’ model. The reasons for this are unclear, but might have to do with the fact that the ‘[P]A only’ relationship was empirically derived by fitting the canopy-scale model used in this study to eddy correlation data from just five sites, all located in poor-nutrient soils. Additionally, even though Vmax parametrizations (equations (2.2) and A 1) were quantitatively similar at many sites (not shown), simulated GP was higher with the ‘[P]A only relationship’ (table 2). This is mostly because the ratio Jmax/Vmax was higher in those simulations (1.92 versus Jmax/Vmax ranging between 1.5 and 1.6), which translates into higher light-limited velocity of photosynthesis, and therefore simulated GP under both sets of parametrizations is not strictly comparable.

Figure 6.

Relationships of simulated Embedded Image against (a) observed stem growth (WP) and (b) basal area growth (ΔB) [6] at 33 sites. Model configurations correspond to assumption of (i) N limitation (using [N]A only) and (ii) P limitation (using [P]A only) across all studied sites. Filled and open symbols correspond to N- and P-limited sites, respectively, according to equation (2.2). The triangle symbols correspond to sites with no available soil phosphorous data.

The ‘min{N : P}’ parametrization was evaluated against single leaf Vcmax and Jmax derived from gas exchange measurements at a site in Tapajos [18]. Results show no statistically significant difference between the mean predicted and observed values (mean values of predicted and observed Vcmax were 0.45 ± 0.09 and 0.39 ± 0.16 µmol g−2 s−1 and mean values of predicted and observed Jmax were 0.69 ± 0.15 and 0.60 ± 0.22 µmol g−2 s−1). Furthermore, estimates of Vcmax derived from gas exchange measurements for tropical rainforest in the Amazon are scarce. Here, we use the only five available data sources that we are aware of to compare against the top of the canopy Vcmax estimated in this study. Our estimates are close to, and within the range of, mean values inferred from gas exchange measurements at MAN-05 [45], CAX-06 [46] TAP-04 in Brazil [47] and LFB-02 in Bolivia (T. Domingues 2007, unpublished data), as shown in figure 7. Estimated Jmax/Vmax in this study ranged between 1.5 and 1.6 using the ‘min{N : P}’ relationship, which agrees well with published values for other tropical rainforest in the Amazon and in Africa. Jmax/Vmax values for upper-canopy leaves at Manaus [45], Caxiuana [46], Tapajos [47] and LSL-02 (T. Domingues 2007, unpublished data) have been reported as 2.27, 1.94, 1.59 and 1.92, respectively. Also, Coste et al. [48] reported a Jmax/Vmax ratio of 2.1 from seedlings of 14 tree species in the tropical forest of French Guiana. Meir et al. [49] obtained Jmax/Vmax ratios of 1.7 from gas exchange measurements at a rainforest in Cameroon with broadleaf, coniferous, shrubs and herbaceous plants.

Figure 7.

Comparison of top of the canopy Vcmax derived from gas exchange measurements at four rainforest sties in the Amazon Basin (crosses, error bars correspond to one standard deviation), and estimated in this study using the ‘min{N : P}’ relationship (circles). Data sources from Brazilian sites are Caxiuana from Vale et al. [46], Manaus from Carswell et al. [45] and Tapajos from Domingues et al. [47]. The Bolivian site LFB-02 is courtesy of T. F. Domingues, unpublished data.

(b) Model evaluation

The model evaluation presented in this study comparing simulated Embedded Image against bottom-up and eddy covariance flux estimates is intended as a point of reference for the model only (i.e. as a general evaluation). This is because simulations correspond to an average over the period 1980–2001 and the available Embedded Image estimates from both type of observations used in this study correspond to a different and shorter period of time, generally a few years at most.

There are errors associated with both methods used to estimate Embedded Image, with both likely to underestimate the fluxes. Underestimation of total net ecosystem exchange using bottom-up approaches for tropical forest can potentially be up to 20 per cent [50]. There are also difficulties in estimating Embedded Image from eddy covariance fluxes. Important sources of error from eddy covariance measurements include data representativity, treatment of data gaps, flux correction for systematic errors and the non-measured night-time fluxes [51]. The latter are associated with a frequent failure of the system to measure night-time fluxes during low wind speed conditions [34,5254]. Therefore, only observations above a defined wind-speed threshold can be used for the purpose of estimating ecosystem respiration, Re which is used together with the daytime net ecosystem exchange measurements to estimate Embedded Image from eddy covariance. Unmeasured night-time fluxes can introduce significant errors into estimates of Re [51]. Additionally, the response of Re to daytime temperatures and a response of leaf respiration to both temperature and light conditions [42] are still unaccounted for in the calculation of Re in the studied sites. However, although these effects might counteract each other, there are no available measurements to quantify them. As for daytime measurements, the total uncertainty estimated for the Manaus K34 and Jaru site are ±12 and ±32 per cent, respectively [55].

Once the difficulties in estimating Embedded Image from bottom-up and eddy covariance measurements are considered, we conclude that our model estimates are in reasonable agreement. Furthermore, the canopy exchange photosynthesis model used in this study showed agreement of between ±10 per cent with eddy covariance estimates of Embedded Image at the five main eddy covariance sites (MAN-C14, MAN-K34, TAP-K67, Jarú and CAX-06, [10]), despite the photosynthesis/nutrient parametrizations. Similarly, our estimates of Embedded Image agree with values reported for tropical humid evergreens in the global dataset of Luyssaert et al. [56] with estimates from the bottom-up approach (36 Mg C ha−1 a−1, from six sites). Nevertheless, the main purpose of this study was to simulate variation in GP across the Amazon Basin and to look for relationships with observed stem wood production (WP) and basal area growth (ΔB) rather than to simulate the absolute value of GP at any single site.

Evaluation of simulated leaf respiration (figure 5 and table 3) shows good model agreement against observations. Analyses with global datasets have shown a consistent coupling between mass-based respiration rates and [N]DW content across species [57]. In this study, we estimate leaf respiration rates at 25°C from Vcmax, but Vcmax is usually constrained by leaf P, according to our models. Although numerically our leaf respiration rates seem to match the scaled-up observations, it remains to be determined if this assumption is valid. Nevertheless, it is worth pointing out that Meir et al. [58] found high correlation between measured leaf respiration rates and foliar [P]A, especially when specific leaf area was included in their analysis for the Jaru tower site.

The modelling framework used in this study does not consider plant-related soil water stress. The stomatal conductance formulation used [41] takes into account drought stress via atmospheric vapour-pressure deficits. However, especially at forest sites with strong seasonality, there is a strong stomatal control of transpiration combined with deep rooting systems that allow water recharge from deep soils [59]. There is little evidence of soil water deficit effects on tropical forest productivity [32,59,60], except for the most peripheral forests close to the forest/savannah transition zone [9] and no such forests were considered in this study.

(c) Relating simulated net carbon uptake to measured stem wood production (WP) and basal area growth (ΔB)

In the analysis presented, we have shown that simulated GP using ‘min{N : P}’ can explain approximately 30 per cent of observed variability in observed wood production including P-limitation alone or when including both N- and P-limitation. Including N-limitation with P-limitation marginally decreases the obtained correlation, and simulated GP under N-limitation alone can only explain 15 per cent of the observed variability in WP. However, from the seven N-limited sites (only six included in calculations of R2), three presented N and P co-limitation. Additionally, our analysis suggests that variability in local meteorological variables given fixed nutrient levels at the various sites does not simulate variation in GP sufficient to explain any of the observed variability of WP and ΔB. Finally, using a species abundance weighting (equation (2.1)) to estimate leaf nutrients at each plot proved to have an impact on simulated Gp and its capability to explain variability of WP and ΔB (adjusted R2 of 26%, p = 0.001 and 60%, p < 0.001 respectively) as opposed to using plot level means which gave lower correlations and significance levels (adjusted R2 of 15%, p = 0.016 and 23%, p = 0.003 respectively, not shown).

Our simulations gave only a limited variability in GP for the 38 sites simulated (with 95% of GP estimates within 10% of our overall basin-wide means). These small variations in simulated GP explained up to 26 and 60 per cent of the observed variability in WP and ΔB, respectively. This suggests not only that variations in WP for the forests examined may indeed be driven to a considerable extent by variations in GP (themselves driven by variations in canopy nutrient concentrations, and especially by variations in phosphorus), but also that as GP increases beyond a certain point, then much of the extra carbohydrate availability is preferentially allocated towards WP. For example, at a low GP of 25 Mg C ha−1 a−1 as simulated by the ‘min{N : P}’ model, WP is estimated at 3.9 Mg C ha−1 a−1 equal to just 0.16 of GP. But at a higher GP of 35 Mg C ha−1 a−1, WP is estimated at 11.3 Mg C ha−1 a−1, equal to 0.32 of GP. This is as suggested by Lloyd & Farquhar [61] and has important implications for any stimulation of GP as might be expected to occur, for example, in response to elevated [CO2] with any stimulation of WP being proportionally much higher than any associated increase in GP.

For ΔB, the differences are even more profound with the ratio ΔB/GP varying more than threefold from 0.012 to 0.037 m2 Mg−1C for the same simulated GP variations. This greater relative sensitivity is due to a consistent decline in tree-wood density with increasing GP. Indeed, based on data presented by Baker et al. [38,62], we have also estimated changes in stand-level wood density (Embedded Image) and found even stronger relationships with GP than for either WP or ΔB with R2 of 0.76 and a slope of around 40 mol µmol kg m−1 s−1 (results not shown). For example, Embedded Image at a low GP of 25 Mg C ha−1 a−1 as simulated by the ‘min{N : P}’ model is estimated at 730 kg m−3, but it is also estimated that Embedded Image would be only 331 kg m−3 for GP of 35 Mg C ha−1 a−1. Thus, higher productivity sites are also characterized by a greater abundance of species characterized by a low Embedded Image capable of high rates of basal area (and presumably height) growth.

Variability in carbon use efficiency (CUE), the ratio of net primary productivity to Embedded Image, across the studied transect owing to variability in autotrophic respiration rates has been proposed as a possible explanation for the observed spatial variability in observed rates of stem wood productivity [1]. This hypothesis suggests lower rates of CUE owing to higher plant respiration at low productivity sites and higher rates of CUE and lower plant respiration rates at the high productivity sites. Ongoing fieldwork campaigns at various sites across the Amazon Basin are currently measuring the individual carbon cycle components. This type of research should lead to a better understanding of (i) variability of main fluxes across the Amazon Basin and, most importantly, (ii) the mechanism behind a possible variability in plant respiration as proposed, and its relationship to nutrient supply, which as shown in the present analysis has a major role in influencing observed spatial variability in stem productivity across the Amazon Basin.

Global vegetation models, which ignore phosphorus availability [63,64], would not predict at least some of the observed spatial variability in stem growth-rates that clearly occurs across the Amazon Basin. This variability appears to be of fundamental ecological importance, because it is closely correlated in space with similar variation in carbon stores (biomass), stand-level species composition, wood density and population dynamics, and with underlying soil properties. The likelihood, therefore, is that these macro-ecological gradients will prove important in determining forest responses to global change drivers such as increasing [CO2] and temperatures, as well as determining how forest ecosystems may respond to other more local drivers. For example, it has been suggested that western Amazon forests are less sensitive than eastern forests to some forms of anthropogenic disturbance, because they are adapted to the naturally much higher turnover rates of these forests compared with the eastern part of the Basin [65,66]. Clearly, global-scale models need to be able to reproduce the present-day spatial variation in tree growth-rates across the Amazon and to include any relevant additional process parametrizations, in our specific case effects of soil fertility on foliar [P] and hence on simulated net carbon uptake. Such advances should lead to an enhanced ability of global models to simulate the response of tropical ecosystems to future changes in climate and atmospheric composition, and to improved quantification of future climate–carbon cycle feedbacks.

5. Summary

We have simulated the sensitivity of GP to parametrizations of photosynthetic capacity, in order to understand spatial dynamics of observed WP. We tested model sensitivity to including nutrient (N and P, only N, and only P) constraints on photosynthesis, and how the simulated photosynthesis relates to variables such as above-ground growth and biomass. In agreement with previous studies showing soil phosphorous as the principal soil fertility factor influencing rates of tree growth across the Amazon Basin, simulations performed in this study suggest that when including a leaf [P] constraint to photosynthetic capacity, simulated GP agrees with the few available observations and can explain up to approximately 30 per cent of the spatial variability in stem growth. Incorporating this improved understanding of the role of soil nutrients in explaining basin-wide variations in GP will lead to improved parametrizations in global carbon cycle models. This should enable improved simulation of the response of tropical ecosystems to future changes in climate and atmospheric composition, and improved quantification of future climate–carbon cycle feedbacks.


Sandra Patiño and Samuel Almeida, two of our co-authors, have left us during the process of the submission and publication of this paper. We would like to dedicate this piece of work to both of them. Financial support for L.M. was provided through the EU Project ‘Carbonsink-LBA Proposal no. EVK 1999–0191’ with most of the modelling work being done at the Max Planck Institute for Biogeochemistry, Jena, Germany. The final analysis was also supported through the UK NERC Tropical Biomes in Transition (TROBIT) consortium, UK Natural Environment Research Council (NERC) grant (NE/A/S/2003/00608/2) and the UK NERC Amazon Integrated Carbon Analysis—(AMAZONICA) consortium grant (NE/F005997/1). O.P. was partially funded through the European Union-funded PANAMAZONIA programme. N.M.F. is supported by a Marie Curie Intra-European Fellowship within the 7th European Community Framework Programme. G.P.W. was supported by the Joint DECC/Defra Met Office Hadley Centre Climate Programme—(GA01101). We also thank Gabriela López-Gonzalez for developing the forest productivity and biomass database, with funding from the Gordon and Betty Moore Foundation (grant to RAINFOR) and the School of Geography, University of Leeds. The RAINFOR programme is supported by the Gordon and Betty Moore Foundation. RAINFOR field campaigns contributing to this study were funded by the Royal Geographical Society, CARBONSINK-LBA and the Max Planck Institute for Biogeochemistry, Jena, Germany. We gratefully acknowledge the support and funding of organizations who have contributed to the establishment and maintenance of individual sites: Natural Environment Research Council, EU Framework V and VI, US National Geographical Society, WWF-US/Garden Club of America, US National Science Foundation, the Nature Conservancy/Mellon Foundation (Ecosystem Function Programme), Conselho Nacional de Desenvolvimento Cientifico e Tecnológico, Museu Goeldi, Estacao Cientifica Ferreira Penna, Andrew W. Mellon Foundation, NASA-LBA Programme, Conservation, Food and Health Foundation, MacArthur Foundation, Fundación Jatun Sacha, Estación Científica Yasuní de la Pontificia Universidad Católica del Ecuador, Estación de Biodiversidad Tiputini, Conservation International, ACEER, Albergue Inkaterra, Explorama Tours S.A., Explorers Inn, IIAP, INRENA, UNAP and UNSAAC. We thank two anonymous reviewers for comments on the manuscript. This paper constitutes Publication No. A/582 of the Royal Society South East Asia Rainforest Research Programme.

Appendix A

In addition to the min{N : P} equation, we investigated the utility of the Vmax versus leaf P (area basis) relationship: the ‘[P]A only’ relationship as obtained in Mercado et al. [10] with the ‘A’ subscript here used to denote values expressed on a leaf area basis. This relationship was obtained by relating the best-fitted top of the canopy Vmax from each of five calibration sites, against measured upper canopy [N]A and [P]A for the same sites. The relationship for [N]A (‘[N]A only’) was very poor (R2 = 0.08), but much better for [P]A (R2 = 0.41). Jmax has been estimated as a constant ratio (r) of Vmax as calibrated in the modelling exercise of Mercado et al. [10],Embedded Image A 1with corresponding values for a, b and r of 386.9 µmol m−2 s−1, 17.9 µmol g−1 s−1 and 1.92, respectively. Respective values for the slope, intercept and r for the ‘[N]A only’ relationship are 6.8 µmol m−2 s−1, 30.8 17.9 µmol g−1 s−1 and 1.92.



View Abstract