Next Article in Journal
Comparative Multi-Criteria Assessment of Hydrological Vulnerability—Case Study: Drainage Basins in the Northeast Region of Romania
Previous Article in Journal
An Improved Transfer Learning Model for Cyanobacterial Bloom Concentration Prediction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Runoff Reduction Effects at Installation of LID Facilities under Different Climate Change Scenarios

1
Chungbuk Research Institute, Cheongju 28517, Korea
2
Department of Agricultural Rural Engineering, Chungbuk National University, Cheongju 28644, Korea
3
Faculty of Agricultural Engineering and Technology, PMAS Arid Agriculture University, Rawalpindi 44000, Pakistan
4
Water Resources Management Research & Planning Department, Korea Water Resources Corporation, Daejeon 34350, Korea
*
Author to whom correspondence should be addressed.
Water 2022, 14(8), 1301; https://doi.org/10.3390/w14081301
Submission received: 12 March 2022 / Revised: 12 April 2022 / Accepted: 14 April 2022 / Published: 16 April 2022
(This article belongs to the Topic Water Management in the Era of Climatic Change)

Abstract

:
In order to prepare for floods and droughts that occur as a result of climate change, various studies in water-related fields are being carried out in various countries around the world. Among them, special attention is being paid to the low-impact development (LID) technique. This study measured the annual maximum daily rainfall data from 2019 to 2100 by climate change scenario, which is the annual maximum daily rainfall series observed for rainfall stations, and tested the hydrological data using statistical analysis. After determining whether the data could be analyzed, the probability distribution was selected, and the parameters of the selected probability distribution were calculated using the L-moment method for each rainfall station. The probabilities of rainfall data were derived using GEV distribution, and the United States Environmental Protection Agency Storm Water Management Model (SWMM), a runoff simulation program, was used to compare and analyze the runoff reduction rate before and after the installation of a permeable pavement as an LID facility. The results of the analysis showed that representative concentration pathway (RCP) 4.5 and RCP 8.5 had the effect of reducing the runoff for more than 100 years at a 30% reduction rate compared with before installation.

1. Introduction

The frequency of typhoons and heavy rains has increased rapidly in Korea due to global warming and abnormal weather conditions, and the scale of flood damage has increased due to climate change and urbanization, resulting in substantial property damage in urban areas. Unlike the flooding of agricultural land, urban flooding amplifies not only economic and human loss, but also the psychological damage and social anxiety of urban residents because, when flooding occurs, many people and a large amount of urban infrastructure are concentrated in a dense space. In Korea, the average temperature recorded in 2016 was the highest since 1973, and there are many signs of a changing climate compared with the past, such as heat waves as well as unexpectedly heavy rains in July and August.
According to the Intergovernmental Panel on Climate Change (IPCC) report [1], the average global temperature will rise by 1 °C in the 2020s, and up to 1.7 billion people worldwide will suffer from water shortages. In the long term, over 30% of coastal areas will be lost due to rising sea level in the 2080s. It has been reported that more than 20% of the world’s population will be at risk from flooding [2].
In preparation for many water-related disasters, such as floods and droughts due to climate change, various studies in water-related fields are being conducted in countries around the world. Among them, much attention is being paid to the low-impact development (LID) technique. LID is a technique that aims to establish water circulation systems similar to those prior to development by installing permeable elements in impervious areas to reduce runoff and improve water quality. Developed countries are responding to the problems caused by climate change and urbanization by applying it to urban areas. LID facilities include bio-retention cells, rain gardens, green roofs, infiltration trenches, permeable pavements, rain barrels, and vegetative swales. Therefore, a permeable pavement is one of the facilities selected to reduce the depletion of groundwater and the impacts of urban flooding.
In this study, by analyzing the changes in the observed rainfall and the probability of rainfall by applying climate change scenarios, by setting the design rainfall for each return period using the probability of rainfall for each scenario, and by simulating the amount of runoff before and after the installation of an LID facility in the object region, we could use the frequency, according to the scenario, to compare and analyze the emission reduction rate. In addition, by comparing the amount of runoff in the return period for each climate change scenario, we sought to show that installing LID facilities has the effect of reducing runoff. Therefore, it is judged that the results analyzed through this study can aid in setting an appropriate design direction when considering an LID facility.

2. Theoretical Analysis

For the case of the current facility design return period, rainfall analysis and return period analysis using rainfall observation data was conducted. However, urban runoff has continued to increase over a short period of time due to the increase in impervious areas, and the return period of damage, accompanied by an irregular climate and locally heavy rainfall, has continued to occur. Therefore, in order to supplement the existing design return period and to design a stable structure, it was necessary to predict the future climate using observational data and utilize the data composed of scenarios, rather than analyze past data only. Therefore, it is suggested that the application of the LID facility in the climate change scenario used in this study be used as the basic premise for disaster prevention standard guidelines, and in measures to respond to climate change in urban areas in the future.

2.1. Climate Change Scenario

For the climate change scenario used in this study, the greenhouse gas concentration is determined as the amount of radiation exerted by human activities on the atmosphere as per the IPCC 5th evaluation report, and the representative concentration pathway (RCP) scenario contains one fixed representative radiative forcing value, with the expression “representative” used in the sense that there can be many socio-economic scenarios. Unlike the existing Special Report on Emission Scenarios (SRES), the RCP scenario reflects the recent trend of changes in greenhouse gas concentrations and has been updated to fit the recent prediction model. The four representative greenhouse gas concentrations in the RCPs are 2.6, 4.5, 6.0, and 8.5. In the process of calculating the greenhouse gas concentrations, the social and economic assumptions were changed from the basis of a future society structure to whether or not to implement climate change response policies.
The Korea Meteorological Administration (KMA) simulated the future climate change scenario by introducing the RCP scenario as forced input data to the Hadley Center Global Environmental Model version 2—Atmosphere and Ocean (HadGEM2-AO) model, a global climate change prediction model, and used the simulated global climate change scenario as input data to the HadGEM3-RA model, a regional climate model. Through epidemiological detailing, regional climate model (RCM) data, a regional climate change scenario that well reflects the effects of complex topography that cannot be expressed by the global model, were calculated. In this study, RCP 4.5 and RCP 8.5 scenarios were selected and analyzed among the four climate change scenarios.

