空间直观森林景观模型空间不确定性的地统计学模拟

doi:10.7523/j.issn.2095-6134.2005.4.007

摘要/Abstract

摘要：

为了用尽可能小的蒙特卡罗模拟样本来反映模型模拟结果中的不确定性,把拉丁超几何体采样引入地统计随机模拟的LU分解算法.首先把拉丁超几何体采样与普通随机采样在LU分解算法中的表现进行比较,然后把基于拉丁超几何体采样的LU分解法应用于空间直观森林景观模型LANDIS的模拟.结果表明,与普通随机采样相比,拉丁超几何体采样能捕获更多的不确定性,特别是在蒙特卡罗模拟次数较少时.LANDIS模型的模拟结果表明,由地统计学随机模拟所引入的不确定性在象元尺度上随模拟时间增加而增加,但是在景观尺度上并没有受很大影响.这表明由地统计学随机模拟所参数化的LANDIS模型可用于大时空尺度森林景观变化的预测

关键词: 不确定性, 克吕格插值, 地统计随机模拟, LU分解, 拉丁超几何体采样, 空间直观森林景观模型

Abstract:

We introduced an effective sampling method (Latin Hypercube sampling) into a stochastic simulation algorithm (LU decomposition simulation). Latin Hypercube sampling is first compared with a common sampling procedure (random simple sampling) in LU decomposition simulation. Then it is applied to the investigation ofuncertainty in the simulation results of a spatially explicit forest model, LANDIS. Results showed that Latin Hypercube sampling can capture more variability in the sample space than simple random sampling especially when the number of simulations is small. Simple as the application is, it gives us general insights about which model results are robust given the uncertainty introduced by interpolation. Application results showed that LANDIS simulation results at the landscape level (species percent area and their spatial pattern measured by an aggregation index) is not sensitive to the uncertainty in species age cohort information at the cell level produced by geostatistical stochastic simulation algorithms. This suggests that LANDIS can be used to predict the forest landscape change at broad spatial and temporal scales even if exhaust species age cohort information at each cell is not available.

Key words: uncertainty, Kriging interpolation, geostatistical stochastic simulation, LU decomposition, LatinHypercube sampling, spatially explicit forest landscape model

中图分类号:

S718.5

徐崇刚, 胡远满, 常禹, 李秀珍, 布仁仓, 贺红士. 空间直观森林景观模型空间不确定性的地统计学模拟[J]. 中国科学院大学学报, 2005, 22(4): 436-446.

XU Chong-Gang, HU Yuan-Man, CHANG Yu, LI Xiu-Zhen, BU Ren-Cang, HE Hong-Shi. Geostatistical Modeling of Spatial Uncertainty in a Spatially Explicit Forest Landscape Model Simulation[J]. , 2005, 22(4): 436-446.

参考文献

[1] Baker WL, Egbert SL, Frazier GF. A spatial model for studying the effects of climatic change on the structure of landscapes subject to large disturbances. Ecol. Model. 1991, 56 : 109～125

[2] Gardner RH, Hargrove WW, Turner MG, et al. Climate change, disturbances and landscape dynamics. In : Walker BH, Steffen WL (Eds.). Global Change and Terrestrial Ecosystems. Cambridge : Cambridge University Press, 1996. 149～172

[3] Gustafson EJ,Crow TR. Modeling the effects of forest harvesting on landscape structure and the spatial distribution of cowbird brood parasitism.Landscape Ecol. 1994, 9 : 237～248

[4] Gustafson EJ, Crow TR. Simulating the effects of alternative forest management strategies on landscape structure. J. Environ. Manage. 1996,46 : 77～94

[5] Gustafson EJ, Crow TR. HARVEST: linking timber harvesting strategies to landscape patterns. In : Klopatek JM, Gardner RH (Eds.). Landscape Ecological Analysis : Issues and Application. New York : Springer2Verlag, 1999. 309～332

