Abstract

The urban blue-green space landscape ecology is an important component that supports the urban landscape ecosystem. In this study, a combination model of the grey system theory (GST), analytic hierarchy process (AHP), and technique for order of preference by similarity to ideal solution (TOPSIS) was used to construct criteria layers from the following four aspects: green space vegetation coverage, ecosystem diversity, air and water quality, social participation, and health. The GST method was employed to screen evaluation indicators, the AHP method was used to calculate the weight of each indicator, and finally, TOPSIS was applied for weighted ranking to verify the scientific and rationality of the evaluation system. The results indicated the following: (1) out of 36 indicator factors, 10 were not selected through the GST method; (2) the weight ranking of the criteria layer was as follows, from highest to lowest: ecosystem diversity 0.4816 >  green space vegetation coverage 0.2608 >  air and water quality 0.1864 >  social participation and health 0.0713; and (3) the relative closeness and ranking of different regions in Jinan city were as follows: Da Ming Hu area > Shan Da Lu area > Gu Cheng area > Shang Bu area. This study can help avoid decision-making errors caused by subjective factors and provide new approaches for the assessment of the ecological health of urban blue-green space landscapes.

1. Introduction

Landscape ecological health (LEH) is an interdisciplinary field that combines landscape ecology and ecosystem health studies. LEH primarily focuses on the issues of landscape ecosystem health that have suffered severe pollution and degradation under the intense human activities [1]. Urban blue-green spaces are composite ecosystems composed of water bodies and green areas [2], shaping the overall landscape structure of urban areas. Simultaneously, they are among the most severely impacted ecosystems by human societal and economic activities [3, 4]. Therefore, the assessment of urban LEH quality holds significant importance.

Urban blue-green spaces are based on blue-green corridors and patches, composed of urban green areas, water bodies, and more, forming an urban ecological spatial network that constitutes the overall ecological framework of the city. It is a complex ecological entity [5]. In China, urban green spaces are often referred to as “urban green areas” and encompass forests, urban green spaces, parks, scenic areas, residential area green spaces, and more [6]. They play a critical role as essential components of urban ecosystems and green infrastructure, fulfilling important ecosystem services (ESC) functions. Urban blue spaces typically refer to water bodies in urban environments, whether natural or artificial, including rivers, lakes, wetlands, waterfalls, and fountains. [79].

In recent years, scholars have increasingly explored urban blue-green space landscapes, with research focuses varying between domestic and international contexts. Internationally, the primary emphasis is on the study of the impact of urban blue-green spaces on residents’ health, socioeconomic factors, environmental justice, and more [7, 10, 11]. In contrast, domestic research primarily centers on the positive effects of urban blue-green spaces in areas such as society, economy, culture, riverbanks, green areas, as well as their effects on urban grids, urban heat islands, sponge cities, spatial scales, river wetlands, and their role as natural barriers [8, 9, 1217]. Currently, there is limited research on the landscape ecological health (LEH) of urban blue-green spaces, and scholars have proposed various indicators such as “landscape accessibility,” “green three-dimensional quantity,” and “even distribution of various ecological facilities” [1820]. These indicators primarily evaluate the effects of green spaces on cities, and the study of blue spaces lags behind that of green spaces, with relatively limited metrics for assessing blue spaces. Therefore, the key to assessing urban blue-green space LEH lies in the scientific selection of evaluation methods. These methods mainly include the analytic hierarchy process (AHP), fuzzy comprehensive evaluation (FCED), grey statistical theory (GST), principal component analysis (PCA), cluster analysis (CA), and TOPSIS-weighted ranking [21]. While these evaluation methods can effectively assess the quality of targets, they often lack a comprehensive understanding of the nature and mechanisms of urban blue-green space LEH. In addition, the selected evaluation indicator factors are often too simplistic, and the ability to select among these factors is lacking, resulting in an evaluation system that lacks accuracy and scientific rigor. Urban spatial layout planning, urban heat environment patterns, and the quality of living environments and physical health all depend on the construction of urban blue-green spaces. It is imperative to establish a scientific, rational, and accurate evaluation system for urban blue-green space landscape LEH under the concept of sustainable development. The objectives of this study are as follows: (1) quantifying the proportions of LEH indicator factors in urban blue-green spaces, (2) analyzing landscape indices using a combination model of GST-AHP-TOPSIS, (3) optimizing blue-green space LEH based on the best evaluation indicator factors, and (4) developing strategies for creating blue-green space landscape LEH under the concept of green development. The results of this research will provide a scientifically sound theoretical guide for urban blue-green space landscape planning, quality healthcare services, and ecosystem development. Figure 1 illustrates the framework of the article.

