Copula-Probabilistic Flood Risk Analysis with an Hourly Flood Monitoring Index

: Floods are a common natural disaster whose severity in terms of duration, water resource volume, peak, and accumulated rainfall-based damage is likely to differ significantly for different geographical regions. In this paper, we first propose a novel hourly flood index ( SWRI 24 − hr − S ) derived from normalising the existing 24-hourly water resources index ( WRI 24 − hr − S ) in the literature to monitor flood risk on an hourly scale. The proposed SWRI 24 − hr − S is adopted to identify a flood situation and derive its characteristics, such as the duration ( D ), volume ( V ), and peak ( Q ). The comprehensive result analysis establishes the practical utility of SWRI 24 − hr − S in identifying flood situations at seven study sites in Fiji between 2014 and 2018 and deriving their characteristics (i.e., D , V , and Q ). Secondly, this study develops a vine copula-probabilistic risk analysis system that models the joint distribution of flood characteristics (i.e., D , V , and Q ) to extract their joint exceedance probability for the seven study sites in Fiji, enabling probabilistic flood risk assessment. The vine copula approach, particularly suited to Fiji’s study sites, introduces a novel probabilistic framework for flood risk assessment. The results show moderate differences in the spatial patterns of joint exceedance probability of flood characteristics in different combination scenarios generated by the proposed vine copula approach. In the worst-case scenario, the probability of any flood event occurring where the flood volume, peak, and duration are likely to exceed the 95th-quantile value (representing an extreme flood event) is found to be less than 5% for all study sites. The proposed hourly flood index and the vine copula approach can be feasible and cost-effective tools for flood risk monitoring and assessment. The methodologies proposed in this study can be applied to other data-scarce regions where only rainfall data are available, offering crucial information for flood risk monitoring and assessment and for the development of effective mitigation strategies.


