Modelling the Potential of Integrated Vegetation Bands ( IVB ) to Retain Stormwater Runoff on Steep Hillslopes of Southeast Queensland , Australia

Rainfall intensity is predicted to increase under a changing climate, leading to increased risks of hillslope erosion, downstream sedimentation and flooding. For many catchments used for grazing and agricultural land uses, it will become increasingly important to maintain ecohydrological functioning despite climatic extremes. One means to achieve this is through strategic reforestation using locally endemic species, in spatial configurations that effectively intercept, retain or and redistribute overland flows. This paper adopts a modelling approach for investigating the potential of one such design termed “integrated vegetation bands” (IVB), to increase the retention of runoff across steep hillslopes, particularly in the sub-tropics where rainstorms are becoming increasingly intense. A spatially distributed simulation model (MIKE-SHE) was applied to a steep, grazed catchment (Maronghi Creek catchment, Southeast Queensland, Australia) to compare stormwater runoff characteristics between: (1) the existing pasture land cover; and (2) a series of hypothetical IVB added across this pasture land. The IVB were approximately 20 m wide, and configured at 5% gradient towards ridgelines. Results for estimates of overland flow depth and infiltration (spatial), and accumulative water balance (temporal), confirm that OPEN ACCESS


Introduction
The reduction of woody and herbaceous vegetation at the hillslope scale is known to increase the volume, velocity and erosive potential of stormwater runoff [1].Many thunderstorms in the sub-tropics can generate rainfall rates well over 100 mm/h [2], and the frequency of severe thunderstorms with intense rainfall is projected to increase under climate change [3].In the study area (mountainous terrain of Southeast Queensland), rainfall intensity is projected to increase by 40% by 2030 [4].This poses a high risk of soil erosion and slope failure for this agriculturally productive region, and is likely to increase the risks of downstream sedimentation and flooding.
Native trees and shrubs modify soil structure and ground surface properties by the production of a surface litter layer, increased soil organic matter and deep root and biopore development.In native woodlands, rain infiltration rate, soil water storage capacity and soil biological activity can vary with changes in accumulated litter and soil organic matter [5,6].While these processes occur naturally in many native ecosystems (e.g., [7,8]), they are often impaired in landscapes that support grazing and cropping land uses.
The structural form of vegetation is also important for determining the ecohydrological function it is able to perform.For example, Duran Zuazo et al. [9] found that the lowest-growing plant covers reduced runoff and soil erosion more effectively than taller sized shrubs on hillslopes of the Sierra Nevada Mountains (SE Spain).Leguedois et al. [10] found that vegetation bands (tree belts) can effectively trap sediment travelling down slope in overland flow by sedimentation processes.Blanco-Canqui et al. [11] also demonstrated that vegetative filter strips applied in cropping landscapes could reduce water velocity and increase sediment and nutrient retention.Filter strips alone, however, are unlikely to generate micro-terraces to the same degree as vegetation bands, as they are missing the woody debris component.Plant litter containing both herbaceous and woody components will form meta-stable micro-scale terraces that directly promote the emergence of backwaters.
Through its effects on soil and surface properties, perennial vegetation has the capacity to reduce surface flow velocity, enhance water infiltration and sediment retention and reduce hillslope runoff [12][13][14].This suggests that by increasing the area of perennial vegetation across hillslopes, the negative impacts associated with extreme rainfall events may be reduced.To achieve these beneficial ecohydrological functions while at the same time minimising competitive effects, perennial vegetation can be used in banded forms that integrate with existing agricultural production (i.e., integrated vegetation bands-IVB).
Vegetation bands can buffer landscapes from intense rainfall by forming semi-permeable through-flow systems that impede and reduce the velocity of overland flows [12,13].Leaves, bark and sticks (plant litter) along with multi-stemmed herbaceous vegetation covering the soil surface, are all efficient at intercepting water, sediment and nutrients [14,15].Litter forms micro-terraces that, in turn, generate "backwaters" [16,17].Backwaters have a particularly intriguing property-they can cause water to spread laterally across slopes [18].This raises an interesting prospect: can IVB be configured to help spread runoff by gravity, so that water velocity is reduced while retention is increased over a greater area of the landscape?
The aim of this desktop study was to simulate how IVB affected resistance to flows down steep hillslopes, and in particular, whether water retention could be increased by passively redistributing the flows across the hillslopes.The MIKE SHE [19] hydrological model was used to simulate the effects of adding IVB to steep grazed hillslopes under intense (up to 63 mm/h) and prolonged rainfall conditions (five days).The IVB simulation was compared directly with the existing pasture land cover in the Maronghi Creek catchment of Southeast Queensland, Australia.The simulations were used to calculate changes in water balance and the spatial distribution of infiltration and overland flow at hillslope and sub-catchment scales.Sensitivity analyses were applied to determine each variables effect on model output.The information gained from these analyses were subsequently used to optimise the design of instrumentation in a paired catchment study that is currently underway in Southeast Queensland (ARC Linkage Project LP100100356).