[6] Mladenoff DJ, Host GE, Boeder J, et al. LANDIS: a spatial model of forest landscape disturbance, succession and management. In : GoodchildMR, Steyaert LT, Parks BO(Eds.). GIS and Environmental Modeling : Progress and Research Issues. Fort Collins : GISWorld Books, 1996. 175～180

[7] Mladenoff DJ, He HS. Design and behavior of LANDIS, an object2oriented model of forest landscape disturbance and succession. In : MlandenoffDJ, Baker WL (Eds.). Advances in Spatial Modeling of Forest Landscape Change : Approaches and Applications. Cambridge : CambridgeUniversity Press, 1999. 1～13

[8] Urban DL, Acevedo MF, Garman SL. Scaling fine2scale processes to large scale patterns using models derived from models : meta2models. In :Mladenoff DJ, Baker WL (Eds.). Spatial Modeling of Forest Landscape Change : Approaches and Applications. Cambridge : Cambridge University Press, 1999. 70～98

[9] Wallin DO, Swanson FJ, Marks B. Landscape pattern response to changes in pattern generation rules : land2use legacies in forestry. Ecol. Appl. 1994, 4 : 569～580

[10] ESRI. Using ArcGISTM Geostatistical Analyst. USA,CA,Redlands :ESRI InC. 2001

[11] Biondi F, Meyers DE, Avery CC. Geostatistically modeling stem size and increment in an old growth forest. Canadian Journal of Forest Research,1994, 24 : 1354～1638

[12] Goovaerts P. Geostatistics for Natural Resources Evaluation. New York : Oxford University Press, 1997. 1～483

[13] Isaaks EH. The application of Monte Carlo methods to the analysis of spatially correlated data : [Ph. D. thesis]. Stanford University, 1990

[14] Pachepsky Y,Acock B. Stochastic imaging of soil parameters to assess variability and uncertainty of crop yield estimates. Geoderma, 1998, 85 :

213～229

[15] Finke PA, Wladis D, Dros J, et al. Quantification and simulation of errors in categorical data for uncertainty analysis of soil acidification modelling. Geoderma, 1999, 93 : 177～194

[16] Goovaerts P. Geostatistical modelling of uncertainty in soil science. Geoderma, 2001, 103 : 3～26

[17] Saito H, Goovaerts P. Accounting for source location and transport direction into geostatistical prediction of contaminants. Environmental Science &Technology, 2001, 35 : 4823～4829

[18] Van Meirvenne M, Goovaerts P. Evaluating the probability of exceeding a site2specific soil cadmium contamination threshold. Geoderma, 2001,102 : 75～100

[19] Viscarra Rossel RA, Goovaerts P, McBratney AB. Assessment of the production and economic risks of site2specific liming using geostatistical uncertainty modelling. Environmetrics, 2001, 12 (8) : 699～711

[20] Xu CG, He HS, Hu YM, et al. Assessing the effect of cell level uncertainty on a forest landscape model simulation in northeastern China. Ecological Modelling,2004,181 (10) :57～72

[21] Deutsch CV. Algorithmically2defined random function models. In : Dimitrakopoulos R (Ed.). Geostatistics for the Next Century. Dordrecht :Kluwer Academic Publishing, 1994. 422～435

[22] Gotway CA, Rutherford BM. Stochastic simulation for imaging spatial uncertainty : comparison and evaluation of available algorithms. In : Amstrong M, Dowd PA (Eds.). Geostatistical Simulations. Dordrecht : Kluwer Academic Publishing, 1994. 1～21

[23] Srivastava MR. An overview of stochastic spatial simulation. In : Mowrer HT, Czapleweski RL, Hamre RH (Eds.). Spatial Accuracy Assessment in Natural Resources and Environmental Sciences : Second International Symposium, Genereal Tchnical Report RM2GTR2277. US department of Agriculture, Forest Service, Fort Collins. 1996. 283～295

[24] Goovaerts P. Geostatistics in soil science : state2of2the2art and perspectives. Geoderma, 1999, 89 : 1～45

