过渡带中岳桦种群生态特征的研究

在长白山北坡halt1525m-1775m的岳桦云冷杉森林过渡带和1990m-2170m的岳桦苔原林线过渡带中，沿海拔高度，采用梯度取样方法，进行岳桦种群生态学研究，分析了岳桦种群存活曲线、种群分布格局和种群竞争指数。同时应用分形分析研究了岳桦种群的分形特征。结果表明：依据岳桦种群的存活曲线，在岳桦云冷杉过渡带中岳桦种群处于由衰退型向稳定型发展的阶段，而在岳桦苔原过渡带中则由稳定型向增长型过渡，从种群分布格局指数分析，岳桦在两个过渡带中的分布格局类型都是杂乱无章的，而且都偏离随机分布，处于由一个生态系统向另一个生态系统转换的相变区，从竞争指数和盒维数来看，在岳桦云冷杉过渡带中，这两项指标都是由小到大迅速递增，而在岳桦苔原过渡 情况刚好相反。     

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.     

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.     

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.     

A geospatial model of forest dynamics with controlled trend surface

This paper proposes a method of controlled trend surface to simultaneously account for large-scale spatial trends and non-spatial local effects. With this method, a geospatial model of forest dynamics was developed for the Alaska boreal forest from 446 constantly monitored permanent sample plots. The geospatial component of this model represented large-scale spatial trends in recruitment, diameter growth, and mortality. The model was tested on two sets of validation plots which represented temporal and spatial extensions of the current sample coverage. The results suggest that the controlled trend surface model was generally more accurate than both the non-spatial and conventional trend surface models. With this model, we mapped the forest dynamics of the entire Alaska boreal region by aggregating predicted stand states across the region. It was predicted that under current conditions of climate and natural disturbances, most of the Alaska boreal forest region may undergo a major shift from deciduous-dominant to conifer-dominant, with an average increase of 0.33 m² ha year⁻¹ in basal area over the Twenty-First Century.     

A 700-year paleoecological record of boreal ecosystem responses to climatic variation from Alaska

Recent observations and model simulations have highlighted the sensitivity of the forest-tundra ecotone to climatic forcing. In contrast, paleoecological studies have not provided evidence of tree-line fluctuations in response to Holocene climatic changes in Alaska, suggesting that the forest-tundra boundary in certain areas may be relatively stable at multicentennial to millennial time scales. We conducted a multiproxy study of sediment cores from an Alaskan lake near the altitudinal limits of key boreal-forest species. Paleoecological data were compared with independent climatic reconstructions to assess ecosystem responses of the forest tundra boundary to Little Ice Age (LIA) climatic fluctuations. Pollen, diatom, charcoal, macrofossil, and magnetic analyses provide the first continuous record of vegetation fire-climate interactions at decadal to centennial time scales during the past 700 years from southern Alaska. Boreal-forest diebacks characterized by declines of Picea mariana, P. glauca, and tree Betula occurred during the LIA (AD 1500-1800), whereas shrubs (Alnus viridis, Betula glandulosa/nana) and herbaceous taxa (Epilobium, Aconitum) expanded. Marked increases in charcoal abundance and changes in magnetic properties suggest increases in fire importance and soil erosion during the same period. In addition, the conspicuous reduction or disappearance of certain aquatic (e.g., Isoetes, Nuphar, Pediastrum) and wetland (Sphagnum) plants and major shifts in diatom assemblages suggest pronounced lake-level fluctuations and rapid ecosystem reorganization in response to LIA climatic deterioration. Our results imply that temperature shifts of 1-2 degrees C, when accompanied by major changes in moisture balance, can greatly alter high-altitudinal terrestrial, wetland, and aquatic ecosystems, including conversion between boreal-forest tree line and tundra. The climatic and ecosystem variations in our study area appear to be coherent with changes in solar irradiance, suggesting that changes in solar activity contributed to the environmental instability of the past 700 years.     

粤北亚热带山地森林土壤有机碳沿海拔梯度的变化

选择位于南岭国家级自然保护区的广东第一峰——石坑崆,从海拔300 m起到山顶部1 900 m范围,每隔100 m高程设置1条10 m×120 m的样带,共17条样带,研究土壤有机碳含量沿海拔梯度的变化规律及其与植被类型和凋落物层厚度的关系。结果表明,0～20 cm和>20～40 cm土层有机碳含量均随海拔梯度变化呈极显著差异(P<0.001),并随林分类型不同而呈高度显著差异(P<0.01),土壤有机碳含量总体上呈随海拔上升而升高的变化趋势。凋落物层厚度仅对0～20 cm土层有机碳含量有显著影响(P<0.05)。可见,海拔梯度变化是影响土壤有机碳含量的综合和主导因素,而最表层土壤有机碳含量还易受林分因子的影响,这些因子反映了土地利用变化及自然保护历史。该研究结果提示,减少人类活动的干扰和保持林地适合的凋落物层厚度有助于增加森林生态系统的碳储量。     

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.     

Experimental evidence for herbivore limitation of the treeline