Conceptual Model
Research on vegetation that influenced the IVB concept include banded vegetation [20], litter micro-terraces [7,8], vegetation filter strips [11,21], and contour bund layouts as used by some innovative farmers [22,23].The particular design applied in this study is adapted from earlier work by J.G. Ryan [24].The IVB simulated in this study were 20 m wide, set at approximately 120 m apart, and configured to follow topography at 5% gradient to the contour angling toward the ridgelines (Figure 1).These parameter choices for IVB are based on a compromise between the land footprint being small to minimise the competition between pasture and tress, whilst still allowing the ecohydrological benefits of the tree bands to influence catchment water balance.Other widths and gradients may offer similar functional benefits to the IVB used in this study.

Study Area
The study area forms part of the headwaters of Maronghi Creek catchment (152°10ʹE, −27°11ʹS), a tributary of the Brisbane River in southeast Queensland, Australia (Figure 2).The average annual rainfall is approximately 709 mm based on Cressbrook Dam (6.8 km to the south), with about 70% falling in the summer months (October to March) as intense thunderstorms or rain depressions.The study area covers 95 ha, varies from 330 to 440 metres in elevation, and has generally steep topography with slopes up to 24%.This catchment was chosen as it was devoid of native perennial vegetation cover, and contained steep grazed land on unstable soils.
The sub-catchment geology is Mesozoic to Proterozoic age igneous rocks termed the Eskdale Granitoids (Triassic), which form an extensive drainage network and steep hills containing little soil development [25,26].The soils are classed as Orthic Tenosols in the Australian Soil Classification system [27], which are broadly similar to Cambisols using the world reference base (WRB) classification [28].
The native vegetation is grassy eucalypt woodland comprising Eucalyptus melanophloia (silver leaved ironbark), E. crebra (narrow leaved ironbark) and Corymbia erythrophloia (gum-topped bloodwood).The landscape has been largely cleared to promote grazing of native Heteropogon contortus (black speargrass) pastures.The parameters used to define the physiology of the trees within the IVB were based on these native woodlands.

