Retrieving leaf area index for coniferous forest in Xingguo County, China with Landsat ETM+ images

Spatial distributions of the leaf area index (LAI) needed for carbon cycle modeling in Xingguo County, China were estimated based on correlations between the field-measurements and vegetation indices (VIs). After making geometric and atmospheric corrections to two Landsat ETM+ images, one in January 2000 and the other in May 2003, three VIs (SR, NDVI, and RSR) were derived, and their separate correlations with ground LAI measurements were established. The correlation with RSR was the highest among the three VIs. The retrieved LAI values for January 2000 were lower than those for May 2003 because of a small seasonal variation in the coniferous forests (predominantly masson pine) and the decrease in the understorey vegetation during winter.     

Spatial scaling between leaf area index maps of different resolutions

We developed algorithms for spatial scaling of leaf area index (LAI) using sub-pixel information. The study area is located near Liping County, Guizhou Province, in China. Methods for LAI spatial scaling were investigated on LAI images with 960 m resolution derived in two ways. LAI from distributed calculation (LAID) was derived using Landsat ETM+ data (30 m), and LAI from lumped calculation (LAIL) was obtained from the coarse (960 m) resolution data derived through resampling the ETM+ data. We found that lumped calculations can be considerably biased compared to the distributed (ETM+) case, suggesting that global and regional LAI maps can be biased if surface heterogeneity within the mapping resolution is ignored. Based on these results, we developed algorithms for removing the biases in lumped LAI maps using sub-pixel land cover-type information, and applied these to correct one coarse resolution LAI product which greatly improved its accuracy.     

Tests of soil organic carbon density modeled by InTEC in China's forest ecosystems

The integrated terrestrial ecosystem C-budget model (InTEC) developed by Chen and co-workers has been used successfully to predict carbon dynamics of forests in Canada. It was tested here for forest soil organic carbon (SOC) density of China's northern temperate zone and southern subtropical zone. The results show that the simulated SOC density is highly correlated and in broad agreement with observations in Liping and in Changbaishan, representing the southern subtropical zone and the northern temperate zone in China, respectively. SOC density ranged from 2.2 to 11.2 kg/m(2) in Liping and from 3.4 to 14.8 kg/m(2) in Changbaishan. The correlation coefficients (r(2)) are 0.63 (N=16) and 0.76 (N=14) between the simulated and measured data in Liping and Changbaishan, respectively. The SOC densities under different vegetation types in Liping decrease in the order of mixed forest, broadleaf forest, Chinese fir, couch grass, and Chinese redpine, and in Changbaishan in the order of mixed forest, silver fir, larch forest, and birch forest.     

Relationships between net primary productivity and stand age for several forest types and their influence on China's carbon balance

Affected by natural and anthropogenic disturbances such as forest fires, insect-induced mortality and harvesting, forest stand age plays an important role in determining the distribution of carbon pools and fluxes in a variety of forest ecosystems. An improved understanding of the relationship between net primary productivity (NPP) and stand age (i.e., age-related increase and decline in forest productivity) is essential for the simulation and prediction of the global carbon cycle at annual, decadal, centurial, or even longer temporal scales. In this paper, we developed functions describing the relationship between national mean NPP and stand age using stand age information derived from forest inventory data and NPP simulated by the BEPS (Boreal Ecosystem Productivity Simulator) model in 2001. Due to differences in ecobiophysical characteristics of different forest types, NPP-age equations were developed for five typical forest ecosystems in China (deciduous needleleaf forest (DNF), evergreen needleleaf forest in tropic and subtropical zones (ENF-S), deciduous broadleaf forest (DBF), evergreen broadleaf forest (EBF), and mixed broadleaf forest (MBF)). For DNF, ENF-S, EBF, and MBF, changes in NPP with age were well fitted with a common non-linear function, with R(2) values equal to 0.90, 0.75, 0.66, and 0.67, respectively. In contrast, a second order polynomial was best suitable for simulating the change of NPP for DBF, with an R(2) value of 0.79. The timing and magnitude of the maximum NPP varied with forest types. DNF, EBF, and MBF reached the peak NPP at the age of 54, 40, and 32 years, respectively, while the NPP of ENF-S maximizes at the age of 13 years. The highest NPP of DBF appeared at 122 years. NPP was generally lower in older stands with the exception of DBF, and this particular finding runs counter to the paradigm of age-related decline in forest growth. Evaluation based on measurements of NPP and stand age at the plot-level demonstrates the reliability and applicability of the fitted NPP-age relationships. These relationships were used to replace the normalized NPP-age relationship used in the original InTEC (Integrated Terrestrial Ecosystem Carbon) model, to improve the accuracy of estimated carbon balance for China's forest ecosystems. With the revised NPP-age relationship, the InTEC model simulated a larger carbon source from 1950-1980 and a larger carbon sink from 1985-2001 for China's forests than the original InTEC model did because of the modification to the age-related carbon dynamics in forests. This finding confirms the importance of considering the dynamics of NPP related to forest age in estimating regional and global terrestrial carbon budgets.     

