Arne Pommerening's Forest Biometrics & Quantitative Ecology Lab

From Arne Pommerening's Forest Biometrics Lab
Revision as of 17:06, 8 June 2018 by Arne Pommerening (Talk | contribs) (Upcoming conferences and workshops)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search
Screenshot of a NetLogo implementation of the individual-based shot-noise model published in Pommerening et al. (2011)‎ Individual-based models are specialised on interactions between organisms and I am using them for examining how marked point patterns of trees evolve through time.
Philodendron International


This website and all associated computer code are part of a research area, which is dedicated to systems analysis of forest ecosystems, an important branch of ecological statistics integrating research on forest structure, tree growth analysis, monitoring and modelling.

All data, that are collected in forest ecosystems, have not only a temporal but also a spatial dimension. The properties of the whole system “forest”, e.g. forest growth and interactions between trees, to a large degree depend on the structure of this system. This has been known in forest science for a long time and has given rise to traditional terms such as “growing space” and “initial spacing”. In the last decades, new methods were developed in the statistical fields of point process statistics, geostatistics and random set statistics. These allow better and more detailed research of the interplay between spatial patterns and ecological processes. Apart from statistics and mathematics also other research fields such as physics and materials science have contributed to structural research. The latter subject area has coined the term of structure-property relationships. According to this term ecological processes not only leave traces as spatial patterns, but the spatial structure of materials or of a forest also determines to a large degree the properties of the system under study. This suggests that the results of forest management, e.g. forest products, but also natural regeneration responding to silvicultural systems, are examples of such system properties, which are largely determined by its structure. Structure-property relationships also apply to ecosystem services such as growth and yield and the provision of habitats and recreation. Competition and survival of trees are such properties as well as the sampling error of forest resource inventories, which is strongly correlated with spatial forest structure. Spatial statistics and inventory research are therefore closely related. Also any impact on forests - whether natural or human-induced - is primarily a change of forest structure. It is crucial for my research vision that the quantification, the understanding and the modelling of (spatial) woodland structure and its temporal evolution play a decisive role for the development of sampling and monitoring designs as well as for forest modelling. Spatial statistics and sampling theory share many overlaps which can be used as synergies. And there is certainly no sustainable silviculture and forest management without a proper understanding of woodland structure as every forest practitioner can tell you. The publications below, particularly those relating to reconstruction (Pommerening, 2006; Pommerening and Stoyan, 2008) prove and illustrate this research vision.


A group photo of my research team at Umeå in 2015: Arne Pommerening, Anton Grafström, Xin Zhao, Anders Muszta, Kenneth Nyström and Jaime Uría Díez (from left to right).

I work as a Professor in Forest Biometrics at the Swedish University of Agricultural Sciences (SLU) in Umeå. I am a theoretical forest scientist and my research areas include woodland structure analysis and modelling, spatio-temporal dynamics of plant point patterns, individual-based modelling with a focus on plant interactions, plant growth analysis, methods of quantifying and monitoring biodiversity and the analysis of human behaviour of selecting trees. Recently my group and I have also developed a keen interest in the question of how marked point patterns based on mapped plant communities evolve through time. The ontogenesis of such patterns can be studied by using individual-based models, which is an intriguing field of research. Much of my research is in quantitative ecology including computer-based simulation experiments. My research is strongly interdisciplinary and international.

My Chair is described in greater detail here. I am maintaining a webblog on current issues in Forest Biometrics and you are welcome to follow the discussions on this address. I am also the director of the Centre for Statistics (Statistics@SLU), a centre for statistical consultation and education at the Swedish University of Agricultural Sciences. In addition I direct the Graduate School in Applied Statistics and Scientific Computing.

Computer-based simulation experiments are one of my main research methods for investigating existing and developing new estimators of statistical characteristics related to forest sciences. For this purpose I often create my own research software by employing higher programming languages in combination with R, SAS or other existing packages.

I disseminate my research results mainly through publications in international peer-reviewed journals. On this websites you find information on some of my research and teaching work. As a service to the international research community I provide computer programs and code that I prepared for teaching and research so that you can explore, test and better understand the methods described in my work. You can find my lab also on [Twitter] and on [GitHub].

