For. Snow Landsc. Res. 78, 1/2: 93–118 (2004) 93

Large-scale biodiversity research in the southern taiga, Northern Mongolia

Michael Mühlenberg1, Hermann Hondong1, Choimaa Dulamsuren1 and Klaus von Gadow2

1 Centre for Nature Conservation, University of Göttingen, von-Siebold-Str. 2, D-37075 Göttingen. [email protected] 2 Institute of Management, University of Göttingen, Büsgenweg 5, D-37077 Göttingen. [email protected]

Abstract The Khentii Mountains of Northern Mongolia, where the Siberian forest belt borders the steppe, represent a unique and greatly untouched ecosystem. Altogether 15 000 km2 of primeval forest and grassland are completely protected by law (Mongolian Ministry for nature and environment 1996). Timber exploitation and use of non-timber forest products are permitted in a buffer zone around the core area. Field research is being conducted since 1996 on an established field station in the western Khentii. The aim of the research is to understand the structure of the spatiotempo- ral heterogeneous forest, its dynamics and the impact of utilisation on the ecology and biodiversi- ty of the forest. To deal with the large map scales, the forest was stratified into different types using a range of vegetation attributes.The types were mapped using a 4 km2 grid within an area of about 150 km2. In order to describe the structure, sample plots were distributed at random within each of the strata. The natural regeneration was studied separately.

Keywords: Mongolia, Siberian taiga, biodiversity, forest structure, forest regeneration, forest man- agement

1 Introduction

This paper describes some results of the Khonin Nuga Large Scale Ecological Research Project in Northern Mongolia. Various research teams representing a variety of scientific disciplines from different countries, operate from the Khonin Nuga research station which has been established in the West-Khentii region by the Centre for Nature Conservation of the University of Goettingen in co-operation with the Faculty of Biology of the National University of Mongolia. The general objective of the Khonin Nuga research project is to understand the spatiotemporal mosaic of the ecosystem, its dynamics and the impact of util- isation on the ecology and biodiversity of the forest. Scientists from many countries, rep- resenting a great variety of disciplines, have been conducting research in the Khonin Nuga project. The project objectives require interdisciplinary and longterm commitment by small groups including one senior scientist and at least one (usually German or Mongolian) junior scientist. Several such teams have been conducting surveys of , insects, small mammals, birds and fish. These surveys and the associated scientific work have been continued by some groups for up to three years. Experts with long-term experience in have visited the Khonin Nuga research station over the years. To become independent from short-term project money and thus ensure long-term survival and continuance the station is run with relatively low cost. 94 Michael Mühlenberg et al.

An important basis for coordinating the different botanical and zoological research activi- ties is a common site classification system. Soil samples were analysed in order to correlate the vegetation types with soil types. Profiles through river valleys were taken by DULAMSUREN (2003) at five different places. These profiles were analysed to understand the natural mosaic of habitats in relation to different forest characteristics. The vegetation classification provided the foundation for a definition of habitat types for the Khentii (MÜHLENBERG et al. 2000). The habitat types allow a more comprehensive description and ecological interpretation than the common community types, by creating compatibility between botanical and zoological studies. For the first few years of investigation it was not our goal to analyse ecosystem functions (e.g. KREBS et al. 2001, reporting about the Kluane Project). Instead, since about half of the flora and fauna of Central is encountered in the region, there is a challenge to simply learn about the ecology of different species by contrasting and comparing the environmental conditions in which they occur. Consequently, the main objectives of the Khonin Nuga research project are:

– To address some fundamental questions of ecology, using reference studies in an environ- ment largely untouched by civilization, including characteristics of natural , the biodiversity of different habitats, landscape heterogeneity, issues of biogeography, phenology, and population biology of selected species. Key hypotheses are: • Low human impact and naturalness are better predictors of species richness than biogeographical factors such as latitude, size of area or regional climate1. • A natural landscape is less fragmented than a cultural landscape and thus facilitates greater mobility and greater niche overlap of the different species. • Habitat selection by species in a natural landscape differs from those found in a land- scape modified by humans, as is the case in Central Europe for example. • The population density of a particular species is often higher in a human-dominated landscape than in a natural landscape2. – To evaluate the conservation value of the region including presence of a near-pristine landscape, occurrence of species which are threatened elsewhere, analysis of communities in primeval habitats as reference for the assessment of anthropogenic impact on species communities in Europe. – To conduct impact assessments in areas where timber has been harvested previously (comparison of faunas between undisturbed and managed taiga, composition of species and dominance of species). – To conduct impact assessments of open-cast gold mining in stream valleys including analysis of water quality, sediment load, and animal communities, up- and downstream of the mining (fish communities and benthos). – To develop ecologically sound natural resource management strategies (forest manage- ment, non-timber-forest-products3).

1 related to the habitat diversity hypothesis presented by GASTON and BLACKBURN (2000) 2 due to the buffer effect (aggregation on few optimal patches), SUTHERLAND (1998) 3 examples are nuts of Pinus sibirica, berries, medicinal plants; deer antlers, … For. Snow Landsc. Res. 78, 1/2 (2004) 95

2 Project area

Mongolia is a country characterized by a rather distinct zonation of vegetation types. The Khentii Mountains of Northern Mongolia, where the Siberian forest belt borders the steppe, represent a unique and greatly untouched ecosystem. The Khentii Mountains have been subdivided in two subprovinces, the Western and Eastern Khentii (SAVIN et al. 1988). In the Western Khentii the and the permafrost starts at a lower altitude than in the East- Khentii (Geokriologicheskie usloviya MNR 1974), resulting in different mountain forest types with different typological structures. The study site is situated in the buffer zone of the Strictly Protected Area of Khan Khentii, in the “forest steppe” which forms the southern extension of the Siberian taiga forest. Natural forest growing on permafrost soils is found on the northern slopes, while the southern slopes receiving greater amounts of solar radiation, are naturally covered with steppe vegetation due to the relatively dry conditions. In this transition zone elements of the boreal forests meet the floristic elements of the Central-Asiatic steppe. That particular region in Mongolia is very much exposed to future development because of its valuable timber, water, and productive pasture resources. A greater number of mesophilic elements are found in the Western Khentii. For example, the Siberian Fir (Abies sibirica) occurs only in the Western Khentii and the Siberian Spruce (Picea obovata) forms forest communities only in the Western Khentii. The herbaceous flora exhibits similar differences in the two subprovinces. According to TSEDENDASCH (1995) and TRETER (1997) the most important factor influencing the formation of closed forest communities in the Khentii is the topo- graphy (slope and exposition). This hypothesis is supported by our observations and our preliminary vegetation assessments. However, important additional factors are precipi- tation, radiation, soil depth and permafrost. The main tree species occurring in the Western Khentii mountains are Siberian Larch (), Siberian ( sibirica), Scots Pine (), Birch species (Betula patyphylla most common, also Betula gmelinii and Betula fruticosa), Siberian Spruce (Picea sibirica). Also found are Poplars (Populus tremula most common, also Populus diversifolia) and Elm (Ulmus pumila). The main vegetation types in the core area of the protected Western Khentii are boreal virgin forest, bog and alpine Tundra.The vegetation classification is given in MÜHLENBERG et al. (2000). Initially, more or less in line with HILBIG and KNAPP (1983), eight vegetation types were identified. However, we found it more useful to develop a hierarchical classification (Table 1). It was not possible to classify forest sites based on the herbaceous flora. Significant differ- ences of the coverage in different habitats (Kruskal-Wallis-Test, P < 0.05) showed only three characteristic plant species for the Pine forest and five species for the Betula–Larix forest, while 108 plant species were indifferent. This result appears to be typical of a heterogenous natural landscape. It would be difficult, if not impossible to develop such a classification using satellite imagery. 96 Michael Mühlenberg et al.

Table 1.The hierarchical classification system of the different habitats around the research station.

