生态模型的灵敏度分析

摘要/Abstract

摘要： 灵敏度分析用于定性或定量地评价模型参数误差对模型结果产生的影响,是模型参数化过程和模型校正过程中的有用工具,具有重要的生态学意义.灵敏度分析包括局部灵敏度分析和全局灵敏度分析.局部灵敏度分析只检验单个参数的变化对模型结果的影响程度;全局灵敏度分析则检验多个参数的变化对模型运行结果总的影响,并分析每一个参数及其参数之间相互作用对模型结果的影响.目前,在对生态模型的灵敏度分析中,越来越倾向于使用全局灵敏度分析的方法.但国内仍多采用局部灵敏度分析方法,很少采用全局灵敏度分析方法.文中详细论述了局部灵敏分析和全局灵敏度分析的主要方法(一次变换法、多元回归法、Morris法、Sobol'法、傅里叶幅度灵敏度检验法和傅里叶幅度灵敏度检验扩展法),希望能为国内生态模型的发展提供一个比较完善的灵敏度分析方法库.结合国内外的灵敏度分析发展现状,指出联合灵敏度研究、灵敏度共性研究及空间直观景观模型的灵敏度分析将为生态模型灵敏度分析研究中的热点和难点.

关键词: 生态模型, 灵敏度分析, 局部灵敏度分析, 全局灵敏度分析, 空间直观景观模型

Abstract: Sensitivity analysis is used to qualitatively or quantitatively apportion the variation of model output to different source of variation.It is a very useful tool in model parameterization and calibration,and has important ecological significance by identifying the governing factors for a certain ecological process simulated.There are two schools of sensitivity analysis, local sensitivity analysis and global sensitivity analysis.The former examines the local response of the output(s) by varying input parameters one at a time,holding other parameters to a central value; and the latter examines the global response (averaged over the variation of all the parameters) of model output(s) by exploring a finite (or even an infinite) region.Since it is very easy to conduct local sensitivity analysis,it is very popular in ecological models.However,local sensitivity analysis is not computationally effective,because it can only get the sensitivity of a single parameter at a time.It can not take into consideration the effect of interaction of different parameters.Additionally,the value of other parameters will affect the sensitivity of the parameter specified.In view of this,global sensitivity analysis is increasingly preferred to local sensitivity in recent years.However,for most of the ecological modeling study published in Chinese,only local sensitivity analysis is conducted.To provide a toolbox of alternative sensitivity analysis algorithms for ecological model development in future study,we reviewed the main methods of both local sensitivity analysis and global sensitivity analysis,including one at a time method,multivariate regression,Morris' method,Sobol's method,Fourier Amplitude Sensitivity Analysis,and Extended Fourier Amplitude Sensitivity Analysis.Based on the state-of-the-art research on sensitivity analysis,the sensitivity of the interaction between two or more than two model parameters,the sensitivity of common model parameters in a set of models,and the sensitivity analysis in spatially explicit landscape model simulation were identified as the key areas and difficulties of future study on sensitivity analysis in ecological modeling.

Key words: Ecological model, Sensitivity analysis, Local sensitivity analysis, Global sensitivity analysis, Spatially explicit landscape model

中图分类号:

徐崇刚, 胡远满, 常禹, 姜艳, 李秀珍, 布仁仓, 贺红士. 生态模型的灵敏度分析[J]. 应用生态学报, 2004, (6): 1056-1062.

XU Chonggang, HU Yuanman, CHANG Yu, JIANG Yan, LI Xiuzhen, BU Renchang, HE Hongshi . Sensitivity analysis in ecological modeling[J]. Chinese Journal of Applied Ecology, 2004, (6): 1056-1062.

参考文献

[1] Battaglia M,Sands P.1998.Application of sensitivity analysis to a model of Eucalyptus globulus plantation productivity.Ecol Model,111:237～259
[2] Beres DL,Hawkins DM.2001.Plackett-Burman technique for sensitivity analysis of many-parametered models.Ecol Model,141:171～183
[3] Campolongo F,Saltelli A.1997.Sensitivity analysis of an environmental model:An application of different analysis methods.Reliab Engin Syst Safety,57:49～69
[4] Carlson DH,Thruow TL.1996.Comprehensive evaluation of improved SPUR model(SPUR～91).Ecol Model,85:229～240
[5] Clark JS,Carpenter SR,Barber M,et al.2001.Ecological forecasts:An emerging imperative.Science,27:657～660
[6] Crosetto M,Ruiz JAM,Crippa B.2001.Uncertainty propagation in models driven by remotely sensed data.Rem Sens Environ,76:373～385
[7] Crosetto M,Tarantola S.2001.Uncertainty and sensitivity analysis:Tools for GIS-based model implementation.INT J Geogr Iinform Sci,15(5):415～437
[8] Crosetto M,Tarantola S,Saltelli A.2000.Sensitivity and uncertainty analysis in spatial modeling based on GIS.Agric Ecosyst Environ,81:71～79
[9] Cryer SA,Havens,PL.1999.Regional sensitivity analysis using a fractional factorial method for the USDA model GLEAMS.Environ Model Software,14:613～624
[10] Cukier RI,Fortuin CM,Shuler KE,et al.1973.Study of the sensitivity of coupled reaction systems to uncertainties in rate coefficients I.Theory.J Chem Phys,59:3873～3878
[11] Cukier RI,Levine HB,Shuler KE.1978.Nonlinear sensitivity analysis of multiparameter model systems.J Comput Phys,26:1～42
[12] Cukier RI,Schaibly JH,Shuler KE.1975.Study of the sensitivity of coupled reaction systems to uncertainties in rate coefficients Ⅲ.Analysis of the approximations.J Chem Phys,63:1140～1149
[13] Downing DJ,Gardner RH,Hoffman FO.1985.An examination of response-surface methodologies for uncertainty analysis in assessment models.Technometrics,27:151～163
[14] Dunning JB,Stewart DJ,Danielson BJ,et al.1995.Spatially explicit population models:Current forms and future uses.Ecol Appl,5:3～11
[15] Fennel K,Losch M,Schroter J,et al.2001.Testing a marine ecosystem model:Sensitivity analysis and parameter optimization.J Mar Syst,28:45～63
[16] Gardner RH,O'Neill RV,Mankin JB,et al.1981.A comparison of sensitivity analysis and error analysis based on a stream ecosystem model.Ecol Model,12:173～190
[17] Gentil S,Blake G.1981.Validation of complex ecosystem models.Ecol Model,14:21～38
[18] Halvorsen E,Pedersen OP,Slagstad D,et al.2001.Microzooplankton and mesozooplankton in an upwelling filament off Galicia:Modelling and sensitivity analysis of the linkages and their impact on the carbon dynamics.Prog Oceanogr,51:499～513
[19] Helton JC.1993.Uncertainty and sensitivity analysis techniques for use in performance assessment for radioactive waste disposal.Reliab Engin Syst Safety,42:327～367
[20] Henderson-Sellers B,Henderson-sellers A.1996.Sensitivity evaluation of environmental models using fractional factorial experimentation.Ecol Model,86:291～295
[21] Homma T,Saltelli A.1996.Importance measures in Global sensitivity analysis of Nonlinear models.Reliab Engin Syst Safety,52:1～17
[22] Iman RL,Hora SC.1990.A robust measure of uncertainty importance for use in fault tree system analysis.Risk Ana,10:401～406
[23] Iman RL.1999.Latin hypercube sampling.Encyclopedia of Statistical Science.Vol.3.New York:Wiley.408～411
[24] Ishigami T,Homma T.1990.An importance quantification technique in uncertainty analysis for computer models.Proceedings of the ISUMA'90,First International Symposium on Uncertainty Modelling and Analysis,University of Maryland,USA.New York:IEEE Computer Society.398～403
[25] K?ihler P,Wirtz KW.2002.Linear understanding of a huge aquatic ecosystem model using a group-collecting sensitivity analysis.Environ Model Software,17:613～625
[26] Koda M,McRae GJ,Seinfeld JH.1979.Automatic sensitivity analysis of kinetic mechanisms.Intern J Chem Kinet,11:427～444
[27] Lu YC,Mohanty S.2001.Sensitivity analysis of a complex,proposed geologic waste disposal system using the Fourier Amplitude sensitivity test method.Reliab Engin Syst Safety,72:275～291
[28] Mahamah DS.1988.Simplified sensitivity analysis applied to a nutrient-biomass model.Ecol Model,42:103～109
[29] Majkowski J,Tidgeway JM,Miller DR.1981.Multiplicative sensitivity analysis and its role in development of simulation models.Ecol Model,12:191～208
[30] McCarthy MA,Burgman MA,Ferson S.1995.Sensitivity analysis for models of population viability.Biol Conserv,73:93～100
[31] McKay MD,Beckman RJ,Conover WJ.1979.A comparison of three methods fro selecting values of input variables in the analysis of output from a computer code.Technometrics,21:239～245
[32] Mladenoff DJ,He HS.1999.Design and behavior of LANDIS,an object-oriented model of forest landscape disturbance and succession.In:Mladenoff DJ,Baker WL,eds.Spatial Modeling of Forest Landscape Change:Approaches and Applications.Cambridge:Cambridge University Press.125～162
[33] Morris MD.1991.Factorial sampling plans for preliminary computational experiments.Technometrics,33:161～174
[34] Nathan R,Safriel UN,Noy-Meir I.2001.Field validation and sensitivity analysis of a mechanistic model for tree seed dispersal by wind.Ecology,82(2):374～388
[35] Omlin M,Brun R,Reichert P.2001.Biogeochemical model of Lake Zurich:Sensitivity,identifiability and uncertainty analysis.Ecol Model,141:105～123
[36] Pacala,SW,Canham CC,Saponara J,et al.1996.Forest models defined by field measurement:Estimation,error analysis and dynamics.Ecol Monogr,66:1～43
[37] Pastres R,Chan K,Solidoro C,et al.1999.Global sensitivity analysis of a shallow water 3-D eutrophication model.Comput Phys Commun,117:62～74
[38] Pastres R,Franco D,Pecenik G,et al.1997.Local sensitivity analysis of a distributed parameters water quality model.Reliab Engin Syst Safety,57:21～30
[39] Ratto M,Tarantola S,Saltelli A.2001.Sensitivity analysis in model calibration:GSA-GLUE approach.Comput Phys Commun,136:212～224
[40] Recknagel F.1984.A comprehensive sensitivity analysis for an ecological simulation model.Ecol Model,26:77～96
[41] Saltelli A,Anderes TH,Homma T.1993.Sensitivity analysis of model output:An investigation of New Techniques.Comput Statist Data Anal,15:211～238
[42] Saltelli A,Chan K,Scott M.2000.eds.Sensitivity Analysis,Probability and Statistics Series.New York:John Wiley & Sons.
[43] Saltelli A,Tarantola S,Chan KPS.1999.A quantitative model-independent method for global sensitivity analysis of model output.Technometrics,41:39～56
[44] Salvador R,Pinol J,Tarantola S,et al.2001.Global sensitivity analysis and scale effects of a fire propagation model used over Mediterranean shrublands.Ecol Model,136:175～189
[45] Sharov AA.1995.Non-linear sensitivity analysis as a tool for explaining population dynamics in a life-system conceptual framework.Comput Electron Agric,13:103～113
[46] Shirley MDF,Rushton PS,Yong AG,et al.2001.Simulating the long-term dynamics of slug populations:A process-based modeling approach for pest control.J Appl Ecol,38:401～411
[47] Sobol' IM.1993.Sensitivity Estimates for nonlinear mathematical models.Math Model Comput Exp,1:407～414
[48] Tarantola S,Giglioli N,Jesinghaus J,et al.2002.Can global sensitivity analysis steer the implementation of models for environmental assessments and decision-making?Stochastic Environ Res Risk Asses,16:63～76
[49] Tomovic ' R.1963.Sensitivity Analysis of Dynamic Systems.New York:McGraw-Hill.142
[50] Waide JB,Webster JR.1976.Engineering systems analysis.In:Patten BC,eds.Systems Analysis and Simulation in Ecology.4.New York:Academic Press.329～371
[51] Zhang H(张宏),Fang Z-L(樊自立).2000.Studies on a NPP model of salinized meadow in the north of Tarim basin.Acta Phytoecol Sin(植物生态学报),24(1):13～17(in Chinese)
[52] Zhu Z-R(祝增荣),Hang C-W(黄次伟),Chen J-A(程家安),et al.1994.Simulation analysis of white-backed flanthopper population dynamics on first season rive in Zhejiang Province,China.Acta Ecol Sin(生态学报),14(2):188～195(in Chinese)