Geography of disservices in urban forests: public participation mapping for closing the loop

ABSTRACT Studies on cultural ecosystem services (CES) in urban forests using Public Participation GIS mapping are a well-established field of research. However, far fewer studies do exist that also illuminate perceived disservices, thus omitting an important part of human-forest relationships. Advancing knowledge on unpleasant places can promote outdoor recreation and help forest managers focus on disservices that matter. This study aims at spatially identifying people’s perceptions in urban forests with focus on unpleasant experiences. We elicited perceived disservices and CES spatially by collecting Public Participation GIS data (PPGIS) from city dwellers in three urban forests in Germany’s Southwest (755 respondents mapped 1552 places of disservices). We investigated relationships between sociodemographic characteristics and the number of mapped places using Spearman’s rank. We also explored the spatial concurrence between disservices and CES using Spearman’s rank. We identified a hierarchy of importance of forest visitors’ perceived dislikes. We revealed that negative perceptions of urban forest visitors originated broadly from people and their traces and rarely from the ecosystem according to existing ecosystem disservices typologies (EDS). We found a significant relationship between age and the number of mapped services. We found clusters of disliked places and correlations between disservices and CES, indicating that hotspots of CES in turn are also hotspots of disservices. We conclude that city foresters may best address disservices by advancing guidance-concepts for spatial movements of urban forest visitors within the forest to avoid conflicting clusters of various stakeholders. Our results may help city foresters to better manage both the forest and visitors’ various interests.


