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.     

Structural changes in lake functioning induced from nutrient loading and climate variability

Climate variability is increasingly recognized as an important regulatory factor, capable of influencing the structural properties of aquatic ecosystems. Lakes appear to be particularly sensitive to the ecological impacts of climate variability, and several long time series have shown a close coupling between climate, lake thermal properties and individual organism physiology, population abundance, community structure, and food web dynamics. Thus, understanding the complex interplay among meteorological forcing, hydrological variability, and ecosystem functioning is essential for improving the credibility of model-based water resources/fisheries management. Our objective herein is to examine the relative importance of the ecological mechanisms underlying plankton seasonal variability in Lake Washington, Washington State (USA), over a 35-year period (1964–1998). Our analysis is founded upon an intermediate complexity plankton model that is used to reproduce the limiting nutrient (phosphate)–phytoplankton–zooplankton–detritus (particulate phosphorus) dynamics in the lake. Model parameterization is based on a Bayesian calibration scheme that offers insights into the degree of information the data contain about model inputs and allows obtaining predictions along with uncertainty bounds for modeled output variables. The model accurately reproduces the key seasonal planktonic patterns in Lake Washington and provides realistic estimates of predictive uncertainty for water quality variables of environmental management interest. A principal component analysis of the annual estimates of the underlying ecological processes highlighted the significant role of the phosphorus recycling stemming from the zooplankton excretion on the planktonic food web variability. We also identified a moderately significant signature of the local climatic conditions (air temperature) on phytoplankton growth (r = 0.41), herbivorous grazing (r = 0.38), and detritus mineralization (r = 0.39). Our study seeks linkages with the conceptual food web model proposed by Hampton et al. [Hampton, S.E., Scheuerell, M.D., Schindler, D.E., 2006b. Coalescence in the Lake Washington story: interaction strengths in a planktonic food web. Limnol. Oceanogr. 51, 2042–2051.] to emphasize the “bottom-up” control of the Lake Washington plankton phenology. The posterior predictive distributions of the plankton model are also used to assess the exceedance frequency and confidence of compliance with total phosphorus (15 μg L⁻¹) and chlorophyll a (4 μg L⁻¹) threshold levels during the summer-stratified period in Lake Washington. Finally, we conclude by underscoring the importance of explicitly acknowledging the uncertainty in ecological forecasts to the management of freshwater ecosystems under a changing global environment.     

Testing hypotheses on shape and distribution of ecological response curves

Niche theory with hypotheses on shape and distribution of ecological response curves is used in the studies of resource sharing of competing plant species. Predictions based on theory should be applicable when, e.g., effects of competing species on the ecological tolerances are assessed or species’ diversity along a resource gradient is evaluated. We studied the ecological response curves of competing plant species along a resource gradient in boreal forests. The study was based on nation-wide soil and vegetation data collected from 455 sample plots on boreal forests in Finland. Species response curves along a soil fertility gradient (in terms of C/N ratio) were estimated using generalized additive models. Distribution of species optima and the relationship of niche width and skewness to the location of the optimum were analyzed with new bootstrap tests. The developed tests can account for the effects of truncation observed in the response curves of several species and for the uneven distribution of observations on the gradient.The estimated response curves of the major field layer species of boreal forests were not evenly distributed along soil C/N gradient. The density of optima peaked with relatively high nitrogen availability. Species with optima at low nitrogen availability had relatively broad realized niches. Niche width was negatively correlated with the density of optima. Species optima were packed and niches were narrow at high resource levels. This result suggests that a greater number of more specialized species can occur and interspecific competition decreases niche widths at high resource levels. Species were packed in the gradient where the C/N ratio was lower than 25, i.e., in conditions where nitrification can take place. This indicates that the majority of the vascular plants of boreal forests are favoured by the availability of NO₃⁻. Those few species thriving at high C/N ratios have broader realized niches.     

Using classification trees to analyze the impact of exotic species on the ecological assessment of polder lakes in Flanders, Belgium

