基于模型数据融合的长白山阔叶红松林碳循环模拟

 充分、有效地利用各种陆地生态系统碳观测数据改善陆地生态系统模型, 是当前我国陆地生态系统碳循环研究领域亟待解决的重要问题之一。该研究以2003~2005年长白山阔叶红松林的6组生物计量观测数据和涡度相关技术测定的碳通量数据为基础, 利用马尔可夫链-蒙特卡罗方法对陆地生态系统模型的关键参数(即碳滞留时间)进行了反演, 进而预测了长白山阔叶红松林生态系统碳库、碳通量及其不确定性。反演结果表明, 长白山阔叶红松林叶凋落物和微生物碳的平均滞留时间最短, 为2~6个月; 其次是叶和细根生物量碳, 二者的平均滞留时间为1~2 a; 慢性土壤有机碳的平均滞留时间为8~16 a; 碳在木质生物量和惰性土壤有机质库中的滞留时间最长, 平均滞留时间分别为77~109 a和409~1 879 a。模拟结果显示, 碳库和累积碳通量模拟值的不确定性将随着模拟时间的延长而增大。当气温升高10%和20%时, 长白山阔叶红松林总初级生产力年总量将分别增加6.5%和9.9%, 净生态系统生产力(NEP)年总量的变化取决于土壤温度的变化。若土壤温度保持不变, NEP年总量将分别增加11.4%~21.9%和17.6%~33.1%; 若土壤温度也相应升高10%和20%, NEP年总量的增幅反而下降甚至低于原来的水平。假设气候和植被保持在2003~2005年的状态, 2020年长白山阔叶红松林NEP年总量为(163±12) g C·m^–2·a^–1, 土壤呼吸年总量为(721±14) g C·m^–2·a^–1。马尔可夫链-蒙特卡罗方法是反演模型参数、优化模拟结果和评估模拟结果不确定性的有效方法, 但今后仍需在惰性土壤碳滞留时间的估计、驱动数据和模型结构的不确定性分析、模型数据融合方法方面进行深入研究, 以进一步提高碳循环模拟的准确性。     

Simulated long‐term effects of harvest and biomass residue removal on soil carbon and nitrogen content and productivity for two Upper Great Lakes forest ecosystems

We used the ecosystem process model Biome‐BGC to simulate the effects of harvest and residue removal management scenarios on soil carbon (C), available soil nitrogen (N), net primary production (NPP), and net ecosystem production (NEP) in jack pine (Pinus banksiana Lamb.) and sugar maple (Acer saccharum Marsh) ecosystems in northern Wisconsin, USA. To assess harvest effects, we simulated short (50‐year) and long (100‐year) harvest intervals, high (clear‐cut) and low (selective) harvest intensities, and three levels of residue retention (15%, 25%, and 35%) over a 500‐year period. The model simulation of NPP, soil C accumulation, and NEP agreed reasonably well with biometric and eddy‐covariance measurements of these two ecosystems. The more intensive (50‐year rotation clear‐cuts with low residue retention) harvest scenarios tended to have the greatest NEP (420 and 678 t C ha^?1 for the 500‐year interval for jack pine and sugar maple, respectively). All the harvest scenarios decreased mineral soil C and available mineral soil N content relative to the no‐harvest scenario for jack pine and sugar maple. The rate of change in mineral soil C decreased the greatest in the most intensive biomass removal scenarios (?0.012 and ?0.072 t C ha^?1 yr^?1 relative to no‐harvest for jack pine and sugar maple, respectively) and the smallest decrease was observed in the least intensive biomass removal scenarios (?0.002 and ?0.009 t C ha^?1 yr^?1 relative to no‐harvest for jack pine and sugar maple, respectively). The more intensive biomass removal harvest scenarios in sugar maple significantly decreased peak productivity (NPP) in the simulation period.     

Carbon storage and fluxes in ponderosa pine forests at different developmental stages

