Near-Bed Monitoring of Suspended Sediment during a Major Flood Event Highlights Deficiencies in Existing Event-Loading Estimates

Rates of fluvial sediment discharge are notoriously difficult to quantify, particularly during major flood events. Measurements are typically undertaken using event stations requiring large capital investment, and the high cost tends to reduce the spatial coverage of monitoring sites. This study aimed to characterise the near-bed suspended sediment dynamics during a major flood event using a low-cost approach. Monitoring nodes consisted of a total suspended sediment (TSS) logger, a single stage sampler, and a time-lapse camera for a total cost of less than US$420. Seven nodes were deployed across an elevation gradient on the stream bank of Laidley Creek, Queensland, Australia, and two of these nodes successfully characterised the near-bed suspended sediment dynamics across a major flood event. Near-bed TSS concentrations were closely related to stream flow, with the contribution of suspended bed material dominating the total suspended load during peak flows. Observed TSS concentrations were orders of magnitude higher than historical monitoring data for this site collected using the State government event station. This difference was attributed to the event station pump inlet screening the suspended bed material prior to sample collection. The ‘first flush’ phenomenon was detected and attributed to a local resuspension of muddy crusts immediately upstream of the study site. This low-cost approach will provide an important addition to the existing monitoring of fluvial sediment discharge during flood events.


Introduction
Quantifying fluvial suspended sediment discharge is critical to understanding the sustainability of agricultural practice, the water quality of downstream systems, and carbon export from terrestrial systems [1,2].Total suspended sediment loads carried during high flow events are a combination of washload and suspended bed material.Washload is usually dominated by silt and clay particles (<63 µm) and its vertical concentration profile is relatively constant.In contrast, suspended bed material is enriched in coarser sand particles and its vertical profiles are characterised by greatly increased concentrations near the sediment bed.To quantify downstream sediment loading rates, it is important to characterise both components of total suspended loads [3].
Erosion of river channels during major flow events represents the dominant erosion process in the cleared catchments of South East Queensland, Australia [4,5].As a result of wide-scale clearing of upper catchments for agricultural land use, increases in fluvial sediment loading are evident in receiving water bodies [6].This is exemplified in the receiving coastal water body of Moreton Bay, Queensland, where areal coverage by terrestrial muds has increased by more than 50% in less than 45 years [7].Quantifying sediment discharge from these catchments is challenging due to the irregular nature of rainfall patterns and resultant flashy, high flow events as well as the large spatial coverage required to effectively monitor the 21,000 km 2 Moreton Bay catchment.
Existing event monitoring approaches focus on capturing samples across the hydrograph and on the post-processing of samples after collection.There is little consideration given to the sampling inlet design and placement, which is a critical consideration given the high gradients in total suspended sediment (TSS) concentrations that occur in the near-bed zone.There is a need to develop monitoring systems to capture the near-bed gradients in TSS concentrations in order to better understand rates of sediment transport in fluvial systems.One major limitation of current event monitoring stations is the high costs associated with the installation and maintenance of systems: the capital costs associated with establishing a new station are in excess of US$75,000.This is primarily due to the need to ensure that sample collection occurs across major flood events and, therefore, monitoring stations require large investment in infrastructure to house the pump inlet and sample collection and storage facilities (Figure 1).
Water 2018, 10, x FOR PEER REVIEW 2 of 17 catchments for agricultural land use, increases in fluvial sediment loading are evident in receiving water bodies [6].This is exemplified in the receiving coastal water body of Moreton Bay, Queensland, where areal coverage by terrestrial muds has increased by more than 50% in less than 45 years [7].
Quantifying sediment discharge from these catchments is challenging due to the irregular nature of rainfall patterns and resultant flashy, high flow events as well as the large spatial coverage required to effectively monitor the 21,000 km 2 Moreton Bay catchment.
Existing event monitoring approaches focus on capturing samples across the hydrograph and on the post-processing of samples after collection.There is little consideration given to the sampling inlet design and placement, which is a critical consideration given the high gradients in total suspended sediment (TSS) concentrations that occur in the near-bed zone.There is a need to develop monitoring systems to capture the near-bed gradients in TSS concentrations in order to better understand rates of sediment transport in fluvial systems.One major limitation of current event monitoring stations is the high costs associated with the installation and maintenance of systems: the capital costs associated with establishing a new station are in excess of US$75,000.This is primarily due to the need to ensure that sample collection occurs across major flood events and, therefore, monitoring stations require large investment in infrastructure to house the pump inlet and sample collection and storage facilities (Figure 1).This study monitored the near-bed suspended sediment concentrations during a major flood event in a channel erosion hotspot.This addresses a gap in the knowledge of existing event-monitoring programs, where it is unclear where the primary area of uncertainty lies in estimating sediment loading during flood events.Does it lie in capturing changes in sediment concentration at a single point across the hydrograph, or does it lie in characterising the near-bed concentration gradient?To provide for a balance between temporal and spatial sampling, we utilized commercially available washing machine turbidity sensors as this allows for a truly low-cost approach.This study monitored the near-bed suspended sediment concentrations during a major flood event in a channel erosion hotspot.This addresses a gap in the knowledge of existing event-monitoring programs, where it is unclear where the primary area of uncertainty lies in estimating sediment loading during flood events.Does it lie in capturing changes in sediment concentration at a single point across the hydrograph, or does it lie in characterising the near-bed concentration gradient?To provide for a balance between temporal and spatial sampling, we utilized commercially available washing machine turbidity sensors as this allows for a truly low-cost approach.

