Simulation of interannual responses of trembling aspen stands to climatic variation and insect defoliation in western Canada

A carbon-based model has been developed to simulate responses of trembling aspen (Populus tremuloides Michx.) stands to interannual climatic variation and insect defoliation. The model is designed for medium time scale (10–100 years) simulations and requires only daily maximum and minimum temperature and precipitation as meteorological inputs. The modelling approach is similar to FOREST-BGC but includes additional processes known to be important in deciduous forests. These include removal of leaf area during outbreaks of forest tent caterpillar (Malacosoma disstria Hbn.), phenological changes in leaf area index, storage and allocation of non-structural carbohydrate and the contribution of understorey vegetation to evapotranspiration. The model was used for simulations of growth and mortality of biomass carbon in two mature aspen forests located in the climatically dry transition zone between the boreal forest and prairie grassland regions of Saskatchewan, Canada. Model inputs of annual defoliation intensity were based on historic records of insect defoliation and the incidence of light-coloured tree rings in disks or cores collected from aspen at each of the two sites. At both sites, moderately good correlations (r²=0.47–0.54) were obtained between modelled interannual changes in stem carbon growth and observed interannual changes in stem basal area increment obtained from tree-ring analysis. Model outputs of stem biomass carbon were found to be highly sensitive to parameters describing seasonal leaf area duration, insect defoliation intensity, photosynthesis and root respiration and carbohydrate allocation to growth versus storage.     

Constitutive or induced defences - how does <Emphasis Type="Italic">Eucalyptus globulus</Emphasis> defend itself from larval feeding?

Summary. Following herbivory, induced responses involving plant secondary metabolites have been reported in a number of tree species. Although a wide range of plant secondary metabolites appear to operate as constitutive plant defences in trees belonging to the Eucalyptus genus, no induced responses have as yet been reported following foliar-chewing insect damage. We empirically tested whether branch defoliation (artificial and larval) of 2-year-old Eucalyptus globulus Labill. trees altered the abundance of specific plant secondary metabolites immediately (3 months after initial larval feeding) and 8 months after the cessation of larval feeding. Metabolites assayed, included essential oils, polyphenolic groups and foliar wax compounds and in all cases their abundance was not significantly altered by defoliation. However, the level of foliar tannins after 3 months of larval feeding did display a trend that suggested elevated levels as the result of defoliation, though this trend was not evident 8 months later, indicating that, if real, the response was a rapid and not a delayed induced response. The level of foliar tannins was also negatively correlated to both average larval survival and average percentage branch defoliation, suggesting that foliar tannins may operate as toxins and/or anti-feedants to M. privata larval feeding.     

Susceptibility of Eucalyptus globulus to Mnesampela privata defoliation in relation to a specific foliar wax compound

Summary. Increased abundance of several aliphatic benzyl and phenylethyl alkanoates were previously associated with reduced defoliation of Eucalyptus globulus Labill. by the autumn gum moth (Mnesampela privata Guenée) in two clonally replicated F₂ families. We examine the robustness of this association by comparing the abundance of the implicated foliar wax compounds with defoliation levels in three plantation sites and a genetics trial in northern Tasmania. At all locations, the aliphatic alkanoate, benzyl n-tetracosanoate (C₂₄), was found in significantly higher abundance in the foliar wax of E. globulus trees that had displayed low M. privata damage (designated resistant) compared to trees that had displayed high damage (designated susceptible). Further, when sites were combined, benzyl n-tetracosanoate (C₂₄) accounted for 16 % of variation in field defoliation. Field defoliation was also positively correlated to M. privata oviposition both in the field and on foliage sprigs in a cage. In the cage bioassay, oviposition accounted for 20 % of variation in field defoliation making it a better predictor of tree defoliation in the field than either field oviposition, which accounted for 9 % of variation, or screening foliage for the abundance of benzyl n-tetracosanoate (C₂₄). Despite both benzyl n-tetracosanoate (C₂₄) and oviposition being related to field defoliation, there was no statistically significant relationship between the abundance of benzyl n-tetracosanoate (C₂₄) and oviposition in the field or cage bioassay. Further work is therefore required to determine if increased levels of this compound in the foliar wax operates as a direct deterrent to M. privataoviposition.     

