Comparison of Regionalisation Techniques for Peak Streamflow Estimation in Small Catchments in the Pilbara, Australia

: Arid and semi-arid regions typically lack high-resolution river gauging data causing dif-ficulties in understanding rainfall-runoff patterns. A common predictive method for discharge estimation within ungauged catchments is regional flood frequency estimation (RFFE), deriving peak discharge estimates from similar, gauged catchments and applying them to the catchment of interest. The majority of RFFE equations are developed for larger catchments where flow events may be larger and of greater interest. We test a series of RFFE methods derived for the Pilbara region, applying them to new ungauged small catchments under 10km 2 . Rainfall values are derived from a guideline Australian design rainfall database, Australian Rainfall and Runoff 2019 (ARR2019) which was recently updated with an additional 30 years of rainfall data. RFFE equations are compared to a direct rainfall model to evaluate their performance within small catchments, identifying key limitations and considerations when modelling small headwater catchments.


Introduction
Flood frequency analysis is commonly undertaken to identify and estimate flood quantiles corresponding to a given return period using the available streamflow observations in a catchment [1]. Where streamflow records are absent, or catchments are widely ungauged, prediction of streamflow involves applying regional flood frequency estimation (RFFE) methods, which are data-driven empirical procedures that attempt to compensate for a lack of temporal data at a given location with spatial data, obtained from other locations within a homogenous region [2]. RFFE approaches are frequently used throughout Australia for the prediction of flood events. Accurately predicting the magnitude of flood events is essential for the planning of water resource systems [3] in addition to adhering to the design standards of engineered structures designed to convey a certain flow [4,5]. A lack of long-term rainfall and streamflow data within arid and semi-arid regions is a major issue for run-off modelling [6] and this challenge is highlighted here within the Pilbara region in Western Australia.
Western Australia accounts for 39% of the global supply of iron ore [7] but has a notable lack of rainfall and streamflow data owing to (a) the high spatial and temporal variability of streamflow, (b) the high cost of establishing dense hydrologic gauging sites (c) the likelihood of gauging sites being destroyed by flash flooding events and (d) the disproportionate interest of flow events in larger river channels resulting in a scarcity of data smaller catchments. This data scarcity results in a fewer opportunities for validation of hydrological models used and therefore there are limited opportunities to demonstrate Citation: Flatley, A.; Rutherfurd, I.

Comparison of Regionalisation
Techniques for Peak Streamflow Estimation in Small Catchments in the Pilbara, Australia. Hydrology that a given-site specific model is capable of making accurate predictions for periods outside a calibration period [8]. Nevertheless, increasing engineering modifications are carried out in small catchments in the Pilbara, such as culvert installation for roads and railways or mining river diversions-such as those found in our study region. River diversions are expensive engineering structures and their design and performance rely on adequately estimating both frequent and rare flow events [9]. There is also increased recognition of the importance to estimate and subsequently preserve natural hydrologic flow conditions resulting in an interest in estimating frequent flow events (e.g., 1EY, 0.5EY and 0.2EY where EY stands for the average number of exceedances per year) particularly in these semi-arid areas where flow is sparse and poorly understood.
Current practices for the design and prediction of peak discharges for river relocation designs are based on several estimation methods such as regionalisation methods, numerical and hydrologic modelling. The outcome of these methodological provides wide-ranging, typically heuristic results. However, the improvement of two-dimensional hydrodynamic models is changing how peak discharges are being estimated. Direct rainfall modelling (or rain-on-grid) has increasingly become a standard approach for predicting design flood behaviour [10,11]. The application of rainfall directly onto a 2D domain allows for the simulation of catchment runoff by applying it directly to the modelling grid [12,13]. Catchment runoff within direct rainfall models is dependent on the grid or mesh cell area, the rainfall depth, grid cell roughness, rainfall losses (IL/CL) and the slope between neighbouring cells. This distributed rainfall approach directly onto the 2D domain can give considerably more detail, particularly in the upper parts of sub-catchments [14] and has been found to provide a better representation of minor overland flowpaths than conventional modelling methods [12]. In the absence of stream gauges, this paper examines the suitability of direct rainfall modelling to test the suitability of RFFE approaches on small ungauged catchments within the Pilbara. The tested RFFE approaches were previously designed for the semi-arid Pilbara region and have previously been validated using gauged flow data within medium to large catchments [15][16][17][18]. Most direct-rainfall models are used as a means of indirect calibration or comparison with traditional hydrological predictive methods and can be used to elucidate discrepancies in other models [19]. Direct rainfall models are useful in the modelling of design flood simulations providing appropriate checks and quality assurance procedures are implemented [13]. There are many sources of uncertainty that can have a strong influence on flood mapping and flow hydrographs, including synthetic rainfall estimation with IDF curves [20], initial soil moisture conditions [21,22] the basin response model, modelling grid sizes [23] and the difference between storm return time and the correspondent flood return time [24] to name a few. This paper tests the suitability of these RFFE approaches within smaller headwater catchments in the Pilbara and evaluates them through a comparison of peak discharge values derived from a 2D hydrodynamic direct rainfall model. This paper provides the first comparative study of RFFE approaches for the Pilbara using updated ARR (2019) [25] values to validate their use within smaller catchments in the same region. To achieve this, we firstly (a) select existing RFFE techniques previously applied and validated in larger catchments, (b) simulate direct rainfall events within a series of smaller catchments (c) use the resulting discharge values to quasi-validate their application within small headwater catchments to provide a range of likely predicted peak discharges for a range of annual exceedance probabilities.

