Comprehensive Assessment of Water Sensitive Urban Design Practices based on Multi-criteria Decision Analysis via a Case Study of the University of Melbourne, Australia

Water sensitive urban design (WSUD), as a typical green stormwater infrastructure (GSI), contains various facilities to decrease the urbanization impacts and enhance the values of amenity, ecosystem, and livability in Australia. Although WSUD has developed over 30 years, existing studies for WSUD performances have sometimes ignored its economic and social benefits, and there is still a lack of an integrated framework to optimize the GSI combinations based on various criteria in a site. This paper aims to utilize “score-rank-select” strategy to comprehensively assess WSUD combination scenarios from functional, economic, social, and environmental aspects, by taking the University of Melbourne (Parkville campus) as a case study. In detail, multi-criteria decision analysis (MCDA) was used for weight determination and scenario comparison. The results showed that scenario 4 with 52% green WSUD facilities had the highest assessment score (0.771) among the five scenarios, while the final score (0.758) of scenario 5 was lower than scenario 4 although its green facility proportion reached 69%. The trade-off relation between the proportion of grey and green WSUD facilities was further demonstrated. Additionally, this paper strongly recommends that the MCDA-based comprehensive assessment framework described here can be generally promoted for the water sector to solve the decision-making problems. The use of such a framework can further promote sustainable development by helping water managers to make informed and inclusive decisions involving a variety of factors.


Introduction
Urbanization in recent years has dramatically increased the amount of impervious surfaces [1] and affected local environmental conditions such as landscape or stream morphology [2], species richness [3,4], air temperature [5], climate change [6], and particularly a region's hydrology [7][8][9][10][11]. When heavy storms hit, a large volume of stormwater with pollutants containing sediment, nutrients, litter, oxygen-demanding waste, and heavy metals runs off the impervious areas, giving rise to flooding risk and water quality concerns [12][13][14], as well as a series of social problems [15]. Despite the fact that the conventional stormwater management approaches including gutters, pipes, and channels can control flooding by conveying runoff directly into receiving points [16][17][18], they have been considered unsuitable for urban sustainability development [16,19,20] because of limited storage [21],

Study Area
The study area (The University of Melbourne: Parkville campus, approximately 580,000 m 2 ) is located at the northern edge of the Elizabeth Street catchment of the City of Melbourne council (Figure 1), which has been significantly affected by flooding in the past. In history, the City of Melbourne and Victoria State Emergency Service recorded many significant flooding issues in this catchment [57,58]. Also, the study area has obvious characteristics of lowland and highly urbanized areas with a variety of land uses [57], meaning the flooding risks will further increase as the urban development process. Compared to other regions of Australia, the University of Melbourne is a representative of highly developed and urbanized area. Therefore, there is a need for WSUD implementation to decrease the urbanization impacts, while the variety of land uses provides sufficient spaces for WSUD construction in the study area. Current stormwater management in the study area is mainly based on traditional drainage pipe systems. It was considered unsuitable for urban sustainable stormwater development, though several WSUD projects have been included in the new campus buildings [59] ( Table 1). It has been found that the existing stormwater tanks are grey-based WSUD facilities, and the green-based WSUD facilities such as bioretention areas and swales are very few. For these reasons, it is far from achieving the sustainable stormwater management goals (Table 2) at the University of Melbourne, and it is necessary to design an optimal WSUD combination for the university to solve the problem of stormwater management and urbanization impacts.
Water 2020, 12, x FOR PEER REVIEW  3 of 35 following section, we briefly introduce the study area and the assessment framework and method. The results, discussions, recommendations, and limitations are shown in Section 3.

Study Area
The study area (The University of Melbourne: Parkville campus, approximately 580,000 m 2 ) is located at the northern edge of the Elizabeth Street catchment of the City of Melbourne council (Figure 1), which has been significantly affected by flooding in the past. In history, the City of Melbourne and Victoria State Emergency Service recorded many significant flooding issues in this catchment [57,58]. Also, the study area has obvious characteristics of lowland and highly urbanized areas with a variety of land uses [57], meaning the flooding risks will further increase as the urban development process. Compared to other regions of Australia, the University of Melbourne is a representative of highly developed and urbanized area. Therefore, there is a need for WSUD implementation to decrease the urbanization impacts, while the variety of land uses provides sufficient spaces for WSUD construction in the study area. Current stormwater management in the study area is mainly based on traditional drainage pipe systems. It was considered unsuitable for urban sustainable stormwater development, though several WSUD projects have been included in the new campus buildings [59] ( Table 1). It has been found that the existing stormwater tanks are grey-based WSUD facilities, and the green-based WSUD facilities such as bioretention areas and swales are very few. For these reasons, it is far from achieving the sustainable stormwater management goals (Table 2) at the University of Melbourne, and it is necessary to design an optimal WSUD combination for the university to solve the problem of stormwater management and urbanization impacts.    Table 2. Sustainable stormwater management goals in the study areas.

No. Sustainable Stormwater Goals Reference
1 Reduce the risks and damage caused by flooding.
[57, 60,62] 3 Improve the health of existing vegetation, urban forest, and ecosystem by water quality improvement.
[57, 60,61,63,64] 4 Increase the provision of open space in the catchment with canopy cover to support a healthy population.
[57, 60,61,63,64] To make better use of WSUD facility combinations in the design, the study area is divided into seven sub-areas by site characteristics. The land use of study area can be simply classified as roof and buildings, road and other impervious areas, lawn and trees, and a few landmarks (e.g., South Lawn, playground, University Square, and Lincoln Square) (See Figure 1). The impervious rates of different land use are summarized in Table 3 and the area size and characteristics of the seven sub-areas are shown in Table 4. It is noted that the land use of University Square is assumed as "urban mix area" because the impervious rate is approximately 50%. The rest of the landmarks are regarded as "lawn and trees" surface type. The rainfall and evapotranspiration data are both from 1 January 2019 to 31 December 2019 in a daily time step and shown in Figure 2. Flemington station (rainfall data) [66] and Melbourne Airport station (evapotranspiration data) [67], as the closest stations to the University of Melbourne, were selected as the data source. The concentration of pollutants of total suspended solid (TSS), total nitrogen (TN), and total phosphorus (TP) were stochastically generated by MUSIC following Gaussian Distribution. Table 5 shows the parameters of these pollutants classified by different surface types.