Polder lakes in Flanders are stagnant waters that were flooded by the sea in the past. Several of these systems are colonized by exotic species, but have hardly been studied until present. The aim of the present study was: (1) to assess the influence of exotic macrobenthic species on the outcome of the Multimetric Macroinvertebrate Index Flanders (MMIF) and (2) to use classification trees for evaluating to what extent physical-chemical characteristics affect the presence of exotic species.In total, 27 mollusc and 10 macro-crustacean species were present in the monitored lakes of which respectively five and four were exotic. The exclusion of the exotic species from the MMIF resulted in a significant decline of this ecological index (−0.03 ± 0.04; p = 0.00). This elimination often resulted into a lower ecological water quality class and more samples were classified into the bad and poor ecological water quality classes.Single-target classification trees for Gammarus tigrinus and Potamopyrgus antipodarum were constructed, relating environmental parameters and ecological status (MMIF) to the occurrence of both exotic invasive species. The major advantages of using single-target classification trees are the transparency of the rule sets and the possibility to use relatively small datasets. However, this classification technique only predicts a single-target attribute and the trees of the different species are often hard to integrate and use for water managers. As a solution, a multi-target approach was used in the present study. Exotic molluscs and crustaceans communities were modelled based on environmental parameters and the ecological status (MMIF) using multi-target classification trees. Multi-target classification trees can be used in management planning and investment decisions as they can lead to integrated decisions for the whole set of exotic species and avoid the construction of many models for each individual species. These trees provide general insights concerning the occurrence patterns of individual crustaceans and molluscs in an integrated way.     

A novel method for mapping reefs and subtidal rocky habitats using artificial neural networks

Reefs and subtidal rocky habitats are sites of high biodiversity and productivity which harbour commercially important species of fish and invertebrates. Although the conservation management of reef associated species has been informed using species distribution models (SDM) and community based approaches, to date their use has been constrained to specific regions where the locality and spatial extent of reefs is well known. Much of the world's subtidal habitats remain either undiscovered or unmapped, including coasts of intense human use. Consequently, to facilitate a stronger understanding of species-environmental relationships there is an urgent need for a cost and time effective standard method to map reefs at fine spatial resolutions across broad geographical extents. We used bathymetric data (∼250 m resolution) to calculate the local slope and curvature of the seabed. We then constructed artificial neural networks (ANNs) to forecast the probability of reef occurrence within grid cells as a function of bathymetric and slope variables. Testing over an independent data set not used in training showed that ANNs were able to accurately predict the location of reefs for 86% of all grid cells (Kappa = 0.63) without over fitting. The ANN with greatest support, combining bathymetric values of the target grid cell with the slope of adjacent grid cells, was used to map inshore reef locations around the Southern Australian coastline (∼250 m resolution). Broadly, our results show that reefs are identifiable from coarse-scale bathymetry data of the seabed. We anticipate that our research technique will strengthen systematic conservation planning tools in many regions of the world, by enabling the identification of rocky substratum and mapping in localities that remain poorly surveyed due to logistics or monetary constraints.     

GIS-based modeling of the geographic distribution of Quercus emoryi Torr. (Fagaceae) in México and identification of significant environmental factors influencing the species’ distribution

The greatest concentration of oak species in the world is believed to be found in Mexico. These species are potentially useful for reforestation because of their capacity to adapt to diverse environments. Knowledge of their geographic distribution and of species–environment relations is essential for decision-making in the management and conservation of natural resources. The objectives of this study were to develop a model of the distribution of Quercus emoryi Torr. in Mexico, using geographic information systems and data layers of climatic and other variables, and to determine the variables that significantly influence the distribution of the species. The study consisted of the following steps: (A) selection of the target species from a botanical scientific collection, (B) characterization of the collecting sites using images with values or categories of the variables, (C) model building with the overlay of images that meet the habitat conditions determined from the characterization of sites, (D) model validation with independent data in order to determine the precision of the model, (E) model calibration through adjustment of the intervals of some variables, and (F) sensitivity analysis using precision and concordance non-parametric statistics applied to pairs of images. Results show that the intervals of the variables that best describe the species’ habitat are the following: altitude from 1650 to 2750 amsl, slope from 0 to 66°; average minimum temperature of January from −12 to −3 °C; mean temperature of June from 11 to 25 °C; mean annual precipitation from 218 to 1225 mm; soil units: lithosol, eutric cambisol, haplic phaeozem, chromic luvisol, rendzina, luvic xerosol, mollic planosol, pellic vertisol, eutric regosol; type of vegetation: oak forest, oak–pine forest, pine forest, pine–oak forest, juniperus forest, low open forest, natural grassland and chaparral. The resulting model of the geographic distribution of Quercus emoryi in Mexico had the following values for non-parametric statistics of precision and agreement: Kappa index of 0.613 and 0.788, overall accuracy of 0.806 and 0.894, sensitivity of 0.650 and 0.825, specificity of 0.963, positive predictive value of 0.945 and 0.957 and negative predictive value of 0.733 and 0.846. Results indicate that the variable average minimum temperature of January, with a maximum value of −3 °C, is an important factor in limiting the species’ distribution.     

