Are small protected habitat patches within boreal production forests effective in conserving species richness, abundance and community composition? A systematic review

In boreal zone forest management is changing and degrading forest habitats, which has caused declines in biodiversity. To mitigate these harmful effects in production forests, small-scale habitats with high biodiversity values have been protected within them. These habitats include woodland key habitats, and other small habitat patches protected by voluntary conservation actions. In this systematic review we synthesize the evidence on the value of small protected habitat patches (SPHP) within production forest landscapes for biodiversity. Review question: Are small protected habitat patches within boreal production forests effective in conserving species richness, abundance, and community composition? Both peer-reviewed and grey literature were searched from bibliographical databases, organizational websites and internet search engines in English, Finnish, Swedish and Russian. Articles were screened at two stages (title/abstract and full text) and the validity of the included studies were assessed. Screening and validity assessment were based on predetermined criteria. After data extraction, narrative and quantitative syntheses were conducted. Influences of effect modifiers were tested, and sensitivity analyses were conducted. During the searches 19,458 articles were found. After duplicate removal and title/abstract screening 336 articles remained. During full text screening 41 articles were included and 35 of them (174 studies) were included in narrative synthesis. 28 articles with 127 studies had suitable data for meta-analysis. SPHPs had significantly higher species richness compared to production forests. When compared to natural forests, there was no significant difference. Forest management in areas surrounding SPHPs did not have impact on species richness of these patches. Individual abundance was significantly higher in SPHPs compared to natural or production forests. There was significantly more dead wood in SPHPs compared to production forests, but when compared to natural forests there was no significant difference. Community composition was different between SPHPs and both production and natural forests. The findings of this review show that small protected patches within production forests are important part of biodiversity conservation. They cannot substitute larger protected areas but supplement the protected area network. However, there were gaps both in geographical distribution of the studies as well as in the selection of target species of the studies. Therefore, generalization of the results must be done carefully.

Background Protecting areas from human usage has traditionally been the main method for biodiversity conservation. However, setting aside areas in productive forest land is expensive and only a relatively small proportion of forests is protected for biodiversity in Northern Europe [1] Globally, the network of protected areas is considered inadequate to maintain species assemblages [2][3][4][5]. Protected areas are often situated in less productive areas with low biodiversity [6,7], and establishing new protected areas is not possible in the magnitude needed for maintaining biodiversity [5]. Therefore, conserving biodiversity in production forests has become increasingly important.
Practices to maintain and increase biodiversity in production forests include green tree retention, prescribed burning, leaving dead wood in forests and creating habitat corridors and buffer strips [7,8]. In addition, preservation of small patches of certain habitats, e.g. woodland key habitats (WKHs), was introduced as a method in biodiversity conservation in early 1990′s [9]. Woodland key habitats are a common concept in Northern Europe first introduced in Sweden [10,11]. Nowadays the concept is used also in Finland, Norway, the Baltic countries, and Russia (mostly in the northern regions of the European part of the country and Siberia [12]) There are differences in the definitions and legal status of the WKHs and their delineation vary between countries [9][10][11][12][13][14][15][16][17][18][19][20]. For example, in Sweden, WKHs are defined by certain structural properties or species present in them and the protection of WKHs is mostly voluntary while in Finland WKHs are protected by the Forest Act [10,14].
In addition to WKHs, other new kinds of conservation policies have been developed. Finnish Forest Biodiversity Program for Southern Finland (METSO) [21], Norwegian Frivillig vern [22] and Swedish Komet programmet [23], are all based on the forest owners´ own initiative to protect their forests. Forest owners offer their forest to be protected by the program, and if the forest has sufficient ecological values, e.g. high volume of dead wood or quantity of large deciduous trees, the owners will get a compensation payment.
In the USA the concept of conservation easements is the main way of voluntary forest conservation. Conservation easements mean that landowners give up their right to develop the forest and they get monetary or taxationbased compensation from the government or a conservation group (land trust) [24]. In Canada and Russia, almost all forests are publicly owned, and therefore comparable voluntary conservation systems do not exist. In the Baltic countries the private forest ownership has been re-established after regaining their independence in 1991 [25]. Environmental values of forests are emphasized both by national forest programs and private forest owners, but most conservation programs are still determined by state authorities [25,26].
Forest management certification systems may also include provisions for conserving certain habitats beyond legal requirements. Both Forest Stewardship Council (FSC) and Program for the Endorsement of Forest Certification (PEFC) systems are widely used in the boreal zone [27,28]. To get certified, forest owners must commit themselves to responsible forest management. National certification standards vary between countries, but they often include preservation of small habitat patches of high conservation value [29][30][31].
But is conserving small habitat patches within production forests an effective way to maintain biodiversity? There has been critique of the small size and scattered distribution of WKHs [32]. It has been suggested that isolated woodland key habitats suffer from extinction debt [33] and that small habitat patches may not be able to maintain species diversity over time [34]. In previous systematic review, WKHs were found to be high in species richness and in the number of red-listed species but their ability to maintain diversity in managed landscape could not be addressed [35]. As the evidence base has not been reviewed recently, it is time for an update. Here we present a systematic review on the importance of small protected habitat patches (SPHPs hereafter) within production forests for maintaining biodiversity. The topic arose from the interest of Finnish forest industry on the effectiveness of conservation actions within production forests to produce beneficial biodiversity outcomes. A stakeholder workshop was held in November 2018 to further discuss the topic. There was a balanced representation of different stakeholder interest groups from government agencies and academia to non-governmental organizations and the private sector. In the workshop review questions initial theory of change ( Fig. 1) and PICO-based search terms were defined and factors potentially creating heterogeneity discussed. Further details of the stakeholder engagement are provided in the protocol [36].