Grasslands (with seven subdivisions) Mountain dry steppe (G1a) and meadow steppe (G1b) G1 Herb meadow on the terrace in the river valley G2 Meadow on the river terrace with shrubs of Padus asiatica, Salix sp. G3 Wet meadow with Salix sp. and Betula fusca shrubs G4 Wet grassland dominated by Carex sp. G5 Peat meadow G6 River bank with Carex sp., Equisetum fluviatile and Calamagrostis purpurea G7 Riparian woodland (with five subdivisions) Dense Betula fusca shrub and Salix sp. in the river valley R1 Salix sp. shrubs on the river bank (R2a) and Salix sp. shrub thickets with Pinus sylvestris, R2 Larix sibirica, Padus asiatica (R2b) Open riparian forest with Larix sibirica and Betula platyphylla with shrub layer R3 Picea obovata–riparian forest R4 Populus laurifolia–riparian forest, mixed with Padus asiatica, Crataegus sanguinea, R5 Cornus alba, Salix sp. Mountain forest (with five subdivisions) Larix sibirica–Betula platyphylla forest with different successional stages F1 Mixed forest with dominant (Pinus sylvestris, Abies sibirica, Picea obovata, F2 Larix sibirica, Pinus sibirica, Betula platyphylla) Pinus sylvestris forest and Populus tremula–stands F3 Picea obovata–Abies sibirica forest F4 Pinus sibirica forest (“dark taiga”) F5

3 Assessment methods

3.1 Biodiversity

Biodiversity is investigated in different taxa: higher plants, birds, small mammals, some groups of insects, and fish. Biodiversity research requires assessment of organisms including their occurrence, abundance and distribution. Even for selected taxa it is obviously impossible to assess all the species within such a big area. Therefore one needs to use sampling techniques adapted to the different kinds of species. Appropriate sampling techniques were used for all groups (Table 2). For. Snow Landsc. Res. 78, 1/2 (2004) 97

Table 2. Sampling methods used at the Khonin Nuga research station. 1 Sherman traps; 2 IMS, it is a joint programme of the three German ornithological research stations Wilhelmshaven, Radolfzell, and Hiddensee using standardised methods (Constant Effort Sites) in order to pool the data from different study sites in Eurasia. 3 Angle count method for assessing basal area/ha after BITTERLICH (1948, using the Dendrometer of KRAMER and AKÇA 1995)

Field of research Methods used in Khonin Nuga Botanical surveys Mapping of plant communities in 10 x 10 m2 plots according to a matrix of different ecological factors (variables: slope exposition, canopy closure, soil depth), applying the Braun-Blanquet method. Entomological Standardised catch of butterflies with the same effort (one hour netting along surveys transects in one study plot of 0.5 ha within the chosen habitat, late morning with sunshine and no wind) in 6 different habitats, twice per month for the whole season (May–September). Two habitat types each with 4 replicates in 2002. Small mammal Standardised catch with live traps1 arranged both in a trapping web of 148 traps surveys in three habitats and in a grid of 100 x 100 m each with 121 traps, operating each for a four-days period monthly. In addition 20 m-ditches with two pitfalls each were established in nine habitats especially to collect shrews (Sorex-species). Ornithological Mistnetting in two habitats (108 m each) according to the integrated monitoring surveys of songbird populations2 for the whole season (May till August); mapping along transects in different habitats in spring time (May until July). Hole-nesting birds were surveyed along transects each 1200 m in length in four forest types, the census was conducted twice in each habitat in May and June. Stream ecological Electro-fishing in the Eröö-river and its tributaries, sample stations of 3–30 m surveys length corresponding to three most typical fish habitats were marked out in advance in the river and fished each 3–5 times. Qualitative samples were taken in addition with gill nets by 12–50 mesh sizes and cast nets; fishing with lines occasionally; measurements of the fishes, investigations of the ectoparasites; survey of benthos-community (assessment of relative abundances, with particular interest in stoneflies); standard measurements of physical and chemical parameters including turbidity by photometer. Forest surveys Stratified random sampling: in four forest types 40–60 points were scattered randomly at which variables of forest structure, including the description of dead wood were evaluated. At each point trees were sampled by plot less method with the help of a dendrometer3. The survey of the variables was prepared by working sheets. In addition dead wood was estimated with methods described by KIRBY et al. (1998).

An assortment of sampling methods are available which allow an estimate of the species richness with affordable sampling effort (COLWELL and CODDINGTON 1994). Cumulative species curves may be used to compare different habitats.The method is known as rarefaction (see COLWELL 1997). To deal with the large map scales, the forest was stratified into different types using a range of vegetation attributes. The types were mapped using a 4 km2 (2 x 2 km) grid4 within an area of about 140 km2. The mapped area covers 35 grid cells. The areas of the different vegetation types are shown in Figure 1.

4 following the 1:50 000 universal transverse mercator (UTM) grid 98 Michael Mühlenberg et al.

800

700

600 a Mountain steppe/meadow steppe (G1) b Peat meadow/wet grassland with Salix 500 (G4–6) research area 2 c Betula fusca + Salix shrub (R1, R2) 400 d Riparian forest types (R3, R4, R5) e Pinus sylvestris forest with Larix sibirica, 300 Betula platyphylla, Populus tremula,

ha on 140 km 200 Abies, Picea, Pinus sibirica (F2, F3) f Larix sibirica–Betula platyphylla forest 100 (F1) g Picea–Abies forest (F4) 0 h Pinus sibirica forest (F5) abcdefgh vegetations types

Fig. 1. Areas of different vegetation types in the 35 grid cells of the study site of the Khonin Nuga region. The dominant vegetation is the Larix–Betula forest with its different successional stages.

Each grid cell was sampled using two parallel transects 500 m apart. The vegetation formation was mapped on the spot and located using a global positioning system (GPS).The natural regeneration was separately assessed.

3.2 Sampling and monitoring forests

Forest composition and structure was investigated at large scales, using stratified random sampling.The forest types were classified according to the dominant tree species, resulting in 4 different strata: Larix–Betula forest with different successional stages (F1), Picea–Abies forest (F4), Pinus sibirica forest (F5), and Populus laurifolia riparian forest (R5). In order to describe the structure, 40 to 60 sample plots were distributed at random within each of the strata. In each plot variables of forest structure were assessed, including the dead wood. In total 184 points were sampled. The dynamics of a forest ecosystem is influenced by tree growth which in turn is a reaction to the specific environmental conditions existing on the site. Tree growth data, obtained in a variety of ways, are essential for predicting the consequences of harvesting decisions. The limited availability of research funds and the increasing complexity of the questions that are being addressed by research, necessitate a continuous evaluation of the optimum design of growth trials. Forest management objectives are continually changing. This requires data that permit prediction of forest growth for any set of site conditions and management objec- tives. Three types of growth trials were established. Permanent plots are established for col- lecting data for a particular silvicultural program. The plots are remeasured, usually at regular intervals, until harvesting. Temporary plots, measured only once, provide age-based information about relevant state variables which is used to construct a yield table, again For. Snow Landsc. Res. 78, 1/2 (2004) 99 assuming normal or representative silviculture. Interval plots are remeasured at least once, thus providing an average rate of change in response to a given set of initial conditions. After each remeasurement, a decision was taken whether to abandon the trial or maintain it for another growth interval.

3.2.1 Permanent plots One of the advantages of a database derived from permanent plots is the potential to describe polymorphic growth patterns by evaluating the data of each plot separately and by expressing the parameters of a growth model as a function of specific site variables. In this way, it is possible to develop non-disjoint polymorphic growth models (CLUTTER et al. 1983; KAHN 1994) and disjoint polymorphic site index equations, which depict the site-specific development of certain forest variables over age. A recent example of a polymorphic height model is presented by JANSEN et al. (1996). Many of the existing yield tables are based on permanent plots (SCHOBER 1987; JANSEN et al. 1996; ROJO and MONTERO 1996). A disadvantage of the permanent plot design is the high maintenance cost of the research infrastructure and the long wait for data. The object of the trial is not always achieved, as plots may be destroyed prematurely by wind or fire, or by unauthorized cutting.

