Modeling seasonal course of carbon fluxes and evapotranspiration in response to low temperature and moisture in a boreal Scots pine ecosystem

Environmental conditions act above and below ground, and regulate carbon fluxes and evapotranspiration. The productivity of boreal forest ecosystems is strongly governed by low temperature and moisture conditions, but the understanding of various feedbacks between vegetation and environmental conditions is still unclear. In order to quantify the seasonal responses of vegetation to environmental factors, the seasonality of carbon and heat fluxes and the corresponding responses for temperature and moisture in air and soil were simulated by merging a process-based model (CoupModel) with detailed measurements representing various components of a forest ecosystem in Hyytiälä, southern Finland. The uncertainties in parameters, model assumptions, and measurements were identified by generalized likelihood uncertainty estimation (GLUE). Seasonal and diurnal courses of sensible and latent heat fluxes and net ecosystem exchange (NEE) of CO₂ were successfully simulated for two contrasting years. Moreover, systematic increases in efficiency of photosynthesis, water uptake, and decomposition occurred from spring to summer, demonstrating the strong coupling between processes. Evapotranspiration and NEE flux both showed a strong response to soil temperature conditions via different direct and indirect ecosystem mechanisms. The rate of photosynthesis was strongly correlated with the corresponding water uptake response and the light use efficiency. With the present data and model assumptions, it was not possible to precisely distinguish the various regulating ecosystem mechanisms. Our approach proved robust for modeling the seasonal course of carbon fluxes and evapotranspiration by combining different independent measurements. It will be highly interesting to continue using long-term series data and to make additional tests of optional stomatal conductance models in order to improve our understanding of the boreal forest ecosystem in response to climate variability and environmental conditions.     

PIXGRO: A model for simulating the ecosystem CO2 exchange and growth of spring barley

A model, PIXGRO, developed by coupling a canopy flux sub-model (PROXEL_NEE; PROcess-based piXEL Net Ecosystem CO₂ Exchange) to a vegetation structure submodel (CGRO), for simulating both net ecosystem CO₂ exchange (NEE) and growth of spring barley is described. PIXGRO is an extension of the stand-level CO₂ and H₂O-flux model PROXEL_NEE, that simulates the NEE on a process basis, but goes further to include the dry matter production, partitioning, and crop development for spring barley. Dry matter partitioned to the leaf was converted to leaf area index (LAI) using relationships for the specific leaf area (SLA). The canopy flux component, PROXEL_NEE was calibrated using information from the literature on C3 plants and was tested using CO₂ flux data from an eddy-covariance (EC) method in Finland with long-term observations. The growth component (CGRO) was calibrated using data from the literature on spring barley as well as data from the Finland site. It was then validated against field data from two sites in Germany and partly via the use of MODIS remotely sensed LAI from the Finland site.Both the diurnal and the seasonal patterns of gross CO₂ uptake were very well simulated (R² = 0.92). A slight seasonal bias may be attributed to leaf ageing. Crop growth was also well simulated; simulated dry matter agreed with field observed data from Germany (R² = 0.90). For LAI, the agreement between the simulated and observed was good (R² = 0.80), giving an indication that functions describing the conversion of fixed CO₂ to dry matter and the subsequent partitioning leaf dry matter and LAI simulation were robust and provided reliable estimates.The MODIS LAI at a resolution of 1000 m agreed poorly (R² = 0.45) with the PIXGRO simulated LAI and the observed LAI at the Finland site in 2001. We attributed this to the coarse resolution of the image and/or the small size of the barley field (about 17 ha or 0.25 km²) at the Finland site. By deriving a regression relation between the observed LAI and NDVI from a higher resolution MODIS (500 m resolution), the MODIS-recalculated LAI agreed better with the PIXGRO-simulated LAI (R² = 0.86).PIXGRO provides a prototype model bridging the disciplines of plant physiology, crop modeling and remote sensing, for use in a spatial context in evaluating carbon balances and plant growth at stand level, landscape, regional, and with some care, continental scales. Since almost 50% of the European land surface is covered by crops, such a model is needed for the dynamic estimation of LAI and NEE of croplands.     