Study Site Description
Laidley Creek is located in South East Queensland, Australia (Figure 2A), a sub-tropical region prone to large, irregular rainfall events [8].The creek has a catchment surface area of 462 km 2 and runs for approximately 60 km, with much of the course through highly erodible alluvium plains.The original vegetation cover was dominated by dry sclerophyll forest and the riparian vegetation cover is typically mixed eucalypt, casuarina, and wattle forest [9].However, the catchment is highly modified, with a large-scale clearing of forest trees in the late 1800s resulting in present day landuse dominated by grazing and high intensity agriculture [4].It is a sub-catchment of the larger Lockyer catchment, and previous studies of both the Laidley Creek and Lockyer Creek systems have clearly demonstrated that the erosion processes are dominated by gulley and channel erosion [10,11].The creek banks are highly degraded, with limited riparian vegetation remaining and very steep channel profiles.During flood events, large sections of the channel are eroded, resulting in the loss of huge quantities of sediment from the upper catchments.The finer sediment fractions (silts and clays) of this material are transported several hundred kilometers downstream to Moreton Bay [4].The Lockyer catchment is part of the Brisbane River system, which drains over 13,000 km 2 into the coastal embayment, Moreton Bay (Figure 2B).The gulley and channel erosion of the upper catchments has resulted in rapid infilling of the Moreton Bay with terrestrial silts and clays [6,7].
Water 2018, 10, x FOR PEER REVIEW 3 of 17

Study Site Description
Laidley Creek is located in South East Queensland, Australia (Figure 2A), a sub-tropical region prone to large, irregular rainfall events [8].The creek has a catchment surface area of 462 km 2 and runs for approximately 60 km, with much of the course through highly erodible alluvium plains.The original vegetation cover was dominated by dry sclerophyll forest and the riparian vegetation cover is typically mixed eucalypt, casuarina, and wattle forest [9].However, the catchment is highly modified, with a large-scale clearing of forest trees in the late 1800s resulting in present day landuse dominated by grazing and high intensity agriculture [4].It is a sub-catchment of the larger Lockyer catchment, and previous studies of both the Laidley Creek and Lockyer Creek systems have clearly demonstrated that the erosion processes are dominated by gulley and channel erosion [10,11].The creek banks are highly degraded, with limited riparian vegetation remaining and very steep channel profiles.During flood events, large sections of the channel are eroded, resulting in the loss of huge quantities of sediment from the upper catchments.The finer sediment fractions (silts and clays) of this material are transported several hundred kilometers downstream to Moreton Bay [4].The Lockyer catchment is part of the Brisbane River system, which drains over 13,000 km 2 into the coastal embayment, Moreton Bay (Figure 2B).The gulley and channel erosion of the upper catchments has resulted in rapid infilling of the Moreton Bay with terrestrial silts and clays [6,7].The deployment of a low-cost monitoring network was undertaken upstream of the Beckman Rd Bridge across Laidley Creek at Mulgowie, Queensland, Australia (Figure 2C).Initial aerial surveys of the study sites were undertaken in July 2016, shortly after extensive streambank works were undertaken to rehabilitate the western bank (Figure S1).The long-term streamflow monitoring at Mulgowie gauging station (1967-present) recorded a median of 0.08 m 3 /s, a mean of 0.83 m 3 /s, and a The deployment of a low-cost monitoring network was undertaken upstream of the Beckman Rd Bridge across Laidley Creek at Mulgowie, Queensland, Australia (Figure 2C).Initial aerial surveys of the study sites were undertaken in July 2016, shortly after extensive streambank works were undertaken to rehabilitate the western bank (Figure S1).The long-term streamflow monitoring at Mulgowie gauging station (1967-present) recorded a median of 0.08 m 3 /s, a mean of 0.83 m 3 /s, and a maximum of 349 m 3 /s.These highly variable flow rates indicate that the system is event-driven, with 90% of the annual flow occurring during the summer rainfall period [12].The low-cost monitoring network was deployed between 29 March 2017 and 1 April 2017 to capture a major flood event resulting from the passage of Ex-Tropical Cyclone Debbie across the study region.On 30 March 2017, over 170 mm of rainfall was recorded in 24 h at Laidley Creek at the Stony Creek Rd monitoring station located 16 km upstream of the study site [13], and the peak water level of 9.35 m was the second-highest on record, just below the maximum of 9.37 m in January 2013 [12].