3.2.2 Temporary plots Temporary plots may provide a quick solution in a situation were nothing is known about forest growth. They are measured only once, but cover a wide range of ages and growing sites. Thus, the sequence of remeasurements in time is substituted by simultaneous point measurements in space. This method has been used extensively during the 19th century (KRAMER 1988, p. 97; ASSMANN 1953; WENK et al. 1990, p. 116)5. Temporary plots are still being used today for constructing growth models in situations where empirical data are not available (BIBER 1996). For this purpose, increment cores may be taken from a reference tree (usually the last five years are evaluated).To explain variations in diameter growth, it is necessary to evaluate the neighbourhood constellation in the immediate vicinity of the tree. The reference tree should be positioned in the centre of a competition area, the size of which depends on the tree density. Temporary plots are often useful for establishing relationships between variables. The main limitation of temporary plots, when increment cores are not used, is the fact that they cannot provide information about the rate of change of a known state variable, thus preventing the use of some contemporary techniques of growth modelling (GARCÍA 1988).

5 During the 19th century, the “Weiserverfahren” and the “Streifenverfahren” were the most popular methods for obtaining growth information rapidly (KRAMER 1988, p. 97). In the approach known as “Weiserverfahren” the growth of single trees was reconstructed using stem analysis techniques. Another method known as the Streifenverfahren was used to gather data in numerous normally stocked temporary plots of different ages and site qualities for developing yield tables (BAUR 1877). 100 Michael Mühlenberg et al.

3.2.3 Interval plots A compromise may be achieved by using a system of growth trials which maintains the advantages of permanent plots, i.e. obtaining rates of change of known initial states, as well as temporary plots, i.e. broad coverage of initial states and minimum wait for data. Interval plots are measured at least twice, the interval between the measurements being sufficiently long to absorb the short-term effects of abnormal climatic extremes.The interval is a period of undisturbed growth. Measurements should coincide with a thinning operation, to obtain data not only about tree growth, but at the same time about the change of state variables resulting from a silvicultural operation. The thinning effects may be assessed at the initial (t1) or at the final (t2) measurement, or at both occasions. The concept is illustrated in Figure 2. GARCÍA (1988) proposed a multi-dimensional system of differential equations, in which the future development of a forest depends solely on the present state.To be able to develop such a model, it is necessary to have data describing initial states as well as the associated changes of the state variables.

W

W2 * a ∆ W * * a b W 1 *

t1 t2 t ∆t

Fig. 2. Two successive measurements for obtaining the change of a state variable W resulting from a) a thinning and b) natural growth.

3.3 Forest regeneration

Natural regeneration is an important element of forest dynamics.Accordingly, the distribution of the density, height and browse damage is often assessed in forest ecological surveys. The method employed in Khonin Nuga involves 10 m2 circular sample plots (KIRCHHOFF 2003). The sapling representing the sample plot with its height and species is the one nearest to the center of the sample plot. The illustration in Figure 3 shows the representative sapling (Abies with a height of 58 cm) plus two saplings within the circular plot which are used to determine the sapling density, which is equal to 3000 plants per ha. The class frequencies derived from the representative trees represent area proportions (STAUPENDAHL 1997). A disadvantage of this otherwise effective method is the difficulty, due to the small plot size, to capture rare species. For. Snow Landsc. Res. 78, 1/2 (2004) 101

Sample point

Abies/58

Circular plot (r = 1.78 m). “Representative” sapling Number of saplings in plot nearest to sample point; defines sapling density. used to determine species and height.

Fig. 3. Schematic representation of a regeneration sample plot showing the tree species and sapling height in cm. The representative sapling in this example is Abies, with a height of 58 cm; the sapling density is 3000 plants/ha.

4 Preliminary results

4.1 Species diversity

More than 1150 plant species characteristic of the steppe ecosystem, 253 bird and more than 50 ungulate species were identified in the protected area of the Khentii. Prominent large mammals are the Maral (Cervus elaphus maral), Moose (Alces alces), Siberian Roe Deer (Capreolus pygargus), Musk Deer (Moschus moschiferus), Wild Boar (Sus scrofa), Brown Bear (Ursos arctos), Wolf (Canis lupus), Lynx (Lynx lynx), Wolverine (Gulo gulo), (Martes zibellina) (READING et al. 1994). Except for the Musk Deer all other mentioned species may occur in Europe.Thus the area under study can serve in some way as a reference area representing natural conditions in Europe. Considering the entire fauna about half of the species encountered in the study area are palearctic. Table 3 presents an overview of the biodiversity found in Khonin Nuga and a comparison with findings from other areas in Europe. 50 percent of the butterfly species and 51 percent of the bird species found in Khonin Nuga are palaearctic and occur also in Central Europe (Fig. 4). The analysis of cumulative species curves for butterflies shows that, in order to evaluate the species diversity in different habitats, it is necessary to capture at least 2000 individuals. Rarefaction curves are suitable for describing differences between habitats. They com- pare species numbers at the same sample size, in our case with the same amount of captured individuals. Figure 5 shows the rarefaction curve of the butterfly community in two habitats. 102 Michael Mühlenberg et al.

Table 3. Comparison of known species numbers in different regions. 1 KARSHOLT and RAZOWSKI 1996, Lepidoptera of Europe; 2 JONSSON 1992, Vögel Europas; 3 HAGE- MEIJER and BLAIR 1997, EBCC Atlas of European Breeding Birds; 4 BfN Rote Liste Deutschlands, 1998; 5 bird species number in Germany inclusive guests: 515 species; 6 DENNIS 1992; 7 GREENWOOD et al. 1993; 8 MÖCHBAYAR 1999; 9 Redkie Zivotnye Mongolii (pozvonocnye), Moskva 1996; 10 DAWAA et al. 1994: Kommentierte Checkliste der Vögel und Säuger der Mongolei; 11 ULZIJCHUTAG 1989; 12 own data 1998–2002. BV = Brutvögel = Breeding bird species.

Region Number of known Number of known Known species Total land butterfly species bird species number of higher area (km2) plants Europe 4681 4692, 500 BV3 12 500 10 531 000 Germany 1854 2883, 260 BV5 2691 357 042 Great Britain 626 215 BV7 1494 241 752 Mongolia (207)8 415 (322 BV)9, 440 282311 1 565 000 (360 BV)10 Khonin Nuga12 146 162, 123BV 553 140

600 Species cumulative curve of butterflies in West Khentei

500 140 s 400 120 pecie s s 100 300 pecie

s 80 200 60 Number of 100 40 Number of

0 20 a i ss olia 0 hina g C weden

Ru 01000 2000 3000 4000 5000 6000 7000 Britain ermany Europe Finland S G Mon Number of individuals

Fig. 4. Species richness and distribution of butterflies. Left: the black part of the columns indicates the number of species shared with Mongolia (only the butterfly assemblage of Khonin Nuga is presented for Mongolia); right: cumulative curve for butterfly species in West Khentei, species pooled from the catch of 2000 and 2001. Broken lines indicate the 95%-confidence limit. For. Snow Landsc. Res. 78, 1/2 (2004) 103

120 Herb Meadow

100

80 Mountain Dry Steppe

60

40 Expected number of species 20

0 0 500 1000 1500 2000

Fig. 5. Rarefaction curves showing the species richness of herb meadow (habitat G2) and mountain dry steppe (habitat G1). Data are pooled from 4 replicates of each habitat in 2002.

Single Linkage Cluster Distance = 1 – Morisita Horn Cluster Diagramm Butterflies 2002 (n = 3027)

HM1

MDS2

MDS1

MDS3

MDS4

HM2

HM3

HM4

0.10 0.15 0.20 0.25 0.30 0.35 Distance

Fig. 6. Cluster diagram of the dissimilarity indices (1 – Morisita Horn) for the butterfly assemblages of the 4 herb meadow plots (HM) and of the 4 mountain dry steppe plots (MDS).