利用数码相片对冬小麦生物量的试算

数码相片属于地面遥感资料的一种,被广泛应用于实地调查中.然而,目前利用数码相片对植被生物量进行推算的研究很少.利用数码相片提取小麦生长期的小麦覆盖度,结合地面实测叶面积指数(LAI)、归一化植被指数(NDVI)和小麦生物量,分析小麦覆盖度与LAI,NDVI,生物量与小麦覆盖度、LAI,NDVI各自的相关系数,建立各自合适的回归模型推算冬小麦的生物量.结果表明,在小麦生长期,小麦覆盖度与LAI,NDVI的相关系数R2分别达到0.765,0.896,生物量与小麦覆盖度、LAI,NDVI的相关系数R2分别为0.774,0.876,0.712,生物量与其小麦覆盖度、LAI,NDVI之间具有很高的相关性,多元回归分析效果更好,系数R2达0.891.结果说明,在不破坏植被生长状态的情况下,基于数码相片可有效地推算小麦的生物量,这为大面积实地测量和遥感监测作物长势及其生物量估算提供了理论参考依据.     

Carbon sinks and sources in China's forests during 1901-2001

This paper reports the annual carbon (C) balance of China's forests during 1901-2001 estimated using the Integrated Terrestrial Ecosystem C-budget model (InTEC). Annual carbon source and sink distributions are simulated for the same period using various spatial datasets including land cover and leaf area index (LAI) obtained from remote sensing, soil texture, climate, forest age, and nitrogen deposition. During 1901-1949, China's forests were a source of 21.0+/-7.8 Tg C yr(-1) due to disturbances (human activities). Its size increased to 122.3+/-25.3 Tg C yr(-1) during 1950-1987 due to intensified human activities in the late 1950s, early 1960s, 1970s and early 1980s. The forests became large sinks of 176.7+/-44.8 Tg C yr(-1) during 1988-2001, owing to large-scale plantation and forest regrowth in previously disturbed areas as well as growth stimulation by nondisturbance factors such as climatic warming, atmospheric CO(2) fertilization, and N deposition. From 1901 to 2001, China's forests were a small carbon source of 3.32 Pg C, about 32.9+/-22.3 Tg C yr(-1). The overall C balance in biomass from InTEC generally agrees with previous results derived from forest inventories of China's forests. InTEC results also include C stock variation in soils and are therefore more comprehensive than previous results. The uncertainty in InTEC results is still large, but it can be reduced if a detailed forest age map becomes available.     

Integration of Satellite Imagery and Forest Inventory in Mapping Dominant and Associated Species at a Regional Scale

To achieve the overall objective of restoring natural environment and sustainable resource usability, each forest management practice effect needs to be predicted using a simulation model. Previous simulation efforts were typically confined to public land. Comprehensive forest management practices entail incorporating interactions between public and private land. To make inclusion of private land into management planning feasible at the regional scale, this study uses a new method of combining Forest Inventory and Analysis (FIA) data with remotely sensed forest group data to retrieve detailed species composition and age information for the Missouri Ozark Highlands. Remote sensed forest group and land form data inferred from topography were integrated to produce distinct combinations (ecotypes). Forest types and size classes were assigned to ecotypes based on their proportions in the FIA data. Then tree species and tree age determined from FIA subplots stratified by forest type and size class were assigned to pixels for the entire study area. The resulting species composition map can improve simulation model performance in that it has spatially explicit and continuous information of dominant and associated species, and tree ages that are unavailable from either satellite imagery or forest inventory data. In addition, the resulting species map revealed that public land and private land in Ozark Highlands differ in species composition and stand size. Shortleaf pine is a co-dominant species in public land, whereas it becomes a minor species in private land. Public forest is older than private forest. Both public and private forests have deviated from historical forest condition in terms of species composition. Based on possible reasons causing the deviation discussed in this study, corresponding management avenues that can assist in restoring natural environment were recommended.     