Visiting researchers

I am delighted to receive research visitors at all levels from MSc students to senior professors. In my view visiting each other is very important for moving research forward, it is most stimulating and excellent for inspiring a good academic culture. If you have your own travelling and subsistence funding visiting me is straightforward, just let me know when you would like to come. If you need funding for the visit I can try to help with good ideas for a funding proposal. I am also happy to co-supervise and host PhD students from other universities. No fees apply at SLU, all I am asking you is to give a research seminar of 20 to 40 minutes. Short or long visits (sabbaticals) are equally fine by me and it would be great to work on a mutual publication during your visit. Assoc. Professor Mari Myllymäki from Natural Resources Institute Finland (LUKE), a recent visitor, commented on her experience: "I felt most welcome to join the group and greatly enjoyed my stay and the scientific discussions that led to many new research ideas that we are now following up together.”

CRANCOD - A Program for the Analysis and Reconstruction of Spatial Forest Structure

In the last few decades an impressive number of structural indices (also referred to as nearest neigbour summary statistics [NNSS]) have been developed to quantify spatial forest structure. Of particular interest in this regard is the development of a family of individual tree neighbourhood-based indices, which are measures of small-scale variations in tree locations, species and dimensions, developed by Gadow and colleagues at Göttingen University (Germany). Especially when expressed as frequency distributions these indices offer valuable information on spatial woodland structure. Forest structure is closely correlated to and an expression of biodiversity at forest stand level (α diversity). Therefore the structural quantities used in CRANCOD play an important role as surrogate measures of biodiversity (Pommerening, 2002).

The CRANCOD program is a virtual laboratory for analysing and experimenting with nearest neighbour summary statistics and second-order characteristics. CRANCOD has been designed for use with large research plots including full enumerations of trees and in addition offers the opportunity to analyse forest inventory data consisting of multiple sample plots of circular or rectangular shape and varying plot size based on a systematic grid. The program has inbuilt flexibility with the user able to select the number of neighbour trees and to choose between six different methods of edge correction. CRANCOD can, of course, also be used to analyse research and sample plots without spatial information (stand analysis).

The integrated sampling simulator ISIS allows the simulation of sampling with systematic selection of sample plots of varying plot geometry.

Saving individual-tree results in addition to the summary files provides the opportunity to carry out individual-tree based follow-up research. A special visualisation tool allows the user to visually explore nearest neighbour summary statistics. Tree species codings and colours can be flexibly edited externally. A number of language options allow optimal adaptation of the program in different countries.

In a recent study (Pommerening and Uria-Diez, 2017)‎ we found our mingling-size hypothesis confirmed that large trees and trees growing at low local densities frequently have a tendency towards high spatial species mingling. This is often a consequence of disturbances in conjunction with Janzen-Connell effects but also of diversity-oriented forest management. MinglingMand.R‎ is the corresponding R script used in this study. In this graph you see the logistic regression curves of twelve forest stands from different parts of the world describing the probaility of high mingling and their dependence on size.

