Localizing and prioritizing roof greening opportunities for urban heat island mitigation: insights from the city of Krefeld, Germany

Climate change may increase the frequency, intensity, and occurrence of urban heat islands (UHI) in cities worldwide, often with harmful impacts on citizens. Strategic planning and implementation of multifunctional green roofs promises to help mitigating UHI effects, but cities often lack up-to-date scientific understanding of best-suited locations. The aim of this paper is to develop and apply a socio-ecological approach to explore and prioritize present and prospective opportunity spaces for roof greening based on remote sensing data to mitigate UHI effects. The city of Krefeld, Germany, serves as a case study. The research design consists of three steps, applied to the conditions of 2019 and a 2030 scenario: (i) Examining residents’ vulnerability to heat, (ii) Assessing existing green roofs and potentials for greening, and (iii) Prioritizing opportunity spaces for roof greening to reduce UHI effects. Findings showed that the area of high vulnerability due to combined high heat exposures and densities of sensitive residents in Krefeld accounts for almost 300 hectares in 2019 and may triple until 2030. More than 90% of evaluated horizontal roofs have no vegetation cover. Highest priority for roof greening is attributed to 59 ha and 113 ha of roofs in 2019 and 2030, respectively. The findings can inform strategic roof greening efforts for climate adaptation, e.g. for the extension of cadasters, and facilitate communication to increase understandings, public and policy support, and implementation.