Incorporating low-resolution historic species location data decreases performance of distribution models

Developing robust species distribution models is important as model outputs are increasingly being incorporated into conservation policy and management decisions. A largely overlooked component of model assessment and refinement is whether to include historic species occurrence data in distribution models to increase the data sample size. Data of different temporal provenance often differ in spatial accuracy and precision. We test the effect of inclusion of historic coarse-resolution occurrence data on distribution model outputs for 187 species of birds in Australian tropical savannas. Models using only recent (after 1990), fine-resolution data had significantly higher model performance scores measured with area under the receiver operating characteristic curve (AUC) than models incorporating both fine- and coarse-resolution data. The drop in AUC score is positively correlated with the total area predicted to be suitable for the species (R² = 0.163-0.187, depending on the environmental predictors in the model), as coarser data generally leads to greater predicted areas. The remaining unexplained variation is likely to be due to the covariate errors resulting from resolution mismatch between species records and environmental predictors. We conclude that decisions regarding data use in species distribution models must be conscious of the variation in predictions that mixed-scale datasets might cause.     

Simulating soil organic matter with CQESTR (v. 2.0): Model description and validation against long-term experiments across North America

Soil carbon (C) models are important tools for examining complex interactions between climate, crop and soil management practices, and to evaluate the long-term effects of management practices on C-storage potential in soils. CQESTR is a process-based carbon balance model that relates crop residue additions and crop and soil management to soil organic matter (SOM) accretion or loss. This model was developed for national use in U.S and calibrated initially in the Pacific Northwest. Our objectives were: (i) to revise the model, making it more applicable for wider geographic areas including potential international application, by modifying the thermal effect and incorporating soil texture and drainage effects, and (ii) to recalibrate and validate it for an extended range of soil properties and climate conditions. The current version of CQESTR (v. 2.0) is presented with the algorithms necessary to simulate SOM at field scale. Input data for SOM calculation include crop rotation, aboveground and belowground biomass additions, tillage, weather, and the nitrogen content of crop residues and any organic amendments. The model was validated with long-term data from across North America. Regression analysis of 306 pairs of predicted and measured SOM data under diverse climate, soil texture and drainage classes, and agronomic practices at 13 agricultural sites having a range of SOM (7.3–57.9 g SOM kg⁻¹), resulted in a linear relationship with an r² of 0.95 (P < 0.0001) and a 95% confidence interval of 4.3 g SOM kg⁻¹. Using the same data the version 1.0 of CQESTR had an r² of 0.71 with a 95% confidence interval of 5.5 g SOM kg⁻¹. The model can be used as a tool to predict and evaluate SOM changes from various management practices and offers the potential to estimate C accretion required for C credits.     

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.     

A general process-based mass-balance model for phosphorus/eutrophication as a tool to estimate historical reference values for key bioindicators,as exemplified using data for the Gulf of Riga