Variation in Vegetation Structure and Soil Properties, and the Relation Between Understory Plants and Environmental Variables Under Different Phyllostachys pubescens Forests in Southeastern China

Biodiversity maintenance and soil improvement are key sustainable forestry objectives. Research on the effects of bamboo forest management on plant diversity and soil properties are therefore necessary in bamboo-growing regions, such as southeastern China’s Shunchang County, that have not been studied from this perspective. We analyzed the effects of different Phyllostachys pubescens proportions in managed forests on vegetation structure and soil properties using pure Cunninghamia lanceolata forests as a contrast, and analyzed the relation between understory plants and environmental variables (i.e., topography, stand and soil characteristics) by canonical correspondence analysis (CCA). The forest with 80% P. pubescens and 20% hardwoods (such as Phoebe bournei, Jatropha curcas, Schima superba) maintained the highest plant diversity and best soil properties, with significantly higher plant diversity than the C. lanceolata forest, and better soil physicochemical and biological properties. The distribution of understory plants is highly related to environmental factors. Silvicultural disturbance strongly influenced the ability of different bamboo forests to maintain biodiversity and soil quality under extensive management, and the forest responses to management were consistent with the intermediate-disturbance hypothesis (i.e., diversity and soil properties were best at intermediate disturbance levels). Our results suggest that biodiversity maintenance and soil improvement are important management goals for sustainable bamboo management. To achieve those objectives, managers should balance the inputs and outputs of nutrients and protect understory plants by using appropriate fertilizer (e.g., organic fertilizer), adjusting stand structure, modifying utilization model and the harvest time, and controlling the intensity of culms and shoots harvests.     

基于遥感数据的生态脆弱区森林景观分类及格局研究——以香格里拉为例

以Landsat-7 ETM＋为数据源，在遥感（RS）和地理信息系统（GIS）技术的支持下，对云南省香格里拉县北部进行森林景观分类及其景观结构分析。结果表明，研究区以针叶林为主要森林景观，阔叶混交林的形状结构较为松散，景观破碎化程度最严重，人为干扰强烈。因此，需加强森林保护，减少人为干扰，防止景观破碎化程度的进一步加深。     

The Vegetation Coverage Dynamic Coupling with Climatic Factors in Northeast China Transect

Based on SPOT-VGT images and meteorological data, this paper applied an integrated method to investigate the vegetation dynamic and its response to climate factors during 1998-2008 in Northeast China Transect, one of 15 ecological transects listed in the International Geosphere-Biosphere Programme. The main findings are as follows: (1) The NDVI time series presented nonlinear patterns that vary with timescales. The series fluctuated greatly at the smallest timescale (20?days), showing no salient trend, whereas a trend manifested itself more and more with the increase of time scale and finally stabilized at the 320-day scale. Little difference was found between vegetation types about the NDVI periodicity, as they occurred on either a 280-day or a 290-day cycle. (2) NDVI exhibited a significant correlation with temperature, precipitation, and sunshine hours. Overall, the correlation between NDVI and temperature was the highest, followed by precipitation, sunshine hours, and relative humidity. For different vegetation types, the correlations between NDVI and climate variables diversified, increasing from desert steppe to typical steppe, meadow steppe, and forest. (3) The periodicity of temperature and precipitation occurred in either a 280-day or 290-day cycle, which was approximately coincident with that of NDVI. This further supported the significant relationship between NDVI and these two climate factors. (4) At all the time scales under examination, NDVI and temperature and precipitation are significantly, positively correlated, especially at the 160-day scale, which can be regarded as the most suitable time scale for investigating the responses of vegetation dynamics to climate factors at most stations.     

The Importance of Temporal and Spatial Vegetation Structure Information in Biotope Mapping Schemes: A Case Study in Helsingborg,Sweden