Introduction
Urban Heat Islands (UHI), i. e. urbanized areas in which higher temperatures are experienced than in the surrounding rural areas, are considered as one of the major societal challenges of the twenty-first century (Rizwan et al. 2008). UHI are extensively researched and documented for many cities across the globe (Santamouris 2020;Tuholske et al. 2021). UHI arise due to shifts in energy fluxes (increase absorptions of solar radiation, increase of sensible heat and decrease of latent heat flow) associated with anthropogenic land cover changes and human activities (Oke et al. 1991;Tzavali et al. 2015). When surveying urban climate and the UHI phenomenon it is important to consider new paradigms as described by Wang (2022) not to oversimplify the complex system determined by interactions and changes along the urbanization gradient. The intensity of UHI varies between 0.4 and 11 °C (Tzavali et al. 2015), and is much intense at night (Azevedo et al. 2016), exposing urban residents to higher thermal stress. Such exposure can have adverse impacts for human comfort and health, causing increased mortality and morbidity, especially amongst low-income and sensitive population groups such as the elderly, children or sick people (Mücke et al. 2013;Huang et al. 2022). For instance, during the European heatwave in 2003, more than 70,000 additional deaths were registered (Robine et al. 2008). UHI-related mortality is expected to increase by climate change and urban growth (Chapman et al. 2017). Projections indicate that urban heat stress intensify further in the future and negative impacts of UHI might exacerbate (Patz et al. 2005;Huang et al. 2019).
A key strategy to mitigate against further increases in temperatures in urban areas and the associated negative impacts on human health is to 'green' urban areas by increasing the abundance and cover of vegetation (Bowler et al. 2010;Huang et al. 2021). Vegetation around and on urban structures, e. g. buildings, can affect its moisture, aerodynamic and thermal properties, and so urban greening could lower the air and surface temperature through different processes such as evapotranspiration, shading drop and resistance to air flows (Manso et al. 2021).
Roof greening, i. e. multi-layered systems that cover the roof of a building with vegetation over a drainage layer of soil (Berardi et al. 2014), presents a popular strategy to moderate the UHI effect as it exploits otherwise unused house roofs, lowers the air temperatures around urban green roofs up to 4.2 °C (Al-Kayiem et al. 2020), reduces building energy use by around 0.7% compared to conventional roofs (Sailor et al. 2011), and provides other multiple cobenefits including the regulation of air quality, carbon regulation, storm water management, habitat provisioning and cultural services (Berardi et al. 2014;Manso et al. 2021). The effectiveness of roof greening to moderate UHI effects depends on the area and greening type of the green roofs (Santamouris 2014;Sun et al. 2016). In general, green roof fraction should be expanded, as the increase of vegetated areas linearly reduces UHI (Li et al. 2014;Sun et al. 2016;Sharma et al. 2018). Especially by increasing the number of green roofs on buildings that hold high potential to reduce temperatures and positively affect surrounding areas, as for example office and industrial buildings that are expected to have a high share of large and flat roofs and low buildings heights and are not greened yet for about 97% in Berlin, UHI can effectively be reduced (Coenradie et al. 2016;Sen-Stadt 2016). Greening types are divided into extensive or intensive roof greening (Pfoser et al. 2013;Brune et al. 2017). Intensive greening generates higher cooling effects then extensive greening of roofs due to larger plants and thicker depth of planting mediums (Brune et al. 2017). However, extensive roof greening is more cost-saving, requires less maintenance and can be installed on a greater number of roof structures due to lighter weight (Getter et al. 2009).
Implementing roof greening is arguably particularly effective when planned at the neighborhood or city-scale (Peng and Jim 2013;Knaus and Haase 2020), however, substantial scientific knowledge gaps still exist. For instance, timely, cost-efficient, and spatially explicit datasets on structural and climatic conditions in urban areas are rare (Shao et al. 2021). This information is essential to derive local needs and potentials for mitigating UHI (Ansel et al. 2015;Wu and Biljecki 2021). Also, there is a lack of cost and time-efficient planning methods for conducting vulnerability analyses and investigating the greening potential of roofs to prioritize opportunity spaces and integrate them into urban planning (Coenradie et al. 2016;Wu and Biljecki 2021). Remote sensing could be a key tool to overcome these challenges (Wellmann et al. 2020). For example, remote sensing approaches have recently been used to investigate greening potentials and mitigate UHI (Karteris et al. 2016;Xu et al. 2020;Shao et al. 2021). Remote sensing approaches provide up-to-date datasets without costly and timeconsuming on-site mapping and facilitate data processing through automation of work steps (Ansel et al. 2015;Shao et al. 2021). Recent approaches use geographic information systems (GIS), remote sensing techniques and high resolution imagery to assess roof greening potentials. Skoryi et al. (2022), for example, estimated roof greening potentials from imagery that were classified in a semi-automatic process using the software of Ansel et al. (2015), instead of manual thresholding. Thresholds were set for reference polygons and used for classification of lager areas across several tiles by the software. However, Skoryi et al. (2022) found difficulties in defining suitable thresholds without misclassifications.
So far, studies research focus either on technical aspects of developing and improving methods for green roof inventory by investigating greening type and spatial distribution using imagery (cf. Hoffert and Lumasegger 2010; Coenradie et al. 2016) or on identifying and evaluating the potential for greening (cf. Karteris et al. 2016;Santos et al. 2016;Xu et al. 2020;Shao et al. 2021), but not on both aspects, although combined approaches are required for decision support and climate adapted urban planning (Santos et al. 2016).
The aim of this paper is to demonstrate and apply a socio-ecological approach to explore and prioritize opportunity spaces for roof greening. Data for the present (2019) and future (2030) situation in the city of Krefeld, Germany was analyzed, Krefeld, was selected as a case study because UHI occur frequently (Möller et al. 2020) and North Rhine-Westphalia was one of the most heat-affected regions in Germany during the 2003 European heatwave (Robine et al. 2008). The following three research questions were addressed: i. Which areas of Krefeld are characterized by high vulnerability to heat and thus require action to mitigate UHI? ii. Which roofs have already been greened or hold the potential for greening within the areas of high vulnerability to heat? iii. Which opportunity spaces for roof greening should be prioritized, taking into account vulnerabilities and cooling effects?
The remainder of this paper is structured as follows. Chapter 2 introduces the approach to prioritize areas for roof greening. In chapter 3 the approach is applied, and the results are presented for the case study of Krefeld. The following chapter 4 contains a critical discussion of the results and the applied approach as well as a reflection on the next steps to advance urban planning practices for an enhanced consideration of roof greening. The final chapter summarizes the findings of this study and provides recommendations for climate adaptation in Krefeld and the usability of the results for implementation of roof greening from a planning perspective.

Study site
The city of Krefeld (51°20′ N, 6°34′ E) is part of the Rhine-Ruhr metropolitan region and is located in the urban agglomeration of Mönchengladbach, Duisburg and Düsseldorf, which has an increasing impairment by heat with significant impact on health and infrastructure (Die Bundesregierung 2020) (Fig. 1). With around an area of 137 km 2 and 234,000 inhabitants, the city is a regional center in the metropolitan area (Stadt Krefeld 2020a).