2. Research Object and Research Methods

2.1. Overview of the Research Area

Jinan is located in the central-western part of Shandong Province, with the Yellow River to the north and Mount Tai to the south. It is situated between approximately 36°01′ to 37°32′ north latitude and 116°13′ to 117°44′ east longitude [22]. The city falls within a warm temperate climate zone, experiencing distinct four seasons and abundant sunshine, with an average annual temperature of 13.6°C and an average annual precipitation of 614 mm [23]. The research area in this study corresponds to the administrative boundaries of Jinan city after the 2019 administrative division adjustment, and its total area has increased to 10,244 square kilometers. In 2022, Jinan had a permanent population of 9.415 million, with a per capita GDP of 128,829 yuan/person and a total regional GDP of 1,202.746 billion yuan [24]. The specific area of focus in this research is a typical central district within the old city of Jinan, including the Shangbu area, the Gucheng area, the Shanda Road area, and the Daming Lake area, constituting the basic control planning units. The total area of this study region is 30.33 square kilometers, extending east to Erhuan East Road, west to Weier Road, south to Jingshi Road, and north to Beiyuan Overpass. This area is considered the urban core of Jinan and encompasses 14 communities within a 15-minute living circle [25]. The choice of this area for the study is based on several considerations: The region is part of the old city, with relatively well-developed urban infrastructure and a substantial amount of blue and green spaces, making it suitable for the analysis of sanitary and ecological elements. It contains significant blue and green space forms such as Daming Lake Park, Baotu Spring Park, Wulongtan Park, and Qianfoshan scenic area. In addition, it boasts well-established green infrastructure and a variety of spatial types. There is existing quantitative research data available for this typical area, making it a good foundation for this study. Figure 2 illustrates the classification of blue-green spaces in Jinan City.

2.2. Research Methods
2.2.1. Grey System Theory (GST)

Grey system theory (GST) is a statistical method that applies whitening functions for mathematical operations and statistics, suitable for dealing with model frameworks involving a significant number of unknown values. Whitening function analysis generates whitening statistics for given values, which describe the level of certainty about the research object. When the greyness is 1, the set’s whiteness is 0, indicating a lack of knowledge about the research object. When the greyness is 0, the whiteness is 1, signifying a complete determination of the object under study. Traditional mathematics discusses sets with a greyness of 0 and whiteness of 1, indicating a known object. Grey sets describe populations containing both known and unknown quantities, meaning a portion is known and a portion is unknown. When evaluating the landscape quality (LEH) of urban blue-green landscapes, there are numerous indicators to consider. To obtain scientifically and reasonably selected indicators, it is necessary to use the grey statistical theory (GST) to select a set of relatively high-importance indicators, thereby obtaining a more comprehensive and authoritative set of indicators. Figure 2 illustrates the classification of blue-green spaces in Jinan city.

2.2.2. Analytic Hierarchy Process (AHP)