Study Area
The Pilbara region of Western Australia is a large arid to semi-arid region with a wide expanse of ungauged catchments with limited streamflow records within its catchments [26,27]. The region in a transitional location between the Eyrean (central desert) and the southern Torresian (tropical) bioclimatic regions [28]. The Pilbara region is classified as either Arid, desert hot (Bwh) and Arid steppe hot (Bsh) by the updated Köppen-Geiger climate classification [29]. Temperatures exceed 30 °C for most of the year, and rarely dip below 10 °C. Because of its positioning, tropical depressions and recurrent cyclonic events comprise most the regions' total rainfall [30,31]. The Pilbara has a low total rainfall averaging between 250-300 mm annually [26,32,33] but the majority of flow events are concentrated in short duration floods of high magnitude [34]. The region is characterized by extreme hydroclimatic conditions, in which the rainfall is highly sporadic [31], driven by infrequent tropical cyclones ( Figure 1) and thunderstorms occurring mainly within the summer months between January to March [26,30,35]. Winter rainfall is typically from low pressure trough systems [26]. For smaller streams, such as those within headwater channels, flow events typically last around 5 days of the year. Rainfall is very localised, causing issues for the correlation of rainfall and runoff. In addition, the Pilbara has very high evaporative losses, the annual potential evaporation is also 10 times higher than the annual rainfall [37]. Runoff is also highly variable, and only 2-13% of mean annual rainfall becomes runoff in the Pilbara [38]. Higher percentages of rainfall almost certainly run off in small headwater catchments, with a lot of runoff infiltrating into streambeds and therefore failing to reach downstream gauging stations [38]. Most gauging sites are located within larger catchments, however most of these initial stream gauges were not operational until 1967 [33,39].
The Pilbara region is also rich in iron ore and open-cut mines resulting in a wide array of engineering structures built to service mine sites, in addition to many watercourse modifications such as culverts, drains and river diversion channels constructed both within large and smaller channels. River diversion channels for mining in the Pilbara region are designed conservatively to convey rare flow events driven by large cyclonic events or infrequent 100 and 1000-year ARI floods (or the 1 to 0.1AEP (%)) [9]. Many river diversions are constructed within smaller catchments that lack gauged rainfall or streamflow data resulting in a poor understanding of the peak flood discharges and more frequent events experienced within these catchments.