Altogether 2114 individuals in 111 species were found in the Herb Meadow (G2) type and 913 individuals in 95 species in the Mountain Dry Steppe (G1), in the sample of the year 2002. Habitat G2 has a higher species richness than habitat G1. The butterfly population was used to test the faunistic similarity between different habitats. As an example the butterfly community of two habitats which appear particularly different is compared: a moist herb meadow with tall grass (G2, G3) and a mountain dry steppe with short grass (G1). Twice per month, between May and August 2002, the butterfly 104 Michael Mühlenberg et al. assemblages were assessed on four sampling areas on each vegetation type using the standard sampling method. The result is presented using a cluster diagram (Fig. 6). The cluster diagram shows the dissimilarity indices (1 – Morisita Horn) for the butterfly assemblages of the four herb meadow plots (HM) and of the four mountain dry steppe plots (MDS). The ANOVA test reveals no overall differences between the plots (R = 4.07, p = 0.089). These preliminary results confirm a high similarity between different habitats in the studied natural landscape. The similarity indices of the small mammal communities between different habitats are surprisingly high as well, indicating the same large overlapping of species in the habitats. These findings support the hypothesis of greater mobility and greater niche overlap in a natural landscape (refer to the Chapter 1).

4.2 Forest spatial structure and diversity

The old-growth forest (Pinus sibirica – taiga, F5) exceeds all other forest types in basal area of living and dead wood, but for cavity-nesting birds the holes in Betula trees are most important. Different forest types (e.g. successional stages of Larix–Betula and Picea–Pinus) can be grouped together due to fire disturbances considering the tree species composition. Riparian woodland (R2, R3, R5) sustains the highest biodiversity but is most restricted in area.

4.2.1 Habitat trees The types of damage investigated were fire, wind and rot.Wind caused breakages, rot hollow trees. Most of fire damages and most of the hollow trees are found in the Betula–Larix forest stands (F1). Most of the tree hollows are provided by Betula platyphylla trees (BAI et al. 2003). Betula is therefore a key species for cavity nesting birds. The Pinus sibirica forest exceeds the other forest types, both in total basal area and number of big-diameter trees. Clustering these samples according to the basal area of the tree species confirms the stratification by methods of vegetation analysis (Fig. 7a, cluster 1). Pinus sibirica forest is most clearly separated from other forest types. The mixture of Betula–Larix forest stands with Picea–Pinus and Populus riparian forest at right hand in the cluster documents the high influence of fire, in all stands Betula platyphylla as a pioneer tree is represented with a rather high basal area. Some plots of Picea–Abies–Pinus forest (green colour) are grouped together like a low disturbed conifer forest (in the cluster 1 right of the Pinus sibirica block). If clustering is done with variables of structure, the picture changes (Fig. 7b, cluster 2): the sample plots are now not grouped in the vegetation formations (the classified 4 forest types). The mixture reflects more the high dynamics in the natural landscape. A relevant feature of interest to biologists is the occurrence of “wildlife trees”, rep- resenting particularly big trees, dead or broken trees, trees with hollows and fire damage. Clustering with these variables leads to the cluster 3 (Fig. 7c). Consider the significant dif- ferent picture of the cluster 3 in comparison with cluster 1. One conclusion is that mapping of vegetation or interpretation of satellite photos according to vegetation classification methods (e.g. dominant tree species) does not necessarily delineate important stands for conservation purposes. It may be concluded that for evaluating conservation values, terrestrial assessment is also needed. For. Snow Landsc. Res. 78, 1/2 (2004) 105

Fig. 7 a.Three clusters created with different sets of variables.The colour assigns the sample point to one of the four stratified forest types. Cluster analysis grouping of 184 point samples according to the basal area of the tree species. Orange = samples in Betula–Larix forests, green = samples in Picea–Abies– Pinus forest, brown = samples of Pinus sibirica forest, blue = samples in riparian forest with Salix and Populus laurifolia.

Fig. 7 b.Three clusters created with different sets of variables.The colour assigns the sample point to one of the four stratified forest types. Cluster analysis grouping of 184 point samples according to the diam- eter classes of the trees. Orange = samples in Betula–Larix forests, green = samples in Picea–Abies– Pinus forest, brown = samples of Pinus sibirica forest, blue = samples in riparian forest with Salix and Populus laurifolia. 106 Michael Mühlenberg et al.

Fig. 7 c.Three clusters created with different sets of variables.The colour assigns the sample point to one of the four stratified forest types. Cluster analysis grouping of 184 point samples according to relevant variables of conservation. The variables are dbh > 50cm, dead-, broken-, hollow-, fire-tree. Orange = samples in Betula–Larix forests, green = samples in Picea–Abies–Pinus forest, brown = samples of Pinus sibirica forest, blue = samples in riparian forest with Salix and Populus laurifolia.

4.2.2 Forest spatial structure in Sangstai Forest The Sangstai Forest, representing of old-growth forest, was studied in greater detail. Old- growth forests are found in places with very low fire frequency. In the Khentii region fire did not affect the remote mountain ridges with wet mossy ground vegetation and shallow soil layers. Another region not affected by fire is situated in the river valley between water bodies where riparian woodland is found. The “structure” of a forest may be defined by the spatial distribution of the tree positions, by the spatial mingling of the different tree species and by the spatial arrangement of the tree dimensions. The spatial structure is one of the characteristic attributes of a forest. The problem is to characterize and describe forests with different spatial characteristics more accurately, using affordable assessment techniques. The Sangstai plot in the Khentii may be used to demonstrate an approach to describe the spatial forest structure and diversity (Fig. 8). L- and Pair correlation functions are useful for describing forest structures, but they require datasets with known tree positions (STOYAN and STOYAN 1992; PRETZSCH 2001; POMMERENING 2002). Such data are hardly ever available in practice and this precludes their use. Aggregate indices, such as the spatial index proposed by CLARK and EVANS (1954), can provide a first general impression of the structure of a particular forest, but they cannot be used to describe the great variety of spatial arrangements (ZENNER and HIBBS 2000).This problem is especially serious in very irregular forests where small-scale structural characteristics are highly variable (ALBERT 1999). For. Snow Landsc. Res. 78, 1/2 (2004) 107

For this reason, three types of neighbourhood-based parameters are used, which are known as Contagion, Mingling and Differentiation. The parameters can be used to provide a comprehensive description of the spatial structure of a forest. Assessment and description may be tree-based or point-based. In the tree-based approach a sample tree closest to a sample point is chosen as reference tree and the attributes of its immediate neighbours (size, species) and the regularity of their positions are related to the reference tree. In the point-based approach, the structural attributes of a neighbourhood group of trees (variation of tree species and sizes; regularity of tree positions) is assessed at each sample point.

Species Trees per ha Basal area per ha Abies sibirica 280 4.22 Larix sibirica 8 1.64 Picea obovata 40 0.64 Pinus sibirica 212 34.69 Total 540 41.20

Fig. 8. Sangstai plot with buffer showing tree positions (left) and the corresponding plot data with trees (N/ha) and basal areas (G/ha) per hectare listed for the four species occurring in the plot (right). Red = samples in Betula–Larix forests, green = samples in Picea–Abies forests, spotted = samples of Pinus sibirica forests, blue = samples in riparian forests with Salix and Populus Laurifolia.

Contagion The variable contagion Wi describes the degree of regularity of the spatial distribution of the 6 four trees nearest to a reference tree i . Wi is based on the classification of the angles αj between these four neighbours. A reference quantity is the standard angle α0, which is expected in a regular point distribution. The binary random variable vj is determined by comparing each αj with the standard angle α0. The Contagion is then defined as the pro- portion of angles αj between the four neighbouring trees which are smaller than the standard angle α0:

4 1 ⎧1, αj < α0 W = ∑v withv = ⎨ and 0 ≤W ≤1 (1) i j j 0, otherwise i 4 j=1 ⎩

6 For details refer to GADOW et al. (1998). Four neighbours have proved to be most suitable based on practical considerations in connection with the field assessment methods (ALBERT 1999; HUI and HU 2001). 108 Michael Mühlenberg et al.

Wi = 0 indicates that the trees in the vicinity of the reference tree are positioned in a regular manner, whereas Wi = 1 points to an irregular or clumped distribution. With four neighbours, there are five possible values that Wi can assume. The estimator for the Contagion of a given forest is W the arithmetic mean of all Wi-values. Although the Contagion mean value W is quite informative for characterizing a point distribution, it is often advisable to study the distribution of the Wi -values which reveals the structural variability in a given forest (Table 4).

Table 4. Distribution of the variable “contagion” which describes the degree of regularity of the spatial distribution of the four trees nearest to a reference tree i.The spatial distribution is random with a small proportion of very clumped neighborhoods.

