Restoration of ecosystems for biodiversity and carbon sequestration: Simulating growth dynamics of brigalow vegetation communities in Australia

Restoration of abandoned and degraded ecosystems through enhanced management of mature remnant patches and naturally regenerating (regrowth) forests is currently being used in the recovery of ecosystems for biodiversity protection and carbon sequestration. Knowledge of long-term dynamics of these ecosystems is often very limited. Vegetation models that examine long-term forest growth and succession of uneven aged, mixed-species forest ecosystems are integral to the planning and assessment of the recovery process of biodiversity values and biomass accumulation. This paper examined the use of the Ecosystem Dynamics Simulator (EDS) in projecting growth dynamics of mature remnant brigalow forest communities and recovery process of regrowth brigalow thickets. We used data from 188 long-term monitored plots of remnant and regrowth forests measured between 1963 and 2010. In this study the model was parameterised for 34 tree and shrub species and tested with independent long-term measurements. The model closely approximated actual development trajectories of mature forests and regrowth thickets but some inaccuracies in estimating regeneration through asexual reproduction and mortality were noted as reflected in stem density projections of remnant plots that had a mean of absolute relative bias of 46.2 (±12.4)%. Changes in species composition in remnant forests were projected with a 10% error. Basal area values observed in all remnant plots ranged from 6 to 29 m² ha⁻¹ and EDS projections between 1966 and 2005 (39 years) were 68.2 (±10.9)% of the observed basal area. Projected live aboveground biomass of remnant plots had a mean of 93.5 (±5.9) t ha⁻¹ compared to a mean of 91.3 (±8.0) t ha⁻¹ observed in the plots. In regrowth thicket, the model produced satisfactory projections of tree density (91%), basal area (89%), height (87%) and aboveground biomass (84%) compared to the observed attributes. Basal area and biomass accumulation in 45-year-old regrowth plots was approximately similar to that in remnant forests but recovery of woody understorey was very slow. The model projected that it would take 95 years for the regrowth to thin down to similar densities observed in original or remnant brigalow forests. These results indicated that EDS can produce relatively accurate projections of growth dynamics of brigalow regrowth forests necessary for informing restoration planning and projecting biomass accumulation.     

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.     

Tradeoffs between forestry resource and conservation values under alternate policy regimes: A spatial analysis of the western Canadian boreal plains

An important element of resource management and conservation is an understanding of the tradeoffs between marketed products, such as timber, and measures of environmental quality, such as biodiversity. In this paper, we develop an integrated economic-ecological spatial optimization model that we then apply to evaluate alternate forest policies on a 560,000 km² study region of managed boreal forest in Alberta and British Columbia, Canada. The integrated model incorporates dynamic forest sector harvesting, current levels of oil and gas sector development, coarse-filter or habitat-based old forest indicators, a set of empirical forest bird abundance models, and statistical models of the natural and current fire regimes. Using our integrated model, economic tradeoff curves, or production possibility frontiers, are developed to illustrate the cost of achieving coarse-filter targets by a set time (50 years) within a 100-year time horizon. We found levels of ecological indicators and economic returns from the timber industry could both be increased if spatial constraints imposed by the current policy environment were relaxed; other factors being equal, this implies current policy should be revised. We explore the production possibility frontier's relationship to the range of natural variation of old forest habitat, and show how this range can be used to guide choices of preferred locations along the frontier. We also show that coarse-filter constraints on the abundance of certain habitat elements are sufficient to satisfy some fine-filter objectives, expressed as the predicted abundances of various species of songbirds.     

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.     

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.     

Dynamics of surface albedo of a boreal forest and its simulation

Surface albedo determines the distribution of solar radiation between the earth's surface and the atmosphere. It affects the global climate directly by altering surface energy balance, and indirectly by controlling ecosystem processes and greenhouse gas exchange. In this study, a land surface albedo model was constructed based on the gap probability approach for ray tracing and the basic optical parameters of ecosystem elements. The model was applied to a boreal deciduous forest and results were compared with field measurements. Results show that seasonal and diurnal albedo dynamics were well simulated by the model. The standard deviation between the simulated and measured reflected radiation was 2.5–5.0 W m⁻² in different seasons. The model also provided an insight into the relationships between surface albedo and radiation components (direct versus diffuse), solar zenith angle, and different wave bands. Model sensitivity analyses show that the surface albedo in winter is very sensitive to the forest wood area index for this boreal aspen forest, suggesting that accurate estimates of wood area index are necessary to improve the accuracy of surface albedo simulation in leafless seasons.     