2.2. Return Period Analysis

There are many different probability distribution types for return period analysis according to climate change scenarios, such as Gumbel (GUM), generalized extreme value (GEV), generalized logistic (GLO), generalized Pareto (GPA), generalized normal (GNO), and Pearson type 3 (PT3). The relationship between the probability weighted moment [3] and the L-moment is described by Hosking [4,5], Maeng et al. [6], Maidment [7], and the World Meteorological Organization [8]. Furthermore, Hosking [4] published that the L-moment, which is the linear combination of the statistical characteristics of the probability weight moment (PWM) based probability distribution, enables efficient and safe parameters in smaller samples more than other moment methods and method of maximum likelihood. In addition, the parameter estimation method by L-moment for each probability distribution has been studied and described by Hosking et al. [4,9] and Maidment [7].

2.3. EPA-SWMM for Application of LID Facility

In this study, the Storm Water Management Model (SWMM) was selected and analyzed from among the models for estimating the amount of flooding caused by rainfall in urban watersheds. In 1971, with the support of the United States Environmental Protection Agency (US EPA), the Metcalf & Eddy Corporation, and the University of Florida and Water Resource Engineering (WRE), the model was developed to simulate the flow rate and water quality of the urban watershed sewage system [10]. In addition, in the 1981 SWMM model, the extended transport (EXTRAN) block, designed to calculate the overflow, drainage, and pressure flow of hand structures, was added to the SWMM model to expand and supplement the transport block. Currently, the EPA-SWMM has been developed up to version 5.1.007 (with LID control).
The LID runoff model using the SWMM was mentioned by Lee [11], and the applied model of the study site was selected and analyzed as a permeable pavement.

3. Study Area and Research Method

3.1. Selection of Study Area and Overview

The watershed of the study area was the Eco-Delta City (EDC) located in Busan city in Korea. The location of the study area is shown in Figure 1. The total area was 105.9 ha (1.059 km2). The West Nakdong River on the west flows from north to south, and the Pyeonggangcheon Stream on the east flows into the West Nakdong River through the Suna sluice gate. The study area was the estuary of the West Nakdong River, which consists of flat land with an elevation of 50 m or less (mostly 0.4 to 4.0 m), and the difference in height between the river and the land is less than 1.0 m, being mostly lowlands. In addition, agricultural water is supplied to the existing irrigation channel in study area through the lower culvert of the Namhae Expressway, but it will be closed due to development of the study area, and there is no external inflow.
Therefore, the study area was selected with regard to the topographical condition of the area and the design stage before the installation of the LID facility, as well as the fact that the permeable pavement will have a large reduction effect on the flat land.

3.2. Collection of Hydrological Data

The rainfall station was selected according to the distance from the study area. Therefore, the Gimhae rainfall station was selected, as it is located close to Gimhae-si, with the data exceeding 30 years. The annual maximum daily rainfall data from the Gimhae rainfall station were used in this study [12]. The location of the Gimhae rainfall station is shown in Figure 2.
In addition, the Korean Peninsula climate change scenario was analyzed using the global climate change scenario data of the HadGEM2-AO model, and was developed in Korea using the HadGEM3-RA model (Hadley Center Global Environment Model version 3 atmosphere regional climate model) at a resolution of 12.5 km.
Production process for climate change scenarios is shown in Figure 3.
In this study, the RCP 4.5 scenario was simulated for the peninsula for 200 years with 12.5 km resolution, and the integral control was based on the RCP 8.5 scenario. The extraction of annual maximum daily rainfall data was based on the past observed time series of daily rainfall (1988–2018) and the climate change scenarios from 2019 to 2100.
The annual maximum daily rainfall for the next 100 years suitable for the EDC, the research study area, was analyzed and used as shown in Table 1.

3.3. Input Data Composition of SWMM Model

The parameters of the SWMM model are watershed area, permeability, Manning’s roughness coefficients, surface storage volume, storage depth, CN value, slope, and sewage pipe data. The rainfall runoff is simulated using these parameters as input data. However, since there is limited observation data to apply all parameters, in this study, watershed area, permeability, Manning’s roughness coefficients, slope, and sewage pipe data were used as parameters.
In Korea, the EPA-SWMM is used for flood analysis of LID facilities, and in the case of LID facility effect analysis, the most accurate and ideal result is obtained by performing monitoring under the same conditions before and after installation of the facility and analyzing the effect. In this section, the parameters and data construction before and after LID facility installation are presented. Watershed delineation is the first step in simulating the physical drainage systems. In the EPA-SWMM model, the subwatershed is assumed to be a rectangle with uniform characteristics (slope, roughness, etc.). The shape of the watershed is defined by factors such as area, watershed width, slope, etc. In the case of the EDC area, it flows out through three outlets, and in the case of each watershed, the model was built by unifying the watershed names A, B, and C, as shown in Figure 4, with all three watersheds flowing into the West Nakdong River. The watersheds were divided into subwatersheds according to the pipe network, and there were 133 subwatersheds in A, 171 in B, and 202 in C. The pipe network, according to watershed area, subwatersheds, outlets, nodes, and links, is presented in Table 2 and Figure 4.
To introduce the LID facilities, variables for each facility were applied. Table 3 shows the designated values of the permeable pavement in this study.

