Applying GIS in Blue-Green Infrastructure Design in Urban Areas for Better Life Quality and Climate Resilience

: The expansion of urban centers and peri-urban zones significantly impacts both the natu-ral world and human well-being, leading to issues such as increased air pollution, the formation of urban heat islands, and challenges in water management. The concept of multifunctional greening serves as a cornerstone, emphasizing the interconnectedness of ecological, social, and health-related factors. This study aimed to identify potential locations for three specific types of blue-green infrastructure (BGI): bioswales, infiltration trenches, and green bus stops. Leveraging geospatial datasets, Geographic Information System (GIS) technology, and remote sensing methodologies, this study conducted a comprehensive analysis and modeling of spatial information. Initial cartographic representations were developed to identify specific locations within Olsztyn, a city in Poland, deemed appropriate for the implementation of the designated blue-green infrastructure (BGI) components. Following this, these models were combined with two additional models created by the researchers: a surface urban heat island (SUHI) model and a demographic model that outlined the age structure of the city’s population. This synergistic approach resulted in the development of a detailed map, which identified potential locations for the implementation of blue-green infrastructure. This was achieved by utilizing vector data acquired with a precision of 1 m. The high level of detail on the map allows for an extremely accurate representation of geographical features and infrastructure layouts, which are essential for precise planning and implementation. This infrastructure is identified as a key strategy for strengthening ecosystem resilience, improving urban livability, and promoting public health and well-being.