Mechanisms of moisture stress in a mid-latitude temperate forest: Implications for feedforward and feedback controls from an irrigation experiment

While it is well established that stomata close during moisture stress, strong correlations among environmental (e.g., vapor pressure deficit, soil moisture, air temperature, radiation) and internal (e.g., leaf water potential, sap flow, root-shoot signaling) variables obscure the identification of causal mechanisms from field experiments. Models of stomatal control fitted to field data therefore suffer from ambiguous parameter identification, with multiple acceptable (i.e., nearly optimal) model structures emphasizing different moisture status indicators and different processes. In an effort to minimize these correlations and improve parameter and process identification, we conducted an irrigation experiment on red maples (Acer rubrum L.) at Harvard Forest (summers of 2005 and 2006). Control and irrigated trees experienced similar radiative and boundary layer forcings, but different soil moisture status, and thus presumably different diurnal cycles of internal leaf water potential. Measured soil moisture and atmospheric forcing were used to drive a transient tree hydraulic model that incorporated a Jarvis-type leaf conductance in a Penman–Monteith framework with a Cowan-type (resistance and capacitance) tree hydraulic representation. The leaf conductance model included dependence on both leaf matric potential, Ψ_L (so-called feedback control) and on vapor pressure deficit, D (so-called feedforward control). Model parameters were estimated by minimizing the error between predicted and measured sap flow. The whole-tree irrigation treatment had the effect of elevating measured transpiration during summer dry-downs, demonstrating the limiting effect that subsurface resistance may have on transpiration during these times of moisture stress. From the best fitted model, we infer that during dry downs, moisture stress manifests itself in an increase of soil resistance with a resulting decrease in Ψ_L, leading to both feedforward and feedback controls in the control trees, but only feedforward control for the irrigated set. Increases in the sum-of-squares error when individual model components were disabled allow us to reject the following three null hypotheses: (1) the f(D) stress is statistically insignificant (p = 0.01); (2) the f(Ψ_L) stress is statistically insignificant (p = 0.07); and (3) plant storage capacitance is independent of moisture status (p = 0.07).     

Inter-annual variability of ecosystem production in boreal jack pine forests (1975-2004) estimated from tree-ring data using CBM-CFS3

We describe and apply a method of using tree-ring data and an ecosystem model to reconstruct past annual rates of ecosystem production. Annual data on merchantable wood volume increment and mortality obtained by dendrochronological stand reconstruction were used as input to the Carbon Budget Model of the Canadian Forest Sector (CBM-CFS3) to estimate net ecosystem production (NEP), net primary production (NPP), and heterotrophic respiration (Rh) annually from 1975 to 2004 at 10 boreal jack pine (Pinus banksiana Lamb.) stands in Saskatchewan and Manitoba, Canada. From 1975 (when sites aged 41-60 years) to 2004 (when they aged 70-89 years), all sites were moderate C sinks except during some warmer than average years where estimated Rh increased. Across all sites and years, estimated annual NEP averaged 57 g Cm⁻² yr⁻¹ (range −31 to 176 g Cm⁻² yr⁻¹), NPP 244 g Cm⁻² yr⁻¹ (147-376 g Cm⁻² yr⁻¹), and Rh 187 g Cm⁻² yr⁻¹ (124-270 g Cm⁻² yr⁻¹). Across all sites, NPP was related to stand age and density, which are proxies for successional changes in leaf area. Regionally, warm spring temperature increased NPP and defoliation by jack pine budworm 1 year previously reduced NPP. Our estimates of NPP, Rh, and NEP were plausible when compared to regional eddy covariance and carbon stock measurements. Inter-annual variability in ecosystem productivity contributes uncertainty to inventory-based assessments of regional forest C budgets that use yield curves predicting averaged growth over time. Our method could expand the spatial and temporal coverage of annual forest productivity estimates, providing additional data for the development of empirical models accounting for factors not presently considered by these models.     