The analytic hierarchy process (AHP) is a multicriteria decision-making method developed by American operations researcher Thomas L. Saaty at the University of Pittsburgh [26] aimed at assisting decision-makers in making rational decisions when faced with complex choices. This method breaks down complex problems into several relatively simple hierarchies. It quantifies subjective qualitative analyses through pairwise comparisons, quantification, and ranking, creating a multilayered analytical model. Compared to other methods, AHP provides a clear mathematical logic and a well-defined hierarchy. It combines decision-makers’ subjective preferences with mathematical logic, making problem-solving systematic, concise, and not limited by mathematical background. AHP effectively transforms multicriteria problems into multilevel single-objective problems. During this process, AHP uses mathematical calculations based on the quantity relationships between elements at the same hierarchy level, avoiding the arbitrariness of a single subjective evaluation [27]. The basic steps of the AHP method include constructing judgment matrices, calculating the maximum eigenvalue and eigenvector, conducting consistency checks, and ranking hierarchical indicators. Due to the numerous factors influencing landscape ecological health (LEH) in urban blue-green spaces and the presence of fuzziness and uncertainty, the AHP method can accurately and effectively determine indicators at various levels, including criteria, objectives, and attributes, leading to a comprehensive analysis of the landscape quality indicators for urban blue-green spaces LEH.

2.2.3. Technique for Order of Preference by Similarity to Ideal Solution (TOPSIS)

TOPSIS is a ranking method initially proposed by Hwang and Yoon to approximate the ideal solution. The basic steps of the TOPSIS method are as follows: establishing an initial decision matrix, normalizing the initial matrix, determining the positive and negative ideal solutions, calculating distances, and ranking. By comparing the landscape’s positive ideal solution and negative ideal solution for various areas in Jinan, if a particular area’s total score is close to the positive ideal solution and far from the negative ideal solution, it is considered the best solution, indicating the highest quality in terms of landscape ecological health (LEH) for blue-green spaces.

3. Construction of the Landscape Ecological Health (LEH) Evaluation Index System for Urban Blue-Green Spaces

3.1. Selection of LEH Evaluation Indicators Using the Grey System Theory (GST)

The grey system theory (GST) [28] can be applied to the abstraction of LEH evaluation indicators for urban blue-green spaces, considering their multicriteria characteristics. This method involves comparing, quantifying, and extracting important indicator factors. The specific steps are as follows.

3.1.1. Preliminary Selection of LEH Evaluation Indicators for Urban Blue-Green Spaces

Through a review of the literature [2938] and consultation with experts, four criterion layers were identified, including visual psychology, ecological function, aesthetic function, and service effect, encompassing a total of 36 indicators as the preliminary selection of LEH evaluation indicators for urban blue-green spaces.

3.1.2. Calculation of Grey Whitening Functions
Questionnaire SurveyA questionnaire survey was distributed to 30 experts in the fields of landscape, geography, and the environment, using the Likert scale method with seven levels to rate the importance of the selected indicator elements. The scale ranged from 1 (very unimportant) to 7 (very important), with 2–6 indicating varying levels of importance. This process generated original data on the importance levels of the 36 selected indicators.Data ProcessingGrey whitening functions were used to process the original data concerning the importance levels of the 36 selected indicators, as obtained from the expert survey. In grey system theory, whitening functions are a mathematical method for dealing with incomplete and uncertain information. The preliminary LEH indicator set for urban blue-green spaces was categorized into three groups (high, medium, and low) based on grey classes. In accordance with the definition of grey correlation functions, high-, medium-, and low-level grey whitening functions were defined for [39] based on the grey system methods.When , the whitening function calculation formula for “high degree of importance” is as follows:When , the whitening function calculation formula for “medium degree of importance” is as follows:When , the whitening function calculation formula for “low degree of importance” is as follows:where a represents the degree of importance, b represents the index number, , represents the whitening function value of the b-th index with an importance level of a, and represents the assigned value corresponding to the importance level of the b-th index as a. Using formulas (1)–(3), the whitening function values for each index can be calculated according to the following three levels of importance: high, medium, and low.
3.1.3. Calculation of Grey Decision Coefficients and Grey Decision Vectors

Grey decision coefficients are used to measure the impact of various factors on decision outcomes. The definition of grey decision coefficient for the b-th indicator in the k-th grey class is as follows: represents the number of experts who assigns an importance value a to the b-th indicator and represents the whitening function value for the b-th indicator when its importance is a. Therefore, is defined as follows:

For each evaluation indicator, the corresponding three-level grey decision coefficients (high, medium, and low) can be calculated using equation (4) as follows:

These three coefficients together form a grey decision vector as follows:

By comparing the grey decision vectors of various evaluation indicators, important indicators with a high level of importance can be selected, completing the process of indicator selection.

3.2. Construction of LEH Evaluation Indicators for Blue-Green Spaces Using the Analytic Hierarchy Process (AHP)

The evaluation indicators selected through the grey system theory (GST) are ranked using the analytic hierarchy process (AHP). The specific steps are as follows.

3.2.1. Construction of the Judgment Matrix

Usually, in the minds of different decision-makers, the importance of different elements varies. Therefore, the analytic hierarchy process often uses pairwise comparisons of decision factors to create pairwise comparison matrices. Commonly, numerical scales are used to quantitatively process the judgments of decision-makers, resulting in a judgment matrix. The scale meanings are shown in Table 1, and these values are determined based on the intuition and judgment of decision-makers during qualitative analysis.

3.2.2. Eigenvalue Method to Obtain Maximum Eigenvalue and Eigenvector (Weight Vector)
Calculate the eigenvalue and corresponding eigenvector for each judgment matrix:Normalize the eigenvectors as follows:Calculate the maximum eigenvalue as follows:
3.2.3. Consistency Test of the Judgment Matrix
Calculate the consistency index (CI) as follows:Determine the corresponding average random consistency index (RI) from the table.The values of the average random consistency index (RI) are shown in Table 2.Based on the order of the judgment matrix, the corresponding average random consistency index RI can be obtained from Table 2.Calculate the consistency ratio CR and make a judgment as follows:

When , it is considered that the judgment matrix is completely consistent; when , it is considered that the consistency of the judgment matrix is acceptable; when , it is considered that the judgment matrix does not meet the consistency requirement, and it needs to be reconstructed and revised until it meets the consistency requirement.

3.2.4. Hierarchy Index Ranking
Individual ranking of hierarchy indices: the ranking of hierarchy indices is determined by the calculated weight vectors.Overall ranking of hierarchy indices: the overall ranking of hierarchy indices is calculated by combining the weight vectors.
3.3. Ranking of LEH Evaluation for Blue-Green Spaces in Different Areas Using the Technique for Order of Preference by Similarity to Ideal Solution (TOPSIS)

The TOPSIS method is a multiobjective decision-making technique that involves ranking based on the distance between objects and the ideal solution [40]. When combined with the AHP method, it allows for the scientific ranking of various evaluation objects through the construction of a weighted matrix. The specific steps are as follows.

Suppose there are m evaluation objects, each with n evaluation indicators.(1)Construct the judgment matrix as follows:where represents the evaluation value of the i-th evaluation object for the j-th indicator. It is specified that . The evaluation values of the i-th evaluation object for various indicators are denoted as follows:The evaluation values of all evaluation objects for various evaluation indicators are represented as the weight vector of that evaluation indicator as follows:Satisfying .(2)Standardize the judgment matrix as follows:where .(3)Construct a weighted normalized judgment matrix as follows:where (4)Determine the best vector and worst vector.Define the positive ideal solution and the negative ideal solution using the maximum and minimum values of each indicator, respectively, as follows:So, is the best vector, and is the worst vector.(5)Calculate the Euclidean distance between each evaluation object and the best vector, resulting in , and between each evaluation object and the worst vector, resulting in :where is the i-th row of the weighted normalized decision matrix .(6)Calculate the relative closeness of each evaluation object to the best vector, denoted as , as follows:If , then ; if , then ; with the condition . Therefore, the closer is to 1, the closer the evaluation object is to the best vector .(7)Rank the evaluation objects: arrange the evaluation objects in descending order based on the calculated values of relative closeness , yielding the ranking results of the evaluation objects.

4. Calculation Result Analysis

4.1. Analysis of GST Indicator Selection Results