W All species Abies sibirica Larix sibirica Picea obovata Pinus sibirica 0.00 0.00 0.00 0.00 0.00 0.00 0.25 0.21 0.16 0.50 0.50 0.23 0.50 0.54 0.57 0.50 0.40 0.53 0.75 0.22 0.24 0.00 0.10 0.23 1.00 0.02 0.01 0.00 0.00 0.02

Based on the work by HUI and GADOW (2002), the spatial distribution may be characterized as random, although about two percent of the trees are situated in a neighbourhood with a very clumped distribution (W = 1.0). Neighbourhoods with a very regular distribution (W = 0) are not encountered.

Species mingling Species diversity has become a very important aspect of forest management and conservation and a number of parameters are available to describe it.An example is the Shannon-Weaver index which has been used in ecological applications by PIELOU (1977, p. 293). We propose to evaluate the species diversity in the vicinity of a reference tree and define mingling as the proportion of the n nearest neighbours that do not belong to the same species as the reference tree (GADOW and FÜLDNER 1992), specifically:

4 = 1 M i ∑v j (2) 4 j=1 ⎧1, neighbour j belongs to the same species as reference tree i with v = ⎨ j ⎩0, otherwise

and 0 ≤ ≤ 1 M i For. Snow Landsc. Res. 78, 1/2 (2004) 109

With four neighbours, the mingling attribute Mi can assume five values. Table 5 presents the mingling distributions for all trees and for each of the four tree species in the Sangstai plot.

Table 5. Distribution of the variable “mingling” which describes the degree of regularity of the spatial distribution of the four trees nearest to a reference tree i.

M All species Abies sibirica Larix sibirica Picea obovata Pinus sibirica 0.00 0.13 0.23 0.00 0.00 0.04 0.25 0.22 0.23 0.00 0.00 0.26 0.50 0.22 0.26 0.00 0.00 0.23 0.75 0.26 0.21 0.00 0.40 0.30 1.00 0.16 0.07 1.00 0.60 0.17

The mingling distribution for all species shows a variety of mingling constellations. 13% of the trees, for example, occur in pure groups and 22% in groups where half of the trees are of the same species. As expected, the rare species (Larix and Picea) have the highest mingling values.The most frequent species Abies and Pinus occur in all the different mingling constel- lations. Abies sibirica is frequently found in pure groups.

Size differentiation and dominance The tree attribute “dominance” of neighbours was proposed by HUI et al. (1998) to relate the relative dominance of a given tree species to the immediate neighbourhood. We define dominance as the proportion of the n nearest neighbours of a given reference tree which are smaller than the reference tree, which is calculated in the same way as the previous tree- based structural parameters:

4 = 1 U i ∑ v j (3) 4 j=1 ⎧1, neighbour j is smaller than reference tree i with v = ⎨ j ⎩0, otherwise

and 0 ≤ ≤ 1 U i

With four neighbours, Ui can assume five values. Low U-values indicate dominance. The dominance criterion is useful if we wish to describe the relative dominance of a particular tree species. The highest value of 1 means that the tree is the smallest one in its immediate neighbourhood. Figure 9 shows the results for the two most common species, Pinus sibirica and Abies sibirica. The high dominance values of Pinus sibirica can be expected as this species is represented with a basal area of almost 35 m2/ha and less trees per ha than Abies sibirica which is rep- resented by a basal area of only 4.22 m2/ha. Pinus sibirica and Larix sibirica are mostly dominant while Abies sibirica and Picea are more subdominant or suppressed. 110 Michael Mühlenberg et al.

0.5

0.4

0.3 P 0.2

0.1

0.0 0.00 0.25 0.50 0.75 1.00 U Fig. 9. Distribution of dominance values. Pinus sibirica (shaded columns) occurs mostly as a dominant tree while Abies sibirica (white columns) is mostly subdominant.

4.3 Forest regeneration

Sustainable development of a natural or managed forest ecosystem depends on the ability of the system to regenerate itself. The recruitment potential is a key factor in the Southern Taiga forests which are regularly affected by sometimes very destructive and large scale wildfires. To evaluate this potential, KIRCHHOFF (2003) made an assessment of the natural regeneration in three different forest environments. Figure 10 presents an impression of three assessment sites and a graph of the terrain features and occurrence of the different tree species.The corresponding distributions of regeneration density classes for the different tree species are also shown.

Sample sites were chosen with three specific questions in mind: 1) What is the capability of the forest to recolonize burnt areas? 2) What is the regeneration potential in the managed forests, which were heavily exploited towards the end of the 20th century and where mostly Larix sibirica was cut? 3) What is the regeneration potential in the virgin dark taiga forests, a very sensitive eco- system which is dominated by Pinus sibirica and Abies sibirica?

This first assessment done by KIRCHHOFF (2003) suggests that forest regeneration is not endangered in any of these three problem sites. For. Snow Landsc. Res. 78, 1/2 (2004) 111

East West

Pinus sylvestris

Larix sibirica

Populus tremula

Betula platyphylla

60 50 40 30

Fig. 10 a. Burnt forest Sharlang percent 20 Altitude 1000–1100 m, slope 20–30°; regeneration in strips parallel with slope; two regeneration age 10 classes are found. 0 no regen. Betula Populus Pinus Larix platyphylla sylvestris sibirica

no regeneration 1000–5000 per ha 5000–10000 per ha 10000–50000 per ha

South North

Betula platyphylla

Larix sibirica

50

40

30

percent 20 Fig. 10 b. Managed forest Hausberg Altitude 900 m, moderate slope; regeneration in 10 clumps; regeneration of Abies sibirica. 0 no regen. Betula Abies Pinus Larix platyphylla sibirica sylvestris sibirica

no regeneration 1000–5000 per ha 5000–10000 per ha 10000–50000 per ha 112 Michael Mühlenberg et al.

West East

Pinus sibirica Abies sibirica Picea obovata

60 50 40 30

Fig. 10 c. Virgin taiga forest Sangstai percent 20 Altitude 1500 m, flat; dense forest with few gaps; 10 regeneration in clumps. 0 no regen. Abies Pinus sibirica sibirica

no regeneration 1000–5000 per ha 5000–10000 per ha 10000–50000 per ha

Fig. 10 a–c. Three examples of areas where natural regeneration was assessed (KIRCHHOFF 2003).

5 Discussion

Biological diversity describes the variety of life at different levels of biological organisation (SPELLERBERG and SAWYER 1999). Inventorying biodiversity involves the surveying, sorting, cataloguing, quantifying and mapping of entities such as genes, individuals, populations, species, habitats, biotopes, ecosystems and landscapes or their components, and the synthesis of the resulting information for the analysis of processes (HEYWOOD 1995)7. Based on practical considerations, assessment and analysis are usually concentrated on the species level. We selected in our study sites some taxa from which we have some knowledge for comparison and the experts available to work in the field. It is not our aim, to approach an “All Taxa Biodiversity Inventory” (JANZEN and HALLWACHS 1994). Half of the species found in our project area are palearctic, i.e. the area under study serves in some way as a reference area with natural conditions for Europe, e.g. the still existing coexistence of all big mammals (top carnivores and big herbivores) in an unchanged land- scape. Our study of biodiversity aims to find answers to the following questions:

7 Global Biodiversity Assessment, (GBA; HEYWOOD 1995) For. Snow Landsc. Res. 78, 1/2 (2004) 113

1) What are the relationships between certain attributes of forest structure and the variables describing biodiversity? Can forest structure attributes be used to predict biodiversity? 2) What is the biogeographic significance of the Khentii region in the international context? 3) What is the conservation value of the Khentii region in the national and international context? 4) Is it possible to use biodiversity indicators to describe human impacts in the region?