Five Scenarios Design
Based on different characteristics and land uses of seven sub-areas, the schematic diagram of five designed scenarios comprising WSUD facilities for the University of Melbourne are shown in Figure 3. In this figure, the dashed arrows indicate the assumed flow paths, and the solid arrow refers to the flow outlet to the municipal pipes. Four WSUD facilities including bioretention areas, pond, rainwater tank, and swale were selected for design, and the size of different WSUD facilities used are summarized in Table 6. The descriptions, functions, parameters, and specifications of different WSUD facilities for designs are provided in Appendix A.

Five Scenarios Design
Based on different characteristics and land uses of seven sub-areas, the schematic diagram of five designed scenarios comprising WSUD facilities for the University of Melbourne are shown in Figure 3. In this figure, the dashed arrows indicate the assumed flow paths, and the solid arrow refers to the flow outlet to the municipal pipes. Four WSUD facilities including bioretention areas, pond, rainwater tank, and swale were selected for design, and the size of different WSUD facilities used are summarized in Table 6. The descriptions, functions, parameters, and specifications of different WSUD facilities for designs are provided in Appendix A. Water 2020, 12, x FOR PEER REVIEW 6 of 35  To explore how the changes of proportion of grey and green WSUD facilities affect the final assessment scores in the study area, five scenarios were generated to be further evaluated. The five scenarios comprised different proportions of grey and green WSUD facilities ( Figure 4). From scenario 1 to scenario 5, the proportion of green-based WSUD facilities was gradually increasing. It is noted that the scenario 1 to 4 contained the existing 225 m 2 stormwater tanks, while scenario 5 assumed that the 750 kiloliters (150 m 2 ) stormwater tank installed in the Glyn Davis Building did not exist. The descriptions of 5 scenarios are presented as follow. To explore how the changes of proportion of grey and green WSUD facilities affect the final assessment scores in the study area, five scenarios were generated to be further evaluated. The five scenarios comprised different proportions of grey and green WSUD facilities ( Figure 4). From scenario 1 to scenario 5, the proportion of green-based WSUD facilities was gradually increasing. It is noted that the scenario 1 to 4 contained the existing 225 m 2 stormwater tanks, while scenario 5 assumed that the 750 kiloliters (150 m 2 ) stormwater tank installed in the Glyn Davis Building did not exist. The descriptions of 5 scenarios are presented as follow.
Scenario 1: "Grey-based" design. This scenario mainly relied on the large size of stormwater tanks and ponds for water harvesting and retention. Its total areas of grey-based facilities reached 675 m 2 while the green-based areas were only 30 m 2 . Scenario 2 and Scenario 3: "Grey-green" design. The two scenarios comprised all the four WSUD technologies with different sizes. Scenario 2 had a higher proportion of bioretention areas and stormwater tanks, and scenario 3 comprised more swales. Scenario 4: "Equal grey-green" design. This scenario had more plant-support facilities including bioretention areas and swales (a total of 320 m 2 ), while its size of grey-based facilities (a total of 300 m 2 ) was much smaller than first three scenarios. The proportions of grey and green facilities were equal. However, the existing 225 m 2 stormwater tanks in this scenario made it difficult to further increase the proportion of green-based facilities.
Scenario 5: "Green-based" design. It designed more bioretention areas and swales (a total of 450 m 2 ) than scenario 4 and assumed that the 750 kiloliters (150 m 2 ) stormwater tank installed in the Glyn Davis Building did not exist.