3.4. Research Methodology

The following procedure was followed to analyze the reduction rate of the runoff of the LID facility installed in selected watersheds.
The extraction of annual maximum daily rainfall data was based on the past observed time series of daily rainfall (1988–2018) and the climate change scenarios from 1988 to 2050 and 1988 to 2100. The statistical tests, such as the mean, standard deviation, coefficient of variation, coefficient of skewness, and coefficient of kurtosis, were applied to hydrological data. After calculating the basic statistical values of the annual maximum daily rainfall series, the time series were tested for independence, homogeneity, and outliers. An appropriate distribution was selected for the GUM, PT3, GEV, GLO, GPA, and GNO distributions through the goodness-of-fit test. The parameters of the appropriate distribution were computed according to the study area and the period of analysis using the L-moment method, and the probability of rainfall was estimated using the selected appropriate distribution. Probability of rainfall was computed according to the climate change scenarios.
A rainfall intensity duration equation was derived, and the regression equation was calculated using the precipitation data from the selected rainfall station. The rainfall distribution of the probability of rainfall was constructed using the calculated regression equation. In the EPA-SWMM, the same pipe network data, parameters, and rainfall data before and after the installation of the LID facility were established, and the parameters for the facility after the installation of the LID facility were additionally configured. The effect analysis was presented by estimating the probability of rainfall for each return period of climate change scenario, such as RCP 4.5 and RCP 8.5, by simulating runoff and comparing the reduction rate.
Figure 5 is a schematic diagram of the probability precipitation calculation considering climate change scenarios and LID facilities. The flow chart used in the study based on the procedure introduced in Figure 5 is shown in Figure 6.

4. Results

4.1. Basic Statistical Analysis and Calculation of Probable Rainfall According to Climate Change Scenario

The mean, standard deviation, coefficient of skewness, coefficient of variation, and coefficient of kurtosis for the annual maximum daily rainfall data from the Gimhae rainfall station were calculated. The mean and standard deviation for each period of the RCP 4.5 scenario ranged from 126.2 m3/day to 135.0 m3/day and 48.502 to 53.361, respectively, the coefficient of skewness and coefficient of variation were from 1.162 to 1.618 and 0.383 to 0.408, respectively, and the coefficient of kurtosis was from 1.616 to 3.769. This indicated overall that the mean was larger than the range of standard deviation, the coefficient of skewness and the coefficient of variation showed positive values and were biased to the right, and the coefficient of kurtosis was larger than 3, the standard value of normal distribution.
In the RCP 8.5 scenario, the mean and standard deviation ranged from 121.2 to 135.0 and 50.727 to 56.491, respectively, the coefficient of skewness and the coefficient of variation were from 1.231 to 2.255 and 0.395 to 0.461, respectively, and the coefficient of kurtosis was from 1.616 to 6.828.
The basic statistics for RCP 4.5 and RCP 8.5 are shown in Table 4 and Table 5.

4.2. Independence, Homogeneity, and Outlier Detection in Annual Maximum Daily Rainfall

The climate change scenario data were checked to determine the existence of independence, homogeneity, and outliers. The Wald–Wolfowitz test was applied to check the independence, and the Mann–Whitney test was applied to observe the homogeneity of the data. In addition, the Grubbs–Beck test was used for the detection of outliers in the time series data.
For the independence test, the Wald–Wolfowitz test [14], which is a non-parametric test that tests the independence of a population, was performed, and it was found that there were no abnormalities in any of them.
The outlier test induces inappropriate statistical parameters in the case of data that appear far above or below the general balanced distribution of hydrological data, resulting in uncertainty in the presentation of the design hydrologic quantity. Therefore, the presence or absence of outliers was tested using the Grubbs–Beck method for hydrological data of the daily maximum flow series for each analysis in the region to which the climate change scenario was applied. It was confirmed that there were no outliers.
Therefore, the annual maximum daily rainfall in the study area was recognized as valid for analysis as hydrological data.

4.3. Estimation of L-Moment Ratio

L-moment ratios were estimated for goodness-of-fit test evaluation of the probability distributions, such as GUM, PT3, GEV, GLO, GPA, and GNO. L-skewness and L-kurtosis were calculated for the 1988–2018, 1988–2050, and 1988–2100 periods using the annual maximum daily rainfall of RCP 4.5 and RCP 8.5.
In the case of RCP 4.5, L-skewness was 0.2535, 0.2136, and 0.2510, and L-kurtosis was 0.2108, 0.1884, and 0.1999. In the case of RCP 8.5, L-skewness was 0.2535, 0.2822, and 0.3321, and L -kurtosis was 0.2108, 0.2019, and 0.2538.

4.4. Goodness-of-Fit Test

For the goodness-of-fit test, L-moment ratio and Kolmogorov–Smirnov (K–S) tests were performed. First, the selection of an appropriate probability distribution for each period was made by plotting the L-moment ratio, using the L-moment ratio calculated in the previous section. The L-moment ratio diagram was plotted to select the best fitted probability distribution from among GUM, PT3, GEV, GLO, GPA, and GNO. The best fitted probability distribution is the one that follows the observed data.
Therefore, Figure 7 shows the average values of the L-moment ratio for the RCP 4.5 and RCP 8.5 scenarios.
The dimensionless L-moment ratios scattered in the diagram are more likely to follow the curve of the GEV distribution, as shown in Figure 7. Therefore, the L-moment method application on the annual maximum daily rainfall for each rainfall station and its plotting in the L-moment ratio diagram showed that the GEV distribution was an appropriate probability distribution compared with GUM, GLO, GPA, GNO and PT3.
Second, the K–S test was applied to choose the best fitted probability distribution for each period of the annual maximum daily rainfall series of the climate change scenario. As a result, at the 5% significance level, data for each period of the annual maximum daily rainfall series were recognized as following the distribution of GEV, GUM, GLO, GPA, GNO, and PT3.
The L-moment ratio diagram and the K–S goodness-of-fit tests showed that the GEV distribution was best fitted probability distribution compared with the other probability distributions. Therefore, the GEV distribution was selected for further analysis.