Author, copyrights and contact: Prof. Arne Pommerening ( or

Software version: 1.4

Licenses: The core version of CRANCOD is a public domain software. However, the program is protected by intellectual property rights and users are expected to acknowledge CRANCOD and its developer when publishing results.

You can find more about the philosophy and objectives of the program along with other details including a downloadable version here.

Selected R scripts and C++ files

Some of the functionality of CRANCOD and additional methods have also been implemented in the R language and in C++. (Computation time can be significantly reduced when combining the advantages of both languages.) When using R scripts/C++ files from the list below for research purposes please acknowledge the author. Instructions on how to integrate R and C++ code are provided here.

Non-spatial analysis

  • Download the script StandAnalysis.R to calculate basal area, trees per hectare, mean diameters, height diameter regression, top height/diameter, diameter distribution and Shannon index along with the Clocaenog 6 sample data. Last updated on 06.07.2012.

Tree growth and modelling

  • Download the script BEM.R and the necessary C++ file BEM.cpp to use the stand development model by Prof. Günter Wenk and colleagues. The model is based on relative growth rates (see my publications). A short model guide will help to get used to the model. Also read Wenk (1994). The file BEMreport.Rnw will allow you to produce Latex/pdf reports of the model outputs. Find an example output here. Last updated on 25.07.2014.
  • Download the R package WenkRegression to model volume, height and diameter development of single stem analysis trees for identifying unusual growth patterns along with sample data from a Sitka spruce tree in Gwydyr forest (North Wales). Unpack the zip file and load the package in R using the package installer. A help file (invoked in R by the command help(WenkRegression)) contains all necessary details. For the modelling theory refer to Wenk (1994) and to my publications in this field. To get started use this example R file. Last updated on 02.12.2015.
  • Download the script AllometricCoefficient.R to calculate the allometric coefficient m quantifying the relationship between diameter and height growth of trees along with sample data from a Sitka spruce tree in Gwydyr forest (North Wales). There is a brief interpretation guide for your convenience. This methodology can also be applied to time series data from forest monitoring with arbitrary survey cycles. Last updated on 20.05.2013.

Diversity indices

Second-order characteristics


  • Download the script Construction.R to simulate a point pattern with a user-defined value of the aggregation index by Clark and Evans (1954) along with a short documentation. Last updated on 18.01.2013.
  • Download a faster version of the script Construction.R along with the necessary C++ code. Last updated on 18.01.2013.
In Pommerening et al. (2018)‎ we found that agreement in tree marking is generally quite low in forest management compared to agreement studies in medicine. We have analysed data from 36 so-called marteloscope experiments from all over Britain and applied the Fleiss kappa characteristic. There was poor to fair agreement in crown thinning experiments (red) and fair to moderate agreement in low thinning experiments (black). Low thinnings are what most British forest managers are used to and crown thinning is a method new to them.

A selection of publications

Pommerening, A., 2002. Approaches to quantifying forest structures. Forestry 75, 305-324. [PdF file]
Pommerening, A. and Murphy, S. T., 2004. A review of the history, definitions and methods of continuous cover forestry with special attention to afforestation and restocking. Forestry. 77, 27–44. [PdF file]
Pommerening, A., 2006. Evaluating structural indices by reversing forest structural analysis. Forest Ecology and Management 224, 266–277. [PdF file]
Pommerening, A. and Stoyan, D., 2006. Edge-correction needs in estimating indices of spatial forest structure. Canadian Journal of Forest Research 36, 1723–1739. [PdF file]
Mason, W. L., Connolly, T., Pommerening, A. and Edwards, C., 2007. Spatial structure of semi-natural and plantation stands of Scots pine (Pinus sylvestris L.) in northern Scotland. Forestry 80, 567-586. [PdF file]
Pommerening, A. and Stoyan, D., 2008. Reconstructing spatial tree point patterns from nearest neighbour summary statistics measured in small subwindows. Canadian Journal of Forest Research 38, 1110–1122. [PdF file]
Pommerening, A., 2008. Analysing & modelling spatial woodland structure. Habilitation thesis BOKU University Vienna. Bangor, 145p. [PdF file]
Davies, O. and Pommerening, A., 2008. The contribution of structural indices to the modelling of Sitka spruce (Picea sitchensis) and birch (Betula spp.) crowns. Forest Ecology and Management 256, 68–77. [PdF file]
Crecente-Campo, F., Pommerening, A. and Rodríguez-Soalleiro, R., 2009. Impacts of thinning on structure, growth and risk of crown fire in a Pinus sylvestris L. plantation in northern Spain. Forest Ecology and Management 257, 1945-1954. [PdF file]
LeMay, V., Pommerening, A. and Marshall, P., 2009. Spatio-temporal structure of multi-storied, multi-aged interior Douglas fir (Pseudotsuga menziesii var glauca) stands. Journal of Ecology 97, 1062-1074. [PdF file]
Motz, K., Sterba, H. and Pommerening, A., 2010. Sampling measures of tree diversity. Forest Ecology and Management 260, 1985–1996. [PdF file]
Murphy, S. T. and Pommerening, A., 2010. Modelling the growth of Sitka spruce (Picea sitchensis (BONG.) CARR.) in Wales using Wenk's model approach. German Journal of Forest Research 181, 35-43. [PdF file]
Pommerening, A., LeMay, V. and Stoyan, D., 2011. Model-based analysis of the influence of ecological processes on forest point pattern formation - A case study. Ecological Modelling 222, 666-678. [PdF file]
Pommerening, A., Gonçalves, A. C. and Rodríguez-Soalleiro, R., 2011. Species mingling and diameter differentiation as second-order characteristics. German Journal of Forest Research 182, 115-129. [PdF file]
Gonçalves, A. C. and Pommerening, A., 2011. Spatial dynamics of cone production in Mediterranean climates: A case study of Pinus pinea L. in Portugal. Forest Ecology and Management 266, 83–93. [PdF file]
Gadow, K. v., Zhang, C. Y., Wehenkel, C., Pommerening, A., Corral-Rivas, J., Korol, M., Myklush, S., Hui, G. Y., Kiviste, A. and Zhao, X. H., 2012. Forest structure and diversity, 29 - 83. In: Pukkala, T. and Gadow, K. v. (Eds.), 2012: Continuous cover forestry. 2nd edition. Managing Forest Ecosystems 23. Springer. Dordrecht, 296p. [PdF file]
Pommerening, A. and Särkkä, A., 2013. What mark variograms tell about spatial plant interactions. Ecological Modelling. 251, 64-72. [PdF file]
Schütz, J.P. and Pommerening, A., 2013. Can Douglas fir (Pseudotsuga menziesii (Mirb.) Franco) sustainably grow in complex forest structures? Forest Ecology and Management. 303, 175-183. [PdF file]
Lilleleht, A., Sims, A. and Pommerening, A., 2014. Spatial forest structure reconstruction as a strategy for mitigating edge-bias in circular monitoring plots. Forest Ecology and Management 316, 47-53. [PdF file]
Hui, G. and Pommerening, A., 2014. Analysing tree species and size diversity patterns in multi-species uneven-aged forests of Northern China. Forest Ecology and Management 316, 125-138. [PdF file]
Pommerening, A. and Maleki, K., 2014. Differences between competition kernels and traditional size-ratio based competition indices used in forest ecology. Forest Ecology and Management 331, 135-143. [PdF file]
Pommerening, A. and Muszta, A., 2015. Methods of modelling relative growth rate. Forest Ecosystems 2, 5. [PdF file]
Abellanas, B., Abellanas, M., Pommerening, A., Lodares, D., Cuadros, S., 2016. A forest simulation approach using weighted Voronoi diagrams. An application to Mediterranean fir Abies pinsapo Boiss stands. Forest Systems 25, e062. [PdF file]
Brzeziecki, B., Pommerening, A., Miścicki, S., Drozdowski, S. and Żybura, H., 2016. A common lack of demographic equilibrium among tree species in Białowieża National Park (NE Poland): evidence from long-term plots. Journal of Vegetation Science 27, 460-469. Doi: 10.1111/jvs.12369. [Link]
Pommerening, A., Brzeziecki, B. and Binkley, D., 2016. Are long-term changes in plant species composition related to asymmetric growth dominance in the pristine Białowieża Forest? Basic and Applied Ecology 17, 408-417. Doi: 10.1016/j.baae.2016.02.002. [Link]
Pommerening, A. and Muszta, A., 2016. Relative plant growth revisited: towards a mathematical standardisation of separate approaches. Ecological Modelling 320, 383-392. [PdF file]
Vítková, L., Ní Dhubháin, Á and Pommerening, A., 2016. Agreement in tree marking: What is the uncertainty of human tree selection in selective forest management? Forest Science 62, 288-296. Doi: 10.5849/forsci.15-133. [PdF file]
Uria-Diez, J. and Pommerening, A., 2017. Crown plasticity in Scots pine (Pinus sylvestris L.) as a strategy of adaptation to competition and environmental factors. Ecological Modelling 356, 117-126. [PdF file]
Pommerening, A. and Uria-Diez, J., 2017. Do large forest trees tend towards high species mingling? Ecological Informatics 42, 139-147. [Link] (R code related to this publication can be found above.)
Stoyan, D., Pommerening, A., Hummel, M. and Kopp-Schneider, A., 2018. Multiple-rater kappas for binary data: Models and interpretation. Biometrical Journal 60, 381-394. [Link]
Kuehne, C., Weiskittel, A., Pommerening, A. and Wagner, R. G., 2018. Evaluation of 10-year temporal and spatial variability in structure and growth across contrasting commercial thinning treatments in spruce-fir forests of northern Maine, USA. Annals of Forest Science 75, 20. [Link]
Pommerening, A., Pallarés Ramos, C., Kędziora, W., Haufe, J. and Stoyan, D., 2018. Rating experiments in forestry: How much agreement is there in tree marking? PLOS ONE 13, e0194747. [Link].
Pommerening, A. and Sánchez Meador, A. J., 2018. Tamm review: Tree interactions between myth and reality. Forest Ecology and Management 428, 164-176. [Link].

Mingling distribution of beech (Fagus sylvatica L., red) and ash (Fraxinus excelsior L., green) in the mixed species woodland Södderich 55B (near Göttingen, Germany). A common situation in woodlands managed for valuable hardwoods: Whilst dominant ash trees mostly have neighbours of different species, the dominated beech trees are arranged in conspecific clusters. The mingling index was calculated with the R script TreeDiversityIndices.R‎.
Gaussian (black) and exponential interaction kernels (red) for trees with a stem diameter at breast height of 60 cm (continuous line) and 20 cm (dashed lines). Interaction kernels are probability density functions describing how ecological processes such as growth, survival and reproduction of an individual depend on the size of and distance to other individuals. A variant of the exponential interaction kernel was applied in the shot-noise model shown at the top of this page (screenshot) and published in Pommerening et al. (2011)‎.

Additional text resources

Visitors of this website who are unfamiliar with the topic of forest structure research and forestry summary characteristics may want to read the following documents:

Feedback and cooperation

The developer of the CRANCOD software and the R scripts is particularly interested in feedback concerning any aspect of the program and of the provided listings to improve their functionality and usefulness in future versions. Please report your feedback to or The author is very open to any kind of cooperation, particularly in terms of joint publications. It would also be possible to maintain CRANCOD and the associated R/C++ files with an international project team in the future.

Graduate School and external teaching

  • R code relating to the NOVA course "Introduction to Scientific Programming and Simulation" at Umeå on 16.06.2017 (Introduction to simulations).
  • Programme and schedule of the Graduate School in Applied Forest Statistics and Scientific Programming (SLU, Umeå) for 2018.
  • New: A free two-days distance learning beginners' course for the software R is available from my GitHub repository.

For SLU courses please visit the SLU website. Do not hesitate to contact Arne Pommerening on or if you are interested in any kind research cooperation or if you simply have a question.

Upcoming conferences and workshops

Cancelled: Hands-on workshop on individual/agent-based modelling, 20-21 August 2018 at Umeå, Sweden.

Research mentoring

It is a simple fact that we all owe a lot to our parents, providence, to our mentors and also to sheer luck. Therefore it is not asking too much to pay back for the good we received by helping others. All academics should help others, particularly students and young researchers, by providing help and support whenever the opportunity arises. As part of this ethos I offer a free, independent research mentoring particularly to those who are at the start of their research careers. The idea is to help with independent and confidential advice on career planning, research paper planning, shaping research visions, promotions, job applications and research questions, conference attendance, proposal writing and academic conflict management outside the normal line management. Meetings can be scheduled during or after working hours. Contact Arne Pommerening on or

I have compiled general research-career and mentoring advice that you find in this document, my Humble guide. Here is a link to Prof. Colin Price's famous guide to academic writing.

A group photo of the NOVA course in Scientific Programming and Simulation at Umeå in June 2017.
A group photo of the course in Applied Spatial Statistics at Umeå in October 2017.


The CRANCOD software, all its libraries and algorithms, the R scripts and C++ files have been prepared with great care. However, the developer can not be made liable in the unlikely case of damage caused to your computer while using CRANCOD, the R scripts/packages and C++ files or of incorrect outputs.