Concept for prioritizing roof greening
A modular approach with three consecutive steps was developed to identify priority areas for roof greening that contribute to lower UHI for the situation in 2019 and in 2030 (Fig. 2). First, the vulnerability of urban areas to heat was analyzed by assessing the heat exposure of urban structures and the heat sensitivity of urban dwellers and facilities. Second, an inventory of green roofs and potentials for greening was compiled by analyzing spectral data from aerial images of roof greening with the Differenced Vegetation Index (DVI). Third, the results from the vulnerability analysis and the inventories were intersected to identify opportunity spaces with high priority for roof greening to reduce UHI in the city of Krefeld.

Assessing heat vulnerability
The vulnerability to heat was assessed for 3,848.6 ha of settlement areas in 42 statistical districts of the city of Krefeld ( Fig. 1) by valuation of the key indicators of heat exposure and sensitivity of the population (Table 1) (Wittig and Schuchardt 2012;Schulze Dieckhoff et al. 2018;HLNUG 2019).
The extent to which the population is exposed to heat depends on urban structures, e.g., type of building, building density, degree of sealing, and urban climatic characteristics such as reflection and heat storage properties, cold air generation and fresh air supply (Gao et al. 2022). Heat exposure was determined by using a climate assessment map that shows the thermal situation and spatial distribution of settlement areas with heat stress in Krefeld modelled by LANUV NRW (2018a, b). The map divided areas into five tiers of heat stress based on their level of heat intensity from low (tier 1 = very favorable thermal situation, physiologically equivalent temperature (PET) < = 35 °C and air temperature < = 17 °C) to high (tier 5 = very unfavorable thermal situation, PET > 35 °C and air temperature > 20 °C). Areas of tiers 1 and 2 were merged as they represent areas with a (very) favorable thermal situation that were of minor relevance for the vulnerability assessment. The resulting tiers were translated to heat exposure from Prioritization concept including three key steps and input data to estimate areas for roof greening very low (former tier 1 and 2 (very) favorable thermal situation) to high (former tier 5 = very unfavorable thermal situation). For the assessment of the prospective exposure, precautionary areas of climate change were used, which represent areas with an increase in heat stress by 2050 according to models of regionalized climate projections of North Rhine-Westphalia assuming a one Kelvin temperature increase between 2021 and 2050 (LANUV NRW 2018a). Precautionary areas of climate change were also grouped and translated into four tiers of heat stress (Table 1).
How sensitive people are to heat depends on sociospatial factors such as demographics and population density of the area. Heat sensitivity was assessed by using four sub-indicators representing population groups that are highly vulnerable to heat according to Bach et al. (2013) and HLNUG (2019). These population groups are the elderly (> 65 years), children (< 6 years), poor people (recipients of social and financial support according to the German Social Code (SGB II 2011/2021; SGB XII 2003/2021)), as well as ill and disabled people. Data from the population registry (Stadt Krefeld 2020a, b) and official spatial maps of settlement areas (Geobasis NRW 2019b) were used to calculate the population density of the heat sensitive population groups in settlement areas for the statistical districts of Krefeld (Table SM1). In accordance with HLNUG (2019), each heat sensitive population group was divided into five classes of population density by calculating the statistical quantiles (Table 1). Population data were not available for ill and disabled people; instead, the spatial locations of hospitals and facilities for people with disabilities (Stadt Krefeld 2019a, b) were considered for the further assessment. Hospitals and facilities for people with disabilities were buffered with a radius of 50 m and valued with high sensitivity based on the recommendations from HLNUG (2019).
To assess future heat sensitivity, population dynamics were modeled for the year 2030, using projections of population development for the elderly and children from Bertelsmann Stiftung (n. d.). The number of poor people in 2030 was calculated using a linear trend derived from data from 2009 to 2019, as no official projections were available. To estimate the population density of heat sensitive populations in settlement areas the same grouping approach was used as for the data of 2019. No changes were assumed for the number and location of the prospective development of hospitals and facilities for disabled people.
Analyzing roof greening and greening potentials Roof greening was analyzed based on the presence of vegetation on roof areas detected from multispectral information of high-resolution (0.1 m ground resolution), digital orthophotos (DOP) Geobasis NRW 2019c) by calculating the DVI for rooftops in areas highly vulnerable to heat in Krefeld in 2019 and 2030. The DVI is a graphical indicator that accentuates the reflectance of electromagnetic waves that are not absorbed by green plants (0.7-1.3 µm bands) by subtracting the visible red reflectance from the near-infrared reflectance. After calculating DVI for each DOP, values were classified to distinguish between areas without vegetation, extensive roof greening that appears as drier and less vital vegetation and intensive roof greening characterized by dense vital vegetation. The definition of thresholds for the classification depends on feature reflectance, which may vary between different imageries due to e.g., time of the day, phenology, incident radiation, cloudiness, shooting angle and differing spectral information of vegetation (Jones and Vaughan 2010;Pafi et al. 2016). Automated approaches for threshold definition are inaccurate (Xie et al. 2010), thus thresholds defining four classes were set individually for each picture by human visual interpretation of each imagery and histogram (Pafi et al. 2016) (Table SM2). DVI values that could not be clearly assigned to one of the three classes were placed in the group "uncertain". To analyze DVI for rooftops, DOPs were overlaid with outlines of rooftops in Krefeld from level of detail 2 (LoD2) 3D building models in city geography markup language (citygml) format (Geobasis NRW 2019a) by using the software FME Desktop 2019 Quick Translator (Safe Software 2021) for further processing in ArcGIS 10.6.1 (Esri 2021). All kinds of buildings with flat roofs and pent roofs were selected as these are not too steep, allowing greening without need for thrust lock and erosion protection (Landeshauptstadt Stuttgart 2010; Pfoser et al. 2013). The rooftop outlines were manually adjusted to reduce misclassification caused by tilting effects of the used DOPs and to fit the visible roof boundaries of the imagery (Jensen 2007;Coenradie et al. 2016;Marconcini and Metz 2016). Neglecting correction of the tilting effects can lead to an overestimation of greened areas by about 15% (Hoffert and Lumasegger 2010; Ansel et al. 2015).
Potentials for roof greening were derived from the results of the analysis of roof greening under the following assumptions: roof areas without vegetation were assigned to a high potential; extensive roof greening and pixels classified as uncertain were classified with medium potential; and intensive roof greening was assigned to a low potential. Estimates of the potential for roof greening in 2030 were made based on the assumptions that no changes in the city's DVI and roof inventory occur.
Prioritization of opportunity spaces for roof greening to reduce UHI For the prioritization of areas for roof greening the spatial results from the heat vulnerability assessment and the analysis of roof greening and greening potentials were intersected and evaluated using Arc-GIS Desktop 10.6.1 (Esri 2021). Roof areas highly affected by heat with a high or medium potential for greening were valued with the highest priority to reduce UHI, while roofs with low potential for greening were valued with medium priority.
Cooling effects of green roofs in Krefeld were assessed using the cooling rates identified by Li et al. (2014) according to the green roof fraction in Krefeld. They used a non-hydrostatic, regional climate model in conjunction with a physically-based urban canopy model to assess cooling impacts of green roof strategies and showed that as the green roof fractions increase, the near-surface temperature (two meter above ground) is reduced almost linearly. Estimated cooling rates on a city-scale varied between 0.0 and 0.59 °C depending on the coverage rate of green roof fractions (Li et al. 2014). To calculate changes of near-surface UHI cooling rates from Li et al. (2014) were attributed to the prioritized roof areas in respective to their share of the total roof area in settlement areas highly affected by heat in the 2019 and 2030 scenarios.

