Effects of residual trees on tree height of 18- and 19-year-old Cryptomeria japonica planted in group selection openings

This study examined the effects of residual trees on the height of 18- and 19-year-old Cryptomeria japonica planted in group selection openings of about 0.1 ha at Yufuin in Oita Prefecture in southern Japan. We first developed a general model expressing variation in the height for all of the planted trees within the openings from indices of both between-cohort competition (the effects of residual trees) and within-cohort competition (among the even-aged planted trees). The between-cohort competition index was calculated both with and without taking into account the direction that residual trees were positioned from subject trees. The two models suggested that residual trees located near the northern edge of the opening were not effective competitors for light for the planted trees within the opening. In addition, our results indicated that both between-cohort competition and within-cohort competition are important for explaining tree height in closed stands within group selection openings. Using the general model, we then estimated the potential loss of height growth for locally dominant trees within the openings that could be attributed only to the influence of residual trees. These results indicated that in addition to the effect of opening size, both the degree of crowding and the direction of residual trees from the planted trees are important factors affecting the height of individual trees planted in group selection openings.     

Effect of topography and geology on the site index of Picea abies in the West Carpathian,Poland

Abstract

The effects of topographic features and the soil parent rock on the site index of Norway spruce [Picea abies (L.) Karst.] were examined. Field data were collected on 347 plots located in the Polish West Carpathian. On these plots, the diameter at breast height and the height of trees were measured and the age of trees was determined. On 64 plots one dominant or codominant tree was felled for stem analysis. In addition, topographic features of each site were measured, including elevation above sea level, aspect and slope. Stem analysis data were used in the development of a local site index system, which was used for calculating the site index for particular plots. In order to interpret site index as functions of elevation, aspect, slope, size of a mountain massif and soil parent rock, the multiple regression model was used, which explained 79.0% of the variation in Norway spruce site index. Elevation above sea level was the major factor responsible for variation in site index, and was negatively correlated (R ² _adj=0.610) with site index. Site index also depended on the soil parent rock. The stands growing on the tertiary Magura sandstones had considerably higher site indices than those located on the geological substratum composed of Istebna and Godula sandstones formed during the Cretaceous period. The site index of Norway spruce is also determined by the aspect, slope and height of mountain massif.     

Predicting the spatial distribution of major species composition in secondary hardwood forests on Mt. Gozu,central Japan,based on environmental factors

Our goal was to predict the spatial distribution of canopy species composition in secondary deciduous hardwood forests at a fine spatial resolution, based on climatic and topographical factors using a geographic information system. We studied secondary forests on Mt. Gozu, Niigata Prefecture, central Japan. Canopy species composition was investigated in 100 sample plots within the study area. A digital elevation model (DEM) was created, and topographical, hydrological, and light factors were calculated using the DEM. Climatic factors were interpolated by kriging. The five major species used as response variables were Fagus crenata, Quercus serrata, Quercus crispula, Magnolia obovata, and Castanea crenata. We prepared three possible explanatory variable sets: climatic variables only, both climatic and topographic variables, and topographic variables only. Multivariate regression trees were derived, and the accuracy of predicting the major species composition was tested. The multivariate regression trees derived from the climatic variable set and from the climatic and topographic variable set had better accuracy than the regression tree derived from the topographic variable set. In the regression tree generated by the climatic and topographic variable set, the warmth index was the principal explanatory variable in classifying forest types, followed by topographic factors. This regression tree would be preferable to the other two regression trees for the prediction of canopy species composition.     

Relationships of stem surface area to other stem dimensions for Japanese cedar (<Emphasis Type="Italic">Cryptomeria japonica</Emphasis> D. Don) and Japanese cypress (<Emphasis Type="Italic">Chamaecyparis obtusa</Emphasis> Endl.) trees