This work describes how a general, process-based mass-balance model (CoastMab) for phosphorus for coastal areas may be used as a tool to estimate realistic values of “natural” or preindustrial reference levels of key bioindicators in coastal science, including the Secchi depth, a standard measure of water clarity, the chlorophyll-a concentration, an operational measure of phytoplankton biomass and the concentration of cyanobacteria, a measure of the concentration of harmful algae. The CoastMab-model is an ecosystem model giving monthly predictions to achieve seasonal variations of basin-wide properties. The selected case-study area, the Gulf of Riga, is sensitive to nutrient loading because of its shallowness and low openness towards the Baltic Proper. The morphometry of any coastal area, as given by the size and form parameters, influences all internal processes, such as sedimentation, resuspension, diffusion in water and from sediments to water, biouptake and retention in biota, stratification, mixing and outflow. There has been no mass-balance modeling for nitrogen (N) in this work because empirical data (from the HELCOM database) clearly indicate that the monthly primary production in the Gulf of Riga is regulated by phosphorus (P) – the mean monthly total-N to total-P ratios are well over 7.2 (the Redfield-ratio) and generally higher than 15 for the data used in this study (from 1992 to 2005). At present anthropogenic loads, the average modeled monthly values for Secchi depth, chlorophyll (Chl), cyanobacteria (CB) and total-P (TP) are 3.2 m, 3.8 μg/l, 78 μg/l and 31.3 μg/l, respectively. If 50% of all anthropogenic sources to the Gulf of Riga via rivers, point sources and diffuse sources were to be removed, these values would be 3.6 m, 3.4 μg Chl/l, 63 μg CB/l and 29.1 μg TP/l. If 60% of the anthropogenic phosphorus fluxes to the Baltic Proper were to be omitted and as well as 75% of all direct anthropogenic sources to the Gulf of Riga, the values would be 4.6 m, 2.7 μg Chl/l, 45 μg CB/l and 25.4 μg TP/l. These values represent the “natural” reference levels and it is not realistic to expect that remedial measures would improve the conditions more than that. Using the CoastWeb-model, similar calculations can be made for any given coastal area and the data necessary for such calculations are discussed in this work.     

贺兰山高山草甸生物多样性和地上生物量的关系

在对贺兰山高山草甸进行群落调查的基础上,研究高山草甸生物多样性和地上生物量与环境因子之间的关系,进而分析生物多样性和地上生物量的关系.结果表明:(1)地上生物量主要与土层深度成正相关关系.(2)海拔高度与生物多样性成负相关关系,而其它影响物种丰富度或Shannon指数的环境因子仅在个别群落类型中起作用.(3)生物多样性与地上生物量主要呈单峰曲线关系.     

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.     

Formulation and testing of a novel river nitrification model

The nitrification process in many river water quality models has been approximated by a simple first order dependency on the water column ammonia concentration, while the benthic contribution has routinely been neglected. In this study a mathematical framework was developed for sediment bed nitrification based on mass transfer theory and Monod bacterial growth kinetics. The model describes ammonia transport across the boundary layer and consumption within the biofilm to quantify the overall nitrification flux. Model results suggest that nitrification is usually controlled by the boundary layer thickness, and we estimated a nitrification velocity range between 0.14 and 0.97 m d⁻¹, assuming typical boundary thicknesses of 0.1–1.0 mm. These ranges compared favorably with reported literature values, including our own measurements. The model was applied to several river systems of different depths where nitrification rates and river depths were available. Assuming that nitrification is exclusively a benthic process, the average velocity of all the rivers evaluated was 0.85 m d⁻¹ (r² = 0.72).     

Assessing spawning ground hydraulic suitability for Chinese sturgeon (Acipenser sinensis) from horizontal mean vorticity in Yangtze River

The Chinese sturgeon, Acipenser sinensis, is an anadromous protected species. Its migration pattern has been blocked since the construction of Gezhouba dam and the length of the natural spawning site reduced to less than 7 km. However, the fish eventually established an alternative spawning ground in the narrow downstream reach of Gezhouba dam. In this article, we applied Delft3D-Flow model to simulate the hydraulic suitability of the spawning ground downstream Gezhouba dam. Horizontal mean vorticity was used to assess the hydraulic environment of spawning ground. The flow field state was determined through model simulation and field-measured data used to validate the model. The computational method was improved by calculating absolute horizontal mean vorticity from estimates the literature. The final vorticity was determined from the simulation output and its distribution pattern retrieved. The horizontal mean vorticity range was 0.71–4.61 10⁻³ s⁻¹ for the entire spawning grounds, with egg mass field upper limit of 1.0 × 10⁻³ s⁻¹. Vorticity strength selection of Chinese sturgeon spawning can enhance our understanding of egg fertilization rate, hence the protection of fertilized eggs. Furthermore, the results of the study would add to existing scientific database for spawning ground hydraulic environment protection, especially in the ecological regulation drive of the Three Gorges Reservoir.     

Estimating soil carbon turnover using radiocarbon data: A case-study for European Russia