General Approach
This research tests a series of RFFE methods to calculate predicted peak discharge (Qpeak) previously used within the Pilbara and applies them specifically to smaller sized catchments using updated ARR2019 IFD rainfall values. The selected RFFE methods were chosen based on their satisfactory performance when applied in larger catchments in the Pilbara. The methods tested include a QRT and PRT method from Taylor et al., 2011 [15], QRT and PRT methods from Rahman et al., 2012 [2], a Regional Flood Frequency Procedure (RFFP2000) from Flavell., 2012 [17], a IFM method from Davies and Yip, 2014 [18] and the ARR 2016 RFFE (Table 1). Further details on the RFFE approaches and their prior performance is provided in Table 1 of the Supplementary Material. Next, the RFFE methods are quasi-validated against a TUFLOW direct rainfall model, acting as the observation of rainfall within the catchments lacking gauging infrastructure. The results yield a range of peak discharge estimates for each AEP. The direct rainfall model (and associated sensitivity analysis) is used to provide validation for the most appropriate RFFE for small catchments. Figure 2 shows the full sequence of steps. Similar approaches have been undertaken to assess the performance of ARR 2016 RFFE using RORB modelling [40] or to incorporate it into a direct rainfall model of complex urban catchments [41]. However, few studies have used direct rainfall modelling approaches to benchmark RFFE methods within small semi-arid headwater channels.  With the largest value from two Q20 equations being adopted for the Q20 value: Q20 = 1.98 × 10 − 23(ASe 0.5 ) 0.79 LAT −15.08 LONG 20.91 Q20 = Q10 = (13.21A 0.61 )/(8.74A 0.60 ) Q50 = Q20 × frequency factor (Q50/Q20) Q100 = Q20 × frequency factor (Q100/Q20) A = catchment area (km 2 ), Se = equivalent uniform slope (m/km) and L = mainstream length (km) Figure 2. Workflow of steps to determine Regional Flood Frequency Estimation (RFFE) selection and for small headwater channels.

Regionalisation Approaches
Regionalisation refers to the process of transferring hydrological information from one catchment to another. RFFE approaches assume a statistical relationship between observable catchment properties and flood discharge characteristics, allowing the construction of flood hydrographs by applying relationships developed for gauged catchments with similar properties [42]. Regionalization approaches are commonly based either on spatial proximity or similar catchment attributes [43]. Within regionalization methods, model parameters are used as an instrument to transfer hydrological information from gauged to ungauged basins [44]. In general, a regional model can be stated in a simplified form defined by Wagener and Wheater (2006) as: where is the estimated hydrological variable of interest at the ungauged site (it can be an estimated model parameter, probability or cumulative distribution function parameter, or hydrological response such as streamflow or flow events), is a functional relation for using a set of catchment attributes-physiographic or meteorological attributes , is a set of regional hydrological variables of interests and is an error term [45]. Regionalisation approaches may be satisfactory if the catchments are similar in some sense, but error prone if they are not [46]. Razavi and Coulibaly (2013) provide a review of methodology for streamflow prediction in ungauged basin using regionalization methods, concluding that most model-dependent methods in arid to warm-temperate climates (e.g., Australia) indicate that physical similarity and spatial proximity appear to be the best approach to estimating streamflow. However, the most regionalisation methods are highly site-specific, and therefore, a comparative study between suitable approaches is suggested before selecting the regionalisation method for a given site or region.

Regional Flood Frequency Estimation (RFFE)
Regional flood frequency estimation (RFFE) is widely used to estimate flood quartiles in ungauged catchments. RFFE approaches provide an alternative method to flood frequency analysis (FFA) where a lack of temporal data is substituted with spatial data to make more accurate flood estimates at ungauged sites [15]. Common RFFE techniques include Probabilistic Rational Method (PRM), Quantile Regression Techniques (QRT) the Index Flood Method (IFM) and a Parameter Regression Technique (PRT). Regression based RFFE methods are more commonly applied to recent studies within Australia. The following section details selected regression based RFFE methods, which are developed from a longer record of data and are considered to give a more reliable estimation of design flows.

Quartile Regression Technique
The quartile regression technique (QRT) is used frequently within ungauged catchments. The method estimates flood quartiles through a multiple regression between recorded streamflow data and a set of climatic and catchment characteristics within a region [47]. The QRT regression technique is expressed as: where B, C and D are catchment and climatic characteristics variables (predictors); a, b, c, d are the regression coefficients and is the flood magnitude with T-year ARI (flood quantile) [48].

Parameter Regression Technique
The parameter regression technique (PRT) is similar to the QRT. However, instead of quartiles, the first three moments of the log-Pearson type 3 (LP3) distribution are taken as dependent variables in regression analysis against catchment characteristics [15]. Let Q be the annual maximum flood series at a site and X = ln(Q), then the mean (M), standard deviation (S) and skew (g) of the X series are taken as dependent variables: where is a flood quantile of T years ARI and is the standardized LP3 frequency factor (which is a function of skew) and can be obtained from ARR or can be approximated [15].

Index Flood Method
The index flood method (IFM) assumes that the exceedance probability distribution of annual peak discharge is identical, except for a site-specific scaling factor called the index flood (average likely flood) [49]. The IFM method is expressed as: where is the flood quantile, is the function basin area, slope and is a regional growth factor (a dimensionless frequency distribution quantity common to all sites within each homogeneous region).