Software
The water balance simulations were modelled with MIKE SHE [19], a fully distributed hydrological water balance model that accounts for spatial heterogeneity using physically based equations [29].MIKE SHE is a part of the MIKE ZERO family of hydrological systems modelling software, as developed by the Danish Hydrological Institute (DHI).DHI's Gold Coast office (Nerang, QLD, Australia) provided the software as an in-kind contribution to the project, and the model development was assisted by their engineers.Simulations used standard MIKE SHE process sub-models for interception, evapotranspiration, overland flow, and unsaturated and saturated flow.These were finite difference for overland flow, simple 2-layer for unsaturated flow, and finite difference for saturated flow.To account for spatial variation in rooting depth and evapotranspiration, differences between pasture and IVB land cover types are determined in reference tables using gravity flow and Richards ET.Richards ET approximates the movement of water in the unsaturated zone of soils due to capillary and gravitational forces.
In the study area, Eskdale Granitoids weather to kaolinite and illite clays and albite [25].Soils are slightly acidic (pH 6.0-6.5),shallow to moderately deep, very dark grey to brown, structure-less "Earthy Sands" [39], that occur between granodiorite rock outcrops, and are often underlain by impermeable sub-soils (shallow clay hardpan and rocks) [40].These clays are relatively stable to erosion but have poor retention of water and cations [41].
Soil parameter were derived from the "ASRIS" database [36,38].Based on the relationship between maximum volumetric water content at saturation equalling total soil pore volume, bulk density values for the surface soil (1525 g/m 3 ) and sub-soil (1733 g/m 3 ) were converted to soil saturated water content (SAT) values of 0.43 and 0.35, respectively.Sivapalan [42] reported SAT at 0.42 for an Orthic Tenosol in a cropping landscape in south-eastern Australia (a more temperate climate than the study area), but given the influences of steep terrain, soil compaction from loss of biopores and organic matter, SAT for IVB was set at 0.38 and the grazed open pastures at 0.24.
Saturated hydraulic conductivity (Ks) of the surface soil was 14 mm/h (pasture) and 19 mm/h (IVB), while the saturated zone horizontal and vertical hydraulic conductivities set to 9.75 × 10 −6 m/s and 1.6 × 10 −8 m/s respectively.These values are consistent with the association of Tenosols with high drainage and low runoff [43].Soil depth was set at 300 mm (pasture) and 650 mm (IVB) to account for differences in compaction and erosion of soil under the pressure of grazing compared to deeper profiles that form in native woodlands under the influence surface plant litter and deep tree roots/biopores.
Surface roughness variability across the hillslopes was represented by a spatial grid containing Manning's M values, also known as "Strickler's value" (Manning's M is the inverse of the more conventional Manning's n).For example, values of n of 0.01 (smooth channel) and 0.10 (thickly vegetated channel) correspond to M values of 100 and 10 [19].Overland flow across hillslopes often reflect much coarser surface roughness than for channels.To reflect this, the M value for pastures was 12.5 [34], and 2.5 for IVB based on the added surface roughness of fallen timber and plant litter under trees [37].These are similar values to those published for other woodlands (e.g., [44]).
Estimates for detention storage (DS) can vary considerably [45], from 5 mm for bare sandy soils and up to 32 mm for some native forests [46].For a sub-humid, semi-arid Eucalypt woodland in northern Queensland, Bonell and Williams [47] estimated DS at 3.6 m (±1.6 mm), whereas Yirdaw and Luukkanen [48] estimated plant litter depths at 5 to 6 mm for young plantation eucalypts.To reflect the contrasting effects of grazing and litter build up under mature eucalypt woodlands in a sub-tropical environment, we used DS values of 3 mm for pasture and 8 mm for IVB.
Physiological differences between grasses and trees, including leaf area index (LAI) [32], rooting depth and evapotranspiration, were included in the model despite the short simulation time (one week).Evapotranspiration was 1.60 mm/day for pasture and 2.60 mm/day for IVB, which equates to approximately 949 mm/yr, similar to the 2.61 mm/day (953 mm/yr) estimate stated by Gordon et al. [49] for average annual evapotranspiration in open eucalypt forests in Australia.Potential evapotranspiration (PET) [33] was based on the long-term average from Cressbrook Dam of 1505 mm/yr or 4.12 mm/day.
This study focused on modelling a representative extreme rainfall event following a dry period.The modelled rainfall event was derived from a six minute pluviograph record collected at the Cressbrook Dam weather station [50] between the 6 and 12 February 1999 (Figure 3).Although more recent heavier rainfall events have occurred in this area (e.g., January 2010), this event was selected as being "representative" of a dry period where soils may have low ground cover that are subsequently exposed to heavy rainfall and gradually become more saturated over several days.Most of the rain fell between the 7 and 11 February (276 mm in total), with a maximum six minute rainfall intensity of 61 mm/h.The average recurrence interval (ARI) for this event was 1:20 years, with a six-minute time-step used for the model due to the extreme variability in rainfall rates that can occur over short timeframes [51].
A 10 m DEM was derived from 5 m contour data and overlain on a rectified 1:5000 orthophoto of the study area to check for spatial consistency.False sinks were identified manually, and eliminated using a post-processing fill operation to provide a hydrologically realistic surface.All spatial data were based on a grid-based or raster GIS format created in ArcGIS 9.1 at 10 m resolution.The existing land cover was mapped from the orthophoto, with IVB added through on-screen digitizing in ArcGIS.
Parameters (as defined above) were set as benchmarks from which sensitivity analyses were performed.Multiple simulations runs were conducted, in which each benchmark estimate was incrementally varied to higher or lower values.Simulated outputs of evapotranspiration (ET), infiltration (INF) and boundary outflow (Q) were compared for each set of parameter estimates and entered into MIKE SHE based on their mean and standard deviation.