We compared carbon storage and fluxes in young and old ponderosa pine stands in Oregon, including plant and soil storage, net primary productivity, respiration fluxes, eddy flux estimates of net ecosystem exchange (NEE), and Biome‐BGC simulations of fluxes. The young forest (Y site) was previously an old‐growth ponderosa pine forest that had been clearcut in 1978, and the old forest (O site), which has never been logged, consists of two primary age classes (50 and 250 years old). Total ecosystem carbon content (vegetation, detritus and soil) of the O forest was about twice that of the Y site (21 vs. 10 kg C m^?2 ground), and significantly more of the total is stored in living vegetation at the O site (61% vs. 15%). Ecosystem respiration (R_e) was higher at the O site (1014 vs. 835 g C m^?2 year^?1), and it was largely from soils at both sites (77% of R_e). The biological data show that above‐ground net primary productivity (ANPP), NPP and net ecosystem production (NEP) were greater at the O site than the Y site. Monte Carlo estimates of NEP show that the young site is a source of CO₂ to the atmosphere, and is significantly lower than NEP(O) by c. 100 g C m^?2 year^?1. Eddy covariance measurements also show that the O site was a stronger sink for CO₂ than the Y site. Across a 15‐km swath in the region, ANPP ranged from 76 g C m^?2 year^?1 at the Y site to 236 g C m^?2 year^?1 (overall mean 158 ± 14 g C m^?2 year^?1). The lowest ANPP values were for the youngest and oldest stands, but there was a large range of ANPP for mature stands. Carbon, water and nitrogen cycle simulations with the Biome‐BGC model suggest that disturbance type and frequency, time since disturbance, age‐dependent changes in below‐ground allocation, and increasing atmospheric concentration of CO₂ all exert significant control on the net ecosystem exchange of carbon at the two sites. Model estimates of major carbon flux components agree with budget‐based observations to within ± 20%, with larger differences for NEP and for several storage terms. Simulations showed the period of regrowth required to replace carbon lost during and after a stand‐replacing fire (O) or a clearcut (Y) to be between 50 and 100 years. In both cases, simulations showed a shift from net carbon source to net sink (on an annual basis) 10–20 years after disturbance. These results suggest that the net ecosystem production of young stands may be low because heterotrophic respiration, particularly from soils, is higher than the NPP of the regrowth. The amount of carbon stored in long‐term pools (biomass and soils) in addition to short‐term fluxes has important implications for management of forests in the Pacific North‐west for carbon sequestration.     

Evaluating the impacts of climate and elevated carbon dioxide on tropical rainforests of the western Amazon basin using ecosystem models and satellite data

Forest inventories from the intact rainforests of the Amazon indicate increasing rates of carbon gain over the past three decades. However, such estimates have been questioned because of the poor spatial representation of the sampling plots and the incomplete understanding of purported mechanisms behind the increases in biomass. Ecosystem models, when used in conjunction with satellite data, are useful in examining the carbon budgets in regions where the observations of carbon flows are sparse. The purpose of this study is to explain observed trends in normalized difference vegetation index (NDVI) using climate observations and ecosystem models of varying complexity in the western Amazon basin for the period of 1984–2002. We first investigated trends in NDVI and found a positive trend during the study period, but the positive trend in NDVI was observed only in the months from August to December. Then, trends in various climate parameters were calculated, and of the climate variables considered, only shortwave radiation was found to have a corresponding significant positive trend. To compare the impact of each climate component, as well as increasing carbon dioxide (CO₂) concentrations, on evergreen forests in the Amazon, we ran three ecosystem models (CASA, Biome‐BGC, and LPJ), and calculated monthly net primary production by changing a climate component selected from the available climate datasets. As expected, CO₂ fertilization effects showed positive trends throughout the year and cannot explain the positive trend in NDVI, which was observed only for the months of August to December. Through these simulations, we demonstrated that the positive trend in shortwave radiation can explain the positive trend in NDVI observed for the period from August to December. We conclude that the positive trend in shortwave radiation is the most likely driver of the increasing trend in NDVI and the corresponding observed increases in forest biomass.     

Disturbance and climate effects on carbon stocks and fluxes across Western Oregon USA