Australian Rainfall and Runoff (ARR) RFFE
The Australian Rainfall and Runoff (ARR) guidelines offer predicted estimates of rainfall intensity, frequency, and duration (IFD) values for Australia. Additionally, ARR also has a regional flood frequency estimation model which is widely used and recommended for design flood estimation [11,50]. The ARR guidelines were updated in 2016 and again in 2019. Five predictor values were adopted for the RFFE technique [51,52]. These predictor values are: catchment area (in km²); design rainfall intensity at catchment centroid (in mm/h) for the 6 h duration and 50% AEP ( 50% I6h); design rainfall intensity at catchment centroid (in mm/h) for the 6 h duration and 2% AEP ( 2% I6h); ratio of design rainfall intensities of AEPs of 2% and 50% for duration of 6 h ( 2% I6h/ 50% I6h); and catchment shape factor (Sf), which is defined as the shortest distance between catchment outlet and centroid divided by the square root of catchment area. The RFFE technique used in ARR is adapted for different regions throughout Australia (Figure 3). Each region is determined based on a Region of Interest (ROI) approach based on geographical proximity of gauging stations, with fringe zones between regions defined by the 500 mm and 400 mm isohyet to delineate between humid and arid/semi-arid regions [53]. The Pilbara region was characterized as an alternative sub-region distinct from the other arid and semi-arid regions of Australia [51]. This was due to (a) concentrations of stream gauging stations in three parts of Western Australia which are separated by long distances (e.g., Kimberley region, Pilbara region and South-West WA) and (b) notable differences in region hydrologic conditions [2]. The approach used for the Pilbara region was modified to an IFM as suggested by Farquharson et al., (1992) [54]. This recommended approach is an IFM with Q10 as an index variable and a dimensionless growth factor (GF) for X% AEP (GFx): A prediction equation was developed for Q10 as a function of catchment characteristics, and regional growth factors were developed based on the estimated at-site flood quartile. A Bayesian parameter estimation procedure with LP3 distribution was used to estimate flood quantiles for each gauged site for AEPs of 50%, 20%, 10%, 5%, 2% and 1%. Rahman et al., (2019) provide further information surrounding this process [52]. The adopted predictive equation for the index variable Q10 has the form: Where b0, b1 and b2 are regression coefficients, estimated using OLS regression, area is the catchment area (in km 2 ) and I6,50 is the design rainfall intensity at catchment centroid for a 6 h duration and 50% AEP. The values of b0, b1 and b2 and the regional growth factors (GFx) are embedded into the RFFE Model 2015.

Direct Rainfall Modelling
Direct Rainfall Modelling (also known as rain-on-grid) was undertaken in TUFLOW HPC, a 2D fixed-grid, adaptive time-step, hydrodynamic solver that uses an explicit finite volume solution [55]. TUFLOW HPC reduces the run time of models. The direct-rainfall approach applied the rainfall hyetograph (mm versus time) uniformly to active cells within the defined grid of the catchment of interest. Each hyetograph value represents the rainfall that fell per increment. The double precision version of TUFLOW Classic was used in initial model set-up to minimize initial model errors (such as deficient or erroneous data) before running greater numbers of simulations in the GPU for faster run-times after model establishment. TUFLOW Classic uses a fixed time step and will highlight any initial errors with the model runs. Instabilities in the model highlight bad data or poor model set up [56]. TUFLOW HPC can hide poor model set up through its adaptive time-stepping to ensure the model remains stable. Model simulation parameters (grid size, time-step) were established to optimise the accuracy, run-time, and stability of the model. To effectively resolve flow events within the channel it is recommended to provide at least 5 grid/mesh elements laterally across the river channel [57].
The time-step for the model runs was maintained at 1/5 of the selected model grid size in meters (a time-step appropriate for TUFLOW classic) [55] with a cell wet/dry depth of 0.0002m to account for the high proportion of shallow flow with a direct rainfall model. Due to the small catchment size, rainfall values are small and the reported IL/CL values for larger Pilbara river channels are scaled to these larger catchment hydrological inputs. Simply, the observed losses in larger catchments are larger than the grid-averaged rainfall inputs in small headwater channels. Therefore, IL/CL values were not applied during the final modelled scenarios. Hall (2015) provides a description of the advantages and disadvantages of direct rainfall modelling through model construction, calibration validation and sensitivity analysis [13].