Temporal and spatial vegetation structure has impact on biodiversity qualities. Yet, current schemes of biotope mapping do only to a limited extend incorporate these factors in the mapping. The purpose of this study is to evaluate the application of a modified biotope mapping scheme that includes temporal and spatial vegetation structure. A refined scheme was developed based on a biotope classification, and applied to a green structure system in Helsingborg city in southern Sweden. It includes four parameters of vegetation structure: continuity of forest cover, age of dominant trees, horizontal structure, and vertical structure. The major green structure sites were determined by interpretation of panchromatic aerial photographs assisted with a field survey. A set of biotope maps was constructed on the basis of each level of modified classification. An evaluation of the scheme included two aspects in particular: comparison of species richness between long-continuity and short-continuity forests based on identification of woodland continuity using ancient woodland indicators (AWI) species and related historical documents, and spatial distribution of animals in the green space in relation to vegetation structure. The results indicate that (1) the relationship between forest continuity: according to verification of historical documents, the richness of AWI species was higher in long-continuity forests; Simpson’s diversity was significantly different between long- and short-continuity forests; the total species richness and Shannon’s diversity were much higher in long-continuity forests shown a very significant difference. (2) The spatial vegetation structure and age of stands influence the richness and abundance of the avian fauna and rabbits, and distance to the nearest tree and shrub was a strong determinant of presence for these animal groups. It is concluded that continuity of forest cover, age of dominant trees, horizontal and vertical structures of vegetation should now be included in urban biotope classifications.     

Soil organic carbon decomposition and carbon pools in temperate and sub-tropical forests in China

Decomposition of soil organic carbon (SOC) is a critical component of the global carbon cycle, and accurate estimates of SOC decomposition are important for forest carbon modeling and ultimately for decision making relative to carbon sequestration and mitigation of global climate change. We determined the major pools of SOC in four sites representing major forest types in China: temperate forests at Changbai Mountain (CBM) and Qilian Mountain (QLM), and sub-tropical forests at Yujiang (YJ) and Liping (LP) counties. A 90-day laboratory incubation was conducted to measure CO(2) evolution from forest soils from each site, and data from the incubation study were fitted to a three-pool first-order model that separated mineralizable soil organic carbon into active (C(a)), slow (C(s)) and resistant (C(r)) carbon pools. Results indicate that: (1) the rate of SOC decomposition in the sub-tropical zone was faster than that in the temperature zone, (2) The C(a) pool comprised approximately 1-3% of SOC with an average mean residence time (MRT) of 219 days. The C(s) pool comprised approximately 25-65% with an average MRT of 78 yr. The C(r) pool accounted for approximately 35-80% of SOC, (3) The YJ site in the sub-tropical zone had the greatest C(a) pool and the lowest MRT, while the QLM in the temperature zone had the greatest MRT for both the C(a) and C(s) pools. The results suggest a higher capacity for long-term C sequestration as SOC in temperature forests than in sub-tropical forests.     

Relationships Between Major Ownerships, Forest Aboveground Biomass Distributions, and Landscape Dynamics in the New England Region of USA

This study utilizes remote sensing derived forest aboveground biomass (AGB) estimates and ownership information obtained from the Protected Areas Database (PAD), combining landscape analyses and GIS techniques to demonstrate how different ownerships (public, regulated private, and other private) relate to the spatial distribution of AGB in New England states of the USA. “Regulated private” lands were dominated by lands in Maine covered by a Land Use Regulatory Commission. The AGB means between all pairs of the identified ownership categories were significantly different (P < 0.05). Mean AGB observed in public lands (156 Mg/ha) was 43% higher than that in regulated private lands (109 Mg/ha), or 30% higher than that of private lands as a whole. Seventy-seven percent of the regional forests (or about 9,300 km²) with AGB >200 Mg/ha were located outside the area designated in the PAD and concentrated in western MA, southern VT, southwestern NH, and northwestern CT. While relatively unfragmented and high-AGB forests (>200 Mg/ha) accounted for about 8% of total forested land, they were unevenly proportioned among the three major ownership groups across the region: 19.6% of the public land, 0.8% of the regulated private land, and 11.0% of the other private land. Mean disturbance rates (in absolute value) between 1992 and 2001 were 16, 66, and 19 percent, respectively, on public, regulated private, and other private land. This indicates that management practices from different ownerships have a strong impact on dynamic changes of landscape structures and AGB distributions. Our results may provide insight information for policy makers on issues regarding forest carbon management, conservation biology, and biodiversity studies at regional level.     

Observation and simulation of net primary productivity in Qilian Mountain, western China