Interpreting spatial heterogeneity of crop yield with a process model and remote sensing

A process-based crop growth model (Vegetation Interface Processes (VIP) model) is used to estimate crop yield with remote sensing over the North China Plain. Spatial pattern of the key parameter—maximum catalytic capacity of Rubisco (V_cmax) for assimilation is retrieved from Normalized Difference of Vegetation Index (NDVI) from Terra-MODIS and statistical yield records. The regional simulation shows that the agreements between the simulated winter wheat yields and census data at county-level are quite well with R² being 0.41-0.50 during 2001-2005. Spatial variability of photosynthetic capacity and yield in irrigated regions depend greatly on nitrogen input. Due to the heavy soil salinity, the photosynthetic capacity and yield in coastal region is less than 50 μmol C m⁻² s⁻¹ and 3000 kg ha⁻¹, respectively, which are much lower than that in non-salinized region, 84.5 μmol C m⁻² s⁻¹ and 5700 kg ha⁻¹. The predicted yield for irrigated wheat ranges from 4000 to 7800 kg ha⁻¹, which is significantly larger than that of rainfed, 1500-3000 kg ha⁻¹. According to the path coefficient analysis, nitrogen significantly affects yield, by which water exerts noticeably indirect influences on yield. The effect of water on yield is regulated, to a certain extent, by crop photosynthetic capacity and nitrogen application. It is believed that photosynthetic parameters retrieved from remote sensing are reliable for regional production prediction with a process-based model.     

浙江千岛湖地区苦槠叶片昆虫取食状类型和取食强度研究

选择千岛湖地区不同大小岛屿为样地,对苦槠(Castanopsis sclerophylla)叶片虫食状类型和取食强度进行了研究.结果表明,千岛湖地区苦槠叶片昆虫取食状类型为16种,但取食频率不同.取食频率和取食部位对叶片危害程度的分析结果显示,缘食状、掏食状、顶食状、阴面食状和叶中脉状对植物叶片的影响更大.16种取食状中,缘食状和顶食状、缘食状和阴面食状之间呈极显著正相关(P＜0.01),缘食状与大孔状、顶食状与大孔状、顶食状与掏食状、阳面食状与阴面食状之间呈显著正相关(P＜0.05).苦槠叶片被昆虫取食的频率为84%,但绝大多数受害程度低.     

Impacts of changing frost regimes on Swedish forests: Incorporating cold hardiness in a regional ecosystem model

Understanding the effects of climate change on boreal forests which hold about 7% of the global terrestrial biomass carbon is a major issue. An important mechanism in boreal tree species is acclimatization to seasonal variations in temperature (cold hardiness) to withstand low temperatures during winter. Temperature drops below the hardiness level may cause frost damage. Increased climate variability under global and regional warming might lead to more severe frost damage events, with consequences for tree individuals, populations and ecosystems. We assessed the potential future impacts of changing frost regimes on Norway spruce (Picea abies L. Karst.) in Sweden. A cold hardiness and frost damage model were incorporated within a dynamic ecosystem model, LPJ-GUESS. The frost tolerance of Norway spruce was calculated based on daily mean temperature fluctuations, corresponding to time and temperature dependent chemical reactions and cellular adjustments. The severity of frost damage was calculated as a growth-reducing factor when the minimum temperature was below the frost tolerance. The hardiness model was linked to the ecosystem model by reducing needle biomass and thereby growth according to the calculated severity of frost damage. A sensitivity analysis of the hardiness model revealed that the severity of frost events was significantly altered by variations in the hardening rate and dehardening rate during current climate conditions. The modelled occurrence and intensity of frost events was related to observed crown defoliation, indicating that 6-12% of the needle loss could be attributed to frost damage. When driving the combined ecosystem-hardiness model with future climate from a regional climate model (RCM), the results suggest a decreasing number and strength of extreme frost events particularly in northern Sweden and strongly increasing productivity for Norway spruce by the end of the 21st century as a result of longer growing seasons and increasing atmospheric CO₂ concentrations. However, according to the model, frost damage might decrease the potential productivity by as much as 25% early in the century.     