Clearly, the answers to these questions may turn out to be dependent on the spatial scale of our work. Until now, research in Khonin Nuga has concentrated on the vascular plants, birds, small mammals, butterflies, grasshoppers, fish and stoneflies. To be able to deal with the questions (1) and (2), it was necessary to develop a classification and description of the vegetation and habitat types. Concerning question (1), relationships have been established between the structural variety of the vegetation and the number of species richness (e.g. KARR and ROTH 1971; WILLSON 1974; MÜHLENBERG 1980; ARNOLD 1983; NILSSON et al. 1988; JEDRZEJEWSKA et al. 1994; KUJAWA 1997; SULLIVAN et al. 2001; LOHR et al. 2002). Animal groups that exploit the environment in three dimensions are most sensitive to plant community structure, which has been shown in the classical study by MACARTHUR and MACARTHUR (1961); MACARTHUR (1964) who found a correlation between foliage height diversity and bird species richness. In our area the Larix–Betula forest (F1) as the most extended forest harbours the richest breeding bird community (41 out of 109 species), followed from the riparian woodland (R5) with 20 species out of 109. 18 breeding bird species are recorded in the Pinus sibirica forest (F5) (WICHMANN 2001). In the successional series of F1 the bird species richness increases from burned area to young forest to old forest, indicating an increase in biodiversity with increasing structure (WICHMANN 2001). BOURSKI (1996) confirmed a highest richness in the breeding bird assemblage for the flood-plain (corresponding to our riparian woodlands), decreasing to the taiga forest and last to the burned areas. The three modes of clustering our sample points of forest structure show us, that basal area of tree species or simple variables like diameter are not useful to predict biodiversity. For that approach a set of specific variables has to be measured on a large scale. The Khentii mountains are part of the “Transbaikal region” (for bird studies see KOZLOVA 1930; GLADKOV and SELIVONIN 1963; BOLD 1984; VASILCHENKO 1987). The species richness of birds and trees is known to be higher in the eastern Palearctic than in the western part due to different histories in the two biogeographical regions (MÖNKKÖNEN 1994; MÖNKKÖNEN and VIRO 1997). Our botanical survey shows that the forested northern slopes are home to more western Eurosiberian and Uralosiberian flora elements whereas on the steppe of the southern slopes East-Asian elements from the Mandshurian-Altaian- Dahurian area are dominating (DULAMSUREN and MÜHLENBERG 2003). The national conservation value of the Khentii region (question 3) is highlighted by map- ping of the plant species out of the Red Data Book of Mongolia (MÜHLENBERG et al. 2000). The international conservation value is documented by the presence of many palearctic species of which the populations in Europe are threatened (MÜHLENBERG et al. 2000; WOYCIECHOWSKI et al. 2001). About three quarters of the palearctic butterflies, half of the palearctic species of birds and a third of the palearctic plant species found in Khonin Nuga have some threat status in Central Europe (BfN 1996, 1998). The overall conservation value of the region exists because of the huge natural landscape itself (>> 20 000 km2) which is not yet altered by humans. The unmanaged forests provide fallen timber and a great amount of woody debris what is generally seen as being of conservation value (JONSELL et al. 1998; KLAUSNITZER 1999; IRMLER et al. 1996; JONSSON and KRUYS 2001; MACNALLY et al. 2002; GÖTMARK and THORELL 2003). 114 Michael Mühlenberg et al.

Question (4) relating to the use of biodiversity indicators to describe human impacts is still under investigation. The importance and potential effects of many proposals for forest habitats has to be ascertained for the Khentii region (e.g. LANDRES et al. 1988; PEARSON and CASSOLA 1992; WEAVER 1995; NILSSON et al. 1995; STORK et al. 1997; DUFRÊNE and LEGENDRE 1997; NIEMELÄ 1997; JONSSON and JONSELL 1999; LINDENMAYER et al. 1999; KERR et al. 2000; MIKUSINSKI et al. 2001; TAYLOR and DORAN 2001; RAINIO and NIEMELÄ 2003). This goes beyond ecological considerations and implies political ones as well. Only large samples can help to understand the correlation between forest structure and biodiver- sity. One of the future challenges is the development of a sustainable management system for the forest resources which does not yet exist in Mongolia. Concerning mature forest habitat and the maintenance of coarse woody debris for biodiversity, European guidelines for sustainable management cannot easily be adapted. Alternative guidelines need to be developed and evaluated for this very unique ecosystem.

Acknowledgements We acknowledge the generous support of the Mongolian project partners, especially the Faculty of Biology of the National University of Mongolia and the German Technical Service (GTZ) in Mongolia. The field data were collected by students from Mongolia, Taiwan and Germany; regen- eration data were collected by B. Kirchhoff; tree structure data by A. Gradel. H. Heydecke helped with data processing using the software developed by Chen BoWang. We are grateful for useful comments received from two anonymous referees.

6 References

ALBERT, M., 1999: Analyse der eingriffsbedingten Strukturveränderung und Durchforstungs- modellierung in Mischbeständen. PhD Diss., Faculty of Forest Sciences, Universität Göttingen, Germany. Hainholz Verlag. 63–68. ARNOLD, G.W., 1983: The influence of ditch and hedgerow structure, length of hedgerows, and area of woodland and garden on bird numbers on farmland. J.Appl. Ecol. 20: 731–750. ASSMANN, E., 1953: Zur Bonitierung süddeutscher Fichtenbestände. Allg. Forst Z. Waldwirtsch. Umweltvorsorge 10: 61–64. BAI, M.-L.; WICHMANN,F.;MÜHLENBERG, M., 2003: Hole abundance and nest-site preference of hole-nesting birds in a primeval boreal forest of Mongolia. Acta Ornithol. 38: 95–102. BAUR, F., 1877: Die Fichte in Bezug auf Ertrag, Zuwachs und Form. Berlin. 103 pp. BIBER,P.,1996: Konstruktion eines einzelbaumorientierten Wachstumssimulators für Fichten- Buchen-Mischbestände im Solling. Diss. Forstwissenschaftliche Fakultät der Ludwig- Maximilians-Universität München. 239 pp. BITTERLICH, W., 1948: Die Winkelzählprobe. Allg. Forstztg. [Wien]. 4–5. BOLD, A., 1984: Die Vögel des Chentej und ihre praktische Bedeutung. Wissenschaftliche Beiträge Universität Halle-Wittenberg 4: 132–151. BOURSKI, O., 1996: Bird population dynamics in relation to habitat quality. In: SETTELE,J.; MARGULES, C.R.; POSCHLOD, P.;HENLE, K. (eds) Species Survival in Fragmented Landscapes, Kluwer Academic Publishers. 52–60. BfN (Bundesamt für Naturschutz), 1996: Rote Liste gefährdeter Pflanzen Deutschlands. Bonn- Bad Godesberg, Schriftenreihe für Vegetationskunde 28: 744 pp. BfN (Bundesamt für Naturschutz), 1998: Rote Liste gefährdeter Tiere Deutschlands. Bonn-Bad Godesberg, Schriftenreihe für Landschaftspflege und Naturschutz 55: 434 pp. For. Snow Landsc. Res. 78, 1/2 (2004) 115

