Global sensitivity analysis of a modified CENTURY model for simulating impacts of harvesting fine woody biomass for bioenergy.

Modeling the long-term effects of intensive forest biomass harvesting scenarios over time, including the more complete removal of tree tops and branches, is a scientific and policy need. Yet, due to our incomplete understanding about complex forest ecosystems, model simulations are to various degrees uncertain. In this study we first modified a well-evaluated and widely used ecosystem model – CENTURY 4.5 – to model management scenarios that retain various sizes and quantities of small-diameter woody material after intensive biomass harvests. Second, we used a global sensitivity analysis approach to evaluate the sensitivity of nine model outputs to 55 parameters, grouped into 17 factors. The values of the parameters were generated with a normal distribution and sampled with the extended Fourier amplitude sensitivity test. Our analysis indicated that within a harvest rotation, the model output sensitivity varied over years in response to different factors. The model was most sensitive to factors consisting of temperature effects on potential production as well as N deposition and non-symbiotic N fixation. In response to the uncertain parameter values, the model simulation revealed that outputs of net N mineralization rates in slow and passive soil organic matter pools had the highest uncertainties. However, due to the very low fraction of the N supplied from these two pools, forest production and other simulations were not strongly affected, ending with overall variations less than 6%. Ultimately, this study exhibits a novel approach in modeling the effects of harvesting fine woody debris for bioenergy on long-term ecosystem C and N cycles, and illustrates that sensitivity testing the most uncertain parameters is crucial for minimizing model uncertainty.

[1]  A. Mäkelä,et al.  Optimal co-allocation of carbon and nitrogen in a forest stand at steady state. , 2008, The New phytologist.

[2]  A. Saltelli,et al.  An alternative way to compute Fourier amplitude sensitivity test (FAST) , 1998 .

[3]  Saltelli Andrea,et al.  Global Sensitivity Analysis: The Primer , 2008 .

[4]  A. Saltelli,et al.  A quantitative model-independent method for global sensitivity analysis of model output , 1999 .

[5]  Igor Linkov,et al.  Exploring vulnerability of coastal habitats to sea level rise through global sensitivity and uncertainty analyses , 2011, Environ. Model. Softw..

[6]  K. Rees,et al.  Contribution of fine roots to ecosystem biomass and net primary production in black spruce, aspen, and jack pine forests in Saskatchewan , 2006 .

[7]  A. Saltelli,et al.  The role of sensitivity analysis in ecological modelling , 2007 .

[8]  Yowhan Son.  Non-symbiotic nitrogen fixation in forest ecosystems , 2001, Ecological Research.

[9]  Robert J. Scholes,et al.  Observations and modeling of biomass and soil organic matter dynamics for the grassland biome worldwide , 1993 .

[10]  K. Wilson,et al.  OAK FOREST CARBON AND WATER SIMULATIONS: MODEL INTERCOMPARISONS AND EVALUATIONS AGAINST INDEPENDENT DATA , 2004 .

[11]  Bert Bredeweg,et al.  Representing and managing uncertainty in qualitative ecological models , 2009, Ecol. Informatics.

[12]  M. Kaufmann,et al.  Effects of water and temperature , 1985 .

[13]  C. Fortuin,et al.  Study of the sensitivity of coupled reaction systems to uncertainties in rate coefficients. I Theory , 1973 .

[14]  Changhui Peng,et al.  Modelling the response of net primary productivity (NPP) of boreal forest ecosystems to changes in climate and fire disturbance regimes , 1999 .

[15]  Ryuichi Hirata,et al.  Simulating carbon and water cycles of larch forests in East Asia by the BIOME-BGC model with AsiaFlux data , 2009 .

[16]  M. G. Ryan,et al.  Factors controlling Eucalyptus productivity: How water availability and stand structure alter production and carbon allocation , 2010 .

[17]  Shuangzhe Liu,et al.  Global Sensitivity Analysis: The Primer by Andrea Saltelli, Marco Ratto, Terry Andres, Francesca Campolongo, Jessica Cariboni, Debora Gatelli, Michaela Saisana, Stefano Tarantola , 2008 .

[18]  Jianguo Wu,et al.  UNCERTAINTY ANALYSIS IN ECOLOGICAL STUDIES:AN OVERVIEW , 2006 .