Modelling the spatial population dynamics of the green oak leaf roller (Tortrix viridana) using density dependent competitive interactions: Effects of herbivore mortality and varying host-plant quality

Cyclic population dynamics of forest insects with periods of more than two generations have been discussed in relation to a variety of extrinsic and intrinsic forces. In the present study, we employed the selection pressure of density dependent competitive interactions according to Witting's equations (Witting, 2000) as driver for a discrete spatiotemporal model of the green oak leaf roller (Tortrix viridana). The model was successfully parameterised to rebuild the cyclic population dynamics of an empirical data set of a 30-year leaf roller monitoring in Russia. Our analysis focussed on the role of herbivore mortality and host plant food quality, which have a significant effect on T. viridana population dynamics. An additional egg or larvae mortality lowers population density and can lead to selection pressures that favour individuals with higher growth rate. This increased population growth rate can not only compensate the additional mortality, but also can lead to higher average moth abundances in subsequent generations. Furthermore, we analysed the effect of inter- and intraspecific variation in host plant quality on herbivore population dynamics and the spatial distribution of abundance and defoliation patterns. We found significant effects of the qualitative composition of a trees neighbourhood on the herbivore population of the respective tree. Also, the patchy damage patterns observable in reality have been reproduced by the present model. The applicability of the model approach and the putative genetic processes underlying Witting's model are discussed.     

Relationships between simulated water stress and mortality and growth rates in underplanted Toona ciliata Roem. in subtropical Argentinean plantations

Toona ciliata Roem. (Australian red cedar) requires a nurse-tree overstory to prevent damage from drought and irradiation in some regions of north-eastern Argentina. T. ciliata was planted in the understory of Pinus taeda L. (625 stems/ha), Pinus elliottii Engelm. × Pinus caribaea Morelet (625 stems/ha), and Grevillea robusta A. Cunn. (833 stems/ha) nurse trees, which were thinned to 0, 25, 50, 75 and 100% of the initial densities. We measured initial T. ciliata mortality and growth as well as Leaf Area Index (LAI) based on light transmission. T. ciliata soil water availability and its effect on early growth and mortality were examined by modelling drought stress using the two-dimensional forest hydrology model ForWaDy. Simulated patterns in T. ciliata water stress for the different overstory treatments were consistent with observed patterns of mortality. Early mortality was lowest with a G. robusta overstory, with corresponding lowest drought stress values and high modelled soil water contents in the top soil layer in intermediate and high overstory densities. Mortality was highest with a P. elliottii × P. caribaea overstory in treatments with the highest modelled drought stress values in the most open treatments. The model supported our field observations by indicating that water stress was an important limitation to T. ciliata survival and growth on our study sites. The linkage between T. ciliata establishment success, early growth and soil water availability as indicated by ForWaDy, leads us to conclude that the model is a suitable stand management tool for guiding establishment of T. ciliata plantations.     

Integrating effects of leaf nitrogen,age, rank,and growth temperature into the photosynthesis-stomatal conductance model LEAFC3-N parameterised for barley (Hordeum vulgare L.)