CLARK, P.J.; EVANS, F.C., 1954: Distance to nearest neighbor as a measure of spatial relationships in populations. Ecology 35, 4: 445–453. CLUTTER, J.L.; FORTSON, J.C.; PIENAAR, L.V.; BRISTER, G.H.; BAILEY, R.L., 1983: Timber man- agement – a quantitative approach. John Wiley: 333 pp. COLWELL, R.K., 1997: EstimateS 5. Statistical Estimation of Species Richness and Shared Species from Samples. Version 5.0.1. [Published online] Available from internet . Department of Ecology and Evolutionary Biology, University of Conneticut Storrs, CT 06269-3042, USA. COLWELL, R.K.; CODDINGTON, J.A., 1994: Estimating terrestrial biodiversity through extra- polation. Philos. Trans. R. Soc. Lond. (Series B) 345: 101–118. DAWAA, N.; BUSCHING, W.-D.; SUMIYAA, D.; BOLD, A.; SAMJAA, R., 1994: Kommentierte Checkliste der Vögel und Säuger der Mongolei (Köthen), Bd. 1, Vögel. DENNIS, R.L.H., 1992: The Ecology of Butterflies in Britain. Oxford University Press. 354 p. DUFRÊNE, M.; LEGENDRE, P., 1997: Species assemblages and individual species: the need for a flexible assymetrical approach. Ecol. Monogr. 67, 3: 345–366. DULAMSUREN, C.; MÜHLENBERG, M., 2003: Additions to the flora of the Khentii, Mongolia. Willdenowia 33: 149–158. DULAMSUREN, C., 2003: Floristische Diversität, Vegetation und Standortbedingungen in der Gebirgstaiga West-Khentii, Mongolia. PhD thesis Georg-August-Universität, Göttingen. DULAMSUREN, C.; SOLONGO,B;MÜHLENBERG, M., submitted: Comment on the red data book of endangered plant species of Mongolia. J. Biol. Sci. Mongolia 1. GADOW,K.V.; FÜLDNER, K., 1992: Bestandesbeschreibung in der Forsteinrichtung. Proc. of the working group Forsteinrichtung. Klieken near Dessau 15.10.92. GADOW,K.V.; HUI, G.Y.; ALBERT, M., 1998: Das Winkelmass – ein Strukturparameter zur Be- schreibung der Individualverteilung in Waldbeständen. Cent.bl. gesamte Forstwes. 115, 1: 1–9. GARCÍA, O., 1988: Experience with an advanced growth modelling methodology. In: EK, A.R.; SHIFLEY, S.R.; BURKE, T.E. (eds) Forest growth modelling and prediction. USDA For. Serv. Gen. Tech. Rep. NC-120: 668–675. GASTON, K.J.; BLACKBURN, T.M., 2000: Pattern and Process in Macroecology. Oxford, Blackwell Science. 376 pp. Geokriologicheskie usloviya MVR, 1974, Moskau. GLADKOV, N.A.; SELIVONIN, J.A., 1963: Die Ornithogeographischen Unterschiede des Östlichen und Westlichen Teiles der Taigazone. Proc. XIII Intern. Ornithol. Congr. 1109–1116. GÖTMARK, F.; M. THORELL, 2003: Size of nature reserves: densities of large trees and dead wood indicate high value of small conservation forests in southern Sweden. Biodiver. Conserv. 12: 1271–1285. GREENWOOD, J.J.D.; BAILLIE, S.R.; CRICK, H.Q.P.; MARCHANT, J.H.; PEACH W.F., 1993: Integrating population monitoring: detecting the effect of diverse changes. In: FURNESS, R.W.; GREENWOOD, J.J.D. (eds) Birds as Monitors of Environmental Changes. London, Chapman and Hall. 293–315. HAGEMEIJER, W.J.M.; BLAIR, M.J., 1997: EBCC Atlas of European Breeding Birds. London, T & A. D. Poyser. 903 pp. HEYWOOD, V.H.H. (ed), 1995: Global biodiversity assessment. Cambridge University Press. 1140 pp. HILBIG, W.; KNAPP, H.D., 1983: Vegetationsmosaik und Florenelemente an der Waldsteppen Grenze im Chentej Gebirge (Mongolei). Flora 174: 1–89. HUI, G.H.; HU, Y.B., 2001: Measuring species spatial segregation in mixed forest. For. Res. 14, 1: 23–27. HUI, G.Y.U.; GADOW,K.V., 2002: Das Winkelmass – Herleitung des optimalen Standardwinkels. Allg. Forst- Jagdztg. 173: 10. HUI, G.Y.;ALBERT, M.; GADOW,K.V., 1998: Das Umgebungsmass als Parameter zur Nachbildung von Bestandsstrukturen. Forstwiss. Cent.bl. 117: 1: 258–266. IRMLER, U.; HELLER, K.; WARNING, J., 1996: Age and tree species as factors influencing the popu- lation of insects living in dead wood (Coleoptera, Diptera: Sciaridae, Mycetophilidae). Pedobiologia 40: 134–148. 116 Michael Mühlenberg et al.

JANSEN, J.J.; SEVENSTER,J.U.; FABER, P.J., 1996: Opbrengsttabellen voor belangrijke boomsoorten in Nederland. Hinkeloord Report Nr. 17, Landbouwuniversiteit, Wageningen. 202 pp. JANZEN, D.H.; HALLWACHS, W., 1994: All Taxa Biodiversity Inventory (ATBI) of Terrestrial Systems. A generic protocol for preparing wildland biodiversity for non-damaging use. Draft Report of National Science Foundation Workshop, April 1993, Philadelphia Pa. 16–18. JEDRZEJEWSKA, B.; OKARMA, H.; JEDRZEJEWSKI, W.; MILKOSKI, L., 1994: Effects of exploitation and protection on forest structure, ungulate density and wolf predation in Bialowieza primeval forest, Poland. J.Appl. Ecol. 31: 664–676. JONSELL, M.; WESHIN, J.; EHNSTRÖM, B., 1998: Substrate requirements of red-listed saproxylic invertebrates in Sweden. Biodiv. Cons. 7, 6: 749–764. JONSSON, B.G.; KRUYS, N. (eds) 2001: Ecology of woody debris in boreal forests. Lund, Sweden, Ecol. Bull. 49: 281 pp. JONSSON, B.G.; JONSELL, M., 1999: Exploring potential biodiversity indicators in boreal forests. Biodiv. Cons. 8: 1417–1433. JONSSON, L., 1992: Die Vögel Europas und des Mittelmeerraumes. Stuttgart, Franckh Kosmos Verlag. 559 pp. KAHN, M., 1994: Modellierung der Höhenentwicklung ausgewählter Baumarten in Abhängigkeit vom Standort. Forstliche Forschungsberichte München, Nr. 141: 204 pp. KARR, J.R.; ROTH, R.R., 1971: Vegetation structure and avian diverstiy in several new world areas. Am. Nat. 105: 423–435. KARSHOLT, O.; RAZOWSKI,J.,1996: Lepidoptera of Europe. A distributional checklist. Apollo Books, Stenstrup, DK. 380 pp. KERR, J.F.; SUGAR, A.; PACKER, L., 2000: Indicator Taxa, rapid biodiversity assessment and nestedness in an endangered ecosystem. Conserv. Biol. 14, 6: 1726–1734. KIRBY, K.J.; REID, C.H.; THOMAS, R.C.; GOLDSMITH, F.B., 1998: Preliminary estimates of fallen dead wood and standing dead trees in managed and unmanaged forests in Britain. J. Appl. Ecol. 35, 1: 148–155. KIRCHHOFF, B., 2003: Erfassung und Beschreibung ausgewählter Waldflächen der Nordmongolei. Master thesis, University of Göttingen. 84 pp. KLAUSNITZER, B., 1999: Wollen wir wirklich Land mit einer reichen Tierwelt, und denken wir dabei auch an die Kleinen, besonders die Insekten? Natur- und Kulturlandschaft 3: 118–131. KOZLOVA, E.V., 1930: Birds of south-east Transbaikalia, North Mongolia and Central Gobi. Materialy Komisii AN SSSR po issledovaniu Mongolskoi I Tuvinskoi Narodnyh Respublik I Buriat-Mongolskoi ASSR, 12: 1–396. KRAMER, H.; AKÇA, A., 1995: Leitfaden zur Waldmesslehre. Frankfurt, J. D. Sauerländer´s Verlag. 266 pp. KRAMER, H., 1988: Waldwachstumslehre. Hamburg und Berlin, Verlag Paul Parey. KREBS, C.J.; BOUTIN, S.; BOONSTRA, R. (eds), 2001: Ecosystem Dynamics of the Boreal Forest. The Kluane Project. Oxford University Press. 511 pp. KUJAWA, K., 1997: Relationships between the structure of mid-field woods and their breeding bird communities. Acta Ornithol. 32, 2: 175–184. LANDRES, P.B.; VERNER, J.; THOMAS, J.W., 1988: Ecological uses of vertebrate indicator species: a critique. Conserv. Biol. 2: 316–328 LEE, W.K.U.; GADOW,K.V., 1996: Iterative Bestimmung der Konkurrenzbäume in Pinus densi- flora Beständen. Allg. Forst- Jagdztg. 168, 3–4: 41–45. LINDENMAYER, D.B.; INCOLL, R.D.; CUNNINGHAM, R.B.; DONELLY, C.F., 1999: Attributes of logs on the floor of Australian Mountain Ash (Eucalyptus regans) forests of different ages. For. Ecol. Manage. 123: 195–203. LOHR, S.M.; GAUTHREAUX, S.A.; KILGO, J.C., 2002: Importance of Coarse Woody Debris to Avian Communities in Loblolly Pine Forests. Conserv. Biol. 16, 3: 767–777. MAC NALLY, R.; HORROCKS,G.;PETTIFER, L., 2002: Experimental evidence for potential benefical effects of fallen timber in forests. Ecol. Appl. 12, 6: 1588–1594. MACARTHUR, R.H.; MACARTHUR, J., 1961: On bird species diversity. Ecology 42: 594–598. MACARTHUR, R.H., 1964: Environmental factors affecting bird species diversity. Am. Nat. 98: 387–397. For. Snow Landsc. Res. 78, 1/2 (2004) 117