Field Monitoring
The low-cost monitoring of TSS concentration was undertaken using custom-built sensors as described in Appendix A, and each sensor has a total parts cost under US$150 (Table A1, Appendix A).Raw sensor output was converted directly to TSS concentration following the methodology in (Figure A5) [14].In addition to the low-cost sensors, single stage samplers were deployed in parallel to collect discrete samples for analysis of suspended sediment particle size distribution and total nutrient concentrations.Time-lapse cameras (Panda UV535 trail camera, UOVision Australia, Oakleigh South, VIC, Australia) were utilized to record changes in water level and colour at 5-min intervals.The total cost of a single monitoring node, consisting of one 1.8 m star picket (US$9), one low-cost sensor (US$148), one single stage sample (US$35), and one time-lapse camera (US$190), was less than US$400.
Monitoring nodes were deployed along an elevation gradient up the streambank to capture changes in sediment loading across different stages of the hydrograph (Figure 3).The deployment configurations were mounted on star pickets driven into the bank and followed the methodology outlined in [3].A total of seven monitoring nodes were deployed across the site with sensor heights ranging from 0.2 m to 9.5 m above local water level, which were then referenced to the State Government gauging station (Figure 2C).Local stream height changes were recorded using a non-vented underwater pressure transducer (Level TROLL 400 Data Logger, In-situ, Fort Collins, CO, USA), which was deployed immediately adjacent to the lowest monitoring node.The pressure transducer was anchored in a 0.3 m water depth using a star picket and the water level was logged at 1-min intervals for the duration of the deployment.During zero flow conditions, surface and sub-surface (10 cm depth) sediment samples were collected within the creek bed for analysis of particle size distribution.The surface extent of bed material deposits was estimated immediately upstream of the study site by visual inspection of a series of holes dug across the deposit at 5-m intervals.
Water 2018, 10, x FOR PEER REVIEW 4 of 17 maximum of 349 m 3 /s.These highly variable flow rates indicate that the system is event-driven, with 90% of the annual flow occurring during the summer rainfall period [12].The low-cost monitoring network was deployed between 29 March 2017 and 1 April 2017 to capture a major flood event resulting from the passage of Ex-Tropical Cyclone Debbie across the study region.On 30 March 2017, over 170 mm of rainfall was recorded in 24 h at Laidley Creek at the Stony Creek Rd monitoring station located 16 km upstream of the study site [13], and the peak water level of 9.35 m was the second-highest on record, just below the maximum of 9.37 m in January 2013 [12].

Field Monitoring
The low-cost monitoring of TSS concentration was undertaken using custom-built sensors as described in Appendix A, and each sensor has a total parts cost under US$150 (Table A1, Appendix A).Raw sensor output was converted directly to TSS concentration following the methodology in (Figure A5) [14].In addition to the low-cost sensors, single stage samplers were deployed in parallel to collect discrete samples for analysis of suspended sediment particle size distribution and total nutrient concentrations.Time-lapse cameras (Panda UV535 trail camera, UOVision Australia, Oakleigh South, VIC, Australia) were utilized to record changes in water level and colour at 5-min intervals.The total cost of a single monitoring node, consisting of one 1.8 m star picket (US$9), one low-cost sensor (US$148), one single stage sample (US$35), and one time-lapse camera (US$190), was less than US$400.
Monitoring nodes were deployed along an elevation gradient up the streambank to capture changes in sediment loading across different stages of the hydrograph (Figure 3).The deployment configurations were mounted on star pickets driven into the bank and followed the methodology outlined in [3].A total of seven monitoring nodes were deployed across the site with sensor heights ranging from 0.2 m to 9.5 m above local water level, which were then referenced to the State Government gauging station (Figure 2C).Local stream height changes were recorded using a nonvented underwater pressure transducer (Level TROLL 400 Data Logger, In-situ, Fort Collins, CO, USA), which was deployed immediately adjacent to the lowest monitoring node.The pressure transducer was anchored in a 0.3 m water depth using a star picket and the water level was logged at 1-min intervals for the duration of the deployment.During zero flow conditions, surface and sub-surface (10 cm depth) sediment samples were collected within the creek bed for analysis of particle size distribution.The surface extent of bed material deposits was estimated immediately upstream of the study site by visual inspection of a series of holes dug across the deposit at 5-m intervals.

Laboratory Analysis
Water samples from single stage samplers were retrieved within 24 h of peak flows.A total nitrogen and total phosphorus analysis was undertaken on these samples using standard methods [15] by the National Association of Testing Authorities, Australia (NATA)-accredited laboratory Environmental Analysis Laboratory at the Southern Cross University.To better understand the relative contribution of washload to downstream nutrient loading, subsamples from two single stage samplers were sieved through a 63 µm sieve to remove the sand fraction and the filtrate was analysed for total phosphorus and total nitrogen.
Samples of suspended and bed material for analysis of particle size distribution were sonicated for one minute prior to measurement to separate aggregated particles, allowing the true particle size to be calculated [16].Particle sizing was undertaken using laser diffraction (Malvern Mastersizer 2000E, Malvern Instruments Ltd., Malvern, UK), and five measurement cycles were performed for each sample.Sediments were classified according to the Wentworth scale, with mud fractions defined as sediment particles with a diameter of <63 µm [17].
Given the extreme TSS concentrations as well as the obvious sand fractions in the single stage samples collected during high flows, a whole-water test was undertaken to determine how representative the TSS particle size distribution estimated by laser diffraction was.The issue of obtaining representative sub-samples for particle size distribution or TSS concentration from samples with significant sand fractions has been well-established [18].The whole-water test followed the methodology outlined in [19], where a 12:10 p.m. single stage sample was wet sieved through a 63 µm micron sieve and the sand and mud fractions were dried at 105 • C. The relative sand and mud fractions were determined gravimetrically and expressed as a percentage.TSS concentrations were determined gravimetrically, where samples were filtered using a dried pre-weighed 0.7 µm filter (Whatman glass microfiber filter GF/F grade) and then dried at 105 • C to a constant weight.Filters were reweighed and the TSS concentration (mg/L) was then calculated from the filter weight difference and normalized to the sample volume.

Imagery Analysis
To obtain high-resolution imagery of the study site before and after a flood event, a series of drone flights were undertaken in July 2016, one day before the flood event (29 March 2017) and one day after the flood event (31 March 2017).The DJI Phantom 3 Pro platform was utilized with the camera (DJI FC300X with a 12 megapixel RGB sensor) in a nadiral orientation for mapping flights.Georeferenced imagery was downloaded from the drone and images selected for photogrammetric data processing.Data processing was undertaken using dedicated processing software (Pix4D Promapper, Pix4D SA, Lausanne, Switzerland), and the resulting mosaics were visualized using the geographical information system ArcGIS (Version 10.3, ESRI Inc., Redlands, CA, USA).