The treeline ecotone divides forest from open alpine or arctic vegetation states. Treelines are generally perceived to be temperature limited. The role of herbivores in limiting the treeline is more controversial, as experimental evidence from relevant large scales is lacking. Here we quantify the impact of different experimentally controlled herbivore densities on the recruitment and survival of birch Betula pubescens tortuosa along an altitudinal gradient in the mountains of southern Norway. After eight years of summer grazing in large-scale enclosures at densities of 0, 25, and 80 sheep/km2, birch recruited within the whole altitudinal range of ungrazed enclosures, but recruitment was rarer in enclosures with low-density sheep and was largely limited to within the treeline in enclosures with high-density sheep. In contrast, the distribution of saplings (birch older than the experiment) did not differ between grazing treatments, suggesting that grazing sheep primarily limit the establishment of new tree recruits rather than decrease the survival of existing individuals. This study provides direct experimental evidence that herbivores can limit the treeline below its potential at the landscape scale and even at low herbivore densities in this climatic zone. Land use changes should thus be considered in addition to climatic changes as potential drivers of ecotone shifts.     

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.     

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.     

Community-scale analysis of the farmland abandonment occurrence process in the mountain region of Ladakh,India

Farmland abandonment has been a significant issue in mountain regions. Thus, many studies have attempted to uncover its driving factors. Most existing studies take a large- or meso-scale view of mountain regions, and there are a few studies that provide community-level analysis. Therefore, it is unclear how ecological and sociopolitical factors are actually combined to influence farmland abandonment at a community level. Thus, we conducted a field survey in the mountain village of Ladakh, northern India, beginning in 2009 and tried to understand the farmland abandonment occurrence process from the viewpoint of traditional altitudinal models of land use. Using data collected during the survey, a relationship between altitude and farmland abandonment was established. However, this relationship was complex, combining traditional forms of land use system in the village. The altitudinal model of land use was very common in mountain regions; thus, discussions of current land use issues should be based on an understanding of such traditional altitudinal systems.     

Bats aloft: variability in echolocation call structure at high altitudes

Bats alter their echolocation in response to changes in ecological and behavioral conditions, but little is known about how they adjust call structure in response to changes in altitude. We examined altitudinal variation in the echolocation of Brazilian free-tailed bats, Tadarida brasiliensis, a species known to fly to altitudes of 3,000 m above the ground. From 50.2 h of recordings, we analyzed 113 high-quality echolocation call sequences recorded from 0 to 862 m above ground level. Bats flying near the ground used shorter, higher-frequency, broader-bandwidth calls compared to bats at higher altitudes, an effect likely due to the greater levels of echo-producing clutter (i.e., vegetation, buildings) found near the ground. When ground-level recordings are excluded, bats continue to shift towards the use of longer-duration, lower-frequency, narrower-bandwidth calls with increasing altitude. We propose that the observed high-altitude changes in call structure are a response to changing acoustic attenuation rates and/or decreasing insect densities at higher altitudes.     

An innovative computer design for modeling forest landscape change in very large spatial extents with fine resolutions

Although forest landscape models (FLMs) have benefited greatly from ongoing advances of computer technology and software engineering, computing capacity remains a bottleneck in the design and development of FLMs. Computer memory overhead and run time efficiency are primary limiting factors when applying forest landscape models to simulate large landscapes with fine spatial resolutions and great vegetation detail. We introduce LANDIS PRO 6.0, a landscape model that simulates forest succession and disturbances on a wide range of spatial and temporal scales. LANDIS PRO 6.0 improves on existing forest landscape models with two new data structures and algorithms (hash table and run-length compression). The innovative computer design enables LANDIS PRO 6.0 to simulate very large (>10⁸ ha) landscapes with a 30-m spatial resolution, which to our knowledge no other raster forest landscape models can do. We demonstrate model behavior and performance through application to five nested forest landscapes with varying sizes (from 1 million to 100 million 0.09-ha cells) in the southern Missouri Ozarks. The simulation results showed significant and variable effects of changing spatial extent on simulated forest succession patterns. Results highlighted the utility of a model like LANDIS PRO 6.0 that is capable of efficiently simulating large landscapes and scaling up forest landscape processes to a common regional scale of analysis. The programming methodology presented here may significantly advance the development of next generation of forest landscape models.     

Impacts of daily weather variability on simulations of the Canadian boreal forest

This study examines the importance of climate variability when simulating forest succession using a process-based model of stand development. The FORSKA-2V forest gap model, originally developed for forcing with monthly mean climate data, was modified to accept daily weather data. The model's performance was compared using different temporal resolutions of forcing along a bioclimatic transect crossing the boreal region of central Canada, including the aspen-parkland and forest-tundra ecotones. Forcing the model with daily weather data improved the simulation of key attributes of present-day forest along the transect, particularly at the ecotones, compared to forcing with monthly data or long term averages. The results support the hypothesis that climatic variability at daily time-scales is an important determinant of present-day boreal forest composition and productivity. To simulate boreal forest response to climatic change it will be necessary to create climatic scenarios that include plausible projections of future daily scale variability.     

Foliage profiles of individual trees determine competition,self-thinning,biomass and NPP of a Cryptomeria japonica forest stand: A simulation study based on a stand-scale process-based forest model