By comparing the grey decision vectors of various evaluation indicators, important indicators with a high level of importance are selected, completing the process of indicator selection. The importance levels of landscape ecological health (LEH) evaluation indicators for urban blue-green spaces, as determined through grey statistical analysis, are presented in Table 3.

4.1.1. Analysis of Green Space Vegetation Coverage

High forest coverage, green space coverage, urban greening rate, and the distribution of green space areas effectively mitigate the impact of urban climate change, improve the quality of life for urban residents, and promote urban sustainability. The temperature regulation capacity, air humidity regulation, and CO2 absorption capacity help alleviate the urban heat island effect and contribute to improving the city’s climate, residents’ quality of life, and urban sustainability.

Three indicators of green space vegetation coverage were not selected. While floral diversity, garden coverage, and lawn coverage contribute to beautifying the urban environment and have a positive impact on the psychological well-being of urban residents, their ecological service functions are limited to a small local area and do not have a significant ecological impact. The evaluation index system for LEH green space vegetation coverage is analyzed as shown in Figure 3.

4.1.2. Analysis of Ecosystem Diversity

The level of wetland conservation helps protect urban wetland ecosystems. The degree of biodiversity and species richness is instrumental in establishing robust ecosystems. The integrity of natural landscapes contributes to urban greening, air quality improvement, and environmental beautification. Maintaining water quality health supports safe drinking water supply, aquatic ecosystems, and recreational opportunities. Ecosystem services can provide various services to the city, promoting the physical and mental well-being of residents and the city’s sustainability.

Two indicators of ecosystem diversity were not selected. The conservation of endangered species and habitat area, due to factors such as data availability, priority, methodology, resource constraints, require accurate data for comprehensive assessment and were, therefore, not included. The analysis of the LEH green and blue space ecosystem diversity evaluation index system is shown in Figure 4.

4.1.3. Air Quality and Water Quality

Air Quality Index, Water Quality Index, pollutant adsorption, wastewater treatment efficiency, water accessibility, and river and lake health directly impact the environmental quality and quality of life for city residents. They also reflect the overall environmental health of the city. Good air and water quality have a significant impact on the lives and health of urban residents.

Two indicators for air quality and water quality were not selected. While waste management efficiency and flood risk are crucial urban environmental issues, their selection was influenced by factors such as prioritization, resources, and policy directions. They are not as urgent as other indicators, which is why they were not included. The analysis of the LEH green and blue space air quality and water quality evaluation index system is shown in Figure 5.

4.1.4. Social Engagement and Health

The Community Health Index is a key factor reflecting the quality of life for community residents. Facility maintenance, timeliness of services, and ecological corridor connectivity contribute to providing city residents with safer, more reliable, and timely services, promoting the ecological sustainability of the city. Improving psychological health, public engagement, and environmental protection provide psychological health services, democratic decision-making, protect various city resources, and enhance the quality of the urban environment.

Three indicators for social engagement and health were not selected. Although environmental education opportunities, public recreational spaces, and information disclosure contribute to raising awareness of environmental issues and providing recreational spaces, these indicators are not highly correlated with urban blue-green spaces and were, therefore, not included. The analysis of the LEH green and blue space social participation and health evaluation index system is shown in Figure 6.

4.2. Analysis of AHP Indicator Weight Results

Based on the evaluation index system, we conducted a questionnaire survey from June 14th to 22nd, 2023, inviting 30 experts from universities and companies such as Shandong University, Jinan University, Shandong University of Art and Design, Shandong Jianzhu University, and Shandong Tongyuan Design Co., Ltd. All experts reside within the research area of this paper. The survey questionnaire consisted of background questions and main questions. The background questions included respondents’ personal basic information and their usage of the research area, such as gender, age, education level, and recreational activities in the area. The main questions asked respondents to quantitatively rate the landscape ecological health of the research area based on their subjective perceptions. The questionnaire was self-administered and consisted of four pages of A4 paper, including four background questions and ratings for 26 evaluation index factors. It could be completed within 5 minutes, and respondents were instructed to select appropriate options and fill in appropriate scores based on their own circumstances.