[25] Gómez2Hernández JJ. Issues on environmental risk assessment. In : Baafi EY, Schofield NA (Eds.). Geostatistics Wollongong’96. Dordrecht :Kluwer Academic Publishing, 1997. 15～26

[26] Alabert FG. The practice of fast conditional simulation through the LU decomposition of the covariance matrix. Mathematical Geology, 1987, 19 :369～386

[27] Davis MW. Production of conditional simulations via the LU decomposition of the covariance matrix. Mathematical geology, 1987, 19 : 91～98

[28] Helton JC, Davis FJ. Latin hypercube sampling and the propagation of uncertainty in analyses of complex systems. Reliability Engineering and System Safety, 2003, 81 : 23～69

[29] McKay MD, Beckman RJ,Conover WJ. A comparison of three methods for selecting values of input variables in the analysis of output from a computer code. Technometrics, 1979, 21 : 239～245

[30] Iman RL. Latin hypercube sampling. Encyclopedia of statistical science update (Vol. 3). New York : Wiley, 1999. 408～411

[31] Wolter PT, Mladenoff FJ, Host GE, et al. Improved forest classification in the Northern Lake States using multi2temporal Landsat imagerPhotogrammetric Engineering and Remote Sensing, 1995, 61 :1129～1143

[32] He HS, Maladenoff DJ, Radeloff VC, et al. Integration of GIS data classified satellite imagery for regional forest assessment. Ecol. Appl. 1998,8 (4) : 1072～1083

[33] He HS, Mladenoff DJ, Boeder J. An object2oriented forest landscape model and its representation of tree species. Ecol. Model. 1999, 119 : 1～19

[34] He HS, Mladenoff DJ. Spatially explicit and stochastic simulation of forest landscape fire disturbance and succession. Ecology, 1999, 80 : 81～99

[35] He HS, Mladenoff DJ. The effects of seed dispersal on the simulation of long2term forest landscape change. Ecosystems, 1999, 2 : 308～319

[36] Gustafson EJ, Shifley SR, Mladenoff DJ, et al. Spatial simulation of forest succession and timber harvesting using LANDIS. Can. J. Forest Res. 2000, 30 : 32～43

[37] He HS, Dezonia B, Mladenoff DJ. An aggregation index (AI) to quantify spatial patterns of landscapes. Landscape Ecol. 2000, 15 : 591～601

[38] Mladenoff DJ, DeZonia B. APACK 2. 14 Users Guide. Department of Forest Ecology and Management. USA,WI : University of Wisconsin2Madison, 2000

[39] Vargas2Guzman JA, Dimitrakopoulos R. Conditional simulation of randomfields by successive residuals. Mathematical geology, 2002, 34 (5) : 597～610

[40] O’Neil, RV. Perspectives in hierarchy and scale. In : Roughgarden J, May RM, Levin SA(Eds.). Perspectives in Ecological Theory. Princeton :Princeton University Press, 1989. 140～156

[41] Wiens JA. Spatial scaling in ecology. Functional Ecology, 1989, 3 : 385～397

[42] Wiens JA, Milne TT. Scaling of“landscapes”in ecology, or landscape ecology from a beetle’s perspective. Landscape Ecol. 1989, 3 : 87～96

[43] Levin SA. The problem of pattern and scale in Ecology. Ecology, 1992, 73 : 1943～1967

[44] Levin SA, Grenfell B, Hastings A, et al. Mathematical and computational challenges in population biology and ecosystems science. Science,1997, 275 : 334～342

[45] Clark JS, Carpenter SR, Barber M, et al. Ecological forecasts : an emerging imperative. Science, 2001, 27 : 657～660

[46] Pacala SW, Canham CD, Saponara J, et al. Forest models defined by field measurements : estimation, error analysis and dynamics. Ecol. Monogr. 1996, 66 : 1～43

[47] Deutschman DH, Levin SA, Pacala SW. Error propagation in a forest succession model :the role of fine2scale heterogeneity in light. Ecology,1999, 80 : 1927～1943