Water and carbon fluxes above European coniferous forests modelled with artificial neural networks

Artificial neural networks are used to select a minimal set of input variables to model water vapour and carbon exchange of coniferous forest ecosystems, independently of tree species and without detailed physiological information. Neural networks are used because of their power to fit highly non-linear relations between input and output-variables. Radiation, temperature, vapour pressure deficit and time of the day showed to be the dynamic input variables that determine ecosystem water fluxes. The same variables, together with projected leaf area index are needed for modelling CO₂-fluxes. The results for the individual sites show that the neural networks found mean water and carbon flux responses to the driving variables valid for all sites. The sensitivity analysis of the derived neural networks shows that the LAI-effect of the CO₂-flux model is overfitted because of the low variability of LAI. However, the predictions of CO₂-fluxes of sites not included in the calibration set indicate that the LAI-response of the network is reliable and that results can be used as a first estimate of the net ecosystem carbon exchange of the forest sites. Independent predictions of forest ecosystem vapour fluxes were equally satisfying as empirical models specifically calibrated for the individual sites. The results indicate that both short term water and carbon fluxes of European coniferous forests can be modelled without using detailed physiological and site specific information.     

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.     

Controls on vegetation structure in Southwestern ponderosa pine forests, 1941 and 2004

Long-term studies can broaden our ecological understanding and are particularly important when examining contingent effects that involve changes to dominance by long-lived species. Such a change occurred during the last century in Southwestern (USA) ponderosa pine (Pinus ponderosa) forests. We used five livestock grazing exclosures established in 1912 to quantify vegetation structure in 1941 and 2004. Our objectives were to (1) assess the effects of historical livestock grazing on overstory structure and age distribution, (2) assess the effects of recent livestock grazing and overstory on understory vegetation, and (3) quantify and explain changes in understory vegetation between 1941 and 2004. In 1941, canopy cover of tree regeneration was significantly higher inside exclosures. In 2004, total tree canopy cover was twice as high, density was three times higher, trees were smaller, and total basal area was 40% higher inside exclosures. Understory species density, herbaceous plant density, and herbaceous cover were negatively correlated with overstory vegetation in both years. Most understory variables did not differ between grazing treatments in 1941 but were lower inside exclosures in 2004. Differences between grazing treatments disappeared once overstory effects were accounted for, indicating that they were due to the differential overstory response to historical livestock grazing practices. Between 1941 and 2004, species density declined by 34%, herbaceous plant density by 37%, shrub cover by 69%, total herbaceous cover by 59%, graminoid cover by 39%, and forb cover by 82%. However, these variables did not differ between grazing treatments or years once overstory effects were accounted for, indicating that the declines were driven by the increased dominance of the overstory during this period. Our results demonstrate that historical livestock grazing practices are an aspect of land-use history that can affect ecosystem development. Grazing history must be considered when extrapolating results from one site to another. In addition, the understory vegetation was more strongly controlled by the ponderosa pine overstory than by recent livestock grazing or by temporal dynamics, indicating that overstory effects must be accounted for when examining understory responses in this ecosystem.     

Identification of important factors for water vapor flux and CO2 exchange in a cropland