Time-Lapse and Drone Imagery
High-resolution drone mosaics before and after the event provided a detailed picture of the impact of this flood event, with clear signs of erosion on the natural eastern bank as well as the levee breach on the rehabilitated western bank (Figure 4A,B; Video S1; Video S2).Erosion of the natural bank occurred along at least 250 m of the creek channel upstream from the study site, with slope failures resulting in up to 3 m (horizontal) and 4 m (vertical) of bank loss (Figure 4B; Video S3).The levee breach on the rehabilitated bank highlighted the extreme nature of the flood event, as this was designed to withstand major floods (Figure 4B; Video S4).Oblique imagery from the site highlighted the contrast between the natural bank and the rehabilitated bank (Figure 5).Despite being on the inside bend, the natural bank experienced substantial erosion, whilst the rehabilitated bank on the outside of the bend was largely stable (Figure 5).This further highlights the potential for streambank rehabilitation to reduce sediment loading from channel erosion, a critical sediment loading source to Moreton Bay [6,11].Oblique imagery from the site highlighted the contrast between the natural bank and the rehabilitated bank (Figure 5).Despite being on the inside bend, the natural bank experienced substantial erosion, whilst the rehabilitated bank on the outside of the bend was largely stable (Figure 5).This further highlights the potential for streambank rehabilitation to reduce sediment loading from channel erosion, a critical sediment loading source to Moreton Bay [6,11].Oblique imagery from the site highlighted the contrast between the natural bank and the rehabilitated bank (Figure 5).Despite being on the inside bend, the natural bank experienced substantial erosion, whilst the rehabilitated bank on the outside of the bend was largely stable (Figure 5).This further highlights the potential for streambank rehabilitation to reduce sediment loading from channel erosion, a critical sediment loading source to Moreton Bay [6,11].

Suspended Sediment Dynamics
The stream flow across this major flood event consisted of a relatively small initial peak of 86 m 3 /s occurring at 1:00 p.m. followed by the far larger peak of over 320 m 3 /s occurring after 8:00 p.m. on 30 March 2017 (Figure 6A).Due to the extremity of the flood event, only two of the seven sensors deployed were recovered (Figure S2).Logged TSS concentrations were strongly correlated to streamflow across the event; however, the timing of peaks in the TSS and streamflow differed slightly.The initial flows were characterised by a peak in the TSS that preceded the peak flow; this was followed by a lag in the TSS peak with the much higher secondary peak in streamflow (Figure 6A).Changes in TSS concentrations were mirrored by both loggers recovered, although logger 2 lost its sensor just prior to the major peak flow.The majority of suspended sediment transport occurred at high flow rates, with TSS concentrations recorded above 80 g/L by both TSS loggers as well as by a single stage sampler.A large discrepancy was observed between the TSS concentration of the 7 p.m. single stage sample (384 mg/L) and that of the logged values (57 and 59 mg/L) recorded at this time (Figure 6A).Overall, these data are in contrast to the long-term event monitoring undertaken by the State Government, where event mean concentrations from 203 samples collected between 2002 and 2015 averaged 0.272 g/L and ranged between 0.094 and 0.762 g/L [20,21].Whilst these data include a number of small flow events, there were at least two major floods (in 2011 and 2013) that occurred in this time period and that recorded similar peaks to this 2017 event.
Water 2018, 10, x FOR PEER REVIEW 7 of 17

Suspended Sediment Dynamics
The stream flow across this major flood event consisted of a relatively small initial peak of 86 m 3 /s occurring at 1:00 p.m. followed by the far larger peak of over 320 m 3 /s occurring after 8:00 p.m. on 30 March 2017 (Figure 6A).Due to the extremity of the flood event, only two of the seven sensors deployed were recovered (Figure S2).Logged TSS concentrations were strongly correlated to streamflow across the event; however, the timing of peaks in the TSS and streamflow differed slightly.The initial flows were characterised by a peak in the TSS that preceded the peak flow; this was followed by a lag in the TSS peak with the much higher secondary peak in streamflow (Figure 6A).Changes in TSS concentrations were mirrored by both loggers recovered, although logger 2 lost its sensor just prior to the major peak flow.The majority of suspended sediment transport occurred at high flow rates, with TSS concentrations recorded above 80 g/L by both TSS loggers as well as by a single stage sampler.A large discrepancy was observed between the TSS concentration of the 7 p.m. single stage sample (384 mg/L) and that of the logged values (57 and 59 mg/L) recorded at this time (Figure 6A).Overall, these data are in contrast to the long-term event monitoring undertaken by the State Government, where event mean concentrations from 203 samples collected between 2002 and 2015 averaged 0.272 g/L and ranged between 0.094 and 0.762 g/L [20,21].Whilst these data include a number of small flow events, there were at least two major floods (in 2011 and 2013) that occurred in this time period and that recorded similar peaks to this 2017 event.The suspended sediment particle size distribution across the event clearly showed an increasing contribution from sand particles during higher flows (Figure 6B).Relatively low flows were almost exclusively composed of fine silt and clay sediment particles (D10% 2.3 µm; D50% 14.8 µm; D90% 50.4 µm).The increasing contribution from sandy particles is observed in the D90% particle diameter, which increased from 50.4 µm (coarse silt) at 11:25 a.m. to 92.5 µm (very fine sand) at 12:10 p.m. and 448 µm at 7:00 p.m. (medium sand).Particle sizing using laser diffraction heavily favoured the fine particles.The suspended sediment particle size distribution across the event clearly showed an increasing contribution from sand particles during higher flows (Figure 6B).Relatively low flows were almost exclusively composed of fine silt and clay sediment particles (D10% 2.3 µm; D50% 14.8 µm; D90% 50.4 µm).The increasing contribution from sandy particles is observed in the D90% particle diameter, which increased from 50.4 µm (coarse silt) at 11:25 a.m. to 92.5 µm (very fine sand) at 12:10 p.m. and 448 µm at 7:00 p.m. (medium sand).Particle sizing using laser diffraction heavily favoured the fine particles.At 12:10 p.m., the mud fraction was estimated to be 80%; however, the sieve analysis of the whole sample estimated the mud fraction to be 40%.This is likely due to the difficulty associated with taking a truly representative sample for laser diffraction with slurries containing such a high sand fraction.The issue of obtaining representative sub-samples applied to quantifying TSS concentrations in the single stage as well, with sub-samples approximately 10% lower than whole-sample estimates.

