Current issue: 56(2)
Under compilation: 56(3)
New mortality models were developed for the purpose of improving long-term growth and yield simulations in Finland, Norway, and Sweden and were based on permanent national forest inventory plots from Sweden and Norway. Mortality was modelled in two steps. The first model predicts the probability of survival, while the second model predicts the proportion of basal area in surviving trees for plots where mortality has occurred. In both models, the logistic function was used. The models incorporate the variation in prediction period length and in plot size. Validation of both models indicated unbiased mortality rates with respect to various stand characteristics such as stand density, average tree diameter, stand age, and the proportion of different tree species, Scots pine (Pinus sylvestris L.), Norway spruce (Picea abies (L.) Karst.), and broadleaves. When testing against an independent dataset of unmanaged spruce-dominated stands in Finland, the models provided unbiased prediction with respect to stand age.
Vaccinium myrtillus L., V. vitis-idaea L. and V. uliginosum L. belong to the genus Vaccinium. These wild species are widely distributed and ecologically important within the Baltic countries but they have not been extensively studied using molecular markers. EST-SSR and cpSSR markers were used to investigate the population structure and genetic diversity of these species to obtain information useful for the development of in situ conservation strategies. Wild Vaccinium species populations are moderately genetically differentiated, with some populations more highly differentiated, but without higher order clustering of groups of populations, indicating that there are no dispersal barriers for these species within the Baltic countries. Genetic diversity of populations growing in protected areas, managed forests and intensively utilised public recreational areas is similar. The results from this study can be utilised for the selection of populations for the in situ conservation of the studied Vaccinium species. In addition, complementary ex situ conservation strategies can be used for the preservation of rare varieties (e.g. V. myrtillus var. leucocarpum).
Expertise in the cost-efficient utilization and treatment of brushwood on forest roadside sites is limited. In the present study, the productivity of brushwood clearing and harvesting on forest roadside sites was defined by creating time-consumption models or parameters for the aforementioned working methods. Compiled time consumption models and parameters for the brushwood clearing and harvesting can be used as a basis for evaluating alternative management practices and to determine when brushwood biomass should be harvested and when it should be left to decay. The harvesting of brushwood was based on the harwarder system and the clearing of brushwood was done with a spiral cutter, which is a novel accessory for cutting roadside vegetation. Based on the study results, the average volume of harvested brushwood and forwarding distance are the key elements that have an effect on harvesting productivity with harwarders. Correspondingly, stump diameter has a strong impact on the clearing productivity of brushwood. The plot-wise productivity of the spiral cutter in brushwood clearings varied in the range of 0.19–0.61 ha per PMh. An increase in stump diameter slowed down the clearing productivity of the spiral cutter and there was a clear step downward in clearing productivity as the average diameter increased from 30 mm to 40 mm. The machinery under study operated well and there were no interruptions due to machine breakdowns.
In Nordic forests, consistent evidence about better seedling survival rate and increased growth due to site preparation have been obtained in numerous studies. Proper site preparation method can reduce costs of the whole regeneration chain through its effects on survival of planted seedlings, abundance of natural regeneration and competition in early stand development. This study compared the natural regeneration of birches (silver birch (Betula pendula Roth) and downy birch (B. pubescens Ehrh.)), amount of exposed mineral soil, and growth of planted seedlings between spot mounding and inverting site preparation methods. Present study was conducted in eight forest stands established in 2012 or 2015. Even though difference was not statistically significant, inverting exposed less mineral soil than spot mounding and thus reduced the natural regeneration of birch seedlings by 6135 seedlings ha–1 compared to spot mounding. However, the variation between regeneration areas was remarkable. There was no difference in seedling mortality or growth between the site preparation methods. In order to achieve high growth of conifers, moderate amount of exposed mineral soil and thus less naturally regenerated birch, inverting should be favored over spot mounding.
The diameter at any point on a stem and tree volume are some of the most important types of information used in forest management planning. One of the methods to predict the diameter at any point on a stem is to develop taper models. Black locust (Robinia pseudoacacia L.) occurs in almost all forests in Poland, with the largest concentration in the western part of the country. Using empirical data obtained from 13 black locust stands (48 felled trees), seven taper models with different numbers of estimated parameters were analysed for section diameters both over and under bark using fixed and mixed-effects modelling approaches. Assuming a lack of additional measurements, the best fitted taper models were used for the prediction of over bark volume using both methods. The predicted volume was compared with the results from different volume equations available for black locust. The variable-form taper model with eight estimated parameters fitted the data the best. The lowest root mean square error for volume prediction was achieved for the elaborated fixed-effects taper model (0.0476), followed by the mixed-effects taper model (0.0489). At the same time, the difference between the volume relative errors achieved based on the taper models does not differ significantly from the results obtained using the volume equations already available for black locust (two of the three analysed).
Fire is a common disturbance in boreal forests causing changes in biological diversity at various spatial scales. In the past 100 years, forest management has limited fire outbreaks, but in the future, the fire-affected forest area is expected to increase in many regions due to climate change. Burned forests are typically salvage-logged, but the effect of this type of management versus natural regeneration on biological diversity is not well understood, particularly the mid-term effect to tree establishment and understory vegetation composition and diversity. Various management methods were used after a large fire in 1992 in a peatland-forest complex and neighbouring managed forests, which created an experimental setup for study of the effect of management after fire in the Sliteres National park, northwestern Latvia. Understory vegetation was described in plots using a design of four forest and three management types: natural regeneration (unmanaged) and managed sites with salvage logging followed by no further human intervention and salvage logging with planting. Post-fire management had different effect in each forest type. Species richness was higher in forest types with salvage logging than in natural regenerated sites on rich wet and rich dry forest types, but not for the poor forest types. Tree regeneration was generally greater in salvage-logged stands, but differed between forest types. Species composition was related to tree regeneration and canopy openness. In contrast to other studies, salvage logging had a positive mid-term effect to ground vegetation diversity and tree establishment in the studied stands, implying potential for concomitant management and conservation of ground cover vegetation in semi-natural stands.
The purpose of this study was to prepare a comprehensive, computerized teak (Tectona grandis L.f) plantation yield model system that can be used to describe the forest dynamics, predict growth and yield and support forest planning and decision-making. Extensive individual tree and permanent sample plot data were used to develop tree-level volume models, taper curve models and stand-level yield models for teak plantations in Panama. Tree volume models were satisfactorily validated against independent measurement data and other published models. Tree height as input parameter improved the stem volume model marginally. Stand level yield models produced comparable harvest volumes with models published in the literature. Stand level volume product outputs were found like actual harvests with an exception that the models marginally underestimate the share of logs in very large diameter classes. The kind of comprehensive model developed in this study and implemented in an easy to use software package provides a very powerful decision support tool. Optimal forest management regimes can be found by simulating different planting densities, thinning regimes and final harvest ages. Forest practitioners can apply growth and yield models in the appropriate stand level inventory data and perform long term harvest scheduling at property level or even at an entire timberland portfolio level. Harvest schedules can be optimized using the applicable financial parameters (silviculture costs, harvesting costs, wood prices and discount rates) and constraints (market size and operational capacity).
In remote sensing-based forest inventories 3D point cloud data, such as acquired from airborne laser scanning, are well suited for estimating the volume of growing stock and stand height, but tree species recognition often requires additional optical imagery. A combination of 3D data and optical imagery can be acquired based on aerial imaging only, by using stereo photogrammetric 3D canopy modeling. The use of aerial imagery is well suited for large-area forest inventories, due to low costs, good area coverage and temporally rapid cycle of data acquisition. Stereo-photogrammetric canopy modeling can also be applied to previously acquired imagery, such as for aerial ortho-mosaic production, assuming that the imagery has sufficient stereo overlap. In this study we compared two stereo-photogrammetric canopy models combined with contemporary satellite imagery in forest inventory. One canopy model was based on standard archived imagery acquired primarily for ortho-mosaic production, and another was based on aerial imagery whose acquisition parameters were better oriented for stereo-photogrammetric canopy modeling, including higher imaging resolution and greater stereo-coverage. Aerial and satellite data were tested in the estimation of growing stock volume, volumes of main tree species, basal area and diameter and height. Despite the better quality of the latter canopy model, the difference of the accuracy of the forest estimates based on the two different data sets was relatively small for most variables (differences in RMSEs were 0–20%, depending on variable). However, the estimates based on stereo-photogrammetrically oriented aerial data retained better the original variation of the forest variables present in the study area.
Forest inventories assisted by wall-to-wall airborne laser scanning (ALS), have become common practice in many countries. One major cost component in these inventories is the measurement of field sample plots used for constructing models relating biophysical forest attributes to metrics derived from ALS data. In areas where ALS-assisted forest inventories are planned, and in which the previous inventories were performed with the same method, reusing previously acquired field data can potentially reduce costs, either by (1) temporally transferring previously constructed models or (2) projecting field reference data using growth models that can serve as field reference data for model construction with up-to-date ALS data. In this study, we analyzed these two approaches of reusing field data acquired 15 years prior to the current ALS acquisition to estimate six up-to-date forest attributes (dominant tree height, mean tree height, stem number, stand basal area, volume, and aboveground biomass). Both approaches were evaluated within small stands with sizes of approximately 0.37 ha, assessing differences between estimates and ground reference values. The estimates were also compared to results from an up-to-date forest inventory relying on concurrent field- and ALS data. The results showed that even though the reuse of historical information has some potential and could be beneficial for forest inventories, systematic errors may appear prominent and need to be overcome to use it operationally. Our study showed systematic trends towards the overestimation of lower-range ground references and underestimation of the upper-range ground references.
The size of Finnish wood harvesting enterprises has grown, and entrepreneurs have become responsible for various additional tasks, resulting in networking with other harvesting enterprises of various sizes and suppliers of supporting services, but the profitability of the wood harvesting sector has remained low. In the present study, the financial performance of 83 wood harvesting companies in Eastern and Northern Finland was evaluated, based on public final account data from a five-year period between 2013 and 2017. The factors underlying economic success were identified based on 19 semi-structured entrepreneur interviews. The Business Model Canvas framework was applied in the analyses. In particular, the smallest companies (with an annual turnover of less than 600 000 €) struggled with profitability. They showed increasing indebtedness, suffered from poor power in negotiations, had typically short-term contracts, and faced difficulties in retaining skilled operators. Most of the small companies were subcontractors of larger wood-harvesting companies. The better economic success of larger companies was likely based on their capacity to provide wood harvesting services in large volumes and supply versatile services, power in negotiations, and more cost-effective operations. The future development of wood harvesting seems to be polarised: larger enterprises are likely to continue growing, while the size of smaller enterprises has stabilised. Enhancing business management skills and practices is required in enterprises of all size groups.
In wintertime, the payload capacity of a timber truck is reduced by snow that accumulates on the structures of the truck. The aim of this study was to quantify the potential payload loss due to snow and winter accessories and to predict the loss with weather variables. Tare weights of eight timber trucks were collected at mill receptions in Finland over a one-year period. Monthly and annual loss of potential payload was estimated using the tare measurements in summer months as a reference. Each load was also connected with weather data at the location and time of delivery and payload loss explained by the weather data with the aid of regression models. The maximum loss of payload varied between 1560 kg and 3100 kg. On a monthly basis, the highest losses occurred in January, when the median values varied between 760 kg and 2180 kg. Over the year, the payload loss ranged between the trucks from 0.5% to 1.5% (from 1.9% and 5.1% in January) of the total number of loads in the study. Payload loss was found to increase with decreasing temperature, increasing relative humidity and increasing precipitation. Although the average payload loss was not very high, the biggest losses occur just during the season of highest capacity utilization. Big differences were also found in the tare weights between the trucks. The results of the study give incentive to develop truck and trailer structures that reduce the adherence of snow.
Scarification is a mechanical site preparation technique designed to create microsites that will favor the growth of planted tree seedlings after clearcutting. However, the positive growth response of black spruce (Picea mariana (Mill.) Britton, Sterns & Poggenb.) to scarification varies across different sites. We hypothesized that this was due to different forms of physiological stress induced by different climates or by the severity of competition from ericaceous shrubs. We thus compared the effects of scarification on black spruce needle gas exchange and other foliar properties, as well as on indices of soil water and nitrogen availability, in relatively warm-dry (Abitibi) vs. cool-humid (Côte-Nord) climates in the province of Québec (Canada). We found a similar positive effect of scarification on tree height in Abitibi and Côte-Nord. Scarification reduced soil moisture in both climatic regions, but increased soil N mineralization in Côte-Nord only. Accordingly, scarification increased the instantaneous water use efficiency in both climate regions, but decreased photosynthetic N use efficiency in Côte-Nord only. In both regions, we found a positive relationship between foliar δ18O and δ13C on scarified plots, providing further evidence that increased growth due to scarification depends on a decrease in stomatal conductance. We conclude that scarification increases total evapotranspiration of trees evenly across the east-to-west moisture gradient in the province of Québec, but also improves long-term soil nutritional quality in a cooler-humid climate.
Ingrowth is an important element of stand dynamics in several silvicultural systems, especially in continuous cover forestry. Earlier predictive models for ingrowth in Finnish forests are few and not based on up-to-date statistical methods. Ingrowth is here defined as the number of trees over 1.3 m entering a plot. This study developed new ingrowth models for Scots pine (Pinus sylvestris L.), Norway spruce (Picea abies (L.) H. Karst.) and birch (Betula pendula Roth and B. pubescens Ehrh.) using data from the permanent sample plots of the Finnish national forest inventory. The data were over-dispersed compared to a Poisson process and had many zeros. Therefore, a zero-inflated negative binomial model was used. The total and species-specific stand basal areas, temperature sum and fertility class were used as predictors in the ingrowth models. Both fixed-effects and mixed-effects models were fitted. The mixed-effects model versions included random plot effects. The mixed-effects models had larger likelihoods but provided biased predictions. Also censored prediction was considered where only a certain maximum number of ingrowth trees were accepted for a plot. The models predicted most pine ingrowth in pine-dominated stands on sub-xeric and xeric sites where stand basal area was low. The predicted amount of spruce ingrowth was maximized when the basal area of spruce was 13 m2 ha–1. Increasing temperature sum increased spruce ingrowth. Predicted birch ingrowth decreased with increasing stand basal area and towards low fertility classes. An admixture of pine increased the predicted amount of spruce ingrowth.
Since the 1990’s, forest resource maps and small area estimates have been produced by combining national forest inventory (NFI) field plot data, optical satellite images and numerical map data using a non-parametric k-nearest neighbour method. In Finland, thematic maps of forest variables have been produced by the means of multi-source NFI (MS-NFI) for eight to ten times depending on the geographical area, but the resulting time series have not been systematically utilized. The objective of this study was to explore the possibilities of the time series for monitoring the key ecosystem condition indicators for forests. To this end, a contextual Mann-Kendall (CMK) test was applied to detect trends in time-series of two decades of thematic maps. The usefulness of the observed trends may depend both on the scale of the phenomena themselves and the uncertainties involved in the maps. Thus, several spatial scales were tested: the MS-NFI maps at 16 × 16 m2 pixel size and units of 240 × 240 m2, 1200 × 1200 m2 and 12 000 × 12 000 m2 aggregated from the MS-NFI map data. The CMK test detected areas of significant increasing trends of mean volume on both study sites and at various unit sizes except for the original thematic map pixel size. For other variables such as the mean volume of tree species groups, the proportion of broadleaved tree species and the stand age, significant trends were mostly found only for the largest unit size, 12 000 × 12 000 m2. The multiple testing corrections decreased the amount of significant p-values from the CMK test strongly. The study showed that significant trends can be detected enabling indicators of ecosystem services to be monitored from a time-series of satellite image-based thematic forest maps.
Despite the importance of spectral properties of woody tree structures, they are seldom represented in research related to forests, remote sensing, and reflectance modeling. This study presents a novel imaging multiangular measurement set-up that utilizes a mobile handheld hyperspectral camera (Specim IQ, 400–1000 nm), and can measure stem bark spectra in a controlled laboratory setting. We measured multiangular reflectance spectra of silver birch (Betula pendula Roth), Scots pine (Pinus sylvestris L.) and Norway spruce (Picea abies (L.) Karst.) stem bark, and demonstrated the potential of using bark spectra in identifying tree species using a Support Vector Machine (SVM) based approach. Intraspecific reflectance variability was the lowest in visible (400–700 nm), and the highest in near-infrared (700–1000 nm) wavelength regions. Interspecific variation was the largest in the red, red-edge and near-infrared spectral bands. Spatial variation of reflectance along the tree height and different sides of the stem (north and south) were found. Both birch and pine had increased reflectance in the forward-scattering directions for visible to near-infrared wavelength regions, whilst spruce displayed the same only for the visible wavelength region. In addition, spruce had increased reflectance in the backward-scattering directions. In spite of the intraspecific variations, SVM could identify tree species with 88.8% overall accuracy when using pixel-specific spectra, and with 97.2% overall accuracy when using mean spectra per image. Based on our results it is possible to identify common boreal tree species based on their stem bark spectra using images from mobile hyperspectral cameras.
Five Scots pine (Pinus sylvestris L.) progeny field trials, each established in different Lithuanian regions of provenance in 1983, were studied. Each progeny field trial consists of 140 half-sib families from seven populations (20 families from each population). The evaluation was carried out in 2012 and 2018 to assess the families resistance to Heterobasidion annosum (Fr.) Bref. An index of resistance in the infected plots was calculated. To verify the accuracy of the method, total phenolic compounds (TPC) was chosen as key parameter to compare with the plant resistance index. During the six years between the two assessments, the percentage of living Scots pine trees in the progeny field trials decreased up to 20 percentage points (range: 4 p.p. to 20 p.p.). In 2018 the area of H. annosum damaged plots (in percentage from total field trial area) varied from 17 to 27%. Tree mortality in the trial correlates with site soil fertility – more fertile soils were distinguished by higher tree percentage loss and vice versa. Using analysis from combined data of all progeny trials, the family variance component reached 13.3 ± 2.2% and family heritability was 0.81. Family heritability estimates for root rot resistance show possibilities of high breeding effectiveness. The correlations between the trials in family resistance estimates were negligible (ranging from 0 to 0.28). The significant high correlation coefficient was determined between the resistance index and TPC concentration (r = 0.77, p = 0.0003). This allows us to assume that plant resistance is directly linked on TPC synthesis. The results indicate that the chosen methods of chemical resistance for identification of root rot-resistant genotypes are applicable for the selection of Scots pine half-sib families in the field trials with higher resistance to pathogens.
The tree stem volume models of Norway spruce, Scots pine and silver and downy birch currently used in Finland are based on data collected during 1968–1972. These models include four different formulations of a volume model, with three different combinations of independent variables: 1) diameter at height of 1.3 m above ground (dbh), 2) dbh and tree height (h) and 3) dbh, h and upper diameter at height of 6 m (d6). In recent National Forest Inventories of Finland, a difference in the mean volume prediction between the models with and without the upper diameter as predictor has been observed. To analyze the causes of this difference, terrestrial laser scanning (TLS) was used to acquire a large dataset in Finland during 2017–2018. Field-measured predictors and volumes predicted using spline functions fitted to the TLS data were used to re-calibrate the current volume models. The trunk form is different in these two datasets. The form height is larger in the new data for all diameter classes, which indicates that the tree trunks are more slender than they used to be. One probable reason for this change is the increase in stand densities, which is at least partly due to changed forest management. In models with both dbh and h as predictors, the volume is smaller a given h class in the data new data than in the old data, and vice versa for the diameter classes. The differences between the old and new models were largest with pine and smallest with birch.
Moose (Alces alces L.) browsing causes severe damage in Scots pine (Pinus sylvestris L.) seedling stands. The effects of this damage on the quality of sawlogs were studied in a long-term controlled experiment. This article reports the stem size and external quality characteristics of Scots pine stems 34 years after artificial moose browsing damage. Damaging the trees by clipping the main stem at the seedling stage reduced the diameter, height, and tree volume of the trees at the end of the experiment. The tree growth reduction was dependent on the severity of clipping. The differences between the damaged and the control trees were more obvious in diameter than in height at the time of final felling. Stem form defects and vertical branches were the most typical externally detectable defects caused by clipping. Defects in the butt logs were detected in 71–89% of the damaged trees, depending on the clipping treatment severity. The stronger the clipping treatment, the more likely the stem form was defected and the more commonly were vertical branches and crooks detected in the stems. The results indicate that both tree dimensions and stem quality suffer from moose browsing. The findings of this controlled experiment more likely underestimate than overestimate the damage in comparison to real moose browsing. Further analyses are required to assess the effects of browsing damage on the internal quality of sawlogs and subsequent economic outcomes.
Using fine-resolution satellite imagery from multiple satellite data products, we assessed the change in forest cover of a state-managed Reserve Forest (RF) located in India’s Eastern Himalaya biodiversity hot-spot. 4.6% of forest cover was lost from Papum RF between 2013 and 2017 at the rate of 8.2 km2 year–1. Three species of hornbills: Great Hornbill Buceros bicornis Linnaeus, 1758, Wreathed Hornbill Rhyticeros undulatus (Shaw, 1811) and Oriental Pied Hornbill Anthracoceros albirostris (Shaw, 1808), that are functionally important are found here with nesting habitat in the areas affected by illegal logging. Therefore, we assessed the habitat loss within a 1 km radius around 29 nest trees. From 2011 to 2019, forest cover declined from 38.55 km2 to 21.94 km2 around these hornbill nest trees. Illegal logging is the main driver that is depleting forest cover within this important bird area. Our results highlight the ongoing threats to biologically-rich forests and the need for urgent measures to halt this loss. We suggest that this study has practical implications for the monitoring and governance of state-managed forests in Arunachal Pradesh.
The use of modern multi-functional forestry machines has already been associated with central nervous system fatigue induced by high mental workload. As these machines are being used under increasingly difficult terrain conditions, further knowledge is required on the expected aggravation of operators’ mental workload, so that suitable work/rest schedules can be developed. Within such a context, the aim of this study was to gauge aggravations of mental workload derived from increasing slope gradient. Measurements of eye activity were obtained from a representative harvester operator working in corridors with the following mean inclinations: 9%, 23% and 47%. The duration, frequency and trajectory of eye movements were used to determine the harvester operator’s mental workload, on the assumption that worsening work conditions would be reflected by increased eyeball activity. The number of fixations during the performance of all tasks increased with the increasing slope gradient. Similarly, fixation duration increased with slope gradient. The mean duration of saccades when working on a 23% slope was 5% shorter compared to work under a 9% gradient. A further significant shortening of saccade duration (~22%) occurred when working on a 47% slope. The good match between eye activity cycles and work cycles, visible especially on steep slopes, indicates that mental workload is related to work conditions. Overall, operating a forest harvester on steep slopes results in a greatly increased mental workload and calls for suitable rest schedules.
In northern Europe, largest part of non-wood forest products (NWFPs) are gathered for recreational purposes and household consumption, but considerable amount of forest berries and mushrooms are sold as well. Retail market, largely invisible for the official statistics, reveals the lifestyle-related aspects of NWFP trade and may help to understand the flows of this ecosystem service when information on wholesale trade is inaccessible. The prices and flows of most common NWFPs – edible berries, mushrooms and tree sap – in the retail market in Latvia in 2017 and 2018 were analysed based on direct interviews with the sellers in marketplaces and telephone interviews with online retailers. The mean retail prices of NWFPs were compared between statistical regions and years and correlated with socio-economic data and forest characteristics. The directions of the NWFP flows were analysed according to the place of origin and place of retail sales. The highest prices were recorded for stinkhorn (Phallus impudicus Pers.) and Boletes spp. among mushrooms, for wild strawberries (Fragaria vesca L.) among berries and for maple (Acer platanoides L.) sap in the product group of tree sap. The retail price of the same products differed between years, most likely due to the product availability, largely caused by meteorological conditions. In more than half of the cases of recorded sales, NWFPs were consumed in the same region as they were gathered. For other cases of sales, the capital, Rīga, was the main service benefitting area of NWFP retail trade, and the largest part of the products originated from the two closest statistical regions.
To preliminary evaluate the potential wood utilization of Betula platyphylla Sukaczev trees naturally regenerated in Mongolia, growth characteristics (stem diameter and tree height), wood properties (annual ring width, basic density, and compressive strength parallel to grain at the green condition) of core samples, and stress-wave velocity in stems were investigated for Betula platyphylla trees grown naturally in three different sites in Selenge, Mongolia. Betula platyphylla trees, naturally grown in Nikko, Japan, were also examined to compare wood properties between the two regions. The mean values of stem diameter, tree height, stress-wave velocity of stems, annual ring width, basic density, and compressive strength parallel to grain at green condition in Mongolian B. platyphylla were 17.6 cm, 14.1 m, 3.50 km s–1, 1.27 mm, 0.51 g cm–3, and 20.4 MPa, respectively. Basic density and compressive strength were decreased first from the pith, and then gradually increased toward the bark. The wood properties of B. platyphylla trees grown naturally in Mongolia were similar to those in B. platyphylla trees grown in Japan. Growth characteristics, especially stem diameter, were positively correlated with the stress-wave velocity of stems and basic density. Early evaluation of basic density in B. platyphylla trees is possible by using wood located 2 cm from the pith. Basic density at the position from the 1st to the 15th annual ring from the pith showed significant between-site differences in Mongolian B. platyphylla. Based on the results, it is concluded that the wood of B. platyphylla trees grown in Mongolia may be used for industrial products as well as those from similar species in other countries.
A major after-use option for former peat harvesting areas has been afforestation. The profitability of afforestation with Scots pine trees (Pinus sylvestris L.) was studied in two 31–32-year old experiments in southern and northern Finland. The stands were established by seeding and planting, and various fertilization treatments and drainage intensities were tested. The financial performance for each plot was assessed in three steps. First, the costs occurred during the measurement time were summed up according to their present value. Then, for the rest of the rotation (i.e., from the age of 31/32 onwards) the stand management was optimized in order to maximize the net present value (MaxNPV). Finally, bare land values (BLVs) were calculated by summing up the present value of costs and the MaxNPV and converting the sum of the series into infinity. The afforestation method did not affect the mean annual increment (MAI; 9.2–9.5 m3 ha–1 a–1) in the southern experiment. In the northern experiment the afforestation method, ditch spacing and fertilization had significant effects on the MAI of the stands. The average MAI of the planted pines was 8.9 m3 ha–1 a–1, and for seeded pines it was 7.5 m3 ha–1 a–1. The BLV at an interest rate of 3% was positive for all stands in both regions. In the northern region afforestation method, ditch spacing and fertilization also had a significant effect on the BLV. When the interest rate was 5%, almost two thirds of the stands had a negative BLV in both regions.
Although crop tree release (CTR) in hardwood stands is an accepted variant of precommercial thinning (PCT), the lack of an affordable and feasible method hinders its adoption. CTR implies selecting between 150 and 500 trees ha–1 when trees are between 7 and 12 m high and cutting only stems competing with the target crop trees. We performed a field trial of a CTR variant of PCT in a 27.8 ha hardwood stand using a backpack mounted chain saw. A detailed time study was performed to document the trial over 13 days. Compared to conventional PCT performed earlier in the life of a stand, precommercial crop tree release required cutting larger stems, which showed to be feasible and productive using a backpack mounted chain saw. Productivity varied between 0.22 to 0.47 ha h–1 during the trial, Although productivity could vary with stand characteristics and worker, this proof of concept trial demonstrates some of the potential uses that this new saw configuration offers and sets the basis for an eventual larger scale deployment of this treatment.
Leaf reflectance and transmittance spectra are essential information in many applications such as developing remote sensing methods, computing shortwave energy balance (albedo) of forest canopies, and monitoring health or stress of trees. Measurement of coniferous needle spectra has usually been carried out with single integrating spheres, which has involved a lot of tedious manual work. A small double integrating sphere would make the measurements considerably faster, because of its ease of operation and small sample sizes required. Here we applied a compact double integrating sphere setup, used previously for measurement of broad leaves, for measurement of coniferous needles. Test measurements with the double integrating sphere showed relative underestimation of needle albedo by 5–39% compared to a well-established single integrating sphere setup. A small part of the bias can be explained by the bias of the single sphere. Yet the observed bias is quite significant if absolute accuracy of measurements is required. For relative measurements, e.g. for monitoring development of needle spectra over time, the double sphere system provides notable improvement. Furthermore, it might be possible to reduce the bias by building an optimized measurement setup that minimizes absorption losses in the sample port. Our study indicates that double spheres, after some technical improvement, may provide a new and fast way to collect extensive spectral libraries of tree species.
The effects of wood ash fertilisation on tree nutrition and growth on forested peatlands has been studied using loose ash, but in practice, ash fertilisation is done almost exclusively with granulated ash. In this study, the effects of granulated ash and loose ash (both 5 Mg ha–1) on the growth and nutrition of Scots pine (Pinus sylvestris L.) stands were compared between a nitrogen-poor and a nitrogen-rich site over 15 years. On the nitrogen-rich site, wood ash application was also compared with commercial PK fertilisation. On the nitrogen-rich site, mean stand volume growth increase over unfertilised control treatment during the 15 year study period using granulated ash and commercial PK fertiliser was of the same magnitude (on average, 2.2–2.3 m3 ha–1 a–1). However, when loose ash was used growth increase over control was higher (3.7 m3 ha–1 a–1). On the nitrogen-poor site, the mean growth increase gained by loose or granulated ash (1.4–1.5 m3 ha–1 a–1) over the unfertilised control treatment was not significant. Fertilisation with loose ash or PK increased foliar P, K and B concentrations already in the first or second growing season, following fertilisation on both sites. Granulated ash increased foliar P concentrations on the nitrogen-rich site less than loose ash. After an initial increase, foliar P, K and B concentrations decreased at the end of study period. On the nitrogen-poor site, foliar P concentrations were below the deficiency limit by the end of the study period.
Growing demand for wood products, combined with efforts to conserve natural forests, have supported a steady increase in the global extent of planted forests. Here, a two-phase sampling strategy for large-scale assessment of the total area and the total wood volume of fast-growing forest tree crops within agricultural land is presented. The first phase is performed using tessellation stratified sampling on high-resolution remotely sensed imagery and is sufficient for estimating the total area of plantations by means of a Monte Carlo integration estimator. The second phase is performed using stratified sampling of the plantations selected in the first phase and is aimed at estimating total wood volume by means of an approximation of the first-phase Horvitz-Thompson estimator. Vegetation indices from Sentinel-2 are exploited as freely available auxiliary information in a linear regression estimator to improve the design-based precision of the estimator based on the sole sample data. Estimators of the totals and of the design-based variances of total estimators are presented. A simulation study is developed in order to check the design-based performance of the two alternative estimators under several artificial distributions supposed for poplar plantations (random, clustered, spatially trended). An application in Northern Italy is also reported. The regression estimator turns out to be invariably better than that based on the sole sample information. Possible integrations of the proposed sampling scheme with conventional national forest inventories adopting tessellation stratified sampling in the first phase are discussed.
Genetic parameters of growth and stem quality traits were estimated for open-pollinated silver birch Betula pendula Roth progenies in Latvia at the age of 10 and 14 years. Tree height and stem volume were found to be under strong genetic control at both inventories (narrow-sense heritabilities varied from 0.41 to 0.66). Mainly low heritabilities were found for stem defects, yet genetic control of branch diameter, stem straightness and overall stem quality varied from low to high depending on study site. High additive genetic coefficient of variation was found for stem volume (25.3–32.5%). Genetic correlations among growth traits were strong and positive (0.90–0.99). Mainly weak genetic correlations between growth and quality traits implied simultaneous improvement. Still, strong negative correlations between branch angle and stem straightness might result in enlarged knot size for straighter logs. The genetic age-age correlations were strong. Weak genotype by environment interaction and stability of best genotypes over different sites was indicated by strong genetic correlations between trials. Each growth or quality trait alone showed substantial improvement in terms of estimated genetic gain (up to 62% over trial mean for stem volume). Therefore, selection index combining both growth and stem quality may be developed.
Cable yarding is a general solution for load handling on sites not accessible to ground-based machinery, and is typically associated with steep terrain. On flat terrain, such conditions can primarily be found on soft or wet soils, most frequently encountered in Central and Northern European countries. Today, changed environmental and market conditions may offer an unprecedented opportunity to the actual implementation of cable yarding on flat terrain in commercial operations. The study goal was to collect cable yarder manufacturers experience regarding the use and adaption of cable yarding technology on flat terrain. European manufacturers of cable yarding technology were interviewed about customer experience, particular challenges, adaptation potential, future potential and main hurdles for the expansion of cable yarding on flat terrain. Almost all manufacturers have received requests for flat-terrain yarding technology solutions, primarily from Germany. Temporal or permanent inaccessibility, regulatory or environmental reasons were the most frequent motivation for considering cable yarding technology. Installation was considered particularly challenging (clearance, stable anchoring). Potential adaptations included higher towers, artificial anchors, mechanized bunching before extraction and un-guyed yarder-systems. An artificial, highly mobile, self-anchoring tail spar was considered the most useful adaptation. While concerned about limited profitability and qualified labour shortage, most manufacturers demonstrated a positive or neutral view concerning the expansion of cable yarding on flat terrain. However, cable yarding is not considered to be cost-competitive wherever ground-based systems can be employed and cable yarding is not subsidized.
In the forest areas of eastern China, there is a change from forest dominated by deciduous broad-leaved trees to forest dominated by evergreen broad-leaved trees as the latitude or altitude decreases. Different life forms have different survival strategies to deal with climate change, and studying the life form dynamics of the tree layers in the mixed forest in eastern China, with increasing temperature, can help us understand how the forest responds. This study was performed in a 1 ha plot in evergreen and deciduous broad-leaved mixed forest in Tianmu Mountain National Nature Reserve. Based on the data from two surveys (1996 and 2017), the changes in life form composition and biodiversity over the past 21 years were analyzed. We obtained the following results: (1) The proportion of evergreen trees increased from 55.0% in 1996 to 67.5% in 2017, and the dominance of evergreen species was enhanced. (2) The diversity of both life forms increased, and the tree species were more abundant. (3) The average annual recruitment rate of the evergreen species was 2.1% greater than their mortality rate, and the average annual recruitment rate of the deciduous species was 0.5% less than their mortality rate. (4) The competition among the trees in the small-diameter class (10 cm ≤ DBH < 20 cm) was fierce for many tree species. The proportion of the evergreen species in the small-diameter class was high. The life forms making up the mixed climax forest community has changed over the past 21 years, with the proportion and dominance of evergreen trees increasing significantly.
This study proposes an original method for tree species classification by satellite remote sensing. The method uses multitemporal multispectral (Landsat OLI) and hyperspectral (Resurs-P) data acquired from determined vegetation periods. The method is based on an original database of spectral features taking into account seasonal variations of tree species spectra. Changes in the spectral signatures of forest classes are analyzed and new spectral–temporal features are created for the classification. Study sites are located in the Czech Republic and northwest (NW) Russia. The differences in spectral reflectance between tree species are shown as statistically significant in the sub-seasons of spring, first half of summer, and main autumn for both study sites. Most of the errors are related to the classification of deciduous species and misclassification of birch as pine (NW Russia site), pine as mixture of pine and spruce, and pine as mixture of spruce and beech (Czech site). Forest species are mapped with accuracy as high as 80% (NW Russia site) and 81% (Czech site). The classification using multitemporal multispectral data has a kappa coefficient 1.7 times higher than does that of classification using a single multispectral image and 1.3 times greater than that of the classification using single hyperspectral images. Potentially, classification accuracy can be improved by the method when applying multitemporal satellite hyperspectral data, such as in using new, near-future products EnMap and/or HyspIRI with high revisit time.
In Nordic countries, tree planting of seedlings is mainly performed during spring and early summer. Interest has increased in extending the planting window throughout the unfrozen growing season. This study compared the success of one-year-old spring, summer and autumn plantings in practical forestry in Norway spruce (Picea abies (L.) Karst.) and Scots pine (Pinus sylvestris L.) in southern and central Finland. Planting success was based on the number of viable seedlings per hectare relative to a species-specific target density. The influence of different factors to poor planting results were determined, including quality of site preparation and planting, and sources of natural damage. Overall, in Norway spruce, 85, 69 and 84% and in Scots pine 53, 55 and 40% of spring, summer and autumn plantings succeeded. In Norway spruce, the planting results were consistent between the southern and central regions, whereas in Scots pine, the success was slightly lower in the south. The poor work quality and a low density of appropriate planting spots, contributed to poor planting results, regardless of planting season, region or tree species. Considering different damages, especially mammal damage contributed to the failure of Scots pine spring plantings, whereas in summer plantings, corresponding single failure reason could not be identified. Based on our findings, extending the planting season of Norway spruce could be recommended in both regions. For Scots pine, there is still significant uncertainty about the success of summer and autumn plantings, partially due to the limited number of plantings available for analyses.
Diameter at breast height (DBH) and height (H) of trees are two important variables used in forest management plans. However, collecting the measurements of H is time-consuming and costly. Instead, the H-DBH relationship is modeled and used to estimate H. But, ignoring the effects of slope, aspect and tree competition on the H-DBH relationship often impedes the improvement of H predictions. In this study, to improve predictions of Cyclobalanopsis glauca (Thunb.) Oerst. tree H in mixed forests, we compared eleven H-DBH models and examined the influence of slope and aspect on the H-DBH relationship using 426 trees. We then improved Hegyi competition index and explored its effect on the H predictions by including it in the selected models. Results showed 1) There were statistically significant effects of slope and aspect on the H-DBH relationship; 2) The log transformation and exponential model performed best for sunny- and shady-steep, respectively, and the Gompertz’s model was optimal for both sunny- and shady-gentle; 3) Compared with the whole dataset, the division of the data into the slope and aspect sub-datasets significantly reduced the RMSE of H predictions; 4) Compared with the selected models without competition index, adding the original Hegyi and improved Hegyi_I into the models improved the H predictions but only the models containing the improved Hegyi_I significantly increased the prediction accuracy at the significant level of 0.1. This study implied that modeling the H-DBH relationship under different slopes and aspects and including the improved Hegyi_I provided the great potential to improve the H predictions.
Accurate estimates of aboveground biomass (AGB) strongly depend on the suitability and precision of allometric models. Diospyros mespiliformis Hochst. ex A. DC. is a key component of most sub-Sahara agroforestry systems and, one of the most economically important trees in Africa. Despite its importance, very few scientific information exists regarding its biomass and carbon storage potential. In this study direct method was used to develop site-specific biomass models for D. mespiliformis tree components in Burkina Faso. Allometric models were developed for stem, branch and leaf biomass using data from 39 tree harvested in Sudanian savannas of Burkina Faso. Diameter at breast height (DBH), tree height, crown diameter (CD) and basal diameter (D20) were regressed on biomass component using non-linear models with DBH alone, and DBH in combination with height and/or CD as predictor variables. Carbon content was estimated for each tree component using the ash method. Allometric models differed between the experimental sites, except for branch biomass models. Site-specific models developed in this study exhibited good model fit and performance, with explained variance of 81–98%. Using models developed from other areas would have underestimated or overestimated biomass by between –72% and +98%. Carbon content in aboveground components of D. mespiliformis in Tiogo, Boulon and Tapoa-Boopo was 55.40% ± 1.50, 55.52% ± 1.06 and 55.63% ± 1.00, respectively, and did not vary significantly (P-value = 0.909). Site-specific models developed in this study are useful tool for estimating carbon stocks and can be used to accurately estimate tree components biomass in vegetation growing under similar conditions.
The time consumption (TC) of pre-commercial thinning (PCT) varies greatly among sites, stands and forest workers. The TC in PCT is usually estimated by field-assessed work difficulty factors. In this study, a linear mixed model for the TC in PCT was prepared by utilizing forest resources data (FRD). The modelling data included 11 848 and validation data included 3035 worksites with TC information recorded by forest workers within the period of 2008–2018. The worksites represented a range of site and stand conditions across a broad geographical area in Finland. Site and stand characteristics and previous management logically explained the TC in PCT. The more fertile the site, the more working time was needed in PCT. On sites of medium fertility, TC in the initial PCT increased with stand age by 0.5 h ha–1 yr–1. Site wetness increased the TC. PCT in summer was more time consuming than in spring. Small areas were more time consuming to PCT per hectare than larger ones. The between-forest worker variation involved in the TC was as high as 35% of the variation unexplained by the TC model. The coefficient of determination in validation data was 19.3%, RMSE 4.75 h ha–1 and bias –1.6%. The TC model based on FRD was slightly less precise than the one based on field-assessed work difficulty factors (removal quantity and type and terrain difficulty): RMSE 4.9 h ha–1 vs. 4.1 h ha–1 (52% vs. 43%). The TC model could be connected to forest information systems where it would facilitate the predictions of the labour costs of PCT without field-assessing work difficulty factors.
We studied variations on different traits of Parkia timoriana (D.C.) Merr. in twelve provenances systematically from their source of origin to a controlled environment where representative seedlings were grown. Among the provenances, P1 gave the best result for seed traits including germination traits, P7 for pod traits and P10 for seedling vigour. Effects of seasonal distribution of rainfall and temperature on seed and pod traits were also determined by computing multiple regression analysis. The results displayed winter rainfall and summer temperature as the most important factor determining pod and seed traits. Latitude also significantly (P < 0.001) affected PWT (r = 0.52), SWP (r = 0.46) and SW (r = 0.50). A common garden study for germination and seedling growth indicated P1 and P10 provenance as the best among all. Seeds drawn from P10 gave the highest seedling vigour with an average growth rate of 0.61 cm/day from 90th to 180th day. Highest broad-sense heritability values (h2) were observed in germination traits, followed by seedling collar diameter. The lowest h2 was observed for seedling height.
We studied the spatial decomposition rates of standardised organic substrates in soils (burned boreal pine-dominated sub-xeric forests in eastern Finland), with respect to charred and non-charred coarse woody debris (CWD). Decomposition rates of rooibos plant litter inside teabags (C:N = 42.870 ± 1.841) and pressed-sheet Nordic hardwood pulp (consisting of mainly alpha-cellulose) were measured at 0.2 m distance from 20 charred (LC0.2) and 40 non-charred logs (LNC0.2). We also measured decomposition at 60 plots located 3–10 m away from downed logs (L3,10). The rooibos decomposition rate constant ‘k’ was 8.4% greater at the LNC0.2 logs than at the L3,10 or LC0.2 logs. Cellulose decomposed more completely in 1 micron mesh bags at LNC0.2 (44% of buried bags had leftover material) than at LC0.2 (76%) or L3,10 (70%). Decomposition of cellulose material was rapid but varied greatly between sampling plots. Our results indicate that decomposition of the standardised organic matter was more rapid close to CWD pieces than further away. However, only the plots located near non-charred logs (LNC0.2) exhibited high decomposition rates, with no corresponding increase observed at the charred logs (LC0.2). This suggests a possible noteworthy indirect effect of forest burning on soil organic matter (SOM) decomposition rates close to charred CWD after forest fires. We urge for more studies on this tentative observation as it may affect the estimates on how fires affect carbon cycling in forests.
Nature-based tourism (NBT) is a growing industry within regions rich in natural amenities worldwide. An important feature of NBT business is the dependence on the quality of surrounding environment. This paper addresses the role of the management of commercial forests owned by the state in Finnish Lapland. The paper explores the NBT entrepreneurs’ willingness to participate in a proposed new landscape and recreational value trading (LRVT) and elaborates the effect of entrepreneur and enterprise characteristics, such as entrepreneurial attitude, venture size, and a variety of services offered to customers, on the experienced and expected growth of NBT enterprise. The survey data on NBT enterprises were analyzed with ordered and binary logit models. The willingness of enterprises to participate in LRVT depended on the venture size, entrepreneurial attitude, and type of activities offered to customers. The results show that relatively young and small-sized enterprises have faced difficulties in developing their business. Entrepreneurial experience, risk-taking and intention to develop new business associate positively with expected increase in turnover.
Trees are particularly susceptible to climate change due to their long lives and slow dispersal. However, trees can adjust the timing of their growing season in response to weather conditions without evolutionary change or long-distance migration. This makes understanding phenological cueing mechanisms a critical task to forecast climate change impacts on forests. Because of slow data accumulation, unconventional and repurposed information is valuable in the study of phenology. Here, I develop and use a framework to interpret what phenological patterns among provenances of a species in a common garden reveal about their leafing cues, and potential climate change responses. Species whose high elevation/latitude provenances leaf first likely have little chilling requirement, or for latitude gradients only, a critical photoperiod cue met relatively early in the season. Species with low latitude/elevation origins leafing first have stronger controls against premature leafing; I argue that these species are likely less phenologically flexible in responding to climate change. Among published studies, the low to high order is predominant among frost-sensitive ring-porous species. Narrow-xylemed species show nearly all possible patterns, sometimes with strong contrasts even within genera for both conifers and angiosperms. Some also show complex patterns, indicating multiple mechanisms at work, and a few are largely undifferentiated across broad latitude gradients, suggesting phenotypic plasticity to a warmer climate. These results provide valuable evidence on which temperate and boreal tree species are most likely to adjust in place to climate change, and provide a framework for interpreting historic or newly-planted common garden studies of phenology.
In the Nordic countries Finland, Norway and Sweden, the most common regeneration method is planting after clearcutting and, often, mechanical site preparation (MSP). The main focus of this study is to review quantitative effects that have been reported for the five main MSP methods in terms of survival and growth of manually planted coniferous seedlings of Norway spruce (Picea abies (L.) Karst.), Scots pine (Pinus sylvestris L.) and lodgepole pine (Pinus contorta var. latifolia Engelm.) in clearcuts in these three countries. Meta analyses are used to compare the effects of MSP methods to control areas where there was no MSP and identify any relationships with temperature sum and number of years after planting. In addition, the area of disturbed soil surface and the emergence of naturally regenerated seedlings are evaluated. The MSP methods considered are patch scarification, disc trenching, mounding, soil inversion and ploughing. Studies performed at sites with predominately mineral soils (with an organic topsoil no thicker than 0.30 m), in boreal, nemo-boreal and nemoral vegetation zones in the three Fenno-Scandinavian countries are included in the review. Data from 26 experimental and five survey studies in total were compiled and evaluated. The results show that survival rates of planted conifers at sites where seedlings are not strongly affected by pine weevil (Hylobius abietis L.) are generally 80–90% after MSP, and 15–20 percent units higher than after planting in non-prepared sites. The experimental data indicated that soil inversion and potentially ploughing (few studies) give marginally greater rates than the other methods in this respect. The effects of MSP on survival seem to be independent of the temperature sum. Below 800 degree days, however, the reported survival rates are more variable. MSP generally results in trees 10–25% taller 10–15 years after planting compared to no MSP. The strength of the growth effect appears to be inversely related to the temperature sum. The compiled data may assist in the design, evaluation and comparison of possible regeneration chains, i.e. analyses of the efficiency and cost-effectiveness of multiple combinations of reforestation measures.
The research was carried out in unmanaged middle-aged (75–85 years) Northern taiga Scots pine (Pinus sylvestris L.) forests in the Kola peninsula. It was established that forests of green moss-lichen and green moss site types are characterised by a predominance (65–70% by stand volume) of moderately and strongly weakened trees. Trees of differing vitality have significant differences in annual increment. Healthy trees had a radial increment (RI) 70–75% greater than that of dying trees, and a basal area increment (BAI) 85–90% greater. The dynamics of the RI and BAI of Scots pine trees for the 70-year period (from 1945 to 2015) is different. The RI of all individuals in the communities studied decreases consistently. The decrease is expressed more strongly in green moss Scots pine forests (80–95% from 1945 to 2015) compared to green moss-lichen forests (60–80%); it manifests itself more in strongly weakened and dying individuals (75–95%) than in healthy and moderately weakened ones (60–80%). Annual basal area increment in green moss Scots pine forests increases by 45–65% from stand establishment until the trees are 25 to 35 years old and subsequently decreases by 50–80% to 70–80 years of age. In green moss-lichen pine forests the BAI of Scots pine remains rather stable in healthy and moderately weakened trees and decreases in strongly weakened and dying individuals by 45% and 75–80%, respectively throughout the studied period.
Forests are affected by climate change in various ways. This includes abiotic factors such as droughts, but also biotic damage by pest insects. There are numerous examples from cases where pest insects have benefitted from longer growing seasons or from warmer summers. Similarly, new pest insects have been able to expand their range due to climatic conditions that have changed from hostile to tolerable. Such seems to be the case with the nun moth (Lymantria monacha), an important defoliator of coniferous trees in Europe. For centuries, the species has had massive outbreaks across Central-Europe, while it has been a rare inhabitant in Northern Europe. Recently, the nun moth population in Finland has not only expanded in range, but also grown more abundant. This research note describes the results from the first years (2018–2019) of a monitoring program that is being conducted with pheromone traps across central and southern Finland. So far, the northernmost individuals were trapped near the 64 N degrees. However, there were more southern locations where no moths were trapped. The species was present in every trapping site below the latitude of 62 N degrees. More importantly, at some sites the abundance of the nun moth suggested that local forest damage may already occur. Given the current climatic scenarios for Fennoscandia, it is likely that the nun moth populations will continue to grow, which is why systematic surveys on their abundance and range expansions will be topical.
One of the main objectives of the International Union of Forest Research Organization’s WP 9.01.06 (Forest Science Publishing) is to clarify the range of forestry-related journals in the context of the overall objective of supporting authors and readers to deepen their knowledge of forestry. This study extended an earlier list of forestry journals using the ISSN database and the CAB abstracts. The study included 451 journal titles, each categorized by ISSN, publisher, language of publication, website, and geographical coverage (i.e. national, regional, international, or three mixed intermediate categories), as well as information on indexing in Web of Science, Scopus, DOAJ, and SherpaRomeo. The included journals are published in 61 countries and in 32 different languages. Those categorized as international are mostly published in English. 21.7% of the journals are indexed in Web of Science and 34.1% in Scopus. 95.6% of the titles are published online or in the print+online publishing model, but only 57.0% of the titles are published in Open Access, of which only 33.7% are indexed in DOAJ.