Catchment DEMs
Headwater channels (first-to-third order) are the areas from which water originates within a channel network and are closely coupled to hillslope processes [58]. This study used a high-resolution dataset around the periphery of the Yandi iron ore mine in the Pilbara, WA (Figure 4). Catchments were selected based on the high-resolution dataset covering the full extent of the catchment. Surface catchment DEMs were derived from a wide area semi-global matching (SGM) survey [59] (horizontal accuracy = 0.5 m and vertical accuracy = 0.25 m) undertaken in the Upper Fortescue catchment. Additionally, the catchments were screened for minimal catchment disturbance, such as mining, agriculture of the hydraulic alteration of the waterway from bridges and culverts. Catchments with engineering infrastructure, such as railways, culverts and main roads were modelled until the upstream contact with these features. Ten catchments (ranging from 0.96 km 2 to 9.23 km 2 ) matched these requirements and were selected to test the RFFE approaches (Table 2).

Rainfall
The Bureau of Meteorology, Australia (BOM) 2016 Design Rainfall Data System [60] was used to provide rainfall intensities for RFFE analysis and as a modelling input for direct rainfall modelling. The BOM 2016 IFD values provide a rainfall intensity (mm/h) or depth (mm) based on the latitude, longitude, and catchment size within a set location ( Figure 5). BOM2016 IFD values replaced the older 1987 IFD and interim 2013 IFD values providing 30 additional years of hydrological data and adjustments to the approach. Additionally, the direct rainfall was applied to the 2D model over the entire catchment as time-series data as mm versus hours. The time-series data has a histogram stair-step shape) therefore rainfall was applied as a stepped approach holding the rainfall constant during the allocated time interval (BMT, 2018) meaning each rainfall value is the amount of rain that fell in mm between the previous time and current time. Rainfall was applied to every active cell within the digitised catchment. Ten varied areal hyetograph patterns were used for each catchment ( Figure 6). Podger et al., (2018) provide a detailed description of the creation of these hyetograph patterns [61].

Roughness
Manning's n values were assigned to account for runoff conditions. The roughness values assigned within a direct rainfall model can affect the timing of runoff. Additionally, constant roughness values may underestimate the effective roughness and ignores the role of spatially varied roughness within the catchment. Catchment floodplains were assigned a Manning's value of n = 0.02 and the channel was assigned a default Manning's value of 0.035. Adjustments of the roughness values were made in the sensitivity analysis of the TUFLOW models.

Discharge
For each catchment a cross section was delineated near the catchment outlet, (Figure 7). At this cross section, a plot output was created in the form of a time-series hydrograph. This time-series provides the flood flow through the catchment during each model scenario and is used to identify peak discharge (Qpeak) values for each rainfall event. Backwater development has been reported during flood modelling of channel confluences within the larger Marillana Catchment [62]. Therefore, cross sections were delineated slightly upstream from the catchment outlet to reduce backwater flow events impacting peak discharge values during flood modelling. In addition to using the plot output for each simulation, maps of water depth (d) and velocity (V) were constructed at each of the model time-steps. The output from this was used to make an independent estimate of discharge using a velocity-area method: where Q is the discharge expressed in cubic meters per second (m 3 s −1 or cumecs), A is the stream cross sectional area (m 2 ) and V is the mean velocity of stream flow (m s −1 ). Measurements of stream water depth (stage) are typically measured at sites within the Pilbara, however continuous flow measurement of river discharge is expensive and logistically unwise given the annual frequency of flow events within the region.