Introduction
The important and urgent need to cope with climate change, coupled with the evolution of urban and peri-urban landscapes, emphasizes the need to seek solutions focused on mitigating and, in retrospect, counteracting the effects of this situation.Consequently, methods for monitoring the environmental impact of human activities on the urbanization processes taking place are increasingly important [1][2][3][4][5][6][7].Multifunctional greening is fundamental to providing a wide range of ecosystem goods and services that benefit the urban population.The motivation for undertaking this research was the prevailing opinion in the literature that there are not many methods and approaches to quantify and monitor multifunctional greening at different spatial scales also with GIS usage [8][9][10].Many of the urban areas that are elements of blue-green infrastructure are undergoing development, thus maintaining a minimal amount of green space, for which the type and area are mainly the result of existing urban planning regulations rather than a necessity arising from minimizing future problems [11][12][13][14][15][16].Currently, geospatial data are the most meaningful type of data and thus a key element for increasing the scope of information about the state of space.Ongoing research on blue-green infrastructure provides examples of the use of GIS tools in spatial analysis [17][18][19][20].
The integration of blue-green infrastructure in urban areas is proposed as a strategy to lessen the impact of escalating temperatures from climate change on health outcomes [11,[21][22][23][24][25][26][27].Highlighting the wealth of research confirmed the positive effects of BGI; the focus has been directed toward elements seamlessly connecting water and natural features, specifically bioswales, infiltration ditches, and green bus stops [28][29][30].
The BGI encompasses design and spatial solutions rooted in the intrinsic qualities of a specific location, often referred to as nature-based solutions (NBSs) [31].The elements of BGI are designed to create a cohesive system where services and spatial configurations enhance and support each other reciprocally [18].In urban environments, BGI serves multiple purposes concurrently.It can store and purify water, enhance the visual appeal of an area [32], absorb carbon dioxide, mitigate air pollution, and counteract the urban heat island effect by regulating the air temperature.Furthermore, these elements provide habitats for urban plants and wildlife, fostering ecological continuity and contributing to environmental education.When properly designed, they mitigate excessive surface runoff and the risk of flooding [33][34][35].BGI contributes to improving mental health and positively influences the well-being of city dwellers by mitigating the negative impacts of climate change with the provision of cohesive ecosystem services [31,36,37].
The technology under discussion ensures the effective collection, storage, and handling of spatial data, enabling the processing and examination of these data to generate novel, dependable spatial information, namely, details associated with a specific spatial location.Some research in this arena focuses on employing the GIS platform to recognize established blue-green infrastructure (BGI) [38,39] as well as whether there are enough BGI elements to ensure at an appropriate level the quality of life of the population and whether this is in line with current legal regulations [15,[40][41][42].According to the literature, there are several difficulties in measuring quality of life, mainly due to the many interpretations of the concept and the lack of a complete, agreed definition [3,43,44].This paper assumes that the concept of quality of life is defined by a model that combines objective and subjective indicators, encompasses different domains of life, and integrates both personal autonomy and individual preferences with actual assessments.This approach suggests that a statistical measurement of quality of life should consider the multifaceted nature of the concept.It proposes that such measurements should include both objective conditions of an individual's life and their subjective experiences or subjective well-being.In conclusion, the term living conditions encompasses a multitude of factors, including material conditions, health, education, economic activity, leisure activities, social relations, personal security, and the quality of the natural environment at the place of living.The multifaceted domains addressed by BGI, as previously mentioned, effectively translate into both objective improvements in living conditions and subjective perceptions of well-being.A holistic approach that takes into account the optimal realization of BGI is in line with the emphasis in the literature on including both tangible and intangible aspects of life in quality of life assessments.Consequently, BGI not only supports the technical functioning of a city but also enriches the lives of its inhabitants by creating healthier, more sustainable, and enjoyable living environments.
The referenced studies demonstrate that Geographic Information System (GIS) tools are crucial for identifying areas in need of further development to ensure sustainability and for effectively selecting components of blue-green infrastructure.Additionally, substantial research focuses on leveraging remote sensing data to analyze land cover change dynamics, especially in scenarios characterized by swift and significant urbanization [45][46][47].Nevertheless, a unifying theme emerges from these studies-the automation of conducted analyses and the capacity to execute them over expansive geographical areas, objectives that this study also espouses.It is also noteworthy that spatial analyses developed using GIS permit the precise identification of phenomena and processes occurring within the context of dynamic suburban development [48][49][50][51][52][53].
The form of urban spaces, as well as areas under pressure from urbanization processes, have been permanently shaped over the years.Thus, it should be noted that it is much easier to design new elements of blue-green infrastructure in open and nature-rich areas than in an area where everything is already occupied, and such spaces are difficult to find.Given the location of the research area, this study aligns with the scientific article [54] in examining information relevant to BGI from the standpoint of its inception, particularly in zones where recreation and tourism are of central importance.The feasibility of performing analyses and spatial data modeling was validated with a GIS, which also facilitated the evaluation of how practicable it was to implement three selected BGI components.Consequently, a map was produced showcasing potential sites for three specified types of blue-green infrastructure within the urban setting, encompassing bioswales, infiltration trenches, and green bus stops.
This research methodology utilized geospatial data, including remote sensing data supported by GIS technology, to identify and optimize the location selection process for BGI (blue-green infrastructure) components.The results obtained confirm the validity of the use of specific spatial datasets as well as the applied spatial analysis methods in improving the location optimization process for different BGI components.

Multifunctional Greening and Blue-Green Infrastructural Components
The implementation analyses of BGI are intended to assist planners and landscape architects in integrating nature-based solutions into urban planning, serving as either an alternative or a complement to traditional infrastructure.These activities, based on the principles of sustainable development, aim to improve the quality of life and protect the environment.It is also essential to conduct research on effective land management, which is crucial for sustainable urban development and the development of rural areas [55,56].
These ecosystems offer critical services, essential for tackling climate change through mitigation and adaptation efforts, notably including [57,58] cooling and insulation; the absorption of CO2; the utilization of low-carbon materials; and the promotion and implementation of sustainable development goals (SDGs) [59].It also aims at increasing the continuity of natural areas within cities, supporting their ecological role and substantially increasing urban biodiversity [60].
Subsequently, the specific BGI components that are most suitable for the optimal use of urbanized areas are delineated, either independently or in combination.The utilization of BGI elements permits the enhancement of environmental and thermal comfort, whereas the deployment of multiple elements would have a more synergistic impact on the urban environment [4].The integration of the elements collectively or individually is contingent upon the specific urban environmental objectives, the availability of space, and the specific climatic conditions of the area.The combined use of the BGI elements analyzed ensures the adequate irrigation of the vegetation and, simultaneously, reduces excessive evaporation through appropriate shading [61].Concurrently, each of these elements contributes to the improvement in water quality through filtration, which has a positive impact on soil quality.Furthermore, the absorption of water directly into the ground serves to reduce the load on urban drainage systems, preventing the flooding of urban areas.At the same time, the introduction of appropriate vegetation and a shading function significantly enhances urban aesthetics, making it more welcoming and attractive for residents.
Table 1 presents the concept and configuration, highlighting the spatial attributes favorable to their installation and offering insights into the proposed locations for the various BGI components [57,[62][63][64][65][66].Concurrently, this study aims to collate a set of spatial features (geodata), which represents a fundamental preliminary step for subsequent analyses to identify optimal locations for these components [57].In the context of this study, the researchers examined the spatial attributes identified in the referenced literature as key to identifying optimal locations for BGI elements [57,62,[64][65][66][67].The analysis conducted unveiled that, with regard to specific features, publicly accessible spatial datasets exhibited variations in accessibility and quality.This allowed for the identification of parameters from the aforementioned table that could contribute to establishing a composite indicator for determining the optimal locations of individual BGI components.It is notable that, in the analysis, the decision was taken to primarily rely on publicly available registries that are distinguished by high quality, timeliness, and standardization.Consequently, the chosen attributes used in the spatial analyses enabled the precise identification of locations and, more critically, ensured that these analyses could be replicated in future years.

Geodata Analysis for Optimum BGI Locations
This study focuses on the city of Olsztyn, covering an area of 88.3 km 2 .Forests constitute 21.3% of this area, with the City Forest alone comprising 18 km 2 .The majority of the 5.6 km 2 of urban green spaces in Olsztyn are parks and cemeteries.The city also boasts eleven lakes and five smaller water bodies, predominantly located in the western part.Four rivers, the Łyna, Wadąg, Kortówka, and Skanda, flow through the city.With a population exceeding 172,000, Olsztyn provides each resident with an average of 139 m 2 of urban forest, 3.5 m 2 of parkland, and 43.5 m 2 of water within the city limits.Despite abundant natural resources, the increasing prevalence of residential and service-oriented development poses a threat to green spaces, often leading to their usurpation or degradation.This trend is confirmed by the heat island map and analysis conducted for this study.The necessity to identify optimal sites for blue-green infrastructure (BGI) deployment is highlighted, aiming to enhance urban sustainability.The developed method assists in selecting sites with favorable spatial characteristics for BGI, ensuring their viability and effective integration.The optimal BGI location map is instrumental in the decision-making process, revealing potential sites that traditional methods might overlook.
The geoinformatic analysis was performed to develop a cartographic model of the optimal location of BGI following the scheme presented in Figure 1.The spatial data were categorized into three sets (Figure 1):

•
Group A-refers to forms of land use that preclude the optimal location of the BGI elements analyzed.Surface waters, buildings, roads, communication zones, railways, cemeteries, sports facilities, and forests were included in this category (Figure 2).The BDOT10k database, which provides the required accuracy and timeliness of data, was used to identify these areas [68].

•
Group B-relates to elements that determine the optimal placement of the analyzed types of BGI.The geospatial features listed in Table 1 and the datasets that allow for their identification that are conducive to the occurrence of bioswales, infiltration trenches, and green bus stops are mentioned later in this paper.
• Group C-spatial datasets that allow for the identification of urban heat islands (SU-HIs) on the analyzed urban area and the analysis of the distribution of residents' places of residence, taking into account their age (DM-Demography Model).The BDOT10k and DTM datasets were used for this analysis.The BDOT10k database, a vector-based repository, contains details of the location and characteristics of topographic features, facilitating the production of 1:10,000 scale maps.The DTM dataset, on the other hand, is a point model of terrain elevation created using aerial laser scanning (ALS).This research used a grid model with a 1 m × 1 m mesh in the PL-KRON86-NH height reference frame to represent the terrain elevation.These datasets were chosen for their clarity, ease of acquisition, and completeness in relation to the area analyzed.Concurrently, this study proposed the utilization of the vector data from OpenStreetMap as a supplement to the topographic data, which may be more up to date in certain urban areas.However, it should be emphasized that currently the BDOT10k dataset is the main source of complete and high-quality data at a national level.10mar The use of these different datasets and their analysis with specialized tools allowed for the identification of methods to determine the optimal locations for the analyzed BGI elements, taking into account the specificities presented in Table 1.

MBGI for Bioswale, Infiltration Trench, and Green Bus Stops
The first BGI element subjected to a geoinformatic analysis was the bioswale.In examining this BGI component, the spatial analysis considered factors such as the areas with slopes of less than 5%, flood-prone areas, the minimum and maximum catchment sizes, and the maximum planned width of the bioswale.In addition, using the favorable characteristics analyzed, the best locations for the bioswale were identified near roads, footpaths, cycle paths, car parks, and public areas.These elements were identified as spatial characteristics that would support the integration of the specified BGI component.
The analyzed facilities required analyses of the downslopes to be first carried out.In accordance with the adopted assumptions concerning BGI (Table 1), areas with slopes of less than 5% were sought.Next, the BDOT10k dataset was used to remove surface water and woodland areas, as well as building and road areas, from areas considered to have the required slope.At the same time, buffer zones of 30 m for bicycle paths and existing roads were established, which would enable the indication of areas with slopes of up to 5% that should be considered in the analysis as the preferred areas for bioswales.The next step assumed, based on the SCALGO and DTM application algorithms, the development of runoffs in vector form for the area under analysis.It was then assumed that the BGI component, i.e., the bioswale, should be established to maximize its function on the existing runoff lines.At the same time, the maximum width of the trench was adopted to be 5 m in accordance with the above findings (Table 1).This was followed by the identification of common areas (intersect) for specific trenches with a width of 5 m and the areas located in the vicinity of roads and paths as well as squares and car parks.In consideration of the fact that bioswales are also designed to impede surface runoff, the analysis included areas that had been flooded with a minimum of 10 mm of rainfall.The analysis, extended to include the possibility of considering the function of minimizing the risk of flooding in the city area, allowed for the location for the BGI component concerned to be made more specific.Finally, based on the analysis concerned, an MBGI cartographic model for bioswales was developed, which helped identify 1615 optimum locations (Figure 3).
Another component of blue-green infrastructure that was analyzed in this study was the infiltration trench.Due to the similar nature, it was adopted that the spatial analysis conducted could also be applied to determine the optimum location of rain gardens in a container (bioretention planter).The evaluation of the spatial characteristics of this BGI component took into account factors such as the largest catchment areas, the planned maximum width of the component, and the local drainage systems.At the same time, the favorable characteristics identified in the analysis were used to identify prime locations in the vicinity of playgrounds, sports complexes, recreational areas or open communal spaces, and parking facilities.These elements were recognized as spatial attributes that facilitate the integration of the chosen BGI component.
The first step assumed the determination of a drainage basin with a maximum overall area of 5 ha within the area of the city under analysis.As regards the infiltration trench, as was performed for the bioswale, the basic analysis was to determine the runoff lines over the entire area under analysis.Next, the drainage lines that passed through areas covered by surface water and forest, as well as areas where buildings and roads were located, were removed from further analysis using the BDOT10k data.At the same time, the maximum width of the infiltration trench was adopted to be 2.5 m in accordance with the adopted findings (Table 1).This was followed by the indication of common areas (intersect) for trenches with a preset width and recreational and sports areas, squares, and car parks, which were extracted from the BDOT10k database.As a result of the analyses, 62 locations for the BGI components, i.e., infiltration trenches, were identified and then juxtaposed with the layer representing a drainage basin with a maximum area of 5 ha.Finally, 39 optimum locations were selected (Figure 4), which enabled the generation of an MBGI cartographic model of the locations of infiltration trenches in the Olsztyn city area.The last blue-green infrastructure component analyzed in this study was the green bus stops.In examining the spatial aspects of this BGI feature, considerations included the location of public transport stops, the maximum size of catchment areas, and the flood risk zones within the area.These considerations were identified as spatial attributes favorable to the deployment of the selected BGI feature.
In order to inventory bus stops in the city area, the BDOT10k and OSM databases were used.However, due to the lack of data or incorrect location of facilities, there was a need for verification and updating, which was carried out using an up-to-date orthophotomap.The above-mentioned measures resulted in the establishment of a database containing information on the location of 369 bus stops.According to the adopted assumptions (Table 1), in relation to the bus stops, the area under analysis covered a maximum area of 60 m 2 .The data compiled in the SCALGO application, representing areas that are subject to flooding following significant daily rainfall of at least 10 mm, were then used.When implementing the additional function of the green bus (Table 1), which concerns the reduction in the risk of local floods and overloading of the storm drainage system, the bus stops, whose drainage basins overlapped spatially with flooded areas, were selected.The analysis resulted in the selection of 39 bus stops, which represented a cartographic model (MBGI) of the optimum locations of green bus stops (Figure 5).

Compilations of Surface Urban Heat Island (SUHI) Model and Demography Model (DM) Applied to the Study Area
Urban heat islands were identified by measuring the land surface temperature (LST) using data from the Landsat 8 mission.Following the recommendations of the United States Geological Survey (USGS), thermal infrared sensor (TIRS) data were used with a spatial resolution of 30 m for channel 10 [70].
The chosen methodology, necessary for accurate measurements, included atmospheric correction and consideration of different emissivity according to the nature of the surface.The single-channel algorithm method [71] was used, which includes atmospheric functions, atmospheric correction, and a threshold method for the normalized vegetation index difference.This method actively incorporates emissivity thresholds using the NDVI in the measurement process [72].Given its proven effectiveness in numerous studies, especially on Landsat-8 OLI/TIRS mission data, the methodology described in this study was considered suitable [71,[73][74][75][76].It is worth noting that, of the currently developed methods, the one by Sobrino et al. [77] stands out for having the most accurate results.
In view of the accuracy and trueness of the analyses performed, historical meteorological data and the maximum cloud cover of 12% were taken into account in the basic conditions for the selection of the available data.The analysis covered the summer periods from the beginning of June to the end of August from the years 2021-2022.Based on the indicators analyzed as part of the procedure for calculating the actual temperature of the area under study, the datasets developed for five days, i.e., 17 June 2021, 12 July 2021, 13 August 2021, 12 June 2022, 28 June 2022, and 14 July 2022, were analyzed (USGS, no date).The next step compared the maximum temperature of a particular day based on archival data [78].
Based on the satellite observations over the examined area, the highest temperature noted at 10:00 was chosen for comparison.Historical meteorological records revealed that the peak temperature, reaching 30.9 °C, occurred on 28 June 2022 [78,79].This resulted in the data collected on that particular day being selected for detailed analysis (Figure 6).In order to analyze the age structure within the study area, demographic tables from the Department of the Olsztyn City Council were used.These data were compared with the boundaries of the residential areas after they had been collected and collated.This approach made it possible to segment the population into defined age groups, 0-10, 11-65, and over 65, across different residential areas.The results of this segmentation are shown in Figure 7, which highlights the regions with a higher concentration of people aged 65 and over who are more susceptible to the adverse effects of elevated temperatures, particularly in the summer months.

Results
In the context of the factors considered when determining the optimal locations of the specific BGI elements, the researchers emphasized the possibility of utilizing diverse data sources, including digital terrain models (DTMs), a topographic database (BDOT10k), and hydrological data.The utilization of the digital terrain model allowed for obtaining information about the terrain's morphology, which is crucial for analyzing terrain slopes, groundwater levels, and surface water runoff directions.On the other hand, the BDOT10k topographic database enabled the identification of both preferred areas, such as parks, areas along roads, and parking lots, suitable for BGI components, as well as unsuitable areas, encompassing forests and areas occupied by existing buildings or infrastructure.The last analyzed aspect was related to the hydrological data and analyses, which resulted in the creation of watershed and surface runoff maps.Determining the flow directions of water and identifying areas prone to flooding supported the selection of the BGI component locations and contributed to flood risk reduction and the mitigation of the negative impact of extreme rainfall events.
In the following step of this investigation, comparisons were made between the MBGI cartographic models, depicting the best locations for the blue-green infrastructure components, and the findings from both the surface urban heat island (SUHI) analysis and the demographic model (DM), as presented in Figure 8.The integration of the different mapping models involved aligns the individual datasets containing information on optimal BGI locations with the temperature data, the age structure of the population, and the geospatial factors that preclude the location of the individual BGI elements.By extending the focus to health challenges, the aim of this phase was to propose targeted solutions in specific zones of the urban study area.As a result, the MBGI mapping models developed became a valuable tool for visualizing, analyzing, and selecting the optimal locations for the BGI components.The aim was to maximize their positive impact on temperature reduction, which translates into more favorable living conditions in the city.It is important to note that the impact of BGI on urban microclimates may vary depending on the specific implementation and environmental context [80].Variability due to geographical and climatic factors should inform the design and implementation of BGI elements, taking into account local conditions.Climatic differences, such as the increased cooling and humidification effects of water bodies in hot and dry climates compared to temperate regions, are important to consider.Similarly, vertical greening may be more effective in urban areas with high solar exposure.The structure and layout of a city also have a significant impact on the effectiveness of BGI.In densely built cities, small enclaves of green infrastructure can provide the necessary relief from heat, while in more dispersed cities, larger and more continuous green spaces may be needed to achieve similar effects.The cultural and social context also plays a role in influencing the use and maintenance of BGI.The utilization of public green spaces exhibits considerable variation between continents, contingent upon the extent of public involvement and municipal support.Economic factors, such as the financial resources available to implement and maintain BGI, also influence its quality, scope, and effectiveness.Cities with greater financial resources are able to implement more comprehensive and technologically advanced BGI systems.Consequently, city planners and city managers should adapt BGI strategies to specific local conditions.It is of the utmost importance to conduct comprehensive climate surveys prior to the implementation of BGI, to analyze community engagement, and to continuously monitor and adjust BGI projects based on their performance and community feedback in order to identify key conditions and optimize projects.Minimizing the negative effects of high temperatures is directly linked to reducing the risk of escalating hazards such as hyperthermia, dehydration, increased chronic diseases, social isolation, and sleep disturbances in vulnerable people, particularly the elderly [81][82][83][84].
Thus, the conducted analyses and developed cartographic models can significantly contribute to health prevention strategies and risk management related to extreme temperatures among the elderly population residing in urban areas.The process of developing MBGI cartographic models represented a pivotal stage in the research, aiming to integrate diverse data and thereby incorporate various factors for identifying optimal locations for BGI components.These developed cartographic models hold practical significance, particularly for urban planning and making informed decisions regarding sustainable city development and improving residents' quality of life.According to the outcomes of this research, proactive measures could be utilized to protect residents during extreme climate change impacts in crisis management, as well as by enhancing air quality, reducing heat islands, and providing suitable recreational spaces that positively impact residents' health.Furthermore, with regard to increasing biodiversity, the analyses pertaining to blue-green infrastructure locations could facilitate the creation of new habitats for various animal and plant species, the preservation of ecological corridors, improvement in water quality, reduction in spatial fragmentation, and the protection of endangered species.
Considering such a multitude of aspects in pinpointing optimal BGI locations allows for the optimal allocation of financial resources and the pursuit of a sustainable urban policy concerning infrastructure and the environment.Moreover, the adaptability and applicability of the models in various contexts and urban areas underscore the universality and practical value of the conducted research.

Discussion
It was stressed that data quality is of paramount significance within the analysis and that data quality is a crucial determinant in this study.Lessons learned from the BDOT10k data revealed prevalent challenges related to data completeness and timeliness within the dataset under review.Despite updating the topographic database every two years, the dynamic nature of spatial development meant that not all components were adequately represented in the dataset.The resolution of issues pertaining to the completeness and timeliness of the topographic database data necessitates a multifaceted approach.Due to financial constraints that preclude more frequent updates of the BDOT10k database, it is prudent to consider the integration of supplementary data sources, particularly those of open access, such as OpenStreetMap.Furthermore, a more intensive utilization of satellite data and remote sensing, in conjunction with machine learning and artificial intelligence techniques, could facilitate the ongoing detection of land use change and data updating.It is also of interest to obtain data for surveys from social networks or via mobile applications and online platforms, where residents can report changes and create spatial data.This solution would allow for the provision of up-to-date and realistic data, which, once verified, would significantly increase the usability and accuracy of topographic data.
The biggest challenge in utilizing the BDOT10k data is the aspect of changes occurring between subsequent updates, leading to discrepancies between the dataset and the actual spatial conditions.As a result, inaccuracies in representing new roads, buildings, or changes in land use may arise, impacting the quality of conducted analyses and potentially leading to erroneous conclusions.It should also be noted that the updating of the topographic dataset relies heavily on reference databases; however, in cases where data are missing from these references, updates are based on the visual identification of changes using orthophotomaps.Consequently, this could lead to the incorrect identification of urban development elements by the person updating the database.Nevertheless, it is important to highlight that urban centers, particularly provincial ones, are updated more frequently than every two years.As a consequence, frequent updates lead to the discovery and rectification of errors within the dataset, thus improving the overall quality of the dataset.
Given the availability of the BDOT10k data in the vector format, there is an opportunity to merge it with additional datasets.Integration provides an opportunity to improve the quality control of individual datasets while promoting a higher degree of timeliness within the combined dataset.As part of the research methodology, it was decided to include OSM data within the BDOT10k dataset.The decentralized data collection process driven by volunteer communities through the OSM database played a significant role in this decision.Volunteers often reside in specific regions and regularly update data based on their accurate knowledge of the updated area.This implies that data updates occur with greater frequency and reflect actual, real-time changes.Consequently, for specific components, such as drainage facilities, support was derived from the OSM database, which does not necessitate an administrative mode of data update.This implies that access to up-to-date data is facilitated in urbanized areas.Certainly, as was the case with the location of the existing bus stops, the dataset used also required a visual verification and supplementation of data based on the current orthophotomap.An example of errors that should have been corrected for the purposes of this study was the lack of road continuity (Figure 9).
Considering the advantages as well as the limitations of each dataset, the decision was made to rely on data integration.On the one hand, BDOT10k, based on reference sources, was used, which allowed for the identification and reliability of technical infrastructure elements.On the other hand, the speed of OSM updates favored the identification of new elements that emerged between the BDOT10k update cycles.The integration of data from disparate sources, as exemplified by the present study, facilitated cross-comparison and cross-validation, thereby enhancing the precision of the analysis.The utilization of multiple data sources not only compensates for individual dataset gaps but also contributes to a more comprehensive and reliable understanding of the studied phenomena.
In this study, the relevant emission coefficients were identified, and atmospheric corrections were applied to obtain a plausible SUHI model.Unfortunately, one of the unavoidable elements of this study was the data acquisition time.A review of the literature indicates that the SUHI phenomenon intensifies during the nighttime hours, while satellite data in Poland are recorded during the daytime, specifically between 09:30 and 10:00 for the selected area.Nevertheless, it was determined that, despite these challenges, the data can still indicate optimal locations for blue-green infrastructure (BGI).
At the same time, the research assumed the use of demographic data, which due to being obtained from official registers required additional activities, including their ordering and aggregation.A significant challenge was the absence of a clearly defined location, which was addressed by the generalization of data pertaining to the age structure of residents to districts.Obviously, acquiring these data was extremely time-consuming, and the presented analysis allows for areas to be indicated through the extent of the residential community to the location of the BGI components.Moreover, given that the implementation of BGI-related facilities is largely overseen by city management policies, it is deemed sufficient to narrow down the location of potential investments.Highlighting the strengths of the applied methodology, it would be important to mention the use of geoinformation systems, enabling the integration and comprehensive, precise analysis of spatial data from various sources.It is also important to highlight the interdisciplinary approach employed in the design of solutions for individual blue-green infrastructure elements.This approach encompasses a number of different disciplines, including geodesy, demography, hydrology, and geography.Application-oriented, the developed models, using appropriate visualization, can be a key role as a tool for decisionmakers and planners, creating solutions that will contribute to improving the quality of life of residents, as well as strengthening individual settlement units on the consequences of climate change.It is postulated that the models developed may serve as a foundation for further analysis, taking into account additional spatial features pertinent to the location of BGI.Important elements from the point of view of the analysis could also be the soil conditions, infiltration capacity of the ground, depth of the groundwater level, and ventilation model of the city.However, the data were excluded from the analyses due to their unavailability, incompleteness, and untimeliness.The utilization of open data exchange formats within geoinformation systems is expected to significantly facilitate this task.
In the context of good practices for the implementation of open data initiatives, and especially for the promotion of local sustainability, the issues of interoperability and standardization of data exchange formats are crucial [85][86][87][88][89][90].The solution is to introduce universal standards, open-source software, and appropriate regulations that legitimize the solutions adopted [91].In the context of data exchange, currently the key format is GML, which is gradually being implemented as the primary spatial data format.Nevertheless, a considerable quantity of data remains in analogue form, necessitating its conversion to the vector format for exchange.Furthermore, facilitating exchange, integration, and verification by merging multiple spatial datasets is the concept of sharing data through web services, such as WMS, WMTS, and WFS.This technology enables wide interoperability between different systems and organizations utilizing geographic information.Independence from specific software ensures the use of these standards.Furthermore, the popularization of GIS applications available under open licenses, such as QGIS, SAGA GIS, gvSIG, as well as ETL (Extract, Transform, Load) software, which is used in data processing, allowing for data to be efficiently extracted from various sources, transformed to fit business needs, and loaded into target systems, such as databases or data warehouses, is also related to this.The realization of these goals will not be possible without the active cooperation of government entities and commercial entities.In order to achieve this, it is necessary to create appropriate regulations to guide the creation, use, and sharing of open spatial datasets.Furthermore, the implementation of such regulations will facilitate the attraction of financial support, which is crucial for the maintenance of infrastructure, data management, and distribution.Furthermore, public institutions should encourage the use and creation of open datasets by organizations, particularly non-profit organizations, through concerted action.
One of the most crucial and practical aspects of BGI localization is the economic costs associated with implementing blue-green infrastructure (BGI) projects, as well as potential sources of funding for these projects.However, these issues are not addressed in this article as they are a broad and separate research problem.The cost of constructing BGI depends largely on the quality and quantity of these elements and the characteristics of the space in which they are to be located.This implies that the costs may vary considerably from one instance to another.Consequently, it is crucial to highlight that the research conducted and the methodology developed are employed to identify the most suitable location and to inform subsequent activities, including their implementation.

Conclusions
This research has shown the remarkable effectiveness of geoinformation systems, ranging from data collection and integration to spatial analysis, during the process of indicating the optimal locations of selected blue-green infrastructure (BGI) elements.The presented approach, using the availability of spatial data, through the use of appropriate formats and services, based on publicly available open-source GIS solutions, allows for significant automation of the performed analyses.The method used enables the precise identification of areas that exhibit an appropriate set of spatial features, which are optimal locations for the arrangement of the analyzed BGI elements.
However, in view of the possibility of errors in the data obtained, it is essential that those responsible for selecting the optimal BGI sites undertake verification based on other data sources and appraise the suitability of the selected locations, including through onsite inspection.This approach offers additional control and enables the improvement in the quality of the geoinformation analyses carried out.
Cartographic models, including the SUHI model and DM, were developed for this study as filters that increase the information layer and thus support decision-making regarding the location of blue-green infrastructure (BGI) elements.The interest in the surveyed aspect of spatial management, which includes improving the living conditions of residents as well as increasing the quantity and diversity of urban vegetation complexes, was also driven by the following aspects:

•
The benefits of reducing pollution, improving air quality in urban areas, and thereby improving the health and well-being of residents.In addition, locating BGI elements within the city allows for an increase in biodiversity, thus providing habitats for a variety of plant and animal species.Water conservation is also an important element in relation to the environment, as the blue-green infrastructure helps to filter rainwater.

•
Improving the aesthetics of the urban landscape by giving the urban space a unique character.Through the creation of new parks, ponds, and accompanying architectural elements, as well as elements creating urban vegetation, the aesthetics of streets and public spaces are improved, offering residents a place for relaxation, recreation, and social integration.

•
Elements of blue-green infrastructure allow for resilience to be built against the effects of climate change, in particular, by providing natural and effective thermal insulation as well as rainwater retention.Green spaces within the city allow for absorption in the event of precipitation, reducing the risk of flooding, as well as the retention of water in the ground in the event of extremely high temperatures.
In future research, it would be advisable to consider the use of alternative sets of spatial data to enhance the quality and accuracy of the data [20,46,91,92].The integration of diverse spatial datasets and advanced data acquisition technologies significantly strengthens the reliability of spatial analyses.Data from public government sources that provide statistics on land use and regulatory-driven changes are crucial for informed and effective spatial planning.In addition, observations of actual user behavior on social media platforms make it possible to assess their reactions to changes and to better understand the direction and dynamics of community action in urban environments.Consequently, the advancement of remote sensing technologies, including drones (UAVs) and LiDAR systems, enables the acquisition of high-resolution data, which are essential for the accurate and ongoing analysis of physical changes in space [93].Satellite data, with its ability to cover large areas, play a pivotal role in environmental and urban monitoring, supporting large-scale strategic planning and providing valuable data for historical analysis [94].The integration of these data sources and the use of modern technology enables planners and analysts to develop comprehensive, multidimensional models.The models presented support effective data-driven decision-making in urban development and other spatial analyses [95].With such tools, it is not only possible to plan accurately but also to adapt quickly to changing conditions and societal needs, resulting in increased quality of life for residents and sustainable spatial development.Furthermore, machine learning techniques and artificial intelligence algorithms could be employed to automate the process of data updates and validation, thereby further enhancing the reliability of the analysis.
The selection of locations for BGI elements should be complemented by taking into account the opinions and needs of the community during the design process of these spaces.Such activities can influence the successful implementation and long-term maintenance of BGI elements.Socially desirable facilities and a sense of collective action positively influence the acceptance of such concepts.Therefore, proposals for the selection of these spaces for such investments should be based on the developed map of the location potential of BGI elements.This should be followed by a field visit and public consultations to justify and strengthen the design and implementation activities.
In accordance with the aforementioned, this study presented here constitutes the initial stage of an analysis aimed at developing a matrix of optimal locational features for the deployment of blue-green infrastructure elements.This research will facilitate a more precise assessment and selection of optimal locations, taking into account a range of social, spatial, and environmental factors.Space management in this aspect is a key element of adaptation strategies to the challenges of today, allowing for the minimization of both the negative effects of urbanization processes as well as climate change.Refining and expanding this methodology is essential for the successful design and execution of blue-green infrastructure.The process involves technical aspects of urban planning and a holistic approach to creating spaces that serve environmental and social purposes.By prioritizing the development of these methodologies, urban planners and developers can ensure that blue-green infrastructure is not only effective in managing natural resources, such as water and green spaces, but also in creating environments that are inherently more welcoming and harmonious for their inhabitants.

Figure 1 .
Figure 1.Diagram for the optimization of a blue-green infrastructure (BGI) location.

Figure 2 .
Figure 2. Data used for identifying areas unsuitable for BGI based on the topographic database BDOT10k.Utilizing a geoinformation application enhances the precision of analyzing the phenomena and processes taking place in space.To ensure the integrity of the data, it was decided to use spatial datasets maintained by the National Land Surveying and Cartographic Resource.The data and their sources used in the analyses are listed below: • Hydrological corrections from SCALGO LIVE [69] (date acquired: 10 March 2023); • Digital Terrain Model (DTM) from GUGiK [68]-the Head Office of Geodesy and Cartography (date acquired: 10 March 2023); • Topographic Reference Database (BDOT10k) surface water from GUGiK [68]-the Head Office of Geodesy and Cartography (date acquired: 15 October 2020); • Topographic Reference Database (BDOT10k) from GUGiK [68]-the Head Office of Geodesy and Cartography (date acquired: 10 March 2023); • OpenStreetMap (OSM) (date acquired: 10 March 2023).

Figure 3 .
Figure 3.A map model (MBGI) showing the optimal location of the bioswale for the city of Olsztyn.

Figure 4 .
Figure 4.A map model (MBGI) showing the optimal location of the infiltration trench for the Olsztyn city.

Figure 5 .
Figure 5.A map model (MBGI) showing the optimal location of the bus stops for the Olsztyn city.

Figure 7 .
Figure 7. DM model for Olsztyn city: population in residential communities.

Figure 8 .
Figure 8.The results of the analysis conducted to determine the optimal locations for the selected BGI features include.

Figure 9 .
Figure 9.A part of the Olsztyn area near Lake Track-a visualization of the errors in the acquired data contained in the BDOT10k database.