4.5. Parameter Estimation of Desired Distribution According to L-Moment Method

The parameters of the GEV distribution were estimated by applying it to the annual maximum daily rainfall of the climate change scenarios of RCP 4.5 and RCP 8.5. The parameters of the GEV distribution consist of shape, scale, and location, estimated using the L-moment method for defined periods. The parameters of the GEV distribution that were estimated for the RCP 4.5 and RCP 8.5 scenarios are shown in Table 6 and Table 7.

4.6. Computation of Return Periods Based on Annual Maximum Daily Rainfall

Return periods were computed through the application of GEV distribution on annual maximum daily rainfall of the climate change scenario using the parameters computed in the section above, as shown in Table 8 and Table 9.

4.7. Analysis of Runoff before and after Application of LID Facilities by Climate Change Scenario

As a result of the analysis of the reduction rate for RCP 4.5 and RCP 8.5 during the 1988–2018 period, it was found that the reduction rate for the total amount of runoff increases as the return period increases.
In addition, as a result of the analysis of the reduction rate from 1988 to 2050 for the RCP 4.5 scenario, the 10-year return period was 45.4%, the 20-year return period was 45.8%, the 50-year return period was 46.2%, the 100-year return period was 46.5%, and the 200-year return period was 46.7%. It was found that the greater the return period, the greater the reduction rate.
Analyzing the RCP 4.5 scenario as a sample, the reduction rate from 1988 to 2100 decreased from 321.8 m3/day to 175.7 m3/day in the 10-year return period of the total runoff, showing a reduction rate of 45.4%. In the 20-year return period, it decreased from 386.8 m3/day to 209.6 m3/day, showing a reduction rate of 45.8%; and in the 50-year return period, it decreased from 481.1 m3/day to 258.4 m3/day, showing a reduction rate of 46.3%. In the case of the 100-year return period, it decreased from 559.9 m3/day to 299.0 m3/day, a reduction rate of 46.6%; and in the 200-year return period, it decreased from 646.2 m3/day to 343.3 m3/day, showing a reduction rate of 46.9%.
Table 10 and Figure 8 show the results of the runoff and reduction rates according to the maximum daily rainfall by return period for the 1988–2100 period analysis of the RCP 4.5 scenario.
As a result of the analysis of the reduction rate for the RCP 8.5 scenario analysis from 2018 to 2050, the 200-year return period decreased from 677.2 m3/day of the total runoff to 359.1 m3/day, indicating the largest reduction rate of 47.0%.
Analyzing the RCP 8.5 scenario as a sample, the reduction rate from 1988 to 2100, during the RCP 8.5 scenario analysis period, for the 10-year return period of the total runoff decreased from 313.1 m3/day to 171.2 m3/day, showing a reduction rate of 45.3%. In the 20-year return period, it decreased from 390.8 m3/day to 211.7 m3/day, showing a reduction rate of 45.8%; and in the 50-year return period, it decreased from 515.3 m3/day to 276.0 m3/day, showing a reduction rate of 46.4%. In the case of the 100-year return period, it decreased from 629.9 m3/day to 334.9 m3/day, a reduction rate of 46.8%; and in the 200-year return period, it decreased from 766.1 m3/day to 404.5 m3/day, showing a reduction rate of 47.2%.
Table 11 and Figure 9 show the results for the runoff and reduction rates according to the maximum daily rainfall by return period for the 1988–2100 period analysis of the RCP 8.5 scenario.

4.8. Analysis of Reduction Rate by Climate Change Scenario

To efficiently manage runoff, it is necessary to determine the rainfall (design rainfall) to estimate the size of the reduction potential facility. A pre-disaster impact assessment is recommended [15]. The usual return period used in the design of temporary and permanent facilities is 30 and 50 years, respectively. As a result, a reduction rate return period analysis based on the 20-year runoff was reported in this study for LID facilities. The return period analysis from 1988 to 2018 was the same for the RCP 4.5 and RCP 8.5 scenarios; hence, it is reported only for the RCP 4.5 scenario.
In the RCP 4.5 scenario, as a result of the comparison and analysis of the total runoff from 1988 to 2018, it was found that the runoff on a 20-year return period of 421.1 m3/day was greater than the 200-year return period of 376.1 m3/day before and after the LID facility installation, which was predicted to reduce the runoff for up to 200 years, as a result of the LID facility installation, for the 20-year return period, when designing the return period.
In addition, a comparison and analysis of the total runoff from 1988 to 2050 was conducted, and the results showed that 380.1 m3/day for the 20-year return period before the LID facility installation was higher than 316.7 m3/day for the 200-year return period after the LID facility installation, which would have the effect of reducing runoff up through the 200-year return period. The results of the comparison and analysis of the total runoff from 1988 to 2100 showed 386.8 m3/day for the 20-year return period before the LID facility installation, and 343.3 m3/day for the 200-year return period after the LID facility installation. This showed a value higher than 343.3 m3/day, and it was predicted that the runoff at the 20-year return period set would be reduced to the 200-year return period as a result of the LID facility installation, when designing the return period.
Table 12 and Figure 10 show the results of the comparison and analysis of the total runoff for the RCP 4.5 scenario by return period before and after the LID facility installation.
In the case of the RCP 8.5 scenario, the analysis of the 20-year return period from 1988 to 2050 before the LID facility installation was based on the comparison and analysis of the total runoff from 1988 to 2050. The return period of 380.0 m3/day was higher than the return period of 359.1 m3/day for the 200 years after LID facility installation. When designing the return period, it was found that the amount of runoff, with a return period of 20 years, would be reduced for around 200 years or more due to the LID facility installation.
The analysis of the 20-year return period from 1988 to 2100 before the LID facility installation was based on the comparison and analysis of the total runoff from 1988 to 2100. The return period of 390.8 m3/day was higher than the return period of 334.9 m3/day for the 100 years after the LID facility installation, and smaller than the return period of 404.5 m3/day for 200 years. When designing the return period, it was found that the amount of runoff with a return period of 20 years would be reduced for around 100 years or more due to the LID facility installation.
Table 13 and Figure 11 show the results of the comparison and analysis of total runoff for the RCP 8.5 scenario by return period before and after the installation of the LID facility.

