**About the Editors**

#### **Akira Kawamura**

Akira Kawamura mandatorily retired from Tokyo Metropolitan University in 2021 after serving as a professor for 16 and half years. Afterward, he was awarded the title of Professor Emeritus. Akira Kawamura graduated with a Doctor of Engineering degree in Hydrology and Water Resources Engineering from Kyushu University in 1985. He was first appointed as a research associate in the Department of Civil Engineering Hydraulics at Kyushu University in 1985 and promoted to associate professor in 1993. He transferred to Tokyo Metropolitan University as a professor in 2004. His main fields of research interest include the real-time prediction of hydro-meteorological time series, urban rainfall–runoff models, integrated flood risk management, the analysis of chaotic hydrological phenomena, the sustainability assessment of water resources, X-band MP radar rainfall data applications, and deep learning models in hydrology, among others. To date, Prof. Kawamura has contributed to 202 peer-reviewed scientific publications, 48 of which are in international journals with high impact factors.

#### **Kei Nakagawa**

Kei Nakagawa is a professor of environmental groundwater science with 30 years of research experience. He was first appointed as an assistant professor in soil science at the Department of Agricultural Chemistry in 1999 at Kyushu University and then promoted to associate professor in water use engineering at the Department of Agricultural Engineering in 2002 at Kagoshima University. In April 2011, he was appointed as a full professor of the Graduate School of Fisheries and Environmental Sciences at Nagasaki University. His main fields of research interest include reactive transport in groundwater, physical and chemical hydrogeology and heterogeneity, saltwater intrusion, the performance evaluation of subsurface dams in the coastal aquifer, groundwater modeling, and the remediation of contaminated soils and groundwater (phytoremediation and electro-kinetic remediation of heavy metal contaminated soils). Since being appointed at Nagasaki University in 2011, his main research topic has been nitrate pollution resulting from agricultural activities.

### *Editorial* **Urban Water-Related Problems**

**Akira Kawamura 1,\* and Kei Nakagawa <sup>2</sup>**


#### **1. Introduction**

Urban areas are considered to be the most vulnerable to water-related problems, which involve a lack or excess of water problems from the perspectives of quantity and quality. These specific phenomena include flash floods and inundation, droughts and water shortages, surface and ground water pollution, tsunamis and storm surges, landslides and mudflows, the degradation of fluvial and aquatic ecosystems, and unsanitary conditions and epidemics, among others. In urban areas, water-related problems cause immense human losses and economic damage. Water-related problems frequently reoccur in urban areas, and are intricately linked with each other, posing major obstacles to the achievement of human security and the sustainable socio-economic development of cities. Thus, it is crucial that they are scientifically and comprehensively discussed, so that they can be better understood, in order to fight against and mitigate these problems.

The special session "Urban Water-Related Problems" has been held in the AOGS (Asia Oceania Geosciences Society) annual meeting since 2017. Thus far, we have discussed a very wide range of urban water-related problems. For this Special Issue published in *Water*, we intend to invite studies on these broad topics and encourage a collective perspective on urban water-related problems.

This Special Issue comprises four review papers and nine research articles, with contributions from 55 authors of six countries. All four of the review papers are related to the problems associated with flash flood phenomena, among the many broad topics mentioned above, but they focus on very different aspects of the phenomena, such as urban runoff modeling in Japan [1], real-time urban flood forecasting systems for Southeast Asia [2], frequency analysis of urban floods [3], and the reduction of non-point pollution as well as flood runoff by porous concrete infiltration [4].

On the other hand, among nine research articles, five articles are related to water quality problems, of which two articles investigate groundwater pollution in Japan and Sweden [5,6], one article deals with lake contamination by phosphorus in Sweden [7], another targets lakes or reservoirs from the perspective of the hydrodynamic response by wind [8], and the last one considers the water pollution problem for a riverbank area in Bangladesh [9]. Three research articles out of the nine are related to problems of urban flood phenomena, in which studies on urban flooding area characteristics in Japan [10], water level forecasting caused by urban floods in Japan [11], and barriers to IFRM (Integrated Flood Risk Management) adaptation in the Philippines [12] are carried out. The last article is related to water resource problems, investigating the urban domestic water supply system from the viewpoint of willingness to pay [13].

Unfortunately, studies dealing with the phenomena related to tsunamis and storm surges, landslides and mudflows, and the degradation of fluvial and aquatic ecosystems are not published in this Special Issue.

**Citation:** Kawamura, A.; Nakagawa, K. Urban Water-Related Problems. *Water* **2023**, *15*, 3280. https:// doi.org/10.3390/w15183280

Received: 28 August 2023 Accepted: 12 September 2023 Published: 16 September 2023

**Copyright:** © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

#### **2. Overview of the Contributions of This Special Issue**

The review paper "Urban Flood Runoff Modeling in Japan: Recent Developments and Future Prospects" [1] summarizes, discusses, and shares key outputs from some of the main research directions in urban flood runoff modeling, significant parts of which have been uniquely developed in Japan and only published in Japanese. In the paper, after a general introduction to urban runoff modeling, the authors present key historical works in Japan, followed by a description of the situation in Japan with respect to observations of precipitation and water level. Then, the storage function model approach is reviewed, including an extension to urban basins, as well as recent experiments with AI-based emulation in Japanese basins. Subsequently, the authors review the prospects of detailed hydrodynamic modeling involving high-resolution, vector-based Geographical Information System (GIS) data for the optimal description of the urban environment with applications in Tokyo. The authors conclude the paper with some future prospects related to urban flood risk modeling and assessment in Japan.

In the paper "Real-Time Urban Flood Forecasting Systems for Southeast Asia—A Review of Present Modelling and Its Future Prospects" [2], the authors reviewed stateof-the-art models of real-time forecasting systems in countries in Southeast Asia, such as Thailand, for urban flash floods. A real-time system basically consists of rainfall forecasting, drainage system modeling, and inundation area mapping. The authors summarized the recent radar data utilization methods for rainfall forecasting, physical-process-based hydraulic models for flood inundation prediction, and data-driven artificial intelligence (AI) models for the real-time forecasting systems. The authors also dealt with available technologies for modeling digital surface models (DSMs) for the finer urban terrain of drainage systems. The review indicated that an obstacle to using process-based hydraulic models was the limited computational resources and shorter lead time for real-time forecasting in many urban areas of tropical Southeast Asia.

The review paper "Frequency Analysis of Hydrological Data for Urban Floods— Review of Traditional Methods and Recent Developments, Especially an Introduction of Japanese Proper Methods" [3] investigates in detail the frequency analysis of hydrological data for urban floods, targeting the Japanese methods in particular. First, the authors introduced well-used Japanese frequency analysis methods, because some techniques that are slightly different from the international standard have been used in Japan for many years. This review emphasized discussions of the parameter estimation of stochastic models and the selection of optimal statistical models, which include the evaluation of goodnessof-fit techniques of statistical models. Based on these results, the authors criticized the Japanese standard procedures recommended by the central government, and indicated that consistency between parameter estimation and the evaluation of goodness-of-fit was necessary. From this perspective, the authors recommended using the maximum likelihood method and AIC. The authors also recommended, when using SLSC method, to apply not only SLSC itself but also SLSC's non-exceedance probability.

In the review paper "Application of Porous Concrete Infiltration Techniques to Street Stormwater Inlets That Simultaneously Mitigate against Non-Point Heavy Metal Pollution and Stormwater Runoff Reduction in Urban Areas: Catchment-Scale Evaluation of the Potential of Discrete and Small-Scale Techniques" [4], first, the literature related to the expansion of pervious areas was introduced, because the expansion of pervious areas was an essential and common concept in mitigating nonpoint pollution runoff in urban areas. Then, the potential application of porous concrete as a medium for constructing the bottom and side walls of street stormwater inlets was investigated. The effectiveness of the medium in reducing (i) the stormwater runoff volume and (ii) the heavy metal pollution runoff loads was assessed. In the paper, sensitivity tests were performed by changing the exfiltration rates. The results indicate that porous concrete used at only the bottom and side walls of the street stormwater inlets was suitable for reducing the runoff volume and removing any heavy metals from stormwater at a catchment scale.

In the paper "Groundwater Quality and Potential Pollution in the Southern Shimabara Peninsula, Japan" [5], the authors aimed to assess the groundwater quality in the southern regions of the Shimabara Peninsula in Nagasaki Prefecture, Japan, for a comprehensive understanding of the overall situation against the background that the peninsula is known for serious groundwater contamination caused by nitrate pollution. To achieve this goal, groundwater samples were collected at 56 locations in Minami-Shimabara City from 28 July to 4 August 2021. The authors showed using a Piper-trilinear diagram that almost 80% of the samples were classified as the alkaline earth carbonate type. From the Stiff diagrams, they showed that most of the locations exhibited a Ca-HCO3 water type, while the Na-HCO3 and Mg-HCO3 types were also observed in coastal areas. Finally, from the principal component analysis, the authors revealed that the first component corresponded to dissolved constituents and denitrification, and the second component indicated the effects of ion exchange and low nitrate pollution.

The paper "PFAS in the Drinking Water Source: Analysis of the Contamination Levels, Origin and Emission Rates" [6] delves into the prevailing groundwater pollution situation in southern Sweden against the background that waterborne pollution caused by polyfluoroalkyl substances (PFAS) has been reported in numerous countries globally, emerging as a novel concern that is capturing attention. The objective of this study was to enhance understanding regarding the historical contamination of water sources due to PFAS and to explore interconnected facets of PFAS-containing aqueous film-forming foam (PFAS-AFFF) application in fire training exercises and equipment trials. To achieve this goal, the authors encompassed the analysis of contamination extents and PFAS compositions, the evaluation of potential PFAS and PFAS-AFFF emission scenarios, and the assessment of contaminant transport conditions.

The paper "Decision Support for Lake Restoration: A Case Study in Swedish Freshwater Bodies" [7] presents techniques aimed at providing decision support for lake restoration through a combination of multi-criteria analysis and decision analysis. The main objective is to identify effective and efficient measures for restoring lakes. The study focuses on evaluating six commonly used lake restoration methods for reducing internal phosphorus loads in two selected lakes, using criteria such as cost, longevity, and effectiveness. The findings indicate that aluminum treatment emerges as the most favorable option due to its superior effectiveness and cost efficiency. The authors suggested that these methodologies could be utilized not only in specific countries but also in the restoration of lakes and ponds worldwide.

The paper "A Detailed Analysis on Hydrodynamic Response of a Highly Stratified Lake to Spatio-Temporally Varying Wind Field" [8] explores the hydrodynamic reactions of a highly stratified lake or reservoir to various inhomogeneous wind conditions through the application of numerical models and integrated analysis. This approach was taken because wind was commonly acknowledged as a vital factor driving transport and mixing processes in stratified, enclosed systems. The authors demonstrated that under the influence of nonuniform wind conditions, often characterized by weak-to-moderate wind with high spatial variances, the lake's hydrodynamic responses differed from those under uniform wind conditions. On the other hand, low spatial variances, often associated with strong wind, resulted in hydrodynamic responses identical to those of uniform wind. Thus, the authors recommended using non-uniform wind in model simulation for achieving realistic results, improving water management in lakes and reservoirs.

The paper "Spatial Distribution and Source Identification of Water Quality Parameters of an Industrial Seaport Riverbank Area in Bangladesh" [9] employed a range of water quality indices, including the Metal Index (MI), Comprehensive Pollution Index (CPI), and Weighted Arithmetic Water Quality Index (WQI), to enhance understanding of pollution dispersion and the underlying processes influencing river water quality. Additionally, multivariate statistical techniques were utilized to assess pollutant loads and sources within the Pasur River system in Bangladesh. The authors indicated that contaminant origins encompass both geogenic and anthropogenic factors, involving untreated or inadequately

treated wastewater from industries and the discharge of urban domestic waste. The effectiveness of the water quality assessment and pollution source identification methodologies introduced by the authors was demonstrated, with potential applicability on a global scale.

In the paper "Topographical Characteristics of Frequent Urban Pluvial Flooding Areas in Osaka and Nagoya Cities, Japan" [10], frequent urban pluvial flooding areas using 20 years of urban pluvial flooding area records during 1993–2012 were identified and analyzed using the principal component analysis of their topographical characteristics in Osaka and Nagoya Cities, Japan. The authors showed that the topographical characteristics of the frequent urban pluvial flooding areas in both cities were different. In Osaka City, not only the topographical characteristics, but also the influence of anthropogenic factors and stormwater drainage improvements were influential. On the other hand, in Nagoya City, the mere presence of dominant structures dammed up the inundated water and caused urban pluvial flooding. The authors quantitatively showed the paradigm shift of urban pluvial flooding factors from topographical characteristics to anthropogenic characteristics by the statistical analysis of newly defined urban pluvial flooding frequency areas.

The paper "Study on a Water-Level-Forecast Method Based on a Time Series Analysis of Urban River Basins—A Case Study of Shibuya River Basin in Tokyo" [11] investigated a vector autoregressive model to develop a water level forecast system that uses observed rainfall and water level. The model was targeted to ensure information conducive to evacuation approximately 20 min in advance without the need to build a physical model. The authors showed that the method based on time series analysis achieved a stable forecast, and indicated that the method can be applied as a water level forecast method for basins with an extremely fast flood arrival times and limited observation data.

In the paper "An Expanded Interpretive Structural Modeling Analysis of the Barriers to Integrated Flood Risk Management Adaptation in Metro Manila" [12], the authors proposed an expanded ISM (Ex-ISM) method to comprehensively analyze the interrelationships between the barriers to integrated flood risk management (IFRM) adaptation in Metro Manila, Philippines. The Ex-ISM enhanced conventional ISM in that the symbolism was modified to explicitly show the contextual interrelationships; then, the hierarchy assignment step was simplified, and a diagram was used to show all of the inter-relationships that allowed a comprehensive analysis. The authors showed that the results obtained using the Ex-ISM method did not deviate from those yielded by the conventional ISM method, but the Ex-ISM method allowed an easy assignment of hierarchy. The authors also showed that the Ex-ISM method was able to draw a diagram incorporating not only the direct but also the indirect interrelationships to provide a comprehensive analysis of the relationships between barriers.

The paper "Willingness to Pay for Improved Urban Domestic Water Supply System: The Case of Hanoi, Vietnam" [13] investigated water users' willingness to pay (WTP) for the improvement of Hanoi's domestic water supply. This was because in Hanoi, the capital of Vietnam, the municipal government is facing a number of difficulties in providing sufficient water in a sustainable manner due to not only the increasing urban population and the serious pollution of water resources, but also a lack of resources to invest in the supply system. In this paper, a contingent valuation process based on a survey of 402 respondents was used to explore citizens' willingness to pay for the improvement of their urban water supply. The authors revealed that Hanoi's urban communities were generally satisfied with the quantity of their water supply, but tended to be dissatisfied with its quality. The developed WTP regression model based on the survey findings showed that the average WTP was 1.4% of the average household income, taken as the affordability level of monthly water payments.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

**Akira Kawamura 1,\*, Hideo Amaguchi 1, Jonas Olsson <sup>2</sup> and Hiroto Tanouchi <sup>3</sup>**


**Abstract:** Since the 20th century, Japan has experienced a period of very rapid urbanization. Cities have experienced substantial densification and expansion, resulting in gradually elevated flood risk. Urban flooding has also occurred in most large cities in Japan, particularly in Tokyo. In response to this growing problem, much effort and resources have been spent on research and development aimed at understanding, simulating, and managing urban flood risk in Japan. The objective of this review is to summarize, discuss, and share key outputs from some of the main research directions in this field, significant parts of which have been uniquely developed in Japan and only published in Japanese. After a general introduction to urban runoff modeling, in the next section, key historical works in Japan are summarized, followed by a description of the situation in Japan with respect to observations of precipitation and water level. Then, the storage function model approach is reviewed, including an extension to urban basins, as well as recent experiments with AI-based emulation in Japanese basins. Subsequently, we review the prospects of detailed hydrodynamic modeling involving high-resolution, vector-based Geographical Information System (GIS) data for the optimal description of the urban environment with applications in Tokyo. We conclude the paper with some future prospects related to urban flood risk modeling and assessment in Japan.

**Citation:** Kawamura, A.; Amaguchi, H.; Olsson, J.; Tanouchi, H. Urban Flood Runoff Modeling in Japan: Recent Developments and Future Prospects. *Water* **2023**, *15*, 2733. https://doi.org/10.3390/ w15152733

Academic Editor: Fernando António Leal Pacheco

Received: 25 June 2023 Revised: 26 July 2023 Accepted: 26 July 2023 Published: 28 July 2023

**Copyright:** © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

**Keywords:** urban flood risk assessment; USF model; TSR model; urban landscape GIS delineation; rainfall data

#### **1. Introduction**

Development in urban areas greatly increases surface runoff, owing to the increase in impervious surfaces, as well as the total runoff volume flowing to the receiving watershed. In addition, the construction of storm sewer systems and the culverting of rivers during the urbanization process increase the runoff rate. The peak velocity of runoff inevitably increases as more runoff is discharged at shorter intervals, resulting in an increased risk of overflow. In conjunction with the concentration of population and property in the watershed, urbanization thus generally makes flood potential higher [1–4].

Since the latter half of the 20th century, urbanization has progressed rapidly in both developed and developing countries. Especially in Japan, urbanization has progressed much more rapidly than in Europe and the United States [5]. Mountains and hills in watershed areas that had been forests and fields were developed one after another, and agricultural land, including paddy fields, was simultaneously converted into residential land [5]. As a result, the watersheds were transformed into urban areas with a high risk of inundation. In an urban watershed, there is a complex distribution of artificially developed impervious areas, such as houses, buildings, paved roads, and parking lots, where rainwater cannot infiltrate. Runoff pathways by storm drainage facilities, such as

sewers and artificial drainage channels, as well as runoff control and flood control facilities (e.g., storage and infiltration facilities), form a very complex urban runoff system that is constantly changing [6]. In Japan, the 1958 Kano River Typhoon Event is considered the first case of urban flooding [7], which caused extensive inundation damage not only in the lowlands but also in previously flood-free plateau areas [7,8].

In terms of flood damage in Japan, there has been a decreasing trend in human losses, although the economic losses have not decreased. The damage per unit area has exhibited an upward trend since the 2000s. Recent floods in urban areas have become more complex than past urban floods. One reason is the increasing frequency of unpredictable localized torrential rains, commonly called "guerrilla downpours" [9–11], which might be considered to be an effect of urbanization [12]. Seino et al. [13] explored how urbanization affects precipitation in the Tokyo metropolitan area. The risk of pluvial flooding caused by surcharge from sewerage facilities has been shown to increase compared to fluvial flooding caused by overflow from river channels [14].

Urban flooding, especially associated with recent record-breaking rainfall events, has become one of the most pressing social issues. Thus, flood mitigation is a crucial water management strategy for reducing destructive damage. For this purpose, rainfallrunoff models are important tools that play a central role, especially in urban watersheds. Generally, the purpose of runoff analysis is to clarify the physical mechanism of the runoff process by applying the amount of precipitation (rainfall, snowfall, etc.) to predict hydrological impacts, such as floods and droughts, in the target watershed. The results are used for the planning and management of rivers and watersheds.

In terms of process description, runoff models can be classified into two main types: lumped and distributed models. The former considers the watershed as a homogeneous unit without spatial variability, and the runoff process is generally expressed as ordinary differential equations of time. Distributed models explicitly consider the spatial variability of watershed characteristics and are usually expressed as partial differential equations of time and place [15,16]. In addition, semi-distributed models also exist, where a watershed is divided into sub-watersheds, each of which is represented by a lumped model. Based on how the response of runoff as output to rainfall (precipitation) as input is described, models may be classified into "response models", which simply use a response relationship, "conceptual models", which conceptualize the response relationship, and "physical models", which model the runoff from the basic equations of physics.

The most commonly used lumped conceptual models in Japan today are the tank model [17] and the storage function (SF) model [18,19]. The tank model was developed by Sugawara in 1956 [17], and the SF model was proposed by Kimura in 1961 [18,20]. The SF model is a lumped model used to simulate individual flood events and has the ability to represent nonlinearity in the runoff process with a relatively simple structural equation and fast calculation. The SF model is commonly used in watersheds with relatively large catchment areas under the jurisdiction of the Ministry of Land, Infrastructure, and Transport (MLIT), and has become the standard runoff model for flood control planning in Japan [21–24].

In contrast to lumped models, distributed models explicitly incorporate the spatial variability of watershed characteristics. With the widespread use of Geographic Information System (GIS) technology, distributed physical models have been widely used in small-tomedium-sized urban watersheds, where the effects of land use change and flood control facility development on runoff characteristics have been evaluated [25–27]. The main advantages of commonly used grid-based distributed models are their simple structure and their use of spatially explicit watershed information, which is generally readily available. Because of these advantages, grid-based distributed models are widely used not only in rural areas but also in urban watersheds. In Japan, Sayama et al. [28] proposed the rainfall-runoff-inundation (RRI) model, which is a 2D model capable of simultaneously analyzing rainfall runoff and flood inundation. The RRI model has been integrated with modeling and visualization tools by the International Centre for Water Hazard and Risk

Management (ICHARM) of the Public Works Research Institute (PWRI) [29] and is used not only in Japan but also in other countries.

Globally, many different basin-scale hydrological models have been developed, including the TOPMODEL [30], Xinanjiang model [31,32], Systeme Hydrologique Europeen (SHE) [33], and the Institute of Hydrology Distributed model (IHDM) [34], and those are also commonly used in Japan. The Hydrologiska Byrån Vattenbalans (HBV) model was developed in Sweden in the 1970s and has been widely used worldwide [35]. More recently, the hydrological processes for the environment (HYPE) model [36,37] was developed, which has been set up from the basin to the global scale and is used for operational forecasting as well as climate change impact assessment. Tanouchi et al. [21] developed the semi-distributed HYPE model to improve the description of urban runoff.

The accurate characterization of impervious surface type and extent is required for the prediction of the corresponding impacts on watershed hydrology [38]. However, in general, the land use information used as watershed characteristics in grid-based models is simply the impervious area ratio and the runoff coefficient. Therefore, the ability to accurately represent the impermeable characteristics of urban watersheds in grid-based models is limited [39]. Many models use raster-type land use data generated from aerial photographs.

With recent technological advances in GIS, digital information has become readily available for urban watersheds. In addition to grid-shaped raster land use information, vector-type polygonal (polygon shape) land use data, which can accurately represent buildings, roads, and other land features, have also been created [40]. However, there are still only a few cases in which polygon-type land use information is directly used for runoff analysis. Sample et al. [41] used GIS to facilitate urban stormwater analysis using land use parcel boundaries. Rodriguez et al. [42] used a vector-based watershed description based on a so-called urban databank, which includes the categories cadastral parcel, building, street, sewer system, and river, to calculate an urban unit hydrograph. Lee et al. [43] developed a model to evaluate low-impact development (LID) by defining a watershed as a hydrological response unit (HRE), which is a sub-watershed for calculating stormwater inflow to the sewer system. However, the area of the HRE is very big at several thousand square meters, and land use information is averaged over the HRE.

In light of this background, we review flood runoff analysis and modeling with a focus on urban watersheds in Japan. Regarding previous research related to urban flood runoff modeling in Japan, this has already been reviewed in detail in the authorized Handbook [44] and Hydroscience Formulary [16,45–47] published by the Japan Society of Hydrology and Water Resources (JSHWR) and the Japan Society of Civil Engineers (JSCE), respectively. Therefore, in this review paper, we mainly focus on recent developments in urban flood runoff modeling. Following this introduction, in Section 2, we briefly describe the history of urban runoff modeling in Japan. Section 3 is devoted to observations of precipitation and water levels, which are required for accurate modeling and the current situation in Japan. Section 4 focuses on lumped modeling, the urban storage function (USF) model is outlined as an example, and the hydrological data required to build the model for urban areas are described together with the prediction accuracy obtained using artificial intelligence (AI). Distributed physical modeling is reviewed in Section 5, including different approaches to describing the urban structure, as well as the Tokyo storm runoff (TSR) model as an example of a distributed model using vector-type GIS data. This paper closes with some future prospects related to urban flood runoff modeling in Japan.

#### **2. Brief History on Urban Runoff Modeling in Japan**

In Japan, Morita [48] characterized an urban watershed as a so-called "small- to medium-sized urban river basin", where most of the watershed is flat land, such as alluvial plains, flood plains, alluvial fans, and hills. Normally, the entire upper, middle, and lower reaches are developed into urban areas, and the watershed area is generally smaller than 100 km2.

As for urban runoff research in Japan, Kinoshita [1] was the first to comprehensively and quantitatively describe and discuss urban flood runoff mechanisms in 1967. This paper was followed by numerous studies on urban flood runoff mechanisms [2,3], which were summarized by Kadoya [4] in 1985. There, the urban flood runoff mechanisms are characterized by (1) an increased runoff rate during floods, (2) a reduced flood concentration time before rainwater flows into rivers, and (3) extremely large peak runoff rates due to highintensity rainfall peaks with a short duration. Generally, most of the flood runoff models were mainly intended for mountainous and natural watersheds and were not developed with urban watersheds in mind, as pointed out in the Hydrology and Water Resources Handbook published by JSHWR in 1997 [49]. In addition, physical runoff processes in urban watersheds are not fundamentally different from those in natural watersheds. As a result, when applying a runoff model to an urban watershed, only the model parameters are usually adjusted [49].

The rational equation [50], which is used in many small-to-medium-sized river and sewage system plans in Japan, was originally proposed as a model for estimating peak flows as a basis for stormwater management in urban watersheds. Parameters, such as the runoff coefficient and flood concentration time, need to be calibrated for the target urban watershed. For example, in 1976, Kadoya and Fukushima [51] proposed an equation relating flood concentration time to flood magnitude based on the equivalent roughness method. The synthesized rational formula [52,53] has also been developed as a model for urban watersheds in 1998. This method creates a hydrograph by dividing an urban watershed into sub-watersheds, each characterized by area, runoff coefficient, and flood concentration time, and combining the peak flows by superimposing the results from the rational equation. The quasi-linear storage model [54] developed in 1977 and the modified RRL method [55] developed in 1972 are also used as flood runoff models for urban areas by setting flood concentration times and runoff coefficients for the urban watershed of interest [56,57].

In 1954, Iwagaki and Sueishi applied the kinematic wave theory to the analysis of hillslope runoff [58,59], and soon after, Lighthill and Whitham [60] developed the kinematic wave theory for river flow tracking. Since then, this theory has been accepted as a standard tool for modeling surface flows and various other elementary hydrological processes [61,62]. The kinematic wave model, in which the watershed is described in terms of a collection of rectangular slopes and river channels, and rainwater flow, which is represented by a hydraulic continuity equation and a kinematic equation, is a widely used distributed physical model in Japan (see, for example, [24,62]). As an example of applications of the kinematic wave model to urban watersheds, in 1991, Ando et al. [63] divided an urban watershed into rectangular sub-watersheds using spatially distributed rainfall and applied the kinematic wave model to each divided sub-watershed. The reproducibility of the hydrograph was improved by using the equivalent roughness value of the urban watershed as an input. In 1995, Fujimura et al. [64] showed that diffusion wave or dynamic wave models are more accurate than kinematic wave models for nonstationary and unsteady river channel flow calculations, although kinematic wave models are used for the slope flow, as well as surface flow, of low-lying watersheds. In 1996, Suzuki et al. [65] developed the Doken model as a distributed model for real-time flood forecasting. The model has been incorporated into the Integrated Flood Analysis System (IFAS) for flood forecasting in watersheds with inadequate water level observations and has been applied in Asian countries, such as Pakistan [66]. In 1997, Chikamori et al. [67] attempted to obtain model parameters, such as equivalent roughness, using GIS information when constructing a kinematic wave model for an urban watershed. In the Technical Standard for River Erosion Control, Survey Edition by MILT [68], the rational equation [50], the SF model [18], the kinematic wave model [24], the quasi-linear storage model [54], the tank model [17], and the Doken model [65] are introduced as runoff models.

Since then, along with the rapid development of computer and IT technology, the analysis of flood runoff in urban watersheds has been improved in several respects. Data collection by, e.g., satellites and radar is a significant advancement, as is the advancement of GIS technology [69]. Furthermore, a range of new methods based on, e.g., optimization methods, network theory, and AI, have been applied to flood runoff analysis and modeling. Shibuo and Furumai [70] recently reviewed urban stormwater management in Japan, focusing on rainfall observations, numerical modeling, and emerging technologies.

#### **3. Rainfall and Water Level Observation System in Japan**

A key component of flood runoff modeling is access to observational data. In Japan, information on rainfall and river levels is mainly provided by the Water Information and National Land Data Management Centre established by the Water Management and Land Conservation Bureau under the MLIT [71]. As of 2023, real-time river observations (water level and streamflow) and rainfall observations are provided online at 7580 and 10,619 stations, respectively [72]. In particular, the recently developed crisis water level gauges have improved river level monitoring during heavy rainfall, with approximately 6000 stations installed from 2018 to 2020. For rainfall, in addition to these, approximately 1300 precipitation monitoring stations have been installed by the Japan Meteorological Agency (JMA) through the Automated Meteorological Data Acquisition System (AMeDAS), with one station per 289 km2 across the country [73].

The X-band multi-parameter (MP) network of weather radars, known as the eXtended RAdar Information Network (XRAIN), was developed in 2008 by MLIT. The network uses 16 high-performance C-band MP radars and 39 X-band MP radars that are specialized for rainfall observation [74]. XRAIN provides high-resolution (250 m), high-frequency (1 min), and near-real-time rainfall information for almost the entire land area of Japan without the need for correction by rain gauges at the ground surface. However, owing to its short wavelength, XRAIN has only a narrow observation area (80 km radius) because of attenuation caused by rainfall, and in some cases, it cannot fully penetrate heavy rainfall areas.

In urban areas, small-to-medium-sized rivers are at high risk of river overflow and flood damage due to flooding caused by short periods of high-intensity rainfall. Therefore, urban runoff analysis requires very high-resolution observations to accurately capture the spatial and temporal distributions of rainfall. Schilling [75] suggested that rainfall data of at least 1–5 min and 1 km resolutions should be used for urban drainage modeling. Fabry et.al. [76] recommended the use of higher-resolution data (1–5 min temporal and 100–500 m spatial) for urban hydrological applications to provide the necessary hydraulic details. In the case of localized torrential rainfall, there is often a delay of only tens of minutes from the start of the rainfall to the peak river water level, and flooding damage associated with such events is frequent in large cities in Japan [77]. For this reason, the validation of flood runoff analysis models for urban areas requires a time resolution of river observations similar to that of the rainfall observations.

In Tokyo, for example, the Tokyo Metropolitan Flood Control Integrated Information System (FCIS) plays a role in flood prevention activities by automatically collecting realtime observations of rainfall and river water levels in Tokyo and providing them to the relevant organizations [78]. The FCIS began observing 10 min rainfall in January 1978 and later increased the resolution to 1 min. Water level observations were initiated to inform residents of the danger of flooding due to the rapid increase in water levels caused by heavy rainfall. A total of 124 water level observation stations were established in Tokyo, and 15 stations (1.67 km2/station) were established in the Kanda River basin. Rainfall observations were introduced to measure high-intensity rainfall with high resolution in time and space, as well as high precision. The number of rainfall observation stations has increased annually, and as of 2022, 93 stations (19 km2/station) have been installed (Figure 1). This is a very large number, considering that there are only 10 AMeDAS stations operated by the JMA [73] in Tokyo (178 km2/station). FCIS data are expected to accurately describe the spatial variability of heavy rainfall.

**Figure 1.** Map of FCIS and AMeDAS stations in Tokyo Metropolis (top) and a zoomed-in map of FCIS stations and XRAIN grid in the upper Kanda River Watershed (bottom).

#### **4. Storage Function Model in Urban Watersheds**

#### *4.1. Conventional SF Models*

SF models are event-based lumped models that are used for short-term hydrograph prediction, characterized by a relationship between storage and discharge. Kimura's original SF model in 1961 is a method to divide effective rainfall into storage and runoff processes by separating the watershed into runoff and infiltration areas [24] and is described as the "total storage function method" in the 1985 edition of the Hydroscience Formulary [45]. In general, the relationship between storage and discharge in a watershed during a flood can be described by a two-valued function, where the relationship between storage and discharge differs between the rising and receding parts of the hydrograph, respectively, and often forms a loop [79,80]. Kimura's SF model with two parameters introduced a delay time to eliminate this loop. In 1967, Prasad [81] presented a three-parameter SF model that directly expressed the bivalence between storage and discharge by adding a term for the loop effect. In 1982, Hoshi and Yamaoka [82] developed a four-parameter SF model by adding another parameter to accurately express the bivalence of the loop effect.

Flood runoff models, including conventional SF models, generally target the direct runoff component of the total runoff, except for tank models that incorporate rainwater loss and storage effects. Thus, general flood runoff models require effective rainfall as their input for the prediction of direct runoff. Hence, at first, the separation of baseflow from the total discharge is needed to calculate the effective rainfall when applied to the actual runoff analysis. There is no definitive method for separating direct runoff, and the calculated direct runoff will be different because of the subjective nature of the separation process [83]. As a result, the identification of model parameters may change significantly, which may cause the stability and reliability of the runoff calculations to decrease. Baba et al. [22] and Hoshi et al. [84] developed an SF model that directly uses the observed rainfall and runoff to solve the problems associated with the estimation of direct runoff. The proposal of a model concept that does not require the separation of baseflow is groundbreaking [69].

The four conventional SF models are the Linear, Kimura [18], Prasad [81], and Hoshi [84] models, and their storage equations are shown in Table 1, where *s* is the storage (mm), *Q* is the observed river discharge (mm/min), *t* is the time (min), *k*1, *k*<sup>2</sup> are model parameters that express the constant of proportionality for *Q* and the time variation of *Q*, respectively, and *p*1, *p*<sup>2</sup> are model parameters that express the exponent related to *Q* itself and the time variation of *Q*, respectively. The associated continuity equations, in which *ql* and *qR* include parameters *k*3, *z,* and *α*, respectively, are also listed, and they are explained in the next Section 4.2. Among these models, Hoshi's model was found to be superior owing to an additional parameter, *p*2, which may be quantified by numerical experiments and can well define the flow characteristics based on kinematic wave theory [85]. Some simplifications in Hoshi's storage model can lead to Prasad's storage model [81]. If we set *p*<sup>2</sup> = 1 in Hoshi's model, we obtain Prasad's storage model. Similarly, if we set *k*<sup>2</sup> = 0 in Prasad's model, the model can be transformed into Kimura's model [18]. The simplest linear model can be obtained by keeping *p*<sup>1</sup> = 1 in Kimura's model. Kimura's SF model with one storage tank, which is widely used as a special case of Kimura's original model with a delay time of zero [85].


**Table 1.** Storage function models with their associated continuity equations (PAR represents a parameter).

#### *4.2. Urban Storage Function (USF) Model*

The runoff characteristics of urban watersheds are significantly different from those of natural and mountainous watersheds because of the high proportion of impervious areas and the increasing prevalence of sewerage systems. Therefore, it is difficult to accurately reproduce the hydrographs of small-to-medium-sized urban rivers with short flood concentration times by using conventional SF models. Takasaki et al. [86] proposed an SF model including a two-valued function that considered the runoff characteristics of combined sewer systems. Later, Takasaki et al. [23] proposed and constructed an urban storage function (USF) model that does not require the calculation of effective rainfall or the separation of baseflow but directly uses observed rainfall and runoff.

Figure 2 shows the schematic diagram of all the possible inflow and outflow components of an urban watershed with a combined sewer system. Combined sewers are common in small-to-medium-sized river watersheds in old cities, such as the Tokyo Metropolitan District, where sewers have been in place since the 1900s [87]. The inflow components in Figure 2 are represented by rainfall, *R* (mm/min), and urban-specific and groundwater inflows from other watersheds, *I* (mm/min). Urban-specific inflows include leakage from water distribution pipes, irrigation flows, etc. The outflow components consist of river discharge, *Q* (mm/min), evapotranspiration *E* (mm/min), storm drainage from the watershed through the combined sewer system *qR* (mm/min), water extraction from the watershed for, e.g., water supply, agricultural needs, etc., *O* (mm/min), and groundwater-related loss *ql* (mm/min). In addition, domestic sewage *qw* is depicted in Figure 2 even though it does not contribute to watershed storage, *s* (mm).

**Figure 2.** Schematic diagram of all inflow and outflow components of an urban watershed with a combined sewer system.

The USF model is an empirical representation of Hoshi's SF model (shown in Table 1), in which the river discharge, *Q*, is replaced by the discharge, including storm drainage, *Q* + *qR*. The associated continuity equations, which do not require the separation of baseflow for both the conventional SF models and the USF model, are listed in Table 1. Here, the groundwater-related loss (*ql*) was defined by considering the infiltration hole height (*z*) by multiplying *k*3, which are the new parameters in the continuity equations in Table 1. Furthermore, the *qR* of the USF is defined by using an additional parameter, α, that expresses the slope of the linear relationship between the total discharge, *Q* + *qR*, and the drainage, *qR* (refer to [23,88] for details).

Thus, seven unknown parameters (*k*1, *k*2, *k*3, *p*1, *p*2, *z,* and *α*) are identified in the USF model. As an optimization method for parameter identification in lumped conceptual models, global search methods using evolutionary computation methods, such as genetic algorithms, have diversified in recent years [89]. Among them, the effectiveness of the shuffled complex evolution University of Arizona (SCE-UA) method developed by Duan et al. [90] has attracted considerable attention [91,92].

#### *4.3. Application of the USF Model in an Urban Watershed*

Takasaki et al. [23] applied the USF model to the upper Kanda River watershed (Figure 1), a typical small-to-medium-sized urban watershed in Tokyo that is entirely urbanized and where urban flooding occurs frequently. The small concentration time implies that river runoff occurs within a short period of time immediately after rainfall, and the use of hydrological data with very short time intervals for rainfall-runoff analysis is required. Therefore, FCIS rainfall and water level data at 1-minute intervals (described in Section 3) were used. They used the SCE-UA method for parameter identification and verified the effectiveness of the USF model, including the much-improved reproducibility of the hydrograph [23].

Padiyedath et al. [88] applied the five SF models shown in Table 1 to the upper Kanda River watershed. They also used the SCE-UA method to identify parameters to verify the performance of the SF models. As for the results, the USF model showed not only the highest hydrograph reproducibility among the five SF models, which is not surprising when considering that the hydrograph reproducibility increases with the number of model parameters, but also the most effective SF model from an Akaike information criterion (AIC) [93] perspective. In addition, the uncertainty of the seven parameters of the USF model was evaluated in detail using a residual-based bootstrap approach [94]. Furthermore, Padiyedath et al. [95] expanded the USF model by considering the spatial distribution of rainfall over watersheds and developed a generalized SF model able to directly estimate the water level [96].

Lumped models, such as the USF model with a relatively small number of model parameters, are suitable for real-time flood runoff forecasting by changing the model parameters sequentially. Parameter identification can be almost automated using global search methods like SCE-UA, and hydrograph reproducibility at the target point is generally higher than that of a distributed model. Takasaki et al. [97,98] applied both the particle filter method and the Kalman filter method to real-time flood runoff forecasting using the USF model in the upper Kanda River watershed and obtained accurate forecasting results.

#### *4.4. Runoff Prediction by USF Using XRAIN Data in an Urban Watershed*

For localized torrential rainfall events in the upper Kanda River watershed (Figure 1), the USF model has also been used to evaluate hydrograph reproducibility using XRAIN rainfall data and the FCIS rain gauge data [99]. At first, 1-minute XRAIN mesh rainfall data were compared to 1-minute FCIS rain gauge data of the same mesh to check the precision of XRAIN [100,101]. The results showed that 1-minute XRAIN values corresponded well with the 1-minute ground observation rainfall values, but the former were observed 1 to 3 minutes earlier than the ground observation, which was considered mainly due to the falling time of the raindrops. Spatiotemporal correlation analysis between XRAIN and FCIS rain gauge data has also been conducted [102,103]. As for the results, for rainfall during typhoons, where the rainfall area moved clearly along with the typhoon, it was shown that radar rainfall was highly correlated with ground rainfall with a 2-minute time delay. For localized torrential rainfall, it has been shown that the spatiotemporal correlation between XRAIN and ground rainfall is much more complicated than in a typhoon case, and the peak correlation does not always occur in the same mesh as the rain gauge but in the surrounding meshes.

Regarding runoff predictability, the best performance was found for FCIS rainfall observations, whereas the performance was lower for XRAIN [99]. This study also showed that it was possible to improve runoff prediction using XRAIN by considering spatiotemporal correlation characteristics. It should be emphasized that there is a spatial uncertainty in XRAIN rainfall, as it is measured at some elevation in the atmosphere and not at ground level. This uncertainty can significantly affect runoff prediction using XRAIN, particularly in the case of localized torrential rainfall occurring near the watershed boundary in small-to-medium-sized watersheds [99,103].

#### *4.5. AI-Based Flood Runoff Modeling in Urban Watersheds*

Machine learning using AI technology has developed remarkably in recent years. In order to utilize machine learning for flood forecasting, many AI-based flood runoff models have been constructed using, e.g., multilayered artificial neural networks (ANN) or deep learning (DL) [69]. Many of these models have been constructed for large river watersheds [104,105] using hourly observation data as input, and some DL applications for small-to-medium-sized mountain river watersheds have used 10-minute data to construct the models [106]. On the other hand, in small-to-medium-sized urban rivers, such as in Tokyo, water level and rainfall observations are carried out at 1-minute intervals, as described in Section 3. It has been concluded that the combination of temporal and volume resolutions in rainfall observations is critical for successful AI modeling. In particular, the hyetographs of 1-minute observations are highly intermittent and differ from the 10-minute and hourly data in large river watersheds.

An attempt was made to emulate the USF model using ANN and DL. Fujizuka et al. [107] tested the performance of a DL model using a virtual hyetograph and a virtual hydrograph and examined the impact of the amount of training data. The DL model was compared with a conventional ANN model on virtual datasets [108]. For both the ANN and DL models, if the number of nodes in the input layer is sufficient, the Nash–Sutcliffe coefficients of the hydrograph of both models are greater than 0.98; therefore, even the conventional ANN model can accurately emulate the USF model. However, when the total number of nodes in the intermediate layer was the same, the accuracy of the DL model was higher than that of the ANN model. Fujizuka et al. [109] also performed a comparison using actual observations and concluded that the DL model was superior when considering the number of parameters in relation to the number of observations.

Recently, hybrid forecasting models that combine an ANN forecasting model with a flood runoff model have been proposed [110]. In conventional ANN models, water level and flow rate are used as the main input data for training and forecasting. In contrast, the hybrid forecasting model incorporates flow rate forecasts from the runoff model into the input layer of the deep ANN model. This is expected to improve the forecasting of unprecedented runoff peaks, which remains a challenge for conventional statistical models.

#### **5. Distributed Physical Models in Urban Watersheds**

#### *5.1. Background*

In recent years, most of the distributed flood runoff models developed in Japan for urban watersheds have focused on surface flooding and inundation analysis, and few have focused on runoff analysis and the generation of runoff volume directly from urban impervious areas. Since 1980, with significant improvements in numerical techniques and hydrological models for slope flow, surface flow hydrological models have become available for predicting inundation, and inundation hazard maps have been prepared as flood damage reduction countermeasures. Most inundation analysis models are based on shallow water equations as applied, e.g., by Iwasa et al. [111]. Toda et al. [112] used a kinematic wave model for runoff from mountainous areas near urban areas and combined it with an inundation model for urban areas that also consider drainage by sewage systems. Subsequently, a new integrated lowland inundation model was developed and used for flood inundation forecasting, considering slope runoff processes, sewer flows, surface inundation, and runoff from manholes. This model configuration has been established as a standard model for flood damage prediction by flood inundation analysis in urban watersheds.

Furthermore, hydraulic models, such as InfoWorks, MOUSE, and the stormwater management model (SWMM), which are mainly intended for the design and flow analysis of sewer pipelines, have been integrated with surface inundation models and are used to analyze flooding due to poor stormwater pipeline drainage, etc. In Japan, although these models are used in practical applications, they are not widely used for flood runoff analysis in urban watersheds because, as mentioned above, the development of flooding analysis models that take sewer pipeline flow into account has been conducted independently. As an example of an integrated sewerage and inundation model, Shibuo et al. [113] developed a comprehensive analysis model for the Tsurumi River watershed using the water energybased distributed hydrological model (WEB-DHM) to calculate runoff from the upper catchment and InfoWorks to analyze the sewers and combined them into an integrated river channel and sewage model. In addition, many inundation analyses of urban watersheds have been conducted overseas using the Hydrologic Engineering Center's river analysis system (HEC-RAS), an open source 2D hydraulic model integrated into GIS [114,115]; this model is not widely used in Japan because there is a need to incorporate aspects such as flow in sewer pipelines and/or flooding processes in subways into the modeling. Detailed explanations regarding these aspects will be provided in the next section.

#### *5.2. Representation of Urban Structures in Distributed Physical Models*

Urban flooding processes can be described by a dual drainage model that combines a 1D hydrodynamic model of the stormwater drainage network model with a 2D hydrodynamic model of overland flow and flooding [116]. In particular, different representations of urban structures have significant effects on surface flood processes [117]. Table 2 summarizes the computational grid shapes used in different studies. The TSR model proposed by the authors modeled direct runoff generation and flood flows from the perspective of urban formation, whereas the other four models were developed to analyze urban flood inundation.


**Table 2.** Characteristics of a distributed urban flood runoff models.


#### **Table 2.** *Cont.*

#### *5.3. Structured Grid*

A structured grid is most commonly used for surface models because of the simplicity of grid formation when analyzing ground surface flows as two-dimensional unsteady flows. The surface flow (inundation analysis) model using a structured grid was developed based on a relatively large computational grid for wide-area inundation. At the beginning of the model development, the effects of structures, road networks, etc., were considered by equivalent roughness coefficients based on land use or synthetic roughness coefficients as a function of building occupancy and inundation depth. This was carried out because it was difficult to take into account the influence of roads and buildings, which are likely to have a significant impact on flood flows in urban areas, in the structured grid. Later, building volume (nonflooded space) and flow disturbance by buildings were considered [118].

Currently, this is the standard method used to prepare municipal flood maps in Japan, and the National Institute for Land and Infrastructure Management (NILIM) of the MLIT provides the analysis program NILIM2.0 [119]. The accuracy of inundation analysis can be improved by using high-resolution elevation data. In recent years, with advances in computer performance, analyses using 5 and 1 m mesh elevation data have been made. Inokawa and Kobayashi carried out an inundation analysis on the surface and in an underground station usinga1m mesh and concluded that it is important to consider slight changes in the inundation level at the station [120].

#### *5.4. Curvilinear Grid*

Curvilinear grids are widely used in the numerical analysis of open channel flows, including bends [121]. When this grid model is applied to an urban area, the attributes of roads and buildings can be classified for each grid if the co-ordinate axes are appropriately set. When compared to a structured grid, grid construction requires much more effort. Since there are only two axes, the direction of the road is restricted [121].

Fukuoka et al. focused on structures, such as houses and road networks, in urban areas and developed a model using a curvilinear grid for the road network, achieving the highly accurate reproduction of water levels and flow velocities on roads and between city blocks of houses [122]. However, in recent studies, this co-ordinate system has rarely been used.

#### *5.5. Unstructured Grid*

Unstructured grids can be created according to the urban structure, for example, by applying detailed grids in regions with complex boundaries and geometries and coarsening the grids in the regions that are not computationally important. However, as with a curvilinear grid, the time and effort required for grid formation are significant for unstructured grids [121].

Akiyama et al. developed a two-dimensional numerical model of flood flows based on the finite volume method (FVM) with an unstructured grid for the floodplain profile to ensure the proper representation of the road network and houses [123,124]. The flood flow in the river and the runoff on the ground were simulated simultaneously using a 2D free surface flow model based on the flux-difference splitting scheme and an unstructured FVM. For the sewer network, a dynamic network model for free-surface pressure flow based on the flux-difference splitting scheme combined with a stormwater drainage network was used. By performing an inundation analysis on Iizuka City and verifying its accuracy against observations of the waterlogged areas, it was shown that the model could estimate waterlogged areas with sufficient accuracy and estimate detailed inundation properties, such as overflow, internal inundation locations, and flood propagation paths.

#### *5.6. Road Network*

Inoue et al. proposed representing roads and buildings as different computational grids and developed an analysis method using a road network model that can represent flooding in real urban areas more faithfully than the method using a structured grid [125]. Road junctions were used as nodes, and the road segments between junctions as links; residential areas were specified as areas surrounded by the network; junctions and residential areas were defined as the analysis grid; and water depths were obtained from the continuity equation. They showed that the road network model could directly parameterize the roughness coefficients at the surface by considering the block walls.

Sekine et al. [126] proposed another flood analysis model using a road network, focusing on flow concentrations on roads and sewers in densely built-up areas. The model is characterized by its ability to faithfully model rainwater runoff from buildings into sewers. The average house in a city block, based on its roof area, was treated as a unit house in the model, and rainwater on the unit house was assumed to flow into the sewers adjacent to the city block using the kinematic wave approximation. Sekine et al. focused on three cities neighboring the 23 wards of Tokyo that were damaged by Typhoon No. 19 in 2019 and performed predictive calculations of inundation in the event of hypothetical maximum rainfall, describing the inundation process in depressed areas and underpasses with significant inundation [127].

#### *5.7. Urban Landscape GIS Delineation*

The authors have developed the Tokyo storm runoff (TSR) model that can simulate urban storm runoff and flood inundation using a vector-based watershed description [6]. The TSR model employs so-called "urban landscape GIS delineation" that realistically describes complicated urban land use features in detail. Urban landscape GIS delineation enables an accurate representation of the infiltration and perviousness of each property, making it possible to calculate the impervious area of a target watershed with very high accuracy. In addition, for pervious land, characteristics such as infiltration capacity can be set specifically for each land use type [40].

In Japan, 10 m resolution grid-based data classified by land use type are generally used in distributed runoff modeling in urban watersheds. Although there are 17 land use categories in the classification, including industrial land, commercial/business land, and general low-rise residential land, the original purpose was to obtain the basic data necessary for residential land-related policies in urban areas, and these land use categories were not compiled with differences in infiltration rates in mind. Since estimating the impervious area ratio (IAR) for each land use category is an important component of accurate runoff analysis, Koga et al. [40] overlaid this landscape GIS and grid on a typical urban watershed in Tokyo and calculated the IAR for each grid. The results revealed that the reference value of the IAR for the land use classification had an error of about 7%.

Since the urban landscape GIS considers individual buildings, roads, parking lots, and other factors, Araki et al. [128] added stormwater runoff control facilities for individual buildings and specific permeable paved roads in the upper Kanda River basin as attributes of their maintenance status and examined their effectiveness in controlling flood runoff. The simulation results show that the developed model allows for the quantitative estimation of the efficiency of infiltration facilities.

The construction of urban landscape GIS data currently requires time-consuming manual work using GIS software (ArcGIS 10) [129]. Therefore, research has been conducted on automating urban landscape GIS delineation [130,131]. In addition, current GIS data only represent the planar structure of ground surface objects, and thus cannot represent the three-dimensional structure of subways, underground road tunnels, and the multilevel intersections of roads and bridges, as well as buildings, especially high-rise buildings. In order to more faithfully reproduce rainfall runoff paths physically, it will be necessary to construct three-dimensional GIS data that also take the height direction into consideration. The construction of a four-dimensional GIS that includes not only three spatial dimensions but also time has been pointed out as a challenging issue.

#### *5.8. TSR Model*

#### 5.8.1. Outline of Tokyo Storm Runoff (TSR) Model

The Tokyo storm runoff (TSR) model is a distributed physical model that utilizes urban landscape GIS data (for details, see Amaguchi et al. [6]). Figure 3 shows the rainfallrunoff process for the TSR model. Rainfall in the watershed is applied to each land use element, and the rainwater that cannot infiltrate is discharged to nearby roads and further to manholes and sewer pipes based on the infiltration and imperviousness characteristics set individually. In a manhole/sewer pipe element, the water level in the manhole section is first calculated based on the inflow/outflow to/from the road and inflow/outflow to/from the connecting pipe. The flow rate is then calculated based on the water level in the manhole section and the cross-sectional characteristics of the pipe. In this calculation process, water in the manhole overflows onto the road when the water level in the manhole exceeds the ground level of the road. The overflowing water flows down the road and eventually returns into a manhole or sewer pipe with sufficient flow capacity. Thus, water entering the sewer pipes is combined with water in the sewer network and eventually flows out into the river channel. In the inundation analysis, the amount of runoff from land use elements and the amount of overflow through manholes in the sewer system is used to calculate the water flows into residential areas when the road water level is above the surrounding ground level. The water in the residential areas flows back into the road when the water level on the road decreases.

Stormwater runoff processes in urban watersheds are highly dependent on local boundary conditions such as buildings, roads, rivers, and surface flows to sewers. Based on these boundary conditions, the TSR model calculates the direct runoff from roads, buildings, and other impervious areas that constitute a city and simulates the transport to the ground surface, sewerage systems, and rivers. When compared to the other four representations of urban structure (grid, curvilinear, unstructured, and road networks), and the associated models, the TSR model is not suitable for analyzing flooding flows with levee failures from large rivers, but it is intended for inundation phenomena, such as pluvial flooding in urban areas.

**Figure 3.** Schematic of the rainfall-runoff process.

#### 5.8.2. Model Application to Urban Watershed

Amaguchi et al. [132,133] used the building element of urban landscape GIS data to set individual values for the site elevation and floor elevation of buildings and made some modifications to the flooding analysis in the TSR model. This made it possible to simulate inundation conditions both inside and outside buildings. The model was set up and evaluated for the upper Kanda catchment in the Tokyo Metropolis, Japan. In the model, the parameters for the foundation height and floor height of each building, as well as the opening ratio below and above the 1st floor, were set for individual buildings based on their structure and purpose (Figure 4). The model was applied to the major flood event of 4 September 2005. There was no flood damage in the target watershed due to the effect of the underground control basins, so the simulation was performed assuming that the control basins were not functioning. The validity of the model was confirmed by comparing the model with river water levels observed at upstream locations not affected by the underground control basins. The validity was confirmed through its ability to overall accurately reproduce the observed river water levels [132]. This model was applied to simulate major flood events and was able to reproduce the differences in the inundation depth of individual buildings, reflecting the additional input parameters. The use of the model to evaluate inundation conditions for specific building properties was also demonstrated. The results show that the suggested approach, which is based on a detailed reproduction of all the relevant elements in an urban watershed, can simulate the inundation of individual buildings (Figure 5).

**Figure 4.** Parameters used for building inundation.

**Figure 5.** Example of inundation analysis results in the upper Kanda River watershed (The arrows indicate north direction).

#### **6. Concluding Remarks and Future Prospects**

In this review, approaches to flood runoff analysis and modeling with a focus on urban rivers in Japan have been reviewed. Relevant essential articles quoted in this paper concerning urban flood modeling in Japan have been summarized chronologically in Table A1 in Appendix A.

The USF model and other lumped conceptual models predict runoff at a given point in a river using rainfall data in the target urban watershed as input. The major advantages of these methods are that the number of parameters is small, the execution time is short, and the reproducibility of the observed hydrograph is generally high.

Concerning distributed modeling, in the TSR model with urban landscape GIS delineation, all land surfaces are represented by polygons divided according to their infiltration characteristics. This concept almost eliminates the problem of the excessive fine graining required in grid-based models to fully represent heterogeneity. However, current GIS data only represent the planar structure of ground surface objects and, thus, cannot represent the three-dimensional structures of subways, underground road tunnels, the multilevel intersections of roads and bridges, and (high-rise) buildings. In order to more faithfully reproduce rainfall runoff paths physically, it will be necessary to construct three-dimensional GIS data that also consider the height direction.

Furthermore, in flood runoff analysis using distributed physical models, it is important to improve observation systems for model calibration and the real-time prediction of flooding and inundation. These observations include spatial real-time rainfall radar information, the water levels in sewers, and the inundation depth during flooding. For this purpose, technology to quickly collect and utilize accurate information is needed. Additionally, runoff and inundation forecasts need to be transmitted in an easy-to-use format and be presented in an easy-to-understand manner to different recipients. Recently, the use of drones to observe early flood damage has been proposed to help local residents decide about evacuation and efficiently promote self-help. Finally, to enable highly accurate real-time flood forecasting with sufficient lead time using any type of runoff model, accurate and reliable high-resolution rainfall forecasts are required at least one hour ahead of an event in the case of urban watersheds.

Especially in lower urban areas, seamless, coupled models will be developed that incorporate not only runoff, channel, and sewerage models but also tidal models, which provide boundary conditions for flooding and inundation, and coastal hydraulic models to simulate storm surges, tidal waves, and tsunamis. In addition, this seamless, coupled modeling will incorporate atmospheric observations and simulations, such as radar rainfall information and land surface schemes. These seamless models, or "Digital Twins", will move beyond the scope of a single urban watershed to encompass a wide range of watersheds, including rural and mountainous watersheds (ultimately, the entire globe).

Finally, it is crucial to recognize that flood runoff models play a significant role in mitigating the impact of destructive floods through the implementation of appropriate protection measures, including flood control planning and policy. As there are very different kinds of urban floods worldwide, they require different protection measures [134]. The implementation of integrated flood risk management (IFRM) is also adopted in both developed and developing countries, yet some countries have encountered barriers to IFRM adaptation [135]. In Japan, a new policy called 'River Basin Disaster Resilience and Sustainability by All' was recently introduced in response to climate change risks and social changes [136]. It takes comprehensive measures, mainly consisting of flood prevention, exposure reduction, and appropriate evacuation, response, and recovery, aiming to strengthen disaster resilience and achieve sustainability through concerted efforts among all stakeholders [136]. Under the present circumstances, various rainfall-runoff models are becoming increasingly important tools, especially in urban watersheds.

**Author Contributions:** Investigation, H.A. and A.K.; Writing—original draft, A.K.; Data curation, H.A., J.O., A.K. and H.T.; Methodology, H.A. and J.O.; Funding acquisition, H.A. and J.O.; Writing—review and editing, H.A., J.O. and A.K.; Supervision, A.K. and H.A. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Japan Society for the Promotion of Science (BR200501) and the Scandinavia-Japan Sasakawa Foundation (GA20-SWE-0065).

**Data Availability Statement:** All data are provided as tables and figures.

**Acknowledgments:** This work was supported by Tokyo Metropolitan Government Advanced Research Grant Number (R4-2).

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **Appendix A**

**Table A1.** Chronological list of essential articles relevant to urban flood modeling in Japan. (Parenthetical Ref. No. is not developed in Japan).



**Table A1.** *Cont.*


#### **Table A1.** *Cont.*


#### **Table A1.** *Cont.*

#### **References**


**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