I investigated the relationships between stem surface area outside bark and other stem dimensions for Japanese cedar (Cryptomeria japonica D. Don) and Japanese cypress (Chamaecyparis obtusa Endl.) trees. The stem dimensions used here were the basal area and the product of diameter at breast height (dbh) and total tree height. The regression equation of the stem surface area s against basal area g was s = 184.216g for the cedar trees and s = 156.878g for the cypress trees. The slope of the equation was significantly different between the two species. For the same dbh, the cedar trees tended to have a higher total tree height than the cypress trees. The cedar trees also had a larger surface area of relative stem form than the cypress trees. These results indicated that the difference in the slope was produced by the differences in both the stem slenderness and tapering between the two species. On the other hand, the regression equation between the stem surface area and the product of dbh and total tree height dh was s = 1.937dh for the cedar trees and s = 1.921dh for the cypress trees, whereas no significant difference in the slope was found. The obtained slopes for the cedar and cypress trees seemed to be in accord with that for other coniferous species reported in earlier studies, suggesting that the variation in the slope among coniferous species would be small. The estimation from the basal area would provide a simpler means for estimating the stem surface area and would be useful in obtaining an approximation of the surface area. By contrast, the estimation from the product of dbh and total tree height would provide a more accurate and precise estimate as well as a wider applicable range, i.e., a parameter for physiological growth models. In conclusion, it could not be judged which regression equation examined in the present study was superior to the other, and thus it was important to select an appropriate equation depending both on the purpose and on the time and labor available.     

Potential of natural and improved fallow using indigenous trees to facilitate cacao replanting in Ghana

Poor establishment, due to loss of soil fertility, weeds and lack of appropriate shade, is a major constraint to replanting cacao on previously used land. Spathodea campanulata, Newbouldia laevis and Ricinodendron heudelotii planted as monospecific improved fallow and Terminalia ivorensis, T. superba and Antiaris toxicaria planted as a multispecies improved fallow and a natural tree fallow were assessed for their potential to facilitate cacao replanting in a randomized complete block design experiment. Simpson and Shannon diversity indices and species richness in the natural tree fallow were 0.6, 1.6 and 20, respectively, at 4 years after trial inception. The Multispecies and the R. heudelotii improved fallows had better height growth, crown development and light transmission characteristics, which are desirable for cacao shade. However, these were not comparable to S. campanulata or the natural tree fallow in terms of improving microsite topsoil pH, % organic carbon and % total nitrogen and site capture. Since optimum fallow period is shortened by growing fast-growing trees, the height growth rate >2.0 m per annum in all the treatments except N. laevis indicates the suitability of these species for improved fallow. The trees species showed different and complementary characteristics and from a standpoint of biodiversity conservation and the future floristic composition of the landscape the natural tree fallow with its diversity of tree species may be recommended as a rehabilitation technique to facilitate the replanting of cacao with a diverse overhead shade.     

Predicting tree regeneration in <Emphasis Type="Italic">Picea abies</Emphasis> snag stands

A bark beetle (Ips typographus) infestation caused the death of almost all Norway spruce (Picea abies) trees in a mountain forest in the Swiss Alps. We developed a tree regeneration model, ‘RegSnag’ (=REGeneration in a SNAG stand), to project the future amount and height of tree regeneration in these snag stands. The model combines a height-class structured tree module with a microsite-based module of snag decay and ground-vegetation succession. Microsite-specific rates of germination, mortality and height growth were modelled for four tree species (Picea abies, Sorbus aucuparia, Acer pseudoplatanus and Betula pendula) in eight height classes (from seedlings to saplings 5 m tall) and on 26 microsite types (e.g. moss, grass). Model tests with independent field data from 8 years after the Picea die-back demonstrated that microsites had a considerable effect on the development of tree regeneration on both the montane and the subalpine level. With microsite-specific parameters, the height and frequency of Picea in each microsite could be simulated more accurately than without considering microsite effects (e.g. bias of 8 vs. 119 saplings ha⁻¹ on the montane level). Results of simulations 40 years into the future suggest that about 330–930 Picea saplings per ha out of those that germinated in 1994 and 1996 will reach a height of 5 m within 30–35 years after Picea die-back. This is due to differences in seed inflow and browsing intensities. Picea and not Betula or Sorbus trees will replace the current herbaceous vegetation in these snag stands.     

Tree species composition,biodiversity and regeneration in response to catena shape and position in a mountain forest

We examined the variation in woody plant species composition, regeneration and biodiversity in relation to topographic heterogeneity, specifically catena shape (V and Λ) and slope position (summit, shoulder, back slope, foot slope and toe slope). All research took place at the TMU Experimental Forest Station, located in northern Iran. Two sets of replicated quadrats were established along transects that encompassed the different catena shapes and positions. We recorded the density, height, and canopy cover of all trees ≥12 cm diameter at breast height within 30 quadrats (400m²), as well as the density of all tree regenerations within 120 quadrats (4m²). In addition, three measurements of biodiversity were calculated for both mature trees and regenerations. Mature tree species composition varied greatly with catena shape and slope position, as did that of regenerating trees. Despite differences in species composition, the various measurements of biodiversity were consistent across the landscape for mature trees, whereas significant variation in evenness was detected for regenerating trees. We concluded that the broad matrix of environmental conditions occurring within this rugged mountainous area generates a variety of microhabitats, which in turn contributes to both the variation in species composition and the consistent biodiversity of these forests.     