Results
The results of the 2019 heat vulnerability assessment showed that almost 90% of the assessed areas were classified as very low. High heat vulnerability was identified in 7% of Krefeld. These highly vulnerable areas are concentrated in the city center (Vier Wälle, Stadtgarten/Drietendorf, Südring, Bleichpfad, Stephanplatz and Schinkenplatz) and are characterized by a (very) unfavorable thermal situation with thermal stress due to high sealing rates and the highest density of elderly people in Krefeld (Fig. 3). Areas with a low heat vulnerability in 2019 were located in 2% of Krefeld, followed by medium vulnerability in 1%. These areas were identified in the southern parts of central Krefeld (Lehmheide and Dießem). The results for the estimation of the vulnerability in 2030 show that the areas affected by UHI will expand by 655 ha. Particularly, areas highly vulnerable to heat will triple in the future due to shifts from areas that have a very low degree of vulnerability in 2019 but are expected to be affected by climate change with temperature increase and consequential thermal stress. Due to changes in the heat sensitive population group of elderly people, with predicted growing numbers, vulnerability also increases (Fig. 3). These areas of high vulnerability will expand further outwards of the city, including districts south of the center and in eastern areas, as well as single districts further out (north-east Elfrath, east Uerdingen-Markt and Linn, as well as the suburb Hüls-Ortskern in the northwest). Areas with a medium level of vulnerability will almost double by 2030, whereas areas with low vulnerability will be reduced to one-third of their state in 2019. Medium and low vulnerability will in the future be found in small patches spread throughout the different districts of the city, except for the center (Fig. 3).
The majority of the investigated roofs, with 94% (57 ha), were not covered by vegetation in 2019 (Fig. 3). Intensive roof greening was found for two hectares in Krefeld. One hectare of the investigated roofs was covered by extensive greening. The results for 2030 showed with 93% (108 ha), the majority of the analyzed roofs remain uncovered by vegetation. Intensive roof greening was identified for a total of four hectares. Extensive greening covers a total of three hectares in Krefeld. The largest areas of green roofs (intensive and extensive) were identified in the Südring district (0.6 and 0.5 ha) in both 2019 and 2030. DVI values that could not be clearly assigned to one of the roof classes, without vegetation, extensive or intensive greening, covered one hectare in 2019 and 2 ha in 2030 in total for the city of Krefeld, with the largest area in the Südring district in 2019 and at Uerdingen-Markt (0.2 ha) in 2030.
The results of the prioritization of areas for roof greening in 2019 suggested that 49% (59 ha) of the total roof area in Krefeld that are located in areas highly vulnerable to heat have a high potential for roof greening and are therefore valued as high priority roof areas. Increasing green roof fraction by 49% would cause a cooling effect on the air temperature two meters above ground of up to approximately 0.2 °C in this area (Fig. 3) and contribute to the reduction of the near-surface urban heat island effect. In 2030, the ratio of opportunity spaces for roof greening remains similar to 2019. 44% (113 ha) of the total roof areas that are located in areas highly vulnerable to heat in Krefeld were classified as high priority roof areas (Fig. 3). Transforming these prioritized roofs into green roofs would change near-surface temperature by 0.2 °C. The largest high priority areas were found in the densely built-up inner-city districts of Südring and Vier Wälle (approximately 15 and 11 ha) in 2019 and 2030, respectively. Areas with medium priority for roof greening cover about 1.5% of the overall roof areas located in areas highly vulnerable to heat in 2019 and approximately 1.6% in 2030. The largest areas of medium priority for urban heat reduction in 2019 and 2030 were found in the Südring district (0.6 ha).