Spatio-temporal modelling of ground vegetation development in mountain spruce forests

More complex models of forest ecosystems are required to understand how land-cover changes can impact vegetation dynamics and spatial pattern. In order to document spatio-temporal modelling abilities, the observations conducted in the declined climax mountain Norway spruce forest during the recovery period (1995-2006) are used for simulation and spatial analysis in the GIS environment. The developed spatio-temporal model is used for simulation of forest vegetation dynamics in a mountain spruce forest in the framework of regeneration processes after stress from air pollution. In order to explore the spatial and temporal phenomena of regeneration processes, the spatio-temporal model is based on a large set of ordinary differential equations that solve dynamic processes in sets of microsites arranged in grids for each ground vegetation species and each age group of Norway spruce seedlings. The spatial extent of the explored site is composed of a set of 50 × 50 microsites. Each microsite is represented by a square with dimensions of 1 m × 1 m. The presented simulation studies are mainly focused on seedlings from the seed year 1992, in order to explore the longest monitored time series of survival. It is based on exponential growth models that are related to the environmental conditions for each microsite. The canopy gaps based on estimates of the local crown projected area, the soil type layer, and the dominant grass density are used to provide case simulation studies. The first case study simulates the influence of microsite positions in relation to the local tree crown projections on the survival of spruce seedlings. It is assumed that the density of the trees is the main factor that determines the light and heat supply to the ground level of the Norway spruce seedlings. The second case study extends the previous study to include terms that determine the growth ratio in dependence on the crown projection area. The third case study provides further extensions in order to simulate growth ratio relations to the local soil type. The fourth case study demonstrates the local influence of the dominant grasses, such as Avenella flexuosa and Calamagrostis villosa, on the natural regeneration of Norway spruce. Starting from the conditions at the sites before the recovery period, the case simulation studies are able to project the short-term succession for a regeneration decade and the approximate long-term development. In addition to the standard simulation procedures based on solution of ordinary differential equations, spatio-temporal modelling in the GIS environment is able to provide spatial data management, analysis and visualization of the data.     

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.     

Spatio-temporal analysis of alpine ecotones: A spatial explicit model targeting altitudinal vegetation shifts

There is general agreement in literature that Alpine vegetation belt ecotones have shown a trend of upward migration in the last few decades. Despite the potential of such shifts as indicators of global change effects in mountain ecosystems, there are relatively few works focused on their assessment in a systematic and spatially explicit way. In this work our aim is to quantify the altitudinal shifts and analyse the spatial pattern dynamics of mountain ecotones. We developed a novel procedure to delineate the current and former state of three characteristic mountain ecotones, which we formalised as forest, tree and tundra lines. Our approach is based on the recognition of altitudinal extreme outposts identified with ecotone locations at a slope scale. The integration of multi-temporal datasets allows the identification and quantification of altitudinal advances and retreats in the outpost locations for a given period. We tested the method in a section of the Italian Alps for the period 1957-2003. Results show a general trend of an increase in altitude for the three ecotones, despite the occurrence of occasional decreases. We estimate decadal altitude increments of 25 m for forest line, 13 m for treeline and 11 m for tundra line. We also identified changes in ecotone spatial morphology between the two dates, with significant implications in connectivity and colonisation dynamics.     

Impact of climate change on soil moisture dynamics in Brandenburg with a focus on nature conservation areas