Evaluating tree competition indices as predictors of basal area increment in western Montana forests

Fire hazard reduction treatments are commonly applied to mixed-species coniferous forests in western Montana, USA, to modify fuels structures and alter the competitive environments of individual trees. An improved understanding of how competition can be measured and how it conditions individual tree growth is needed for projecting the development of these forests, with and without treatment. Numerous studies have evaluated how competition affects tree growth and many indices have been developed to quantify the competition an individual tree experiences. These studies suggest that no single competition index or a single class of indices is universally superior; indices perform differently according to forest type and forest conditions. We chose several widely used distance-independent and distance-dependent competition indices, and also derived anisotropic distance-dependent indices from estimates of light interception by tree crowns. We evaluated the effectiveness of these competition measures for predicting basal area increment (BAI) of Pinus ponderosa, Pseudotsuga menziesii, and Larix occidentalis in western Montana. The best distance-dependent competition indices explained a larger proportion of growth variation than the best distance-independent indices (64% vs. 56%). This result indicates that competition is an important growth determinant in these forests and that competition varies locally, with variable tree densities and relatively complex stand structures creating heterogeneous neighborhood conditions. Competition indices derived from light interception models were only weakly correlated with other indices and performed poorly in terms of predicting tree growth. This result accords with previous observations that competition for light is not the primarily growth limitation for trees in the semi-arid conditions of western Montana. More sophisticated light availability models could be used to better assess variability in light interception and its marginal contribution to predictive accuracy of radial tree growth. Diameter and distance-dependent BAI models were developed for growth prediction at the species level and for all species combined.     

Estimating individual tree heights of sugi (Cryptomeria japonica D. Don) plantations in mountainous areas using small-footprint airborne LiDAR

Recently, it was shown that individual tree heights could be accurately estimated using small-footprint airborne light detection and ranging (LiDAR) remote sensing. Because most of the areas studied previously were limited to flat terrain, we investigated the accuracy of LiDAR-derived individual tree height estimates for different types of topographical features in mountainous forests with a steeper and more complex topography. Several middle-aged (40–50 years old) sugi (Cryptomeria japonica D. Don) plantations are found in the mountainous regions in Japan; hence, we chose 48-year-old sugi plantations to investigate the accuracy of these estimates. The surveyed area was divided into three types of topographical features; steep slope (mean slope ± SD; 37.6° ± 5.8°), gentle slope (15.6° ± 3.7°), and gentle yet rough terrain (16.8° ± 7.8°). Before estimating tree heights, the number of detected trees within each topographical feature was researched. In each of these terrains, the percentage of trees detected correctly was 74%, 86%, and 92%; the average error between LiDAR-derived and field-measured tree heights was 0.227m, –0.473m, and –0.183m; and the accuracy of the LiDAR-derived tree height estimates, given as root mean square error (RMSE), was 0.901m, 0.846m, and 0.576m, respectively. Consequently, the procedure presented in this study could detect most canopy trees and estimate individual tree heights with an accuracy better than 1m, even in a forest with a mean slope angle of approximately 38°; thus, indicating that small-footprint airborne LiDAR will be a useful tool for accurately estimating the heights of individual canopy trees in sugi plantations in mountainous areas.     

Individual-based measurement and analysis of root system development: case studies for <Emphasis Type="Italic">Larix gmelinii</Emphasis> trees growing on the permafrost region in Siberia