First Flush Phenomenon
The first flush is often defined as the disproportionate increase in concentration of soluble and particulate materials during the rising stage of a flow event [21].The initial flows observed in this study underwent a first flush phenomenon, as water colour rapidly changed from transparent to highly coloured with a very small rise (approximately 20 cm) in water level (Figure 7; Video S4).The cause of this rapid colour change was a large increase in suspended clay and silt particles (Figure 6B).The water discoloration persisted with the subsequent flows across the whole event (Video S5).Together with the rapid increase in TSS preceding the initial peak (Figure 6A), these data suggest the importance of this 'first flush' phenomenon and that even low streamflow events may be important contributors to downstream sediment loading.
Water 2018, 10, x FOR PEER REVIEW 8 of 17 At 12:10 p.m., the mud fraction was estimated to be 80%; however, the sieve analysis of the whole sample estimated the mud fraction to be 40%.This is likely due to the difficulty associated with taking a truly representative sample for laser diffraction with slurries containing such a high sand fraction.The issue of obtaining representative sub-samples applied to quantifying TSS concentrations in the single stage as well, with sub-samples approximately 10% lower than whole-sample estimates.

First Flush Phenomenon
The first flush is often defined as the disproportionate increase in concentration of soluble and particulate materials during the rising stage of a flow event [21].The initial flows observed in this study underwent a first flush phenomenon, as water colour rapidly changed from transparent to highly coloured with a very small rise (approximately 20 cm) in water level (Figure 7; Video S4).The cause of this rapid colour change was a large increase in suspended clay and silt particles (Figure 6B).The water discoloration persisted with the subsequent flows across the whole event (Video S5).Together with the rapid increase in TSS preceding the initial peak (Figure 6A), these data suggest the importance of this 'first flush' phenomenon and that even low streamflow events may be important contributors to downstream sediment loading.A large sandy deposit was observed within the creek channel immediately upstream of the study site with a surface area of approximately 174 m 2 and a minimum deposit depth of 70 cm (Figure 8A).The entire surface of this deposit was covered in a thin crust between 1 and 3 cm deep, which was greatly enriched in silt and clay fractions (50% mud) relative to the underlying sandy bed material (10% mud) (Figure 8B).
Total nutrient concentrations were highly elevated in both the single stage samples collected at relatively low flows and relatively high flows (Table 1).In addition, there was relatively little change between total nutrient concentrations of total and sieved samples, which suggests that much of the nitrogen and phosphorus species was associated with the washload (Table 1).The association of phosphorus with fine sediment particles in washload has been observed during high flow events in this region [22].Total nitrogen concentrations during the high flows at 7:00 p.m. were lower concentrations after sieving (Table 1).In this region, total nitrogen is dominated by the particulate organic nitrogen fractions during high flows [23], and larger fragments may have been trapped by the sample sieving, reducing the concentration of total nitrogen.Long-term event monitoring total nutrient concentrations from the State Government event station were approximately one order of magnitude lower than those found in the single stage samples (Table 1).A large sandy deposit was observed within the creek channel immediately upstream of the study site with a surface area of approximately 174 m 2 and a minimum deposit depth of 70 cm (Figure 8A).The entire surface of this deposit was covered in a thin crust between 1 and 3 cm deep, which was greatly enriched in silt and clay fractions (50% mud) relative to the underlying sandy bed material (10% mud) (Figure 8B).
Total nutrient concentrations were highly elevated in both the single stage samples collected at relatively low flows and relatively high flows (Table 1).In addition, there was relatively little change between total nutrient concentrations of total and sieved samples, which suggests that much of the nitrogen and phosphorus species was associated with the washload (Table 1).The association of phosphorus with fine sediment particles in washload has been observed during high flow events in this region [22].Total nitrogen concentrations during the high flows at 7:00 p.m. were lower concentrations after sieving (Table 1).In this region, total nitrogen is dominated by the particulate organic nitrogen fractions during high flows [23], and larger fragments may have been trapped by the sample sieving, reducing the concentration of total nitrogen.Long-term event monitoring total nutrient concentrations from the State Government event station were approximately one order of magnitude lower than those found in the single stage samples (Table 1).Table 1.Total nutrient concentrations from two single stage samplers during the flow event collected at 12:10 p.m. and 7:00 p.m. on 30 March 2017.The washload refers to samples that passed through a 63micron sieve.Total phosphorus (TP) and total nitrogen (TN) concentration units in mg/L.