Objective of the review
The objective of this article is to systematically review and synthesize the impacts of conservation of small habitat patches within production forests on biodiversity. Fundamentally, biodiversity is a complex measure at genetics, species and ecosystem level, but here we concentrate on species level diversity, yet still calling it biodiversity. Legally designated woodland key habitats as well as any small-scale voluntary conservation areas are included. The review focuses on terrestrial biodiversity in boreal forests. We follow the definition of boreal zone by Keenan et al. [37] but will also include Baltic countries. The geographical scope was determined at the stakeholder meeting.
Review question: Are small protected habitat patches within boreal production forests effective in conserving species richness, abundance, and community composition?
The question components are: P: Boreal forests (Canada, USA, Norway, Sweden, Finland, Russia, Estonia, Latvia, Lithuania). I: Small habitat patches set aside for conservation within production forests (SPHPs).
C: Managed forests, natural forests. O: Species diversity, abundance, assemblage, dead wood volume and diversity.

Methods
This systematic review follows the guidelines of Collaboration for Environmental Evidence and complies with the ROSES reporting standards. The ROSES form is included as an Additional file 1. All the deviations from the a priori systematic review protocol [36] are reported and explained below.

Deviations from the protocol Search string
After the publication of our systematic review protocol [36] we noticed some mistakes in the search string and few changes were made. The retention term was  1 Simplified theory of change on the contribution of different actions to maintain or increase biodiversity in production forests to counterbalance the effects of forest harvesting mistakenly included in the search string when the protocol was prepared for publication after the peer review.

Internet searches
In internet searches, one addition to the search method was made compared to the protocol [36], as no citations or patents were considered in Google Scholar.

Validity assessment
Because of limited time resources, all studies were assessed by one person and their decision checked by another person resulting in a small deviation of the protocol [36] as double-blind assessment was not used.

Potential effect modifiers and sources of heterogeneity
Two effect modifiers were added to the list in the protocol [36]: altitude and differences in the forest management surrounding SPHP, and one was removed: size and disposition of retention trees because this information was merged with another effect modifier: differences in forest management (comparator). These effect modifiers were added during the data extraction when differences concerning these potential effect modifiers were noted between included studies.

Searching for articles
Articles were searched from 1990 onwards because small scale habitat protection within production forests was integrated in production forestry in the 1990s.

Search terms
Based on the discussions at the stakeholder meeting, a search string was formulated. The performance of the search string was tested in the Web of Science and Scopus using a test list of 20 articles collected from previous reviews and from experts. The search string used for testing in the protocol [36] mistakenly included term TS = (retention*) and, also the test list in the protocol mistakenly included three articles dealing with retention as originally a broader scope was proposed for the review. These mistakes were caused by human error when modifying the manuscript after peer review. Since retention was not included in the final manuscript, the three articles about retention were removed from the test list and the performance of the final (corrected) search string was re-tested. The corrected test list is in Additional file 2 and the results from testing of the search string in Additional file 3. The final search string (Web of Science format) was: #1 TS = ((Boreal NEAR/5 (forest* OR zone OR tree*)) OR taiga OR spruce* OR picea OR pine* OR pinus OR birch* OR aspen* OR populus) #2 TS = (Finland OR Finnish OR Swed* OR Norw* OR Russia* OR Estonia* OR Latvia* OR Lithuania* OR Fennoscan* OR Scandin* OR Baltic OR "North* Europ*" OR Canad* OR "North* Ameri*" OR Siber* OR Alaska OR "United States" OR USA) and TS = (forest* OR tree*) #3 = #1 OR #2 #4 TS = ("key habitat*") OR TS = ("forest act habitat*") OR TS = (reserve* NEAR/5 (forest* OR OR privat* OR area* OR patch* OR habitat*)) OR TS = ("private* protected area*") OR TS = (voluntar* NEAR/5 (conservation* or set-aside*)) OR TS = (METSO NEAR/5 program*) OR TS = (Komet NEAR/5 program*) OR TS = (conservation NEAR/5 easement*) OR TS = (connectiv*) #5 = #3 AND #4 The search string was translated to Russian. A simplified version of the search string containing only key terms was used when the search interface did not have the capacity to handle the whole search string. The simplified search strings were translated to Swedish and Finnish when needed. Boolean operators were used where possible. Final search strings used are reported in Additional files 4, 5. To detect articles that were published before the data synthesis was started, a search alert was set in two bibliographic databases (Russian Science Citation Index on the Web of Science (RSCI) and Scopus) and Web of Science Core Collection (WoS)), which consists of six databases: Science Citation Index Expanded (1945-present), Social Sciences Citation Index (1956-present), Arts & Humanities Citation Index (1975-present), Conference Proceedings Citation Index-Science (1990-present), Conference Proceedings Citation Index-Social Science & Humanities (1990-present), Emerging Sources Citation Index (2015-present) (hereafter treated as one for simplicity).

Languages
This systematic review includes studies published in English, Finnish, Swedish and Russian. The language selection is based on the geographical scope of the systematic review and is limited by the language skills of the review team. Organizational websites and bibliographic databases were searched in English, except of Finnish, Swedish and Russian ones that were searched in the primary language the website was published. In addition, if the publications section / library catalogue included studies published in other of the review languages (e.g. main website language is Swedish but there are also unique publications in English), the search was conducted in those languages as well. The searches in search engines were conducted in all four languages.