Water vapor flux and carbon dioxide (CO₂) exchange in croplands are crucial to water and carbon cycle research as well as to global warming evaluation. In this study, a standard three-layer feed-forward back propagation neural network technique associated with the Bayesian technique of automatic relevance determination (ARD) was employed to investigate water vapor and CO₂ exchange between the canopy of summer maize and atmosphere in responses to variations of environmental and physiological factors. These factors, namely the photosynthetically active radiation (PAR), air temperature (T), vapor pressure deficient (VPD), leaf-area index (LAI), soil water content in root zone (W), and friction velocity (U*), were used as inputs in neural network analysis. Results showed that PAR, VPD, T and LAI were the primary factors regulating both water vapor and CO₂ fluxes with VPD and W more critical to water vapor flux and PAR and T more crucial to CO₂ exchange. Furthermore, two time variables “day of the year (DOY)” and “time of the day (TOD)” could also improve the simulation results of neural network analysis. The important factors identified by the neural network technique used in this study were in the order of PAR > T > VPD > LAI > U* > TOD for water vapor flux and in the order of VPD > W > LAI > T > PAR > DOY for CO₂ exchange. This study suggests that neural network technique associated with ARD could be a useful tool for identifying important factors regulating water vapor and CO₂ fluxes in terrestrial ecosystem.     

Effects of boreal forest management practices on the climate impact of CO2 emissions from bioenergy

In Life Cycle Assessment (LCA), carbon dioxide (CO₂) emissions from biomass combustion are traditionally assumed climate neutral if the bioenergy system is CO₂ flux neutral, i.e. the quantity of CO₂ released approximately equals the amount of CO₂ sequestered in biomass. This convention is a plausible assumption for fast growing biomass species, but is inappropriate for slower growing biomass, like forests. In this case, the climate impact from biomass combustion can be potentially underestimated if CO₂ emissions are ignored, or overestimated, if biogenic CO₂ is considered equal to anthropogenic CO₂. The estimation of the effective climate impact should take into account how the CO₂ fluxes are distributed over time: the emission of CO₂ from bioenergy approximately occurs at a single point in time, while the absorption by the new trees is spread over several decades. Our research target is to include this dynamic time dimension in unit-based impact analysis, using a boreal forest stand as case study. The boreal forest growth is modelled with an appropriate function, and is investigated under different forestry regimes (affecting the growth rate and the year of harvest). Specific atmospheric decay functions for biomass-derived CO₂ are then elaborated for selected combinations of forest management options. The contribution to global warming is finally quantified using the GWP_bio index as climate metric. Results estimates the effects of these practices on the characterization factor used for the global warming potential of CO₂ from bioenergy, and point out the key role played by the selected time horizon.     

Leaf-area index and light attenuation in rapidly expanding shrub thickets

There is increasing interest in the changes in ecosystem services that accompany the conversion of grasslands to shrub-dominated communities. Shrub structure and associated effects on the light environment may be especially important in affecting productivity and diversity. Leaf-area index (LAI) and understory light levels of Morella cerifera shrub thickets were assessed on Hog Island, Virginia, USA, at four sites along a soil chronosequence. LAI was estimated from annual leaf litter, with allometric models relating stem diameter to leaf area, with a portable integrating radiometer (LI-COR LAI-2000), and from photosynthetically active radiation (PAR) using the Beer-Lambert law. For the two youngest thickets, LAI estimates from leaf litter (approximately 10.0) approached levels often associated with tropical rain forest. Allometric models estimated LAI values at 9.8 and 12.5 for the same thickets. High LAI in thickets also results in high light attenuation. Light levels within thickets were as low as 0.7% of above-canopy PAR in the youngest thicket. These data suggest that M. cerifera shrub thickets have a very high potential for annual net primary production. Furthermore, extreme modification of the light environment, coupled with heavy shrub litter fall, may exclude potential competitors during thicket establishment and rapidly alter community structure and ecosystem function.     

More asymmetric tree competition brings about more evapotranspiration and less runoff from the forest ecosystems: A simulation study