Criteria and Indicator Selection
According to general stormwater goals (See Table 2), GSI cases [68][69][70][71], and functional specifications of selected WSUD (See Appendix A), 10 sub-criteria and several indicators from functional, economic, social, and environmental variables were determined and shown in Figure 5. The rationale for each sub-criterion and indicator is as follows. Scenario 1: "Grey-based" design. This scenario mainly relied on the large size of stormwater tanks and ponds for water harvesting and retention. Its total areas of grey-based facilities reached 675 m 2 while the green-based areas were only 30 m 2 . Scenario 2 and Scenario 3: "Grey-green" design. The two scenarios comprised all the four WSUD technologies with different sizes. Scenario 2 had a higher proportion of bioretention areas and stormwater tanks, and scenario 3 comprised more swales. Scenario 4: "Equal grey-green" design. This scenario had more plant-support facilities including bioretention areas and swales (a total of 320 m 2 ), while its size of grey-based facilities (a total of 300 m 2 ) was much smaller than first three scenarios. The proportions of grey and green facilities were equal. However, the existing 225 m 2 stormwater tanks in this scenario made it difficult to further increase the proportion of green-based facilities.
Scenario 5: "Green-based" design. It designed more bioretention areas and swales (a total of 450 m 2 ) than scenario 4 and assumed that the 750 kiloliters (150 m 2 ) stormwater tank installed in the Glyn Davis Building did not exist. According to general stormwater goals (See Table 2), GSI cases [68][69][70][71], and functional specifications of selected WSUD (See Appendix A), 10 sub-criteria and several indicators from functional, economic, social, and environmental variables were determined and shown in Figure 5. The rationale for each sub-criterion and indicator is as follows. Functional aspect. The functions of WSUD mainly comprised flooding control and water quality improvement [38,39,72,73]. For stormwater flooding control, total water quantity reduction and surface runoff, as important factors affecting runoff control efficiency, were selected as two indicators based on related cases [36,42]. WSUD facilities can promote evapotranspiration to reduce the total water quantity and temporarily retain water to reduce surface runoff [13,74]. For water quality improvement, the major pollutants included total suspended solid (TSS), total nitrogen (TN), and total phosphorus (TP) [75]. The WSUD's performance of TSS, TN, and TP treatment has been demonstrated by Parker [76], Hatt et al. [38] and Passeport et al. [72]. Hence, TSS, TN, and TP were selected as three indicators of water quality improvement.
Economic aspect. Compared to traditional approaches, the application of WSUD facilities or general GSI can reduce the capital costs and operation and maintenance (O&M) costs spent on stormwater treatment [77,78]. Additionally, land use cost is also a consideration in WSUD's economic performance because the land used for WSUD could have been used for other purposes [29,47]. The more land WSUD facilities occupied, the greater the land cost incurred, while the conventional pipe systems are generally buried under the ground. Thus, capital cost, O&M cost, and land use cost were selected as economic subcriteria.
Social aspect. Livability, water reuse functions, and water retaining functions were selected as three social sub-criteria because they are not only important factors for people who live around the study area, but also are identified as significant goals and targets in stormwater management (see Table 2). The functions of water reuse and water retaining provided by WSUD have been broadly demonstrated in the guidelines [32]. Livability comprises recreational values, people's wellbeing, and aesthetic values, which are closely related to people themselves [44]. Further, Leonard et al. [44], Moores and Batstone [79], and Functional aspect. The functions of WSUD mainly comprised flooding control and water quality improvement [38,39,72,73]. For stormwater flooding control, total water quantity reduction and surface runoff, as important factors affecting runoff control efficiency, were selected as two indicators based on related cases [36,42]. WSUD facilities can promote evapotranspiration to reduce the total water quantity and temporarily retain water to reduce surface runoff [13,74]. For water quality improvement, the major pollutants included total suspended solid (TSS), total nitrogen (TN), and total phosphorus (TP) [75]. The WSUD's performance of TSS, TN, and TP treatment has been demonstrated by Parker [76], Hatt et al. [38] and Passeport et al. [72]. Hence, TSS, TN, and TP were selected as three indicators of water quality improvement.
Economic aspect. Compared to traditional approaches, the application of WSUD facilities or general GSI can reduce the capital costs and operation and maintenance (O&M) costs spent on stormwater treatment [77,78]. Additionally, land use cost is also a consideration in WSUD's economic performance because the land used for WSUD could have been used for other purposes [29,47]. The more land WSUD facilities occupied, the greater the land cost incurred, while the conventional pipe systems are generally buried under the ground. Thus, capital cost, O&M cost, and land use cost were selected as economic sub-criteria. Social aspect. Livability, water reuse functions, and water retaining functions were selected as three social sub-criteria because they are not only important factors for people who live around the study area, but also are identified as significant goals and targets in stormwater management (see Table 2). The functions of water reuse and water retaining provided by WSUD have been broadly demonstrated in the guidelines [32]. Livability comprises recreational values, people's wellbeing, and aesthetic values, which are closely related to people themselves [44]. Further, Leonard et al. [44], Moores and Batstone [79], and Bowen et al. [80] reported that WSUD or other GSI can effectively improve the human-related values. Because of the different performance of water retaining functions for grey and green facilities [30,81], the two types of facilities were assessed respectively under the sub-criteria of water retaining function.
Environmental aspect. Many cases reported that carbon emission is an essential environmental indicator in stormwater infrastructure development [40,82] and it was demonstrated that WSUD can provide less carbon emission than conventional approaches [40,41]. Also, this indicator is consistent with the target of "net-zero emission" and "climate change adaptation" in the water sector (See Table 2). Meanwhile, ecosystem value is another environmental benefit that the selected WSUD facilities can provide [83], which is helpful to improve the water's health and biodiversity of the study area.

Weight Determination
The AHP method, developed by Saaty, is mainly used to determine the weights of sub-criteria and indicators [55]. The term consistency index (CI) is used to ensure the accuracy of AHP method. If CI result is less than 0.1, the weights determined by AHP method are considered reliable and the further assessment can be conducted [84]. The initial pair-wise comparison matrix of AHP is determined by seven procedures and shown in Figure 6. The designed questionnaire template is provided in Appendix B and the questionnaire results are enclosed in the Supplementary Materials. Also, the procedure of questionnaire assessment is to select effective questionnaires. In detail, the results provided by the experts who are not very relevant to this study and the results whose CI values are less than 0.1 will not be accepted. Taking the subjectivity of these experts into account, the procedure of peer review allows more experts to assess the rationality of these values to ensure the accuracy of the assessment. In this research, four experts from China University of Geoscience filled out the questionnaires and three of them were accepted by questionnaire assessment. Two experts from the University of Melbourne participated in the peer review. Because of the different performance of water retaining functions for grey and green facilities [30,81], the two types of facilities were assessed respectively under the sub-criteria of water retaining function. Environmental aspect. Many cases reported that carbon emission is an essential environmental indicator in stormwater infrastructure development [40,82] and it was demonstrated that WSUD can provide less carbon emission than conventional approaches [40,41]. Also, this indicator is consistent with the target of "net-zero emission" and "climate change adaptation" in the water sector (See Table 2). Meanwhile, ecosystem value is another environmental benefit that the selected WSUD facilities can provide [83], which is helpful to improve the water's health and biodiversity of the study area.