A simulation study was carried out to investigate simultaneously the effects of eco-physiological parameters on competitive asymmetry, self-thinning, stand biomass and NPP in a temperate forest using an atmosphere–vegetation dynamics interactive model (MINoSGI). In this study, we selected three eco-physiological relevant parameters as foliage profiles (i.e. vertical distribution of leaf area density) of individual trees (distribution pattern is described by the parameter η), biomass allocation pattern in individual tree growth (χ) and the maximum carboxylation velocity (V_max). The position of the maximal leaf area density shifts upward in the canopy with increasing η. For scenarios with η < 4 (foliage concentrated in the lowest canopy layer) or η > 12 (foliage concentrated in the uppermost canopy layer), a low degree of competitive asymmetry was produced. These scenarios resulted in the survival of subordinate trees due to a brighter lower canopy environment when η < 4 or the generation of spatially separated foliage profiles between dominant and subordinate trees when η > 12. In contrast, competition between trees was most asymmetric when 4 ≤ η ≤ 12 (vertically widespread foliage profile in the canopy), especially when η = 8. In such cases, vertically widespread foliage of dominant trees lowered the opportunity of light acquisition for subordinate trees and reduced their carbon gain. The resulting reduction in carbon gain of subordinate trees yielded a higher degree of competitive asymmetry and ultimately higher mortality of subordinate trees. It was also shown that 4 ≤ η ≤ 12 generated higher self-thinning speed, smaller accumulated NPP, litter-fall and potential stand biomass as compared with the scenarios with η < 4 or η > 12. In contrast, our simulation revealed small effects of χ or V_max on the above-mentioned variables as compared with those of η. In particular, it is notable that greater V_max would not produce greater potential stand biomass and accumulated NPP although it has been thought that physiological parameters relevant to photosynthesis such as V_max influence dynamic changes in forest stand biomass and NPP (e.g. the greater the V_max, the greater the NPP). Overall, it is suggested that foliage profiles rather than biomass allocation or maximum carboxylation velocity greatly govern forest dynamics, stand biomass, NPP and litter-fall.     

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.     

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.     

Soil-geomorphic heterogeneity governs patchy vegetation dynamics at an arid ecotone

Soil properties are well known to affect vegetation, but the role of soil heterogeneity in the patterning of vegetation dynamics is poorly documented. We asked whether the location of an ecotone separating grass-dominated and sparsely vegetated areas reflected only historical variation in degradation or was related to variation in inherent soil properties. We then asked whether changes in the cover and spatial organization of vegetated and bare patches assessed using repeat aerial photography reflected self-organizing dynamics unrelated to soil variation or the stable patterning of soil variation. We found that the present-day ecotone was related to a shift from more weakly to more strongly developed soils. Parts of the ecotone were stable over a 60-year period, but shifts between bare and vegetated states, as well as persistently vegetated and bare states, occurred largely in small (<40 m2) patches throughout the study area. The probability that patches were presently vegetated or bare, as well as the probability that vegetation persisted and/or established over the 60-year period, was negatively related to surface calcium carbonate and positively related to subsurface clay content. Thus, only a fraction of the landscape was susceptible to vegetation change, and the sparsely vegetated area probably featured a higher frequency of susceptible soil patches. Patch dynamics and self-organizing processes can be constrained by subtle (and often unrecognized) soil heterogeneity.     

The effects of forest harvest intensity in combination with wind disturbance on carbon dynamics in Lake States Mesic Forests

Total forest carbon (C) storage is determined by succession, disturbances, climate, and the edaphic properties of a site or region. Forest harvesting substantially affects C dynamics; these effects may be amplified if forest harvesting is intensified to provide biofuel feedstock. We tested the effects of harvest intensity on landscape C using a simulation modeling approach that included C dynamics, multiple disturbances, and successional changes in composition. We developed a new extension for the LANDIS-II forest landscape disturbance and succession model that incorporates belowground soil C dynamics derived from the CENTURY soil model. The extension was parameterized and calibrated using data from an experimental forest in northeastern Wisconsin, USA. We simulated a 9800 ha forested landscape over 400 years with wind disturbance combined with no harvesting, harvesting with residual slash left on site (‘standard harvest’), and whole-tree harvesting. We also simulated landscapes without wind disturbance and without eastern hemlock (Tsuga canadensis) to examine the effects of detrital quantity and quality on C dynamics. We estimated changes in live C, detrital C, soil organic C, total C, and forest composition. Overall, the simulations without harvesting had substantially greater total C and continued to sequester C. Standard harvest simulations had more C than the whole tree harvest simulations. Under both harvest regimes, C accrual was not evident after 150 years. Without hemlock, SOC was reduced due to a decline in detritus and a shift in detrital chemistry. In conclusion, if the intensity of harvesting increases we can expect a corresponding reduction in potential C storage. Compositional changes due to historic circumstances (loss of hemlock) may also affect forest C although to a lesser degree than harvesting. The modeling approach presented enabled us to consider multiple, interacting drivers of landscape change and the subsequent changes in forest C.