The present paper reports how stand size-structure dynamics due to competition between different-sized trees affect long-term forested water balance in Japanese cool-temperate planted stands (evergreen coniferous Cryptomeria japonica and deciduous coniferous Larix kaempferi stands) using a fully coupled multi-layered meteorological surface physics—terrestrial ecosystems model. The simulation captured the well-known annual variation in leaf area index (LAI) accurately with stand age in monocultured and even-aged stands; the occurrence of maximum LAI during the early growth stage and then a gradual decline followed by a steady state after the maximum LAI. The simulations also detected a high dependency of annual evapotranspiration (AETr) on LAI with stand age that is well known by prior observational researches. In the C. japonica (shade-tolerant late-successional species) stand, the relationship between annual net primary productivity of an individual tree (NPP_ind) and individual tree mass (w) changed from linear to a convex curve during self-thinning, indicating that the degree of asymmetric tree competition intensified with forest stand development. The higher degree of competitive asymmetry characterized by the convex-shaped NPPind-w relationship produced greater size inequality, i.e., the formation of trees stratified by height. Under such conditions, AETr and annual transpiration (ATr) were mainly regulated by larger trees. On the other hand, the NPPind-w relationships in the L. kaempferi (shade-intolerant early-successional species) stand were linear throughout the simulated period, indicating the lower degree of competitive asymmetry. Under such conditions, the growth of intermediate-sized trees was enhanced and these trees became a dominant source of AETr (and also ATr) during self-thinning. Furthermore, the sensitivity analysis of the effects of ecophysiological parameters such as foliage profile (i.e., vertical distribution of leaf area density) of an individual tree (distribution pattern is described by the parameter η), the maximum carboxylation velocity (V_cmax0) and biomass allocation pattern of individual plant growth (μ₁) on AETr, ATr and annual runoff (AR_off) showed that the temporal trends of AETr, ATr, AR_off and NPPind-w relationships were completely the same as those in the control simulations. However, the NPPind-w relationship during self-thinning indicated higher degrees of competitive asymmetry when η or V_cmax0 were greater than those in the control simulation and generated greater AETr and ATr and thus smaller AR_off. We found that more asymmetric tree competition brings about greater size inequality between different-sized trees and thus more evapotranspiration and less runoff in a forest stand. Overall, our simulation approach revealed that not only LAI dynamics but also plant competition, and thus size-structure dynamics, in a forest ecosystem are essential to long-term future projections of forested water balance.     

Long-term aboveground and belowground consequences of red wood ant exclusion in boreal forest

Despite their ubiquity, the role of ants in driving ecosystem processes both aboveground and belowground has been seldom explored, except within the nest. During 1995 we established 16 ant exclusion plots of approximately 1.1 x 1.1 m, together with paired control plots, in the understory layer of a boreal forest ecosystem in northern Sweden that supports high densities of the mound-forming ant Formica aquilonia, a red wood ant species of the Formica rufa group. Aboveground and belowground measurements were then made on destructively sampled subplots in 2001 and 2008, i.e., 6 and 13 years after set-up. While ant exclusion had no effect on total understory plant biomass, it did greatly increase the relative contribution of herbaceous species, most likely through preventing ants from removing their seeds. This in turn led to higher quality resources entering the belowground subsystem, which in turn stimulated soil microbial biomass and activity and the rates of loss of mass and carbon (C) and nitrogen (N) from litter in litterbags placed in the plots. This was accompanied by losses of approximately 15% of N and C stored in the humus on a per area basis. Ant exclusion also had some effects on foliar stable isotope ratios for both C and N, most probably as a consequence of greater soil fertility. Further, exclusion of ants had multitrophic effects on a microbe-nematode soil food web with three consumer trophic levels and after six years promoted the bacterial-based relative to the fungal-based energy channel in this food web. Our results point to a major role of red wood ants in determining forest floor vegetation and thereby exerting wide-ranging effects on belowground properties and processes. Given that the boreal forest occupies 11% of the Earth's terrestrial surface and stores more C than any other forest biome, our results suggest that this role of ants could potentially be of widespread significance for biogeochemical nutrient cycling, soil nutrient capital, and sequestration of belowground carbon.     

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.     

A simplified approach to implement forest eco-hydrological properties in regional hydrological modelling