[19]  Stefano Tarantola,et al.  Can global sensitivity analysis steer the implementation of models for environmental assessments and decision-making? , 2002 .

[20]  J. Probst,et al.  Aspen ecosystem properties in the Upper Great Lakes , 1991 .

[21]  J. Yarie,et al.  Carbon balance of the taiga forest within Alaska: present and future , 2002 .

[22]  W. Parton,et al.  Dynamics of C, N, P and S in grassland soils: a model , 1988 .

[23]  Andrew D. Richardson,et al.  Phenology of a northern hardwood forest canopy , 2006 .

[24]  S. Pacala,et al.  Evolutionarily Stable Strategy Carbon Allocation to Foliage, Wood, and Fine Roots in Trees Competing for Light and Nitrogen: An Analytically Tractable, Individual-Based Model and Quantitative Comparisons to Data , 2011, The American Naturalist.

[25]  D. Mladenoff.  LANDIS and forest landscape models , 2004 .

[26]  D. J. Lodge,et al.  Hurricane effects on soil organic matter dynamics and forest production in the Luquillo experimental forest, Puerto Rico : results of simulation modeling , 1991 .

[27]  Y. Xu,et al.  Projecting Climate Change Effects on Forest Net Primary Productivity in Subtropical Louisiana, USA , 2011, AMBIO.

[28]  Nicholas Skowronski,et al.  Simulation and sensitivity analysis of carbon storage and fluxes in the New Jersey Pinelands , 2011, Environ. Model. Softw..

[29]  Maria K. Janowiak,et al.  Promoting Ecological Sustainability in Woody Biomass Harvesting , 2010, Journal of Forestry.

[30]  S. Fraver,et al.  Forest management for mitigation and adaptation to climate change: Insights from long-term silviculture experiments , 2011 .

[31]  S. T. Gower,et al.  Simulated long‐term effects of harvest and biomass residue removal on soil carbon and nitrogen content and productivity for two Upper Great Lakes forest ecosystems , 2011 .

[32]  W. Parton,et al.  Analysis of factors controlling soil organic matter levels in Great Plains grasslands , 1987 .

[33]  Vemap Participants.  Vegetation/ecosystem modeling and analysis project: Comparing biogeography and biogeochemistry models in a continental-scale study of terrestrial ecosystem responses to climate change and CO2 doubling , 1995 .

[34]  N. Debyle,et al.  Aspen : Ecology and Management in the Western United States , 1985 .

[35]  Kristofer D. Johnson,et al.  Short- and long-term responses of total soil organic carbon to harvesting in a northern hardwood forest , 2010 .

[36]  Juan A. Blanco,et al.  Forecasting Forest Futures: A Hybrid Modelling Approach to the Assessment of Sustainability of Forest Ecosystems and their Values , 2010 .

[37]  James S. Clark.  Uncertainty in Ecological Inference and Forecasting1 , 2003 .

[38]  J. A. Trofymow,et al.  Estimating branch production in trembling aspen, Douglas-fir, jack pine, black spruce, and balsam fir , 2007 .

[39]  Orie L. Loucks,et al.  Scaling and uncertainty analysis in ecology : methods and applications , 2006 .

[40]  K. Pregitzer,et al.  Soil respiration, root biomass, and root turnover following long-term exposure of northern forests to elevated atmospheric CO2 and tropospheric O3. , 2008, The New phytologist.

[41]  Stefano Tarantola,et al.  Sensitivity Analysis in Practice: A Guide to Assessing Scientific Models , 2004 .

[42]  Elizabeth B. Harper,et al.  Global sensitivity analysis for complex ecological models: a case study of riparian cottonwood population dynamics. , 2011, Ecological applications : a publication of the Ecological Society of America.

[43]  R. Birdsey,et al.  The effects of forest harvest intensity in combination with wind disturbance on carbon dynamics in Lake States Mesic Forests , 2011 .

[44]  R. Katz.  Techniques for estimating uncertainty in climate change scenarios and impact studies , 2002 .

[45]  Tom Andersen,et al.  MyLake—A multi-year lake simulation model code suitable for uncertainty and sensitivity analysis simulations , 2007 .