We used a spatially nested hierarchy of field and remote‐sensing observations and a process model, Biome‐BGC, to produce a carbon budget for the forested region of Oregon, and to determine the relative influence of differences in climate and disturbance among the ecoregions on carbon stocks and fluxes. The simulations suggest that annual net uptake (net ecosystem production (NEP)) for the whole forested region (8.2 million hectares) was 13.8 Tg C (168 g C m^?2 yr^?1), with the highest mean uptake in the Coast Range ecoregion (226 g C m^?2 yr^?1), and the lowest mean NEP in the East Cascades (EC) ecoregion (88 g C m^?2 yr^?1). Carbon stocks totaled 2765 Tg C (33 700 g C m^?2), with wide variability among ecoregions in the mean stock and in the partitioning above‐ and belowground. The flux of carbon from the land to the atmosphere that is driven by wildfire was relatively low during the late 1990s (～0.1 Tg C yr^?1), however, wildfires in 2002 generated a much larger C source (～4.1 Tg C). Annual harvest removals from the study area over the period 1995–2000 were ～5.5 Tg C yr^?1. The removals were disproportionately from the Coast Range, which is heavily managed for timber production (approximately 50% of all of Oregon's forest land has been managed for timber in the past 5 years). The estimate for the annual increase in C stored in long‐lived forest products and land fills was 1.4 Tg C yr^?1. Net biome production (NBP) on the land, the net effect of NEP, harvest removals, and wildfire emissions indicates that the study area was a sink (8.2 Tg C yr^?1). NBP of the study area, which is the more heavily forested half of the state, compensated for ～52% of Oregon's fossil carbon dioxide emissions of 15.6 Tg C yr^?1 in 2000. The Biscuit Fire in 2002 reduced NBP dramatically, exacerbating net emissions that year. The regional total reflects the strong east–west gradient in potential productivity associated with the climatic gradient, and a disturbance regime that has been dominated in recent decades by commercial forestry.     

Across-model spread and shrinking in predicting peatland carbon dynamics under global change

Large across-model spread in simulating land carbon (C) dynamics has been ubiquitously demonstrated in model intercomparison projects (MIPs), and became a major impediment in advancing climate change prediction. Thus, it is imperative to identify underlying sources of the spread. Here, we used a novel matrix approach to analytically pin down the sources of across-model spread in transient peatland C dynamics in response to a factorial combination of two atmospheric CO₂ levels and five temperature levels. We developed a matrix-based MIP by converting the C cycle module of eight land models (i.e., TEM, CENTURY4, DALEC2, TECO, FBDC, CASA, CLM4.5 and ORCHIDEE) into eight matrix models. While the model average of ecosystem C storage was comparable to the measurement, the simulation differed largely among models, mainly due to inter-model difference in baseline C residence time. Models generally overestimated net ecosystem production (NEP), with a large spread that was mainly attributed to inter-model difference in environmental scalar. Based on the sources of spreads identified, we sequentially standardized model parameters to shrink simulated ecosystem C storage and NEP to almost none. Models generally captured the observed negative response of NEP to warming, but differed largely in the magnitude of response, due to differences in baseline C residence time and temperature sensitivity of decomposition. While there was a lack of response of NEP to elevated CO₂ (eCO₂) concentrations in the measurements, simulated NEP responded positively to eCO₂ concentrations in most models, due to the positive responses of simulated net primary production. Our study used one case study in Minnesota peatland to demonstrate that the sources of across-model spreads in simulating transient C dynamics can be precisely traced to model structures and parameters, regardless of their complexity, given the protocol that all the matrix models were driven by the same gross primary production and environmental variables.     

Ecosystem carbon pools,fluxes, and balances within mature tallgrass prairie restorations

Tallgrass prairie restorations can quickly accrue organic C in soil and biomass, but the rate of C accumulation diminishes through time and is highly variable among more mature prairies. Long‐term soil organic carbon (SOC) accumulation in prairies has been linked to edaphic factors such as soil texture, soil moisture, and SOC content, but it is unclear how these factors affect the ecosystem processes that are responsible for observed differences in C accumulation rates in older prairies. We measured belowground plant and SOC pools and fluxes within 27–36‐year‐old restored tallgrass prairies in order to quantify total C storage, determine the net ecosystem production of C (NEP‐C), and explore which edaphic factors influence the ecosystem processes responsible for divergent NEP‐C. We found that 11% of organic C was stored in biomass, and we estimate that one‐third of post‐restoration C sequestration has occurred in biomass, thereby highlighting biomass as a large but often overlooked C pool. Belowground biomass and soil C pools were notably smaller than those reported for remnant prairie, suggesting that future belowground C accumulation could still occur. During this study, the prairies appeared to be a net source of C, although the range of NEP‐C values encompassed zero. Sand content positively affected NEP‐C via increased belowground biomass production‐C inputs, and SOC negatively affected NEP‐C due to increased soil respiration C outputs. However, soil moisture had a smaller negative effect on soil respiration, indicating that both SOC and soil moisture play important roles in determining prairie C balance.     

Identification of vegetation and soil carbon pools out of equilibrium in a process model via eddy covariance and biometric constraints

Assumptions of steady‐state conditions in biogeochemical modelling are often invoked because knowledge on the development status of the modelling domain is generally unavailable. Here, we investigate the role of vegetation pool sizes on nonequilibrium conditions through model‐data integration approaches for a set of sites using eddy covariance CO₂ flux data. The study is based on the Carnegie–Ames–Stanford Approach (CASA) model, modified (CASA_G) in order to evaluate the sensitivity of simulated net ecosystem production (NEP) fluxes to vegetation pool sizes. The experimental design is based on the inverse model optimization of different parameter vectors performed at the measurement site level. Each parameter vector prescribes different simulation dynamics that embody different model structural assumptions concerning (non)steady‐state conditions in vegetation and soil carbon pools. We further explore the potential of assimilating biometric constraints through the cost function for sites where in situ information on aboveground biomass or wood pools is available. The integration of biometric data yields marked improvements in the simulation of vegetation C pools compared to single constraints with eddy flux data. Overall, it is necessary to relax both vegetation and soil carbon pools for consistency with the observed data streams. Multiple constraints approaches also leads to variable model performance among the different experimental setups and model structures. We identify and assess the limitations of various model structures and the role of multiple constraints approaches for tackling issues of equifinality. These studies emphasize the need for establishing consistent data sets of fluxes and biometric data for successful model‐data fusion.     

Partitioning direct and indirect human-induced effects on carbon sequestration of managed coniferous forests using model simulations and forest inventories

Temperate forest ecosystems have recently been identified as an important net sink in the global carbon budget. The factors responsible for the strength of the sinks and their permanence, however, are less evident. In this paper, we quantify the present carbon sequestration in Thuringian managed coniferous forests. We quantify the effects of indirect human‐induced environmental changes (increasing temperature, increasing atmospheric CO₂ concentration and nitrogen fertilization), during the last century using BIOME‐BGC, as well as the legacy effect of the current age‐class distribution (forest inventories and BIOME‐BGC). We focused on coniferous forests because these forests represent a large area of central European forests and detailed forest inventories were available. The model indicates that environmental changes induced an increase in biomass C accumulation for all age classes during the last 20 years (1982–2001). Young and old stands had the highest changes in the biomass C accumulation during this period. During the last century mature stands (older than 80 years) turned from being almost carbon neutral to carbon sinks. In high elevations nitrogen deposition explained most of the increase of net ecosystem production (NEP) of forests. CO₂ fertilization was the main factor increasing NEP of forests in the middle and low elevations. According to the model, at present, total biomass C accumulation in coniferous forests of Thuringia was estimated at 1.51 t C ha^?1 yr^?1 with an averaged annual NEP of 1.42 t C ha^?1 yr^?1 and total net biome production of 1.03 t C ha^?1 yr^?1 (accounting for harvest). The annual averaged biomass carbon balance (BCB: biomass accumulation rate‐harvest) was 1.12 t C ha^?1 yr^?1 (not including soil respiration), and was close to BCB from forest inventories (1.15 t C ha^?1 yr^?1). Indirect human impact resulted in 33% increase in modeled biomass carbon accumulation in coniferous forests in Thuringia during the last century. From the forest inventory data we estimated the legacy effect of the age‐class distribution to account for 17% of the inventory‐based sink. Isolating the environmental change effects showed that these effects can be large in a long‐term, managed conifer forest.     

Estimating carbon fixation of plant organs for afforestation monitoring using a process‐based ecosystem model and ecophysiological parameter optimization

Afforestation projects for mitigating CO₂ emissions require to monitor the carbon fixation and plant growth as key indicators. We proposed a monitoring method for predicting carbon fixation in afforestation projects, combining a process‐based ecosystem model and field data and addressed the uncertainty of predicted carbon fixation and ecophysiological characteristics with plant growth. Carbon pools were simulated using the Biome‐BGC model tuned by parameter optimization using measured carbon density of biomass pools on an 11‐year‐old Eucommia ulmoides plantation on Loess Plateau, China. The allocation parameters fine root carbon to leaf carbon (FRC:LC) and stem carbon to leaf carbon (SC:LC), along with specific leaf area (SLA) and maximum stomatal conductance (g_smax) strongly affected aboveground woody (AC) and leaf carbon (LC) density in sensitivity analysis and were selected as adjusting parameters. We assessed the uncertainty of carbon fixation and plant growth predictions by modeling three growth phases with corresponding parameters: (i) before afforestation using default parameters, (ii) early monitoring using parameters optimized with data from years 1 to 5, and (iii) updated monitoring at year 11 using parameters optimized with 11‐year data. The predicted carbon fixation and optimized parameters differed in the three phases. Overall, 30‐year average carbon fixation rate in plantation (AC, LC, belowground woody parts and soil pools) was ranged 0.14–0.35 kg‐C m^?2 y^?1 in simulations using parameters of phases (i)–(iii). Updating parameters by periodic field surveys reduced the uncertainty and revealed changes in ecophysiological characteristics with plant growth. This monitoring method should support management of afforestation projects by carbon fixation estimation adapting to observation gap, noncommon species and variable growing conditions such as climate change, land use change.     

Carbon cycling and storage in world forests: biome patterns related to forest age

Forest age, which is affected by stand‐replacing ecosystem disturbances (such as forest fires, harvesting, or insects), plays a distinguishing role in determining the distribution of carbon (C) pools and fluxes in different forested ecosystems. In this synthesis, net primary productivity (NPP), net ecosystem productivity (NEP), and five pools of C (living biomass, coarse woody debris, organic soil horizons, soil, and total ecosystem) are summarized by age class for tropical, temperate, and boreal forest biomes. Estimates of variability in NPP, NEP, and C pools are provided for each biome‐age class combination and the sources of variability are discussed. Aggregated biome‐level estimates of NPP and NEP were higher in intermediate‐aged forests (e.g., 30–120 years), while older forests (e.g., >120 years) were generally less productive. The mean NEP in the youngest forests (0–10 years) was negative (source to the atmosphere) in both boreal and temperate biomes (?0.1 and –1.9 Mg C ha^?1 yr^?1, respectively). Forest age is a highly significant source of variability in NEP at the biome scale; for example, mean temperate forest NEP was ?1.9, 4.5, 2.4, 1.9 and 1.7 Mg C ha^?1 yr^?1 across five age classes (0–10, 11–30, 31–70, 71–120, 121–200 years, respectively). In general, median NPP and NEP are strongly correlated (R²=0.83) across all biomes and age classes, with the exception of the youngest temperate forests. Using the information gained from calculating the summary statistics for NPP and NEP, we calculated heterotrophic soil respiration (R_h) for each age class in each biome. The mean R_h was high in the youngest temperate age class (9.7 Mg C ha^?1 yr^?1) and declined with age, implying that forest ecosystem respiration peaks when forests are young, not old. With notable exceptions, carbon pool sizes increased with age in all biomes, including soil C. Age trends in C cycling and storage are very apparent in all three biomes and it is clear that a better understanding of how forest age and disturbance history interact will greatly improve our fundamental knowledge of the terrestrial C cycle.     

Modeling analysis of primary controls on net ecosystem productivity of seven boreal and temperate coniferous forests across a continental transect

Process‐based models are effective tools to synthesize and/or extrapolate measured carbon (C) exchanges from individual sites to large scales. In this study, we used a C‐ and nitrogen (N)‐cycle coupled ecosystem model named CN‐CLASS (Carbon Nitrogen‐Canadian Land Surface Scheme) to study the role of primary climatic controls and site‐specific C stocks on the net ecosystem productivity (NEP) of seven intermediate‐aged to mature coniferous forest sites across an east–west continental transect in Canada. The model was parameterized using a common set of parameters, except for two used in empirical canopy conductance–assimilation, and leaf area–sapwood relationships, and then validated using observed eddy covariance flux data. Leaf Rubisco‐N dynamics that are associated with soil–plant N cycling, and depend on canopy temperature, enabled the model to simulate site‐specific gross ecosystem productivity (GEP) reasonably well for all seven sites. Overall GEP simulations had relatively smaller differences compared with observations vs. ecosystem respiration (RE), which was the sum of many plant and soil components with larger variability and/or uncertainty associated with them. Both observed and simulated data showed that, on an annual basis, boreal forest sites were either carbon‐neutral or a weak C sink, ranging from 30 to 180 g C m^?2 yr^?1; while temperate forests were either a medium or strong C sink, ranging from 150 to 500 g C m^?2 yr^?1, depending on forest age and climatic regime. Model sensitivity tests illustrated that air temperature, among climate variables, and aboveground biomass, among major C stocks, were dominant factors impacting annual NEP. Vegetation biomass effects on annual GEP, RE and NEP showed similar patterns of variability at four boreal and three temperate forests. Air temperature showed different impacts on GEP and RE, and the response varied considerably from site to site. Higher solar radiation enhanced GEP, while precipitation differences had a minor effect. Magnitude of forest litter content and soil organic matter (SOM) affected RE. SOM also affected GEP, but only at low levels of SOM, because of low N mineralization that limited soil nutrient (N) availability. The results of this study will help to evaluate the impact of future climatic changes and/or forest C stock variations on C uptake and loss in forest ecosystems growing in diverse environments.     

基于数据同化的哈佛森林地区水、碳通量模拟

模型模拟和站点观测是陆地生态系统水、碳循环研究最主要的两种手段,但各有优势和不足,若二者相互结合,则能更准确地反映生态系统水、碳通量的动态变化.数据同化为模型与观测结合提供了一条有效的途径.本文采用哈佛森林环境监测站相关数据,利用集合卡曼滤波同化算法,将实测叶面积指数(LAI)和遥感LAI同化进入Biome BGC模型中,对该地区水、碳通量进行模拟.结果表明：与未同化模拟相比,将1998、1999和2006年实测LAI数据同化后,模型模拟碳通量(NEE)与通量观测NEE的决定系数(R²)平均提升8.4%;蒸散发(ET)的R²平均提升10.6%;NEE的绝对误差和(SAE)和均方根误差(RMSE)平均下降17.7%和21.2%,ET的SAE和RMSE平均下降26.8%和28.3%.将2000-2004年MODIS LAI 产品与模型同化后,NEE、ET模拟值与观测值间的R²分别提升7.8%和4.7%;NEE的SAE和 RMSE分别下降21.9%和26.3%,ET的SAE和 RMSE分别下降24.5%和25.5%.无论实测LAI还是遥感观测LAI,同化进入模型都能不同程度地提高水碳通量的模拟精度.     

Biomass and carbon accumulation in a fire chronosequence of a seasonally dry tropical forest

Seasonally dry tropical forests (SDTF) are a widely distributed vegetation type in the tropics, characterized by seasonal rainfall with several months of drought when they are subject to fire. This study is one of the first attempts to quantify above- and belowground biomass (AGB and BGB) and above- and belowground carbon (AGC and BGC) pools to calculate their recovery after fire, using a chronosequence approach (six forests that ranged form 1 to 29 years after fire and mature forest). We quantified AGB and AGC pools of trees, lianas, palms, and seedlings, and BGB and BGC pools (Oi, Oe, Oa soil horizons, and fine roots). Total AGC ranged from 0.05 to nearly 72 Mg C ha⁻¹, BGC from 21.6 to nearly 85 Mg C ha⁻¹, and total ecosystem carbon from 21.7 to 153.5 Mg C ha⁻¹; all these pools increased with forest age. Nearly 50% of the total ecosystem carbon was stored in the Oa horizon of mature forests, and up to 90% was stored in the Oa-horizon of early successional SDTF stands. The soils were shallow with a depth of <20 cm at the study site. To recover values similar to mature forests, BGC and BGB required <19 years with accumulation rates greater than 20 Mg C ha⁻¹ yr⁻¹, while AGB required 80 years with accumulation rates nearly 2.5 Mg C ha⁻¹ yr⁻¹. Total ecosystem biomass and carbon required 70 and 50 years, respectively, to recover values similar to mature forests. When belowground pools are not included in the calculation of total ecosystem biomass or carbon recovery, we estimated an overestimation of 10 and 30 years, respectively.     

Ecosystem carbon storage does not vary with mean annual temperature in Hawaiian tropical montane wet forests

Theory and experiment agree that climate warming will increase carbon fluxes between terrestrial ecosystems and the atmosphere. The effect of this increased exchange on terrestrial carbon storage is less predictable, with important implications for potential feedbacks to the climate system. We quantified how increased mean annual temperature (MAT) affects ecosystem carbon storage in above‐ and belowground live biomass and detritus across a well‐constrained 5.2 °C MAT gradient in tropical montane wet forests on the Island of Hawaii. This gradient does not systematically vary in biotic or abiotic factors other than MAT (i.e. dominant vegetation, substrate type and age, soil water balance, and disturbance history), allowing us to isolate the impact of MAT on ecosystem carbon storage. Live biomass carbon did not vary predictably as a function of MAT, while detrital carbon declined by ~14 Mg of carbon ha^?1 for each 1 °C rise in temperature – a trend driven entirely by coarse woody debris and litter. The largest detrital pool, soil organic carbon, was the most stable with MAT and averaged 48% of total ecosystem carbon across the MAT gradient. Total ecosystem carbon did not vary significantly with MAT, and the distribution of ecosystem carbon between live biomass and detritus remained relatively constant across the MAT gradient at ~44% and ~56%, respectively. These findings suggest that in the absence of alterations to precipitation or disturbance regimes, the size and distribution of carbon pools in tropical montane wet forests will be less sensitive to rising MAT than predicted by ecosystem models. This article also provides needed detail on how individual carbon pools and ecosystem‐level carbon storage will respond to future warming.     

Understanding global teleconnections of climate to regional model estimates of Amazon ecosystem carbon fluxes

We have investigated global teleconnections of climate to regional satellite‐driven observations for prediction of Amazon ecosystem production, in the form of monthly estimates of net carbon exchange over the period 1982–1998 from the NASA–CASA (Carnegie–Ames–Stanford) biosphere model. This model is driven by observed surface climate and monthly estimates of vegetation leaf area index (LAI) and fraction of absorbed PAR (fraction of photosynthetically active radiation, FPAR) generated from the NOAA satellite advanced very high‐resolution radiometer (AVHRR) and similar sensors. Land surface AVHRR data processing using modified moderate‐resolution imaging spectroradiometer radiative transfer algorithms includes improved calibration for intra‐ and intersensor variations, partial atmospheric correction for gaseous absorption and scattering, and correction for stratospheric aerosol effects associated with volcanic eruptions. Results from our analysis suggest that anomalies of net primary production and net ecosystem production predicted from the NASA–CASA model over large areas of the Amazon region east of 60°W longitude are strongly correlated with the Southern Oscillation index. Extensive areas of the south‐central Amazon show strong linkages of the FPAR and the NASA–CASA anomaly record to the Arctic Oscillation index, which help confirm a strong relation to southern Atlantic climate anomalies, with associated impacts on Amazon rainfall patterns. Processes are investigated for these teleconnections of global climate to Amazon ecosystem carbon fluxes and regional land surface climate.     

Canopy duration has little influence on annual carbon storage in the deciduous broad leaf forest

Vegetation phenology, the study of the timing and length of the terrestrial growing season and its connection to climate, is increasingly important in integrated Earth system science. Phenological variability is an excellent barometer of short‐ and long‐term climatic variability, strongly influences surface meteorology, and may influence the carbon cycle. Here, using the 1895–1993 Vegetation/Ecosystem Modelling and Analysis dataset and the Biome‐BGC terrestrial ecosystem model, we investigated the relationship between phenological metrics and annual net ecosystem exchange (NEE) of carbon. For the 1167 deciduous broad leaf forest pixels, we found that NEE was extremely weakly related to canopy duration (days from leaf appearance to complete leaf fall). Longer canopy duration, did, however, sequester more carbon if warm season precipitation was above average. Carbon uptake period (number of days with net CO₂ uptake from the atmosphere), which integrates the influence of all ecosystem states and processes, was strongly related to NEE. Results from the Harvard Forest eddy‐covariance site supported our findings. Such dramatically different results from two definitions of ‘growing season length’ highlight the potential for confusion among the many disciplines engaged in phenological research.     

Changes in carbon storage and fluxes in a chronosequence of ponderosa pine

Forest development following stand‐replacing disturbance influences a variety of ecosystem processes including carbon exchange with the atmosphere. On a series of ponderosa pine (Pinius ponderosa var. Laws.) stands ranging from 9 to> 300 years in central Oregon, USA, we used biological measurements to estimate carbon storage in vegetation and soil pools, net primary productivity (NPP) and net ecosystem productivity (NEP) to examine variation with stand age. Measurements were made on plots representing four age classes with three replications: initiation (I, 9–23 years), young (Y, 56–89 years), mature (M, 95–106 years), and old (O, 190–316 years) stands typical of the forest type in the region. Net ecosystem productivity was lowest in the I stands (?124 g C m^?2 yr^?1), moderate in Y stands (118 g C m^?2 yr^?1), highest in M stands (170 g C m^?2 yr^?1), and low in the O stands (35 g C m^?2 yr^?1). Net primary productivity followed similar trends, but did not decline as much in the O stands. The ratio of fine root to foliage carbon was highest in the I stands, which is likely necessary for establishment in the semiarid environment, where forests are subject to drought during the growing season (300–800 mm precipitation per year). Carbon storage in live mass was the highest in the O stands (mean 17.6 kg C m^?2). Total ecosystem carbon storage and the fraction of ecosystem carbon in aboveground wood mass increased rapidly until 150–200 years, and did not decline in older stands. Forest inventory data on 950 ponderosa pine plots in Oregon show that the greatest proportion of plots exist in stands ～ 100 years old, indicating that a majority of stands are approaching maximum carbon storage and net carbon uptake. Our data suggests that NEP averages ～ 70 g C m^?2 year^?1 for ponderosa pine forests in Oregon. About 85% of the total carbon storage in biomass on the survey plots exists in stands greater than 100 years, which has implications for managing forests for carbon sequestration. To investigate variation in carbon storage and fluxes with disturbance, simulation with process models requires a dynamic parameterization for biomass allocation that depends on stand age, and should include a representation of competition between multiple plant functional types for space, water, and nutrients.     

DROUGHT EFFECTS ON CARBON EXCHANGE IN A SUBTROPICAL CONIFEROUS PLANTATION IN CHINA

 干旱对陆地生态系统的影响已成为全球变化研究的焦点问题之一。该研究基于生态系统过程模型——CEVSA2, 结合涡度相关通量观测, 分析了不同程度干旱对亚热带人工针叶林碳交换的影响及其关键控制因素。结果表明: 1)干旱使生态系统碳交换显著下降, 2003和2004年的干旱使得年净生态系统生产力(Net ecosystem production, NEP)相比无干旱影响情景的模拟结果分别减少了63%和47%; 2)光合和呼吸对干旱具有不同的响应, 干旱时光合的下降比呼吸更为显著, 这导致了NEP的显著下降; 3)当饱和水气压差(Vapor pressure deficit, VPD)达到1.5 kPa以上时, 生态系统的光合、呼吸和净碳吸收均开始下降, 当VPD大于2.5 kPa、土壤相对含水量(土壤含水量/土壤饱和含水量)(Relative soil water content, RSW)低于40%时, 生态系统的碳收支由碳汇转为碳源; 4)土壤干旱是造成碳交换下降的主要驱动因素, 对年NEP下降的平均贡献率为46%, 而大气干旱的贡献率仅为4%。     