Weight Determination
The AHP method, developed by Saaty, is mainly used to determine the weights of sub-criteria and indicators [55]. The term consistency index (CI) is used to ensure the accuracy of AHP method. If CI result is less than 0.1, the weights determined by AHP method are considered reliable and the further assessment can be conducted [84]. The initial pair-wise comparison matrix of AHP is determined by seven procedures and shown in Figure 6. The designed questionnaire template is provided in Appendix B and the questionnaire results are enclosed in the Supplementary Materials. Also, the procedure of questionnaire assessment is to select effective questionnaires. In detail, the results provided by the experts who are not very relevant to this study and the results whose CI values are less than 0.1 will not be accepted. Taking the subjectivity of these experts into account, the procedure of peer review allows more experts to assess the rationality of these values to ensure the accuracy of the assessment. In this research, four experts from China University of Geoscience filled out the questionnaires and three of them were accepted by questionnaire assessment. Two experts from the University of Melbourne participated in the peer review.   Weight determination procedures based on AHP method are shown in Figure 7 and the detailed calculation results are provided in Appendix C. The weights of indicators under a sub-criterion were qualitatively determined from literature (water retaining function) or equal weighting systems (flooding control and water quality improvement).
Weight determination procedures based on AHP method are shown in Figure 7 and the detailed calculation results are provided in Appendix C. The weights of indicators under a sub-criterion were qualitatively determined from literature (water retaining function) or equal weighting systems (flooding control and water quality improvement).   Table 7 shows, each sub-criteria and indicator were respectively assessed qualitatively or quantitatively. The functional benefits including flooding control and water quality improvement were evaluated by MUSIC. The rainfall-runoff model in MUSIC is based on the research of Chiew and McMahon [85]. In this model, the rainfall can go into impervious store, pervious store, and ground water. The outflow can be divided into surface runoff, baseflow, and deep loss. In this paper, the models of scenarios disregarded deep loss and only considered surface runoff and baseflow.  Table 8 summarizes the important parameters and coefficients for sub-criteria calculation. The data sources of initial capital cost and annual O&M cost of each WSUD facility were extracted from Melbourne Water [86] and Jayasuriya and Khastagir [87]. The indicators' values of capital cost and O&M cost can be easily obtained by simple cost calculations. The land use cost, water reuse, and water retaining function were estimated by the size of corresponding WSUD facilities in different scenarios. It is noted that the water retaining functions of grey-based facility (tank and pond) and green-based facility (bioretention area and swale) were calculated separately because of different performance features. As a human-related indicator, the livability of each scenario was assessed by using hedonic pricing method. It utilizes the home price fluctuation to estimate the values related to people's willingness to pay. Ira [88] provided the proportion of house price increase around each WSUD facility. The carbon footprint was used to calculate the carbon emissions. The emissions contribute to global warming, which can be defined as the total amount of greenhouse gases (GHG) released directly or indirectly into the atmosphere by human activities [89]. The raw data of the carbon footprint of different WSUD facilities were provided by Moore and Hunt's research [40] based on a 30-year life cycle. The ecosystem value was assessed by qualitative analysis. It was assigned an identical score with the water quality indicator because of the close relationships between the two sub-criteria.
As some sub-criteria comprised several indicators (see Figure 5), weighted sum model (WSM) (Equation (1)) was utilized to calculate the values of flooding control, water quality improvement, and water retaining function.
where: A i is the score of sub-criteria in scenario i, w i a ij is the performance of alternative A i when it is evaluated in terms of indicator j, n is the number of indicators under one sub-criterion, and a ij is the linear normalized results for indicator j, which is given by Equation (2) where: B ij is the initial value of indicator j in scenario i, and N refers to the number of WSUD scenarios. The indicator values comprised benefit value and non-benefit value, and the non-benefit values should be transferred into benefit forms by Equation (3).
where: b ij is the benefit form value after transformation, B max(ij) is the maximum value of indicator j in scenario i, and B min(ij) is the minimum values of indicator j in scenario i.