We modeled net primary productivity (NPP) at high spatial resolution using an advanced spaceborne thermal emission and reflection radiometer (ASTER) image of a Qilian Mountain study area using the boreal ecosystem productivity simulator (BEPS). Two key driving variables of the model, leaf area index (LAI) and land cover type, were derived from ASTER and moderate resolution imaging spectroradiometer (MODIS) data. Other spatially explicit inputs included daily meteorological data (radiation, precipitation, temperature, humidity), available soil water holding capacity (AWC), and forest biomass. NPP was estimated for coniferous forests and other land cover types in the study area. The result showed that NPP of coniferous forests in the study area was about 4.4 tCha(-1)y(-1). The correlation coefficient between the modeled NPP and ground measurements was 0.84, with a mean relative error of about 13.9%.     

Effect of tree stand and site variables on alpha diversity of ground vegetation in the forests of Northern Karelia

Generalised linear models (GLMs) were used to estimate the alpha diversity of the ground vegetation of boreal managed forests. Alpha diversity was described using Simpson's reciprocal, which was employed as a response variable in modelling, the explanatory variables being the tree and site characteristics. The data were from Northern Karelia, Eastern Finland. Several combinations of error and link functions were used to examine the patterns in the data, the best results being obtained with a model containing the normal error and identity link function, and the second power of the explanatory variable and one containing the gamma error and reciprocal link function. Site type, the presence of pubescent birch and the basal area of the stand were the significant variables in these models.It was concluded that GLMs are suitable tools for estimating biodiversity within the scope of the present data. In addition, it was found that the highest ground vegetation diversity exists within fertile stands with a mixture of birch and a low stand basal area. Models of this kind would be useful in forest management planning in order to study the effect of management on vegetational diversity using the values of different tree stand and site variables before and after treatment.     

LAI inversion algorithm based on directional reflectance kernels

Leaf area index (LAI) is an important ecological and environmental parameter. A new LAI algorithm is developed using the principles of ground LAI measurements based on canopy gap fraction. First, the relationship between LAI and gap fraction at various zenith angles is derived from the definition of LAI. Then, the directional gap fraction is acquired from a remote sensing bidirectional reflectance distribution function (BRDF) product. This acquisition is obtained by using a kernel driven model and a large-scale directional gap fraction algorithm. The algorithm has been applied to estimate a LAI distribution in China in mid-July 2002. The ground data acquired from two field experiments in Changbai Mountain and Qilian Mountain were used to validate the algorithm. To resolve the scale discrepancy between high resolution ground observations and low resolution remote sensing data, two TM images with a resolution approaching the size of ground plots were used to relate the coarse resolution LAI map to ground measurements. First, an empirical relationship between the measured LAI and a vegetation index was established. Next, a high resolution LAI map was generated using the relationship. The LAI value of a low resolution pixel was calculated from the area-weighted sum of high resolution LAIs composing the low resolution pixel. The results of this comparison showed that the inversion algorithm has an accuracy of 82%. Factors that may influence the accuracy are also discussed in this paper.     

CLIMATE CHANGE IMPACTS ON THE HYDROLOGY AND PRODUCTIVITY OF A PINE PLANTATION1

ABSTRACT: There are increasing concerns in the forestry community about global climate change and variability associated with elevated atmospheric CO₂. Changes in precipitation and increases in air temperature could impose additional stress on forests during the next century. For a study site in Carteret County, North Carolina, the General Circulation Model, HADCM2, predicts that by the year 2099, maximum air temperature will increase 1.6 to 1.9°C, minimum temperature will increase 2.5 to 2.8°C, and precipitation will increase 0 to 10 percent compared to the mid‐1990s. These changes vary from season to season. We utilized a forest ecosystem process model, PnET‐II, for studying the potential effects of climate change on drainage outflow, evapotranspiration, leaf area index (LAI) and forest Net Primary Productivity (NPP). This model was first validated with long term drainage and LAI data collected at a 25‐ha mature loblolly pine (Pinus taeda L.) experimental watershed located in the North Carolina lower coastal plain. The site is flat with poorly drained soils and high groundwater table. Therefore, a high field capacity of 20 cm was used in the simulation to account for the topographic effects. This modeling study suggested that future climate change would cause a significant increase of drainage (6 percent) and forest productivity (2.5 percent). Future studies should consider the biological feedback (i.e., stomata conductance and water use efficiency) to air temperature change.     

Water‐Yield Reduction After Afforestation and Related Processes in the Semiarid Liupan Mountains,Northwest China1