Bibliographic searches
The searches in RSCI, Scopus and WoS were conducted using institutional subscriptions of the University of Helsinki and the full search strings in English and Russian were used. For the rest of the databases, no subscriptions were needed, and simplified search strings adjusted for the functionality of the search field was used. For the full search details see Additional file 5.

Search engines
Searches with internet search engines were conducted in 'private' mode to prevent the influence of previous browsing history and location on search results. The searches were conducted in English, Finnish, Russian and Swedish with simplified search strings (Additional file 4). Results were organized by relevance. After the first 50 hits, results were checked until relevant articles were no longer retrieved as advised in Livoreil et al. [38]. To safeguard against finishing the search too early, a hundred irrelevant hits were allowed. If no relevant hits appeared after a hundred irrelevant hits, the search was terminated. The maximum amount of hits shown in Google and Google Scholar is 1000, so this was also the maximum amount of screened hits per search. For the full search details see Additional file 5. The following search engines were used: Google Scholar (https ://schol ar.googl e.com/) Google (https ://www.googl e.com/)

Organizational websites
Besides known sources of potential literature, internet was searched extensively to find additional sources. Organizational websites dealing with related issues were chosen and scoping for possibly relevant material was conducted before decision on which organizational websites will be searched was made. Searches were conducted in the publications section of the website, if available. If not, the search was conducted using the "search" function of the front page. In cases where the search option was missing or it was obviously not working correctly, the websites were searched manually. Also, all Russian webpages were searched manually. The search strings/terms were adjusted to each website according to its functionality. The final list of organizational websites and the full search details in each website are given in the Additional file 5.

Supplementary searches
To supplement the search, citation chasing was undertaken in the relevant review articles included at title/ abstract stage but excluded at full text stage (the reviews, where citation chasing was conducted are identified in Additional file 6). Also, a call for unpublished data was published on the website of the Evidence-Based Forestry in Finland project (http://npmet sa.fi/en/front page/) and in ResearchGate and sent directly to stakeholder organizations that may have unpublished data on the topic.

Screening process
Articles were screened in two stages by three people: first by title and abstract (MH, AJ, SS) and then by full text (MH, AJ). At the first stage a random set of 100 articles was independently screened by all three persons. If their inclusion decisions did not vary more than 5%, the screening process could commence. If the screening decisions differed more than 5%, a second set of 100 articles were independently screened by all three persons. This procedure continued until at least 95% consensus was achieved. All discrepancies in screening decisions after one set of 100 articles were discussed to facilitate consistency in the screening process. If a reviewer was unsure whether to include or exclude an article on the title/abstract stage, the article was included.
At the full text stage, the articles were independently screened by three people (MH, AJ, SS) All the articles that a reviewer was unsure about were discussed with other research group members. Articles in Russian were screened by only one person (AJ). As they were screened after testing consistency of reviewers' decisions, the Russian speaker talked through the decision process of the articles she was unsure about. In most cases it was clear whether the article fulfilled the inclusion criteria or not.
Authors of this review who had also authored papers considered in this review did not take part on the article screening.
Full texts of articles included at the title/abstract stage were searched in the internet, Helsinki University Library, Jyväskylä University Library and by interlibrary loans. In the case an article was not found through these searches, it was classified as unretrievable. The inclusion of articles published by the authors of the review at the screening and critical appraisal stage was jointly determined by the other review group members in accordance with the eligibility and appraisal criteria.
Search records from WoS and Scopus were exported into the reference management software EndNote. The files were merged, and duplicates were removed. To conduct screening, articles were exported into Colandr, a machine-learning assisted application for conducting systematic reviews [39]. Search records from other databases were collated and the articles were screened manually.

Eligibility criteria
To be included in this review, screened articles had to pass eligibility criteria based on PECO-components, study design and language ( Table 1). The criteria were specified more in detail after publication of the research protocol [36].

Study validity assessment
All studies included in the full text stage were critically appraised and categorized. Critical appraisal was based on study design, sampling, accounting for potential effect modifiers and data analysis methods. These were considered by the authors of this review to be key variables related to generalisability and reliability of study findings. The criteria related to risk of bias and confounding factors (internal validity) and generalisability of the results (external validity). They include, for example, appropriate and representative replication, suitable outcome measuring methods, possibility of researchers to control baseline differences, measures taken to address potential baseline differences (e.g. accounting for spatial heterogeneity), and information about potential confounding factors) ( Table 2). Following the critical appraisal, studies were categorised as having 'low' , Control-intervention field studies of small habitat patches Control-intervention studies of small habitat patches where no field study was conducted but the data analysed are originally from a field study Excluded Simulation/modelling studies (even partly simulated/modelled) Habitat selection studies Studies first mapping/finding species and afterwards checking whether the place of occurrence was located in SPHP or not

Language
Included English Finnish Swedish Russian 'medium' , or 'high' risk of bias. If information was inadequate to make an assessment, the study was labelled as 'unclear' . Studies in 'low' and 'medium' categories were considered to have sufficient quality to provide reliable evidence base for quantitative synthesis and hence, the categorisation was used mainly in the narrative synthesis.
The studies were categorized based on their lowest performance, i.e. if any of the factors falls in the category 'high risk' , the study will be categorized as high risk of bias. No studies were excluded based on validity assessment, but during the meta-analysis the studies were weighted according to their category (see data synthesis and presentation section for testing robustness of results with sensitivity analysis). The validity assessment was conducted at the same time with data extraction by two reviewers (MH, AJ). Any inconsistencies or uncertainties were discussed with other research group members.

Data coding and extraction strategy
When there were multiple articles from one study site (i.e. linked articles), they were appraised as a group to avoid inclusion of duplicate data following Frampton et al. [40]. When there were results from several taxa in one article, these were treated as separate studies in data extraction. Data from included studies were extracted and saved in an Excel spreadsheet and are available as Additional file 7. To ensure consistency during data extraction, MH and AJ extracted data from five articles together before beginning the extraction independently.
Basic publication details and information about intervention, comparator, outcome, study subject (species) and study design as well as brief description of the sampling method were extracted. For quantitative analysis, outcome means, standard deviations (SDs), standard errors (SEs), and sample sizes were extracted from text, tables, and graphs. Image analysis tool WebPlotDigitizer [41] was used where necessary. In cases where information on means, SDs or SEs were not available, data on test statistics that could be used in conversions were collected. For analysis of assemblage data, total number of species (intervention + comparator) and number of shared species between intervention and comparator or other description of the community similarity were extracted. Furthermore, data on effect modifiers and potential sources of heterogeneity were extracted whenever available to enable statistical analysis of the relationships between outcomes and sources of heterogeneity. The management of the intervention surroundings was determined as intensive (clear cut, sapling stand after clear cut), moderate (middle-aged or mature even-aged production forest stand) or light (uneven-aged production forest or near natural stand).
In many articles there were data on all species and rare species, i.e. red-listed / indicator / other rare species. In those cases, data was extracted only on rare species as the purpose of small protected habitat patches is to protect especially rare, endangered species. By extracting data specifically on those species, conclusions can be drawn how well small protected habitat patches fulfil their purpose. No data on single species were extracted, but data on separate species groups (bryophytes, lichens, vascular plants etc.) were extracted when available. Two kinds of dead-wood data were extracted: volume and diversity. Whenever given, the total volume was extracted.
If an article had studies from various locations, the data on results combining these locations were extracted. If combined data were not available, data on all locations were extracted as different studies. However, the locations had to be far enough from each other based on the dispersal ability of the taxa in question to ensure results were independent. If there were studies from multiple years, the data on results combining these years were extracted. If combined data was not available, only data from the latest year were extracted (to avoid non-independence of the data).
Data were extracted by two persons (MH, AJ). Thus, to ensure consistency, a set of five studies were first coded together. Any uncertainties with data extraction were discussed among group members. To retrieve missing information or data, authors of the articles were contacted via email.

Potential effects modifiers and sources of heterogeneity
To understand better possible variation in the effects of the studies, possible effect modifiers were extracted from the studies. The list below was compiled based on the authors' experience and consultation at the stakeholder meeting and it was peer reviewed and published in the protocol of this review [36].

Data synthesis and presentation
A narrative synthesis of data from all the individual studies included was produced. All the information described in the narrative synthesis can be found in Additional files 6 and 7.
In addition to the narrative synthesis, a quantitative synthesis, i.e. meta-analysis, was conducted to assess the effects of protection of SPHPs on biodiversity. Studies that provided quantitative data on outcome measures were included in the meta-analysis. Studies that had incomplete or missing information that could not be retrieved were excluded from the meta-analysis.

Meta-analysis
We used the standardised mean difference (Hedges' d) as a measure of the effect size for species richness, abundance, and deadwood volume: where X̄i and X̄c were the means of the intervention (SPHP) and control groups (production and natural forest), s is the pooled standard deviation and J is a correction term mitigating for small sample size bias. The pooled standard deviation was calculated as where n i and n c are the sample sizes of the intervention and control groups and SD is the standard deviation. Correction term J [42] was calculated as Variance for Hedges' d was calculated as If there were non-independent cases within a study (i.e. SPHP was compared to both production and natural forest), we calculated a corrected overall sample size, N corrected , to avoid double counting the intervention sample size following the method described in Gleser and Olkin [43]: where n i is the sample size of the intervention and n c are sample sizes for the controls.
N corrected was then used to calculate corrected S pooled , J, Hedges' d and its variance for studies with non-independent cases: If information on means, SDs or SEs was missing, we used available test statistics to either calculate SDs or to convert available test statistics to Hedges' d. In the deadwood data set, we used imputation to calculate missing SDs for five studies following Lajeunesse [44]: where X j , is the observed mean of the study with missing information, and K is the number of jth studies with complete information (i). The imputation was done separately for intervention and comparator arms of the study and per comparator forest type to minimise biases.
Meta-analyses were conducted in R [45] using the rma.mv function in the metafor package [46]. The rma. mv function allows fitting of a multi-level model with a nested structure with or without moderators. In our models, studies were nested within articles as there were articles, which reported multiple outcomes. We used profile likelihood plots to check that both variance components were identifiable by the model. For deadwood, study level variance was estimated to be zero in the model, which means that the model corresponds to a model with one level only. We used inversevariance weights in the model, i.e. weights were equal to w i = (v i + τ 2 ) −1 . Cochran's Q-test [47] was used to test whether the variability in the observed effect sizes is larger than it would be expected based on sampling variability alone. A significant test indicates that the true effects are heterogeneous. When moderators were included in the model, the Q-test tests for residual heterogeneity, i.e., remaining variability in effect sizes after variability due to moderators is accounted for.
For each outcome variable (species richness, abundance, and deadwood) we first calculated the overall mean effect size. If significant heterogeneity in the effect sizes was found, we explored the influence of effect modifiers thought to be biologically significant in separate models. Following models were performed for species richness: 1. comparator, taxa and their interaction, 2. influence of forest owner, 3. intensity of management of the area surrounding SPHPs, and 4. age of the comparator forest. For the models 2-4 subsets of data were used as not all the studies reported same information and combining them in same model would have resulted in a significant loss of power. For abundance, we tested the influence of taxa and comparator only as the number of studies per comparator forest category was limited (n = 15). For deadwood, we tested only whether age of the comparator forest influenced effect sizes as the dataset was very small (n = 7).
In addition to the models above, we explored publication bias in the datasets by producing funnel plots and conducting trim and fill-tests. If outliers were found, we tested their influence on the results. We also tested whether effect sizes were influenced by study characteristics: the type of publication (MSc thesis, journal article, report), country or year when data was collected. We conducted sensitivity analyses by excluding studies that had corrected sample sizes, imputed SDs or were studies with 'high risk' of bias. The script used to run models in R is provided in Additional file 8 and the data used in these models is provided in Additional file 7.

Analysis of the assemblage data
A variety of different methods used to examine changes in species composition makes it difficult to quantitatively assess the effects of habitat modification on species composition. Hence, to have a standardized measure to assess changes in species composition, we used a simple averaging method following Nichols et al. [48]. For each study, we first extracted or calculated the number of shared species between SPHPs and comparator forests. Then we divided the number of shared species by the total number of species recorded in SPHPs. To know the overall response, we calculated the mean of all the studies and its 95% confidence intervals. The overall response was considered significant when the confidence interval did not include one.

Review descriptive statistics Search and screening of articles
Searches in Scopus and WoS were conducted from December 2nd to 3rd, 2019 resulting in 8,303 hits. Other database searches were conducted December 12th to 17th 2019 producing 4 787 hits. Search engine searches resulted 4250 hits between February 4th to 17th 2020 and searches in organisational websites (not Russian) resulted in 1914 hits during February 7th to 12th 2020. The search alerts were on from December 2nd, 2019 to March 19th, 2020 producing 201 hits. The manual searches in Russian databases produced 3 articles. Altogether the searches resulted in 19,458 hits. Duplicates were first removed from the records imported to EndNote (hits from Scopus and WoS). Because articles from other sources could not be imported to EndNote, they were screened manually, and duplicates were removed during title/abstract screening. After duplicate removal and screening by title and abstract, 336 articles remained (Fig. 2). No articles were obtained through call for unpublished data. One set of 100 randomly chosen articles was needed to achieve 95% consistency in the screeners' (MH, AJ, SS) decisions at the title and abstract stage.
Thirteen full texts were not found which means that a total of 323 full texts were read. After full text screening 38 articles were included. In addition, two articles were added during citation chasing, and one article was found by chance encounter outside the predetermined sources increasing the total number of articles included by full text to 41 (Fig. 1). The most common reason for exclusion was that there was no eligible intervention (no small-scale protected habitat patches or these patches were combined with other forest types). Also, the lack of suitable comparator was a common reason for exclusion (for example articles, where all the studied areas were SPHPs). The list of all articles included and excluded by full text with the reasons for exclusion is given in Additional file 6 and a description of the search of each database and the articles included in each screening stage is given in Additional file 5.
Of the 41 articles included at the full text stage, 20 belonged to some group (articles from the same study site formed a group). There were eight separate groups with 2-4 articles in each. Despite the substantial number of articles belonging to a group, a total of 35 articles had independent data. Therefore, the number of articles and studies included in narrative synthesis was 35 and 174, respectively. Some studies included in narrative synthesis did not have suitable data for further data synthesis and eventually 28 articles with 127 studies were included for quantitative analysis (Fig. 1).
Of the articles included at full text stage, three were authored by an author of this review (MM). Two of these had independent data and were also included in the quantitative analysis. The screening and critical appraisal for these articles were conducted by MH, SS and AJ (Fig. 2).

Sources of the articles included in narrative synthesis
The majority of the included articles (22 of the 35 articles, 63%) were found through searches in Scopus and WoS. Three articles (9%) were found in other bibliographic databases, six (17%) in Google Scholar and two (6%) in organisational websites. Citation chasing and other sources both resulted in one article.
There were three types of articles: journal articles, master's theses, and reports. 25 articles were journal articles, six were master's theses and four were reports (25 peer reviewed and 10 grey literature articles). In addition to WoS and Scopus, journal articles were found in Google Scholar (Finnish search), Doria (Finnish search) and from other sources, one from each source.
Most of the articles were written in English (27 out

Narrative synthesis including validity assessment Locations
Most of the studies were conducted in Sweden, followed by Finland, Norway, and Latvia. No studies from other eligible countries were included (Fig. 3).

Interventions
Almost all the studies considered were WKHs. There were only five studies from two articles where the areas were not WKHs but protected by METSO program [50] or by certification [51]. In most cases, it was not written what kind of WKHs were studied but stream-side forests [52][53][54][55][56] and herb-rich forests [57] were mentioned.

Comparators
In 113 studies (65%) the comparator was production forest, whereas in 61 (35%) it was natural forest. Production forest comparators were mostly described as 'ordinary managed forest' which usually means even-aged forest where thinning has been done and which consists of only few tree species. Age of the production forest was not mentioned in majority of the articles but where it was, it varied between 50 [53] and 120 years [50]. In one article [58] production forest comparator consisted of a mix of clear cuts and older forests. The natural forest comparator mostly referred to protected natural reserves and national parks. Six articles (with 42 studies) had both production forest and natural forest comparators [51,[59][60][61][62][63].

Outcomes
Species diversity (i.e. richness) was the outcome in 61 studies (35.1% of the studies). The number of studies for other outcomes were: assemblage 44 (25.3%), individual abundance 39 (22.4%), dead wood volume 27 (15.5%) and dead wood diversity 3 (1.7%). Studied subjects could be divided into eight groups (Fig. 4). The most studied  [49] showing literature sources and inclusion/exclusion process. Note that duplicate removal was conducted at two stages. Duplicates were first removed from the Scopus and WoS databases after the records were exported to EndNote. Because records from other searches could not be exported to EndNote, duplicate removal was conducted at title/abstract stage subject was fungi (mostly polypores) and vascular plants followed by dead wood and bryophytes. 100 studies (57%) focused on rare (red-listed) species, indicator/signal species or species associated with dead wood.

Potential effects modifiers and sources of heterogeneity
Information about potential effects modifiers and sources of heterogeneity was mostly insufficient (Additional file 7). Geographic location, study year and forest type were reported in most studies. Climatic conditions, soil type, altitude, category and area of SPHP, age of the forest, certification status and owner of the forest were reported occasionally.

Study validity assessment
Almost all studies (168 of 174) were categorised as medium risk of bias during the critical appraisal and the rest were assessed as having a high risk of bias because  Fig. 4 Number of studies by study subject no effect modifiers were identified [64,65] or because the intervention and comparator data were collected with different methods [66]. The studies with high risk of bias were not excluded from the quantitative analysis but their influence on results was tested by conducting a sensitivity analysis. The reason for absence of studies with low risk of bias was the lack of experimental studies. The critical appraisals for all studies are reported in Additional file 7.

Data synthesis Description of studies included in the quantitative synthesis
127 studies (from 28 articles) had suitable data for quantitative analysis. These evidences represent 73% of all the studies included in the narrative synthesis. The outcome was species richness in 45, abundance in 30, assemblage in 26 and dead wood volume in 26 studies. There were not enough data to conduct meta-analysis on dead wood diversity. All the studies where intervention was not WKH (but METSO habitat or protected by certification) had dead wood volume as outcome.
75 studies compared SPHP to production forest and 52 studies were comparing SPHP to natural forest. 63 of the studies were conducted in Sweden, 49 in Finland, 15 in Norway and none in Latvia. Vascular plants and fungi were the most studied taxa (31 and 25 studies, respectively). 72 studies were published in journal articles, 49 in master's theses and 6 in reports. Three studies were assessed as having a high and the rest as having a medium risk of bias. In all three high risk studies the outcome was dead wood volume.

Species richness
The overall mean effect size was significantly different from zero (d = 0.812, p < 0.0001, 95% CI 0.407, 1.217) meaning that SPHPs had significantly more species than comparator forests. There was significant heterogeneity in the effect sizes (Q = 298.592, p < 0.0001, n = 45). Forest plot showed that there was a clear outlier (article id 7 with production forest comparator) (Fig. 5)

Fig. 5
Forest plot for the overall species richness between SPHPs and comparator forest types (production and natural forest). Effect sizes on the right side of zero mean that SPHP has more species than comparator forest. The grand mean noted by a diamond shows the summary effect of all the individual effect sizes shown on the right side of the figure. The error bars represent 95% confidence intervals by excluding studies that had multiple treatments and thus, corrected sample sizes. Exclusion had no effect on the results: species richness remained significantly higher in SPHPs than in comparator forests. (d = 0.816, p = 0.0003, 95% CI 0.371, 1.26, n = 33). Similarly, heterogeneity remained significant (Q = 190.691, p < 0.0001, n = 33). Type of publication (grey or peer-reviewed literature), country or year when data was collected had no influence on effect sizes (QM = 1.762, p = 0.779, n = 45). Funnel plot did not show any obvious asymmetry (Additional file 9). Further analysis showed that species richness in SPHPs does not differ significantly from natural forest but they are significantly more diverse than production forests (d = 1.0807, p = 0.04 0, 95% CI 0.252, 1.891, n = 45). Taxa had no effect on species richness (p > 0.05 for all). When interaction between forest type and taxa was studied, we found that indicator species (lichens and fungi) were significantly more common in SPHPs than in production forests but this was based only on one study (d = 10.813, p = 0.0001, 95% CI 7.085, 14.541). When we removed the outlier study (id 7) and reran the analysis, no impact of comparator, taxa or their interaction was found (QM = 18.267, p = 0.076, n = 44). Subsets of data were used to test the influence of age of the comparator forest, forest owner (private or company; there was only one state owned forest, so it was excluded from the analysis), and intensity of management in the area surrounding the SPHP. There was no impact of the management intensity of the area surrounding the SPHPs on the effect sizes (QM = 0.494, p = 0.482, n = 26). Also, the owner of the forest (QM = 2.773, p = 0.096, n = 19) nor the age of the comparator forest (QM = 2.634, p = 0.104, n = 16) influenced species richness.

Individual abundance
The overall mean effect size for abundance showed that SPHPs have significantly more individuals than comparator areas (d = 1.911, p = 0.0003, 95% CI 0.868, 2.954, n = 30) (Fig. 6) but there was significant heterogeneity in the results (Q = 254.813, p < 0.0001). The effect remained significant after removal of studies that had multiple treatments and thus, corrected sample sizes (d = 2.101, p = 0.0007, 95% CI 0.891, 3.311, n = 22). The influence of year when data was collected, country and publication type (journal article or master's thesis) on effect sizes was   Fig. 6 Forest plot for individual abundance between SPHPs and comparator forest types (production and natural forest). Effect sizes on the right side of zero mean that SPHP has more individuals than comparator forest. The grand mean noted by a diamond shows the summary effect of all the individual effect sizes shown on the right side of the figure. The error bars represent 95% confidence intervals tested but none of these were significant (QM = 1.486, p = 0.829, n = 30). Visual inspection of the funnel plot did not show an asymmetrical pattern. In addition, based on the results of the trim and fill-method no publication bias was detected (Additional file 9). When the influence of comparator type (natural or production forest) and taxa were tested, no differences in species abundance were found (QM = 7.216, p = 0.205, n = 30).

Deadwood
There was more deadwood overall in SPHPs than in controls (d = 0.625, p = 0.025, 95% CI: 0.079, 1.171, n = 26) (Fig. 7). There was significant heterogeneity in the effect sizes (Q = 882.700, p < 0.0001). When studies with imputed SDs were excluded, there was more deadwood in SPHPs than in controls, but the overall effect size was only marginally significant ( Effect sizes were not influenced by country, publication type and data collection year when all the studies were included (QM = 4.054, p = 0.399, n = 26). When studies with imputed SDs were removed, study year and publication type influenced the results (QM = 9.532, p = 0.049, n = 24). This was caused mainly by one study that was published in a peer-reviewed journal and had a relatively large positive effect size (d = 4.822).
No publication bias was visually detected, and this was confirmed by the results of the trim and fill-test (Additional file 9).
Comparator had significant influence on effect sizes (QM = 9.901, p = 0.002, n = 26)) but heterogeneity remained significant (Q = 232.258, p < 0.0001) indicating that there are other moderators influencing effect sizes. Based on the analysis, deadwood volumes were significantly higher in SPHPs compared to production forests (d = 1.564, p = 0.002, 95% CI 0.590, 2.539, n = 26). The influence of age of production forest was tested but no effect was found (QM = 0.954, p = 0.323, n = 9). When studies with imputed SDs were excluded, the results   Fig. 7 Forest plot for deadwood volumes between SPHPs and comparator forest types (production and natural forest). Effect sizes on the right side of zero mean that volume of deadwood is greater in the SPHP than in comparator forest. The grand mean noted by a diamond shows the summary effect of all the individual effect sizes shown on the right side of the figure. The error bars represent 95% confidence intervals remained similar. Deadwood volumes were significantly higher in SPHPs compared to production forests (d = 1.491, p = 0.0005, 95% CI 0.647, 2.336, n = 24). The excluded studies did not record the age of production forest.

Studies with high risk of bias
When studies with high risk of bias (n = 3) were excluded from the dataset, the overall effect size was only marginally significant (d = 0.612, p = 0.06, n = 23). Also, age of the comparator forest had marginally significant impact on the effect sizes (d = − 0.009, p = 0.061, 95% CI − 0.019, 0.0004, n = 7). Otherwise, the results were similar to the analyses where studies with high risk were not removed. Results without high-risk studies are included as Additional file 10.

Assemblage
Similarity of species composition could be statistically assessed in 26 studies. Further 18 studies provided some information about species composition (Additional file 11). A mean of shared species between SPHPs and controls was calculated to compare similarity of species composition. There were 14 studies on production forests and 12 on natural forests that provided suitable data for comparison. On average, production forests and SPHPs shared approximately 58% of species and natural forests and SPHPs 48% of species after the values were standardized (Fig. 8). This represents significant differences in community composition regarding both forest types.

Evidence of effects
Overall, the evidence supports using SPHPs, specifically WKHs, to maintain and increase biodiversity in production forests. There were significantly more species in SPHPs than in production forests concurring with earlier systematic review on WKHs [35]. Also, dead wood volumes were higher in SPHPs than in production forests. The larger volume of dead wood in SPHPs is not surprising since removal of dead wood in SPHPs is usually forbidden, while in production forests dying and dead wood have mainly been harvested in the past. However, nowadays it is often recommended to leave some deadwood also in production forests (for example [67,68]). Increasing the amounts of deadwood have been shown to have a positive effect on saproxylic insects and fungi (i.e. species dependent on dead wood) [69], which indicates positive reinforcement of biodiversity benefits of SPHPs.
Contrary to expectations, management intensity of the area surrounding SPHPs did not have an influence on species richness. There have been studies indicating that even in the case of protected areas, the management of the surroundings can affect species diversity (for example [70] on birds). However, these effects can vary between species groups, and most of the studies in this review were about plants and arthropods. In addition, it must be noted that in many studies the management intensity of the area surrounding SPHPs was not reported, and none of the studies reported light management intensity. Hence, the lack of effect can also mean that taxa react similarly to both moderate and high intensity management. It is also possible that SPHPs are suffering from extinction debt, i.e. even though the species are still present in SPHPs, they can become extinct in the future because isolation of the patch can lead to lack of gene flow and viable populations in long term [33,34]. However, long-term studies are needed to explore the phenomenon. Owner of SPHP did not have influence on species richness.
No differences in species richness between SPHPs and natural forest were noticed concurring with the results of Timonen et al. [9]. When species assemblages were studied, differences between SPHPs and production forests as well as SPHPs and natural forests were found. This means that although SPHPs harbour similar number of species than natural forests, they are not 'mini natural forests' in terms of species composition. This may result partly from the definition criteria of WKHs (all SPHPs analysed here were WKHs). In some cases, the criteria include indicator species, which means these species inevitably exist in the WKH. In addition, WKHs are already by definition so called biodiversity hotspots where large number of species occur within small area. In production forests and natural forests, species are distributed within a larger area. The same species may not occur in the randomly chosen study plot in the comparator forest, but this does not mean that the species could not occur in some other  Fig. 8 The mean proportion of shared species between SPHPs and control areas with 95% confidence intervals. Data were standardized by the total number of species recorded in SPHPs (the number of species = 1) location within the same forest. Also, the differences in species assemblages probably reflect differences in habitat requirements, e.g. in terms of light, humidity and exposure to edge effects. However, more studies are needed, especially regarding species assemblages in SPHPs and neighbouring production forests, to have conclusive evidence on the causal links that lead to certain composition of species. Species abundance in SPHPs was higher overall than in natural forests and production forests but not when the effect was separated per comparator forest type. The explanation can be related to special attributes of SPHPs. As the evidence shows, they can harbour more resources, like dead wood, and therefore they can host more individuals dependent on those resources.

Magnitude of effects
Where SPHPs had significant or even marginally significant impact on species richness or abundance, the effect sizes were in most cases large. Usually, effect size of 0.2 is considered a small effect, d = 0.5 an intermediate effect and d = 0.8 a large effect [71]. In our results, effect sizes were commonly above one, and large effect sizes indicate a strong response from the studied groups. Considering that most of the observational studies had aimed to minimise bias from environmental variation across study sites, we are confident that the large effects found in this systematic review are representative of true effects in nature for the studied species groups.

Review limitations
Even though different sources were searched widely, it is likely that some relevant articles were not found. Because the languages and geographical scope of the review did not fully overlap, there may be relevant articles written in Norwegian, Estonian, Latvian, Lithuanian, or French (Canada) that were not included in this review. Searching every organisational website or database that may include relevant material is impossible, but the search in this review was comprehensive and no publication bias was detected. It is worth noticing that almost 40% of the articles included in this review were found outside the citation databases (Scopus and WoS) and almost 30% were grey literature.
Studies included in this review were not 100% comparable with each other but restricting the eligibility criteria even more could have resulted in a situation where quantitative analysis was not possible. For example, the comparators within one comparator class (natural or production forest) differed. The age of production forests varied, and natural forests were both protected and unprotected. However, these differences reflect variability in nature, and hence natural sources of heterogeneity.
Comparability of production forests was enhanced by excluding articles where the production forest was already, before the study, known to have low or high conservation values. These were, for example, areas not accepted in the METSO program and areas receiving environmental subsidies.
Geographical distribution of the articles included in this review was limited. There were many articles from Finland and Sweden, few from Norway and only one from Latvia. Since the concept of WKH is widely used in all the Baltic Countries and to some extent also in Western Russia, the lack of studies from these countries is unexpected. In Latvia there has been a lot of research on WKHs, but the focus has been on describing the species in WKHs or comparing different types of WKHs (for example [63,64]). Also, many forests in Latvia are temperate and studies conducted in non-boreal forests were excluded from this review.
The results can also be affected by the studied taxa. Most of the studies included in this review were concentrating on small, sedentary, and poorly distributing organisms like insects, plants, and fungi. The importance of SPHPs for these small and local organisms is clear but based on this review the importance of SPHPs for example for mammals and birds is not known.
Not all studies included in narrative synthesis had suitable data for quantitative analysis and hence the quantitative synthesis could not be based on the whole knowledge base. Reporting on the effect modifiers was also inadequate in many articles. This naturally affected metaanalyses when the influence of effect modifiers could not be analysed. Especially age of the comparator forest and the type of management of the area surrounding SPHPs would have been interesting effect moderators to analyse with more data. The lack of studies comparing WKHs with differently managed surroundings was highlighted already in a former systematic review by Timonen et al. [35].

Review conclusions Implication for policy/management
There were more species, individuals, and dead wood in WKHs than in production forests which indicates that setting aside SPHPs helps to protect biodiversity in production forests. However, species assemblages differed between SPHPs and natural forests meaning that SPHPs cannot be considered as a substitute for large natural forests like national parks and reserves. But SPHPs do complement the protected area network and may increase connectivity between larger reserves by providing habitats for species with restricted dispersal abilities.
Most of the studies concentrated on rare/indicator/ signal species and data extraction were also targeted towards these species. Therefore, the results cannot be straightforwardly generalized to all species, not even to all species of the most studied species groups. However, because the purpose of setting aside small habitat patches is to protect especially rare species, the results indicate the success of achieving this target.
Also, clear majority of the studies were about WKHs and therefore, the results cannot be generalized directly to other types of small protected habitat patches when they do not share similar attributes with WKHs. The geographical distribution of the included studies should also be considered when generalizing the results. Most of the studies were conducted in Finland and Sweden, and even between these two countries the definitions, sizes, and protection status of WKHs differs.

Implication for research
Based on this systematic review, there are knowledge gaps regarding geographical distribution, interventions and studied organisms. More research on SPHPs in general is needed, and especially research comparing SPHPs and production forests or natural forests is needed in Baltic Countries and Russia, and to some extent also in Norway. Alongside WKHs, future research should concentrate on other small protected habitat patches, such as voluntarily protected areas.
The results of this review indicate that importance of SPHPs to mammals and birds is not known and more studies concentrating on these species is needed. In addition, replicating past studies in the future would offer more information about the ability of SPHPs to maintain biodiversity in the long run. To enhance reliability of the individual studies and meta-analyses in the future, it is important to report more in detail about potential effect modifiers and reasons for heterogeneity. Especially reporting the management intensity of the area surrounding the SPHPs is important for achieving information on the ability of SPHPs to maintain biodiversity within different landscape matrices.