Discussion
This study successfully captured changes in near-bed TSS concentrations across a major flood event using a low-cost approach.Changes in suspended sediment particle size distribution clearly showed the increasing contribution of sandy fractions to TSS as the flood event developed.Total nitrogen and phosphorus fractions were closely associated with washload, suggesting that the majority of nutrient load was exported from this system to downstream receiving water bodies.The near-bed TSS concentrations were in excess of that required to hinder settling (5 to 10 g/L) of suspended sand Table 1.Total nutrient concentrations from two single stage samplers during the flow event collected at 12:10 p.m. and 7:00 p.m. on 30 March 2017.The washload refers to samples that passed through a 63-micron sieve.Total phosphorus (TP) and total nitrogen (TN) concentration units in mg/L.

Discussion
This study successfully captured changes in near-bed TSS concentrations across a major flood event using a low-cost approach.Changes in suspended sediment particle size distribution clearly showed the increasing contribution of sandy fractions to TSS as the flood event developed.Total nitrogen and phosphorus fractions were closely associated with washload, suggesting that the majority of nutrient load was exported from this system to downstream receiving water bodies.
The near-bed TSS concentrations were in excess of that required to hinder settling (5 to 10 g/L) of suspended sand particles [24]; this would reduce the settling rates of these particles and, therefore, increase their time in suspension and subsequent transport distances downstream [25].
We suggest that the first flush phenomenon observed in this study is a result of local resuspension of muddy crusts immediately upstream of the monitoring site.Whilst this resulted in relatively high TSS concentrations in the near-bed, the higher flows experienced subsequently entrained the underlying sandy bed, which resulted in the extreme TSS concentrations during peak flows.The widespread channel erosion occurring throughout the creek system would likely result in a continuous sediment supply across the event.The finer mud fractions would settle more slowly than the sandy particles as the flood event recedes, reforming the muddy crusts observed in the stream bed.This behavior has been clearly demonstrated in laboratory experiments examining sand-mud suspensions [26].
The discrepancy in TSS concentration between the 7 p.m. single stage sample and the logged values was likely due to the shroud structure surrounding the low-cost sensors (Figure A1).This shroud would shield the sensor face from direct water flows, slowing the near-field water velocities and allowing larger sand particles to rapidly settle in this zone [27].The net effect would reduce the concentration of sand particles passing the sensor face and greatly reduce the observed TSS concentrations [28].The lower flows during the initial stage of the flood event were dominated by silt and clay fractions which would not settle as quickly and showed good agreement between single stage measurements and logged values.In addition, the extremely high TSS concentrations during peak flows resulted in the sensors approaching their upper detection limits (Figure A5) [14].These limitations will form the basis of future research to improve the shroud design and sensor signal-to-noise ratio at extreme TSS concentrations.
The TSS and total nutrient concentrations observed in this study were greatly elevated compared to the long-term monitoring data from the State Government event station located approximately 250 m downstream.The pump inlet for the event station was located on the stream bed and would have been exposed to similarly high concentrations of suspended bed material and washload as the low-cost sensors.One explanation for this discrepancy lies in the protective cover surrounding the pump inlet (Figure 1C) as well as the 500-µm screen on the inlet pipe.This would greatly reduce water flow in the immediate vicinity of the inlet as well as physically filter out much of the coarser fraction of the suspended bed material prior to collection.Sediment particle fractions that were able to pass through these would likely be limited to very fine particles, which would lead to significant undersampling of the actual sediment load [27].This highlights an area of uncertainty in current event-monitoring practices; standardisation of the inlet structure and its elevation above the bed should be considered in standard operating procedures in order to ensure comparable results are obtained.For instance, sample collection systems that integrate over the entire water column as described in [29] should be considered.
These findings have direct relevance to understanding and analysing sediment budgets, e.g., the contribution of catchment-wide sediment yield, erosion from upstream channel regions or floodplain sources, and changes to the volume and residence time within the system [30].Particularly in larger catchment areas, the quantification of these processes is complex [31].Careful and extensive monitoring is required [32,33], which poses a logistical and financial challenge in many river systems from mountainous streams to the ocean.In addition, the variation of suspended sediment concentration during flood events as reported here and in numerous other studies [34][35][36][37][38][39] adds to the difficulty and uncertainty of an efficient and robust monitoring strategy.
The development of a truly low-cost approach to monitoring sediment dynamics during flood events will provide an opportunity for spatially intense regional monitoring programs as well as a cost-effective method to quantify vertical TSS gradients at single sites.This is of particular importance to guide catchment rehabilitation management actions where it will be critical to separate the effects of regional landuse changes from local stream bed effects.In addition, as climate-induced changes in rainfall patterns increase the likelihood of extreme flood events [40], this low-cost approach will provide an important addition to the existing monitoring of these events.The simplicity of the approach will also provide opportunities globally to deploy these sensors in areas currently lacking any formal monitoring system.