Abstract: The increase of coverage of forest/vegetation is imperative to improve the environment in dry‐land areas of China, especially for protecting soil against serious erosion and sandstorms. However, inherent severe water shortages, drought stresses, and increasing water use competition greatly restrict the reforestation. Notably, the water‐yield reduction after afforestation generates intense debate about the correct approach to afforestation and forest management in dry‐land areas. However, most studies on water‐yield reduction of forests have been at catchment scales, and there are few studies of the response of total evapotranspiration (ET) and its partitioning to vegetation structure change. This motivates us to learn the linkage between hydrological processes and vegetation structure in slope ecosystems. Therefore, an ecohydrological study was carried out by measuring the individual items of water balance on sloping plots covered by different vegetation types in the semiarid Liupan Mountains of northwest China. The ratio of precipitation consumed as ET was about 60% for grassland, 93% for shrubs, and >95% for forestland. Thus, the water yield was very low, site‐specific, and sensitive to vegetation change. Conversion of grassland to forest decreased the annual water yield from slope by 50‐100 mm. In certain periods, the plantations at lower slopes even consumed the runon from upper slopes. Reducing the density of forests and shrubs by thinning was not an efficient approach to minimize water use. Leaf area index was a better indicator than plant density to relate ET to vegetation structure and to evaluate the soil water carrying capacity for vegetation (i.e., the maximum amount of vegetation that can be supported by the available soil water for an extended time). Selecting proper vegetation types and plant species, based on site soil water condition, may be more effective than the forest density regulation to minimize water‐yield reduction by vegetation coverage increase and notably by reforestation. Finally, the focuses in future research to improve the forest‐water relations in dry‐land areas are recommended as follows: vegetation growth dynamics driven by environment especially water conditions, coupling of ecological and hydrological processes, further development of distributed ecohydrological models, quantitative relation of eco‐water quota of ecosystems with vegetation structures, multi‐scaled evaluation of soil water carrying capacity for vegetation, and the development of widely applicable decision support tools.     

Silvicultural management in maintaining biodiversity and resistance of forests in Europe-boreal zone: case Finland

The majority of untouched natural boreal forests have been regenerated through large catastrophes, occurring by intervals between 50 and 100 years. Storm and fire will open the landscape, result in a huge amount of dead or dying trees and let the pioneer tree species germinate. These processes are the guideline for Finnish forest management today.The main focus by maintaining the biodiversity in Finnish boreal forest zone is directed to managed forests. Nature-orientated silviculture on stand level is practised. The site type classification, a reflection of the modern concept of biodiversity and developed by Cajander early in 1900s, on the basis of natural vegetation composition of the site, has the central role by choosing tree species, regeneration methods and thinning procedure, and reflects also on the site productivity. The small size of stands, the abundance of natural seedlings in planted stands and the popularity of mixed stands have a positive impact on biodiversity of forests. The protection of small-sized valuable habitats in commercially managed stands, the leaving of retention trees standing and lying in the forest in all phases of the rotation, are activities made for biodiversity. Many insects and fungi are adapted to catastrophes and so they can survive in single stems left on regeneration areas. Maintaining the biodiversity in multifunctional forests is also supported by the new forest legislation and by the criteria of Finnish Forest Certification System.     

Can intensive management increase carbon storage in forests?

A possible response to increasing atmospheric CO₂ concentration is to attempt to increase the amount of carbon stored in terrestrial vegetation. One approach to increasing the size of the terrestrial carbon sink is to increase the growth of forests by utilizing intensive forest management practices. This article uses data from the literature and from forest growth and yield models to analyze the impact of three management practices on carbon storage: thinning, fertilization, and control of competing vegetation. Using Douglas-fir (Pseudotsuga menziesii) and loblolly pine (Pinus taeda) as example species, results from experiments with computer simulation models suggest that, for these two species, thinning generally does not increase carbon storage and may actually cause a decrease. The exception is thinning of very dense young stands. Fertilization generally increases carbon storage, although the response can be quite variable. The largest gains in carbon storage are likely to come from fertilizing lower-quality sites and from fertilizing thinned or less dense stands. Forests usually show increased growth in response to fertilization over a wide range of ages. Simulation of the growth of loblolly pine indicates that controlling competing vegetation at an early age helps to maximize stand growth and carbon storage. The research described in this article has been funded by the US Environmental Protection Agency. This document has been prepared at the EPA Environmental Research Laboratory in Corvallis, Oregon, through contract number 68-C8-0006 to NSI Technology Inc. It has been subjected to the agency’s peer and administrative review and approved for publication. Mention of trade names or commercial products does not constitute endorsement or recommendation for use.