A total of 30 questionnaires were distributed, all of which were returned, resulting in 30 valid questionnaires, with an effective rate of 100%. The respondents covered the following three age groups: youth (18–40 years old), middle aged (40–65 years old), and elderly (over 65 years old), accounting for 36.16%, 50.44%, and 13.04%, respectively. There were 22 males and 8 females; educational backgrounds ranged from undergraduate to doctoral degrees, with a majority holding master’s degrees, accounting for 27.21%, 40.07%, and 20.12%, respectively. Regarding the frequency of recreational activities in the area, 1 person visited for the first time, accounting for approximately 1.8% of the total respondents, 5 people visited for the second time, accounting for about 19.6%, and the remaining 24 people visited three times or more, accounting for approximately 78.6%. The 30 experts conducted pairwise comparisons of the 26 different evaluation indicators based on their relative importance and determined the results of the 30 judgments through weighted averaging, upon which five judgment matrices were constructed AD, , , , and . Using the eigenvalue method to calculate the eigenvectors, which represent the weight values of urban blue-green space LEH evaluation indicators (Table 4), and conducting a consistency check, the calculation results indicate that the CR values for the judgment matrices AD, , , , and are 0.0518, 0.0448, 0.0548, 0.0329, and 0.0261, respectively. All CR values are less than 0.1, passing the consistency test.

4.2.1. Analysis of Evaluation Indicator Weight Calculation Results

Through the analytic hierarchy process, a comprehensive evaluation analysis of the LEH effectiveness of hygienic city blue-green spaces was conducted, and the weight ranking results for the criterion layer are as follows: ecosystem diversity 0.4816 >  green area vegetation coverage 0.2608 >  air and water quality 0.1864 >  social engagement and health 0.0713. As hygienic city blue-green spaces primarily aim to provide a healthy and livable environment for urban residents, they receive the highest weighting. Green area vegetation coverage in urban blue-green spaces has a higher weight than air and water quality, and its impact is also greater than ecosystem diversity, so it receives a higher weight than social engagement and health.

In the total weight ranking of the 26 indicator levels, the following conclusions were drawn: ecosystem services () has the highest weight, followed by natural landscape integrity (). These two indicators have a significant impact on the comprehensive evaluation of hygienic city blue-green space LEH. The data also show that forest coverage (), diversity index (), and air quality index () are similarly crucial for hygienic city blue-green space LEH. The weights of other indicators are as follows, in descending order: green area coverage rate () > species richness () > Water Quality Index () > temperature regulation capability () > wetland conservation level () > air humidity regulation () > river and lake health () > Community Health Index () > wastewater treatment efficiency () > CO2 absorption capacity () > water quality health () > psychological health improvement () > distribution of green areas () > pollutant adsorption () > public participation () > urban greening rate () > water accessibility () > environmental protection () > timeliness of services () > ecological corridor connectivity () > facility maintenance and upkeep ().

4.3. TOPSIS Indicator Weight Analysis
4.3.1. Questionnaire Survey

Satisfaction surveys were conducted based on 36 indicators within four major categories that revolve around factors influencing the usage of hygienic city blue-green spaces, as shown in Table 5. The Likert five-point scale was employed, with a total of 25 questions. Respondents were required to provide a satisfaction rating for each question, ranging from “very dissatisfied,” “dissatisfied,” “neutral,” “satisfied,” to “very satisfied,” which corresponded to scores of 1, 2, 3, 4, and 5, respectively.

4.3.2. Construction of the Judgment Matrix

The evaluation of the use of urban blue-green space LEH includes four segments to be evaluated, and each segment has 26 evaluation indicators. Therefore, the judgment matrix is as follows:

Since all 26 indicators studied in this paper are positive indicators, there is no need for further normalization.

4.3.3. Standardized Judgment Matrix

Using the range standardization method, the original data were subjected to linear transformation to obtain the standardized judgment matrix, which represents the normalized data for the Shangbu District, Gucheng District, Shanda Road District, and Daming Lake District in Jinan, as shown in Table 6.