We present results of individual-based root system measurement and analysis applied for Larix gmelinii trees growing on the continuous permafrost region of central Siberia. The data of root excavation taken from the three stands were used for the analyses; young (26 years old), mature (105 years old), and uneven-aged over-mature stand (220 years old). In this article, we highlight two topics: (1) factors affecting spatio-temporal pattern of root system development, and (2) interactions between aboveground (i.e., crown) and belowground (i.e., root) competition. For the first topic, the detailed observation of lateral roots was applied to one sample tree of the overmature stand. The tree constructed a superficial (<30 cm in depth) and rather asymmetric root system, and each lateral root expanded mainly into elevated mounds rather than depressed troughs. This indicated that spatial development of an individual root system was largely affected by microtopography (i.e., earth hummocks). For these lateral roots, elongation growth curves were reconstructed using annual-ring data, and annual growth rates and patterns were compared among them. The comparison suggested that temporal root system development is associated with differences in carbon allocation among the lateral roots. For the second topic, we examined relationships between individual crown projection area (CA) and horizontal rooting area (RA) for the sample trees of each stand. RA was almost equal to CA in the young stand, while RA was much larger (three or four times) than CA in the mature and overmature stands. Two measures of stand-level space occupation, crown area index (aboveground: CAI; sum of CAs per unit land area) and rooting area index (belowground: RAI; sum of RAs), were estimated in each stand. The estimates of RAI (1.3–1.8 m² m₋₂) exceeded unity in all stands. In contrast, CAI exceeded unity (1.3 m² m₋₂) only in the young stand, and was much smaller (<0.3 m² m₋₂) in the two older stands. These between-stand differences in RAI–CAI relationships suggest that intertree competition for both aboveground and belowground spaces occurred in the young stand, but only belowground competition still occurred in the two older stands. Based on this finding, we hypothesized that competition below the ground may become predominant as a stand ages in L. gmelinii forests. Methodological limitations of our analysis are also discussed, especially for the analysis using the two indices of space occupation (CAI, RAI).     

Nondestructive sampling of Eucalyptus globulus and E. nitens for wood properties. III. Predicted pulp yield using Near Infrared Reflectance Analysis

  Within-tree variation in kraft pulp yield, predicted using near infrared reflectance analysis, was studied in thirty trees of E. globulus and fifty trees of E. nitens to develop a non-destructive sampling strategy. Trees, aged 5 to 9 years, were sampled across a range of sites in southern Australia. Simulated core samples were removed at six fixed heights easily accessible from the ground (0.5, 0.7, ... 1.5 m) and at seven percentage heights (0, 20, 30, ... 70%). Whole-tree values, calculated from percentage height data, were correlated with the core data to determine the optimal sampling height. Core samples were found to be good predictors of whole-tree pulp yield for E. globulus, with simulated cores taken from the recommended sampling height (1.1 m) explaining more than 50% of variation in whole-tree pulp yield. Results for E. nitens were variable with large site differences apparent. On high quality sites, core samples from the recommended sampling height (0.9 m) were good predictors of whole-tree pulp yield, explaining around 60% of the variation. On poor quality sites, cores were poor predictors of whole-tree pulp yield. Radial orientation of cores was not important and predicted pulp yield was not related to tree size, basic density or fibre length. To estimate stand mean pulp yield to an accuracy of ±1% would require sampling 6 trees of E. globulus and 4 trees for E. nitens using either multiple discs or core samples. A single sampling height (1.1 m) is recommended for sampling for basic density, fibre length, fibre coarseness and predicted pulp yield in E. globulus. For E. nitens the recommended sampling height for basic density and fibre length is 0.7 m and 0.9 m is recommended for predicted pulp yield on good quality sites. Received 17 September 1998     

帽儿山实验林场白桦人工幼林适宜微立地研究

在相同立地上的不同类型的微立地环境中,测定白桦幼树生长指标和土壤含水率,并以邻近微环境外的4株白桦作为对照,分析微立地环境对白桦幼树生长的影响。研究结果表明：由伐根和倒木等粗大木质残体构成的微立地对白桦幼树的生长均产生明显的促进作用。凸凹地形微立地对幼树生长的作用因所处立地条件而异,凸地形、凹地形分别在下坡位和上坡位有利于白桦幼树生长。灰烬堆微立地在上坡位不利于白桦的生长,中、下坡有利。白桦人工造林时,根据立地条件科学选择微立地作为栽植点可以有效改善幼林的生长表现。     

Phenotypic plasticity of two invasive alien plant species inside a deciduous forest in a strict nature reserve in Italy

ABSTRACT

Invasive alien plant species (IAPs) represent one of the main biological threats to biodiversity worldwide. Information about their phenotypic plasticity are needed to increase awareness about their future invasive potential. A study about phenotypic plasticity in response to contrasting light regimes and its quantification by a plasticity index (PI) of two IAPs (Ailanthus altissima and Robinia pseudoacacia) inside a Strict Nature Reserve was conducted. R. pseudoacacia showed a 70% higher PI, with a strongly greater value at morphological leaf level, associated with a greater ability to survive and grow in forest understory, explaining its greater widespread. Otherwise, A. altissima showed its highest PI at physiological level, which was associated with the ability to colonize and grow in environments with high-light regimes. Based on these results, the conservative management has limited the presence of A. altissima by its lower ability to grow in forest understory. In fact, the small-scale gaps in the forest infrastructure, that could allow its recruitment, are originated only from the death of a single tree or small group of trees. Regarding R. pseudoacacia, it is critical to maintain this type of management because any disturbances resulting in large openings could further promote its presence inside the Reserve.     