A crucial challenge for including biophysical photosynthesis–transpiration models into complex crop growth models is to integrate the plasticity of photosynthetic processes that is related to factors like nitrogen (N) content, age, and rank of leaves, or to the adaptation of plants to growth temperature (T_g). Here we present a new version of the combined photosynthesis-stomatal conductance model LEAFC3-N [Müller, J., Wernecke, P., Diepenbrock, W., 2005. LEAFC3-N: a nitrogen sensitive extension of the CO₂ and H₂O gas exchange model LEAFC3 parameterised and tested for winter wheat (Triticum aestivum L.). Ecological Modelling 183, 183–210.] that was revised, extended and completely re-parameterised for barley (Hordeum vulgare L.) with special regard for these factors to facilitate the use of the model in ecophysiological studies and in crop modelling. The analysis is based on novel comprehensive data on photosynthetic CO₂ and light response curves measured at two oxygen concentrations and different temperatures on leaves of barley (H. vulgare L.) differing in leaf N and chlorophyll content. Plants were grown in climatic chambers or in the field at different N and T_g.We thoroughly revised the existing and introduced new nitrogen relations for key model parameters that account for a linear increase with leaf N of V_max, J_max, T_p, and R_dmax (maximum rates of carboxylation, electron transport, triose phosphate export, and mitochondrial respiration), a saturation-type increase of φ (quantum yield of electron transport), and a non-linear decrease of θ and m (curvature of the light dependence of electron transport rate, scaling factor of the stomata model). The adaptation of photosynthetic characteristics to T_g was included into the model by linear relations that were observed between T_g and the activation energy ΔH_a of the temperature response characteristics of V_max, J_max, and T_p as well as of the nitrogen dependency of these characteristics. Based on an analysis of diurnal time courses of gas exchange rates it was found necessary including not only the relation between leaf water potential (Ψ) and stomatal conductance as used originally in LEAFC3, but additional effects on V_max and J_max. With the above-listed extensions, the model was capable to reproduce the observed plasticity and the recorded diurnal time courses of gas exchange rates fairly well. Thus, we conclude that the new model version can be used under a broad range of conditions, both for ecophysiological studies and as a submodel of crop growth models. The results presented here for barley will facilitate adapting photosynthesis models like LEAFC3-N to other C₃-species as well. The modelling of the effects of drought stress should be further elaborated in future based on more specific experiments.     

Development of TRIPLEX-Management model for simulating the response of forest growth to pre-commercial thinning

In order to simulate forest growth response to pre-commercial thinning (PCT), TRIPLEX1.0 - a process-based model designed to predict forest growth as well as carbon (C) and nitrogen (N) dynamics - was modified and improved to also simulate managed forest ecosystem thinning practices. A three-parameter Weibull distribution model was integrated to simulate thinning treatments within the newly developed TRIPLEX-Management model. The thinning intensity component within the model allows users to simulate thinning treatments by applying basal area, stand density and volume to quantify thinning intensity. Natural mortality decreased following thinning due to an increase in growing space for residual stems. Predicted litterfall pools also increased after thinning events took place. The TRIPLEX-Management model was tested against published observational data for Jack Pine (Pinus banksiana Lamb.) stands subjected to PCT in Northwestern Ontario, Canada. The coefficients of determination (R²) between the predicted and observed variables including stand density, mean DBH (diameter at breast height), the quadratic mean DBH, total volume and merchantable volume as well as belowground, aboveground, and total biomass ranged from 0.50 to 0.88 (n = 20, P < 0.001) with the exception of mean tree height (R² = 0.25, n = 20, P < 0.05). Overall, the Willmott index of agreement between predicted and observed variables ranged from 0.97 to 1.00. Results show that the TRIPLEX-Management model is generally capable of simulating growth response to PCT for Jack Pine stands.     

Mathematical analysis of change in forest carbon use efficiency with stand development: A case study on Abies veitchii Lindl.