4.3.4. Construction of the Weighted Normalized Judgment Matrix

The weight values of each factor in the comprehensive evaluation indicator system for the LEH effect of urban blue-green spaces obtained through the AHP method are multiplied by the corresponding values in the standardized judgment matrix to construct the weighted normalized judgment matrix. This matrix represents the weighted values for Shangbu District, Gucheng District, Shanda Road District, and Daming Lake District in Jinan, as shown in Table 7.

4.3.5. Calculating the Relative Closeness of Each District

Based on the magnitude of the relative closeness (), Jinan’s Shangbu District, Gucheng District, Shanda Road District, and Daming Lake District are ranked to provide decision-making criteria. Using the calculated weighted values for each district, we first construct the best vector ( positive ideal solution) and the worst vector ( negative ideal solution). Then, we calculate the Euclidean distances between each district and the positive and negative ideal solutions, denoted as and , respectively. Using the formula for relative closeness , a higher value indicates a higher level of satisfaction for that district. Therefore, arranging the relative closeness in the descending order provides the satisfaction ranking results for Shangbu District, Gucheng District, Shanda Road District, and Daming Lake District, as shown in Table 8.

5. Urban Blue-Green Space LEH Landscape Quality Optimization Strategies

5.1. Building a “Blue-Green Integration, City-Water Unity” Blue-Green Space Layout

Under the LEH concept, by implementing “interconnected river systems, integrated green corridor landscapes, and complementary park and square green spaces” [41], it is possible to create a convergence of urban water networks, urban ecological networks, urban landscape networks, and urban public space networks. This fosters a livable, healthy, and sustainable urban environment, helping to maintain urban ecological balance, alleviate the urban heat island effect, improve the quality of life for urban residents, and enhance resource efficiency [42]. Through high-performance core nodes, the optimization of blue-green spaces is expanded from point to area, transforming Jinan into a three-level blue-green space system: “regional blue-green space, community blue-green space, and node blue-green space.” This system is anchored by core optimization frameworks, such as Daming Lake, Baotu Spring, and Qianfoshan, and integrates urban blue-green space networks within the city and edge ecological green belts, creating an “urban space structure of mountains, springs, lakes, rivers, and the city” [43].

5.2. Constructing an Ecological Grid Structure to Enhance Urban Space Quality

Focusing on the integration of urban attributes related to “nature-society-economy-culture,” rational planning of blue-green spaces within the city’s boundaries is undertaken to create an ecological grid system with continuous connectivity and integrity, including mountains, rivers, lakes, and green spaces [44, 45]. To better harness the coupled security effects of blue-green spaces, the main measures include (1) connecting rivers and lakes through new construction, identifying low-power water systems, strengthening linear water body planning, restructuring the urban water grid structure, and preventing water disasters. (2) Building a blue-green space grid in the urban main and subsidiary rivers, large urban green belts, and urban main roads, among other natural and traffic buffer zones, to truly mitigate urban pollution effects and enhance connectivity of urban ecological sectors [46]. (3) Linking natural ecological resources such as mountains, water bodies, wetlands, parks, and green spaces with human-made resources such as farmland and protective forests to interconnect and improve the ecological network, forming an interconnected and urban ecosystem.

5.3. Constructing Human-Centric Blue-Green Spaces to Enhance Public Service Capabilities

The key to constructing the human-made environment in urban blue-green spaces lies in establishing an urban park and green space system that seamlessly integrates with the city. One of the goals of the ecological health of urban blue-green spaces is to create core cultural elements of urban greenery and water. To enhance the service capabilities of urban blue-green spaces, it is necessary to increase service level indicators, including per capita blue space, accessibility, and the water view rate of surrounding residents, to supplement traditional park and green space systems [47, 48]. This involves unearthing historical and cultural elements, including material spatial culture and nonmaterial spiritual cultural activities, integrating point-like cultural spaces within urban blue-green spaces to meet various production and service activities in the city and enhance the city’s economic and cultural level [49].