Sensitivity Analysis
Sensitivity analysis or calibration methods are critical steps in rainfall-runoff and developing useful models of complicated hydrologic systems [63,64]. Models can be calibrated to observed data to demonstrate that the model can produce an observed flow time series with an acceptable level of accuracy [64]. Alternatively, a model may be available that has been previously calibrated for a catchment as part of another study. Sensitive model parameters should be recognised and appropriately evaluated to ensure they are constrained within acceptable ranges. Prediction in ungauged basins is challenging to validate owing to the data limitations within the area. When models are not able to be calibrated to measurements sensitivity testing should at least be carried out to assess the sensitivity of the model to variations in the main model parameters [11]. Model evaluation may not be limited to how accurately model predictions match historical observations, but how well the model represents the hydrological system. In this paper, we use a direct rainfall model to quasi-validate and select the most appropriate RFFE procedure for headwater, small size catchments. All RFFE procedures used in this investigation have been previously calibrated on larger catchments in the Pilbara.
There is no standard method for estimating uncertainty in streamflow in ungauged basins using regionalisation techniques [44]. Uncertainty is estimated here in the sensitivity analysis of the direct rainfall modelling approach by adjusting model parameters on three catchments. (Figure 8). Sensitivity analysis was carried out to determine the optimum model running conditions and explore the suitability of direct rainfall modelling within these small headwater catchments. Catchments were selected on the basis of channel gradient to encompass a range of catchment types within the analysis. Catchment 5 was the steepest, with greater expected areas of supercritical flow and areas of complex terrain which would challenge model performance [56]. Other catchments included catchment 1; a small, shallow, unconfined channel with multiple flow paths and 10; a larger catchment with a predominantly unconfined single-thread channel representative of many of the headwater channels in the vicinity. Existing direct rainfall models have been found to be most sensitive to Manning's roughness and rainfall [13] and therefore sensitivity analysis was carried out to address the parameters that have larger uncertainties within the model. Sensitivity analysis was carried out to address: rainfall hyetograph shape, Manning's roughness (adjusted to +/− 20% or the upper and lower bounds for characteristic minor natural streams [65] and grid sizing, assessed with a 1 m, 2 m, 5 m and 10 m spacing.

Procedure for Evaluating RFFE Approaches
To evaluate the RFFE approaches, we use the result of the direct rainfall model to act as a measured series in this "quasi-validation". This quasi-validation uses joint plots of the RFFE output and simulated rain-on-grid discharge to compare the output of predictive methods. The ARR RFFE approach provides confidence limits of 5% Lower confidence limit (LCL) and a 95% Upper confidence limit (UCL) which are used as the absolute cut off for RFFE values. The evaluation procedure was as follows: (1) Was the RFFE output within the 5-95% confidence limit range for the ARR approach? (2) Did the approach demonstrate appropriate hydrological scaling across space and flood return period? (i.e., did the results increase with an increase in rainfall magnitude or catchment area) and lastly, (3) Did RFFE results display good agreement with findings from the rain-on-grid model? If not, did any approaches align with other tested RFFE methods.

Estimates of Peak Discharge for RFFE Methods
Rainfall values for RFFE approaches were obtained from BOM Design Rainfall Data System (2016) using 2016 Intensity, Frequency and Duration (IFD) values. The rainfall output showed a relatively homogenous depth for catchments ranging between 0.96 km 2 and 2 km 2 . Larger catchments (sized 5.99 and 3.23 km 2 ) were predicted to have higher rainfall depths for the 12 h storm across the range of AEP (Table 3). A wide range of Qpp values were obtained using the RFFE equations ( Figure 9). The Flavell RFFP2000 procedure produced higher Qpp values for Q2 events. Between Q5 and Q100, the ARR values were highest, with greatest agreement with the IFM and Flavell RFFP2000. The PRT and QRT approach from Taylor et al., (2011) show nearly identical values with no increase in peak discharge values for an increase in catchment size. The fixed region PRT (Rahman et al., 2012) provided the lowest estimates of predicted peak discharge. Both QRT and PRT approaches produced low discharge values across Q10-Q100 flow events.  9. Comparison of small flood predicted peak discharges using RFFE approaches for Q2, Q5, Q10, Q20, Q50 and Q100 return intervals. Table 4 shows the output values from direct rainfall modelling from the 10 catchments. Figure 10 show the results of direct rainfall modelling. These values have also been plotted next to the RFFE approaches with similar output (ARR RFFE (2019), IFM (Davies and Yip, 2014) and the Flavell RFFP (Flavell, (2012)) selected from deductive reasoning, omitting PRT and QRT approaches that yielded uncharacteristically small flood magnitudes across all annual exceedance probability.  Figure 10. IFM, ARR RFFE and RFFP2000 approaches compared with results from direct rainfall modelling from TUFLOW for the modelled headwater catchments.

TUFLOW Sensitivity Analysis
The results of the sensitivity analysis are found in Figure 11. The grid size of the model caused large variability in peak discharge values. Changes in peak velocity values were sensitive to this change in grid size, resulting in higher velocities (as one would expect) for a higher resolution DEM (e.g., 2 m) compared to the 10 m grid size. This increase in velocity is attributed to the more detailed topographic representation of channel constrictions and channel bed heterogeneity in the higher resolution grid size ( Figure 12). C5 was the steepest and most topographically varied catchment, with steep confining areas. Within this catchment there was greater variability in the peak discharge from adjustments to Manning's n value, hyetograph shape, (and most prominently) grid size; where Qpp values were doubled between a grid size increment increase of 5 m to 10 m. Final model scenarios used a grid size of 2 m, rainfall pattern 6, which resulted in the highest peak discharge in sensitivity analysis and a channel-wide Manning's roughness value of 0.035. The range of output values for the sensitivity analysis are found in Table 5.

Sensitivity Analysis
Sensitivity analysis of the direct rainfall model has indicated that the steeper, confined C5 is more sensitive to alterations in grid size, Manning's n and hyetograph shape. Further research and integration of IL/CL values would improve the model's representation of headwater channel processes. A grid size of 2 m was selected to avoid excessive computation time but to integrate detailed topographic forms such as smaller flowpaths and to capture the heterogeneity of the natural environment. Higher within a catchment, flowpaths become smaller and may be poorly represented by the model if they exist on a sub-grid scale, affecting the timing of runoff routing within the catchment [66] (ARR, 2012). This is likely a key contributor to the increased Qpp values for coarser resolution model grid sizes (5-10 m). Mesh (grid-size) resolution has been shown to have a high impact on a models output flow volume in other studies conducting direct rainfall in small catchments [23]. Additionally, the slope between neighbouring cells has a key influence on the behaviour of catchment runoff, and therefore peak discharge values. Care should be taken to adequately understand the potential variability in peak discharge estimation using direct rainfall models. Additionally, explicit consideration of local variation in rainfall patterns in peak flood discharge modelling efforts is advised.

Applicability of Direct Rainfall Modelling in Headwater Catchments
Direct rainfall models are subject to higher levels of mass error when using the direct rainfall approach particularly where the model has areas of steep, complex flow or the model is located at a high elevation above sea level and experience relatively small inflows [56]. The catchments modelled here have many of these challenges. They are small catchments with small inflows, and experience sudden changes in slope conditions. Drainage pathways are not clearly defined and for some catchments, the hydrological boundaries between channel and floodplain are difficult to model. Flow is routed along the floodplain or within a claypan environment culminating in diffusive and complex drainage patterns. For these catchments, there were unexpected peak discharge values which may not fully represent what would occur in these catchments outside of where flow is rapidly dispersed away from channel setting (e.g., catchment 10).
The catchment outlet of several small headwater channels was prone to backwater effects. Such backwater effects have been reported for other channels in the Pilbara region from constrictions within the channel [34]. Backwaters can decrease velocities, raise the water surface elevation (WSE) and extend upstream [34]. Peak discharge values obtained from the catchment outlet were based on a cross section beyond the extent of backwater effects from outlet constrictions or tributary junctions. Care was taken to inspect the mapped output to interpret the water surface elevation at catchment cross sections.
Despite these warnings, the use of direct rainfall modelling is useful within small catchments with well-defined drainage pathways to quasi-calibrate selected RFFE approaches in the complete absence of rainfall and streamflow data. In doing so, an approximation of likely flood flows within these ungauged settings can be made. However, direct rainfall modelling with these catchment characteristics should not replace a hydrologic modellers perspicacity in dealing with flow estimation within dryland headwater environments. Even when calibration is properly done, models tend to have greater predictive strengths over shorter timeframes than longer timeframes as the system over short time scales is more similar to the one which it was calibrated for. Therefore, greater trust can be placed on the more frequent flow predictions than longterm flow dynamics within these headwater catchments.

RFFE Evaluation
The evaluation procedure was as follows: (1) Was the RFFE output within the 5-95% confidence limit range for the ARR approach. (2) Did the approach demonstrate appropriate hydrological scaling across space and flood return period (i.e., did the results increase with an increase in rainfall magnitude or catchment area) and lastly, (3) Did RFFE results display good agreement with findings from the rain-on-grid model and if not, did any approaches align with other tested RFFE methods. Many of the RFFE approaches for the Pilbara region are created using far larger catchment sizes. For example, the Flavell RFFP2000 approach is derived from an average catchment size of 5570 km 2 and this research tests the application of this RFFE in catchments with an average size of 2.16 km 2 . Despite this magnitude variation in catchment size, these RFFE approaches are used in these smaller catchments for design discharge calculations. However, the direct modelling highlights that catchment representativeness remains an issue. This issue is exacerbated by sparse gauging coverage and continues to be a high priority for future research [66].
The RFFP2000 (Flavell, 2012), IFM (Davies and Yip, 2014) and ARR RFFE (2019) methods produce higher estimates of peak discharge. These three methods diverge from the other approaches in estimating higher peak discharge estimates and this deviation is more pronounced in less frequent flood flows (Q20, Q50 and Q100). These three methods are also more sensitive to changes in catchment area with larger catchments (above 5 km 2 ) predicted to have peak discharges double that of smaller catchments (around 2 km 2 ).
Both PRT methods estimate the lowest peak discharge values across all ARI, with little increase in discharge across less frequent flood flow events. These estimation techniques are less sensitive to catchment size and may overlook flood flow magnitudes within smaller catchments as they were created with a dataset comprised of larger catchments. Additionally, parameter regression techniques give emphasis to the mean and therefore frequent flows [67]. These methods provide a lower estimate of peak discharges, with emphasis on frequent flow events and are likely not suited to smaller headwater catchments such as those in this study The QRT methods show slightly higher predicted peak discharges but still do not provide convincing flood discharges at higher return intervals. PRT and QRT methods should be used with caution on smaller headwater channels. The IFM (Davies and Yip, 2014) method uses greater weighted frequency factors which results in the larger range of peak flood values across return intervals. This greater range of frequency factors for the IFM method compared to the PRT methods is likely to better represent higher magnitude flow events within these headwater catchments.
The ARR RFFE (2019) [25] approach is industry standard for larger catchments and is useful for small headwater catchments even if it is likely to provide a conservative estimate. The ARR RFFE model is noted to have large uncertainty with mean relative errors of 50-60% [66]. Alternative methods developed by Flavell (RFFP2000 and RFFP2006) and by Davies and Yip (2014) have been suggested as viable locally developed replacements to ARR2015 methods for the Pilbara region [66]. The Flavell RFFE approach is relatively complicated in comparison to the IFM and ARR (RFFE) 2019 approaches. The Flavell RFFP approach requires some catchment analysis to provide equation values (such as slope, latitude, longitude, length of catchment) but it is likely these values would be easy to populate for any catchment of interest. The Davies and Yip (2014) IFM approach uses a wide spatial distribution and catchment range in developing design equations making it applicable throughout the Indian Ocean drainage division [18]. The RFFP2000, IFM and ARR RFFE (2019) are therefore suggested for the estimate of peak discharge values within small ungauged headwater channels.

Conclusions
The role of this research was to identify peak discharge values from a direct rainfall model, and to compare these findings with RFFE approaches conducted within the region.
Within the TUFLOW direct rainfall model simulations, the model was predominantly sensitive to grid size and Manning's n, with the rainfall pattern having a smaller impact on Qpeak. It is worth noting however that even the variation in the shape of the hyetograph produced up to 50% variation in the final Qpeak value within the sensitivity analysis. Therefore, awareness of the limitations of direct rainfall modelling within these steeper catchments is necessary when using this approach or selecting an RFFE approach. Despite this caution, the Flavell RFFP2000 and Davies and Yip, (2014) IFM appear the most reasonable estimates of peak discharge within the channels. These approaches provide lower estimates than the ARR2019 RFFE models, which are suggested to provide higher flow rates and are commonly used as part of a conservative approach to waterway design within the Pilbara. Whilst the direct rainfall modelling approach has many barriers when applied to headwater channels, this approach is useful to provide a frame of reference when selecting an RFFE approach for ungauged headwater catchments. It is common for a RFFE approach to be applied without first considering the weaknesses and limitations of the derived equation, or its suitability within small headwater catchments. Here, we show the potential for direct rainfall modelling and its applicability in an environment with limited opportunities to validate peak flood discharge estimation methods. The lack of gauged real-world measurements within these catchments continues to present the main obstacle in improving our understanding of flow conditions and rainfall values within the Pilbara region.
Supplementary Materials: The following supporting information can be downloaded at: www.mdpi.com/article/10.3390/hydrology9100165/s1, Table S1: Supplementary material highlighting key method details for the RFFE estimation methods applied to headwater catchments in the Pilbara, WA.
Author Contributions: Conceptualization, A.F. and I.R.; methodology, A.F. and I.R.; software, A.F., writing-original draft preparation, A.F.; writing-review and editing, A.F. and I.R. All authors have read and agreed to the published version of the manuscript.