Supplementary Materials:
The following are available online at http://www.mdpi.com/2073-4441/10/2/34/s1, Figure S1 2. Measure the diameter of the waterproof cable glands threaded section (≈19 mm) and drill the corresponding size hole in the centre of the third 50 mm PVC Cap.Clean the edges of the hole, ensuring there is no burr or raised markings.Secure the cable gland in the drilled hole (Figure A2B).Hot glue can be used on the nut of the gland to provide additional support when tightening. 3. Cut a 75 mm length of 50 mm diameter PVC pipe, clean the edges of the cut.4. Cut the appropriate length of cable for the proposed deployment, noting that the length should be 500 mm greater than required.5. Strip 50 mm of outer black coating from both ends of the sensor cable.Then, strip a small amount (≈10 mm) from the red, yellow, and internal black wire covers.Cut the green cable off to the same level at the outer black coating (Figure A3A).6. Twist the ends of the exposed wires separately and apply lead-free solder.7. Remove the black cap from the Amphenol Sensor and then remove the sensor board from the clear housing.Cut the white plastic plug off the sensor board, and then solder the sensor cable wires to the uncovered pins (Figure A3A), covering connections in heat-shrink or electrical tape.Place the sensor board back into the housing and reattach the black cap, ensuring that the sensor cable exits through the hole.Seal the cable gap with hot glue.8. Place the sensor into the glued caps assembled in Step 1, ensuring that the lower section of the sensor is the on the side of shroud.Hot glue the sensor in this position, and slide the pipe down the cable.Glue the lower caps and the pipe together with PVC cement/glue (Figure A3D).

2.
Measure the diameter of the waterproof cable glands threaded section (≈19 mm) and drill the corresponding size hole in the centre of the third 50 mm PVC Cap.Clean the edges of the hole, ensuring there is no burr or raised markings.Secure the cable gland in the drilled hole (Figure A2B).Hot glue can be used on the nut of the gland to provide additional support when tightening.

3.
Cut a 75 mm length of 50 mm diameter PVC pipe, clean the edges of the cut.4.
Cut the appropriate length of cable for the proposed deployment, noting that the length should be 500 mm greater than required.5.
Strip 50 mm of outer black coating from both ends of the sensor cable.Then, strip a small amount (≈10 mm) from the red, yellow, and internal black wire covers.Cut the green cable off to the same level at the outer black coating (Figure A3A).6.
Twist the ends of the exposed wires separately and apply lead-free solder.7.
Remove the black cap from the Amphenol Sensor and then remove the sensor board from the clear housing.Cut the white plastic plug off the sensor board, and then solder the sensor cable wires to the uncovered pins (Figure A3A), covering connections in heat-shrink or electrical tape.Place the sensor board back into the housing and reattach the black cap, ensuring that the sensor cable exits through the hole.Seal the cable gap with hot glue.2. Measure the diameter of the waterproof cable glands threaded section (≈19 mm) and drill the corresponding size hole in the centre of the third 50 mm PVC Cap.Clean the edges of the hole, ensuring there is no burr or raised markings.Secure the cable gland in the drilled hole (Figure A2B).Hot glue can be used on the nut of the gland to provide additional support when tightening. 3. Cut a 75 mm length of 50 mm diameter PVC pipe, clean the edges of the cut.4. Cut the appropriate length of cable for the proposed deployment, noting that the length should be 500 mm greater than required.5. Strip 50 mm of outer black coating from both ends of the sensor cable.Then, strip a small amount (≈10 mm) from the red, yellow, and internal black wire covers.Cut the green cable off to the same level at the outer black coating (Figure A3A).6. Twist the ends of the exposed wires separately and apply lead-free solder.7. Remove the black cap from the Amphenol Sensor and then remove the sensor board from the clear housing.Cut the white plastic plug off the sensor board, and then solder the sensor cable wires to the uncovered pins (Figure A3A), covering connections in heat-shrink or electrical tape.Place the sensor board back into the housing and reattach the black cap, ensuring that the sensor cable exits through the hole.Seal the cable gap with hot glue.8. Place the sensor into the glued caps assembled in Step 1, ensuring that the lower section of the sensor is the on the side of shroud.Hot glue the sensor in this position, and slide the pipe down the cable.Glue the lower caps and the pipe together with PVC cement/glue (Figure A3D).

Figure 1 .
Figure 1.Typical arrangement for current event monitoring stations in the South East Queensland network: (A) Monitoring equipment housed in a large infrastructure facility located above a streambank; (B) Inset showing the location of a pump inlet for water collection; (C) Inset showing the protective cover of the pump inlet relative to the stream bed.The pump inlet has an additional screen (500 µm) within this cover.

Figure 1 .
Figure 1.Typical arrangement for current event monitoring stations in the South East Queensland network: (A) Monitoring equipment housed in a large infrastructure facility located above a streambank; (B) Inset showing the location of a pump inlet for water collection; (C) Inset showing the protective cover of the pump inlet relative to the stream bed.The pump inlet has an additional screen (500 µm) within this cover.

Figure 2 .
Figure 2. Study site location relative to (A) Queensland, Australia; (B) Brisbane River catchment in South East Queensland (sub-catchments indicated in red outlines) and (C) long-term State Government monitoring station: Laidley Creek at Mulgowie gauging station (No. 143209B).Elevation referenced to Australian Height Datum (AHD).

Figure 2 .
Figure 2. Study site location relative to (A) Queensland, Australia; (B) Brisbane River catchment in South East Queensland (sub-catchments indicated in red outlines) and (C) long-term State Government monitoring station: Laidley Creek at Mulgowie gauging station (No. 143209B).Elevation referenced to Australian Height Datum (AHD).

Figure 3 .
Figure 3. Field installation of low-cost total suspended sediment (TSS) sensors and single stage samplers on star pickets.(A) Installation directly adjacent to stream.(B) Upslope installation on a rehabilitated bank 7 m above the stream bed; note that this installation has a time-lapse camera on top of the star picket.