Hillslope Scale
Two hillslope profiles were taken from either side of a central ridge within the catchment to compare the response of overland flow (OF), and its distribution (depth) for pastures (P) and IVB (Figure 4a,b).These profiles correspond to transects A-B and A-C in Figure 2. The two profiles begin at a common point (A), and are orientated down the steepest gradient to points B and C respectively.Two time-steps are shown in each profile, one time-step (TS362) during a very heavy rain period at the beginning of the simulation, the other time-step (TS639) when the catchments soils were fully saturated.At TS362 for the pasture land cover, OF depth was 6 mm and 6 mm for transect A-B and 6 mm and 7 mm for A-C.Under IVB at TS362, OF depth was 16 mm and 14 mm in transect A-B and 14 mm and 14 mm in A-C.This represents up to a 266% increase in OF depth for transect A-B and 200% for A-C under IVB compared to the pasture land cover at the corresponding time and location.Under saturated conditions (TS639), the pasture OF depth was 6 mm and 7 mm (A-B) and 7 mm and Close observation of Figure 4b (transect A-C) also highlights an interesting shift in the location of the peak in OF depth where it interacts with the IVB.During the period of heavy rainfall (TS362), OF is at its deepest at the top edge of the IVB which indicates a backwater has formed, whereas during the saturated conditions (TS639), OF is deeper within the IVB themselves, suggesting the extra surface roughness of the IVB (i.e., plant litter) helps to hold more OF at such times.In addition, during the heavy rainfall period where the backwater formed before the IVB, OF depth concurrently decreased immediately down-slope (5 mm) of the IVB compared to the corresponding location in the pasture land cover (7 mm).Flow directions and velocity results show that during the maximum rain rate on the same section of the main ridgeline (Figure 5), water depth was reduced downslope of the IVB but slightly increased in the backwaters before the bands.The direction of the vectors also moved towards the ridge at these locations, increasing slightly in velocity before the IVB and then decreasing within the bands themselves.This supports the conceptual model of IVB functioning-runoff creates backwaters behind vegetation bands, which in turn, moves runoff laterally across the hillslope.These modifications to runoff resulted in OF depth on the main ridge (Figure 2, inset) being +28% (0.009 m, S.E.1.46 × 10 −4 ) in IVB simulations compared to pasture (0.007 m, S.E.5.01 × 10 −5 ), which was significantly different (p < 0.0001, n = 896).Changes in mean infiltration rates and distribution were also evident.Based on a time-step when rainfall was at its most intense and the soils reasonably dry, the difference between IVB and pasture land covers was clear (Figure 6).During this period, infiltration was 7% and 10% greater at the sub-catchment and hillslope scales respectively.Infiltration slightly increased above the IVB in the backwater area, significantly increased directly beneath them, and slightly decreased immediately down slope of them.Considering the spatial extent of grid cells where both overland flow depth and infiltration were modified, it was apparent that the area of hillslopes retaining > 10 mm/day more runoff increased by 22% under IVB compared to the pasture land use.Given that the IVB covered an area equal to 11% of the catchment area, this reflects an increase of 11% in the area where runoff was retained outside of the IVB themselves.

Sub-Catchment Scale
At the scale of the sub-catchment, a grid based time-series of maximum depth in overland flow (OF) illustrates a greater water depth under IVB than for pasture, but little difference in water depth in the channels (Figure 7).The time-series were based on 24 h time-steps beginning 12:06 pm 7 February and ending 12:00 pm 11 February 1999.One exception to these results: the extent of saturation on footslope areas was less in the IVB simulations.
Changes to overland flow depth recorded during the rainfall event showed significant variability through time, as well as between hillslope and sub-catchment scales (Table 2).Mean and maximum overland flow depths were greater in the IVB at the hillslope scale, although at the sub-catchment scale this trend was mostly absent.This implies that the greater water depth under the IVB was temporary during heavy rainfall, and after infiltration across the hillslopes, discharge at the sub-catchment scale was lower.However, towards the end of the simulation after rainfall had ceased, the IVB contributed to a slight increase in water depth in the channel compared to the pasture land cover.This had the effect of slightly extending the temporal period of baseflows within the stream.

Model Water Balance Predictions
The greatest changes in water balance between the pasture and IVB simulations at the sub-catchment scale were 23% greater maximum infiltration (INF) (54 to 66 mm) and 7% lower maximum boundary outflow (Q) (211 to 196 mm) (Table 3).Mean values were also significant for these two estimates, whereas the remaining water balance relationships showed only minor statistical differences.The main reason behind the small differences overall is that IVB cover only a small footprint of land, with pasture land cover dominant in both simulations.

Soil Water Content at Saturation (SAT)
Changes in estimates of soil water content at saturation (SAT) had a significant bearing on simulated water balance outputs (Figure 8a).As SAT increased INF also increased, whereas Q decreased.Variability in water balance outputs for the pasture land cover was much greater (ET σ 1.34; INF σ 29.33; Q σ 34.87) than for comparable IVB simulation runs (ET σ 0.05; INF σ 5.71; Q σ 6.79).This implies that SAT has a major bearing on the retention of runoff on hillslopes when it decreases through land use management (i.e., grazing on steep slopes).