The topic of this paper is a simplified model for simulating the hydrological properties of forest stands based on a robust computation of the temporal LAI (leaf area index) dynamics. The approach allows the simulation of all hydrologically relevant processes. It includes interception of precipitation and transpiration of forest stands with and without groundwater in the rooting zone. The model also considers phenology, mortality and simple management practice. It was implemented as a module in the eco-hydrological model SWIM (Soil and Water Integrated Model). The approach was tested on Scots pine (Pinus sylvestris) and common oak (Quercus robur and Q. petraea).The results demonstrate a good simulation of annual biomass increase and LAI and satisfactory simulation of litter production (annual mean value). A comparison of the date of May sprout for Scots pine and leaf unfolding for Oak (1980–1990) with observed data of the DWD (German Weather Service) shows a good reproduction of the temporal dynamic. The daily simulation of transpiration shows an excellent correlation of r = 0.81 for the year 1998 but only r = 0.65 for 1999. The interception losses were also simulated and compared with weekly observed data showing satisfactory results in the vegetation periods and annual sums, but worse agreement in autumn and spring time. A regional assessment study was done in the federal state of Brandenburg (Germany) to test the applicability and multi-criteria evaluation capabilities of the approach on the landscape and catchments scale using forest data, daily river discharge and regional water balance.     

川西高山森林生态系统林下生物量及其随林窗的变化特征

作为森林生态系统的重要组成部分,林下植被及其残体的分布受到林冠层的影响,但迄今有关林窗对林下植被和残体生物量的影响尚无研究报道.于2013 年8 月2 日至20 日,以海拔3 600 m 的川西岷江冷杉原始林林下植被为研究对象,根据区域内的坡向和林分组成等因素设置3 个100 m×100 m 的典型样地,调查其生物量及其随林窗的变化特征.在每个样地内选择3 个大林窗,在林窗、林缘和林下分别设置3 个20 m×20 m 的样方,调查粗木质残体长度或高度、大小头直径、枯立木记录胸径、腐烂等级等;在林窗、林缘和林下分别设置3 个5 m×5 m 的样方,采用“收获法”收集样方内直径在2.5-10 cm之间的细木质残体和灌木生物量;在林窗、林缘和林下分别设置3 个1 m×1 m 的样方来调查凋落物储量和草本生物量;在1m×1 m 的样方内随机选择1 个20 cm×20 cm 的小样方来调查地被植物生物量.结果表明,（1）川西高山森林生态系统总生物量为72.75 t·hm^-2,其中林下生物量为67.92 t·hm^-2,占生态系统生物量的95.17%.活体植被以灌木为主,其生物量为9.81t·hm^-2;残体部分以粗木质残体为主,其储量为53.00 t·hm^-2;（2）林窗对灌木、草本、地被植物的影响各不相同,且不同物种的灌木生物量表现出不同的分布规律;草本生物量表现出明显的“边缘效应”,在林缘显著高于林下;林窗和林缘的地被植物生物量相对较低;（3）粗木质残体储量从林下到林窗呈现减小的趋势,但总体储量仍然较大,林窗和林缘的细木质残体储量高于林下.这些结果为认识高山森林生态系统林下生物量及其格局,以及林窗在森林生态系统的重要作用提供了基础理论依据.     

Beurteilung des Ozonrisikos für die Waldregionen Bayerns am Beispiel des Jahres 2002 und des Extremtrockenjahres 2003 auf der Basis der externen Ozonexposition und der internen Ozonaufnahme