Discussion
This paper identified opportunity spaces for greening potentials of rooftops at the city-scale and suggested priority areas for further planning and implementation of roof greening measures with regard to a climate sensitive urban development. Although Krefeld is known as a garden city (Kuttler et al. 2003) with many green spaces (Fachbereich Stadtplanung 2015), the vulnerability to heat was already high in most areas of the city center in 2019. This is cause for concern as current assessments in Germany show that 92% of flat roof areas built in 2020 are not covered with vegetation (BuGG 2021). Furthermore, global phenomena such as urbanization (Kahlenborn et al. 2021), global warming (Die Bundesregierung 2020; Kahlenborn et al. 2021), and demographic change (Bunz and Mücke 2017) are driving local impacts, causing the expansion of UHI into the urban periphery and an increase of areas with high exposure and sensitivity to heat, as shown in the vulnerability assessment for Krefeld in 2030. The analysis resulted in 59 and 113 hectares (for 2019 and 2030, respectively) of opportunity spaces for roof greening in Krefeld, which is the amount of almost all flat roofs and pent roofs investigated in areas highly vulnerable to heat. These opportunity spaces can now be used in more detailed analyses of roof greening potentials, for example taking into account local experts' knowledge and information beyond available geodata.
The analysis shows a potential cooling effect of 0.2 °C for temperature in two meter above ground, if prioritized opportunity spaces would be greened. This decrease of near-surface UHI temperature could reduce heat-related mortality, potentially by up to 5% (Santamouris and Osmond 2020). However, these values represent modelled approximations of the reality and depend on highly dynamic feedbacks between, for instance, the ambient air temperatures and humidities, the surface energy balance as well as individual behaviour and health status of residents (Givoni 1991;Tong et al. 2021;Gao et al. 2022). Quantitative models of urban microclimate conditions and effects on human health are affected by uncertainties (Evola et al. 2017). Simulations on cooling effects due to roof greening and impacts on human health vary broadly by the selection of a model and associated set-up of parameters and indicators, as well as quality of input data and scenario design (Krayenhoff et al. 2021). Thus, results from the different methods used for simulations are hard to compare and should be interpreted carefully (Brune et al. 2017).
The proposed approach was designed for application in spatial planning practice and thereby needed to be relatively simple but robust enough for decision-making (Barbier et al. 2008). Several options for method improvement in each step of the presented approach were identified: Critique of methods applied for the assessment of vulnerability to heat The vulnerability assessment considered the location of settlement areas, but ignored the adaptive capacities of heat-sensitive population groups (i.e., the ability to prepare for, respond to, and cope with heat events), given the lack of suitable tested methods and indicators (cf. Heiland et al. 2012;Hartz et al. 2013;UBA 2021). Future applications could try to advance this, for example by drawing on insights from valuation approaches of people's adaptive capacity as those by Siders (2019), Buth et al. (2015) and Hemberger and Utz (2014). The analysis of heat sensitivity in 2030 was conducted under simplified assumptions such as a linear projection of the population development of poor people in Krefeld as well as a constant number and location of the prospective development of hospitals and facilities for disabled people. The estimation of prospective, structural urban changes can partly be improved by considering official landuse plans that define the development of future building areas. The projection of population development of poor people in Krefeld could be optimized, for example by using more complex models to achieve higher accuracy of forecasts, but would require more data and capacities (George et al. 2004). The vulnerability assessment could further be improved by also considering greening potentials of buildings' roofs in industrial areas. In particular, industrial buildings with large and predominantly flat roof areas at low building heights should be considered as these are promising in terms of cooling effects (Rittel et al. 2011;Pfoser et al. 2013;Santamouris 2014;SenStadt 2016).