Figure 3 .
Figure 3. Field installation of low-cost total suspended sediment (TSS) sensors and single stage samplers on star pickets.(A) Installation directly adjacent to stream.(B) Upslope installation on a rehabilitated bank 7 m above the stream bed; note that this installation has a time-lapse camera on top of the star picket.

Figure 4 .
Figure 4. High-resolution mosaic imagery (pixel resolution <2 cm) of the study area collected (A) prior to the major flood event on 16 July 2016 and (B) immediately after the event, 31 March 2017, showing the widespread channel erosion on the natural bank and the levee breach on the rehabilitated bank.

Figure 5 .
Figure 5. Drone image facing upstream from the study site clearly showing channel degradation along the far (natural) bank and the relatively intact nature of the near (rehabilitated) bank after the major flood event.

Figure 4 .
Figure 4. High-resolution mosaic imagery (pixel resolution <2 cm) of the study area collected (A) prior to the major flood event on 16 July 2016 and (B) immediately after the event, 31 March 2017, showing the widespread channel erosion on the natural bank and the levee breach on the rehabilitated bank.

Figure 4 .
Figure 4. High-resolution mosaic imagery (pixel resolution <2 cm) of the study area collected (A) prior to the major flood event on 16 July 2016 and (B) immediately after the event, 31 March 2017, showing the widespread channel erosion on the natural bank and the levee breach on the rehabilitated bank.

Figure 5 .
Figure 5. Drone image facing upstream from the study site clearly showing channel degradation along the far (natural) bank and the relatively intact nature of the near (rehabilitated) bank after the major flood event.

Figure 5 .
Figure 5. Drone image facing upstream from the study site clearly showing channel degradation along the far (natural) bank and the relatively intact nature of the near (rehabilitated) bank after the major flood event.

Figure 6 .
Figure 6.(A) TSS measurements recorded from two low-cost sensors and stream flow during the high flow event proceeding the passage of Ex-Tropical Cyclone Debbie.Discrete single stage sample timing and concentrations are indicated with arrows.(B) Particle size distribution for samples collected by single stage samplers on 30 March 2017 during the initial phases of the flood event.

Figure 6 .
Figure 6.(A) TSS measurements recorded from two low-cost sensors and stream flow during the high flow event proceeding the passage of Ex-Tropical Cyclone Debbie.Discrete single stage sample timing and concentrations are indicated with arrows.(B) Particle size distribution for samples collected by single stage samplers on 30 March 2017 during the initial phases of the flood event.

Figure 7 .
Figure 7. Selected time-lapse images recorded during the initial phase of the flood event showing rapid discoloration associated with a 20 cm rise in stream height.

Figure 7 .
Figure 7. Selected time-lapse images recorded during the initial phase of the flood event showing rapid discoloration associated with a 20 cm rise in stream height.

Figure 8 .
Figure 8. (A) Image of creek bed facing upstream of the study site showing muddy crust overlying the large sandy bed deposit.(B) Particle size distribution of the muddy crust as well as the underlying sandy bed.

Figure 8 .
Figure 8. (A) Image of creek bed facing upstream of the study site showing muddy crust overlying the large sandy bed deposit.(B) Particle size distribution of the muddy crust as well as the underlying sandy bed.

:Figure A2 .
Figure A2.(A) Amphenol sensor.(B) Third PVC cap and cable gland.(C) Two caps produced in step 1 for forming the bottom of the sensor assembly.

Figure A3 .
Figure A3.(A) Illustrating the cable soldered on the sensor board, ensuring the solder does not bridge between the sensors.The wires connected as; Red to A/C, Black to K, and Yellow to E. (B) Epoxy-filled sensor housing.The harder epoxy is first to provide a dense base, and hold the cable firmly in place to prevent movement of the soldered joint.(C) Silicon filling is used in the top to provide a more flexible housing fill so the cable can be positioned for the top cap.(D) Completed sensor assembly.

Figure A2 .
Figure A2.(A) Amphenol sensor.(B) Third PVC cap and cable gland.(C) Two caps produced in step 1 for forming the bottom of the sensor assembly.

Figure A2 .
Figure A2.(A) Amphenol sensor.(B) Third PVC cap and cable gland.(C) Two caps produced in step 1 for forming the bottom of the sensor assembly.

Figure A3 .
Figure A3.(A) Illustrating the cable soldered on the sensor board, ensuring the solder does not bridge between the sensors.The wires connected as; Red to A/C, Black to K, and Yellow to E. (B) Epoxy-filled sensor housing.The harder epoxy is first to provide a dense base, and hold the cable firmly in place to prevent movement of the soldered joint.(C) Silicon filling is used in the top to provide a more flexible housing fill so the cable can be positioned for the top cap.(D) Completed sensor assembly.

Figure A3 .
Figure A3.(A) Illustrating the cable soldered on the sensor board, ensuring the solder does not bridge between the sensors.The wires connected as; Red to A/C, Black to K, and Yellow to E. (B) Epoxy-filled sensor housing.The harder epoxy is first to provide a dense base, and hold the cable firmly in place to prevent movement of the soldered joint.(C) Silicon filling is used in the top to provide a more flexible housing fill so the cable can be positioned for the top cap.(D) Completed sensor assembly.
[19][19]data summary from the State water quality monitoring program.Data is from 194 samples collected during events at the Mulgowie event monitoring station from 2002 to 2015.Values refer to the 50th percentile with the 20th to 80th in brackets.
[19]f.[19]data summary from the State water quality monitoring program.Data is from 194 samples collected during events at the Mulgowie event monitoring station from 2002 to 2015.Values refer to the 50th percentile with the 20th to 80th in brackets. Note