Background, aim, and scope Increasing background concentrations of ground-level tropospheric ozone and more frequent and prolonged summer drought incidences due to climate change are supposed to increase the stress on Bavarian forests. For such scenarios growth reduction and yield losses are predicted. Sustainable forest management in Bavaria aims to significantly increase the proportion of beech (Fagus sylvatica L.) because of its broad ecological amplitude. In our regional study different approaches for calculating ozone impact were used to estimate the risks for Bavarian forests in the average climatic, rather moist year 2002 and the extremely dry year 2003.Materials and methods Measurements were conducted for eleven forest ecosystem sites and two forest research sites representing typical Bavarian forest stands under different climatic conditions and situated in different altitudes. For risk assessment currently used approaches were applied either based on the calculation of the cumulative ozone exposure (external dose; MPOC maximal permitted ozone concentration; critical level AOT40_phen? accumulated ozone exposure over a threshold of 40 nl [O₃] l^–1, for the effective phenolgy of beech) or based on the calculation of the phytomedically relevant ozone flux into the stomata (internal dose, critical level AF_st>1,6, accumulated stomatal flux above a flux threshold of 1.6 nmol O₃?m^–2 PLA; PLA = projected leaf area). For calculations continuously recorded ozone concentrations and meteorological and phenological data from nearby rural open field background measuring stations from the national air pollution control and from forested sites were used. Additionally ozone induced leaf symptoms were assessed.Results The exposure-based indices AOT40_phen and MPOC as well as the flux-based index AF_st>1.6suggest that Bavarian forests are at risk from O₃ during a rather moist average year concerning climate conditions (2002) as well as in an extreme dry year (2003). Thus, growth reductions of 5?% are predicted when thresholds are exceeded. Threshold exceedance occurred in both years at all plots, mostly already at the beginning of the growing season and often even many times over. Ozone induced leaf symptoms could be detected only on a few plots in a very slight occurrence.Discussion The results for the applied critical level indices differed depending on climatic conditions during the growing seasons: Regarding exposure-based indices, the highest degree of threshold exceedance occurred in the dry year of 2003 at all plots; the flux-based approach indicated the highest stomatal ozone uptake and thus an increased risk at moist sites or during humid years, whereas the risk was decreasing at dry sites with prolonged water limitation. Hence, soil and accordingly plant water availability was the decisive factor for the flux-modelled internal ozone uptake via stomata. Drought and increased ozone impact can generate synergistic, but also antagonistic effects for forest trees. At water limited rather dry forest sites restricted transpiration and thus production, but concurrently lower ozone uptake and reduced risk for damage can be expected.Conclusions, recommendations, and perspectives For realistic site-specific risk assessment in forest stands the determination of the internal ozone dose via modeling flux based internal stomatal ozone uptake is more appropriate than the calculation of the external ozone dose. The predicted 5?% growth reductions are in discrepancy with the frequently observed increment increase during the last decades in forest stands. Comprehensive and significant statistical verification for ozone induced forest growth reduction as well as the systematic validation of thresholds for ozone in the field is still lacking. However, a multiplicity of different specific new and retrospective growth analysis data should allow closing the gap. Moreover, the determination of canopy transpiration with sap flow measurements is a novel approach to provide cause-effect related, site specific results for the effective internal ozone dose as well as for canopy water supply and consecutively for regional risk estimation. A further future objective is the refinement of O₃ flux modelling by further consideration of soil/water budget characteristics and the above mentioned improved estimations of crown and canopy transpiration. Further, the introduction of threshold ranges for forest trees in view of their specific regional climatic conditions and their validation in real forest stands is necessary for developing meaningful ozone risk predictions for forests.     

Effects of in situ CO2 enrichment on the structural and chemical characteristics of the seagrass Thalassia testudinum

Seagrasses commonly display carbon-limited photosynthetic rates. Thus, increases in atmospheric pCO₂, and consequentially oceanic CO_2(aq) concentrations, may prove beneficial. While addressed in mesocosms, these hypotheses have not been tested in the field with manipulative experimentation. This study examines the effects of in situ CO_2(aq) enrichment on the structural and chemical characteristics of the tropical seagrass, Thalassia testudinum. CO_2(aq) availability was manipulated for 6 months in clear, open-top chambers within a shallow seagrass meadow in the Florida Keys (USA), reproducing forecasts for the year 2100. Structural characteristics (leaf area, leaf growth, shoot mass, and shoot density) were unresponsive to CO_2(aq) enrichment. However, leaf nitrogen and phosphorus content declined on average by 11 and 21 %, respectively. Belowground, non-structural carbohydrates increased by 29 %. These results indicate that increased CO_2(aq) availability may primarily alter the chemical composition of seagrasses, influencing both the nutrient status and resilience of these systems.     

Estimating California ecosystem carbon change using process model and land cover disturbance data: 1951-2000

Land use change, natural disturbance, and climate change directly alter ecosystem productivity and carbon stock level. The estimation of ecosystem carbon dynamics depends on the quality of land cover change data and the effectiveness of the ecosystem models that represent the vegetation growth processes and disturbance effects. We used the Integrated Biosphere Simulator (IBIS) and a set of 30- to 60-m resolution fire and land cover change data to examine the carbon changes of California's forests, shrublands, and grasslands. Simulation results indicate that during 1951-2000, the net primary productivity (NPP) increased by 7%, from 72.2 to 77.1 Tg C yr⁻¹ (1 teragram = 10¹² g), mainly due to CO₂ fertilization, since the climate hardly changed during this period. Similarly, heterotrophic respiration increased by 5%, from 69.4 to 73.1 Tg C yr⁻¹, mainly due to increased forest soil carbon and temperature. Net ecosystem production (NEP) was highly variable in the 50-year period but on average equalled 3.0 Tg C yr⁻¹ (total of 149 Tg C). As with NEP, the net biome production (NBP) was also highly variable but averaged −0.55 Tg C yr⁻¹ (total of -27.3 Tg C) because NBP in the 1980s was very low (-5.34 Tg C yr⁻¹). During the study period, a total of 126 Tg carbon were removed by logging and land use change, and 50 Tg carbon were directly removed by wildland fires. For carbon pools, the estimated total living upper canopy (tree) biomass decreased from 928 to 834 Tg C, and the understory (including shrub and grass) biomass increased from 59 to 63 Tg C. Soil carbon and dead biomass carbon increased from 1136 to 1197 Tg C.Our analyses suggest that both natural and human processes have significant influence on the carbon change in California. During 1951-2000, climate interannual variability was the key driving force for the large interannual changes of ecosystem carbon source and sink at the state level, while logging and fire were the dominant driving forces for carbon balances in several specific ecoregions. From a long-term perspective, CO₂ fertilization plays a key role in maintaining higher NPP. However, our study shows that the increase in C sequestration by CO₂ fertilization is largely offset by logging/land use change and wildland fires.     

Application of a three-dimensional model for assessing effects of small clear-cuttings on radiation and soil temperature

A three-dimensional model Mixfor-3D of soil–vegetation–atmosphere transfer (SVAT) was developed and applied to estimate possible effects of tree clear-cutting on radiation and soil temperature regimes of a forest ecosystem. The Mixfor-3D model consists of several closely coupled 3D sub-models describing: forest stand structure; radiative transfer in a forest canopy; turbulent transfer of sensible heat, H₂O and CO₂ between ground surface and the atmospheric surface layer; evapotranspiration of ground surface vegetation and soil; heat and moisture transfer in soil. The model operates with the horizontal grid resolution, 2 m × 2 m; vertical resolution, 1 m and primary time step, 1 h.     

Evaluating weather effects on interannual variation in net ecosystem productivity of a coastal temperate forest landscape: A model intercomparison