6. Discussion and Conclusion

6.1. Discussion

Integrating the grey system theory (GST), analytic hierarchy process (AHP), and technique for order of preference by similarity to ideal solution (TOPSIS) in assessing the landscape ecological health (LEH) of urban blue-green spaces offers a comprehensive framework for evaluating the intricate relationships among ecological, social, and economic factors. Our analysis revealed several key points as follows: (1) methodological approach: adopting the GST-AHP-TOPSIS combination model enabled systematic evaluation of diverse criteria layers, facilitating the identification of influential factors impacting urban blue-green space LEH. This methodological approach ensured a comprehensive understanding of landscape health and offered a structured framework for decision-making. (2) Indicator selection and weighting: the selection and weighting of evaluation indicators were critical steps in our assessment process. Through expert consultations and surveys, we identified and prioritized key indicators such as ecosystem diversity, green space vegetation coverage, air and water quality, and social participation and health. The AHP method facilitated the calculation of indicator weights, offering insights into the relative importance of each criterion layer. (3) Regional disparities: while areas like Daming Lake District demonstrated higher landscape health levels, others, such as the Shangbu District, encountered more significant challenges. These results highlight the necessity for targeted interventions and policy actions to tackle specific environmental issues and enhance overall landscape quality. (4) Sensitivity analysis: in order to ensure the reliability of our conclusions, we performed a sensitivity analysis to gauge the impact of potential variations in input parameters on our results. This analysis involved altering key parameters, such as the weights assigned to different evaluation indicators and the criteria used for district ranking. Remarkably, the results of the sensitivity analysis indicated that despite fluctuations in these parameters, the overall district ranking remained relatively consistent. This underscores the robustness of our conclusions regarding the relative performance of districts in terms of urban blue-green space LEH, as they are not substantially affected by variations in input parameters. Moving forward, we aim to conduct more extensive sensitivity analyses to investigate the influence of additional factors on our findings. By doing so, we will further enhance the credibility and applicability of our evaluation framework.

6.2. Conclusion

In conclusion, our study highlights the significance of adopting an integrated approach to assess the landscape ecological health (LEH) of urban blue-green spaces. By leveraging the GST-AHP-TOPSIS combination model, we have developed a systematic and scientifically grounded evaluation framework that can inform decision-making processes and guide urban planning initiatives. Our findings emphasize the need for comprehensive strategies aimed at enhancing the quality and resilience of urban landscapes, thereby promoting the well-being of residents and fostering sustainable urban development. Moving forward, continued research in this area is essential to address emerging environmental challenges and ensure the long-term health and vitality of urban blue-green spaces.

This study has several limitations, including data limitations, methodological limitations, regional limitations, and temporal limitations. First, regarding data, the data used in this study primarily came from literature and expert surveys, which may suffer from issues of incompleteness or lack of comprehensiveness, potentially affecting the accuracy and reliability of the research results. Second, concerning methodology, although this study utilized the GST-AHP-TOPSIS combination model for evaluation, different evaluation methods may yield varying results, introducing a certain degree of subjectivity and bias. In addition, in terms of region, this study primarily focused on evaluating the blue-green spaces of specific cities, and the actual situations in different regions may exhibit significant differences, thereby limiting the generalizability and applicability of the research findings due to regional constraints. Lastly, regarding time, the data and information involved in this study are concentrated within specific time periods, and future changes may impact the stability of the research results. Regular updates of data and information are necessary to ensure the timeliness and reliability of the research findings. In future research, efforts can be made to diversify and broaden the sources of data, explore and compare the effectiveness of different evaluation methods, expand the scope of research, and engage in long-term tracking and monitoring of the development and changes in blue-green spaces. These endeavors aim to better understand and address the various limitations inherent in the evaluation of urban blue-green space health, thus advancing the continuous development and progress of this field.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors wish to thank all participants of this study for their selfless participation and dedication. The authors also acknowledge National Social Science Fund Art Project, 2018 (project no.: 18BA009).