Global warming impacts the water cycle not only by changing regional precipitation levels and temporal variability, but also by affecting water flows and soil moisture dynamics. In Brandenburg, increasing average annual temperature and decreasing precipitation in summer have already been observed. For this study, past trends and future effects of climate change on soil moisture dynamics in Brandenburg were investigated, considering regional and specific spatial impacts. Special Areas of Conservation (SACs) were focused on in particular. A decreasing trend in soil water content was shown for the past by analyzing simulation results from 1951 to 2003 using the integrated ecohydrological model SWIM [Krysanova, V., Müller-Wohlfeil, D.-I., Becker, A., 1998. Development and test of a spatially distributed hydrological/water quality model for mesoscale watersheds. Ecol. Model. 106, 261–289]. The trend was statistically significant for some areas, but not for the entire region. Simulated soil water content was particularly low in the extremely dry year 2003. Comparisons of simulated trends in soil moisture dynamics with trends in the average annual Palmer Drought Severity Index for the region showed largely congruent patterns, though the modeled soil moisture trends are characterized by a much higher spatial resolution. Regionally downscaled climate change projections representing the range between wetter and drier realizations were used to evaluate future trends of available soil water. A further decrease of average available soil water ranging from −4% to −15% was projected for all climate realizations up to the middle of the 21st century. An average decrease of more than 25 mm was simulated for 34% of the total area in the dry realization. Available soil water contents in SACs were generally higher and trends in soil moisture dynamics were lower mainly due to their favorable edaphic conditions. Stronger absolute and relative changes in the simulated trends for the past and future were shown for SACs within Brandenburg than for the state as a whole, indicating a high level of risk for many wetland areas. Nonetheless, soil water content in SACs is expected to remain higher than average under climate change conditions as well, and SACs therefore have an important buffer function under the projected climate change. They are thus essential for local climate and water regulation and their status as protected areas in Brandenburg should be preserved.     

The sensitivity of carbon sequestration to harvesting and climate conditions in a temperate cypress forest: Observations and modeling

The role of disturbance and climate factors in determining the forest carbon balance was investigated at a Japanese cypress forest in central Japan with eddy flux measurements, tree-ring analyses, and a terrestrial biosphere model. The forest was established as a plantation after intermittent harvesting and replanting between 1959 and 1977, and acted as a strong carbon sink of approximately 500 g C m⁻² year⁻¹ for the measurement years between 2001 and 2007. A terrestrial biosphere model, BIOME-BGC, was validated using the eddy flux data at daily to interannual timescales, and the tree-ring width data at interannual to decadal timescales. According to the model simulation, during the observation period 270 ± 55 g C m⁻² year⁻¹ was additionally sequestered due to the indirect effects of the harvesting and planting, whereas the increase of CO₂ concentration and the change in climate increased the sink of 110 ± 40 and 30 ± 80 g C m⁻² year⁻¹, respectively. The model simulation shows that the forest is now recovering from harvesting, and that harvesting is a more important determinant of the current carbon sink than either interannual climate anomalies or increased atmospheric CO₂ concentration. We found that harvesting with long rotation length could be effective management activity in order to increase carbon sequestration, if the harvested timber is converted into products with long lifecycles.     

The application of artificial neural networks to flow and phosphorus dynamics in small streams on the Boreal Plain,with emphasis on the role of wetlands

Historically, management strategies in Canada's boreal forest have focused on forest polygons and terrestrial biodiversity to address ecological considerations in forest management. The Forest Watershed and Riparian Disturbance (FORWARD) project examines the problem from a watershed perspective rather than a forest polygon viewpoint. The main objective of this study was to devise an artificial neural network (ANN) modeling tool that can predict flow and total phosphorus (TP) concentration for ungauged watersheds (where daily flow is not monitored). This dictates that all inputs should be easily accessed via a public domain database, like the Environment Canada weather database, without the need to install flow gauges in each modeled watershed. Daily flow and TP concentration for two of the project watersheds were modeled using ANNs. The two watersheds (1A Creek, 5.1 km² and Willow Creek, 15.6 km²) were chosen to reflect variations in wetland area and composition in the study area. Flow was modeled with a feed-forward multilayer perceptron ANN trained with the error back-propagation algorithm. Simulated values for flow were then used, as inputs, to model TP concentration using the same neural networks algorithm. One hidden layer with three slabs; each operating with a different activation function was utilized to simulate the conceptual differences between base flow, snowmelt, and storm events. Time domain analysis was conducted to identify possible model time-lagged inputs reflecting the time dependency of the modeled variables. Spectral analysis was used to address data hystereses. Our results highlight the capabilities of ANN in modeling complex ecosystems and highly correlated variables. Results also indicated that more research towards the phosphorus dynamics in wetlands is required to better represent the impact of wetland area and composition on the water-phase phosphorus in ANN modeling.     