The contribution of structural indices to the modelling of Sitka spruce (Picea sitchensis) and birch (Betula spp.) crowns

Crown dimensions are important for the quantification of tree interactions in some growth models. This study investigates the potential for structural indices and other spatial measures to improve the prediction of crown radius and crown length for birch (Betula spp.) and Sitka spruce (Picea sitchensis (Bong.) Carr.) in forests in Wales. Crown dimensions were measured for 125 birch and 154 spruce in six fully stem-mapped research plots. These data were used to test the performance of a crown radius model and a crown length model which estimated crown dimensions on the basis of allometric relationships with stem dimensions. Spatial data from the six plots were used to calculate the structural indices mean directional index, diameter correlation index, species mingling, dbh and height dominance, and dbh differentiation, as well as the Hegyi competition index, and basal area of neighbours and larger neighbours, for each crown measurement sample tree, using various numbers of nearest neighbours. Two non-spatial indices, BAL and BALMOD, were also calculated for all sample trees for comparison. These spatial and non-spatial variables were then incorporated into modified crown dimension models. Model performances, in terms of efficiency and relative bias, were compared to determine whether the inclusion of spatial or non-spatial variables resulted in any improvements over models using tree dimensions alone. Crown length and radius were found to be correlated with most of the spatial measures studied. Models incorporating spatial variables gave improvements in performance over allometric models for every data set, and performed more consistently than models containing non-spatial variables. The greatest improvements were achieved for suppressed birch in unthinned forests which had irregularly shaped and strongly displaced crowns. The spatial variable contributing to the most efficient model for each data set varied widely. This points to the complexity of tree spatial interactions and indicates that there is a great deal of scope for investigating other structural indices and crown dimension model forms.     

Nondestructive sampling of Eucalyptus globulus and E. nitens for wood properties. I. Basic density

A non-destructive sampling strategy for basic density, based on removing 12 mm bark-to-bark cores, was developed in E. globulus and E. nitens. Fifty trees of each species, aged 5 to 9 years, were sampled across a range of sites. Core samples were removed on both a north-south and an east-west axis from 6 fixed heights in the base of the tree (0.5 m, 0.7 m, …, 1.5 m). Whole-tree values were calculated from disc samples removed at eight percentage heights (0, 10, 20, …, 70%) and correlations between the cores and whole-tree values were used to determine the optimal sampling height. Core samples were found to be reliable predictors of whole-tree density, explaining between 84% and 89% of the variation between trees. Core sampling of E. globulus and E. nitens to estimate basic density of whole-trees and stands is feasible; cores from trees at all E. globulus sites gave high correlations with whole-tree values. For E. nitens, site differences were apparent, and it is recommended that a small destructive sampling program should be undertaken prior to commencing a major sampling program. Recommended optimal sampling heights are 1.1 m for E. globulus and 0.7 m for E. nitens. Core orientation was not important and density was not related to tree size. Six whole-tree samples or eight core samples are required for estimating the mean density of a stand at a specific site to an accuracy of ±20 kg m⁻³ with a 95% confidence interval. Received 17 September 1998     

Nondestructive sampling of Eucalyptus globulus and E. nitens for wood properties; II. Fibre length and coarseness

Within-tree variation in fibre length and coarseness was studied in fifty trees of E. globulus and E. nitens to develop a non-destructive sampling strategy. Trees, aged 5 to 9 years, were sampled across a range of sites in southern Australia. Simulated core samples were removed at six fixed heights easily accessible from the ground (0.5, 0.7, ... 1.5 m) and at eight percentage heights (0, 10, 20, ... 70%). Whole-tree values, calculated from percentage height data, were correlated with the core data to determine the optimal sampling height. Core samples were found to be reliable predictors of whole-tree fibre length, but results were variable for fibre coarseness. Simulated cores taken from the recommended sampling heights explained 87% and 71% of variation in whole-tree fibre length for E. globulus and E. nitens respectively and 54% and 45% of the variation in whole-tree fibre coarseness. Fibre length at all fixed heights showed good correlations with whole-tree values at all sites for E. globulus. For E. nitens the correlations were slightly lower and variable across sites. Results for fibre coarseness varied across sampling heights and sites for both species. The recommended sampling height for fibre length is 1.5 m for both species, whilst for fibre coarseness, the recommended sampling heights are 0.9 and 1.1 m for E. globulus, and 0.9 and 1.3 m for E. nitens. Radial orientation of cores was not important and neither fibre length nor coarseness were related to tree size or basic density. To estimate stand mean fibre length to an accuracy of ±5% would require sampling 9 whole trees or taking cores from 13 trees for E. globulus and 4 whole trees or cores from 8 trees for E. nitens. For estimating stand mean fibre coarseness, 10 whole trees of E. globulus and 7 whole trees are needed for E. nitens. Core sampling for stand mean coarseness would require more trees: 13 to 21 for E. globulus and 11 to 16 trees for E. nitens. Received 17 September 1998     