5. Conclusions

The runoff reduction rate by return period according to the scenario was compared and decreased in this study by setting the rainfall distribution by return period via an analysis of the change in the probability of rainfall. This was achieved by applying the climate change scenario and the probability of rainfall for each scenario, and simulating the amount of runoff before and after the LID facility installation in the object region.
The above analysis results are summarized as follows:
  • After composing the annual maximum daily rainfall data based on the previously observed annual maximum daily rainfall series (1988–2018) and climate change scenarios from 1988 to 2050 and 1988 to 2100 (provided by the Korea Meteorological Administration), statistical analysis was undertaken to determine its suitability as hydrological data.
  • The GEV distribution was found to be more appropriate than the other five probability distributions applied to the series of annual maximum daily rainfalls per scenario, according to the analysis after performing the L-moment ratio and K–S tests, which are goodness-of-fit tests.
  • The reduction rate was analyzed by constructing a scenario that considered the same parameters and rainfall data before and after LID facility installation, and the LID parameters after installation. Assuming the same cost aspect in installing the LID, we focused on analyzing and selecting the rainwater reduction rate. Regarding RCP 4.5 and RCP 8.5, we observed a maximum reduction rate of 46.9% for the 200-year return period in the 1988–2100 analysis period for RCP 4.5, and the highest reduction rate of 47.2% for the 200-year return period in the 1988–2100 period for RCP 8.5.
  • As a result of analyzing the reduction rate return period for each scenario, in the case of scenario RCP 4.5, in the analysis of the 1988–2018, 1988–2050 and 1988–2100 periods, the 20-year return period before installation of the LID facility was the same as the amount of runoff after installation of the LID facility. In the analysis of the 1988–2050 period, it was found that the effect of reducing the runoff return period would remain for around 200 years, and that the runoff return period of 20 years before the installation of the LID facility would show the effect of reducing runoff for more than 100 years after the installation.
  • Even in the case of scenario RCP 8.5, in the 1988–2050 analysis period, the amount of runoff with a return period of 20 years before the installation of the LID facility had the effect of reducing the amount of runoff for up to 200 years after installation. In the 1988–2100 analysis period, the amount of runoff with a return period of 20 years before the installation of the LID facility was effective in reducing the amount of runoff for 100 years or more.
This study presented the results of analyzing the reduction rate by return period for the design of LID facilities. In the case of Korea, it is difficult to improve efficiency after the installation of LID facilities due to a lack of basic data on LID specifications and an analysis of their effects. In comparison with other regions, there are many things to consider, such as differences in parameters, rainfall amounts, and land use for each region, and there are limits to construction based only on the related literature. Therefore, analysis through reduction research and actual monitoring data should be continuously undertaken. In addition, this study was limited in that the rate of change in the runoff reduction according to the timing of the permeable pavement and the selected LID facility was not considered.
However, as described in Section 2, the existing studies analyzed only the observational data, but this study considered a new analysis in that the LID facilities used data on climate change scenarios.

Author Contributions