Critique of methods applied for analysis of roof greening and greening potentials
Presence and absence of green roofs in Krefeld was determined by calculating DVI. There are alternative indices that can be used, such as the normalized difference vegetation index (NDVI), the color index of vegetation extraction (CIVE), and the chlorophyll vegetation index (CVI). However, the DVI turned out to be most useful since it is less influenced by atmospheric distortion, shadow, and background noise (Roujean and Breon 1995;Jones and Vaughan 2010) and is suitable to identify even sparsely vegetated areas (Broge and Leblanc 2001;Barati et al. 2011;Xue and Su 2017). The inventory of roof greening was done for a prespecified investigation area, which was classified as highly vulnerable to heat in the first step of the approach. In reality, the boundaries of vulnerability to heat are fluid and not strictly delineated (LANUV NRW 2018a). Results of the vulnerability analysis of Krefeld show clear boundaries that come from the input datasets and Fig. 3 Spatial changes of vulnerability to heat and prioritization of opportunity spaces for roof greening in 2019 and 2030. The alluvial diagram (top) shows the spatial changes (size of ribbons) of the vulnerability to heat in Krefeld (colored code of four classes: very low to high, in hectares). The maps show the results of the vulnerability assessment (polygons) and the analysis of roof greening and greening potentials (circles plotted scaled by 1.5) (middle map), and the results of the prioritization of areas for roof greening to reduce UHI in areas highly vulnerable to heat (bottom map). Background maps: statistical districts and DTK50 (Stadt Krefeld 2019c; Geobasis NRW 2021) ◂ were needed to select buildings for further spatial analysis and prioritization using GIS (Kuttler et al. 2013). Only buildings located completely within highly vulnerable areas were considered in the subsequent steps of analyzing roof greening and greening potentials as well as prioritizing opportunity spaces for roof greening. At the time of the study, only classic DOP were available for the identification and localization of green roofs. To achieve a higher accuracy in the localization process, true DOP should be used as these are corrected for image displacement caused by terrain relief and camera tilts as well as building tilting (Jensen 2007), thus reduce the time needed for correction of data from building outlines. The provision of true DOP in North Rhine-Westphalia is planned for 2021 (Bezirksregierung Köln 2020) and should be used in future research. A final error margin of green roof estimation was the accuracy of building polygons. For example, circa 0.4 ha of false polygon outlines of buildings that include vegetated inner courtyards were found. To reduce errors of building outlines in further studies, most recent data should be used and be proved by visual data inspection.
Furthermore, when analyzing roof greening additional roof types besides flat and pent roofs could be considered that are also suitable for greening. Those include, for instance, offset pent roofs, mixed roof types, or also more steeply inclined areas such as (flat) gable roofs. Most of the studies that investigated (potentials for) roof greening tend to underestimate the greening potential or the spatial extent of existing green roofs, because they only consider flat roofs (Ansel et al. 2015). However, the technical approaches today are more sophisticated enabling roof greening with slopes of up to 45°, if drainage and security measures are installed (cf. Pfoser et al. 2013;Brune et al. 2017). A further source of error related to inaccurate estimations of green roof area are DOPs distorted by different lighting conditions, e. g. gable roofs that are half shaded (Coenradie et al. 2016;Marconcini and Metz 2016). Additional steps are necessary to delineate and separately classify these shaded areas, for example through solar irradiation analysis as demonstrated by IP SYSCON (2020). The consideration of structural properties of buildings and rooftops could also be improved, for example regarding constructive and static characteristics of the buildings, accessibility of the roof areas, the overall building condition, and the duration of its planned maintenance and financial scope (Ansel et al. 2015;Cascone et al. 2018). Regarding the analysis of building statics for retrofitted greening, Ansel et al. (2015) provided a manual on how to estimate the bearing load of flat roofs from aerial photographs and identified roofs with potential for extensive greening.
Critique of methods applied for prioritization of opportunity spaces for roof greening to reduce UHI Prioritizing opportunity spaces could include additional spatial indicators that influence the effectiveness of UHI mitigation, such as building height and urban design in terms of the ratio of building height to street width (Ng et al. 2012;Peng and Jim 2013;Santamouris 2014;Sun et al. 2016). These indicators can be derived from LoD2 building models, and used as an input to enhance estimations of potential local cooling effects from roof greening on the neighborhood or city-scale (cf. Rosenzweig et al. 2006;Peng and Jim 2013;Santamouris 2014;Sun et al. 2016). Also, a more integrative analysis of the multifunctional effects of green roofs may influence the prioritization of opportunity spaces (cf. Semeraro et al. 2021). Green roofs provide co-benefits that contribute to address multiple development goals, but also disservices that affect human well-being negatively such as allergenic pollen from roof plants or deteriorated water quality (van Seters et al. 2009;Gregoire and Clausen 2011). A framework to assess the co-benefits, trade-offs, and unintended consequences of green roofs is presented by Wang (2021). As a next step for integrating the study's results into planning and practical implementation of greening measures, pixelbased results should be aggregated to entire roof areas of buildings or building parts (cf. Ansel et al. 2015). Therefore, the classified pixels of a roof need to be reclassified to the vegetation class that holds the largest share within the respective building outline. From a practical point of view, it is furthermore recommended to set a threshold for minimum areas suitable for roof greening, as done by Ansel et al. (2015). The generated information should also be fed into planning discourses by utilizing multi-criteria decision analysis, as suggested -among others -by Langemeyer et al. (2020) and Venter et al. (2021) to evaluate potentials and prioritize areas for roof greening.
The transferability of the three-step approach depends on the availability and quality of the input data and would benefit from more automation of individual steps. For instance, the analysis of aerial images to identify green roofs should be performed with automated processes rather than human interpretation. To do so, automatic or semi-automatic classification methods should be tested to reduce the time required for manual classification through e.g. machine learning approaches (Shao et al. 2021). For example, Chen et al. (2022) proposed a segmentation decoder for the High-Resolution Network (HRNet) model that accurately detects urban green on a fine scale with improved performance on shaded green.

Conclusion
The approach to identify roof areas with potential for greening and priority for UHI reduction presented in this paper was developed through the modular combination of socio-ecological and remote sensing geodata in three consecutive steps. The results show that the approach can be applied to estimate opportunity spaces of roof greening. The opportunity spaces can subsequently be further assessed with more detailed data and input from local experts. The opportunity spaces also form a basis for deliberation in support of the creation of urban climate adaptation strategies. The Krefeld case study highlighted the urgency to respond to the negative impacts of urban heat, which will become even more critical in the future. While the settlement area identified as highly vulnerable to heat in 2019 was about 273 ha, the results of the 2030 scenario showed a tripling in size of the area of high vulnerability, and an expansion towards the outskirts. The majority of the investigated roofs (> 90%) were not covered with any vegetation in areas with population groups highly vulnerable to heat. Roofs located in those areas should therefore be considered in more detail for the implementation of greening measures, as they hold substantial potential for greening and lowering the UHI. According to the estimates and based on the assumption that all rooftops identified with high priority in Krefeld were greened, the near-surface UHI temperature could be mitigated by 0.2 °C and the heat-related mortality reduced by 5% (in both 2019 and 2030). The application of the presented approach provides valuable information that complements the data basis in green roof cadasters, as these were already established in the form of online GIS-applications in urban and rural areas in Germany to provide information for planners, decisionmakers and citizens (Die Senatorin für Klimaschutz, Umwelt, Mobilität, Stadtentwicklung und Wohnungsbau 2019; Landratsamt Schweinfurt 2019; LANUV NRW 2021). The information provided in green roof cadasters varies in quality and quantity of the data. For example, in the recently published cadaster of the federal state of North Rhine-Westphalia (LANUV NRW 2021), information on potential for roof greening, storm water retention, fixation of carbon dioxide and dust binding is available. However, a documentation of existing green roofs and population's sensitivity to heat is not yet considered in the estimated potentials for roof greening of the cadaster. Complementing the cadaster with place-based socio-ecological information from the proposed approach can contribute to increase the relevance for practical usage. A more integrative cadaster may better educate citizens regarding UHI and promote greening strategies to cope with negative impacts of urban heat. The presentation of further ecosystem services provided by green roofs and their co-benefits highlight the added value of greening measurements. Identified and presented synergies can contribute to the realisation of multiple development goals from city scale (e.g. citys' climate concept as e.g. "KrefeldKlima 2030") to national scale (e.g. National Climate Adaptation Strategy) (GERTEC 2018;The Federal Government 2008). By adding the information on synergies from roof greening on a city-scale in the cadaster can support the visibility, communication, promotion, and ultimately the implementation of measurements. This could be supported, in particular, by complementing cadasters with monitoring of roof greening. Implementing roof greening on public buildings can present insightful best practice examples that may serve as role models for rapid outscaling (MULNV 2015;difu 2018). Developing and implementing strategies for roof greening in practice needs to consider the advantages and limitations of such approaches. As discussed in this paper, substantial spatial opportunities exist to harness green roofs for climate adaptation. However, planners and decision-makers need to consider the limited direct effects of green roofs on the thermal environment at the pedestrian level as experienced by citizens. In response, planners should consider the broad spectrum and diverging requirements, cooling effects, and co-benefits (Gunawardena et al. 2017;Han et al. 2022) of green and blue infrastructure actions available to maximize cooling such as tree planting, artificial wetlands, and green walls within integrated development strategies. If planned and implemented comprehensively across local to city scales, such coordinated efforts could yield substantially positive health effects (Iungman et al. 2023).
Future research may focus on further advancing the proposed approach with more detailed methods and automated interpretations. It could also explore the applicability of the approach in other data, climate, and governance contexts. In addition, further research could try to assess the effects of roof greening on other benefits and costs, for example in terms of impacts on biodiversity and recreation as well as on reduction of energy consumption for buildings cooling and heating. Future assessments of greater detail and broader scope could thereby provide more comprehensive accounts of the opportunities and implications of roof greening and allow citizens and decision-makers to make better-informed decisions.
Author contributions All authors contributed to the study conception and design. Methodology was developed by JB with guidance from SS. Material preparation, data collection and analysis were performed by JB. The original draft was written by JB and supervision, writing -reviewing and editing of the manuscript was done by SS and CA. All authors read and approved the final manuscript.
Funding Open Access funding enabled and organized by Projekt DEAL. The work was realized within the PlanSmart Research Group funded by Grant 01UU1601 A and B from the German Federal Ministry for Education and Research (Bundesministerium für Bildung und Forschung-BMBF).

Competing interests
The authors have no relevant financial or non-financial interests to disclose.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.