Introduction
Since 2003, when the Millennium Ecosystem Assessment typology was agreed upon by many researchers around the globe, numerous studies have provided insights into urban nature's invaluable benefits for human health and well-being (Bieling et al. 2014;Pietilä et al. 2015;Völker et al. 2018). Urban nature offers besides provisioning-, and regulating-, also cultural ecosystem services (CES). CES represent non-material nature values for city dwellers, such as spiritual and religious values, educational values, social relations, cultural heritage, and recreation (Millennium Ecosystem Assessment 2005). Urban forests represent a substantial component of nature in cities and encompass urban parks, single trees, urban woodlands, and woods with closed canopy covers (Konijnendijk et al. 2005). People experience CES in urban forests while recreating, enjoying the beauty of nature, or being socially engaged (Pulighe et al. 2016). The global COVID-19 pandemic in 2020 has highlighted the growing importance of urban forests as critical infrastructure for CES delivery since visitor numbers have increased significantly (Derks et al. 2020;Beckmann-Wübbelt et al. 2021).
People perceive urban forests as a symbolic environment, created by a human act of conferring meaning on nature and on the environment (Greider and Garkovich 1994). Focussing on people, the perception of urban forests is consequently socially constructed (Stedman 2003;Winthrop 2014). In fact, benefits from ecosystems are not produced by ecosystems independently of humans, but arise because of people's interactions with the ecosystem: ecosystem services are thus a co-production of human-nature interactions (Fischer and Eastwood 2016;Palomo et al. 2016;Raymond et al. 2017;Torralba et al. 2018). Quantifying and evaluating CES is challenging because of the complexity of interactions (Fish et al. 2016b) and the 'intangible nature' of CES Hernández-Morcillo et al. 2013). Nevertheless, the study of CES of urban green infrastructure is now a well-established field of research. Numerous publications have explored the linkages between ecosystem characteristics (supply) and people's characteristics and needs (demand) for delivery of CES (Hegetschweiler et al. 2017a). Participatory mapping approaches (using Public Participation Geographic Information Systems, PPGIS) are a sound method being used to assess preferences and values directly to quantify demand (Palomo et al. 2013;. PPGIS can capture the heterogeneity of human desires, values, and preferences, which is essential for understanding human demand (Wolff et al. 2015). In fact, existing studies identified the linkages between CES hotspots (demand) and physical landscape features (supply) and revealed a hierarchy of importance of various types of CES, also with special focus on urban forests (Agimass et al. 2017;Baumeister et al. 2020;Gerstenberg et al. 2020).
However, these studies did rarely encompass negative perceptions of urban forest visitors and far less attention has been given to equally address negative effects of natural environments. Such knowledge would be needed to provide a complete picture on people's preferences: it would illuminate two sides of the same coin and could manage to close the loop. In fact, nature may account for costs in terms of maintenance and management (Mullaney et al. 2015;McPherson et al. 2016;Saunders and Lemasson 2020). Moreover, nature may have negative effects on people's health (Douglas 2012): Wildlife and animals may be dangerous (Bixler et al. 1994;Butler et al. 2003;Destefano and Deblinger 2005) or may cause danger to humans through transmission of zoonotic pathogens (Dunn et al. 2010;Ostfeld and Keesing 2017;Ehrmann et al. 2018). Furthermore, urban trees may cause allergic reactions to humans (D'Amato 2000; Baro et al. 2014;Livesley et al. 2016), emit unpleasant odours (Shackleton and Mograbi 2019), or create nuisances by dropping leaves in autumn (Flannigan 2012;Camacho-Cervantes et al. 2014). Likewise, urban trees may be seen as potential hazards in storm events (Conway and Yip 2016) or disliked because of blocking views (Schroeder and Anderson 1984;Jorgensen et al. 2002;Herzog and Kropscott 2016). Finally, urban forests may be perceived unpleasantly because of missing aesthetics or poor maintenance Rodríguez-Morales et al. 2020).
The negative effects of nature have been referred to as ecosystem disservices (EDS) in the literature. The term EDS describes functions or properties of ecosystems that cause effects that are perceived as harmful, unpleasant, or unwanted (Lyytimäki 2015). However, there is ambiguity in the term 'EDS' since some negative effects of nature, e.g. wood damage by wood-boring beetles, could, in turn, support greater functional diversity of beneficial species such as cavity-nesting bees (Sydenham et al. 2016). Moreover, given many studies that focused on forest preferences (Gundersen et al. 2015;Ciesielski and Sterenczak 2017) it is apparent that some perceived services such as shading of trees may be in turn perceived as disservices because of blocked views (Lyytimäki et al. 2008). A growing number of researchers argue that ecosystem services research should incorporate both services and disservices equally (Shapiro and Báldi 2014;Schaubroeck 2017) to provide a full picture of how people benefit or suffer from nature (Lyytimäki and Sipilä 2009;Sandbrook and Burgess 2015;Shackleton et al. 2016;Blanco et al. 2019) Typologies of EDS have been developed (Lyytimäki 2014;Shackleton et al. 2016;Vaz et al. 2017), including some with focus on urban nature (Gómez-Baggethun and Barton 2013;von Döhren and Haase 2015). Urban trees and forest-specific typologies encompass disservices causing social problems (e.g. fears of trees or crime), health problems (e.g. allergy and animal attacks), visual and aesthetic problems (e. g. darkness and blocked view), environmental problems (e.g. cracked sidewalks, volatile organic compounds), and costs (e.g. maintenance, administration) (Escobedo et al. 2011;Roy et al. 2012;Kirkpatrick et al. 2013).
Existing studies focusing on EDS do not reflect the full picture of negative aspects that people may experience when visiting urban forests. There are additional disservices that would not be encompassed by existing typologies for EDS since they are of a more anthropogenic nature. These additional disservices encompass for instance fear of crime (Weinborn et al. 2017;Mak and Jim 2018;Williams et al. 2020) and anthropogenic noise (Liu et al. 2018a(Liu et al. , 2018b(Liu et al. , 2017Song et al. 2018). Furthermore, forest visitors may perceive discomfort because of garbage dropped by other people (Agbenyega et al. 2009;Fish et al. 2016a), bad forest road conditions (Rodríguez-Morales et al. 2020), or forest management that impairs aesthetic values (Bethmann et al. 2018). Finally, there may be other disservices perceived by city dwellers that are still unknown. In this study, we understand the term disservices encompassing both EDS and disservices of anthropogenic origin.
To better understand city dwellers' individual perceptions of CES and disservices in urban forests, it is vital to look at their knowledge about the forest as well as sociodemographic characteristics (Arnberger 2006;Plieninger et al. 2013;Spangenberg et al. 2014;Hegetschweiler et al. 2017a;Meyer et al. 2019). Previous studies indicate that economic activity status, age, number of children, ethnicity, educational attainment, and marital status significantly influence urban forest visit-related preferences (Mak and Jim 2019). For example, Hegetschweiler et al. (2017b) found that frequency of visits positively affects visual attractiveness of forests. Meyer et al. (2019) found that age, education, household income, family status, and recreational activity affect forest perception. Similarly, Tyrväinen et al. (2003) found effects of gender, age, housing style, and education on the evaluation of forests' scenic beauty. However, several studies found no effects of sociodemographic characteristics on visual attractiveness of forests (Eriksson et al. 2012;Chen et al. 2016;Hegetschweiler et al. 2020). In summary, there is a knowledge gap on linkages between urban forest visitors' socio-demographic characteristics and visit-related preferences. This information may aid in visitor-group specific urban forest management and decision-making.
Finally, we want to explore whether there are spatial bundles of disservices and perceived CES. Existing studies found that various types of CES spatially overlap, indicating that people do not clearly separate one cultural service category from the other. This can be understood as evidence of the interlinked nature of CES (Daniel et al. 2012;Bieling and Plieninger 2013) and the bundled occurrence of CES (Raudsepp-Hearne et al. 2010). Furthermore, existing studies also indicate a spatial congruence between CES and disservices , but this needs further testing.
In summary, little research has been done that spatially examined a broad range of negative perceptions of urban forest visitors in addition to CES. In Germany, forest authorities have sound inventory data to assess ecosystem services supply. Additionally, forest function maps help balancing economic, ecological, and social objectives in forest management. Forest function mapping as a planning approach is different from the ecosystem services framework (Meyer and Schulz 2017;Tiemann and Ring 2018) and based upon a regionalscale expert assessment of forest agencies that incorporate the expertise of other sectorial authorities (for example, agencies for water resource management). The public is involved either in participatory workshops or through public displays of forestry plans. Forest function maps, however, do not illustrate perceived urban forest values, and hotspots of visitor frequencies are not finely graduated -they only display heavy and moderate recreational use intensity. Concerning disservices, there are no maps available to forest services that account for disliked places -neither spatially nor qualitatively to indicate where and to what extent urban forest visitors are concerned. However, forest managers need sound data on the spatial distribution of unpopular places in urban forests to ameliorate disservices in urban forests on-site and focussing their capacities and monetary endeavors on disservices that matter. This study therefore aims at exploring perceived disservices and CES (demand) in three urban forests in Germany's Southwest. Consequently, we raise the following research questions and hypothesize as follows: (1) How important are disservices for urban forest visitors? (2) Do sociodemographic characteristics of respondents have an effect on perceived services and disservices? (3) Are mapped disservices spatially clustered? (4) Are there spatial bundles of both CES and disservices and, in turn, are there cold spots with little public interests within urban forests? We hypothesize that mapped disservices are clustered and spatially congruent with CES.

Study forest
We selected three urban forests in the federal state of Baden-Württemberg in Germany's Southwest to gather city dwellers' perceived disservices and CES.
The urban forests we selected as case studies represent varying degrees of peri-urbanity and urbanity (Konijnendijk et al. 2005). The forests -located in the districts of Stuttgart, Karlsruhe, and Heidelbergprovide important green infrastructure for outdoor recreation. Selected urban forests are characterized by broadly closed canopy covers and thus form contiguous pieces of forest. The urban forest of Karlsruhe touches the inner-city-center. The urban forests of Stuttgart and Heidelberg have a more peri-urban character, being adjacent to urban areas. The forest in the district of Heidelberg (Schwetzinger Hardt) is surrounded by eight middle-sized municipalities (275,000 inhabitants). The city of Heidelberg is the largest city in the district, and downtown is roughly 7 km away northeast. The urban forest in the district of Karlsruhe (Karlsruher Hardt) is adjacent to the city center and ranges about 15 km north from downtown. It is also surrounded by five middle-size municipalities located north of the city of Karlsruhe (375,000 inhabitants). The urban forest in the district of Stuttgart (Rotwildpark) is located west of the city and touches two larger cities (810,000 inhabitants) ( Figure 1). The three urban forests were selected because of their similarity in terms of forest coverage (Stuttgart 81%, Karlsruhe 84%, and Heidelberg 80%) and size (Stuttgart 51 km 2 , Karlsruhe 53 km 2 , and Heidelberg 45 km 2 ). Furthermore, all three forests are owned by the federal state and subject to similar management approaches. The three urban forests are managed according to multipurpose forestry objectives: social, ecological, and economic aims are equally addressed. In the urban forests, there are areas dedicated for recreation and conservation as well as areas for water, climate, noise, and air pollution protection. Logging activities take place in all three study areas.

Survey
We used a map-based online questionnaire tool to assess the perception of CES and disservices in urban forests, following standard PPGIS protocols (Brown and Fagerholm 2015). We expected to reach more participants with broader representativeness using web-based approaches instead of paper-based mapping on locations. Existing studies indicate that an internet-based PPGIS system is relatively easy to use ) and spatial results between paper-and web-based mapping approaches typically do not differ (Pocewicz et al. 2012). The used survey-tool (Mapita, www.maptionnaire.com, 2017) allows for combining survey questions with an online mapping interface to collect spatially explicit data. The survey consisted of 16 web-pages and started with descriptions of objectives, technical user instructions including an explanatory video, and the privacy policy. Subsequently, we asked participants about their access to and familiarity with the urban forest. In the following section, we asked for sociodemographic background characteristics, including age, gender, education degree (including school, university, and vocational education), marital status, housing, and work situation, as well as income. Furthermore, using continuous scales, we asked participants to rate their life satisfaction, mental and physical health, place attachment, expertise in urban forestry, and importance of urban forest design. Also, we asked how important provisioning, regulating, and supporting ecosystem services were perceived using examples from the urban forests. Finally, we asked for intensity of recreational use of the urban forests and how participants were advertised for participating in the survey. In the main part of the survey, we asked participants to map CES and disservices in the respective forests. We labelled point markers to be placed within the forest. In terms of CES, we predefined categories based on the Millennium Ecosystem Assessment (2005) and according to existing studies (Rall et al. 2017). The results of cultural ecosystem services perceptions have been published in Baumeister et al. (2020). Concerning disservices, we referred to suggested categories in existing literature (Escobedo et al. 2011;Roy et al. 2012;Kirkpatrick et al. 2013;von Döhren and Haase 2015;Vaz et al. 2017). Each pre-defined marker included a short description of the respective category to operationalize the intangible CES-and disservices framework for participants (c.f. Table 1).
The supplementary category labelled 'place of discomfort: places I find uncomfortable for another reason (e.g. forest management, anthropogenic noise, garbage, missing accessibility, dogs or other bothering visitors)' was added to elicit additional disservices. Using this marker was followed by the open question: 'What makes this place uncomfortable for you?'. Two types of maps were provided to enhance orientation for participants in the urban forest: an orthophoto-map (Microsoft® BingTM Maps 2018) and a street map (OpenStreetMap contributors 2018). Respondents were able to zoom in and out and the default map displayed the center of the respective urban forest. We neither limited the number of markers per respondent nor the places where markers could be placed spatially within the urban forest to account for individual differences in perception of urban forests. Likewise, we did not ask participants to use every individual marker being provided. We pretested the survey with students and elderly people not familiar with the topic. Furthermore, we asked colleagues and friends to take part in the survey in regard to another urban forest not presented in this paper. Subsequently, we adjusted the text of the survey by making some technical terms more tangible for people not acquainted with the topic. Also, we decided to produce a short technical video in which we described the mapping procedure to make mapping easier to accomplishespecially for the elderly.

Sampling approach
We advertised study participation using various approaches. We drew a random sample of households, addressing city dwellers between 16 and 75 years of age. The gender ratio for the random sample was 50/ 50. We sent flyers to advertise the study to 9327 homes -located in cities that were adjacent to the urban forests (district of Stuttgart 3274, district of Karlsruhe 3194, district of Heidelberg 2859). Addresses were provided by city administrations. Additionally, we sent 1189 letters to affiliated people with the forestry administration of the state of Baden-Württemberg (district of Stuttgart 414, district of Karlsruhe 423, district of Heidelberg 352), such as timber and venison customers, hunters, and sawmill industry representatives located in the targeted cities. In addition, we advertised our study by sending emails to schools and local associations. We were also present in radio and newspapers. Advertisement of the study was in German. The survey was online from May to September 2017.

Identifying the relative importance of disservices
We first used inductive category formation (Mayring 2014) to analyze the responses to the open question 'What makes this place uncomfortable for you?', which followed after using the 'place of discomfort' marker. This involves first defining categories according to given answers and second assigning each answer according to this scheme. Statements that did not match any category and were mentioned only once were subsumed to the category 'other'. Statements that involved several meanings were assigned to each addressed category (e.g. 'place of bad smell and loud place' were assigned to 'odour pollution' as well as 'anthropogenic noise'). This resulted in a differentiation of disservices.

Investigating the impact of individual background characteristics on perception of services and disservices
We used Kruskal-Wallis tests and Wilcoxon ranksum tests to investigate potential effects of gender, education degree (including school and university education), marital status, housing and work situation, and income, as well as expertise in urban forestry on the total of mapped liked and disliked places (across all services and disservices categories) per person. Furthermore, we used Spearman's rank correlation to examine the relationship of age, perceived life satisfaction, mental and physical health, place attachment, familiarity with the respective forest, and rated importance of urban forest design to the number of mapped popular and unpopular places. We used R statistical software (R Core Team 2020, version 4.0.2) for these data analyses.

Visualizing hotspots of disservices
We visualized the spatial distribution of perceived disservices within the urban forest using hotspot maps. We calculated the average distance from every mapped point to its nearest neighbour based on feature centroids to further analyse use intensity patterns. Therefore, we used the spatial statistics tool 'average nearest neighbour' in ArcGIS version 10.7.1 (ESRI 2019). Subsequently, we visualized spatial density of mapped unpopular places for each urban forest using the Kernel Density tool for point features (ESRI 2019). The tool calculates the density of point features around each output raster cell and fits a smoothly curved surface over each point (Mitchell 1999). To create this density surface, a search radius around each cell center needs to be defined as it specifies the level of smoothing. Small values produce peaky estimates, whereas large values produce smooth functions (Bivand et al. 2013). The Kernel Density tool then calculates the density of each raster cell output by adding the number of features of all kernel surfaces where they overlay the raster cell center (Mitchell 1999). We used the default search radius -which is calculated based on the spatial configuration and number of input points to correct for spatial outliers -and set the output cell size to ten meters. The map outputs were visualized using standard deviation (interval size 1/4 SD) with 13-17 classes depending on the number of mapped points. For each urban forest, we used these settings to visualize intensity for each disservice-category individually (n = 7). We removed mapped places of discomfort not being placed within the urban forest. Additionally, we removed mapped disservices being left on the default position since we assumed them having no spatial meaning.

Analyzing the relationship between cultural ecosystem services and disservices in urban forests
In order to find spatial bundles of CES and disservices, we created grids for three urban forests with cell sizes of 100 m × 100 m, resulting in a total of 16713 grids (Ungaro et al. 2016). Subsequently, we counted the number of markers per grid being placed for disservices and CES separately. We used Spearman's rank correlation to analyze the relationship between CES and disservices (Bagstad et al. 2018;Campagne et al. 2018;Nesbitt et al. 2019). We conducted all statistical analyses using R version 4.0.2 (R Core Team 2020). For visualization, we used circle-symbols be correlated in size with the number of mapped points per grid. We normalized values, ranging from zero to one for a better comparison of the total high number of services and lower number of disservices in one map. One outlier value (city of Stuttgart) has been visualized differently (grey circle) since referring all other values to the outlier would have strongly decreased visibility of all mapped positive and negative markers.

Respondents characteristics by study district
The extrapolated response rates of participants who had been invited via postal mail were 8% (district of Stuttgart), 10% (district of Karlsruhe), and 7% (district of Heidelberg), respectively. Given the expenses of rounded EUR 4000, for mail advertisement (shipping and printing), we had an imputed expenditure of EUR 5 for each response that had been acquired by mail. Overall, 1849 residents accomplished placing markers in the online-survey map (Stuttgart n = 798, Karlsruhe n = 630, and Heidelberg n = 421), and 755 of them also mapped places of disservices (n = 1552). 63% (n = 1173) of these participants completed all questions about their sociodemographic background. The average respondent spent 19 minutes to take part in the survey. The number of placed markers increased with the time they were engaged (CES, rs = 0.582; disservices, rs = 0.44). The mean age of our respondents was 47.2 years (SD = 15.17, range = 15-85), and 37% were female. Within our sample, the educational background was high, since 80% held a higher secondary school certificate. 74% of our participants completed a vocational training, and 62% held a university degree. For each of our study districts, most respondents were married (53% to 64%), employed (61% to 64%), lived in multi-person households (79% to 86%), and had monthly net household incomes between 2600 and 7500 Euro (52% to 54%) (Table A1).

The relative importance of cultural ecosystem services and disservices in urban forests
In the three urban forests, respondents mapped places of discomfort, making use of all given pre-defined categories. The marker 'place of discomfort' has most often been used (53%), followed by 'insecurity' (21%), 'blocked view' (6%), 'invasive species' (7%), 'restricted access' (5%), 'inconvenient plants' (4%), and 'allergy' (4%). Qualitative answers of respondents to the open question 'places I find uncomfortable for another reason' are summarized in Figure 2.

The impact of individual background characteristics on perception of disservices
The total number of positive and negative attributes did not vary significantly by either education degree, marital status, or income (Table 2). Similarly, the total of mapped attributes did not differ significantly by housing situation (single-person household versus multi-person household), work situation (working versus non-working), and expertise in urban forestry (expert versus layperson; Table 3). However, the total of mapped attributes varied significantly by gender. Female respondents mapped significantly more disservices than male respondents (on Table 3. Results of Wilcoxon rank-sum tests of sociodemographic effects on the total of mapped services and disservices.    (Table 3). Finally, using the middle 90% of the data, the time respondents spent for taking part in the survey varied by gender: female respondents spent significantly more time with the questionnaire than male respondents (Mdn = 1300 seconds versus Mdn = 1204 seconds, respectively; mean = 1932 seconds versus mean = 1739 seconds, respectively). However, the effect size was weak (W = 299,217, r = −0.050, p = 0.021). Spearman's rank correlation analyses revealed no significant relationships of age, perceived physical health, perceived psychological health, and place attachment to the total of mapped disservices per respondent (Table 4). In terms of mapped services, we found significant relationships of age and place attachment to the total of mapped services. Furthermore, we found significant positive relationships between familiarity with the respective urban forest as well as rated importance of urban forest design with the total of mapped services and disservices per respondent. Furthermore, we found a significant negative relationship of weak effect size between perceived life satisfaction and the total of mapped disservices per respondent (Table 4).

Spatial patterns of disservices
Hotspot maps of disservices derived from the kernel density estimates differed in their patterns for our three urban forests (Figure 3). With respect to the urban forest of Stuttgart, clustering of disservices was greatest there and hotspots were located at the center of the forest close to a lake that is being frequently visited within the urban forest. For the urban forest of Karlsruhe and Heidelberg, the hotspots were located towards the forest edges and closer to build infrastructure near the city centers (in Karlsruhe) and close to the surrounding communities (urban forest of Heidelberg) (Figure 1). Furthermore, hotspots were more dispersed in the forests of Karlsruhe and Heidelberg. The calculated Nearest Neighbor Ratios for the three urban forests were 0.548 in Stuttgart (z = −23.059, p = <0.001), 0.662 in Karlsruhe (z = −13.201, p = <0.001), and 0.661 in Heidelberg (z = −13.383, p = <0.001). In summary, mapped disservices were spatially clustered and not evenly distributed over the study area. The analysis of various types of disservices showed several spatial differences in the patterns of hotspots: Allergy (Stuttgart, Heidelberg), restricted access (Stuttgart, Karlsruhe), and blocked view (Heidelberg and Karlsruhe) were more concentrated on fewer hotspots in the urban forest and less distributed within the urban forest. In contrast, hotspots of insecure place and place of discomfort were more numerous and more distributed all over the study area in the three urban forests (Appendix A Figures A1-A3). Concerning the spatial distribution of services, human-made infrastructure as well as bluespaces both were most important predictors for the supply of services, whereas natural features and biophysical characteristics (habitat) were less important ).

Discussion
This study aimed at spatially identifying disservices that matter to urban forest visitors. We expected that the use of a PPGIS approach would reveal a more Table 4. Results of Spearman's rank correlation analyses of relationships between individual background characteristics and the total of mapped services and disservices. Additionally, we examined the spatial distribution of disservices to identify hot and cold spots of public interests within the study area. Finally, we looked at the congruence of cultural ecosystem services and disservices in urban forests.

Respondent characteristics
With regard to age, our sample is representative for the German general population, but it is not for the gender ratio (mean Germany = 44.4 years, 51% female) (Federal Statistical Office 2019, 2020b). Within our sample, the educational background was higher than the one compared with the general population of Germany, since 80% held a higher secondary school certificate (versus 34% in the German general population) (Federal Statistical Office 2020a). Similarly, 74% of our participants had completed a vocational training and 62% had held a university degree (versus 58% and 14%, respectively, in the German general population) (Federal Statistical Office 2020a). The higher educational status of respondents may be explained with the used method (webbased survey) as suggested by . Furthermore, people with higher education are more likely to visit urban forests (Mak and Jim 2019). For each of our study districts, most of our respondents were married, employed, lived in multi-person households, and had monthly net household incomes between 2600 and 7500 Euro. In comparison, the average monthly net household income of Baden-Württemberg's inhabitants is 3525 Euro (Eisenreich 2015). In the German general population, 43% are married (Federal Statistical Office 2020b), 58% live in a multi-person household, and the average monthly net household income is 3399 Euro (Federal Statistical Office 2019).

The relative importance of disservices in urban forests
According to our results, disservices that fit into existing EDS typologies are of lower relevance compared with other disservices resulting from people and their traces: people found anthropogenic noise being the most important disservice in the three urban forests. This is consistent with a previous study from Germany . Traffic noise -emitted nearby but not within the urban forest -was mapped more frequently than all other anthropogenic noises resulting from people within the forest (e.g. by social gathering or celebrating), most likely caused by large roads close to the urban forests (Stuttgart, Heidelberg) and the nearby city center in Karlsruhe with high traffic volume. High noise emissions may further stem from a racetrack close to the urban forest in Heidelberg. Traffic noise has previously been identified as important disservice in urban park settings (Liu et al. 2017(Liu et al. , 2018bSong et al. 2018). Garbage as well as dogs and dog owners were the second and third most stated disservice in our study. Garbage is frequently left in the urban forests, for instance, around officially declared barbecue places that are intensively used by the public especially during the summer. Dogs may be perceived disservice, as dog walkers follow a spatially clustered pattern that is similar to that of other user groups (Gerstenberg et al. 2020). This finding is also consistent with previous research (Agbenyega et al. 2009;Fish et al. 2016a). Finally, road conditions were ranked as the fourth disservice in our urban forests. Forest roads are typically broad, unpaved, and composed of loose gravel; they are also used for logging activities. Our findings concerning forest road dislikes are consistent with other research that focused on single trees (Delshammar et al. 2015) and peri-urban community forests (Rodríguez-Morales et al. 2020). Ultimately, places of crime and fear as well as health problems were ranked of minor importance in our study, which is inconsistent with previous research that emphasized scary places (Madge 1997;Skår 2010;Shackleton et al. 2015) and negative effects of nature on human health (Lohr et al. 2004; Campagne et al. 2018; Larson et al. 2019) as the most important disservices. Inconsistent results in terms of crime may be explained by the fact that surrounding cities in Germany have lower crime rates than surroundings of other presented study areas (Leicester, UK (Madge 1997) and Bela Bela and Tzaneen in Limpopo Province, SA (Shackleton et al. 2015)

Impact of individual background characteristics on perception of disservices
In line with some previous findings (Eriksson et al. 2012;Chen et al. 2016;Hegetschweiler et al. 2020), our study revealed that sociodemographic background characteristics do rarely affect perceptions of urban forest disservices. However, some other studies had identified age (Tyrväinen et al. 2003;Meyer et al. 2019), gender, and educational level (Tyrväinen et al. 2003) being influential. These findings are consistent with our study in terms of gender (services and disservices) as well as age (services). In fact, elder respondents mapped services more frequently than younger participants of our study. This is consistent with other studies (Ode Sang et al. 2016). Our finding that female respondents mapped significantly more disservices than male respondents may be associated with the fact that women spent more time with the questionnaire than men. Differences in terms of sociodemographic characteristics may be explained by different methodological approaches and different cultural and geographical contexts. Finally, the more familiar the respondents are with the forest, the more they mapped liked and disliked places.

Patterns of disservices and relations between cultural ecosystem services and disservices
A PPGIS approach may better satisfy people's demand for public participation in environmental decision-making and may give voice to the 'silent majority' represented by the general public (Brown and Reed 2009). Existing studies highlight that forestry planning maps are overall spatially well congruent with maps using PPGIS approaches, but the former are typically too general to guide place-specific forest planning needs (Brown and Donovan 2013). Maps such as those generated in our study may thus assist in place-specific forest management decisions. A less obvious but equally important benefit may be the restoration of public trust in state forest planning and decision processes (Brown and Reed 2009). Moreover, spatially visualized disliked places could encourage forest managers to either mitigate disservices or to promote CES on-site. The maps of unpopular places in urban forests that we present in this paper have the potential to advance existing planning instruments, such as forest function maps. Furthermore, presented maps of disservices can help local management to focus on disservices that matter. For instance, noise emitted by traffic and industry is in practice being addressed in forest function maps using exemplary measured and subsequently calculated sound emissions (Wirth et al. 2015;Palm and Wirth 2020). Our findings have the potential to complement this good practice by integrating spatially disliked soundscapes and in doing so focussing on people's perceptions as opposed to measured sound intensities.
In terms of spatial distribution of mapped places of disservices within the urban forests, we found hotspots of unpopular places. This is consistent with previous research that used transferable indicators instead of survey data to derive places of discomfort and calculate spatial congruence (Dobbs et al. 2014). Finally, we found CES and disservices to be spatially congruent. In fact, outcomes demonstrate overlapping clusters of perceived CES and disservices (bundles): locations where people referred the most positive meanings to were also the places where the most negative disservices were perceived. This is consistent with previous studies that also found statistically significant correlation of definite services with disservices Campagne et al. 2018;Rodríguez-Morales et al. 2020). In fact, the spatial congruence of perceived CES and disservices underpin the competing nature of both: the delivery of CES brings along disservices caused by recreational use intensity, for instance, garbage and noise caused by urban forest visitors. In turn, there were existing cold spots of CES and disservices with little public interest within the urban forest. This fact may indicate that people merely mapped places they are familiar with and as a result of this often visit. However, given that CES-and consequently also disservicesarise because of people's interactions with the ecosystem as described in section 1 introduction, there is no CES nor disservices without people (Potschin-Young et al. 2017).
In terms of identified cold spots in the urban forest, forestry practitioners may be encouraged to give priority to other than recreational objectives on these areas. For example, disturbance-sensitive wildlife could be promoted there or these sites could alternatively be managed with higher priority on timber production. Furthermore, our results enable forest managers to look at specific individual dislikes to consider them in planning and management. Finally, the high importance of disservices originating from people and their traces may indicate concentrating on the public -instead of focussing on natural environments. City foresters could, for instance, establish visitor group-specific zones in urban forests to ameliorate places of discomfort.

Limitations of the study and implications for future research
Elicited PPGIS-data using points instead of polygons may not account for the spatial extent of perceived disservices. For instance, poisonous plants have typically a low spatial extent, whereas noise can be perceived across large distances from the emitting source. Future research could develop specific buffers for different types of disservices, as suggested by a recent study (Koprowska et al. 2018). Furthermore, future studies could enable study participants to map the spatial extent by incorporating polygon features in the survey. However, existing studies that explored the validity of point versus polygon features for eliciting spatial perceptions of respondents indicate a strong degree of overlap of both when the number of collected point data is high enough (Brown and Pullar 2012). Our mapping approach required respondents to be familiar with maps as well as with the used software tool. Although the method has been often applied and validity of spatial data has been found to be moderate to high (Brown et al. 2017), mapped locations may lack accuracy. We faced this challenge allowing respondents to zoom in for mapping to gain detailed information on maps. Future research could develop approaches that overcome limitations in terms of orientation by providing small-scale maps using a definite zoom-in as a mandatory prerequisite. Revealed CES and disservices may not reflect the full range of perceived CES and disservices of respondents, since the number of placed markers increased with the time participants spent for taking part in the survey and they might have had lack of time. Future research could better represent people with lower educational and different cultural backgrounds by making more effort inviting them to participate in the study. This could be accomplished by other methodical approaches like face-to-face interviews or mapping on physical maps, since web-based survey methods were reported to be cognitively challenging ) and addressing younger participants with higher educational status (Pocewicz et al. 2012). Likewise, future research in other parts of the world may study disservices with other environmental settings.
The presented share of EDS as opposed to all perceived disservices may be too high since our seven markers included five pre-defined categories that are based upon existing EDS typologies and only two other disservices categories. Because of the timing of the survey (spring to summer), it is likely that perceived disservices in winter are under-represented. Consequently, effects of the daytime on perception of disservices might have influenced results of studies. Finally, survey data were collected before the global COVID-19 pandemic. Hence, one future research question may be the potential impact of the pandemic on perceptions of ecosystem services and disservices-such as perceived risks of vectorborne diseases.

Conclusions
In this study, we revealed that perceived dislikes in urban forests are rarely caused by ecosystem disservices in a narrow sense. Human-caused disservices are in contrast much more dominant. This finding may put the discussion of integrating disservices research into the ecosystem services framework in perspective. It also implies that city foresters should focus more on managing disservices caused by forest visitors instead of those provided by natural environments. Our visualization of hotspots of unpopular places in urban forests can complement existing planning instruments, such as forest function maps, as they are based on more reliable spatially explicit data. In most cases, forestry authorities do not have spatial data available that show places where people feel unpleasant. Known locations of disliked places may help forest managers to anticipate potential conflicts. We also found that cultural ecosystem services were much more important than disservices, which implies that urban forest visitors have largely positive perceptions when recreating in the forest. Our spatial analysis indicated a strong overlap of sites perceived for their services and their disservices: The delivery of cultural ecosystem services comes along with the drawback of perceived disservices on-site. Urban forest managers may thus guide the spatial movements of urban forest visitors within the forest more actively to avoid conflicts among various stakeholders. Finally, we identified cold spots of little public interest within the urban forest. Urban forestry practitioners may be encouraged to dedicate these sites to other than recreational objectives, such as timber production or conservation.

Disclosure statement
No potential conflict of interest was reported by the author (s).

Geolocation information
Urban Forest in the district of Stuttgart, Germany (WGS84): 48.765708, 9.099172 Urban Forest in the district of Karlsruhe, Germany (WGS84): 49. 084120, 8.433979 Urban Forest in the district of Heidelberg, Germany (WGS84): 49.337616, 8.584177 Table A1. Sociodemographic characteristics of study participants who placed markers of CES and disservices.  Figure A1. Hotspot maps of individual urban forest disservices (urban forest of Stuttgart), using kernel density estimates for point features (dark blue represents hotspots; point densities decrease from dark blue to bright blue; disservices n = 7, output cell size = 10 m, classification using standard deviation ¼, classes 13-15). Figure A2. Hotspot maps of individual urban forest disservices (urban forest of Karlsruhe), using kernel density estimates for point features (dark blue represents hotspots; point densities decrease from dark blue to bright blue; disservices n = 7, output cell size = 10 m, classification using standard deviation ¼, classes 15-16). Figure A3. Hotspot maps of individual urban forest disservices (urban forest of Heidelberg), using kernel density estimates for point features (dark blue represents hotspots; point densities decrease from dark blue to bright blue; disservices n = 7, output cell size = 10 m, classification using standard deviation ¼, classes 16-17).