MIKUSINSKI, G.; GROMADZKI, M.; CHYLARECKI, P., 2001: Woodpeckers as indicators of forest bird diversity. Conserv. Biol. 15, 1: 208–217. MÖNCHBAYAR, H., 1999: Some results of the species composition and ecology of Lepidopterous (Rhopalocera) of the Bogd Khan Uul. Magisterarbeit, Ulan Bator, Mongolei. Mongolian Ministry for nature and environment, 1996: United Nations development program, Global environment facility, World wide fund for nature: Mongolias Wild Heritage, Ulaan- baatar. MÖNKKÖNEN, M., 1994: Diversity pattern in palearctic and nearctic forest bird assemblages. J. Biogeogr. 21: 183–195. MÖNKKÖNEN, M.; VIRO, P.,1997: Taxonomic diversity of the terrestrial bird and mammal fauna in temperate and boreal biomes of the northern hemisphere. J. Biogeogr. 24: 603–612. MÜHLENBERG, M., 1980: Kann man mit Hilfe von Vegetationsstrukturmessungen Veränderungen von Tiergemeinschaften bestimmen? Vaduz, Bericht des Internationalen Symposiums der Internationalen Vereinigung für Vegetationskunde, Rinteln 1979, J. Kramer. 309–327. MÜHLENBERG, M.; SLOWIK, J.; SAMIYA, R.; DULAMSUREN, C.; GANTIGMAA, C.; WOYCIE- CHOSWKI, M., 2000: The Conservation Value of West-Khentii, North Mongolia: Evaluation of Plant and Butterfly Communities. Fragm. Flor. Geobot. 45, 1–2: 63–90. National Atlas of the Mongolian People’s Rebublik, Ulan Bator 1990: Mongol Ulsyn Undesnii Atlas. Edited in Ulan Bator and Moscow, printed in Minsk, cartography (in mongolian). NIEMELÄ, J., 1997: Invertebrates and Boreal Forest Management. Conserv. Biol. 11, 3: 601–610. NILSSON, C.; GRELSSON, G.; JOHANSSON, M.; SPERENS U., 1988: Can rarity and diversity be predicted along river banks? Biol. Conserv. 44: 201–212. NILSSON, S.G.; ARUP, U.; BARANOWSKI, R.; EKMAN, S., 1995: Tree dependent lichens and beetles as indicators in conservation forests. Conserv. Biol. 9, 5: 1208–1215. PEARSON, D.L.; CASSOLA, F., 1992: World wide species richness patterns of Tiger Beetles (Coleoptera: Cincidelidae): Indicator taxon for biodiversity and conservation studies. Conserv. Biol. 6, 3: 376–391. PIELOU, E.C., 1977: Mathematical Ecology. New York, John Wiley & Sons. 384 pp. POMMERENING, A., 2002: Approaches to quantifying forest structures. Forestry 75, 3: 305–324. PRETZSCH, H., 2001: Modellierung des Waldwachstums. Berlin, Parey. 341 pp. RAINIO, J.; NIEMELÄ, J., 2003: Ground beetles (Coleoptera: Carabidae) as bioindicators. Biodiv. Conserv. 12: 487–506. Redkie Zivotnye Mongolii (pozvonocnye), 1996: Moskva, 182 pp. ROJO, A.U.; MONTERO, G., 1996: El Pino Silvestre en la Sierra de Guadarrama. Ministerio de Agriculture, Madrid, Pesca y Alimentacion. 293 pp. SAVIN, E.N.; MILYITIN, L.I.; KRASNOSCEKOV, Y.N.; KOROTKOV, I.A.; SUNLOV, A.V.; DUGARZZAV, C.; DASZEVEG, T.S.; COGOO, Z.; DORJSUREN, C.; ZHAM’JANSUREN, S.; GOMBO SUREN,N., 1988: Wald der Mongolei (Lärchenwälder im Ostchentej). Moskva “Nauka” 176. SCHOBER, R., 1987: Ertragstafeln wichtiger Baumarten. Frankfurt a. M., J.D. Sauerländer´s Verlag. 154 pp. SPELLERBERG, I.F.; SAWYER, W. D., 1999: An introduction to applied biogeography. Cambridge University Press. 243 pp. STAUPENDAHL, K., 1997: Ein neues Stichprobenverfahren zur Erfassung und Beschreibung von Naturverjüngung. In: Jahrestagung der Sektion Forstliche Biometrie und Informatik des Deutschen Verbandes Forstlicher Forschungsanstalten in Freiburg i. Br., 24.–26.9. 1997: 32–49. STORK, N.E.; BOYLE, T.J.B.; DALE, V.; EELEY, H.; FINEGAN, B.; LAWES, M.; MANOKARAN,N.; PRABHU, R.; SOBERON, J., 1997: Criteria and indicators for assessing the sustainability of forest management: Conservation of Biodiversity. CIFOR Working Paper 17:1–29. STOYAN, D.; STOYAN, H., 1992: Fraktale Formen Punktfelder. Methoden der Geometrie-Statistik. Berlin, Akademie-Verlag. 394 pp. SULLIVAN, T.P.; SULLIVAN, P.S.; LINDGREN, P.M., 2001: Stand structure and small mammals in young lodgepole pine forest: 10-year research after thinning. Ecol. Appl. 11, 4: 1151–1173. SUTHERLAND, W.J., 1998: The importance of behavioural studies in conservation biology. Anim. Behav. 56: 801–809. 118 Michael Mühlenberg et al.

TAYLOR, R.J.; DORAN, N., 2001: Use of terrestrial invertebrates as indicators of the ecological sustainability of forest management under the Montreal Process. J. Ins. Conserv. 5: 221–231. TRETER, U., 1997: Gebirgs-Waldsteppe in der Mongolei. Geogr. Rundsch. 5: 655–661. TSEDENDASCH, G., 1995: Khentejn huruuny oj – urgamalschil (Forest Vegetation of the Khentii Mountains) PhD thesis National University of Mongolia, Ulan Bator (in mongolian). ULZIJCHUTAG, N., 1989: Mongol ornj urgamljn aimgijn tojm, Ulan Bator. VASILCHENKO, A.A., 1987: Birds of Khamar-Daban. Novosibirsk, Nauka Press, Siberian dep., 103 pp. (in russian). WEAVER, J.C., 1995: Indicator Species and Scale of Observation. Conserv. Biol. 9, 4: 939–942. WENK, G.; ANTANAITIS, V.U.; SMELKO, S., 1990: Waldertragslehre. Berlin, Deutscher Landwirt- schaftsverlag. 448 pp. WICHMANN,F.,2001: Analyse der Vogelgemeinschaft im Habitatmosaik einer Naturlandschaft im Norden der Mongolei, Westchentie. Diplomarbeit Univeristät Göttingen. WILLSON, M.F., 1974: Avian Community organization and habitat structure. Ecology 55: 1017–1029. WOYCIECHOWSKI, M.; SLOWIK, J.; MÜHLENBERG, M., 2001: Red listed species beyond the European border: Maculinea species of Mongolia. Poster at the International Symposium for European Butterflies, ANL Laufen, December 2001. ZENNER, E.K.; HIBBS, D.E., 2000: A new method for modeling the heterogeneity of forest struc- ture. For. Ecol. Manage 129, 1: 75–87.

Accepted April 9, 2004