Modelling forest management within a global vegetation model—Part 1: Model structure and general behaviour

This article describes a new forest management module (FMM) that explicitly simulates forest stand growth and management within a process-based global vegetation model (GVM) called ORCHIDEE. The net primary productivity simulated by ORCHIDEE is used as an input to the FMM. The FMM then calculates stand and management characteristics such as stand density, tree size distribution, tree growth, the timing and intensity of thinnings and clear-cuts, wood extraction and litter generated after thinning. Some of these variables are then fed back to ORCHIDEE. These computations are made possible with a distribution-based modelling of individual tree size. The model derives natural mortality from the relative density index (rdi), a competition index based on tree size and stand density. Based on the common forestry management principle of avoiding natural mortality, a set of rules is defined to calculate the recurrent intensity and frequency of forestry operations during the stand lifetime. The new-coupled model is called ORCHIDEE-FM (forest management).The general behaviour of ORCHIDEE-FM is analysed for a broadleaf forest in north-eastern France. Flux simulation throughout a forest rotation compare well with the literature values, both in absolute values and dynamics.Results from ORCHIDEE-FM highlight the impact of forest management on ecosystem C-cycling, both in terms of carbon fluxes and stocks. In particular, the average net ecosystem productivity (NEP) of 225 gC m⁻² year⁻¹ is close to the biome average of 311 gC m⁻² year⁻¹. The NEP of the “unmanaged” case is 40% lower, leading us to conclude that management explains 40% of the cumulated carbon sink over 150 years. A sensitivity analysis reveals 4 major avenues for improvement: a better determination of initial conditions, an improved allocation scheme to explain age-related decline in productivity, and an increased specificity of both the self-thinning curve and the biomass-diameter allometry.     

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.     

Response of tree growth to a changing climate in boreal central Canada: A comparison of empirical,process-based,and hybrid modelling approaches

The impact of 2 × CO₂ driven climate change on radial growth of boreal tree species Pinus banksiana Lamb., Populus tremuloides Michx. and Picea mariana (Mill.) BSP growing in the Duck Mountain Provincial Forest of Manitoba (DMPF), Canada, is simulated using empirical and process-based model approaches. First, empirical relationships between growth and climate are developed. Stepwise multiple-regression models are conducted between tree-ring growth increments (TRGI) and monthly drought, precipitation and temperature series. Predictive skills are tested using a calibration–verification scheme. The established relationships are then transferred to climates driven by 1× and 2 × CO₂ scenarios using outputs from the Canadian second-generation coupled global climate model. Second, empirical results are contrasted with process-based projections of net primary productivity allocated to stem development (NPP_s). At the finest scale, a leaf-level model of photosynthesis is used to simulate canopy properties per species and their interaction with the variability in radiation, temperature and vapour pressure deficit. Then, a top-down plot-level model of forest productivity is used to simulate landscape-level productivity by capturing the between-stand variability in forest cover. Results show that the predicted TRGI from the empirical models account for up to 56.3% of the variance in the observed TRGI over the period 1912–1999. Under a 2 × CO₂ scenario, the predicted impact of climate change is a radial growth decline for all three species under study. However, projections obtained from the process-based model suggest that an increasing growing season length in a changing climate could counteract and potentially overwhelm the negative influence of increased drought stress. The divergence between TRGI and NPP_s simulations likely resulted, among others, from assumptions about soil water holding capacity and from calibration of variables affecting gross primary productivity. An attempt was therefore made to bridge the gap between the two modelling approaches by using physiological variables as TRGI predictors. Results obtained in this manner are similar to those obtained using climate variables, and suggest that the positive effect of increasing growing season length would be counteracted by increasing summer temperatures. Notwithstanding uncertainties in these simulations (CO₂ fertilization effect, feedback from disturbance regimes, phenology of species, and uncertainties in future CO₂ emissions), a decrease in forest productivity with climate change should be considered as a plausible scenario in sustainable forest management planning of the DMPF.     

Spatially distributed modeling of the long-term carbon balance of a boreal landscape