Forest productivity is strongly affected by seasonal weather patterns and by natural or anthropogenic disturbances. However weather effects on forest productivity are not currently represented in inventory-based models such as CBM-CFS3 used in national forest C accounting programs. To evaluate different approaches to modelling these effects, a model intercomparison was conducted among CBM-CFS3 and four process models (ecosys, CN-CLASS, Can-IBIS and 3PG) over a 2500 ha landscape in the Oyster River (OR) area of British Columbia, Canada. The process models used local weather data to simulate net primary productivity (NPP), net ecosystem productivity (NEP) and net biome productivity (NBP) from 1920 to 2005. Other inputs used by the process and inventory models were generated from soil, land cover and disturbance records. During a period of intense disturbance from 1928 to 1943, simulated NBP diverged considerably among the models. This divergence was attributed to differences among models in the sizes of detrital and humus C stocks in different soil layers to which a uniform set of soil C transformation coefficients was applied during disturbances. After the disturbance period, divergence in modelled NBP among models was much smaller, and attributed mainly to differences in simulated NPP caused by different approaches to modelling weather effects on productivity. In spite of these differences, age-detrended variation in annual NPP and NEP of closed canopy forest stands was negatively correlated with mean daily maximum air temperature during July-September (T_amax) in all process models (R² = 0.4-0.6), indicating that these correlations were robust. The negative correlation between T_amax and NEP was attributed to different processes in different models, which were tested by comparing CO₂ fluxes from these models with those measured by eddy covariance (EC) under contrasting air temperatures (T_a). The general agreement in sensitivity of annual NPP to T_amax among the process models led to the development of a generalized algorithm for weather effects on NPP of coastal temperate coniferous forests for use in inventory-based models such as CBM-CFS3: NPP′ = NPP − 57.1 (T_amax − 18.6), where NPP and NPP′ are the current and temperature-adjusted annual NPP estimates from the inventory-based model, 18.6 is the long-term mean daily maximum air temperature during July-September, and T_amax is the mean value for the current year. Our analysis indicated that the sensitivity of NPP to T_amax was nonlinear, so that this algorithm should not be extrapolated beyond the conditions of this study. However the process-based methodology to estimate weather effects on NPP and NEP developed in this study is widely applicable to other forest types and may be adopted for other inventory based forest carbon cycle models.     

Quality control and flux gap filling strategy for Bowen ratio method: revisiting the Priestley–Taylor evaporation model

Micrometeorological measurements were made over an irrigated rice paddy during the Huaihe River Basin Experiment (HUBEX) in 1999. This study addresses the quality control and gap filling strategy for the heat fluxes with the Bowen ratio energy balance (BREB) method. It also endeavors to benefit future studies through comparing five methodologies to estimate the net water exchange. First, a three-step quality control strategy is constructed. Its first two steps guarantee the correct flux directions and reject suspicious data, respectively. The third step forbids supersaturation by considering the Bowen ratio ranges for different flux combinations (termed as “scenarios”). The quality-controlled latent heat (LE) and sensible heat (H) fluxes fall in three scenarios, namely I (LE > 0, H > 0), II (LE > 0, H < 0), and IV (LE < 0, H < 0). Second, the Priestley–Taylor evaporation model (PTEM) is applied to fill the gaps of LE, while the energy balance relationship is used to gap-fill H, namely H = R _N − G − LE (R _N : net radiation; G: soil heat flux). Central to the success of this strategy is the idea to calibrate the Priestley–Taylor parameter (α) in a scenario-specific manner. On average, α values are calibrated as 1.20 and 1.35 for scenario I and II, respectively. For scenario IV, most α values lie in a narrow range, namely from 0.9 to 1.0. Then, α is calibrated as 0.97 to extend the applicability of the PTEM to condensation (negative LE). The scenario-specific treatment explicitly explains the diurnal variation of α derived without distinguishing the scenarios. Third, five methodologies are compared in the calculation of net water exchange, including PTEM-gap-filled BREB method (M1), energy-balance-based eddy covariance method (M2), and three simplified models to estimate LE, respectively as R _N − H, R _N − G, and R _N . A major finding is that G and H have a similar effect of about 3% in the net water exchange. Thus, if either is neglected, the net evaporation should be slightly lowered to avoid overestimate, by 3% as a rule of thumb from this work.