Effects of topography and planted trees on the distribution of naturally regenerated broad-leaved trees in a 140-year-old <Emphasis Type="Italic">Cryptomeria japonica</Emphasis> plantation in northern Kyushu,Japan

We investigated factors affecting the distribution of naturally regenerated broad-leaved trees in a 140-year-old Cryptomeria japonica plantation. We used path analysis to examine the relationship among microtopography, the biomass of planted trees, and the biomasses of canopy and subcanopy trees of broad-leaved species. The study plot was divided into three topographic types (ridge, slope, and valley), and we discuss how the different topographic types are affected. For all topographic types, the biomass of canopy trees of broad-leaved species decreased with convexity. For slope and valley topographies, the biomass of subcanopy trees of broad-leaved species also decreased with convexity. For ridge topography, the biomass of subcanopy trees of broad-leaved species increased with the biomass of planted trees, and decreased with the biomass of canopy trees of broad-leaved species. These results suggest the effects of microtopography on the biomass of subcanopy trees were much larger than the effects of canopy trees for slope and valley topographies, while the effects of microtopography were smaller for ridge topography.     

Factors influencing tree species diversity and Betula alleghaniensis establishment in silvicultural openings

We examined the short-term effects of group-selection harvestingwith seed tree retention on the diversity, abundance and establishmentof tree seedlings in a northern hardwood forest in the UpperPeninsula of Michigan (49 openings, 20 closed canopy referencesites). Three opening sizes were examined – opening radius0.5 x canopy height (267 ± 62 m², n = 16), 0.75 x canopyheight (642 ± 85 m², n = 17) and 1.0 x canopy height(1192 ± 155 m², n = 16) (canopy height = 22 m). A singleyellow birch (Betula alleghaniensis Britt.) seed tree was retainedin the centre of each opening. Tree seedling density was significantlyhigher in the largest group-selection opening than at the closedcanopy reference sites (P < 0.05), the main factor for thiswas the increased proportion of yellow birch, red maple (Acerrubrum L.) and other minor species. Nevertheless, yellow birchwas still a minor component of the developing gap cohort, comprising5.9 per cent of the seedlings and 1.1 per cent of the saplings.Within openings, microsite variables, such as per cent coversof bare soil and coarse woody debris, were the best predictorsof yellow birch occurrence and density. Our results suggestthat microsite limitations and competing vegetation may greatlyreduce the efficacy of openings for ensuring the maintenanceof mid-tolerants.     

Modelling trees' effect on maize in the Grevillea Robusta + maize system in Central Kenya

In the Central Highlands on Kenya, Grevillearobusta (A. Cunn.) trees are commonly planted on maize (Zeamays L.) fields. This type of production system covers about 750 000ha in the highland area around Mt Kenya. Knowing the influence oftrees on maize yield would help to regulate the density and size distribution ofthe tree cover. The spatial distribution of trees varies from rather uniform toextremely aggregated patterns, calling for a spatial modelling approach. Thisstudy employed three different spatial competition index types to model theeffect of trees on the maize yield at a given location. The maize yield wasexpressed as a function of competition index, which was calculated from thediameters and distances of trees. The data were collected from 14 maize fields,and included measurements on 857 trees and 2514 one-m² maize plots. The models indicated that high competition by treesdecreases maize yields considerably, but at the field level the decrease issmall with normal tree stockings (about 200 trees per hectare). Characteristicto the models was low degree of explained variance, which was partly due to highsampling error in maize yield measurement and low variation in competition(places with high competition by trees were seldom planted with maize). One ofthe maize yield models was combined with an earlier simulation program forG. robusta. The simulator allows one to predict thetemporal development of any G. robusta – maize fieldin Central Kenya. This revised version was published online in June 2006 with corrections to the Cover Date.