Spatially and temporally distributed information on the sizes of biomass carbon (C) pools (BCPs) and soil C pools (SCPs) is vital for improving our understanding of biosphere-atmosphere C fluxes. Because the sizes of C pools result from the integrated effects of primary production, age-effects, changes in climate, atmospheric CO₂ concentration, N deposition, and disturbances, a modeling scheme that interactively considers these processes is important. We used the InTEC model, driven by various spatio-temporal datasets to simulate the long-term C-balance in a boreal landscape in eastern Canada. Our results suggested that in this boreal landscape, mature coniferous stands had stabilized their productivity and fluctuated as a weak C-sink or C-source depending on the interannual variations in hydrometeorological factors. Disturbed deciduous stands were larger C-sinks (NEP₂₀₀₄ = 150 gC m⁻² yr⁻¹) than undisturbed coniferous stands (e.g. NEP₂₀₀₄ = 8 gC m⁻² yr⁻¹). Wetlands had lower NPP but showed temporally consistent C accumulation patterns. The simulated spatio-temporal patterns of BCPs and SCPs were unique and reflected the integrated effects of climate, plant growth and atmospheric chemistry besides the inherent properties of the C pool themselves. The simulated BCPs and SCPs generally compared well with the biometric estimates (BCPs: r = 0.86, SCPs: r = 0.84). The largest BCP biases were found in recently disturbed stands and the largest SCP biases were seen in locations where moss necro-masses were abundant. Reconstructing C pools and C fluxes in the ecosystem in such a spatio-temporal manner could help reduce the uncertainties in our understanding of terrestrial C-cycle.     

Uncertainty in future soil carbon trends at a central U.S. site under an ensemble of GCM scenario climates

This study uses DAYCENT model to investigate the sensitivity of soil organic carbon (SOC) at an intensely cultivated site in the U.S. Midwest under an ensemble of scenario climates predicted by IPCC models. The model ensemble includes three IPCC models (Canadian, French, German), three emission scenarios (B1, A1B, A2) and three time periods (late 20th, mid-21st, late 21st century). DAYCENT shows that SOC at the site would decline by 0.3-2.6 kg m⁻² (5-35%) depending on the models and scenarios from late 20th to mid-21st century despite a larger increase of future net primary production (NPP) than respiration. The future SOC decrease is mostly attributable to harvest loss. The wide spread in future SOC decline rates are in part because SOC decrease (by respiration) is directly proportional to SOC itself. Any uncertainty in absolute SOC in DAYCENT would translate directly into its trend, unlike other variables such as temperature whose trends are independent of their values themselves, contrasting the reliability of SOC trend with temperature change.     

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.     

Estimating the potential impact of vegetation on the water cycle requires accurate soil water parameter estimation

It is well known that vegetation dynamics at the catchment scale depends on the prevailing weather and soil moisture conditions. Soil moisture, however, is not equally distributed in space due to differences in topography, weather patterns, soil properties and the type and amount of vegetation cover. To elucidate the complex interaction between vegetation and soil moisture, the dynamic vegetation model LPJ-GUESS (Smith et al., 2001), which provides estimations of vegetation dynamics, but does not consider lateral water fluxes was coupled with the hydrological TOPMODEL (cf. Beven, 2001) in order to be able to evaluate the importance of these lateral fluxes. The new model LG-TM was calibrated and validated in two climatically different mountain catchments. The estimations of runoff were good, when monthly and weekly time scales were considered, although the low flow periods at winter time were somewhat underestimated. The uncertainty in the climate induced change vegetation carbon storage caused by the uncertainty in soil parameters was up to 3-5 kg C m⁻² (depending on elevation and catchment), compared to the total change in vegetation carbon storage of 5-9 kg C m⁻². Therefore accurate estimates of the parameters influencing the water holding capacity of the soil, for example depth and porosity, are necessary when estimating future changes in vegetation carbon storage. Similarly, changes in plant transpiration due to climatic changes could be almost double as high (88 mm m⁻²) in the not calibrated model compared to the new model version (ca 50 mm m⁻² transpiration change). The uncertainties in these soil properties were found to be more important than the lateral water exchange between grid cells, even in steep topography at least for the temporal and spatial resolution used here.     

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.