Turnover rates of soil carbon for 20 soil types typical for a 3.7 million km² area of European Russia were estimated based on ¹⁴C data. The rates are corrected for bomb radiocarbon which strongly affects the topsoil ¹⁴C balance. The approach is applied for carbon stored in the organic and mineral layers of the upper 1 m of the soil profile. The turnover rates of carbon in the upper 20 cm are relatively high for forest soils (0.16–0.78% year⁻¹), intermediate for tundra soils (0.25% year⁻¹), and low for grassland soils (0.02–0.08% year⁻¹) with the exception of southern Chernozems (0.32% year⁻¹). In the soil layer of 20–100 cm depth, the turnover rates were much lower for all soil types (0.01–0.06% year⁻¹) except for peat bog soils of the southern taiga (0.14% year⁻¹). Combined with a map of soil type distribution and a dataset of several hundred soil carbon profiles, the method provides annual fluxes for the slowest components of soil carbon assuming that the latter is in equilibrium with climate and vegetation cover. The estimated carbon flux from the soil is highest for forest soils (12–147 gC/(m² year)), intermediate for tundra soils (33 gC/(m² year)), and lowest for grassland soils (1–26 gC/(m² year)). The approach does not distinguish active and recalcitrant carbon fractions and this explains the low turnover rates in the top layer. Since changes in soil types will follow changes in climate and land cover, we suggest that pedogenesis is an important factor influencing the future dynamics of soil carbon fluxes. Up to now, both the effect of soil type changes and the clear evidence from ¹⁴C measurements that most soil organic carbon has a millennial time scale, are basically neglected in the global carbon cycle models used for projections of atmospheric CO₂ in 21st century and beyond.     

Convexity in projection matrices: Projection to a calibration problem

Convexity, as a fundamental property of sets and functions defined on convex sets, plays an important role in many mathematical and applied disciplines, including extremal and optimal-control problems. We prove the set of all feasible projection matrices in a general class of matrix models for stage-structured population dynamics to be convex and the dominant eigenvalue (λ₁) of any projection 2 × 2 matrix to be either a convex, or a concave function on a simplex of the matrix first-row entries (i.e., stage-specific reproduction rates). The latter is also conjectured for the general n × n case. Though looking far from practical needs of matrix population models, this mathematical result has appeared to be quite useful in solving a practical problem to calibrate the projection matrix, i.e., to estimate all the stage-specific vital rates, from empirical data. The data from monitoring of individual life histories of marked plants on permanent sample plots during successive years enable direct calculation of the stage-specific survival and ontogenetic transition rates, but the rates of reproduction do remain uncertain as far as the parent plants can hardly be determined for the (not yet marked!) recruitment.     

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.     

Surface modelling of global terrestrial ecosystems under three climate change scenarios

A high accuracy and speed method (HASM) of surface modelling is developed to find a solution for error problem and to improve computation speed. A digital elevation model (DEM) is established on spatial resolution of 13.5 km × 13.5 km. Regression formulations among temperature, elevation and latitude are simulated in terms of data from 2766 weather observation stations scattered over the world by using the 13.5 km × 13.5 km DEM as auxiliary data. Three climate scenarios of HadCM3 are refined from spatial resolution of 405 km × 270 km to 13.5 km × 13.5 km in terms of the regression formulations. HASM is employed to simulate surfaces of mean annual bio-temperature, mean annual precipitation and potential evapotranspiration ratio during the periods from 1961 to 1990 (T₁), from 2010 to 2039 (T₂), from 2040 to 2069 (T₃), and from 2070 to 2099 (T₄) on spatial resolution of 13.5 km × 13.5 km. Three scenarios of terrestrial ecosystems on global level are finally developed on the basis of the simulated climate surfaces. The scenarios show that all polar/nival, subpolar/alpine and cold ecosystem types would continuously shrink and all tropical types, except tropical rain forest in scenario A1Fi, would expand because of the climate warming. Especially at least 80% of moist tundra and 22% of nival area might disappear in period T₄ comparing with the ones in the period T₁. Tropical thorn woodland might increase by more than 97%. Subpolar/alpine moist tundra would be the most sensitive ecosystem type because its area would have the rapidest decreasing rate and its mean center would shift the longest distance towards west. Subpolar/alpine moist tundra might be able to serve as an indicator of climatic change. In general, climate change would lead to a continuous reduction of ecological diversity.