Changes in carbon use efficiency (CUE), which is defined as the ratio of net primary production (NPP) to gross primary production (GPP), were analyzed for Abies veitchii Lindl. forests with respect to stand development by developing a simple mathematical model incorporating data on physiological variables and leaf mass ratio. A decrease in CUE with stand development was successfully expressed as a function of stand biomass (y) based on the following three assumptions: (1) a power-law relationship between mean respiration and mean individual tree mass, (2) a power-functional relationship between mean gross primary production and mean individual tree mass, and (3) self-thinning relationship between stand biomass and density. Based on this model, a parameter of CUE–y relationship was defined, and it was clarified that CUE decrease with stand development is caused not by the ratio of specific respiration rate to specific gross photosynthetic rate, but by leaf mass ratio. Since CUE is high in young forests, helpful information on selecting woody species when planting seedlings was provided from the viewpoints of reducing CO₂ in the atmosphere and global warming.     

Dealing with double trouble: consequences of single and double herbivory in Brassica juncea

In their natural environment, plants are often attacked simultaneously by many insect species. The specificity of induced plant responses that is reported after single herbivore attacks may be compromised under double herbivory and this may influence later arriving herbivores. The present study focuses on the dynamics of induced plant responses induced by single and double herbivory, and their effects on successive herbivores. Morphological (leaf length, area and trichome density) and chemical changes (leaf alkenyl and indole glucosinolates) in Brassica juncea were evaluated 4, 10, 14 and 20 days after damage by the specialist Plutella xylostella alone, or together with the generalist Spodoptera litura. To assess the biological effect of the plant’s responses, the preference and performance of both herbivores on previously induced plants were measured. We found that alkenyl glucosinolates were induced 20 days after damage by P. xylostella alone, whereas their levels were elevated as early as 4 days after double herbivory. Trichome density was increased in both treatments, but was higher after double herbivory. Interestingly, there was an overall decrease in indole glucosinolates and an increase in leaf size due to damage by P. xylostella, which was not observed during double damage. S. litura preferred and performed better on undamaged plants, whereas P. xylostella preferred damaged plants and performed better on plants damaged 14 and 10 days after single and double herbivory, respectively. Our results suggest that temporal studies involving single versus multiple attacker situations are necessary to comprehend the role of induced plant responses in plant–herbivore interactions.     

Modeling net primary production of a fast-growing forest using a light use efficiency model

As interest grows in the quantification of global carbon cycles, Light Use Efficiency (LUE) model predictions of the forest net primary production (NPP) are being developed at an accelerating rate. Such models can provide useful predictions at large scales, but evaluating their performance has been difficult. In this study, a remote sensing-based LUE model was established to estimate forest NPP. Using the forest inventory data (FID) from the regional forest inventory survey in China and established allometric biomass equations, we calculated the biomass, the biomass increment, and the NPP of Eucalyptus urophylla (E. urophylla) plantation plots in the forestry jurisdiction of the Leizhou Forestry Bureau, Southern China. The FID-based NPP and the NPP from LUE model predictions were then compared to each other. Results show that the NPP from model predictions at a spatial resolution of 30 m × 30 m varied from 0 to 265 gC/(m² month) and showed regional differences. In addition, the stand age had variable effects on the average individual biomass of the E. urophylla plantation plots. The average individual biomass of the young and mid-age forests increased exponentially and logarithmically with the stand age (R² = 0.9178 and R² = 0.8683), respectively. For young and mid-age E. urophylla plantation plots, the LUE model-predicted NPP was fairly consistent with the FID-based NPP, but the model predictions of the NPP were higher than the estimates from FID. Through the analysis of the causes of uncertainty and the possible reasons for the discrepancy between the model-based NPP and FID-based NPP, the FID-derived estimates provided a foundation for model evaluation.     

Individual-based modeling of flooding and salinity effects on a coastal swamp forest