Soil Depth (ACQ)
In the MIKE SHE model, soil depth was limited by the depth to an impermeable layer or aquifer (ACQ).Varying the depth to this layer had a significant effect on model outputs when less than 0.3 m, otherwise it had little effect (Figure 8b).This is evident in the standard deviation across all depths (ET σ 1.68; INF σ 18.26; Q σ 19.86) compared to the standard deviation when the 0.3 m depth was excluded (ET σ 0.06; INF σ 1.35; Q σ 0.04).The most noticeable effect when depth was at 0.3 m was a higher estimate of Q (40 mm), whereas INF was reduced by a similar margin (−35 mm).There was little change in ET for any depth.

Detention Storage (DS)
Detention storage had a moderate effect on ET, INF, and Q (Figure 8c).Increasing the estimate of detention storage slightly increased ET and slightly decreased Q, with little effect on INF.The total variation in water balance estimates across both pasture and IVB model runs combined was low (ET σ 0.76; INF σ 0.16; Q σ 1.02).Estimates varied more in pasture simulation runs (σ ET 1.16; INF σ 0.25; Q σ 1.54) than in IVB simulation runs (σ ET 0.13; INF σ 0.03; Q σ 0.45).This implies that detention storage is a short-term process, important during and for a short time after rainfall occurs, after which the detained store of water evaporates or infiltrates.

Surface Roughness (M)
Modifying estimates of surface roughness (M) had a slight effect on Q (Figure 8d).As estimates of M were increased (i.e., less rough) water balance outputs indicated that Q increased, with little change in ET or INF.This was evident in the standard deviation (σ) of ET across all estimates (ET σ 0.04; INF σ 0.01; Q σ 0.79).Variability in water balance estimates were again greater for pasture simulation runs (σ ET 0.05; INF σ 0.003; Q σ 1.05) than for IVB simulation runs (σ ET 0.01; INF σ 0.02; Q σ 0.44).

Hydraulic Conductivity of Surface Soil (Ks)
Varying estimates for infiltration rate in the surface soil layer (Ks) had only a slight effect on Q, and a small effect on INF or ET (Figure 8e).The reason behind the small effect on infiltration is that soils were saturated and unable to "infiltrate" more water.Including estimates from both pasture and IVB simulation runs combined, the variation in outputs were minimal (ET σ 0.13; INF σ 0.14; Q σ 0.55).Model outputs for the pasture land cover (σ ET 0.18; INF σ 0.21; Q σ 0.83) varied more than for vegetation bands simulation runs (σ ET 0.04; INF σ 0.05; Q σ 0.07).

Major Findings and Comparisons to Other Vegetation Systems
The simulations and sensitivity analyses indicate that changes to soil water holding capacity and soil depth are the major drivers of increased water retention in IVB compared to pasturelands without bands.Other surface related factors such as detention storage, surface roughness and rain infiltration rate have a slight effect on water retention individually, whereas their combined effect does have a noticeable effect on overland flow distribution and longevity of discharge in the stream.
Water flow was detained at greater depths in the IVB, causing water to spread and infiltrate across the slope.The actual area of hillslope where more than 10 mm of runoff was retained was increased by 22% in IVB compared to the pasture land cover, or an 11% increase of the hillslope area outside of the land footprint of IVB themselves.The extra water stored in IVB also contributed to a slight lengthening of base flow in the main drainage line.This suggests that where IVB are applied over a much larger area (many properties), some extension of base flow duration in ephemeral streams may occur.The results also show that peak discharge is lower under IVB, which may also reduce the potential for damage to the downstream fluvial network following heavy rainstorms.
The results are consistent with known effects of vegetation on catchment hydrology.For example, catchments with native cover of woody perennials often have greater water retention than for comparable, but deforested, catchments [52].Planted bands of perennial vegetation (tree belts) were also found by Ellis et al. [16,53] to retain significant runoff under moderate rainfall rates in a temperate agricultural landscape of southern Australia.Using a field based rainfall simulator situated on a Red Chromosol (duplex) soil of 12% slope, they found that tree belts could intercept 100% of runoff at rainfall intensities of 45 mm/h over 13 minutes duration, 50% at 45 mm/h over 30 minutes duration, but became breached by runoff at 75 mm/h over 30 minutes duration.
The spatial configuration of IVB applied in this study had a direct effect on runoff retention and distribution.Other research has also shown this for grazed semi-arid landscapes.For example, Bartley et al. [54] found hillslopes of the semi-arid Burdekin River catchment in northern Queensland, Australia, had up to nine times more runoff and 60 times more sediment loss under a heterogeneous mix of vegetation cover and bare ground compared to hillslopes with the same mean vegetation cover but no large bare soil patches.Similarly, Boer and Puigdefabregas [6] found that water and sediment yields from hillslopes in semi-arid SE Spain, were greater when vegetation was heterogeneously distributed compared to hillslopes with more uniform distributions.In semi-arid regions, patches, bands or isolated clumps form naturally depending on runoff and run-on characteristics [55].In semi-arid regions of Australia, banded vegetation patterns have been suggested to capture up to 8% more runoff than other vegetation patterns [55].
The orientation of the IVB modelled in this study helped to retain runoff and spread water more than the homogenous pasture land cover.Infiltration was directly affected, with slightly more infiltration directly beneath and upslope of the IVB and slightly less down slope of the bands themselves.The upslope increase is indicative of the effects of backwaters forming due to resistance to runoff created by the IVB.