Conceptualization, S.L.; methodology, S.L., S.M.; software, S.L., M.A.; validation, S.L., D.K., M.A.; formal analysis, D.K., B.L.; investigation, D.K., B.L.; resources, M.A., D.K., B.L.; writing—original draft preparation, S.L., S.M.; writing—review and editing, M.A., D.K.; visualization, S.L., D.K.; supervision, S.L., S.M.; project administration, S.L., S.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Korea Agency for Infrastructure Technology Advancement (KAIA) grant funded by the Ministry of Land, Infrastructure and Transport (Grant number 21RITD-C161134-01.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. National Institute of Meteorological Sciences (NIMS), South Korea. Global Climate Change Report 2012 for IPCC 5th Assessment Report: Climate Change Projections According to RCP 2.6/4.5/6.0/8.5; National Institute of Meteorological Sciences: Seoul, Korea, 2012. (In Korean)
  2. Ministry of Land, Infrastructure and Transport, South Korea. Plan of Securing for National Water Security against Climate Change, 1st ed.; Ministry of Construction and Transportation: Seoul, Korea, 2008. (In Korean)
  3. Greenwood, J.A.; Landwehr, J.M.; Matalas, N.C.; Wallis, J.R. Probability weighted moments: Definition and relation to parameters of several distributions expressed in inverse form. Water Resour. Res. 1979, 15, 1049–1054. [Google Scholar] [CrossRef] [Green Version]
  4. Hosking, J.R.M. The Theory of Probability Weighted Moments; IBM Research Center Yorktown: Heights, NY, USA, 1986. [Google Scholar]
  5. Hosking, J.R.M. FORTRAN Routines for Use with the Method of L-moments; IBM Research Center: Yorktown Height, NY, USA, 1996; pp. 1–43. [Google Scholar]
  6. Maeng, S.J.; Lee, S.H. Estimation of Design Flood by Determination of Best Fitting Order of LH-Moments(II). J. Korean Soc. Agric. Eng. 2003, 45, 33–44. [Google Scholar]
  7. Maidment, D.R. Handbook of Hydrology, 1st ed.McGraw-Hill Inc.: New York, NY, USA, 1992; pp. 1–55. [Google Scholar]
  8. Cunnane, C. Statistical Distributions for Flood Frequency Analysis; World Meteorological Organization (WMO): Geneva, Switzerland, 1989; pp. A4.1–14. [Google Scholar]
  9. Hosking, J.R.M.; Wallis, J.R. Regional Frequency Analysis; Cambridge University Press: Cambridge, UK, 1997; p. 240. [Google Scholar]
  10. Environmental Protection Agency (US-EPA). Protecting Water Quality Urban Runoff; US-EPA: Washington, DC, USA, 2015. Available online: https://www3.epa.gov/npdes/pubs/nps_urban-facts_final.pdf (accessed on 6 May 2019).
  11. Lee, D.C. Runoff Reduction by LID of A Newtown Considering Soil Characteristic. Doctoral Dissertation, Chung-Ang University, Seoul, Korea, 2017. [Google Scholar]
  12. Korea Meteorological Administration. Meteorological Data Open Portal. Available online: https://data.kma.go.kr (accessed on 20 February 2019).
  13. Korea Meteorological Administration. Climate Information Portal. Available online: http://www.climate.go.kr (accessed on 7 March 2019).
  14. Wald, A.; Wolfowitz, J. On a test whether two samples are from the same population. Ann. Math. Stat. 1940, 11, 147–162. [Google Scholar] [CrossRef]
  15. National Emergency Management Agency (NEMA), South Korea. Guideline for Preliminary Disaster Inspection; National Emergency Management Agency: Seoul, Korea, 2012. (In Korean)
Figure 1. Study area.
Figure 1. Study area.
Water 14 01301 g001
Figure 2. Selection of rainfall station.
Figure 2. Selection of rainfall station.
Water 14 01301 g002
Figure 3. Production process for climate change scenarios [13].
Figure 3. Production process for climate change scenarios [13].
Water 14 01301 g003
Figure 4. Pipe network of watersheds. (A) Outlet of section A; (B) Outlet of section B; (C) Outlet of section C.
Figure 4. Pipe network of watersheds. (A) Outlet of section A; (B) Outlet of section B; (C) Outlet of section C.
Water 14 01301 g004
Figure 5. Estimation of probability precipitation by climate change scenario.
Figure 5. Estimation of probability precipitation by climate change scenario.
Water 14 01301 g005
Figure 6. Study flow chart.
Figure 6. Study flow chart.
Water 14 01301 g006
Figure 7. Identification of distributions by L-moment.
Figure 7. Identification of distributions by L-moment.
Water 14 01301 g007
Figure 8. Runoff by return period according to LID facility installation (RCP 4.5, 1988–2100): (a) 10-year return period; (b) 20-year return period; (c) 50-year return period; (d) 100-year return period; (e) 200-year return period.
Figure 8. Runoff by return period according to LID facility installation (RCP 4.5, 1988–2100): (a) 10-year return period; (b) 20-year return period; (c) 50-year return period; (d) 100-year return period; (e) 200-year return period.
Water 14 01301 g008aWater 14 01301 g008b
Figure 9. Runoff by return period according to LID facility installation (RCP 8.5, 1988–2100): (a) 10-year return period; (b) 20-year return period; (c) 50-year return period; (d) 100-year return period; (e) 200-year return period.
Figure 9. Runoff by return period according to LID facility installation (RCP 8.5, 1988–2100): (a) 10-year return period; (b) 20-year return period; (c) 50-year return period; (d) 100-year return period; (e) 200-year return period.
Water 14 01301 g009aWater 14 01301 g009bWater 14 01301 g009c
Figure 10. Change in runoff by return period in RCP 4.5 scenario: (a) 1988–2018 (year); (b) 1988–2050 (year); (c) 1988–2100 (year).
Figure 10. Change in runoff by return period in RCP 4.5 scenario: (a) 1988–2018 (year); (b) 1988–2050 (year); (c) 1988–2100 (year).
Water 14 01301 g010aWater 14 01301 g010b
Figure 11. Change in runoff by return period in RCP 8.5 scenario: (a) 1988–2050 (year); (b) 1988–2100 (year).
Figure 11. Change in runoff by return period in RCP 8.5 scenario: (a) 1988–2050 (year); (b) 1988–2100 (year).
Water 14 01301 g011
Table 1. Annual maximum daily rainfall according to climate change scenarios.
Table 1. Annual maximum daily rainfall according to climate change scenarios.
YearRCP ScenariosYearRCP ScenariosYearRCP Scenarios
RCP 4.5RCP 8.5RCP 4.5RCP 8.5RCP 4.5RCP 8.5
198883.083.02026145.0 205.8 2064145.687.7
1989296.0296.0202793.7 70.0 2065128.2106.6
1990126.0126.02028134.0 117.3 2066219.9109.7
1991218.0218.0202985.7 123.9 206798.765.7
1992164.0164.02030117.0 68.0 206891.5112.0
1993140.0140.02031113.0 212.1 206984.9101.5
199497.097.0203250.9 95.4 207065.0134.7
199593.093.0203366.3 256.1 2071116.591.9
1996103.0103.02034189.7 94.0 2072264.4120.2
1997119.0119.02035124.1 105.8 2073157.4187.8
1998109.0109.02036168.7 108.7 2074102.268.3
1999131.0131.0203757.5 124.8 207579.071.7
200083.083.0203882.4 83.4 207665.697.5
200173.073.02039229.4 76.1 207783.8163.5
2002160.0160.0204085.2 121.6 207886.580.9
2003108.0108.02041113.7 63.5 2079129.582.2
2004127.0127.02042133.6 71.9 2080108.6123.7
2005179.0179.0204380.1 121.1 208195.4128.6
2006131.0131.02044130.2 82.7 2082156.3127.3
2007100.0100.0204580.9 112.8 2083116.7171.2
2008111.0111.02046216.8 131.9 2084129.9101.4
2009218.0218.0204786.7 96.9 2085249.089.2
2010119.0119.02048156.3 82.0 208679.572.0
2011189.0189.0 204993.2 60.0 2087107.5122.0
2012119.0119.0 2050110.4 69.1 2088137.3403.6
2013102.0102.0 2051111.5 147.9 208979.983.3
2014238.0238.0 2052157.8 334.9 209075.297.9
201553.053.0 2053176.7 111.8 2091162.381.1
2016140.0140.0 2054114.2 113.3 2092112.9269.1
2017168.0168.0 2055140.5 162.6 209369.2142.0
201888.088.0 205699.1 63.4 2094150.4121.1
201988.9121.8 2057355.9 105.8 2095152.4124.4
2020121.181.9 2058108.3 101.7 2096123.0141.5
2021157.2151.2 2059189.6 69.5 209796.4181.9
2022131.685.5 206075.3 91.7 2098132.4151.2
2023133.297.4 2061152.4 110.9 2099100.3105.6
202484.785.2 206281.2 87.0 210065.774.6
2025127.673.9 2063103.8 104.6
Table 2. Pipe network.
Table 2. Pipe network.
Watershed Area
(ha)
Subwatershed
(pcs)
Outlet
(pcs)
Node (pcs)Link
(pcs)
A29.413314343
B32.417115252
C44.120219393
Total105.95063188188
Table 3. Parameter values of permeable pavement.
Table 3. Parameter values of permeable pavement.
DivisionAttribute Value
SurfaceBerm height (mm)0
Vegetation volume fraction0
Surface roughness (Manning n)0.012
Surface slope (%)1
PavementThickness (mm)60
Void ratio (Voids/Solids)0.15
Impervious surface fraction0.95
Permeability (mm/h)360
Clogging factor0
StorageHeight (mm)300
Void ratio (Voids/Solids)0.25
Conductivity (mm/h)25.9
Clogging factor0
UnderdrainDrain coefficient (mm/h)0
Drain exponent0
Drain offset height0
In addition, the drainage area of EDC is 52.95 ha, which is 50% of the total watershed area of 105.90 ha. The roughness coefficients among the conduit input data, 0.030 for the permeable area and 0.014 for the impervious area, were applied.
Table 4. Basic statistics for Gimhae rainfall station using RCP 4.5 scenario.
Table 4. Basic statistics for Gimhae rainfall station using RCP 4.5 scenario.
Duration of
Observed Data
(Years)
Mean ( x ¯ )
(m3/s)
Standard Deviation
(S)
Coefficient
of Skewness
( C s )
Coefficient
of Variation
( C v )
Coefficient
of Kurtosis
(CA)
1988–2018135.053.3611.2310.3951.616
1988–2050126.648.5021.1620.3831.646
1988–2100126.251.4281.6180.4083.769
Table 5. Basic statistics for Gimhae rainfall station using RCP 8.5 scenario.
Table 5. Basic statistics for Gimhae rainfall station using RCP 8.5 scenario.
Duration of
Observed Data
(Years)
Mean ( x ¯ )
(m3/s)
Standard Deviation
(S)
Coefficient
of Skewness
( C s )
Coefficient
of Variation
( C v )
Coefficient
of Kurtosis
(CA)
1988–2018135.0 53.3611.2310.3951.616
1988–2050121.2 50.7271.4210.4181.947
1988–2100122.4 56.4912.2550.4616.828
Table 6. Parameters calculated by GEV distribution using L-moment method considering RCP 4.5 scenario at Gimhae rainfall station.
Table 6. Parameters calculated by GEV distribution using L-moment method considering RCP 4.5 scenario at Gimhae rainfall station.
Duration of
Observed Data
(Years)
Parameter
ξ α β
1988–2018108.244296.9893−0.1260
1988–2050103.495535.6241−0.0668
1988–2100101.838834.1054−0.1224
Table 7. Parameters calculated by GEV distribution using L-moment method considering RCP 8.5 scenario at Gimhae rainfall station.
Table 7. Parameters calculated by GEV distribution using L-moment method considering RCP 8.5 scenario at Gimhae rainfall station.
Duration of
Observed Data
(Years)
Parameter
ξ α β
1988–2018108.424296.9893−0.1260
1988–205096.226832.2964−0.1675
1988–210095.866330.1751−0.2377
Table 8. Probable annual maximum daily rainfall at Gimhae rainfall station using RCP 4.5 scenario (unit: mm/day).
Table 8. Probable annual maximum daily rainfall at Gimhae rainfall station using RCP 4.5 scenario (unit: mm/day).
Duration of
Observed Data
(Years)
Return Period (Years)
102050100200
1988–2018204.66241.67294.83338.97386.98
1988–2050190.00220.54262.30295.35329.85
1988–2100190.20224.00272.41312.49355.98
Table 9. Probable annual maximum daily rainfall at Gimhae rainfall station using RCP 8.5 scenario (unit: mm/day).
Table 9. Probable annual maximum daily rainfall at Gimhae rainfall station using RCP 8.5 scenario (unit: mm/day).
Duration of
Observed Data
(Years)
Return Period (Years)
102050100200
1988–2018204.66241.67294.83338.97386.98
1988–2050184.50220.51274.07320.05371.53
1988–2100185.65226.10289.85347.79415.92
Table 10. Runoff and reduction rate by return period of maximum daily rainfall (RCP 4.5, 1988–2100).
Table 10. Runoff and reduction rate by return period of maximum daily rainfall (RCP 4.5, 1988–2100).
Return
Period
(Years)
Watershed
Area
Peak Rate of Runoff
(m3/s)
Reduction
Efficiency
(%)
Runoff
(m3/Day)
Reduction
Efficiency
(%)
BeforeAfterBeforeAfter
10A0.084 0.052 38.5 105.2 54.6 48.1
B0.079 0.050 37.1 101.4 57.4 43.4
C0.088 0.058 34.6 115.2 63.7 44.7
Total0.252 0.159 36.7 321.8 175.7 45.4
20A0.104 0.063 39.4 126.0 64.9 48.5
B0.097 0.061 37.7 122.2 68.7 43.8
C0.109 0.070 35.4 138.5 76.0 45.2
Total0.310 0.194 37.4 386.8 209.6 45.8
50A0.133 0.080 40.3 156.1 79.7 49.0
B0.124 0.077 38.3 152.6 85.1 44.2
C0.140 0.089 36.3 172.3 93.6 45.7
Total0.397 0.245 38.3 481.1 258.4 46.3
100A0.158 0.093 41.0 181.2 92.0 49.2
B0.147 0.090 38.7 178.2 98.8 44.5
C0.166 0.105 36.9 200.6 108.2 46.0
Total0.471 0.288 38.8 559.9 299.0 46.6
200A0.186 0.101 45.7 208.6 105.4 49.5
B0.172 0.105 39.0 206.2 113.7 44.8
C0.196 0.122 37.5 231.4 124.2 46.3
Total0.554 0.329 40.6 646.2 343.3 46.9
Table 11. Runoff and reduction rate by return period of maximum daily rainfall (RCP 8.5, 1988–2100).
Table 11. Runoff and reduction rate by return period of maximum daily rainfall (RCP 8.5, 1988–2100).
Return
Period
(Years)
Watershed
Area
Peak Rate of Runoff
(m3/s)
Reduction
Efficiency
(%)
Runoff
(m3/day)
Reduction
Efficiency
(%)
BeforeAfterBeforeAfter
10A0.082 0.05038.4102.453.248.1
B0.077 0.04837.098.655.943.3
C0.085 0.056 34.5112.162.144.6
Total0.244 0.155 36.6313.1171.245.3
20A0.105 0.064 39.4127.365.548.5
B0.099 0.061 37.7123.569.443.8
C0.110 0.071 35.4140.076.745.2
Total0.314 0.196 37.5390.8211.745.8
50A0.144 0.086 40.6167.085.149.1
B0.134 0.083 38.5163.791.044.4
C0.151 0.096 36.6184.699.945.9
Total0.429 0.264 38.5515.3276.046.4
100A0.181 0.106 41.4203.4102.949.4
B0.168 0.102 39.0200.9110.944.8
C0.190 0.119 37.4225.6121.246.3
Total0.538 0.327 39.3629.9334.946.8
200A0.225 0.130 42.2246.5123.849.8
B0.208 0.126 39.4245.3134.545.2
C0.238 0.147 38.2274.2146.246.7
Total0.671 0.403 39.9766.1404.547.2
Table 12. Runoff by return period using RCP 4.5 scenario (unit: m3/day).
Table 12. Runoff by return period using RCP 4.5 scenario (unit: m3/day).
Year of AnalysisDevelopmentReturn Period
10-Year20-Year50-Year100-Year200-Year
1988–2018Before349.5 421.1 525.1 612.4 708.1
After190.2 227.4 281.1 326.0 376.1
1988–2050Before321.4 380.1 461.3 526.1 594.3
After175.5 206.1 248.2 281.6 316.7
1988–2100Before321.8 386.8 481.1 559.9 646.2
After175.7 209.6 258.4 299.0 343.3
Table 13. Runoff by return period using RCP 8.5 scenario (unit: m3/day).
Table 13. Runoff by return period using RCP 8.5 scenario (unit: m3/day).
Duration of Analysis
(Year)
DevelopmentReturn Period
10-Year20-Year50-Year100-Year200-Year
1988–2050Before310.9 380.0 484.3 574.9 677.2
After170.0 206.1 260.1 306.7 359.1
1988–2100Before313.1 390.8 515.3 629.9 766.1
After171.2 211.7 276.0 334.9 404.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lee, S.; Kim, D.; Maeng, S.; Azam, M.; Lee, B. Runoff Reduction Effects at Installation of LID Facilities under Different Climate Change Scenarios. Water 2022, 14, 1301. https://doi.org/10.3390/w14081301

AMA Style

Lee S, Kim D, Maeng S, Azam M, Lee B. Runoff Reduction Effects at Installation of LID Facilities under Different Climate Change Scenarios. Water. 2022; 14(8):1301. https://doi.org/10.3390/w14081301

Chicago/Turabian Style

Lee, Seungwook, Daye Kim, Seungjin Maeng, Muhammad Azam, and Bongguk Lee. 2022. "Runoff Reduction Effects at Installation of LID Facilities under Different Climate Change Scenarios" Water 14, no. 8: 1301. https://doi.org/10.3390/w14081301

APA Style

Lee, S., Kim, D., Maeng, S., Azam, M., & Lee, B. (2022). Runoff Reduction Effects at Installation of LID Facilities under Different Climate Change Scenarios. Water, 14(8), 1301. https://doi.org/10.3390/w14081301

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop