1. Introduction
Floods are a frequent and damaging natural disaster that negatively impacts the socioeconomic development and lives of millions of people worldwide [
1,
2,
3,
4,
5]. It is estimated that between 1990 and 2016, worldwide losses from flood damage amounted to USD 723 billion [
6]. Due to population growth and climate change, urban areas in particular are expected to be seriously threatened by the effects of increased flood intensity and frequency [
7,
8,
9,
10]. Approximately 40% of the world’s cities will be in areas under high flood risk by 2030, thereby affecting about 54 million people [
11]. Furthermore, 82% of urban areas in Southeastern Asia will be in high-frequency flood zones by 2030 [
12]. Understanding the complex relationships between urban growth and floods is essential in developing risk management strategies for sustainable land use planning [
13].
The flood system is very complicated and has the characteristics of spatial-temporal dynamics, with uncertainties and integration of different challenges in a system generating complex phenomena [
14]. Research is an essential tool in flood risk management since it plays an essential role in predicting flood hazard and improving societal comprehension of the complex environmental and socioeconomic components of flood risk [
15]. The objective of this study is to develop a comprehensive approach to estimating the evolution in flood risk between 1995 and 2040 in a rapidly changing land cover context in Vietnam.
Flood risk is a combination of hazard, exposure, and vulnerability [
16], so risk can be managed by reducing the flood intensity and/or damage [
17,
18]. Over the past few decades, with the development of science, the world has brought many changes in the approach to minimize flood effects. Traditional methods of flood control (i.e., structural measures such as dikes, embankments, etc.) are gradually being replaced by more comprehensive flood risk management models [
19,
20]. These approaches take into account the probabilities and potential consequences of flood events based on risk assessment studies where flood hazard and exposure/vulnerability factors are quantified [
21]. Several methods have been applied to compute these indicators. The hazard approaches include: measurement-based, field surveys [
22,
23], hydrodynamic models [
24,
25], and GIS and Remote Sensing [
26,
27] in linear modeling of flood risk through overlaying component layers with associated analytical hierarchical process (AHP)-based computed weights. Land cover indicators can be grouped into two categories: (1) traditional terrestrial mapping [
28]; and (2) land cover classification based mainly on satellite observations [
29,
30]. The emergence of satellite sensors like Landsat, Satellite Pour Observation de la Terre (SPOT), and Sentinel 2, among others, has greatly facilitated the rapid classification of land cover and its evolution. In addition, remote sensing has the advantage of rapid data acquisition with lower costs over field survey methods [
31].
Flood risk assessment is an essential step in defining appropriate management strategies [
32]. In recent decades, several studies have focused on developing flood risk assessment methods at different scales and with different objectives. Mishra et al. [
33] developed a flood risk index for the Kosi River of India based on hazard (geomorphologic, distance to the active channel, slope, and rainfall) and socioeconomic vulnerability (population, household, and female densities; literacy rates; land cover and use changes; road–river intersections; road density). Chinh Luu et al. [
34] studied the temporal changes of flood risk that integrated hazard, exposure, and vulnerability to better understand the evolving dynamics and formulate appropriate mitigation strategies. Dang et al. [
35] outlined the essential roles in improving flood risk assessment methods to support decision-making processes. The authors classified flood risk indices into three components: social–economic, physical, and environmental. Kron [
36] constructed flood risk indices based on the likelihood of flood occurrence, social–economic vulnerability, the environment, and flood consequences. Begun et al. [
37] integrated floods-related damage with the probability of their occurrence. Several clues or methods of assessing flood risk have been developed in different areas. However, they are limited in a comprehensive framework that can support decision-makers to understand the aggravating risk causes. In addition, these studies focus on assessing the flood risk at a specific time while, according to Penning-Rowell et al. [
38], flood risk reduction strategies are most effective when they are evaluated continuously. Jhong et al. [
39] emphasize that understanding hazard, exposure, and vulnerability at different dates is an essential task in flood risk reduction that allows land managers to better see the spatial and temporal trends likely to arise in the future [
40]. To fill the gaps identified, we integrate hydraulic modeling with land cover change analysis and prediction to assess flood risk at three dates: 1995, 2019, and 2040.
This research is different from the previous studies cited because it provides a novel and comprehensive approach to flood risk assessment based on state-of-the-art remote sensing and modeling techniques and assesses both historical and future trends. The initial hypothesis tested in this study is that flooding risk has increased substantially in the study area due to an increase in the population living in the flood plain. Flood hazard, exposure, and vulnerability were all accounted for in order to assess their potential interactions and compensatory effects. While this study explicitly examines flood risk in a catchment located in Vietnam, the findings are of importance to other rapidly evolving countries affected by floods and experiencing urban growth.
3. Data and Methods
In this study, flood risk is determined by a combination of flood hazard, exposure, and vulnerability. Hazard is a physical phenomenon—natural and uncontrollable with intense occurrence [
41]. It is an elementary notion that expresses the probability of a situation, an event, or some causality—a flood in this case [
36,
42]. It is directly proportional to the intensity of the phenomenon’s occurrence [
43], since the combination of flood depth and flood peak velocity represents the capacity to destroy objects in the area through which the flood passes, directly affecting houses, buildings, and people’s lives and health. Exposure is understood as values present at the location of a possible flood. These values can be goods, infrastructure, cultural heritage, people, or agricultural areas. This can be interpreted as the presence or availability of assets and people in flood risk areas. The level of exposure depends on the frequency of flood occurrence, the flood’s intensity, and the value of the available properties and people [
42,
44]. In this study, land cover and population density were selected to analyze the level of exposure. The land cover categories are detailed below and include Cropland, Forest, Urban area, and Water. Potential flood impacts in urban areas are considered greater since the infrastructure is well developed and they concentrate more value than other land cover types; they also have the greatest population densities. Agricultural land follows urban areas as it provides the livelihood for many of the watershed inhabitants. Population density is a critical variable of flood exposure analysis; because it is linked directly with humans, it is proportional with the level of flood vulnerability [
34]. Vulnerability represents the lack of individuals’ or groups’ ability or capacity to anticipate, counter, and resist floods [
45]. The International Plant Protection Convention (IPPC) emphasizes that these capacities depend on the prevailing economic situation and social and political characteristics. Vulnerability increases in direct proportion to the flood level and decreases socioeconomic wellbeing during floods [
46]. In this study, the poverty rates and number of hospitals were selected to build the vulnerability maps. The wealthy-to-poor ratio characterizes a community’s resilience; the wealthy have sturdier homes, have the ability to access information about danger through modern media, and quick resumption of their normal lives is enhanced by their good economic standing [
47]. Hospitals play a critical role during and after events. Patients place a heavy burden on communities and emergency management services in the event of an evacuation [
48]. A flowchart of the methodological workflow is displayed in
Figure 2. The quantification of hazards, vulnerability, and exposure is described in detail in the following sections.
3.1. Flood Hazard Estimation
Flow depth and velocity were used to estimate flood hazard. The hydraulic modeling tool MIKE FLOOD, which combines the MIKE 11 and MIKE 21 models, was used to build the depth and velocity map for the 2013 Quang Ngai historic flood.
3.1.1. Establishing the 1-Dimensional (1D) River Network
The 1D network is necessary for modeling water flow through the channels into the flood plains and was established for the Ve and Tra Khuc rivers with 46 cross-sections along 27.1 km of the Ve River (average distance about 0.59 km/sections), and 21 cross-sections along 51.6 km of the Tra Khuc (average distance about 2.6 km/sections). The cross-sections were measured directly in the field using total station (
Figure 1). The discharge observation data at An Chi and Son Giang gauges were applied to the upstream boundary of Ve and Tra Khuc, respectively, and the tide level was applied to the downstream boundary using tidal sea level harmonic constants. MIKE NAM is the rainfall runoff model part of the MIKE 11 module developed by the Danish Hydraulic Institute (DHI), Denmark. This model provided discharge for sub-basins without streamflow data, and sub-catchment boundaries are presented in
Figure 1. MIKE NAM was integrated into the MIKE 11 model and values were applied to the inflow boundaries for the MIKE 11 hydraulic model.
3.1.2. Establishing the River Network in the 2-Dimensional (2D) Hydraulic and FLOOD Models
The 2D channel network is necessary for estimating the spatial extent of the flooded areas when the river depth exceeds bank height. The area modeled covered an area of 507.2 km2 based on the surface elevation data and observed flood inundation area defined by flood marks. The 1:10,000 topographic map published in 2010 by the Ministry of Water Resources and Environment was used for defining the 2D domain’s elevation in MIKE-21. The 2D computational mesh was generated by discretizing the computational domain into 48,200 elements, with the size of elements ranging from 457 to 2500 m2.
MIKE FLOOD is used to link MIKE 11 and MIKE 21 through the Lateral link which connects the tops of riverbanks in the 1D model to the 2D model’s mesh elements. The water levels and time series of two flood events—one in November 2013 and the other in November 2017 at Tra Khuc gauge—served to calibrate and validate the model. Parameters for the 2D model were calibrated based on flood trace measurements taken in November 2013 (50 locations) and November 2017 (32 locations). Nash–Sutcliffe efficiency (NSE) [
49], flood peak error, and coefficient of determination (R
2) were used to measure the model’s reliability.
3.2. Indicators of Flood Exposure
Exposure was quantified based on property value and the population [
22,
24]; land-use categories and population density were selected to estimate these variables.
3.2.1. Land Cover Mapping in 1995, 2019, and 2040
Land Cover Mapping in 1995 and 2019
Land cover maps were created using satellite imagery—SPOT-3 (12/03/1995) with a 20 m spatial resolution and Sentinel-2 (27/02/2019) with a 10 m resolution, for 1995 and 2019, respectively. The data were projected to a common datum and coordinate system (WGS84/UTM Zone 48N). Both images were acquired with cloud cover less than 5%. Images underwent radiometric and geometric corrections before re-sampling all bands to a 10 m spatial resolution using the bilinear interpolation sampling method.
The eCognition Developer of Trimble was used to perform object-oriented classification on the SPOT-3 and Sentinel-2 images; a sample output of the classification process is displayed in
Figure 3. First, image objects were segmented by integrating similar pixels; then, each segment was assigned to land cover object layers [
50,
51]. In order to optimize the object-oriented classification, it was necessary to select the appropriate segmentation parameters as each segment must be homogeneous and separate from its neighbors [
52].
Values for the 3 segmentation parameters were 400 for Scale, 0.2 for Shape, and 0.99 for Compactness. After segmentation, the objective variables were selected using random formation points, and the information related to the objects was extracted. Different values of objective characteristics were derived from spectral indices such as the Normalized Difference Vegetation Index (NDVI) and the Enhanced Vegetation Index (EVI). The samples and extracted objective information were then exported to a calculation table. The object-oriented classification was based on the training polygons with 5 land-use categories: Cropland, Developed (built area), Bare soil, Forest, and Water bodies.
The land cover map in 1995 was compared to panchromatic aerial photographs collected from the Ministry of Natural Resources and Environment (MONRE 1995), while the 2019 classification was compared to the 2020 land use map (MONRE 2020). Validation points were randomized with n = 205 for the map in 1995 and n = 262 for the map in 2019; the number of sampling points for each class depended on the class coverage.
Change Prediction in 2040
The initial map had five categories, but Bare soil was mostly sand bars in the river or exposed beach areas along the coast, so it was integrated into the Water body category to facilitate land cover change modeling.
Future land cover change was modeled using the Land Change Modeler (LCM) module from Terrset [
53] using the default Multi-Layer Perceptron (MLP) neural network option for predicting future change. LCM facilitates the quantification and mapping of historical land cover changes and predicts future changes based on past transition rates. Quantity of change from one land cover to another is based on Markov chain analysis, and spatial allocation depends on transition potentials derived from historical trends and explanatory variables such as topographic or distance variables [
53,
54]. Driver variables generally belong to one of three categories: accessibility, suitability, and zoning [
55]. Accessibility is frequently modeled as the distance from roads or existing urban areas, and suitability for urban areas is strongly dependent on topography. There were no zoning constraints or initiatives in the study area. The initial explanatory variables tested were the following: altitude (10 m DEM), slope inclination (%), and distances from roads, developed area in 1995, and from water (river or sea). However, since slope inclination and altitude are closely correlated in the study area, only altitude was retained. The contribution of explanatory variables to predicting land cover change was estimated from Cramer’s V values and from accuracy rates produced during the creation of potential transition layers [
53]. Accuracy rates are calculated by running calibrating transitions on a sub-sample of cells and then evaluating the accuracy of predictions on the remaining cells. A land cover map for 2040 was produced from the transition potential maps.
3.2.2. Population Density
Population density is considered one of the critical indices for flood exposure. Population density by municipality was obtained from the municipalities for 1995 and 2019, and an estimate for 2040 by the municipalities, was obtained from the statistics and general offices of six districts in the study area (Quang Ngai, Son Tinh, Mo Duc, Tu Nghia, Nghia Hanh, and Binh Son).
3.3. Indicators of Flood Vulnerability
Poverty data by municipality and number of hospitals in the flood zone were obtained from six Departments of Statistics (Quang Ngai, Son Tinh, Mo Duc, Tu Nghia, Nghia Hanh, Binh Son) to build the vulnerability map as described below.
3.4. Assigning Normalized Weights for Hazard, Exposure and Vulnerability Using Analytical Hierarchical Process (AHP)
The components of flood risk are hazard, exposure, and vulnerability [
36,
56,
57]. Each component was divided into sub-components and the AHP method was used to quantify the weight of each sub-component [
35].
AHP is a quantitative method used to organize decision options and choose an option that satisfies a given criterion. Developed by Thomas L. Saaty in the 1970s, it is considered an efficient and flexible method for analyzing multi-criteria decisions [
58]. It operates by setting priorities for multi-criteria rankings that are judged by experts involved in the decision-making process to derive the best decision possible. In this study, assessments were based on the authors’ group experience and from previous similar studies [
35,
59,
60,
61,
62,
63]. The weights of flood risks in 4 steps:
- (1)
Constructing component and sub-component hierarchies: hazard, exposure, and vulnerability. Hazard is divided into two sub-components (depth and velocity), exposure into two sub-components (land cover and population density), and vulnerability into two sub-components (poverty rate and number of hospitals).
- (2)
Establishing priorities: The relative importance of the elements in each pair of sub-components was evaluated subjectively and assigned a value ranging from 1 to 9 as per the Saaty scale [
64].
- (3)
After preparing a pairwise comparison matrix (
Table 1), each column was divided by the corresponding sums to obtain the priority factors. This process is a normalized Eigenvector of the matrix (
Table 2); the average values of the row were used as the priority of the subcomponents to calculate flood risk.
- (4)
Estimation of the Consistency Ratio (CR).
The Consistency Ratio (CR) is used to ensure consistency in the experts’ judgement throughout the application and is applied as follows:
CR = CI/RI, where Consistency Index (CI) can be calculated as per the equation CI = (λmax − n)/(n − 1); λmax is the comparison matrix’s eigenvalue obtained by multiplying each parameter’s correlation matrix column’s total summation by the normalized value of corresponding parameters; and n is the number of elements compared in pairs during a calculation.
RI must be defined in [
64]
After calculating the relative importance and determining the weight of each factor in the hierarchy, flood hazard, exposure, and vulnerability were calculated using the following equations:
Flood hazard = 0.75 flood depth + 0.25 velocity
Flood exposure = 0.66 population density + 0.34 land cover
Vulnerability = 0.66 poverty rates + 0.34 number of hospitals
The output values were normalized, ranging from 0 to 1, and reclassified into five classes using the natural break method: Very low, Low, Moderate, High, and Very high.
5. Discussion
Flood risk analysis is a critical step towards sustainable economic development and protection against flood damage [
67,
68]. This study introduces a comprehensive flood risk assessment method that integrates hydraulic modeling, land cover change analysis and prediction, and socioeconomic trends. The flood risk maps of 1995, 2019, and 2020 have five risk levels (very low, low, moderate, high, and very high). These levels were classified by the natural break method, based on applying Jenk’s optimization formula to minimize the variability of each category, having the advantage of automatically defining the final classes, underlining disparities in the best way possible [
69,
70]. Using different breakpoints would perhaps have divided the flood zone area slightly differently but the overall trend of diminishing risk over time would have remained the same.
The final results were unexpected as our initial hypothesis was that flood risk was increasing with increasing population density in the flood zone and with the transition from Cropland to Built-up area. This trend is typical of most cities undergoing rapid urbanization in a flood plain. Huu Duy et al., 2018 [
71] reported that rapid urbanization growth, in addition to poor planning in the Gianh River watershed in Vietnam, has resulted in many populations becoming more vulnerable to floods. Areu-Rangel et al. [
72] found urbanization was an important factor in increasing flood risk. The work by Waghwala [
73] in the case of Surat City hammers in this fact—changes in LULC caused by urban expansion increase the flood risk. Bahrawi et al., 2020 [
74] analyzed the flood risk in East-Ern Jeddah, Saudi Arabia. The authors pointed out that the rapid growth of urbanization is radically changing the characteristics of flooding and increasing the flood risk in this region. These results were confirmed by Zhang et al., 2018 [
75]; according to the authors, the exacerbation of urbanization is not the only one to pose difficulties to the response to the floods, but also the total precipitation during the storms. Handayani et al., 2020 [
76] reported that poorly planned urbanization has increased the surface area exposed to floods. Mustafa et al., 2018 [
77] indicates that urbanization increases the risk of floods in the future due to the growth of population and infrastructure in the flood-prone area. This was confirmed by Neumann et al., 2015 [
78]; the authors point out that population growth and urbanization increase vulnerability and risk. In 2060, Egypt, Nigeria, China, India, Bangladesh, Indonesia, and Vietnam are considered to be the countries most vulnerable to floods due to population growth. However, in this study, the reduction in poverty between 1995 and 2019 substantially increased the population’s capacity to resist and bounce back from floods, therefore vulnerability, and consequently flood risk, diminished over time. This has been justified by research in other countries such as that by Rayhan et al., 2010 [
79]. Tasnuva et al., 2020 [
80] highlighted that poverty is considered one of the main factors of vulnerability. Thus, our study brings additional evidence for integrating socioeconomic status for flood risk assessment in the context of urbanization.
The trend identified here raises two major questions that must be addressed with regard to this specific study: variable weightings in the AHP method and the consequences of an extreme flood. Variable weightings were carried out based on the authors’ experience and the scientific literature; the weightings have a subjective component that is difficult to evaluate. Flood hazard modeling gave particularly reliable results based on the calibration and validation event statistics; therefore, flow velocity and depth values are considered robust. However, the weightings attributed to these parameters in the AHP remain subjective, both for within hazard weighting (
Table 1) and between risk component variables (
Table 2). The same can be said for both exposure and vulnerability. In the end, the change in poverty level was the dominant factor accounting for an improvement in the flood risk levels mapped in
Figure 11, and two further comments can be made about this. First, as cited above, a reduction in poverty has been shown to improve flood resilience in other studies, so this is not entirely surprising. However, the weights attributed to poverty level in
Table 1 and
Table 2 are difficult to justify objectively since an increase in wealth is translated into a range of considerations that are beyond the limits of this study: home improvements that make houses more resistant to floods, insurance coverage that minimizes losses for personal homeowners and businesses, savings and other investments that allow homeowners and businesses to renovate or rebuild quickly, municipal wealth that can reestablish and/or rebuild municipal infrastructures quickly. These elements could not be quantified here and have been little explored in the scientific literature, particularly for Asian countries, but they represent a critical field of investigation in flood risk management, especially in countries where poverty rates are evolving quickly.
A second aspect related to the weight of poverty rate in determining the final flood risk map is related to the unexplored role of spatial scale in the data. Flood hazard and land cover were mapped at 10 m spatial resolutions. Population density and poverty levels were provided according to municipal boundaries where large areas were attributed identical values. The spatial limits of the municipal boundaries extend beyond the flood zone limitations; therefore, the results suggest that increases in population density were perhaps not as great in the flood zone as they were elsewhere in the study area. Hence, the important change in population density noted in the text above (overall increase of 47.3% between 1995 and 2040) is representative of municipalities within the watershed but not necessarily within the flood zone where population density may have changed less. A few municipalities even experienced minor decreases in population density, but it was not possible to find population density at a spatial scale compatible with the 10 m spatial resolution used in the study. Harmonizing socioeconomic data with physical (DEM) or remotely sensed raster data remains a permanent challenge in land cover change modeling.
In the case of an extreme event, like the November 2013 reference event, we could expect flood damage to be greater over time despite the decrease in flood risk. The flood zone is increasingly occupied by more people and land covers with greater economic value, so the damage from a given event would be expected to be of greater financial value. However, “flood risk” in our study is not limited to expected damage but to the overall capacity of the system to bounce back from that damage—its resilience. As noted by Fox et al. (2012) [
81] in SE France, improvements in the river channel networks largely compensated the impacts of peri-urbanization on runoff in the catchment; there was a net increase in exposure, but the flood risk decreased. Despite this, in the case of an extreme event, more people and buildings would be affected. Therefore, although flood risk decreased in the Tra Khuc watershed thanks to improved socioeconomic conditions, flood prevention measures, largely ignored to date, must become a priority for the catchment.
The remaining question of the study is whether this novel approach may resolve the problems of previous studies related to sustainable land use planning. This question is very important because the “urban” status of the territory is the subject of an administrative decision, and other studies from countries experiencing economic transitions connect decisions with the socioeconomic dynamics [
82,
83,
84]. Other studies have used approaches similar to ours to assess flood risk. Huu Xuan Nguyen et al., (2020) [
85] used GIS-Based Fuzzy AHP–TOPSIS for assessing flood haz-ards along the South-Central Coast of Vietnam. Chinh Luu et al., in 2020 [
59] used 300 flood marks in 2013 and a 5 m DEM to build a flood hazard map using the AHP method. Although these methods have some advantages, such as easy access to data, our study is more comprehensive compared to previous studies thanks to the use of hydrodynamic and land cover modeling to assess historical and future flood hazard and thanks to the integration of socioeconomic data such as poverty rate. Integrating socioeconomic changes to estimate flood risk was an essential part of this study, and it confirms that appropriate attention must be given to changes in resilience over time in estimating flood risk despite the uncertainties around AHP weights. Over the past 25 years, Vietnam has experienced rapid economic transition from an economy dependent on natural resources such as agriculture, fishing, and forestry to one based on industry and services [
71]. The role of a continuous flood risk assessment from the past to the present is very important. This is a missing strategy in flood risk management in general and especially in developing countries like Vietnam.
The economic transition policy since 1986 has led to rapid urbanization. The coastal plains in Vietnam in general and in the study area are significantly affected by urbanization. This has been justified by several studies [
71,
86]. Zaninetti [
86] emphasizes that people do not hesitate to occupy and prosper in areas afflicted by natural hazards such as floods and storms. Apprehension towards these dangers is offset by the daily usefulness of the land, fertility of the soil, presence of aquaculture and natural resources such as building materials, or simply the availability of flat, easy-to-build-on land. This transition will have serious consequences on people and assets, especially when no proper assessment strategy is in place. Therefore, the research presented in this paper is an appropriate flood mitigation strategy in the context of rapid urbanization. Our results are a useful contribution to the theoretical advancement of the field because they can assess flood risk at different times, including future probable trends. Based on this, policymakers can identify priority regions in need of spatial planning. Although this study was carried out in Vietnam, the results apply to other countries with flood challenges, particularly to regions with rapid socioeconomic changes and urbanization.
Finally, land cover change studies for flood risk assessment require accurate historical land cover maps [
87]. A long-term, consistent database of satellite images offers researchers the opportunity to analyze phenomena from a historical perspective, and it is possible to assess long-term changes in local natural parameters. In addition, it provides a synoptic view of large areas [
88,
89,
90]. Technically, the selection of an effective method to detect changes in land cover is considered an important step. Object-Based Image Analysis (OBIA) was selected to analyze the satellite images in this study. This method gives better results when compared to traditional pixel-based techniques for image classification, so it was used in this study. Pixel-based techniques often produce more “noise” (often referred to as a “salt-and-pepper” effect) than object-based approaches due to differences in reflection between adjacent pixels. The supervised OBIA technique helps researchers to be proactive about the size of the subjects they study (e.g., fine, medium-coarse), thereby choosing the appropriate level of segmentation for their research. In addition to multi-scale measures, OBIA can incorporate other thematic data (e.g., DEM, thematic maps, etc.) into the classification process to help filter out noise. Hence, object-based classification results will bring more accurate results.
This study is subject to the general limitations characteristic of the use of topography data in hydrodynamic modeling. Several previous studies have shown that the accuracy of hydrodynamic modeling depends on topography resolution. This study used the DEM extracted from the 1:10,000 topography map; however, future research would benefit more when more detailed information is used, e.g., LiDAR, because it presents information on land cover such as dike networks or vegetation to compute the friction coefficients [
91,
92]. In this study, several methods were used to limit these disadvantages, for example the artificial surface like dike networks or the main roads were measured and integrated into the 1D hydrodynamic modeling.