Simulating flash floods at hourly time-step using the SWAT model

: Flash ﬂoods are natural phenomena with environmental, social and economic impacts. To date, few numerical models are able to simulate hydrological processes at catchment scale at a reasonable time scale to describe ﬂash events with accurate details. Considering a ~810 km 2 Mediterranean river coastal basin (southwestern France) as a study case, the objective of the present study was to assess the ability of the sub-daily module of the lumped Soil and Water Assessment Tool (SWAT) to simulate discharge (1) time-continuously, by testing two sub-basin delineation schemes, two catchment sizes, and two output time-steps; and (2) at ﬂood time-scale, by comparing the performances of SWAT to the performances of the event-based fully distributed MARINE model when simulating ﬂash ﬂood events. We showed that there was no beneﬁt of decreasing the size of the minimum drainage area (e.g., from ~15 km 2 down to ~1 km 2 ) when delineating sub-basins in SWAT. We also showed that both the MARINE and SWAT models were equally able to reproduce peak discharge, ﬂood timing and volume, and that they were both limited by rainfall and soil data. Hence, the SWAT model appears to be a reliable modelling tool to predict discharge over long periods of time in large ﬂash-ﬂood-prone basins.


Introduction
Floods are natural phenomena with environmental, social and economic impacts [1][2][3].Impacts of floods may be magnified by some human activities, such as increasing settlements and economic assets in floodplains, or the reduction of the soil's natural water retention under changing land uses.Floods are also responsible for the transport of organic and inorganic contaminants in both dissolved and sorbed phases [4][5][6][7][8][9].Global change is expected to increase the flooding risk in the next decades, with more people settling in floodplains, whereas climate change is likely to increase the frequency of intense rainfall events and consequent floods, especially across the prominent climate response hot-spot Mediterranean basin [10][11][12].
Among floods, flash floods are generally triggered by quasi-stationary convective systems, with local and intense rainfall events [13,14].The term "flash" refers to the rapid response, with water Water 2017, 9, 929 2 of 25 depths in the drainage network reaching peak levels within minutes to a few hours after the onset of the rain event [15].In particular, typical Mediterranean shallow soils with sparse land cover and steep slopes in coastal basins' uplands produce instantaneous and violent runoffs because of limited water infiltration, especially when the soil is already saturated at the beginning of the storm [16].Indeed, flash floods in the northwestern Mediterranean coastal basins are of major environmental and socio-economic concern because of the basins' distinctive topography and geographical locations.The region experiences yearly heavy rainfall events and flash floods that result in up to hundreds of millions of euros in damages and often casualties [2,13,[17][18][19].
To date, few numerical models are able to simulate hydrological processes at basin-scale at a reasonable time scale to describe these flash events with accurate details [20][21][22].For flash flood prediction, an accurate simulation of the flood peak timing, flood peak magnitude, and total exported water during floods is required.Among basin-scale hydrological models, one can distinguish between lumped or distributed models, and between conceptual or process-oriented models.Event-based models simulate fine time periods with intervals as small as several seconds.These models are generally formulated with mechanistic process-oriented equations, and the solutions are numerically approximated.A generally recognized drawback of event-based models is their need for initial conditions (e.g., initial soil water content).As a result, they are typically used only for small catchments [20,23].Distributed models such as the MARINE event-based model [22], are able to simulate hydrological processes at the flood time-scale in flash-flood-prone areas such as the Mediterranean coastal basins [16,[24][25][26], but they disregard inter-event processes such as river discharge during low flow, evapotranspiration, and deep groundwater flows.On the other hand, lumped models such as the conceptual Soil and Water Assessment Tool (SWAT) [27,28] often include more processes and are able to simulate long periods of time, including low flow and high flow periods.The ease of implementing lumped models and their relatively small computational costs are also an advantage.
The capability of simulating individual floods is important for watershed models to adequately capture hydrological processes between short intervals (e.g., flood peak discharge), while time-continuous simulation is required to investigate, e.g., inter-flood processes such as groundwater contribution to river flow and contaminant exportation, or long-term impacts of changes on watershed hydrology, such as land use change or climate change [21].Therefore, the time-continuous, lumped SWAT model has recently been upgraded to sub-daily time-step calculations [20,21].However, the SWAT sub-daily module has to date mostly been tested in small catchments (~1 km 2 ) and mostly urbanized areas [20,21,[29][30][31][32] and rarely over large basins (5800 km 2 , [33]).In addition, even when operated at the sub-daily time step, SWAT has mostly been used for daily runoff and daily stream flow simulation [29,30,32,33].Hence, the upgraded SWAT model has not yet been tested to simulate sub-daily discharge in large basins (~1000 km 2 ), including large non-urbanized areas, such as the reactive Mediterranean coastal basins.
Considering the Têt Mediterranean river basin (810 km 2 , southwestern France) as a study case, the objective of the present study was to assess the ability of the sub-daily module of the SWAT model to simulate discharge (1) time-continuously, by testing two sub-basin delineation schemes (1 and 15 km 2 minimum drainage areas), two catchment sizes (90 and 810 km 2 ), and two output time-steps (daily and hourly); and (2) at flood time-scale, by comparing the performance of SWAT to the performance of the process-oriented, event-based, fully-distributed MARINE model when simulating flash flood events.

Study Design
We ran SWAT at sub-daily time-step for the following cases: (1) comparison of daily and hourly time-step discharge simulation outputs with daily and hourly time-step discharge measurements, respectively; (2) comparison of discharge simulation outputs using small and large sub-basin delineation areas (1 and 15 km 2 ); and (3) application to a smaller sub-basin with no infrastructure Water 2017, 9, 929 3 of 25 (90 km 2 ).We also compared the performance of the time-continuous SWAT model to the performance of the event-based MARINE model by looking at the discharge statistics, flood timing and magnitude, and volume of water exported during 6 flash flood events, and soil water content before and after the floods.Both models were calibrated with the same input forcing, including sub-daily rainfall.

Hydrological Models Description
Major similarities and differences between the SWAT and the MARINE models are listed in Appendix A.

The SWAT Model
The SWAT model is a basin-scale, time-continuous model developed to evaluate the short-term and long-term impacts of management practices on water, sediment and agricultural chemical yields in ungauged basins [27,28].Key processes simulated by the model include hydrology, erosion, plant growth, nutrient and pesticide cycles, and land management.SWAT is a lumped (i.e., semi-distributed) model where the catchment is divided into sub-basins, which are further divided into hydrological response units (HRUs), which are homogeneous units of land use, soil type and slope.
The sub-daily runoff module of SWAT is based on the Green and Ampt [34] excess rainfall method, which was further modified by Mein and Larson [35] to calculate infiltration as a function of wetting front matric potential and effective hydraulic conductivity.The Green-Ampt-Mein-Larson (GAML) equation is described as: where f is the infiltration rate at time t (mm h −1 ), K e is the effective hydraulic conductivity (mm h −1 ) in which the impact of land cover is incorporated, Ψ is the wetting front matric potential (mm), ∆θ is the change in volumetric moisture content across the wetting front (mm mm −1 ), and F is the cumulative infiltration at time t (mm H 2 O).Therefore, the parameters from SWAT that are involved in the GAML equation are the saturated hydraulic conductivity (SOL_K), the curve number (CN2), the soil texture (SOL_CLAY, SOL_SAND, SOL_BD) and the soil available water content (SOL_AWC).
The GAML equation is a physically-based model that allows continuous simulation of the infiltration process assuming the soil profile is homogeneous and antecedent moisture is uniformly distributed in the profile [21].The GAML method requires sub-daily precipitation data.
SWAT then calculates the peak runoff rate with a modified rational method [36].The lateral subsurface flow in the soil profile is determined for each soil layer with the kinematic storage routing model [37], which is calculated simultaneously with percolation.The groundwater flow contribution to the total streamflow is simulated by creating shallow aquifer storage [38], where the percolation from the bottom of the root zone is considered as recharge to the shallow aquifer.
In the river channel, the flow is routed using either the variable storage coefficient method [39] or the Muskingum routing method [40].In this study, we used the variable storage coefficient method because it was adapted to the calculation of very small discharges during low flow periods (<1 m 3 s −1 along the river) over the study basin.The potential evapotranspiration can be estimated with various methods, including the Penman-Monteith method, which was used in this study.In this study, we used the revision 635 executable of the SWAT model.

The MARINE Model
MARINE is a distributed, parsimonious, event-based, rainfall-runoff model.It was developed for the simulation of flash floods, including in poorly-gauged basins, and can run in an operational mode for real-time forecasting.Key physical processes simulated by the model are soil saturation dynamics (including infiltration, subsurface runoff, exfiltration) and surface water flow (in runoff and open channels).Inter-event processes such as evapotranspiration and deep aquifer percolation are deemed negligible with respect to the time scale of a flash flood and therefore not modelled.The input data describing the catchment used by the process-oriented MARINE model are topographical data (slope and downhill directions), soil data (thickness and texture) and land cover data.No description of the river is required, as MARINE infers the drainage network as well as the geometric characteristics of the network reaches based on inputted topographical data and geomorphologic considerations [41,42].Initial conditions are specified by soil water content at the beginning of the event, while the base flow is set to zero.The use of distributed input, as for rainfall measurements and soil properties, enables the model to take advantage of the spatial information.The hydrological model operates over a regular grid of interconnected square cells (a 500 m grid resolution is typical for a catchment with a drainage area above 500 km 2 ).In each time step, MARINE calculates the water balance for every cell according to its location whether it is on the overland or in the drainage network.Local infiltration is computed using the Green and Ampt model [34].Subsurface flow is computed using an approximation of Darcy's law giving the flow per unit width q: where T 0 is the local transmissivity of fully saturated soil (m 2 s −1 ), θ s and θ are the saturated and local water contents (m 3 m −3 ) respectively, m is the transmissivity decay parameter (-), and β is the local slope angle (assumed to be equal to the slope of the water table in the saturated zone, in radians).
The value of m is set to 0.08 according to previous studies [43]; all other parameter values are derived from input data.Soil water may exfiltrate in two cases: when the soil water content θ exceeds saturated water content θ s , or when the soil water reaches the drainage network.Overland and channel flows are calculated using a 1-dimentsion kinematic wave approximation of the Saint-Venant equations with the Manning friction law.For the overland flow, the equation translates to: where h is the water depth (m), t is time (s), S is the bed slope (m m −1 ), n is the Manning's roughness coefficient (s m −1/3 ), x is the space variable (m), r is the rainfall rate (m s −1 ), i is the infiltration rate (m s −1 )-all values being derived from input data.The forcing function on the right-hand side of Equation (3) expresses the rainfall excess resulting from the difference between the rainfall rate and the soil infiltration rate (handled by the Green and Ampt equation as explained above).Flow for the drainage network is represented using the previous kinematic wave equation, taking into account the cross-section of the network reach.For a complete description of the MARINE model the reader can refer to [22].
To initialize the event-based model with soil water contents, we used the operational SIM (Safran-Isba-Modcou) application from Météo-France [44], which is traditionally used for MARINE initialization [22,45,46].

The Têt Watershed
The case study is the Mediterranean coastal Têt River basin (1400 km 2 ), located in the South of France (Figure 1a).The river length is about 120 km and the average discharge at the downstream Perpignan gauging station is 10 m 3 s −1 .The altitudes range from 2800 m a.s.l. in the Pyrenees Mountains down to sea level (Figure 1b,c).In the upper part of the basin the elevation above 2500 m a.s.l.makes the area regularly snow-covered during winter [47,48].Average annual rainfall across the basin is 780 mm [49] and long dry periods alternate with short and intense storm events.The main components of the basin's bedrock are granite and primary era formations, i.e., mainly schists but also locally highly karstified limestone (source: BD Million-Géol, BRGM) and altered bedrocks may lead to water losses across the basin [24].The basin is mostly covered by range land (35.9%) and forests (36.8%).Agricultural land covers 19.9% of the basin's area and includes vine, orchards and annual crops (Figure 1d).Soil hydraulic conductivities are low (~6 mm h −1 ) and soils are mostly shallow (~1.4 m) (Figure 1e,f).The basin embeds two reservoirs (Figure 1b): one in the upper part of the basin, for hydropower purpose (Bouillouses reservoir, 17 × 10 6 m 3 ), and another one in the lower part of the basin (Vinça reservoir, 24 × 10 6 m 3 ) for the purpose of irrigation and flood peak clipping, in order to protect the lower and flat alluvial area of the Têt River basin (Plaine du Roussillon).Irrigation in the Têt River basin has been developed since the ninth century.Many irrigation channels exist across the basin and the water can bypass gauging stations at several places [47].No information was available about the flows in irrigation channels in the upper part of the Têt River [50], which may especially impact the simulation of low flow discharge.

Selected Flash Flood Events
We selected all of the extreme events that occurred during the last decade (2006-2016), i.e., events where maximal discharges were higher than the two-year return period of 59 m 3 s −1 at Marquixanes gauging station.We considered eight flash flood events for MARINE calibration, six of which were used for the comparison with SWAT (Table 1).The basin is mostly covered by range land (35.9%) and forests (36.8%).Agricultural land covers 19.9% of the basin's area and includes vine, orchards and annual crops (Figure 1d).Soil hydraulic conductivities are low (~6 mm h −1 ) and soils are mostly shallow (~1.4 m) (Figure 1e,f).The basin embeds two reservoirs (Figure 1b): one in the upper part of the basin, for hydropower purpose (Bouillouses reservoir, 17 × 10 6 m 3 ), and another one in the lower part of the basin (Vinça reservoir, 24 × 10 6 m 3 ) for the purpose of irrigation and flood peak clipping, in order to protect the lower and flat alluvial area of the Têt River basin (Plaine du Roussillon).Irrigation in the Têt River basin has been developed since the ninth century.Many irrigation channels exist across the basin and the water can bypass gauging stations at several places [47].No information was available about the flows in irrigation channels in the upper part of the Têt River [50], which may especially impact the simulation of low flow discharge.
Hence, with steep upland catchments and sparse soil, the Têt River responds quickly to weather variations such as storms: it is a typical flash-flood-prone Mediterranean river [13].In this study, we focused on two embedded sub-basins of the Têt River basin (Figure 1b): the sub-basin gauged at the Marquixanes station (station Y0444010, 42 • 38 37.54 N, 2 • 29 9.49 E, 810 km 2 drainage area), and the sub-basin gauged at the Catllar station (Y0445010, 42 • 37 55.81 N, 2 • 25 18.54 E, 90 km 2 ).The latter is the outlet of the Castellane Têt tributary.Its discharge is not altered by any infrastructure.

Selected Flash Flood Events
We selected all of the extreme events that occurred during the last decade (2006-2016), i.e., events where maximal discharges were higher than the two-year return period of 59 m 3 s −1 at Marquixanes gauging station.We considered eight flash flood events for MARINE calibration, six of which were used for the comparison with SWAT (Table 1).In order to be able to compare the performance of the two models, the same input data were used whenever possible (Figure 1).

•
Topography was described for both models with the SRTM 90 m digital elevation model [ No outflow data was available for the Bouillouses reservoir.In SWAT, an inlet was thus created ~12 km downstream at the Mont-Louis gauging station (Y0404010, 42 • 32 19.99"N, 2 • 3 21.39E, 39 km 2 ) and daily discharge was fed with January 2005-December 2013 records from the regional flood forecasting service for the Languedoc Roussillon region (Service de Prévision des Crues Méditerranée Ouest-SPCMO).All measured aperiodic discharge data from 2006 to 2014 at the Marquixanes and Catllar gauging stations were provided by SPCMO.To compare SWAT periodic outputs to observations, we resampled observed aperiodic records to periodic ones by allocating the nearest recorded value to each hour.If no flow value was recorded in a period lasting longer than two hours, no value was allocated.
Hourly rainfall from January 2005 to December 2014 was collected from the operational rain gauge network for flood monitoring purposes provided by the SPCMO (Figure 1b).For the time-continuous SWAT model, additional daily temperature, humidity, wind, and radiation data were gathered from the SAFRAN (Système d'Analyse Fournissant des Renseignements Adaptés à la Nivologie) meteorological model, which data are organized in an 8 km × 8 km grid [53][54][55].At the time of the study, SAFRAN outputs were available from January 2005 to July 2014, hence the last complete year for the SWAT run was 2013.
Also, in the specific case of SWAT modelling, we used a typical SWAT watershed delineation scheme (15 km 2 minimum drainage area, Figure S1) and a very detailed watershed delineation (1 km 2 minimum drainage area, Figure S1) in order to assess the effect of the sub-basins' delineation on the quality of hourly time-step simulation results and to determine the effect of the SWAT delineation scheme on the comparison with MARINE results: on one hand, a finer delineation may capture the meteorological information from more climate stations and meteorological information may be better spatialized, since sub-basins are smaller (see Appendix A for the rainfall stations-sub-basins association); on the other hand, a finer delineation may more correctly represent the transfer dynamics at the sub-daily time-step and may approach the 500 × 500 m resolution of the MARINE calculation grid.In SWAT, the 1 km 2 delineation led to 671 sub-basins and 2290 HRUs whereas the 15 km 2 delineation led to 64 sub-basins and 536 HRUs.
With both delineations, five elevation bands were set up in sub-basins where the difference between the maximal and minimum elevation was over 700 m.The default setting of one elevation band was kept when the difference was under 700 m.

Sensitivity Analysis and Model Calibration and Validation
Calibration and validation of the two models were performed by comparing the simulated discharge to the measured discharge.For the SWAT and SWAT-CUP SUFI-2 modelling tools, we resampled the aperiodic records to hourly data to fit the precipitation input time-step (see also Section 2.3.3).For both SWAT time-continuous calibration and event-scale MARINE calibration, we used the Nash-Sutcliffe model efficiency (NS) as an objective function.
2.4.1.Procedure for the SWAT Model SUFI-2 [56][57][58] was used for calibration, validation, sensitivity, and uncertainty analysis of the SWAT model.SUFI-2 is an iterative algorithm.It maps all model uncertainties on the parameters, which are described by a multivariate uniform distribution in a parameter hypercube.The output uncertainty is quantified by the 95% prediction uncertainty (95PPU) calculated at the 2.5% and 97.5% levels of the cumulative distribution function of the output variables.Latin hypercube sampling is used to draw independent parameter sets, which lead to the calculation of 95PPU for a given output variable.Two statistics quantify the goodness-of-fit and model output uncertainty.These are the p-factor, which is the percentage of measured data being bracketed by the 95PPU, and the r-factor, which is the average thickness of the 95PPU band.The p-factor has a highest value of 1, while the lowest value for the r-factor is 0. For flow, a value of 0.6-0.8 for the p-factor and a value of around 1 for the r-factor is suggested [56,57].In this definition, (1 − p-factor) can be thought of as the model error, which is the percentage of the measured points not accounted for by the model.
As SUFI-2 is iterative, it does not need a large number of simulations in each iteration [57].We chose 500 simulations in each iteration to balance the requirement of the algorithm with the time needed to run the model.As the computation time for the 1 km 2 delineation was quite long (20 min per simulation), we ran the model 500 times (about one week) to obtain the sensitivity results.
Parameter sensitivity was determined by using a multiple regression system, which regresses parameters against the objective function values derived by the following equation [56,57]: where g is the goal function and b's are the parameters.A t-test is then used to identify the relative significance of each parameter b i .The sensitivities given above are estimates of the average changes in the objective function resulting from changes in each parameter, while all other parameters are changing.This gives relative sensitivities based on a linear approximation and, hence, only provides partial information about the sensitivity of the objective function to model parameters.
To validate the model, we tested the calibrated parameters' ranges over the validation period.We then compared the goodness-of-fit indices from the validation period to the goodness-of-fit indices of the calibration period.As stated in Section 2.3.3,SWAT time-continuous inputs were available from 2005 to 2013.To minimize the uncertainty stemming from initial conditions, we used a four-year warm-up period prior to each simulation (2005)(2006)(2007)(2008).The calibration and the validation periods were a priori chosen as 2009-2011 and 2012-2013, respectively.We then calculated the statistics on annual precipitation and daily discharge to check the possible biases in discharge patterns over these two periods.Both time periods alternate dry and wet years (2009: 530 mm; 2010: 721 mm; 2011: 839 mm; 2012: 532 mm; 2013: 794 mm).The discharge was about 10% higher during the validation period with respect to the calibration period (Table S1).
The 23 parameters (Table 2) were included in the automatic calibration process (500 simulations per iteration, and up to three iterations, until a satisfactory balance was found between p-factor and r-factor).We calibrated and validated the model at the Marquixanes gauging station for three combinations of time-steps and watershed delineations: The rationale behind the latter choices of time-steps and watershed delineation schemes is our focus on the ability of SWAT to simulate discharge at an hourly time-step on a large, highly-responsive catchment.The purpose of the daily time-step simulation was to be able to compare our modelling results with the very large body of literature dealing with SWAT hydrological simulations at a daily time-step, as summarized by Moriasi et al. [59].Therefore, we did not test the case with 1 km 2 minimum drainage area and daily time-step at Marquixanes.Similarly, we did not test the case of the daily time-step at Catllar.Also, as for Catllar, we wanted to check if we could get better results by calibrating only this sub-catchment, without consideration for the others, and using the typical delineation scheme of 15 km 2 minimal drainage area.
To evaluate the goodness-of-fit of the SWAT time-continuous simulations we calculated statistics: the NS coefficient, the coefficient of determination R 2 and the percentage of bias (PBIAS).The statistics were calculated daily for the daily simulation and hourly for three hourly simulations.We then discussed the performance of the SWAT model at both the daily and hourly time-steps in the light of the guidelines given by Moriasi et al. [59] at the daily (R 2 and NS) and monthly (PBIAS) time-steps.We also calculated hourly NS coefficients at flood event scale to compare with the efficiency of the MARINE model when simulating floods.
Table 2. List of the parameters in SWAT changed for sensitivity analysis and for the calibration of discharge at an hourly time-step at the Marquixanes gauging station with a 15 km 2 minimum drainage area watershed delineation: sensitivity ranks (1 being the most sensitive, 23 the least sensitive), initial ranges, and final calibrated ranges.Significantly sensitive parameters (p-values p < 0.001) are marked with * (H 0 : β = 0, reject if p < 0.001).In SWAT-CUP SUFI-2, two types of changes were used: "R" means the existing parameter value is multiplied by (1 + a given value); "V" means the existing parameter value is to be replaced by a given value.The calibration procedure in MARINE consists of estimating, for a given basin, the value of three coefficients of correction applied to distributed parameter maps: one for saturated hydraulic conductivity, named C K , a second for lateral subsurface flow transmissivity, C KSS , and the third for soil thickness, C Z .The Strickler roughness coefficients of the main channel K D1 and of the overbank of the drainage network K D2 are also calibrated.
The total set of eight events from 2006 to 2014 (Table 1) was used for sensitivity analysis and for calibrating and validating the MARINE model at the Marquixanes gauging station.We divided the events into calibration and validation groups of events, based on the results of the sensitivity analysis, as described by Garambois et al. [24].
This methodology is based on a global sensitivity analysis-generalized likelihood uncertainty estimation (GSA-GLUE) approach [60]: we conducted a generalized sensitivity analysis for the five aforementioned model parameters and for all events.We calculated Nash-Sutcliffe efficiencies (NS) for the 5000 Monte Carlo simulations generated from uniform distributions over the variation ranges of the parameters (Table 3).We classified the 250 best simulations as behavioral to obtain a higher threshold value between behavioral and non-behavioral simulations, as recommended by several studies on the uncertainty assessment of hydrological models [61,62].For each event, the cumulative probability distributions were subjected to the Kolmogorov-Smirnov test, to rank the parameters according to the model's sensitivity.Following the same methodology, events for which the model's behavior was consistently atypical were disregarded.We compared, for a given parameter, the cumulative distribution shapes associated with each event as a heuristic for selecting candidate subsets of calibration (and related validation) events.Indeed, previous results [24] showed that events with a similar sensitivity to C Z are likely to present similar average behaviors in terms of rainfall-to-runoff volume conservation.Therefore, an event is discarded from the calibration process if it presents very dissimilar C Z posterior distribution functions with respect to the other events.For a given catchment, large dissimilarities between events for C Z sensitivity and unusual C Z values can be related to questionable quantitative precipitation estimation or to model structure inadequacy (e.g., snowmelt is not simulated).Of course, the events discarded from calibration are kept for validation.
Following the methodology established in previous studies pertaining to the Mediterranean region [22,24], we performed the estimation of the model parameters for the purpose of this study through a local search method, maximizing an NS coefficient among a number of solutions, the NS coefficient being aggregated over a subset of calibration events.In order to avoid the local-optima problem, we repeated local searches five times with different initial conditions.Finally, we evaluated the calibrated parameters against validation events.

Parameters' Sensitivity and Calibration
The sensitive parameters (p-value < 0.001) at the daily time-step for the Marquixanes gauging station (15 km 2 delineation scheme) were precipitation laps rate (PLAPS), deep aquifer percolation fraction (RCHRG_DP), and water level threshold in the shallow aquifer allowing groundwater contribution to the main channel to occur (GWQMN).The best fitted values for PLAPS and RCHRG_DP were 142 mm H 2 O km −1 and 0.46 (-), respectively.The relative small best-fitted value of 1035 mm H 2 O for GWQMN reflects the highly responsive drainage network within the Têt River basin, with a high contribution of the shallow groundwater flow to the main channel flow.
The sensitive parameters at the hourly time-step (15 km 2 delineation scheme, Table 2) were the Manning's roughness coefficient in the main channel (CH_N2), precipitation laps rate (PLAPS), curve number (CN2), deep aquifer percolation fraction (RCHRG_DP) and effective hydraulic conductivity in the main channel alluvium (CH_K2).The CH_N2 parameter defines the rate and the velocity of flow.Hence, its high sensitivity may be related to the torrential bed of the Têt River and of its tributaries upstream the Marquixanes gauging station.Its best fitted value of 0.12 (m −1/3 s) also reflects the number of big stones and brushes within the stream bed.The PLAPS parameter drives the orographic effect on precipitation when defining elevation bands in sub-basins.In our model, some sub-basins have differences of up to 2367 m between their least-elevated area and their most-elevated area.Not surprisingly, the PLAPS with the best fitted value of 138 mm H 2 O km −1 is thus a key factor to adjust precipitation within the river basin.The CN2 parameter drives the infiltration rate within the GAML equation.The best CN2 value was, however, found to be the default value and the parameter was not further adjusted.The RCHRG_DP and the CH_K2 parameters both drive water loss to the groundwater, respectively, along the soil profile and within the stream bed.Their respective best-fitted values of 0.52 (-) and 52 mm h −1 account for the pervious bedrock of the basin (including karst), but they also reflect the water abstraction for irrigation purposes: diverted water is lost for the channel.
At the hourly time-step and 1 km 2 delineation scheme, none of the SWAT parameters were significantly sensitive when simulating discharge.This may indicate that the SWAT hourly 1 km 2 model is not adequate to simulating the system behavior.
At the Catllar gauging station and at the hourly time-step (15 km 2 delineation scheme), the sensitive parameters were the precipitation laps rate (PLAPS), Manning's roughness coefficient in the main channel (CH_N2), deep aquifer percolation fraction (RCHRG_DP), water level threshold in the shallow aquifer allowing groundwater contribution to the main channel to occur (GWQMN) and effective hydraulic conductivity in the main channel alluvium (CH_K2).The fitted value of PLAPS was 27 mm H 2 O km −1 , which reflects the lower impact of elevation amplitude on this sub-catchment (1187-1797 m) located in the eastern part of the whole basin, where the influence of the mountain is smaller.The fitted value of CH_N2 is 0.15 (m −1/3 s) and reflects the roughness of the torrential stream bed.The best-fitted value of 72 mm h −1 for CH_K2 may reflect a more pervious bedrock in the sub-catchment compared to the overall basin.The fitted values of 1265 mm H 2 O for GWQMN and of 0.05 (-) for RCHRG_DP may reflect the highly-responsive drainage network within the sub-catchment.
As previously highlighted by Jeong et al.
[21], the sensitivity of the SWAT parameters was significantly influenced by the model operational time step.The parameters related to channel routing (CH_N2, CH_K2) were relatively more sensitive at the hourly time step, whereas the significance of groundwater flow parameters (RCHRG_DP, GWQMN) was relatively higher at the daily time interval.For instance, at the hourly time step, the sensitivity of the riverbed roughness coefficient CH_N2 is consistent with the fact that this parameter drives the flow velocity and, therefore, the time of peak and the shape of the hydrograph.At the daily time step, the time of peak and the shape of the hydrograph are integrated among other processes and the sensitivity of CH_N2 decreases.
It should be emphasized that our results in this study are conditioned on the data and procedures used for this study.We are aware that input data [63], discretization of the region of study [64,65], regionalization of the parameters [66], method of calibration and the choice of objective function [67] all affect final parameter ranges and their sensitivities.

Discharge Calibration and Validation Results in the Main Channel
Overall, the agreement between observations and simulations (Figure 2) is consistent between the calibration (2009-2011) and the validation (2012-2013) periods (Table 4), suggesting that the choice of the calibration and validation periods was relevant and that the model is robust regardless of the gauging station, the time-step, and the delineation scheme.
and the shape of the hydrograph are integrated among other processes and the sensitivity of CH_N2 decreases.
It should be emphasized that our results in this study are conditioned on the data and procedures used for this study.We are aware that input data [63], discretization of the region of study [64,65], regionalization of the parameters [66], method of calibration and the choice of objective function [67] all affect final parameter ranges and their sensitivities.

Discharge Calibration and Validation Results in the Main Channel
Overall, the agreement between observations and simulations (Figure 2) is consistent between the calibration (2009-2011) and the validation (2012-2013) periods (Table 4), suggesting that the choice of the calibration and validation periods was relevant and that the model is robust regardless of the gauging station, the time-step, and the delineation scheme.Considering the 15 km 2 delineation scheme at Marquixanes, the best values of NS, R 2 and PBIAS were obtained at the daily time-step, which is consistent with the fact that the daily outputs integrate the variability of smaller time-steps.Considering the guidelines suggested by Moriasi et al. [59] at the daily (R 2 and NS) and monthly (PBIAS) time-steps, the daily time-step R 2 , NS and PBIAS values were deemed satisfactory for both the calibration and validation periods.Similarly, the hourly time-step R 2 , NS and PBIAS values were also deemed very satisfactory for both the calibration and validation periods.
The hourly time-step NS and PBIAS values for the 1 km 2 delineation scheme at Marquixanes resulted in poor performance, although the PBIAS was small (Table 4).Discharge was much noisier with the 1 km 2 delineation scheme than with the 15 km 2 delineation scheme, although discharge peaks were correctly captured (Figure 2).
Finally, the hourly time-step NS and PBIAS values for the 15 km 2 delineation scheme at Catllar resulted in intermediate performance, suggesting that the rain gauge network and rainfall spatialization may be limiting the simulation quality in this small and non-anthropogenized sub-catchment.Also, small catchments such as this 90 km 2 catchment may not fully represent the variability of the spatial inputs (e.g., soil, land use, rainfall) of the 810 km 2 large catchment.
The good agreement between observations and simulations is also shown by the flow duration curves (Figure 3).SWAT simulated a lower flow magnitude during baseflow (>80% exceedance) at Marquixanes at the daily time-step.For both the 15 and the 1 km 2 delineation schemes at Marquixanes, the hourly time-step SWAT simulated a higher flow magnitude for both extreme high flows and extreme low flows.At Catllar, SWAT simulated a higher flow magnitude during baseflow (>80% exceedance, considering that the river is intermittent and that the lower limit of discharge record at Catllar station was 0.01 m 3 s −1 ).
The improvement in the model's performance at the sub-daily time-step is mostly due to the improvement in predicting high flows [21].Indeed, in our study, the peak discharge is better simulated at the hourly time-step with the 15 km 2 delineation scheme.Decreasing the minimum drainage area did not improve the simulation [68,69] but introduced noise instead.The channel routing scheme described in SWAT may be too simplified for the finest spatial resolution, i.e., the 1 km 2 delineation scheme, and thus contribute to the noise in the outputs.Also, in the Têt basin, the number of available rain gauges and their spatialization, together with the resolution of the input maps, may be factors limiting the improvement of the model performances when decreasing the minimum drainage area.

Parameters' Sensitivity and Calibration
As stated in previous studies [25], the sensitivity ranking of parameters varied as a function of the event.However, lateral subsurface flow transmissivity CKSS and soil thickness CZ came up as the two most sensitive parameters for the Têt catchment.These two parameters control the runoff production simulated by the model (CZ) and the subsurface water dynamics (CKSS).The distribution functions of behavioral simulations for these two parameters is illustrated in Figure 4.As explained in Section 2.4.2, the 250 best simulations have been classified as behavioral, i.e., representative of the observed behavior of the system as defined by [70].

Parameters' Sensitivity and Calibration
As stated in previous studies [25], the sensitivity ranking of parameters varied as a function of the event.However, lateral subsurface flow transmissivity C KSS and soil thickness C Z came up as the two most sensitive parameters for the Têt catchment.These two parameters control the runoff production simulated by the model (C Z ) and the subsurface water dynamics (C KSS ).The distribution functions of behavioral simulations for these two parameters is illustrated in Figure 4.As explained in Section 2.4.2, the 250 best simulations have been classified as behavioral, i.e., representative of the observed behavior of the system as defined by [70].Concerning the winter event 20060127, the shape of the PDF of CZ was similar to the shape of the autumn events, whereas the shape of the PDF of CKSS was similar to the shape of the spring events.
Under our methodology, large PDF dissimilarities between events can be attributed to questionable quantitative precipitation estimation.This is in part corroborated by the dissimilarities in the PDF of the winter event: indeed for 20060127, part of the precipitation is likely to have been in the form of snow, which is well known to be difficult to measure.
As a result, the fact that the model behaves differently during spring or autumn events led us to choose a calibration strategy that mixes both spring and autumn events (Table 5).We thus tested several combinations of two spring and two autumn events while excluding the two events exhibiting behaviors highly dissimilar to the average spring and autumn behaviors.The resulting calibration strategy is highlighted is Table 5.
Table 5. Goodness-of-fit indices (Nash-Sutcliffe efficiencies) for the six flood events simulated with both the MARINE and the SWAT models at Marquixanes gauging station.Numbers in bold highlight Nash-Sutcliffe efficiencies higher than 0.5.Events marked with * are calibration events for MARINE.5).As expected, the two events excluded from the calibration process gave low performances in validation.Surprisingly, another event (20130428) also gave poor performance in validation with an NS coefficient of 0.38.This may be due to a high snow melt contribution during this event, which occurred in the spring, as the snow melt process is not described in the MARINE model.This is corroborated by the simulated The shape of the PDF of C KSS was similar to the shape of the spring events.Concerning the winter event 20060127, the shape of the PDF of C Z was similar to the shape of the autumn events, whereas the shape of the PDF of C KSS was similar to the shape of the spring events.
Under our methodology, large PDF dissimilarities between events can be attributed to questionable quantitative precipitation estimation.This is in part corroborated by the dissimilarities in the PDF of the winter event: indeed for 20060127, part of the precipitation is likely to have been in the form of snow, which is well known to be difficult to measure.
As a result, the fact that the model behaves differently during spring or autumn events led us to choose a calibration strategy that mixes both spring and autumn events (Table 5).We thus tested several combinations of two spring and two autumn events while excluding the two events exhibiting behaviors highly dissimilar to the average spring and autumn behaviors.The resulting calibration strategy is highlighted is Table 5.
Table 5. Goodness-of-fit indices (Nash-Sutcliffe efficiencies) for the six flood events simulated with both the MARINE and the SWAT models at Marquixanes gauging station.Numbers in bold highlight Nash-Sutcliffe efficiencies higher than 0.5.Events marked with * are calibration events for MARINE.5).As expected, the two events excluded from the calibration process gave low performances in validation.Surprisingly, Water 2017, 9, 929 another event (20130428) also gave poor performance in validation with an NS coefficient of 0.38.This may be due to a high snow melt contribution during this event, which occurred in the spring, as the snow melt process is not described in the MARINE model.This is corroborated by the simulated discharges (Figure 5e), where MARINE significantly underestimated the flow peak.A complementary explanation can be that the input data exhibited a lower spatial variability than the data traditionally used as input for the MARINE model.Yet it has already been shown that the spatial variability has a great impact on the simulation of the catchment response [25].Fitted parameter values are listed in Table 3.
Water 2017, 9, 929 17 of 25 discharges (Figure 5e), where MARINE significantly underestimated the flow peak.A complementary explanation can be that the input data exhibited a lower spatial variability than the data traditionally used as input for the MARINE model.Yet it has already been shown that the spatial variability has a great impact on the simulation of the catchment response [25].Fitted parameter values are listed in Table 3.

Discharge Simulation
Figure 5 compares observed discharges with both SWAT (15 km 2 delineation scheme, hourly simulation) and MARINE simulated discharges at the Marquixanes gauging station for the six selected flood events.The performance of the MARINE model were better than the SWAT model, considering NS value (Table 5).For these same six events of interest, MARINE was slightly more accurate in simulating the magnitude of the flood peaks (see events 20101007, 20110312, 20130303 and 20131116) and the timing of the simulated peak discharge (20101007, 20110312, 20130303 and 20130428), whereas SWAT was slightly more accurate in simulating the overall exported volume of water during a given event (20101007, 20110312, 20111103 and 20131116) (Table 1).

Discharge Simulation
Figure 5 compares observed discharges with both SWAT (15 km 2 delineation scheme, hourly simulation) and MARINE simulated discharges at the Marquixanes gauging station for the six selected flood events.The performance of the MARINE model were better than the SWAT model, considering NS value (Table 5).For these same six events of interest, MARINE was slightly more accurate in simulating the magnitude of the flood peaks (see events 20101007, 20110312, 20130303 and 20131116) and the timing of the simulated peak discharge (20101007, 20110312, 20130303 and 20130428), whereas SWAT was slightly more accurate in simulating the overall exported volume of water during a given event (20101007, 20110312, 20111103 and 20131116) (Table 1).
Several previous studies have compared the performances of lumped models with the performances of fully-distributed models when modelling flood events [71,72].They found that the distributed modelling approaches were often better than the lumped modelling approaches, which is consistent with the spatially uniform assumptions that underpin the model structure [73].In this study, the MARINE model performed better than the SWAT model in terms of NS coefficients, but it is likely that the NS coefficients of SWAT would have been higher if the model had been calibrated and validated over flood periods only.
However, the floods that were poorly simulated in SWAT also performed poorly in MARINE.This may be explained by the spatial and temporal accuracy of the rainfall measurements, as previously mentioned (Sections 3.1 and 3.2) and discussed by e.g., Grusson et al. [74].Indeed, quantitative precipitation estimation over a large area is a difficult task.Garambois et al. [24] emphasized that in the case of rainfall-runoff modelling, a key issue is that the rain gauge network may have not been adequately designed to accurately estimate the overall spatial pattern of the rainfall amount, location and timing, and thus to accurately explain the hydrological response of the catchment.For instance, the sensitivity analysis of the MARINE model highlighted the questionable precipitation estimation for the flood of November 2013.In the case of SWAT at the daily time-step, and to a lesser extent the hourly time-step, it is also worth noting that autumn floods were better simulated than spring floods.This may be due to the limited snow-melt description in SWAT [75], possibly affecting spring discharge.However, the uncertainty 95PPU bands of both models mostly enveloped the observed peak discharge.The SWAT simulations, including 95PPU bands, were noisier than MARINE's, but SWAT better simulated the flood recessions.For each flood, simulated peak discharge, peak timing, exported water volume and recession discharge depend on the soil water content at the beginning and the end of the rainfall event.Soil water content in both modelling approaches are discussed in the next section.

Soil Water Content
For both the SWAT and SIM-calculated initial soil water content for MARINE, the average soil water saturation before the 6 floods of interest was about 40-60% (Figure 6).In SWAT, the northern edge of the basin had a lower saturation related to the higher soil hydraulic conductivity in these areas (Figure 1e).The spatial distribution of the saturation states was quite different between SWAT and SIM before the floods and between SWAT and MARINE after the floods.MARINE simulation resulted in a higher saturation at the end of the events (Figure 6), especially in the middle upper part of the catchment where the soil depth is low (100 mm) as is the hydraulic conductivity (6 mm h −1 ). Figure 6 also shows the role of the drainage network in the evolution of the saturation state for the MARINE model.The differences between the models are of course related to the description of the infiltration and the subsurface flow dynamics, but also to the spatial resolution of each model.As the flood events last several days after the rainfall peak (Figure 5) and the spatial resolution of the rainfall data is quite coarse (only nine rain gauges in this 810 km 2 area), this spatial resolution does not seem to have a major impact at the end of the simulations.All these results emphasize the significance of the impact of the soil properties on simulated soil saturation dynamics, as already mentioned in previous studies [76,77].Measures describing the spatial distribution of the saturation state would be valuable to evaluate the performance of the models.

Sensitive Parameters
In our study, the sensitive parameters were different depending on the modelling approach (Tables 2 and 3).The most sensitive parameters in MARINE were soil depth and soil lateral

Sensitive Parameters
In our study, the sensitive parameters were different depending on the modelling approach (Tables 2 and 3).The most sensitive parameters in MARINE were soil depth and soil lateral transmissivity.In SWAT, the soil depth was not a calibrated parameter.In comparison to the MARINE parameters, the effective hydraulic conductivity in the main channel alluvium was a highly sensitive parameter in SWAT.Finally, the Manning's roughness coefficient in the main channel was one of the second most sensitive parameters in SWAT, whereas they were the least sensitive in MARINE.To sum up, the most sensitive parameters in MARINE were the parameters that control the runoff production amount (soil depth) and the subsurface water dynamics (lateral transmissivity) whereas the most sensitive parameters in SWAT were the parameters that control the stream-to-aquifer flow (hydraulic conductivity in the main channel alluvium) and the runoff dynamics (roughness coefficient).
The difference in the rank of the sensitive parameters may be explained by the processes respectively described in both models.In addition, the difference in the sensitivity ranks can also be explained by the calibration procedures for both models: SWAT was calibrated over a continuous three-year period of time (2009-2011), thus including both low flow and high flow periods, whereas MARINE was calibrated on four flood events.Also, for both models, the calibration procedure induced a compensation effect, e.g., in SWAT the high values of the effective hydraulic conductivity in the main channel alluvium may partly compensate for the flow diverted for irrigation purposes.Higher infiltration in the channel could also compensate for a too low water infiltration capacity in the catchment.Indeed, the soil depth was not calibrated in SWAT, but the calibration factor C Z in MARINE has a fitted value of 5.1 (Table 3), which means that MARINE needs to enhance the catchment storage capacity to be able to correctly simulate discharge.The input data and especially their low spatial resolution are actually likely to have an impact on the results of the sensitivity analysis.Previous studies using soil data with a higher spatial resolution confirmed that the MARINE model is very sensitive to the soil depth and lateral subsurface transfer, which explains more than 80% of model output variance when the hydrographs are peaking and during the slow-declining limb, respectively [26].However, they also showed an important sensitivity to the catchment soil infiltration capacity (hydraulic conductivity, [24]) that does not appear in the present study.

Conclusions
In this study, we tested the sub-daily module of the SWAT model when applied to a catchment of ~1000 km 2 .To our knowledge, this is the first time that SWAT has been run at a sub-daily time step over such a catchment and calibrated and validated with sub-daily-measured discharge data.We compared simulated discharges with observed discharges at two gauging stations (Marquixanes and Catllar).We compared the time-continuous hourly discharges simulated with SWAT (with two sub-basin delineation schemes) to the event-scale, sub-daily discharges simulated with the MARINE model.We also compared the soil water contents before the flood and after the flood.To our knowledge, this study is also the first comparing the performance of the conceptual SWAT model at the sub-daily time-step with a process-oriented model, namely the MARINE model, over a large river catchment.
In this study, we showed the following: • During flash floods, the fully-distributed MARINE model better simulated peak discharge and peak timing, whereas the lumped SWAT model better simulated recession discharge and exported water volume; • Overall, the latter suggests that both SWAT and MARINE performed equally satisfactorily during floods, especially considering that SWAT was calibrated and validated time-continuously over both low and high flow periods, whereas MARINE calibration and validation was focused on floods only.This also suggests that the simulation performance of the two modelling approaches was not limited by the described processes and by the different calibration procedures, but rather by the rainfall and soil data; • Due to the limitations in rainfall and soil data, the performance of the SWAT model at the hourly time-step were not necessarily improved by the decreasing in the size of the minimum drainage area when delineating sub-basins.They were also not improved when separately calibrating a smaller sub-catchment.
Hence, the SWAT model appears to be a reliable modelling tool to predict discharge over long periods of time in large flash-flood-prone basins.The correct simulation of exported water volumes by SWAT during floods makes it a promising tool for future investigations of more hydrological processes simulated at the sub-daily time-step, such as suspended sediment transport.

Supplementary Materials:
The following are available online at www.mdpi.com/2073-4441/9/12/929/s1, Figure S1: Sub-basins delineations within the Têt River basin with (a) 1 km 2 and (b) 15 km 2 minimum drainage areas for SWAT simulations, Table S1: Statistics calculated on measured daily discharge records at Marquixanes gauging station for the SWAT model calibration and the validation periods.

•
Model delineated with 15 km 2 minimum drainage area and simulation at a daily time-step • Model delineated with 15 km 2 minimum drainage area and simulation at an hourly time-step • Model delineated with 1 km 2 minimum drainage area and simulation at an hourly time-step We also performed additional calibration (2009-2011) and validation (2012-2013) of the SWAT model at hourly time-step and with 15 km 2 minimum drainage area at the Catllar gauging station (500 simulations).

Figure 2 .
Figure 2. Observed and SWAT simulated discharges (m 3 s −1 ) with 95% confidence intervals (95PPU band) for both calibration (2009-2011) and validation (2012-2013) periods: (a) Simulated and observed discharge at Marquixanes at the daily time-step, using a minimum drainage area of 15 km 2 ; (b) Simulated and observed discharge at Marquixanes at the hourly time-step, using a minimum drainage area of 15 km 2 ; (c) Simulated and observed discharge at Marquixanes at the hourly time-step, using a minimum drainage area of 1 km 2 ; (d) Simulated and observed discharge at Catllar at the hourly timestep, using a minimum drainage area of 15 km 2 .For better readability, the Y axes of the hourly figures were limited above the maximal observed hourly discharge, i.e., 161 m 3 s −1 in Marquixanes and 50 m 3 s −1 in Catllar, respectively.

Figure 2 .
Figure 2. Observed and SWAT simulated discharges (m 3 s −1 ) with 95% confidence intervals (95PPU band) for both calibration (2009-2011) and validation (2012-2013) periods: (a) Simulated and observed discharge at Marquixanes at the daily time-step, using a minimum drainage area of 15 km 2 ; (b) Simulated and observed discharge at Marquixanes at the hourly time-step, using a minimum drainage area of 15 km 2 ; (c) Simulated and observed discharge at Marquixanes at the hourly time-step, using a minimum drainage area of 1 km 2 ; (d) Simulated and observed discharge at Catllar at the hourly time-step, using a minimum drainage area of 15 km 2 .For better readability, the Y axes of the hourly figures were limited above the maximal observed hourly discharge, i.e., 161 m 3 s −1 in Marquixanes and 50 m 3 s −1 in Catllar, respectively.

Figure 3 .
Figure 3. Observed and SWAT-simulated flow duration curves (m 3 s −1 ): (a) Simulated and observed discharge at Marquixanes at the daily time-step, using a minimum drainage area of 15 km 2 ; (b) Simulated and observed discharge at Marquixanes at the hourly time-step, using a minimum drainage area of 15 km 2 ; (c) Simulated and observed discharge at Marquixanes at the hourly timestep, using a minimum drainage area of 1 km 2 ; (d) Simulated and observed discharge at Catllar at the hourly time-step, using a minimum drainage area of 15 km 2 .

Figure 3 .
Figure 3. Observed and SWAT-simulated flow duration curves (m 3 s −1 ): (a) Simulated and observed discharge at Marquixanes at the daily time-step, using a minimum drainage area of 15 km 2 ; (b) Simulated and observed discharge at Marquixanes at the hourly time-step, using a minimum drainage area of 15 km 2 ; (c) Simulated and observed discharge at Marquixanes at the hourly time-step, using a minimum drainage area of 1 km 2 ; (d) Simulated and observed discharge at Catllar at the hourly time-step, using a minimum drainage area of 15 km 2 .

Figure 4 .
Figure 4. Distribution functions of behavioral simulations per event for (a) the CKSS (-) parameter and (b) the CZ (-) parameter.

Figure 4 .
Figure 4. Distribution functions of behavioral simulations per event for (a) the C KSS (-) parameter and (b) the C Z (-) parameter.

Figure 5 .
Figure 5. Sub-daily discharges (m 3 s −1 ) simulated with both the SWAT (best simulation at hourly timestep) and MARINE (multi-event calibration at variable time-step) models and observed data for the six flash flood events that occurred in the Têt River basin (Marquixanes gauging station) from 2009 to 2013: (a) event from 7 October 7 to 27 October 2010; (b) event from 12 March to 20 March 2011; (c) event from 3 November to 15 November 2011; (d) event from 3 March to 15 March 2013; (e) event from 28 April to 7 May 2013; (f) event from 16 November to 23 November 2013.Maximal 95PPU upper limit is 806 m 3 s −1 .

Figure 5 .
Figure 5. Sub-daily discharges (m 3 s −1 ) simulated with both the SWAT (best simulation at hourly time-step) and MARINE (multi-event calibration at variable time-step) models and observed data for the six flash flood events that occurred in the Têt River basin (Marquixanes gauging station) from 2009 to 2013: (a) event from 7 October 7 to 27 October 2010; (b) event from 12 March to 20 March 2011; (c) event from 3 November to 15 November 2011; (d) event from 3 March to 15 March 2013; (e) event from 28 April to 7 May 2013; (f) event from 16 November to 23 November 2013.Maximal 95PPU upper limit is 806 m 3 s −1 .

Figure 6 .
Figure 6.Soil water content (%) before and after the floods, and cumulative rainfall (mm) during the six major flood events in the Têt River basin from 2009 to 2013.Initial soil water content (before flood) was simulated with SWAT and SIM and final soil water content (after flood) was simulated with SWAT and MARINE.Cumulative rainfall was spatialized from rain gauge records with Thiessen polygons.

Figure 6 .
Figure 6.Soil water content (%) before and after the floods, and cumulative rainfall (mm) during the six major flood events in the Têt River basin from 2009 to 2013.Initial soil water content (before flood) was simulated with SWAT and SIM and final soil water content (after flood) was simulated with SWAT and MARINE.Cumulative rainfall was spatialized from rain gauge records with Thiessen polygons.

Table 1 .
List of the flood events of interest for MARINE model calibration and validation from 2006 to 2014 at Marquixanes gauging station, and for comparison with the Soil and Water Assessment Tool (SWAT) at flood event-scale (hourly outputs, 15 km 2 delineation scheme): peak discharge, time of peak and water exported volumes.
2.3.3.Data Selection for Models' Set up and Comparison 51] • Land use maps were derived for both models from the 100 m 2006 Corine Land Cover map (Source: European Environmental Agency) • A common map of soil depth and soil texture was computed based on the Food and Agriculture Organization (FAO) Digital Soil Map of the World (Source: Land and Water Development Division, FAO, Rome), with associated soil properties from the French National Institute of Agronomic Research (INRA).MARINE assumes homogeneous soil profiles and estimates soil classes according to pedotransfer function [52], each soil class being associated to saturated water content, suction and saturated hydraulic conductivity.

Table 3 .
List of the catchment parameters and parameter ranges in MARINE.Sensitivity ranks are computed according to the maximum distance between the behavioral and non-behavioral distributions (1 is the most sensitive, 5 is the least sensitive).Fitted values are obtained with the best subset of calibration events.

Table 4 .
[59]ness-of-fit indices for calibration(2009)(2010)(2011)and validation (2012-2013) of the SWAT model at both the hourly and daily time-steps, for two minimum drainage area watershed delineations, and at two gauging stations.Numbers in bold highlight indices higher than the "satisfactory" thresholds suggested by Moriasi et al.[59]at daily (R 2 and NS) and monthly (PBIAS) time-steps.