Patterns of NPP,GPP, respiration,and NEP during boreal forest succession

We combined year‐round eddy covariance with biometry and biomass harvests along a chronosequence of boreal forest stands that were 1, 6, 15, 23, 40, ～74, and ～154 years old to understand how ecosystem production and carbon stocks change during recovery from stand‐replacing crown fire. Live biomass (C_live) was low in the 1‐ and 6‐year‐old stands, and increased following a logistic pattern to high levels in the 74‐ and 154‐year‐old stands. Carbon stocks in the forest floor (C_{forest floor}) and coarse woody debris (C_CWD) were comparatively high in the 1‐year‐old stand, reduced in the 6‐ through 40‐year‐old stands, and highest in the 74‐ and 154‐year‐old stands. Total net primary production (TNPP) was reduced in the 1‐ and 6‐year‐old stands, highest in the 23‐ through 74‐year‐old stands and somewhat reduced in the 154‐year‐old stand. The NPP decline at the 154‐year‐old stand was related to increased autotrophic respiration rather than decreased gross primary production (GPP). Net ecosystem production (NEP), calculated by integrated eddy covariance, indicated the 1‐ and 6‐year‐old stands were losing carbon, the 15‐year‐old stand was gaining a small amount of carbon, the 23‐ and 74‐year‐old stands were gaining considerable carbon, and the 40‐ and 154‐year‐old stands were gaining modest amounts of carbon. The recovery from fire was rapid; a linear fit through the NEP observations at the 6‐ and 15‐year‐old stands indicated the transition from carbon source to sink occurred within 11–12 years. The NEP decline at the 154‐year‐old stand appears related to increased losses from C_live by tree mortality and possibly from C_{forest floor} by decomposition. Our findings support the idea that NPP, carbon production efficiency (NPP/GPP), NEP, and carbon storage efficiency (NEP/TNPP) all decrease in old boreal stands.