Comprehensive Assessment
After determining the values of sub-criteria, the MAUT method, as an extension of Multi-Attribute Value Theory (MAVT) [51], was adopted to obtain the comprehensive results. It is the closest to the public acceptance MCDA method [55]. It assumes that the preferences of decision-makers can be presented as an analytical utility function and effectively combined the risk preferences and uncertainty in decision-making [90]. Generally, MAUT's functions are divided into three main types: additive, multiplicative, and reference point [91]. In this paper, the additive method was selected for comprehensive assessment and calculation follows Equation (4).
where: U(v) is the overall utility of an option, a(i) is the magnitude of the option on the i th attribute, u(i) is the importance of the option on the i th attribute, and n is the number of attributes. One step of MAUT is determining the benefit and non-benefit indicators for the following normalization. In this paper, the beneficial indicators and non-beneficial indicators are normalized by Equations (5) and (6), respectively.
where A n1 is the normalized value of indicator j in scenario i for beneficial values, A n2 is the normalized value of indicator j in scenario i for non-beneficial values, X ij is the values of indicator j in scenario i, Max X j is the maximum value of indicator j in all the scenarios, and Min X j is the minimum value of indicator j in all the scenarios. Because MAUT procedures comprised normalizing the evaluating matrix, the normalized results can be used to respectively obtain the scores of functional, economic, social, and environmental aspects of different scenarios. Based on AHP, MAUT, and several qualitative and quantitative methods, the comprehensive assessment framework used is summarized in Figure 8. This framework is made of four main phases: preparation, design, assessment, and recommendation.
Water 2020, 12, x FOR PEER REVIEW 13 of 35 Figure 8. The framework of comprehensive assessment for five WSUD scenarios. Table 9 shows the weights of different criteria, sub-criteria, and indicators calculated by AHP method. The higher weight of a sub-criteria represents the higher importance for stormwater management in the study area. The CI determined by AHP is 0.0171 (less than 0.1). This means that the AHP results for weight determination are reasonable (see Appendix C).  Table 9 shows the weights of different criteria, sub-criteria, and indicators calculated by AHP method. The higher weight of a sub-criteria represents the higher importance for stormwater management in the study area. The CI determined by AHP is 0.0171 (less than 0.1). This means that the AHP results for weight determination are reasonable (see Appendix C). The social aspect is assigned the highest weight (43.56%) and all the three sub-criteria of the social aspect have relative high weights (15.31% for water reuse function, 13.78% for livability, and 14.47% for water retaining function). This is reasonable because many development goals and strategies of stormwater management in the study area emphasized the improvement of the social values such as water reuse and water retaining. Melbourne Water [83] and Wong [30] reported that the performance of tank and pond is better than bioretention areas and swales in water retaining function. Therefore, their weights are quantitively assigned as 8.68% and 5.79%, respectively. In contrast, the economic aspect is ranked fourth with a 11.91% weight, and the weights of O&M cost and capital cost are very low (1.91% and 2.60%, respectively). Considering the importance of open space for other uses, the weight of land use cost (7.40%) is higher than capital cost and O&M cost. Meanwhile, the weights of the functional aspect (22.73%) and the environmental aspect (21.80%) are similar. Flooding control (12.44%) and water quality improvement (10.29%) are the functions that all stormwater infrastructures must consider regardless of grey or green facilities. The former is given a higher weight because of the severe flooding risks in the area. The carbon emission is assigned 14.78%, which is higher than ecosystem values (7.03%) because "zero-emission" in the water sector has been emphasized many times by stakeholders [92]. It is an essential environmental indicator in the infrastructure construction. Table 10 presents the initial results and normalized results of the functional aspect simulated by MUSIC. Surface runoff, as a non-beneficial indicator, has been transferred to beneficial form before normalization by Equation (3). As the MUSIC software cannot directly calculate the surface runoff, its values were taken from the amount of weir out flow of the South Lawn treatment node. The flooding control and water quality improvement for five scenarios can be further calculated by Equation (1) and are shown in Table 11. It was found that the differences of five scenarios in flooding control are not obvious, except that scenario 2 is slightly better than other alternatives. For the water quality aspect, scenario 5 obtained the highest score in water quality improvement because the large proportion of bioretention areas and swales have pollution removal mechanisms and play a considerable role in improving stormwater quality. Similarly, as an equal grey-green design, the performance of water quality improvement for scenario 4 is also outstanding due to the large size of green-based facilities. On the country, the grey-based scenario's assessment score in water quality improvement is not satisfactory. Scenario 2 performs better in water quality improvement than scenario 3, although the green facility proportion of the latter is higher. This may be because bioretention areas are more efficient than swales in improving stormwater quality. Table 12 summarizes the individual results of eight sub-criteria of economic, social, and environmental aspects and the calculation procedures are provided in the Supplementary Materials. In the calculation, the normal level of stormwater tank and pond is assumed as 1 m. Scenario 1 and 2 have higher capital cost because the cost of stormwater tank is quite high and up to $1000 per kiloliter. Additionally, the two scenarios also perform well in water reuse and retaining function because of the higher performance of these grey facilities in the two indicators, while the bioretention areas and swales cannot contribute to the water reuse. Scenario 3 is relatively average in all the sub-criteria, and there is no one criterion that performs best or worst. Scenario 4 has the highest score among the three economic sub-criteria and also performs satisfactorily in other factors, particularly in carbon emission and ecosystem value. As a green-based option, the pros and cons of scenario 5 are very clear. Its performances in carbon emission, livability improvement, and ecological value are the best among the five alternatives, but the scores of water reuse and retaining function are not satisfactory.

Comprehensive Results
The values of 10 sub-criteria have been determined qualitatively and quantitatively, as shown in Tables 11 and 12. Calculated using Equations (5) and (6), the normalized results of these sub-criteria as well as their weights are presented by spider diagrams (Figure 9). The calculation procedures are enclosed in Appendix D and the Supplementary Materials.
The scores of functional, economic, social, and environm cvental aspects can be calculated respectively by Equation (4), and the overall scores of five scenarios can be further calculated (Table 13). After the comprehensive assessment based on MAUT, scenario 4 (0.771 out of 1), as the equal grey-green strategy, is the optimal option for stormwater management in the University of Melbourne. It provides maximum comprehensive benefits, although it does not perform best in the functional, social, and environmental aspects. It had the highest economic benefits (0.119) because of the smallest areas of WSUD facilities and the relative few stormwater tanks in design. The optimal scenario comprised 250 m 2 stormwater tanks (including 225 m 2 existing tanks), 50 m 2 pond, 170 m 2 bioretention areas, and 150 m 2 swales. The proportion of grey and green facilities are approximately half (52% green facilities). research of Damodram et al. [95] proposed to combine LIDs and BMPs to achieve sustainability goals given limited resources, and its essence is also seeking an optimal combination of grey and green stormwater facilities. Therefore, despite the fact that adopting GSI to manage stormwater has become a worldwide trend, only when a suitable combination of these GSI facilities was applied could their performances be maximized to achieve sustainable development.  In contrast, scenario 1 (0.693 out of 1), as a grey-based strategy with 4.3% green facilities, had a considerable gap for achieving the sustainable stormwater management goals. Its social benefits are excellent but cannot counteract the deficits in terms of functionality and the environmental aspect. Scenario 2 (26.7% green facilities) and scenario 3 (31.7% green facilities) had the approximate 0.75 score but less than scenario 4, which showed the significance of increasing the green WSUD facilities in the study area. If only considering the first four scenarios, we can preliminarily conclude that the overall score increases as the ratio of green facilities gradually increases. Scenario 5 also had the highest environmental and functional benefits and its green facility areas have reached 69.2%, which was mainly attributed by green facilities' effectiveness in water quality improvement and lower carbon emissions. However, it showed poor social benefits and, thus, caused slightly lower overall score than scenario 4 (0.758 and 0.771). This drag-down is due to its poor social benefits that are strongly related to the "water reuse" and "water retaining" services featured in grey facilities. We must notice herein that scenario 5 has the smallest proportion of grey facilities. Therefore, we can conclude that to keep increasing the green facility is not equal to a better GSI combination design, but there is a trade-off relation we must consider between the proportion of green and grey facilities and also among social, economic, environmental, and functional aspects. A proper combination of grey and green WSUD facilities is the best response to the stormwater management in the study area to achieve the sustainable stormwater management goals.
The results are consistent with the research of Alves et al. [93], Gallo et al. [50], Bakhshipour et al. [94], and Damodram et al. [95]. Alves et al. [93] applied MCDA methods to assess the flooding management in three study areas (Marbella, Ayutthaya, and Sukhumvit) and concluded that the green and grey combination measures provided more advantages compared with previously developed methods. Gallo et al. [50] developed a modeling framework that encompasses green and grey stormwater facilities and was applied into the Berkeley neighborhood. They also reported that the optimal stormwater solution in the study area is a mix of green and grey combination. Bakhshipour et al. [94] illustrated that the hybrid green-blue-grey stormwater infrastructures can economically compete with conventional grey-only pipe networks and the green parts can effectively increase the sustainability and environmental friendliness. The research of Damodram et al. [95] proposed to combine LIDs and BMPs to achieve sustainability goals given limited resources, and its essence is also seeking an optimal combination of grey and green stormwater facilities. Therefore, despite the fact that adopting GSI to manage stormwater has become a worldwide trend, only when a suitable combination of these GSI facilities was applied could their performances be maximized to achieve sustainable development.