Introduction
Flooding is a catastrophic natural disaster progressively increasing in frequency and severity, primarily attributed to climate change-induced phenomena such as increased rainfall intensity.Generally, there are three prevalent flood types: fluvial or river floods, pluvial or flash floods, and coastal floods [1,2].A flash flood is a sudden and severe local inundation often resulting from high-intensity rainfall (e.g., tropical cyclones, slow-moving tropical depressions, or thunderstorms) within a short period (usually less than six hours) and/or may also be caused by sudden discharge of impounded water (e.g., dam or levee failures, ice jam release, or a glacier lake outburst) [3,4].Flash floods can affect a range of locations, including river plains, valleys, and areas with steep terrain, elevated surface runoff rates, constrained stream channels, and persistent heavy convective rainfall [3].They often necessitate prompt action to mitigate their severe impact, typically relying on expeditious decision-making and emergency response.Flash floods make up about 85% of all floods, resulting in more than 5000 deaths annually and causing severe social, economic, and environmental impacts [5].The repercussions of flood disasters are more devastating in developing countries such as Fiji [6], where this study is focused.Therefore, developing a real-time flood risk monitoring tool remains an ongoing research motivation to enable an assessment of flood occurrences for early warning systems in Fiji.
Fiji experiences regular flooding events arising from orographic rainfall due to the topography of its larger islands, including Viti Levu and Vanua Levu, which have a maximum elevation of 1300 m above sea level, along with the impact of prevailing southeast trade winds [7].For Fiji, about 90% of its population resides in coastal areas susceptible to floods [8].Between 1970 and 2016, Fiji faced 44 major flood events, impacting approximately 563,310 people and resulting in 103 fatalities [9].The most catastrophic floods occurred in 2004, 2009, 2012 (including the January and March flood events), and 2014.The 2009 and 2012 events, considered among the worst in the nation's history, resulted in over 200 million FJD in damages and losses, causing 15 fatalities and directly affecting more than 160,000 people [9].For Fiji, the estimated average annual flood losses exceed 400 million FJD, equivalent to 4.2% of Fiji's Gross Domestic Product (GDP) [9].These are substantial losses for a small island nation with a population of less than a million and a GDP of less than 5 billion USD [10].Under the assumption that climate change conditions will significantly increase rainfall, the annual flood-related asset losses could exceed 5% of Fiji's GDP by 2050 without adaptation measures [9].As a result, it is imperative to develop reliable methods for accurately monitoring flood risk on near real-time (e.g., hourly) timescales to mitigate the severe impacts of flooding.
Intense and/or prolonged precipitation is among the primary causes of floods.However, to better understand flood risk, multiple factors must be considered.These include the combined impact of various weather or climate events such as temperature and precipitation (while these events may not individually reach extreme levels, their cumulative effect can result in severe impacts) as well as the vulnerability and exposure of the affected area [11,12].Climate change also induces changes in various flood-related factors, including precipitation, soil moisture content, sea level, and glacial lake conditions, potentially changing flood characteristics [12].Other factors, such as land use and cover, catchment size and shape, drainage networks, and dam or levee construction, can also influence flood dynamics.Hence, an index that integrates all these factors and the accumulative impacts of other weather or climate characteristics is essential for a comprehensive flood risk assessment.
In many developing nations where flood monitoring resources, hydro-meteorological datasets, and risk monitoring facilities are underdeveloped, applying a mathematically derived flood index utilising only rainfall data provides a key strategy for assessing an impending flood risk situation.Some of the key mathematical indices used previously in flood risk monitoring include the Standardised Precipitation Index (SPI), the Available Water Resources Index (AWRI), the Weighted Average of Precipitation (WAP), the Standardised WAP (SWAP), the Flood Index (I F ), and the Standardised Antecedent Precipitation Index (SAPI) [13][14][15][16][17][18].The flood indices, such as AWRI, SWAP, I F , and SAPI, are robust as they are designed to account for changes in antecedent or immediate past rainfall by employing a suitable time-dependent reduction function that accounts for the depletion of water resources through various hydrological processes.For example, the daily flood index, I F , applied in Australia [17,19,20], Iran [21], Bangladesh [22,23], and Fiji [6], has shown good performance in monitoring flood events on a daily scale.Despite its benefits, one primary weakness of I F and other indices, such as SPI, lies in their utilisation of daily, monthly, or annual accumulated rainfall data, which represent much longer timescales than what is required in a flash flood monitoring system.Consequently, these indices fail to adequately represent the flood risk caused by bursts of high-intensity rainfall and rapid responses leading to flash flood events.
The present study draws relevance from a pilot study conducted by Deo et al. [24] that proposed a 24-h water resources index (WRI 24−hr−S ) based on a concept similar to the AWRI, which was applied in two study locations, Australia and South Korea, to monitor the flash flood risk in sustained extreme rainfall periods.The WRI 24−hr−S monitors flood risk by considering the contribution of accumulated rainfall in the past 24 h, whereby the rainfall contribution from the preceding hours is subjected to the time-dependent reduction function that accounts for the depletion of water resources through various hydrological processes such as evaporation, percolation, seepage, runoff, and drainage.However, unlike SAPI and I F , which are normalised values derived from the Antecedent Precipitation Index (API) and the AWRI, respectively, the identification of flood events and the computation of their characteristics are not achievable with the current form of WRI 24−hr−S .This is primarily because this index is unnormalised and does not enable objective assessment of flood risk across geographically diverse study sites.
To enhance the understanding of flood risks, it is essential to calculate the flood characteristics, including the volume (V), peak (Q), and duration (D) that concurrently result in major collateral damage.As the flood characteristics such as the D, V, and Q are mostly interrelated, we envisage that these characteristics should be jointly considered in a multivariate analysis model to estimate the actual probability of a flood occurrence [18,25].Importantly, any model representing the joint distribution of D, V, Q, and other crucial flood characteristics, such as the onset and withdrawal of a flood event, can provide significant insights into the relative severity of any flood event.After the initial study of Sklar [26], copula-based models became attractive in modelling interrelated variables, albeit using a multivariate approach.As such, copulas can jointly model the distribution of flood characteristics such as the D, V, and Q, regardless of their marginal distributions or whether their dependence structure is linear or non-linear [18].Vine copulas have recently shown superior capabilities in modelling flood characteristics compared to traditional Archimedean and elliptical copulas [18,[27][28][29].Therefore, this research follows a recent study by Nguyen-Huy et al. [18], which used vine copulas to model the joint distribution of extreme flood characteristics derived using the SAPI in Myanmar.Their study has provided interesting insights into the probabilistic flood risk analysis.To the best of our knowledge, no prior research has applied vine copulas to analyse the probabilistic flood risk in Fiji despite floods being a catastrophic phenomenon on this small island nation.
The scientific contributions of this paper, with significant implications for flood risk monitoring and assessment, are threefold.Firstly, the paper advances the concept of the 24 h water resources index pioneered by Deo et al. [24] and formulates a novel hourly flood index (SWRI 24−hr−S ) (a normalised metric) by normalising the 24-hourly water resources index in such a way that enables the objective assessment of flood risk across geographically diverse study sites.Secondly, the present study adopts the SWRI 24−hr−S , which is computed using real-time hourly rainfall data from 2014 to 2018 obtained from the Fiji Meteorological Services (FMS), to evaluate its practical utility in identifying flood situations and computing their associated flood characteristics (i.e., D, V, and Q) for seven different study sites in Fiji.Thirdly, the present study develops the vine copula approach for the first time to model the joint distribution of D, V, and Q derived from the SWRI 24−hr−S for specific cases of Fiji's flood events to extract their joint exceedance probabilities for probabilistic flood risk assessment.Fiji is a Pacific Small Island Developing State (PSIDS) that frequently experiences recurrent flooding.The lack of advanced infrastructure and necessary data in Fiji makes continuous flood risk monitoring and assessment challenging.The main aim of the present study is, therefore, to develop a novel hourly flood index, SWRI 24−hr−S , using only hourly rainfall data (which are readily available for the present study sites) and to conduct a probabilistic flood risk assessment by modelling the joint distribution of extreme flood characteristics derived from the SWRI 24−hr−S for Fiji's case studies.Hence, the methodologies presented in this study aim to enhance and contribute to existing monitoring and early warning systems for flash floods in Fiji.Moreover, these approaches may also be applied in other flood-prone regions globally, particularly those facing similar data scarcity challenges.By adapting these methodologies, vulnerable communities can benefit from improved flood preparedness and mitigation strategies.
The rest of the paper is structured as follows.Section 2 provides information on the study area, the dataset used, and data pre-processing steps.It also encompasses the mathematical methodology for computing the hourly flood index and the flood characteristics.Additionally, it provides details on Vine copula models and equations used for computing the joint exceedance probability of flood characteristics.Section 3 provides the results and discussion.Section 4 highlights the key findings, the study's limitations, and insights for future research.

Study Area
The proposed copula-probabilistic flood risk analysis system based on the hourly flood index was applied to geographically diverse sites in Fiji.Fiji is located in the South Pacific Ocean at a latitude of 15°S to 22°S and a longitude of 177°W to 174°E (Figure 1), with a tropical maritime climate characterised by warm temperatures throughout the year [30,31].The nation comprises an archipelago of 332 islands, 111 of which are permanently inhabited, with a total land area of about 18,333 km 2 [30].Viti Levu (10,400 km 2 ) and Vanua Levu (5540 km 2 ) are two large mountainous islands covering about 87% of the total land area [31].Fiji has a distinct dry season (May-October) and wet season (November-April).This seasonal variation is mainly attributed to the South Pacific Convergence Zone (SPCZ), the primary rainfall-producing system for the region, which lies typically over Fiji during the wet season [31,32].The rivers and stream basins in Fiji are predominantly small in size and flow from steep mountainous terrain, resulting in rapid shifts in water levels during periods of intense rainfall, which can lead to flash floods within a few hours [33].This study included sites only in Viti Levu due to the lack of rainfall data in other locations.These sites are areas in Fiji prone to recurrent and severe flooding events.Figure 1 shows a map of the study area and the corresponding study sites.

Dataset
The rainfall data for the Lautoka, Sigatoka, Nasinu, Rakiraki, Navua, Nadi, Ba, and Tavua sites from 1 January 2014 to 31 December 2018 (5 years) were successfully acquired from Fiji Meteorological Services.The rainfall data were provided in 5 min intervals for Tavua, Rakiraki, Nasinu, Sigatoka, and Lautoka and in 10 min intervals for the Ba, Nadi, and Navua sites.During the data pre-processing phase, the rainfall data for each site were aggregated to obtain the hourly rainfall needed for this study.If at least 66.67% of the data points (i.e., at least 4 out of 6 data points for a 10 min interval or at least 8 out of 12 data points for a 5 min interval) were available for a particular hour, they were summed to determine the total rainfall for that hour.Otherwise, the rainfall value for that hour was marked as missing.This approach was adopted to maximise the recovery of data values.
Table 1 summarises the hourly rainfall datasets and geographic settings of the study sites.The Ba site, which had a high percentage of missing values, was excluded from this study.The remaining sites had less than 5% missing values; therefore, any gap-filling method could fill in the missing values [34].Based on the study by Oriani et al. [35], the Iterative K-nearest Neighbour (IKNN) technique was used to fill in all the missing data.The data from 1 January 2014 to 31 December 2018 were used for all the computations.However, WRI 24−hr−S followed by SWRI 24−hr−S were calculated from 2 January 2014 as antecedent rainfall of 24 h (the hourly rainfall data for 1 January 2014), which was necessary to allow the calculations of these metrics.This research formulates a novel hourly flood index (SWRI 24−hr−S ), which is a normalised version of the 24-hourly water resources index (WRI 24−hr−S ) proposed by Deo et al. [24].Applying this index to Fiji is significantly advantageous because it is mathematically derived using only hourly rainfall data, which are readily available for the present study sites.The proposed SWRI 24−hr−S is implemented using the Python programming language.
The following steps are taken to obtain the SWRI 24−hr−S .The first step is calculating the WRI 24−hr−S .The WRI 24−hr−S for the current (i th ) hour is given by the following equation [24]: where P 1 is the total rainfall recorded an hour before, P 2 is the total rainfall recorded 2 h before, and so on; W is the time-reduction weighting factor (W ≡ 1 + 1/2 + 1/3 + ... + 1/24 ≈ 3.8) that incorporates the contributions of accumulated rainfall in the latest 24 h [24].This weighting factor ensures that the decay of accumulated rainfall or its potential impact on a flood event depends on several hydrological effects such as evapotranspiration, percolation, seepage, runoff, drainage, etc., in accordance with earlier works [15,24].The substitution of W = 3.8 into Equation (1) yields the following: Notably, the WRI 24−hr−S for a current (ith) hour is expected to accumulate 100% of rainfall received an hour before, ≈74% of that received two hours before, ≈61% of that received three hours before, and eventually ≈2% of that received 24 h before.
After calculating WRI 24−hr−S for any study period, the mathematical form of SWRI 24−hr−S for a current (ith) hour is calculated as a normalised version of WRI 24−hr−S : where WRI max 24−hr−S is the mean monthly maximum values of WRI 24−hr−S for the respective study period and σ(WRI max 24−hr−S ) is the standard deviation of the monthly maximum values of WRI 24−hr−S for the respective study period.
For the purpose of this paper, we follow the notion that if the magnitude of SWRI 24−hr−S for the current (ith) hour is greater than zero (or that the water resources are higher than normal), it is regarded as a flood situation.In this paper, we defined flood characteristics using the running-sum methodology of Yevjevich [36], which has also been used in several other studies [6,17,18].With reference to Figure 2, the flood duration, D, is estimated as the number of hours between the start of a flood, t onset , and the end of a flood, t end , derived from the SWRI 24−hr−S time series.In accordance with the onset of a flood, the volume of the flood, V, is calculated as the sum of all values of SWRI 24−hr−S between t onset and t end of a flood situation, and the peak of the flood, Q, is determined as the maximum SWRI 24−hr−S during any flood situation.Equations ( 4)- (6) show the mathematical equations used to calculate the flood characteristics before developing the copula-probabilistic flood risk analysis model.
where SWRI 24-hr-S > 0 D = t end − t onset (hours) where t onset and t end are the onset and end timestamps of a flood situation, respectively.
To demonstrate the practical use of SWRI 24−hr−S for hourly risk flood monitoring, Figure 2 illustrates the SWRI 24−hr−S applied to identify flood events in January 2014 at the Nadi site located in the western division of Fiji.As illustrated in Figure 2, the onset timestamp of the flood situation, i.e., the exact hour when the magnitude of SWRI 24−hr−S starts to rise above zero, was on 29 January 2014 at 8 a.m.To verify this particular flood situation, we now refer to the report from the FMS [37], which showed indeed that an active trough that caused widespread rain across Fiji was noticeable from 29 to 30 January 2014 and resulted in flooding, particularly in the western division of Fiji, where this study site is located.Thus, this verification confirms that the proposed SWRI 24−hr−S has correctly identified this flood event, demonstrating its practicality in identifying a flood situation at an hourly scale.
To further verify the potency of SWRI 24−hr−S for hourly flood risk monitoring, in Figure 3, WRI 24−hr−S is plotted alongside the hourly rainfall data for the same study site during the same period.Compared to SWRI 24−hr−S or the raw hourly rainfall data, SWRI 24−hr−S simplifies the process of identifying a flood situation.This is because a simple criterion, SWRI 24−hr−S > 0, provides a good indicator of flood risk, which is impossible when using WRI 24−hr−S and the raw hourly rainfall values.

Joint Exceedance Probability between Flood Characteristics
For a comprehensive flood risk assessment, this study follows the original approach of Nguyen-Huy et al. [18] to develop vine copula-based joint exceedance probability models.This task entails developing a multivariate analysis system of flood characteristics that considers the joint exceedance probability of a flood duration D, volume V, and peak Q for the present study sites.This study specifically aimed to estimate the probability that the duration, volume, and peak were concurrently greater than or equal to some threshold scenarios, as presented below: Equation ( 7) requires the modelling of a joint distribution function of three variables, F(x d , x v , x q ).Thus, in this study, we have developed a copula-based model, described in the following section, to estimate the joint exceedance probability of the flood characteristics, i.e., D, V, and Q to perform a probabilistic flood risk analysis. [26]:

Copula Analytical
with the corresponding joint density distribution function (JPDF) in terms of marginal and copula probability density functions [26]: where f i (x i ) and c(.) are the corresponding marginal and copula PDFs, respectively.When the marginal distributions are continuous, a unique copula exists.Equations ( 8) and ( 9) demonstrate an advanced capability of copulas, allowing a JCDF of random variables to be constructed through two separate processes: (i) modelling a copula function that captures the dependence structure among correlated variables and (ii) modelling the univariate marginal distributions.This aspect of copulas presents a more flexible approach for choosing suitable univariate distribution functions to fit the observed data in practical applications.From Equation ( 8), the joint distribution of duration, volume, and peak given in Equation ( 7) can be written using copulas as follows: Different copula families, such as empirical, Archimedean, extreme value, elliptical, vine, and entropy copulas, can be used to model the copula function given in Equation (10) [38,39].Vine copulas, among other copulas, can be used to achieve the utmost flexibility in constructing the JCDF and JPDF, given in Equations ( 8) and ( 9), respectively.Vine copulas have been applied in recent studies across various fields, such as weather and climate risk in agriculture [40,41], hydrology and water resources [27,[42][43][44][45][46][47][48], and finance and insurance [49][50][51].The following section provides more details on vine copulas.

Vine Copulas
The vine copula was first introduced by Joe [52], whose concept was to decompose the JPDF into a cascade of iteratively conditioned bivariate copulas, also called pair copulas.While this decomposition is not unique, all possible decompositions can be organised into a graphical model called a regular vine (R-vine) [53].
Within the R-vine framework, two main types of vine copula decomposition exist: the canonical (C-vine) and drawable (D-vine) distributions.These modes determine the parametric construction of an R-vine.The D-vine copula offers higher flexibility than the C-vine copula, especially when considering all mutual inter-correlations between targeted random variables one after another [47].However, the C-and D-vine frameworks are the same when considering a 3-dimensional (3D) (or tri-variate) joint distribution framework [18,47].
In this study, we have focused on tri-variate cases to model the joint distribution of flood duration, volume, and peak for a detailed probabilistic risk analysis of any flood event.Figure 4 shows a graphical representation of D-and C-vine copulas in the form of trees, edges, and nodes.In the tri-variate case where D, V, and Q are modelled simultaneously, the C-vine copula is the D-vine copula with a specified centre variable, as previously mentioned.For instance, when the flood duration (D) variable serves as the centre variable, the D-vine copula depicted in Figure 4b is identical to the C-vine copula shown in Figure 4d.The edges are linked to bivariate copulas, such as the edge DV associated with the bivariate copula C DV , which captures the dependence structure between D and V.
To fit the univariate marginal distribution functions, we employed the univariate local-polynomial likelihood kernel density estimation method capable of handling discrete (duration) and continuous (volume and peak) data [54].Additionally, the following bivariate copula families were utilised to develop the 3D vine copula models in this study: independence, parametric (elliptical, Archimedean, and their rotated versions), and nonparametric (transformation kernel) families [54-56] (Table A1).
To estimate the parameters of bivariate copulas, we employed maximum likelihood for parametric models and local-likelihood approximations for non-parametric models.Moreover, the modified vine copula Bayesian information criteria (mBICv) [57] was utilised to select the bivariate copulas, and Kendall's tau (τ) was adopted to select tree sequences [18].The mBICv can address the issues of the Bayesian information criterion (BIC), which assumes that the number of possible parameters grows sufficiently slowly with the sample size n and that all models are equally likely.Additionally, the mBICv was explicitly tailored to vine copula models [57].The vine copula models were developed using the R programming language utilising the 'rvinecopulib' library package [54].

Results and Discussion
We now provide a detailed appraisal of the hourly flood index SWRI 24−hr−S for detecting hourly flood possibility in terms of the onset and the end timestamps, duration, peak, volume, and total accumulated rainfall during any flood situation for seven study sites in Fiji over the study period (2014)(2015)(2016)(2017)(2018).We also provide joint distribution model results using the newly proposed vine copulas to provide a probabilistic risk analysis framework for flash floods.

Application of the Hourly Flood Index for Flood Event Analysis
The WRI 24−hr−S , followed by the SWRI 24−hr−S for the study period (2014-2018), were successfully computed for each of the seven study sites.The practicality of SWRI 24−hr−S for determining a flood situation has already been demonstrated in Figure 2. Similarly, the flood events between 3 and 6 April 2016 were quantified.This was first done for the Tavua site as it was one of the severely flooding-impacted areas in the western division of Fiji [58].
Our analysis identified four flood events using the criterion SWRI 24−hr−S > 0 to indicate a flood situation (Figure 5) for the Tavua site.Of four flood events, the two significant events were predominantly caused by heavy rain in the past 24 h.The first major flood event started on 3 April at 5 p.m. and ended on 4 April at 4 p.m. with a total duration of 23 h, a volume of 20.37, and a peak of 2.36.Approximately 154 mm of rainfall was recorded during this event.The second major flood event started on the 6th of April at 3 a.m. and lasted for 14 h.This flood event recorded a total volume of 8.77, with a peak of 1.03, while approximately 69.50 mm of rainfall was recorded for this event.The combined volume of all four flood events for the Tavua site between 3 April 2016 and 6 April 2016 (96 h) was 30.59.The flood events for the same period were also determined for the other six sites.The flood characteristics, i.e., D, V, and Q of the floods, varied among these sites, as shown in Figure 6.The results showed that areas in the western division of Fiji were severely affected by flooding, as was reported by FMS [58] [Tavua (V ≈ 30.59),Lautoka (V ≈ 25.63), Sigatoka (V ≈ 19.45), Nadi (V ≈ 9.12), and Rakiraki (V ≈ 8.79)] compared to the areas in the central division [Navua (V ≈ 0.12) (minor flood event), and Nasinu (no floods)].These results demonstrate the utility of SWRI 24−hr−S in identifying flood situations and determining their characteristics.Consequently, the proposed SWRI 24−hr−S can be considered a feasible and cost-effective tool to monitor the flash flood risk in Fiji.The variation in flood characteristics among our study sites demonstrates the importance of flood risk assessments for each site separately, despite the proximity of these sites, as also highlighted in previous work [6].Figure 7 depicts the occurrence of floods, encompassing minor events with minimal volume and potentially insignificant impacts at each of the seven study sites over 2014-2018.Over this five-year study period, a slight fluctuation in flood frequency was observed across the study sites, as illustrated in Figure 7. Notably, the Tavua study site exhibited the highest frequency, while the Nasinu study site recorded the lowest frequency of flood events.The occurrence of frequent severe weather events, such as tropical cyclones and depressions, results in significant flood events in Fiji, and this usually occurs during the wet season (November-April) and occasionally in the dry season (May-October), especially in La Niña years [33].This is evident in Figure 8, which indicates the wet season, including May and October, experiencing high rainfall, leading to a higher frequency of flood events and greater flood volume (Figure 9) compared to the other months.This emphasises the need for Fiji's National Disaster Management Office (NDMO) and other relevant stakeholders to implement comprehensive flood mitigation and resilience measures.Public education on flood safety and preparedness for the wet season is also crucial, particularly for those residing in flood-prone areas.The annual rainfall and the occurrence of flood events across seven sites from 2014 to 2018 are illustrated in Figure 10.This figure shows that the year 2015 had the lowest rainfall among all years examined.According to FMS [59], the rainfall trends in 2015 exhibited a typical El Niño pattern, characterised by below-average rainfall at most of the study stations.Consequently, there were fewer flood events (Figure 10) and a lower flood volume (Figure 11) in 2015 compared to the other years in the present study.Table 2 lists the five severe floods at each of the seven sites during the study period.The flood severity was determined by ranking the flood events at each site based on their volume, with 1 indicating the most severe and 5 indicating the least severe.For each of the seven sites, the table displays the onset time, duration, volume, peak, total WRI 24−hr−S , total rain, and maximum WRI 24−hr−S for each flood event.Statistics such as these may aid relevant organisations in understanding past flood events at these sites, which will facilitate future flood mitigation decisions to minimise the severity of floods at these locations.
A brief analysis of the most severe flood event at each study site was performed and validated using the annual climate summaries published by the Fiji Meteorological Services to ensure that the SWRI 24−hr−S accurately identified them.The analysis of floods in Nadi (Table 2a) showed that the most severe flood started on 29 January 2014 at 8 a.m. and recorded a volume of 157.28.This flood lasted 49 h and reached a peak of 6.80, making it the most severe flood event with the most prolonged duration among all the study sites during the 5-year study period.During this flood event, about 1590 mm of rainfall was recorded.This flood event was described in Figure 2 and validated using Fiji's annual climate summary for 2014 [37].
Table 2b shows that the most severe flood event in Lautoka recorded a total volume of 25.05.This flood event started on 14 January 2018 at 2 p.m. and lasted 19 h, reaching a peak of 3.29.During this flood event, about 175 mm of rainfall was recorded.This flood event was verified using Fiji's climate summaries 2018, which stated that heavy rainfall occurred from 13 to 15 January 2018 due to an active trough of low pressure, resulting in widespread flooding [60,61].According to Table 2c, the most severe flood in Nasinu started on 27 February 2014 at 9 a.m. and lasted 18 h.This flood had a volume of 24.90 and reached a peak of 2.83.
Approximately 173 mm of rainfall was recorded during this flood event.As reported by FMS [37], the tropical depressions TD14F and TD15F caused heavy rainfall in Fiji's central and eastern divisions between 25 and 27 February 2014.As a result, parts of Fiji, particularly the major river systems in the central division (where this study site is located), experienced flooding during this period.
The most severe flood, both in Navua (Table 2d) and Rakiraki (Table 2e), occurred in December 2016.For Navua, this event started on 15 December 2016 at 6 a.m. and lasted for 28 h, during which it recorded a volume of 56.98 and reached a peak of 4.29.For Rakiraki, it started on the 19th of December at 5 a.m. and lasted for 21 h, during which it recorded a volume of 33.99 and reached a peak of 4.28.Approximately 392.50 mm and 265.50 mm of rainfall were recorded during this flood event for Navua and Rakiraki, respectively.The most severe flood event that occurred at the Navua and Rakiraki sites was validated using Fiji's climate summaries 2016, which stated that the trough of low pressure and active rain bands associated with the tropical depression TD04F resulted in heavy rainfall that caused severe flooding in some parts of the country's main island of Viti Levu (where these study sites are located) from 12 to 20 December 2016 [58].
Based on Table 2f, the most severe flood event in Sigatoka started on 30 January 2014 at 10 a.m. and lasted for 16 h.This flood had a volume of 23.32 and reached a peak of 2.84.Approximately 121 mm of rainfall was recorded during this flood event.As mentioned earlier, an active trough that caused widespread rain across Fiji from the 29 to the 30 January 2014 resulted in flooding, particularly in the western division of Fiji, where this site is located [37].
Lastly, the analysis of floods in Tavua (Table 2g) showed that the most severe flood started on 8 February 2017 at 10 a.m. and recorded a volume of 45.88.This flood lasted 23 h and reached a peak of 4.04.During this flood event, about 238 mm of rainfall was recorded.As per FMS [62], the tropical depression TD09F affected the country between 6 and 8 February 2017 and led to flooding in parts of the western division of Fiji, where this study site is located.

Application of the Vine Copula Model for Probabilistic Flood Risk Analysis
The frequency of flood events at each study site is demonstrated in Figure 7. Similarly, the flood characteristics, i.e., D, V, and Q, for each flood event were calculated for all study sites.Table 3 shows the five-number summary, including the mean, standard deviation, skewness, and kurtosis for each flood characteristic at each study site.
Moreover, as shown in Table 3, the minimum flood duration was 1 h at all study sites.The maximum flood duration, volume, and peak were recorded at the Nadi site (this flood event is described in Figure 2).The skewness and kurtosis information of each flood characteristic, which describes the shape and distribution of a dataset, were more than +1 and +3, respectively, for all study sites, indicating that their distribution is highly right-skewed.This means the flood characteristics dataset for all study sites contains extreme flood duration, volume, and peak values.
The results in Table 3 also show that flood characteristics exhibit high variability across all study sites in terms of their median and inter-quartile range (IQR).The median flood duration for all study sites was 3 h, while the median volume and peak ranged from approximately 0.52 to 0.89 and 0.24 to 0.42, respectively.The IQR for flood duration, volume, and peak varied from 4 to 7 h, 1.78 to 2.94, and 0.52 to 0.91, respectively.The high spatiotemporal variation in flood characteristics highlights the importance of modelling these characteristics simultaneously, and employing a robust model like the copulas used in this study is crucial for accurately capturing the dependence among these flood characteristics.Conducting a comprehensive correlation analysis among flood characteristics and understanding the relationship between each characteristic pair is another crucial step in modelling their joint distribution [18].In this study, both parametric measures-Pearson's correlation coefficient (r)-and nonparametric rank-based correlation measures-Spearman's correlation coefficient (ρ) and Kendall's tau (τ)-were employed to examine the relationship between each pair of flood characteristics.Additionally, Mutual Information (MI) was utilised to examine the degree of dependence between each pair of flood characteristics.The results obtained are presented in Table 4.The correlation coefficients, i.e., r, ρ, and τ, obtained between each pair of flood characteristics, were statistically significant at the 1% level of significance.Overall, a strong positive dependency is observed between each pair of flood characteristics across all study sites, as shown in Table 4.However, there would be cases where the peak value is extremely high while the duration is low or the peak is moderate while the duration is long, and the volumes would also be different.For example, as shown in Table 4, at the Tavua site, the linear correlation measured by Pearson's correlation coefficient between D and Q (r D&Q ≈ 0.859) is higher than that between D and V (r D&Q ≈ 0.838).However, the rank correlation measured by Spearman's and Kendall's tau correlation coefficients between D and Q (ρ D&Q ≈ 0.902, τ D&Q ≈ 0.698) are lower than those between D and V (τ D&V ≈ 0.946, τ D&V ≈ 0.836).This implies complex and non-linear relationships among flood characteristics.Subsequently, the risks of different flood events are different.Therefore, this requires a robust model like copulas used in this study to capture the full dependence among flood characteristics.
Additionally, it must be noted that across all study sites, the strongest dependencies exist between D − V and V − Q compared to D − Q.This is particularly evident based on Kendall's tau (τ) coefficient utilised in this study to select the most optimal structure of the vine copula model at each study site.This implies that the flood volume, V, can be positioned between the other two flood characteristics (i.e., D and Q), as illustrated in Figure 4a, to model the joint distribution of D, V, and Q using the 3D D-vine copula.Table 4.The statistical correlation in terms of the Pearson's correlation coefficient (r), Spearman's rank correlation coefficient (ρ), Kendall's tau (τ), and Mutual Information (MI) computed between the pairs of flood characteristics, i.e., Duration (D, in hours), Volume (V), and Peak (Q) for each study site.Table 5 shows the results obtained when the 3D D-vine copula is fitted to the flood characteristics data at each study site.As depicted in Table 5, the results confirm that the D-vine structure illustrated in Figure 4a with flood volume (V) as the conditioning variable is the most appropriate to model the joint distribution of flood characteristics (i.e., D, V, and Q) across all study sites.The table also shows the best-fitted bivariate copula function and its associated parameters at each tree level for each study site.For instance, at the Sigatoka site, in the first tree (Tree 1), the Frank (C DV ) and Survival Gumbel (C VQ ) copulas are selected between D − V and V − Q, respectively.In the second tree (Tree 2), the Frank copula is chosen as the most parsimonious for identifying the bivariate copula (C DQ|V ).

Site
To derive the joint exceedance probability of the flood event characteristics (i.e., D, V, and Q) for different combination scenarios using the best-fitted 3D D-vine copula selected at each study site, we first quantify the probability that the flood duration, volume, and peak exceed specific thresholds simultaneously (Equation ( 7)).The thresholds were selected at the 50th quantile (median), 75th quantile (moderate), and 95th quantile (extreme).The quantile values of each flood characteristic were computed and subsequently averaged for all study sites, as presented in Table 6.As seen in Table 6, for example, the averaged 50th-quantile value of duration is q D (0.5) = 3 h.Similarly, the averaged 75th-quantile value of duration is q D (0.75) = 6 h, and the averaged 95th-quantile value of duration is q D (0.95) = 15 h.As for the spatial pattern, Table 6 demonstrates a moderate variation in the quantile values of each flood characteristic across all study sites.By applying the vine copula probabilistic model, we show the joint exceedance probabilities of the duration, volume, and peak in different combination scenarios for each study site in Figures 12-14.From a flood risk analysis perspective, the present results clearly demonstrate a moderate yet notable difference in spatial patterns of the joint exceedance probability of flood event characteristics in different combination scenarios.As shown in Figure 12a, the probabilities of a flood event occurring where both the volume and peak exceed the 50th-quantile (median) values (i.e., V ≥ q (0.50) = 0.622 and Q ≥ q (0.50) = 0.352) and the duration (D) exceeds the median (i.e., D ≥ q (0.50) = 3 h), moderate (i.e., D ≥ q (0.75) = 6 h), and extreme (i.e., D ≥ q (0.95) = 15 h) values are approximately 50-59%, 23-39%, and 4-10% across all study sites, respectively.
In general, the probability of a flood event with a volume exceeding the 50th-quantile (median) or 75th-quantile (moderate) values and both the peak and duration exceeding the 95th-quantile (extreme) value was less than 5% across all study sites.In the case when both the flood volume and duration exceeded the 95th-quantile (extreme) value, the probability of a flood event with the peak exceeding the 50th-quantile (median) or 75th-quantile (moderate) values was less than 6% across all study sites.
In the worst-case scenario, when the flood risk could be more severe, we found that the probability of a flood event occurring where the volume, peak, and duration exceeded the extreme values (i.e., V ≥ q (0.95) = 12.799, Q ≥ q (0.95) = 1.915, and D ≥ q (0.95) = 15 h) was less than 5% at all study sites (Figure 14c).These findings imply a moderate probability of a flood event characterised by median (i.e., 50th-quantile) duration, volume, and peak values across all study sites.The results also suggest that the likelihood of a flood event characterised by extreme duration, volume, and peak is exceptionally low across all study sites.

Conclusions, Limitations of the Study and Future Research Directions
This study has made novel contributions to flood risk monitoring and assessment by developing a mathematically convenient hourly flood index, SWRI 24−hr−S , and testing its practical use in identifying flood situations over 2014-2018 at seven different study sites in Fiji while jointly modelling flood characteristics such as flood duration, volume, and peak using a vine copula model for probabilistic flood risk assessment.
The results have unambiguously established the practical use of the newly proposed SWRI 24−hr−S as a potent indicator to identify the flood situation at an hourly scale while computing the associated flood characteristics that were impossible with a 24-hourly water resources index used in the literature.The results also showed that Fiji mainly experienced high rainfall during the wet/cyclone season (November to April), including May and October.Consequently, the number of flood events was higher in these months than in the other months.This highlights the critical importance of implementing comprehensive and well-structured flood preparedness and risk mitigation strategies tailored explicitly for these months characterised by increased rainfall and flood events, thus ensuring the safety and security of communities and their properties.This study also presented the flood characteristics and water-intensive properties of five severe flood events for each study site.Relevant organisations, such as Fiji's NDMO, are expected to use these findings to understand the attributes of past flood events at these study sites.This, in turn, can support future decision-making on flood mitigation, ultimately reducing the severe impacts of such events.
The results also demonstrated a strong positive dependency between each pair of flood characteristics across all study sites.The D-vine structure with flood volume (V) as the conditioning variable was identified as the most appropriate for modelling the joint distribution of flood characteristics (i.e., D, V, and Q) across all study sites.Therefore, it was utilised to model the joint distributions among the extreme flood characteristics to extract their joint exceedance probability, providing crucial information for probabilistic flood risk assessment at each study site.The findings revealed moderate variations in the spatial patterns of joint exceedance probability of extreme flood event characteristics across different combination scenarios, underscoring exceptionally low probabilities of floods with extreme duration, volume, and peak at all study sites.
Despite the merits of the present study, a primary limitation of this research was the unavailability of rainfall data required for many of the flood-prone sites across Fiji.Consequently, this research was confined to selected sites within the western and central divisions of the country.As a result, this study could not perform a comparative analysis across all four divisions (i.e., the western, central, northern, and eastern divisions), which could have provided valuable insights into extreme flood risk areas in the nation.It is important to note that the Ba study site in the western division of Fiji, a frequent flooding zone, had to be excluded due to a high percentage of missing rainfall data.Therefore, in future research, our methodology can be improved with SWRI 24−hr−S derived from satellite-based rainfall products covering a wider area, including major towns and cities, following the recent approach for Myanmar [18].
Another limitation of the present study was using a prior/fixed time-dependent reduction function with the weighting factor, W = 3.8, derived in an earlier study [24] to determine the contributions of accumulated rainfall in the latest 24 h.As discussed, the proposed SWRI 24−hr−S is a normalised version of an existing WRI 24−hr−S that used a suitable time-dependent reduction function to account for the depletion of water resources through various hydrological processes.It must be noted that the results of this study are sensitive to the value of W. A change in W will alter the contributions of accumulated rainfall in the latest 24 h in Equation (2).Consequently, both WRI 24−hr−S and the proposed SWRI 24−hr−S will change accordingly.
In the future, further studies can test the correctness of this weighting factor (W) more comprehensively for study sites where the topography may vary considerably.This could require a major correlation of this weighting factor against rainfall-runoff and other physical models to capture more accurately the actual value of the decay of accumulated rainfall and its impacts on a flood event [15,24].Several tests with hydrological parameters, including evapotranspiration rates, percolation, seepage, surface runoff, and drainage conditions, etc., may be required to ascertain the time-dependent reduction function for SWRI 24−hr−S .In regions with different decay rates of rainfall-accumulated water volume, it is crucial to incorporate them when formulating SWRI 24−hr−S .The proposed SWRI 24−hr−S must also be verified for its broader adoption as an index-based risk monitoring tool.Therefore, its feasibility is expected to be demonstrated in other flood-prone regions globally in future studies, contingent on the availability of well-documented flood records for validation and hourly rainfall data.While doing so, a different technique to normalise the existing WRI 24−hr−S may be selected, depending on how the normalised WRI 24−hr−S index represents flood risk situations in those climatic conditions.
This study has undertaken a purely mathematical-based approach to monitoring flood risk, so in future studies, it is anticipated that the proposed SWRI 24−hr−S , in conjunction with additional data such as the catchment hydrology, drainage information, and river flows, will be utilised to develop hourly hydrographs for various sites.This approach will further cement the accuracy of flood characteristic estimation and the monitoring of flash flood events.There is also the potential to develop an innovative SWRI 24−hr−S -based forecasting system with sufficient lead time, presenting a novel approach for early flash flood warnings in Fiji and other regions.
A key advantage of SWRI 24−hr−S , as an hourly flood risk monitoring tool, is its simplistic mathematical formula that is easy to compute, analyse, and interpret for non-expert audiences compared to physical or hydrological models, including rainfall-runoff models for flood risk monitoring that involves complex development.However, in the future, especially in varied hydrological and topographic settings, it is crucial to comprehensively compare the proposed SWRI 24−hr−S with other established flood monitoring systems, including the Flash Flood Guidance System (FFGS).
Despite these limitations, using SWRI 24−hr−S has demonstrated acceptable accuracy in detecting flood situations on an hourly scale.Therefore, our proposed methodology can be considered a feasible and cost-effective tool for hourly flood risk monitoring in Fiji and perhaps other similar geographical locations.Applying the proposed probabilistic flood risk analysis using vine copulas can enhance the nation's overall flood risk assessment and mitigation strategies.

Figure 1 .
Figure 1.The geographical map shows Fiji's location in the South Pacific region and a detailed inset map highlighting various study sites within Fiji.

Figure 2 .
Figure 2. The SWRI 24−hr−S developed to identify a flood event in January 2014 at the Nadi study site, demonstrating its capability to determine the duration, volume, and peak of any flood event.The flood volume, representing the accumulated water resources, is the cumulative SWRI 24−hr−S under the curve closed by the onset and end of a flood situation and the zero line.

Figure 4 .
Figure 4. Schematic diagram of different ways of constructing the 3D D-vine copula structure in this study.(a) Case 1: the flood volume (V) as a conditioning variable, (b) Case 2: the flood duration (D) as a conditioning variable, (c) Case 3: the flood peak (Q) as a conditioning variable, and (d) an example of the C-vine copula structure.Source: Adapted from Nguyen-Huy et al. [18].

Figure 5 .
Figure 5.The SWRI 24−hr−S applied to identify the flood events in April 2016 at the Tavua site (96 h).

Figure 6 .
Figure 6.The SWRI 24−hr−S applied to identify the flood events in April 2016 at the other study sites (96 h).

Figure 7 .
Figure 7. Geographic analysis of flood frequency between 2014 and 2018.

Figure 8 .
Figure 8. Temporal (monthly) analysis of flood frequency and total monthly rainfall aggregated from 2014 to 2018.

Figure 9 .
Figure 9. Temporal (monthly) analysis of the combined volume of flood events across 7 study sites from 2014 to 2018.

Figure 10 .
Figure 10.Frequency of floods and total rainfall across 7 study sites aggregated from 2014 to 2018.

Figure 11 .
Figure 11.Yearly combined volume of flood events across 7 study sites from 2014 to 2018.

Table 1 .
Overview of the hourly rainfall datasets for the 8 sites in Fiji.(Note: The hourly rainfall spans from 1 January 2014 to 31 December 2018, with 43,824 expected observations.)

Table 3 .
Descriptive statistics of flood characteristics at each study site.