Coastal swamps are among the rapidly vanishing wetland habitats in Louisiana. Increased flooding, nutrient and sediment deprivation, and salt-water intrusion have been implicated as probable causes of the decline of coastal swamps. We developed a two-species individual-based forest succession model to compare the growth and composition of a cypress-tupelo swamp under various combinations of flooding intensity and salinity levels, using historical time-series of stage and salinity data as inputs. Our model simulates forest succession over 500 years by representing the growth, mortality, and reproduction of individual Taxodium distichum (baldcypress) and Nyssa aquatica (water tupelo) trees in a 1-km² spatial grid of 10 m × 10 m cells that vary in water levels and salinity through differences in elevation. We independently adjusted the elevations of each cell to obtain different grid-wide mean elevations and standard deviations of elevation; this affected the temporal and spatial pattern of flooding. We calibrated the model by adjusting selected parameters until averaged basal area, stem density and wood production rates under two different mean elevations (partially versus highly flooded) were qualitatively similar to comparable values reported for swamps in the literature. Corroboration involved comparing model predictions to four well-monitored contrasting habitat sites within the Maurepas Basin, Louisiana, USA. Model predictions of both species combined showed the same patterns among sites as the data, but the model overestimated wood production and the dominance of T. distichum. Exploratory simulations predicted that increased flooding leads to swamps with reduced basal areas and stem densities, while increased salinity resulted in lower basal areas at low salinity concentration (∼1-3 psu) and complete tree mortality at higher salinity concentrations (∼2-6 psu). Our model can provide insight into the succession dynamics of coastal swamps and information for the effective design of restoration actions.     

Estimation of annual potential evapotranspiration at regional scale based on the effect of moisture on soil respiration

Potential evapotranspiration (PET) is an important component of water cycle. For traditional models derived from the principle of aerodynamics and the surface energy balance, its calculation always includes many parameters, such as net radiation, water vapor pressure, air temperature and wind speed. We found that it can be acquired in an easier way in specific regions. In this study, a new PET model (PETP model) derived from two empirical models of soil respiration was evaluated using the Penman-Monteith equation as a standard method. The results indicate that the PETP model estimation concur with the Penman-Monteith equation in sites where annual precipitation ranges from 717.71 mm to 1727.37 mm (R² = 0.68, p = 0.0002), but show large discrepancies in all sites (R² = 0.07, p = 0.1280). Then we applied our PETP model at the global scale to the regions with precipitation higher than 700 mm using 2.5° CMAP data to obtain the annual PET for 2006. As expected, the spatial pattern is satisfactory overall, with the highest PET values distributed in the lower latitudes or coastal regions, and with an average of 1292.60 ± 540.15 mm year⁻¹. This PETP model provides a convenient approach to estimate PET at regional scales.     

Water budget and the consequent duration of canopy carbon gain in a teak plantation in a dry tropical region: Analysis using a soil–plant–air continuum multilayer model

A soil–plant–air continuum multilayer model was used to numerically simulate canopy net assimilation (A_n), evapotranspiration (ET), and soil moisture in a deciduous teak plantation in a dry tropical climate of northern Thailand to examine the influence of soil drought on A_n. The timings of leaf flush and the end of the canopy duration period (CDP) were also investigated from the perspective of the temporal positive carbon gain. Two numerical experiments with different seasonal patterns of leaf area index (LAI) were carried out using above-canopy hydrometeorological data as input data. The first experiment involved seasonally varying LAI estimated based on time-series of radiative transmittance through the canopy, and the second experiment applied an annually constant LAI. The first simulation captured the measured seasonal changes in soil surface moisture; the simulated transpiration agreed with seasonal changes in heat pulse velocity, corresponding to the water use of individual trees, and the simulated A_n became slightly negative. However, in the second simulation, A_n became negative in the dry season because the decline in stomatal conductance due to severe soil drought limited the assimilation, and the simultaneous increase in leaf temperature increased dark respiration. Thus, these experiments revealed that the leaflessness in the dry season is reasonable for carbon gain and emphasized the unfavorable soil water status for carbon gain in the dry season. Examining the duration of positive A_n (DPA) in the second simulation showed that the start of the longest DPA (LDPA) in a year approached the timing of leaf flush in the teak plantation after the spring equinox. On the other hand, the end appeared earlier than that of all CDPs. This result is consistent with the sap flow stopping earlier than the complete leaf fall, implying that the carbon assimilation period ends before the completion of defoliation. The model sensitivity analysis in the second simulation suggests that a smaller LAI and slower maximum rate of carboxylation likely extend the LDPA because soil water from the surface to rooting depth is maintained longer at levels adequate for carbon gain by decreased canopy transpiration. The experiments also suggest that lower soil hydraulic conductivity and deeper rooting depth can postpone the end of the LDPA by increasing soil water retention and the soil water capacity, respectively.     