Insights and Practical Significance of the Case Study
Reviewing the entire case study, it was found that setting the sustainable goals for the specified study area can facilitate the effective WSUD design if water managers or local administration can provide a clear statement of commitment towards stormwater water management. These sustainable goals underpin the proposal and determination of criteria and weights by an assessment framework that serves for WSUD scenarios comparison. We must notice herein that there could be large flexibilities between goals and determination of criteria, indicators, and weights. This paper used the methods of questionnaires, discussion with experts, and literature review to transform the qualitative goals into measurable indicators. These methods are worth recommending but entail local water managers having a good understanding of WSUD benefits, sustainable development goals, and GSI principles. In a real case, water managers should not only tightly follow the local stormwater management goals, but also choose sufficient, suitable, and comprehensive criteria and indicators truly matched with the real situation (e.g., from functional, economic, social, and environmental perspectives).
There are many proposed MCDA methods in previous studies, and the combination of two or more MCDA methods can more easily and practically meet the requirements of decision-makers. In this study, the AHP method helps to determine the weights of different sub-criteria, while MAUT standardizes all the criteria and indicators with or without units. It is anticipated that this combined MCDA assessment tool is more practical and convenient for water managers. In terms of this case study, the framework (Figure 8) well contained the stormwater goals, five different scenarios, combined MCDA methods, and sufficient criteria and indicators. This framework should be promoted for the water sector to solve complex decision-making problems.
Since the proportion of green facilities of scenario 4 is 52% (approximately 50%), the "equal green grey" WSUD design philosophy could be empirically applied in some Melbourne campus. This is mainly because these areas have similar external factors such as land use, stormwater goals, and meteorological data, and the optimal results could be similar. However, for an area with different land uses, policies, and meteorological conditions, it is necessary to conduct a comprehensive assessment as described in this paper. By adjusting the proportions of green-grey facilities to generate a number of scenarios, the model can simulate possible WSUD combinations as alternatives. The method of scenario designs can also be used in many other similar cases where both proportions of green and grey facilities play an important role.

Recommendations and Prospects
The problem that the University of Melbourne urgently needs to solve is related to the implementation of more green stormwater facilities, while the 900 kiloliters existing stormwater tanks basically can meet people's demands of water reusing and retaining. According the results analysis, scenario 4, as "equal grey-green" design, can be applied as a priority design for the University of Melbourne. By comparing scenario 4 and 5, a higher proportion of green WSUD facilities will lower the overall profits. It gives an insight that the performance of WSUD combinations are often not as simple linear relation as the proportional relation, but a more complex mechanism fundamentally associated with criteria, indicators, and weighting systems. Therefore, more scenarios with different green facility proportion are recommended to design for further optimization assessment in the University of Melbourne. In this way, a more accurate optimal result can be further obtained.
As illustrated in Section 3.4, different land uses, policies, and meteorological conditions could lead to various assessment results. In fact, there are many other factors such as stormwater goals, weighting systems, site characteristics, public engagement, and decision-makers' perspectives that should be considered. Taking decision-makers' perspectives as an example, the decision-makers can be engineers, urban planners, and environmentalists who can directly or indirectly impact the assessment of WSUD or GSI [96]. Jayasooriya et al. [47] reported that the priority of the engineers, urban planners, and environmentalists for GSI strategies may be cost-effectiveness, amenity values, and environmental impacts mitigation, respectively. Similarly, for a site with different policies, financial support, location characteristics, public awareness and engagement, etc., the optimal GSI combination may be entirely different. These diversified factors make it a challenge for decision-makers to select the most appropriate stormwater strategy for a site [97].
For the challenge of diversified factors, this paper strongly recommends that a general comprehensive assessment framework can be used to optimize different scenarios based on MCDA methods in Australia or a global scale, rather than the specific case study. In fact, there are several good attempts of combining stormwater management and MCDA to score, rank, and select the optimal option by a more generalizable framework. At the theoretical level, Jayasooriya et al. [47] Sapkota et al. [98], Sapkota et al. [46], and Wu et al. [71] reported that the establishment of a standard comprehensive assessment framework for GSI is significant and it can further assist decision-makers to make a more reasonable, inclusive, and well-informed decision. At the practical level, Morales-Torres et al. [48] reported a decision supported tool called E2STORMED based on MCDA method to assess the multi-benefits of different GSI, and Kuller et al. [31] also developed a rapid GSI-based MCDA tool for WSUD assets (called SSANTO), to assist stakeholders engaged in the urban planning in focusing on opportunities and needs of WSUD. Additionally, related cases regarding comprehensive assessments for GSI have also been published in recent years [42,[99][100][101][102]. Importantly, such a framework can well integrate current or future sustainable strategies by questionnaire or other methods, thereby promoting sustainable development in a site [103]. Specifically, when the sustainable policies and strategies change, the goals to be achieved and the indicators to be assessed will change accordingly [104]. Under various external conditions, the framework will continuously help water managers to make informed and inclusive decisions and play a role in achieving sustainable environmental development involving a variety of factors.
In 2019, Dandy et al. [105] proposed a general comprehensive assessment framework for the assessment of the performance of different stormwater harvesting alternatives based on MCDA (Figure 10), which is similar but more generalizable with the framework described in Figure 8. It has guiding significance for the establishment and promotion of general WSUD or GSI assessment framework for stormwater management. Four main steps shown in Dandy et al.'s framework [105] include problem definition, option design, performance evaluation, and option selection, which could be recommended for GSI assessment framework establishment. The future work will continue to develop the assessment framework on the basis of the four-step theory. As the knowledge in GSI and MCDA progresses, the future research for developing the framework and promoting sustainable stormwater management may include the following aspects. These aspects are also shown and illustrated in a variety of publications.

•
Propose a clear statement of commitment towards sustainable stormwater water management based on different site characteristics [104].

•
Develop a better methodology of transforming the qualitative goals into measurable indicators except for questionnaire [106].

•
Develop more reasonable combined scenarios of GSI or WSUD [107].

•
Develop more accurate and indicator calculation methods [106].

•
Consider the perspectives of all the stakeholders including the public engagement rather than experts only [108,109]. • Select more appropriate MCDA methods for water managers in decision-making processes [110].
assessment framework can be used to optimize different scenarios based on MCDA methods in Australia or a global scale, rather than the specific case study. In fact, there are several good attempts of combining stormwater management and MCDA to score, rank, and select the optimal option by a more generalizable framework. At the theoretical level, Jayasooriya et al. [47] Sapkota et al. [98], Sapkota et al. [46], and Wu et al. [71] reported that the establishment of a standard comprehensive assessment framework for GSI is significant and it can further assist decision-makers to make a more reasonable, inclusive, and wellinformed decision. At the practical level, Morales-Torres et al. [48] reported a decision supported tool called E2STORMED based on MCDA method to assess the multi-benefits of different GSI, and Kuller et al. [31] also developed a rapid GSI-based MCDA tool for WSUD assets (called SSANTO), to assist stakeholders engaged in the urban planning in focusing on opportunities and needs of WSUD. Additionally, related cases regarding comprehensive assessments for GSI have also been published in recent years [42,[99][100][101][102]. Importantly, such a framework can well integrate current or future sustainable strategies by questionnaire or other methods, thereby promoting sustainable development in a site [103]. Specifically, when the sustainable policies and strategies change, the goals to be achieved and the indicators to be assessed will change accordingly [104]. Under various external conditions, the framework will continuously help water managers to make informed and inclusive decisions and play a role in achieving sustainable environmental development involving a variety of factors.
In 2019, Dandy et al. [105] proposed a general comprehensive assessment framework for the assessment of the performance of different stormwater harvesting alternatives based on MCDA (Figure 10), which is similar but more generalizable with the framework described in Figure 8. It has guiding significance for the establishment and promotion of general WSUD or GSI assessment framework for stormwater management. Four main steps shown in Dandy et al.'s framework [105] include problem definition, option design, performance evaluation, and option selection, which could be recommended for GSI assessment framework establishment. The future work will continue to develop the assessment framework on the basis of the four-step theory. As the knowledge in GSI and MCDA progresses, the future research for developing the framework and promoting sustainable stormwater management may include the following aspects. These aspects are also shown and illustrated in a variety of publications.


Propose a clear statement of commitment towards sustainable stormwater water management based on different site characteristics [104].  Develop a better methodology of transforming the qualitative goals into measurable indicators except for questionnaire [106].  Develop more reasonable combined scenarios of GSI or WSUD [107].  Develop more accurate and indicator calculation methods [106].  Consider the perspectives of all the stakeholders including the public engagement rather than experts only [108,109].  Select more appropriate MCDA methods for water managers in decision-making processes [110].

Limitations
The limitations of this paper included MUSIC modeling restrictions, subjectivity of experts in the AHP method, and uncertainty of indicator values. For the MUSIC modeling design, firstly, the campus was divided into seven sub-areas to build the model in MUSIC as shown in Figure 1, where it was assumed that these separate areas are independent. In fact, the seven sub-areas constitute the whole campus and they can influence each other. Additionally, the rainfall events used in the model is of a one-day interval, while a long-term and high intensity of rainfall data (e.g., one-hour interval) are preferable for assessment. Meanwhile, the results produced from MUSIC do not consider the underground pipes and it assumed that these existing conventional approaches have no influence on the facilities performance. Finally, MUSIC is a software only applicable to Australia and New Zealand, and the regional limitations will make it difficult to extend to the global scale [71]. Subjectivity of experts is another limitation of the study. Despite the fact that this study invited a total of six experts to conduct the questionnaire and peer view, this sample size is not sufficient to eliminate the subjectivity of influence. When conditions permit, it would more accurate if more numbers of relevant experts and scholars can participate in the weight determination. For the indicator values, many unit values of indicators such as carbon emission, construction cost, O&M cost, and livability are extracted from previous research and may be out of date. Also, the qualitative assessment method of indicators will also cause some errors.

Conclusions
It is a worldwide trend for city planners to select a GSI philosophy rather than "rapid-draining"-based traditional approaches in stormwater management. In Australia, WSUD facilities are designed to accelerate economic, social, and environmental development while dealing with a series of stormwater functional problems. Current stormwater management assessments sometimes ignore the economic and social benefits and there is a lack of a comprehensive selection framework to optimize the stormwater facility combinations based on various criteria. This paper focuses on assessing and optimizing WSUD facility combinations based on AHP and MAUT through a case study of Melbourne University's Parkville Campus. Based on the results, the key findings and conclusions are summarized as follows.
Scenario 4 (equal grey-green design) containing 52% green WSUD facilities obtained the highest scores (0.771) among five designed scenarios. It provides maximum comprehensive benefits although it does not perform best in the functional, social, and environmental aspects. It can be applied as a priority design for achieving the goals of stormwater management in the University of Melbourne. Implementing more green stormwater facilities is an urgent issue in the study area. Meanwhile, considering the similar land uses, stormwater goals and meteorological data in most of Melbourne campuses, the "equal green grey" WSUD design philosophy could be empirically applied in these areas.
Scenario 5 came second with a score of 0.758 despite the fact that the green facility proportion reached 69%. It indicates that to keep increasing the green facility is not equal to a better GSI combination design, but there is a trade-off relation we must consider between the proportion of green and grey facilities and also among social, economic, environmental, and functional aspects. Therefore, more scenarios with different green facility proportion are recommended for further optimization assessment to obtain a more accurate result.
The methods of questionnaires, discussion with experts, and literature review can transform the qualitative goals into measurable indicators, but it is inevitable that the subjectivity of experts will cause some errors. It would more accurate if more numbers of relevant experts and scholars can participate in the weight determination. Also, the combination of two or more MCDA methods can more easily and practically meet the requirements of decision-makers. It is anticipated that this combined MCDA assessment tool is more practical and convenient for water managers.
It is highly possible that a general comprehensive WSUD or GSI assessment framework can be used to optimize different scenarios based on MCDA methods in Australia or a global scale. In recent years, there are several good attempts of assessments combining stormwater management and MCDA in a framework from the theoretical level to the practical level. Four steps including problem definition, option design, performance evaluation, and option selection are recommended to form the GSI assessment framework structure. The future work will continue to develop GSI assessment framework based on the four-step theory. When such a more applicable framework is broadly promoted, it can further promote sustainable development by helping water managers to make informed and inclusive decisions involving a variety of factors.

Appendix A.2 Pond
Pond is a simple basin for open water body, which is quite useful for water temporary storage, reuse, and extended detention [113]. At the same time, pond removes particles in inflow by gravity [114]. Figure A2 presents the scheme of a pond and

. Pond
Pond is a simple basin for open water body, which is quite useful for water temporary storage, reuse, and extended detention [113]. At the same time, pond removes particles in inflow by gravity [114]. Figure  A2 presents the scheme of a pond and Table A2 summarizes the design parameters of two types of ponds in MUSIC.

. Bioretention Areas
Bioretention areas treat stormwater by vegetated filter media to remove pollutants, and at the same time store stormwater and allow infiltration into the ground [115]. The vegetation types and media properties can clearly impact bioretention areas' performance [116]. A typical example of a bioretention system is a rain garden. In MUSIC modeling, three types of bioretention areas with different sizes were

Appendix A.3 Bioretention Areas
Bioretention areas treat stormwater by vegetated filter media to remove pollutants, and at the same time store stormwater and allow infiltration into the ground [115]. The vegetation types and media properties can clearly impact bioretention areas' performance [116]. A typical example of a bioretention system is a rain garden. In MUSIC modeling, three types of bioretention areas with different sizes were selected. We assume that vegetation and media properties are the same and all bioretention systems are square. Figure A3 and Table A3 show the schematic diagram and design parameters of bioretention areas, respectively.

Appendix A.4 Swales
Vegetated swales are open water channels covered with vegetation. Swales can remove suspended solids effectively by sedimentation and vegetation trapping [117]. The performance of vegetated swales can be impacted by vegetation density and heights. In MUSIC modeling, the cross section of swales is assumed to be trapezoid. Two types of swales including 50 m 2 size and 100 m 2 size were selected for design. Figure A4 and Table A4 show the schematic diagram and design parameters of swales, respectively. selected. We assume that vegetation and media properties are the same and all bioretention systems are square. Figure A3 and Table A3 show the schematic diagram and design parameters of bioretention areas, respectively.  Vegetated swales are open water channels covered with vegetation. Swales can remove suspended solids effectively by sedimentation and vegetation trapping [117]. The performance of vegetated swales can be impacted by vegetation density and heights. In MUSIC modeling, the cross section of swales is assumed to be trapezoid. Two types of swales including 50 m 2 size and 100 m 2 size were selected for design. Figure  A4 and Table A4 show the schematic diagram and design parameters of swales, respectively.    This questionnaire is to determine the relative importance for 10 sub-criteria in terms of Melbourne University's stormwater management in the future. The questionnaire results will be used as the raw data of AHP (analytic hierarchy process) method for Hanxiang Xiong, Yafei Sun, and Xinwei Ren's paper. Please refer to the examples below, compare the sub-criteria's relative importance of table's each row and column, and use numbers (from 1 to 9) to represent different relative importance (see Table A5) in terms of stormwater management development in the study area (Fractions are ALLOWED to use to express more accurately). More information about site's stormwater goals and strategies are summarized as follows. Table A5. Scale of relative importance for 1 to 9.

Number
Scale of Relative Importance If you think flooding control is of extreme importance than capital cost, write the number "9" in the form as shown in Table A6. Table A6. Writing example (flooding control is of extreme importance than capital cost).

ATTRIBUTE Flooding Control Capital Cost
Flooding control 1 9 Capital cost 1 If you think capital cost is of between equal importance and moderate importance than flooding control, write the number "2" or other fractions in the form as shown in the Table A7. As the value of CI is less than 0.1, the initial pair-wise comparison matrix determined my questionnaires are reasonable. Therefore, the following comprehensive assessment can be further conducted.