Model Performance
The sensitivity analyses indicated that water balance outputs were most sensitive to incremental changes in parameters related to soil water holding capacity at saturation and soil depth.Other parameters that had an influence on model outputs, but to a lesser degree, were detention storage, surface roughness and infiltration rate.These variables had more effect on runoff distribution.The responses of the model to changes in parameters provide a valuable indicator of the causal effects instigated by IVB, which help to optimise the monitoring program for the full field-trial.The remaining parameters had little effect on model outputs and are not reported on here.

Saturated Soil Water Content (SAT)
Modifying estimates of soil water content at saturation (SAT) had a noticeable influence on water balance outputs.In general, higher estimates of SAT resulted in greater infiltration (INF) and decreased discharge (Q).In a study on steep semi-arid hillslopes of the Tigra region, Ethiopia, Descheemaeker et al. [56] found that water holding capacity, saturated hydraulic conductivity and bulk density were parameters that affected the volume of runoff.More specific to SAT estimates within MIKE SHE, Christiaens and Feyen [57] also found that their model was sensitive soil water content.Similarly, Oliver and Smettem [58] found that SAT and the soil water retention curve for a "sandy soil" had a significant effect on their modelling output.

Depth to Impermeable Layer (ACQ)
Soil depth (depth to hardpan) had a noticeable effect on water balance outputs.Scherrer et al. [59] consider that estimates of the depth to the impervious layer can greatly affect the potential storage volume of the soil, and in turn, determine how much water is needed to saturate the soil.However, Sahoo et al. [60] found drainage depth had little effect on MIKE SHE outputs in their modelling of a steep, high rainfall Hawaiian catchment.Over longer simulation periods the depth to an impermeable layer or aquifer may become more important because impermeable layers can limit rooting depth [61].Limits to rooting depth, in turn, can modify both surface and sub-surface lateral flow characteristics [62].This suggests that longer simulation periods would be needed to quantify how tree root and biopore development modify soil depth, depth to an impermeable layer or aquifer, and possible effects on hydraulic conductivity of the saturated zone.

Detention Storage
Variations to detention storage estimates had a less noticeable effect on model outputs.The increased detention storage under IVB slightly increased ET whereas Q slightly decreased, with little effect on INF.The surface detention storage estimates defined in MIKE SHE are representative of the development of a litter layer beneath IVB.This increase of detention storage reflects additional moisture trapped behind layers or micro-terraces of litter, which may then be available for infiltration.Guevara-Escobar et al. [63] found that the maximum storage capacity of litter from Populus nigra increased as rainfall intensity increased.These authors also cited work from Huber and Oyarzun [64] who found litter intercepted up to 8.9% of annual rainfall in a native evergreen forest in Chile, whereas the study of Balazs [65] found the litter of Larix decidua intercepted up to 18% of rainfall.