Influence of sediment nitrogen-availability on carbon and nitrogen dynamics in the seagrass Thalassia testudinum

Variations in tissue carbon (C), nitrogen (N) content, and non-structural carbohydrate (NSC) reserves in the turtle grass Thalassia testudinum Banks ex König were examined in relation to changes in sediment-N availability in Corpus Christi Bay (CCB) and lower Laguna Madre (LLM), Texas, USA, from May to October 1997. Under natural conditions, sediment pore-water NH⁺ ₄-concentrations were higher at CCB (100?μM) than at LLM (30?μM); this difference was reflected in a significantly higher leaf and rhizome N-content at CCB than at LLM. However, sediment NH⁺ ₄-enrichment using a commercial fertilizer resulted in significantly higher tissue N-content relative to controls at both sites. N enrichment also influenced plant carbon metabolism, as reflected by distinct increases in leaf C-content at both sites. Significant decreases in rhizome NSC-content was recorded during the first two months of the experiment, suggesting that C was reallocated from rhizomes to leaves to support stimulated leaf growth at both sites. At LLM, leaf growth-rates increased and leaf turnover-time decreased as a result of sediment NH⁺ ₄-enrichment. With respect to chlorophyll, concentrations did not change significantly at CCB, but increased steadily at LLM after the first month following fertilization. In general however, chlorophyll concentrations in control plots were significantly higher at CCB than that at LLM. These observations suggest that leaf function related to C-fixation is enhanced under higher sediment N conditions, as reflected in higher leaf growth-rates and increased blade chlorophyll-content. In contrast, under low-N conditions, below-ground tissue production is enhanced at the expense of the above-ground shoots and leaves, resulting in the high below: above-ground biomass ratios often observed in seagrass beds of oligotrophic environments.     

A process-based population dynamics model to explore target and non-target impacts of a biological control agent

The risks and benefits associated with efforts to control invasive alien species using classical biological control are being subjected to increasing scrutiny. A process-based population dynamics model was developed to explore the interactions between a folivorous biological control agent, Cleopus japonicus, and its plant host Buddleja davidii. The model revealed that climate could have a significant impact upon the interactions between B. davidii and C. japonicus. At the coolest sites, the impact of C. japonicus on B. davidii was slowed, but it was still eventually capable of controlling populations of B. davidii. At the warmer sites where both B. davidii and C. japonicus grew faster, B. davidii succumbed rapidly to weevil damage. We hypothesise that barring an encounter with a natural enemy, C. japonicus will eventually be able to provide sustained control B. davidii throughout the North Island of New Zealand. The model scenarios illustrate the potential for the C. japonicus population to attain high densities rapidly, and to defoliate patches of B. davidii, creating the potential for spill-over feeding on non-target plants. The potential magnitude of this threat will depend partly on the climate suitability for C. japonicus, the pattern by which it migrates in response to a reduction in the available leaf resource, and the suitability of non-target plants as hosts. In all migration scenarios considered, the pattern of population growth and resource consumption by C. japonicus was exponential, with a strong tendency toward complete utilisation of resource patches more quickly at the warmer compared to colder sites. In addition to providing some useful hypotheses about the effects of climate on the biological control system, and the non-target risks, it also provides some insight into the mechanisms by which climate affects the system.