Surface Roughness (Manning's M)
The result of increasing surface roughness was to decrease discharge (Q).Jinkang Du et al. [66] also found that runoff on sub-tropical hillslopes in the Jiaokou catchment (259 km 2 ), China, was sensitive to changes in surface roughness (M).Similarly, Sahoo et al. [60] found that flood peaks decreased by nearly 1 m 3 /s when M was reduced from 60 to 10 (i.e., smooth to rough) in a mountainous catchment of Hawaii, USA.All of these studies (including this one) were based in steep terrain.It is likely that as slope gradient decreases, surface roughness coefficients would have a greater effect on runoff and water balance outputs.

Saturated Hydraulic Conductivity (Ks)
As infiltration increased under the pasture land cover, a slight increase in discharge (Q) was apparent, but there was little change in Q when infiltration rates were increased in the IVB.Although seemingly counterintuitive, as increased infiltration is generally expected to result in a decrease of Q, it is likely that more water was stored within the IVB themselves.Carroll et al. [67], for example, found that infiltration rates were up to 60 times higher under young planted trees than in adjacent pastures in a Welsh grazing landscape.The effect of increased infiltration and enhanced detention storage within the shelterbelts planted across this landscape resulted in a 23% reduction in runoff [68].Infiltration rates can also vary with topography, soil depth and soil hydraulic conductivity [35,69].Without knowledge of these parameters from field measurement or a spatial coverage at suitable resolution on which to base simulation models, it is not possible to estimate the degree to which infiltration modifies the partitioning of water to ET and Q.

Benefits of IVB
The major benefit of applying IVB in a landscape is that they can help to retain more stormwater runoff over greater areas of a hillslope.The changes in soil properties and increased resistance to flow from surface roughness are vital.A significant effect of IVB is the restoration of SAT to levels similar to native woodlands.Plant litter is an important pool of organic matter and nutrients available to be recycled to the soil, which also helps reverse soil compaction [70,71].The loss of organic matter enhances the effects of compaction on soil that, in turn, reduces the size and number of biopores and the spaces available for water to be retained in a soil [43].A loss of organic matter and biopores decreases SAT and increases dispersibility when saturated [72], whereas amendments of organic matter can help restore soil water holding capacity [73].If the species that comprise IVB are biologically diverse, than they also produce a great diversity of litter types which functions to increase the food available to a broad array of soil biota, which also form biopores.
A particularly important function of IVB is the increased surface roughness and detention storage emerging from the accumulation of plant litter.Plant litter and microtopography create rough surfaces that increase surface detention storage, rates of infiltration, hydraulic conductivity and soil water contents [74].Conversely, if erosion or fire removes the litter layer, there are fewer obstructions to overland flow which reduces soil water-holding capacity, while runoff and sediment removal is increased [75,76].Together, surface roughness and plant litter attributes increase stormwater retention by slowing the speed of runoff, and increasing detention storage and infiltration.Other ecosystem services that may accrue from IVB are enhanced soil ecosystem processes, water quality, carbon retention and biodiversity conservation [77,78].

Limitations of IVB
The major limitations to installing IVB is the complexity of designing them to match terrain and soil, their cost, competition with pasture/crops, stock trampling, breach failure, drought, fire risk and weeds/feral animals.Any potential functional gain needs to be weighed against these limitations.Trees can compete with crops and pastures [79], although this can be offset by root pruning which has been shown to create higher soil moisture compared to unpruned vegetation [80].Both fire [81] and heavy grazing can destroy the litter layer and micro-terraces within IVB, although grazing can be controlled with fencing.IVB may pose a potential source area for weeds and feral animals, although due to their compact nature, these should be more easily managed than for more homogenous layouts.Where other functions such as erosion control, shade and shelter, timber products, carbon storage or biodiversity habitat are integrated into IVB, their costs/benefits ratio are likely to increase substantially.

Research Limitations and Questions
Distributed parameter models such as MIKE SHE can simulate the spatial variability in land cover and changes in ecohydrological processes quiet effectively.Whereas the sensitivity analyses provide some guide to the validity of the simulations, model calibration is complex due to the multiple feedbacks between model parameters and their distributed and temporally variable outputs [82].As rainfall and hillslope hydrological responses are inherently non-linear [83], data must be mapped at hillslope scales in a manner that accurately captures their spatial heterogeneity.However, model development and validation are often constrained by the availability of data and parameter approximations [84].
The temporal resolution of data will also affect the response of the catchment to changes in rainfall intensities, including hydrological responses to land change and configuration [85,86].The simulations presented here were based on a week long period of thunderstorms preceded by a dry period.The reasoning behind this was that regional projections suggest increased drought and up to 40% greater rainfall intensity by 2030 [4].This may have affected some hydrological responses, as infiltration capacity and macropore flow connectivity at hillslope scales display strong nonlinear relationships based on antecedent wetness [87,88].
The depth to an impermeable layer is a significant variable determining the proportioning of overland flow in this study, so this factor ideally would need to be accurately mapped across a hillslope.In steeper terrain this would include mapping of rock fragments, as these also affect infiltration [89].Soil water contents at the upper and lower limits also need careful approximation, as they are likely to be highly heterogeneous at hillslope scales.Along with SAT, data on the type and depth of plant litter is needed to generate estimates of detention storage and surface roughness.Simulations may also be improved by adding bypass flows, such as soil biopores and preferential flow pathways [90].All of these data, however, are costly to obtain by field sampling.
Further improvements to the design of IVB can be gained through modelling or applied experimentation.This might include different configurations, the use of multiple species with different rooting architecture and water use efficiencies, lower competitiveness with pasture or crop species, and the production of different litter types and biomass which will effect detention storage.Simulating alternate configurations would ascertain if the 5% slope gradient is more important than their orientation in relation to topography.It is likely that for alternate slope gradients, soils and vegetation combinations, the ecohydrological response would show considerable variation.In the case of the IVB modelled in this study, the findings and knowledge generated have been directly applied to optimise a paired catchment study that is currently underway in steep grazing lands of Southeast Queensland, Australia.

Conclusions
The model results presented in this simulation study underlies the potential for IVB to retain stormwater runoff in high-slope landscapes and under moderate to high intensity rainfall events.The results reveal four key findings:  At the hillslope scale, simulations predicted that water depth increased beneath IVB but decreased in inter-band spaces.The actual area of hillslope where more than 10 mm of runoff was retained was increased by 22% in IVB compared to the pasture land cover, or an 11% increase in hillslope area external to IVB where extra runoff was retained. At the sub-catchment scale, simulations predicted that landscapes with IVB would have 23% more infiltration and 7% less discharge than landscapes without these bands.There was a slight lengthening of base flow in the channel. Higher SAT levels and greater soil depth were the primary drivers for differences in ecohydrological functioning between pasture and IVB, with secondary functional differences including increased surface roughness, detention storage and infiltration rate. Information on the effects of each parameter was used to optimise the selection of instrumentation for a current paired catchment field trial in Southeast Queensland.The findings may also serve as a guide for establishing similar adaptations of IVB in other regions where protection of hillslopes and catchment systems against greater rainfall intensity may be required.

Figure 1 .
Figure 1.Conceptual model of the integrated vegetation bands (IVB) design and general ecohydrological effects.The inset shows the hypothesised difference in partitioning of overland flow on hillslopes between (A) pasture, or (B) IVB.

Figure 2 .
Figure 2. Location of Maronghi Creek sub-catchment showing contours, fluvial system and the existing land cover onto which IVB were applied.Note the inset which is comprised of a central ridge was selected for closer analysis of hillslope dynamics, including the two transects denoted as A-B and A-C (Coordinates are GDA 1994 MGA Zone 56).
8 mm (A-C), whereas IVB held 20 mm and 19 mm (A-B) and 18 mm and 19 mm (A-C).These values show IVB holding up to 286% greater OF depth in transect A-B and 238% in A-C when soils are saturated.

Figure 4 .
Figure 4. Two hillslope profiles showing overland flow (OF) depth for two transects A-B (a) and A-C (b) at two time-steps (T362 and T639).The two transects start at a common cell (A) located on a central ridge of the sub-catchment (see Figure 2), with the trajectory of each at right angles down the steepest gradient to points B and C. TS362 was at the beginning of the first thunderstorm event when the soil was dry (12:15 PM 7 February 1999), whereas TS639 was taken when the catchment was saturated (4.15 PM 8 February 1999).The location of IVB along the transects are displayed as shaded bars.The most obvious indicator is the much greater depth of overland flow before and within the IVB, and the slight decrease in depth immediately down-slope.This indicates increased water retention in the IVB.

Figure 5 .
Figure 5. Raster model showing OF distribution between pasture (top) and IVB (bottom) simulations.The direction and velocity are related to vector orientation and length respectively.Note the significant reduction in overland flow depth beneath the IVB, and the slight redistribution toward the central ridge compared to the pasture simulation.

Figure 6 .
Figure 6.Raster model showing change in mean infiltration (mm/day) under IVB (NB: the location is the same as Figure 5).

Figure 7 .
Figure 7.A time-series OF for pasture (left) and IVB (right) models.The time-steps are in series from top to bottom showing changes in overland flow every 24 h beginning at 4:15 pm 7 February 1999.The major changes are increased flow depth in IVB themselves, as well as a broader spatial distribution and longer duration of flow in the IVB model.

Figure 8 .
Figure 8. Scatter plots showing the difference (mm) between pasture and IVB simulation models in terms of evapotranspiration (ET), infiltration (INF) and boundary outflow (Q) water balance estimates.Five parameters had varying degrees of influence on water balance estimates, in order of significance these were: (a) soil water content at saturation (SAT); (b) depth to an impermeable layer (ACQ); (c) detention storage (DS); (d) surface roughness (M); and (e) infiltration rate (Ks).

Table 1 .
Simulation variables and parameter estimates for two land covers.
Notes: in this series day 1 begins at 12:00 pm 7 February 1999, and each day thereafter is taken at 12:00 pm.* denotes change is negligible.