Next Article in Journal
Low-Frequency Noise Characteristics of (Al, Ga)As and Ga(As, Bi) Quantum Well Structures for NIR Laser Diodes
Previous Article in Journal
MmWave Physical Layer Network Modeling and Planning for Fixed Wireless Access Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Resolution Estimation of Suspended Solids and Particulate Phosphorus Using Acoustic Devices in a Hydrologically Managed Canal in South Florida, USA

by
Ikechukwu S. Onwuka
1,2,
Leonard J. Scinto
1,2,* and
David C. Fugate
3
1
Department of Earth and Environment, Florida International University, Miami, FL 33199, USA
2
Institute of Environment, Florida International University, Miami, FL 33199, USA
3
Department of Marine and Earth Sciences, Florida Gulf Coast University, 10501 FGCU Blvd. South, Fort Myers, FL 33965, USA
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(4), 2281; https://doi.org/10.3390/s23042281
Submission received: 26 January 2023 / Revised: 13 February 2023 / Accepted: 15 February 2023 / Published: 17 February 2023
(This article belongs to the Section Physical Sensors)

Abstract

:
Conventional methods of measuring total suspended sediments (TSS) and total particulate phosphorus (TPP) are typically low-resolution and miss critical processes that impact their exports in aquatic environments. To create high-resolution TSS and TPP estimates, echo intensity (EI), a biproduct of velocity measurements from acoustic devices, was utilized. An acoustic Doppler velocimeter (ADV) and an acoustic Doppler current profiler (ADCP) were deployed in three locations in the L-29 Canal in South Florida, USA, to obtain estimates near the canal bed and in the water column, respectively. Corrections for transmission losses from the ADCP proved unnecessary due to the low vertical variability in the measured EI. EI calibrations were performed using artificially created TSS obtained from bed sediments (ADV) and gravimetrically measured TSS from water samples that matched the depths and times of the ADCP deployments. The measured TSS values were then analyzed for total phosphorus and converted to TPP estimates. The results showed that high TSS and TPP were caused by the rapid discharge releases typical of managed canals. This work demonstrates that high-resolution estimates are imperative for assessing the effects of such swift hydrologic changes on the potential export of sediments and nutrients to delicate ecosystems downstream.

1. Introduction

The movement of suspended sediments impacts the quality of aquatic systems, including freshwater environments, because they transport particulate nutrients [1,2,3,4]. Phosphorus (P) is one such nutrient [5,6,7], and excess P is responsible for adverse water quality effects, such as eutrophication and harmful algal blooms [8,9]. Due to hydrologic and biogeochemical processes, suspended sediments and their associated nutrients are dynamic in space and time. Therefore, high-resolution measurements are necessary to fully evaluate the conditions that can increase sediment and nutrient loadings to enable the restoration and preservation of impaired aquatic systems [3,10,11,12].
Total suspended sediment (TSS) concentrations are generally determined in two ways, namely through direct and surrogate methods. Direct methods involve the conventional collection of water samples via bottles and pumps, which are then analyzed gravimetrically in the laboratory [10,13,14]. This method is usually expensive and labor-intensive, often requiring costly field sampling techniques and laboratory protocols; intrusive, as delicate flocculated aggregates can easily be broken down during sample retrieval or handling; and, in some instances, sampling is physically impossible or dangerous for field personnel [15,16,17]. Furthermore, the collected samples may not be representative of the range of conditions in aquatic systems due to the low temporal and spatial resolutions of the sampling regimes employed [3,13,16].
Surrogate methods are indirect ways of estimating TSS concentrations without a reliance on water sample collection and can thus be implemented autonomously [13]. Examples of surrogate methods include optical backscatter sensors, transmissometers, laser diffraction, and acoustic devices, which all make nonintrusive high-temporal-resolution estimates of TSS [10,13,18]. Acoustic devices are resistant to biofouling, unlike optical backscatter sensors [19,20], and can also make high-resolution spatial measurements [13,17].
Commonly used acoustic devices include acoustic Doppler velocimeters (ADVs) and acoustic Doppler current profilers (ADCPs). These devices are primarily velocity instruments that measure the frequency shift (Doppler shift) of transmitted sound waves (back) scattered by suspended materials in the water and convert the sound waves or echoes to three-dimensional (north/south, east/west, and vertical) flow velocities [10,13,21]. Acoustic Doppler velocimeters (ADVs) make single-point velocity measurements for a small water volume (measurement volume) at a point approximately 5–10 cm away from the device [13]. Acoustic Doppler current profilers (ADCPs) can further provide velocity by splitting reflected sound signals into segments (range gating) so that velocities of water current can be determined at preset or user-predetermined intervals along the acoustic path (called cells or bins) [10]. ADCPs have the added advantage of providing data that reflect conditions across large areas in the water column [13].
Backscattered sound waves (acoustic backscatter) from ADV and ADCP velocity measurements have been used to estimate TSS in aquatic systems [10,14,18,22,23,24,25,26]. The fundamental theory is that acoustic waves moving through water containing sediments will scatter and attenuate depending on the characteristics of the fluid, sediment (concentration, shape size), and acoustic device [25]. Furthermore, adjustments (corrections) are usually introduced to account for the influence of the fluid and device characteristics [13,25]. Estimates of TSS are obtained by calibrating the acoustic backscatter with actual measured TSS from the deployment site using regressions [10,23]. This calibration can then be used to provide a time series of TSS estimates during the deployment of the acoustic sensor at the site.
Although high-resolution estimates of TSS now exist, there is still a need to develop similar estimates for sediment-associated constituents such as particulate phosphorus to monitor and prevent excessive nutrient loading into delicate ecosystems. Studies have found higher concentrations of TSS and phosphorus during high discharge events in rivers and canals, with sediment resuspension and desorption of sediment P highlighted as potential reasons for such high concentrations [27,28].
The human-made canals of South Florida (home of the Everglades) were constructed to meet the agricultural, urban, and environmental needs of the region [29,30]. Discharges in these canals are highly managed through inflow structures [30], and it is not uncommon to have large masses of water discharged into canals over relatively short periods. Therefore, there is a need to understand how rapid discharge releases can impact the export of sediments and phosphorus into the receiving oligotrophic Everglades wetlands. The objectives of this study were to (1) use acoustic backscatter to develop high-resolution TSS and TPP estimates in a major South Florida canal and (2) evaluate how canal discharges can impact the downstream export of TSS and TPP.

2. Materials and Methods

2.1. Instrumentation

2.1.1. Characteristics and Limitations of Acoustic Devices

Acoustic devices collect many types of data, including velocity, signal amplitude/signal-to-noise ratio (SNR), correlation coefficient, temperature, and pressure (Table 1). Velocity is typically the primary parameter of interest, while the signal amplitude and correlation are used to provide data quality information [31]. The main purpose of the signal amplitude is to determine if there are enough particles in the water, making it an excellent illustration of sediment dynamics and fluctuations [32]. This signal amplitude constitutes the acoustic backscatter parameter used as a surrogate for total suspended solids. The signal amplitude data are usually accessed from the device as either signal amplitude (units of counts) or the SNR. For SONTEK instruments, the SNR is obtained from the signal amplitude by subtracting the ambient, background electronic noise level and converting it to decibel units (dB) (multiplying by 0.43) [32]. To obtain accurate velocity measurements, SNR values must be at least 5 dB (ADV) and 3 dB (ADCP) (Table 1). An SNR value of 0 dB indicates that the water is too clear and that there is no distinction between the signal (acoustic energy) and the ambient noise level [31,32]. The ADV uses the correlation coefficient as a second data quality parameter, and it is expressed as a percentage. A correlation score of 100% means reliable low-noise velocity measurements, which occur at an SNR level of about 15 dB, while a 0% correlation reflects velocity measurements dominated by noise [31,32]. Interferences that can lead to low correlation values include high velocities, turbulence, and aerated water. A score of 70% or higher indicates accurate velocity data [32].

2.1.2. Acoustic Method

The conversion of the acoustic backscatter to TSS estimates (TSS) requires a series of steps that are summarized in this equation:
TSS = 10 [ A + B RB ]
where A—intercept, B—slope, and RB—relative acoustic backscatter.
The exponent of Equation (1) contains a term for the relative acoustic backscatter, RB, measured by the acoustic device, and the terms for the equation coefficients—intercept, A, and slope, B. These coefficients are determined by regressions of backscatter with artificially created TSS concentrations or gravimetrically (differences in weights) measured TSS from water samples [10]. This equation can also be expressed as:
Log   TSS = A + B RB
RB is composed of the measured acoustic backscatter and a correction for transmission losses [22] in units of counts or decibels (dB), as shown below:
  RB = RL + 2 TL
where RL is the reverberation level (the measured backscatter) and 2TL is the two-way transmission loss, assumed to be the same in each direction. The RL can be used directly as the actual reported values from the device (or the reported values minus the noise level) if the noise level is relatively constant [13]. This RL is multiplied by a manufacturer-supplied scale factor to convert it from count to decibels (dB). After the conversion, the RL can be referred to as the echo intensity (EI).
The two-way transmission loss, 2TL, is introduced to compensate for losses due to the spherical spreading of the acoustic beam from the acoustic device and losses due to absorption in the water, as shown below:
2 TL   = 20   log 10 [ R   Ψ ] + 2 α R
where R—Range, Ψ—near-field correction, and α—absorption correction.
The first part of the equation is the correction for beam spreading, and the second part is the absorption correction. R is the slant range from the transducer head (ADCP) to the measured bin (m) given as:
R = r + D 4
where r is the slant distance from the transducer head to the center of the bin in meters and D is the bin size in meters [24]. The next term, Ψ, is the transducer near-field correction [37] that accounts for the non-spherical spreading of acoustic energy close to the transducer and is defined as:
Ψ   = [ 1 + 1.35 z + ( 2.5 z ) 3.2 ] / [ 1.35 z + ( 2.5 z ) 3.2 ]
and
z = R   λ   π a t 2
where at is the transducer radius in meters and λ is the acoustic wavelength in meters. However, Ψ is empirically derived and generally not as reliable as the other parts of the transmission loss equation, and could, thus, introduce some level of uncertainty [25].
The next term in the 2TL equation is the α coefficient, which describes the absorption of energy by water (αw) and attenuation due to suspended sediments (αs)—that is, α = αw + αs (all in dB m−1). The coefficient αw describes the absorption of energy by water and depends on the salinity, temperature, and pressure of the water and the sound frequency of the acoustic device [13,25]. Pressure does not have a significant effect for shallow water environments (depth ≤ 20 m), and salinity is not applicable in freshwater systems [13,24], leaving temperature and frequency as the determining variables in canals, rivers, and lakes. Next, the coefficient αs, the sediment attenuation coefficient, is controlled by the sediment characteristics for a given acoustic frequency and is negligible in water with a TSS concentration less than 300 mg L−1 [38].
The correction for transmission losses is more significant for ADCPs, where the range, R, can extend many meters, while ADVs have much shorter ranges (e.g., 10 cm), and any losses are comparatively negligible [13].

2.2. Study Area

The Tamiami (L-29) Canal in the Everglades region in South Florida was created by means of the excavation of limestone used to construct the Tamiami Trail roadway (Figure 1). The roadway acted as a barrier to natural overland flow and restricted freshwater inflows to the Northeast Shark River Slough (NESRS) within the Everglades National Park (ENP). To remediate the limitations on freshwater inflows, the Modified Deliveries Project was created to develop the necessary infrastructure to allow more freshwater delivery into the ENP and restore a “more natural” hydrologic condition [39,40]. The infrastructure includes the construction of approximately 5.8 km of bridges on the Tamiami Trail to enable more flow into the NESRS through the S333 water control (inflow) structure (spillway) along L-29 [41]. A section of the L-29 Canal east of the spillway was used in this study.

2.3. Deployment in the Tamiami (L-29) Canal and Backscatter Processing

Two acoustic devices, a 10 MHz side-looking ADV and a 1500 kHz Argonaut ADCP, both manufactured by SonTek (SonTek—a Xylem Brand, San Diego, CA, USA), were deployed at three sites in the Tamiami (L-29) Canal, namely upstream (UP—25°45′40.8″ N, 80°40′18.8″ W), downstream (DS—25°45′40.6″ N, 80°39′06.3″ W), and interior (INT—25°45′39.1″ N, 80°32′03.5″ W), at different periods between June and December 2021. The ADV and ADCP were deployed during the same period at DS (June) and INT (August), while at UP, the ADV was deployed in November, and the ADCP was deployed in September. UP was closest to the inflow structure S333, while DS and INT were located before the 4.2 km and 1.6 km bridges, respectively. These sites were strategically chosen to determine how discharge releases into the L-29 Canal can transport suspended sediments and particulate phosphorus into the ENP via the bridges.

2.3.1. Acoustic Doppler Velocimeter (ADV)

The ADV was deployed near the canal bed (0.14 m–0.5 m above the bed) to measure the velocity and backscatter every 2 h. The ADV had three receivers and made measurements 10 cm away from the probe (sampling volume). The backscatter from the three receivers was averaged and converted from signal amplitude (counts) into echo intensity (decibels) by multiplying by a 0.43 scaling factor [32] before being used for calibration.

2.3.2. Acoustic Doppler Current Profiler (ADCP)

The three-beamed Argonaut ADCP was deployed at the canal bottom, in an upward facing direction, and was configured to make temporal and vertical velocity and backscatter measurements with a vertical resolution (bin size) that was either set at 0.5 m or 1 m (Figure 2). During deployments, the water depths in the canal ranged from 3.6 to 5.0 m and provided between 4 and 7 measurement bins. During longer deployments (1–3 weeks), the ADCP was programmed to make either 120- or 300-s measurements (averaging interval) every 1200 or 3000 s, respectively (sampling interval). For shorter deployments, the ADCP was programmed to make 10-s measurements every 100 s. Additionally, the blanking distance—the region directly above the ADCP where no measurements can be made—ranged from 0.5 m to 1 m.

2.4. Acoustic Backscatter Conversion and Correction

For the ADV, laboratory experiments were conducted to develop calibration curves between TSS and echo intensity (EI). Sediment cores from the surface of the canal bed were first taken from each site. Then, portions of the cores were weighed (wet weight), dried to a constant weight and reweighed to determine the water content of the sediments. Next, more fresh wet surface sediments were weighed, and the water content was accounted for so that the true (dry) weights of the sediments could be calculated. These true weights were then diluted in measured volumes of tap water in a bucket to artificially create TSS between 5 mg L−1 and 150 mg L−1 in a large bucket. Next, a pump was inserted into the bucket to create velocity and keep the sediments in suspension [26]. During agitation, the ADV was introduced into the bucket and programmed to make 15-s velocity and backscatter measurements for approximately 2 min. The agitation and measurements were repeated for all the created TSS (between 5 mg L−1 and 150 mg L−1), and the corresponding acoustic backscatter from the ADV was recorded, converted to echo intensity (EI) by multiplying by 0.43, and used to generate a calibration curve for each site.
To develop time series profiles of suspended sediments from ADCP backscatter, a series of steps were followed. First, the backscatter was downloaded from the ADCP and averaged (the three beams of the ADCP produced three backscatter values per measurement). In cases when the backscatter from one beam was much higher than the other, or when the difference in SNR values was greater than 10 dB [42], the averages of the other two beams were used. Next, the backscatter was converted to EI in dB at each bin. Since the noise level was relatively constant during deployment, it was not subtracted from the backscatter before conversion to EI. Then, the EI was corrected for transmission loss from beam spreading and attenuation due to water to give the final RB (Equation (3)). The αw values were based on the instrument frequency and water temperature since the study site was a freshwater canal with a water depth much less than 20 m. Published αw values provided by [25] were used. Attenuation due to sediment was ignored because of the low TSS measured in the water column. Data and graphical analyses were conducted using R [43], Microsoft Excel and PowerPoint, and MATLAB [44]. An alpha (p) value to define statistical significance was set at 0.05 for calibration regressions.

2.5. Water Sampling for TSS and TPP for Calibrations and Estimations

To provide TSS measurements for the calibration of the ADCP backscatter, water samples were collected simultaneously during the ADCP deployments (at UP, a calibration performed in November was applied to the September deployment). To perform the calibrations, water samples were collected from different depths in the water column that matched the bins of the ADCP (Figure 2). A peristaltic pump was used to collect the water samples in prerinsed 1 L sampling bottles, which were kept on ice and transported to the lab where they were analyzed gravimetrically for TSS. In the laboratory, predried (105 °C for 1 h) and preweighed 47 mm GFF Whatman filters were used to filter 1 L (volume) of the collected water samples (Wpre). After filtration, the filter and particles (residue) were dried again (105 °C for 1 h), cooled in a desiccator, and weighed (Wpost) [14,45]. TSS was then calculated according to the following equation:
TSS   ( mgL 1 ) = Wpost ( g ) Wpre   ( g ) 1000   V ( L )
where W—weight and V—volume.
During the November calibration, filtering was performed in the field using predried and preweighed filters in a container. After filtering, the filters were kept, sealed, and labeled in 47 mm Petri dishes.
After analysis for TSS, the filters together with the residues were digested and analyzed to obtain the concentration of total phosphorus on the particles—total particulate phosphorus (TPP)—using the ascorbic acid method [46]. This TPP (termed TPP in TSS, μg mg−1) was multiplied by the TSS estimates (mg L−1) that were derived from applying the calibration curve to obtain the total particulate phosphorus (TPP in μg L−1) at each depth in the water column for every deployment. The equation is shown below:
TPP   estimate in   water ( μ g   L 1 ) = TSS   estimate   ( mg   L 1 ) TPP   in   TSS   ( μ g   mg 1 )

2.6. Discharge Data

Instantaneous discharges at UP (25°45′40.8″ N, 80°40′18.8″ W) were downloaded from the S333 inflow structure available at South Florida Water Management District’s online data repository—DBHYDRO (https://www.sfwmd.gov/science-data/dbhydro, accessed on 12 November 2022). Instantaneous discharges at DS (25°45′41″ N, 80°39′17″ W) and INT (25°45′41.0″ N, 80°32′12.0″ W) were downloaded from the United States Geological Survey-National Water Information System (https://waterdata.usgs.gov/nwis, accessed on 26 October 2022). The discharges at INT had both positive and negative values. The negative discharges were taken as flows in the reverse direction (east to west) caused by factors including gate operations at a pump at the eastern end of the canal and wind [47]. Such negative values were converted to positive values for analyses.

3. Results

3.1. Acoustic Backscatter Processing for the Acoustic Doppler Current Profiler (ADCP)

After converting the RL in counts to EI in decibels (dB) and correcting for beam spreading and water absorption, the relative acoustic backscatter (RB, Equation (3)) from the ADCP was obtained. The RB was not clearly distinguishable between depths in the water column (Figure 3a,b). Using the RB values to develop calibration curves resulted in a very narrow range of backscatter values, for instance, 75.5–77.5 dB at UP (not shown). Applying the equation coefficients (slope and intercept) of the calibration curves to obtain the TSS estimates resulted in values that ranged from 0 to over 1000 mg L−1 at UP, which was unrealistic in this relatively low-particulate-level canal (Table 2). Similarly, converting the estimated TSS into TPP resulted in higher concentrations in the surface layer of the water compared to the deeper layer, which conflicted with the measured results that were used for the calibrations. Plots using only the EI yielded better results, with the values being more distinguishable between depths (Figure 3c,d).

3.2. Calibration Curves

There were statistically significant positive correlations between echo intensity and log TSS for ADV and ADCP across all three sites, as shown in the calibration curves (Figure 4a,b). Additionally, the uncorrected EI produced better regression for calibration because there was more spread in the EI: 60–70 dB at UP, 50–70 dB at DS, and 55–65 dB at INT (Figure 4b).

3.3. Measured Total Suspended Solids (TSS) and Total Particulate Phosphorus (TPP) Used in Calibrations

The gravimetrically analyzed TSS (mg L−1) and the resulting calculated TPP in the water column (μg L−1) revealed that the top layers of the water column generally had higher concentrations than the deeper layers across the three sites (Table 2). Conversely, the phosphorus fraction of the analyzed sediments (TPP in TSS, μg mg−1) generally decreased from the top to the deeper layers. Although the sample collections and analyses were performed at different times of the year across sites, the TPP in TSS was lowest at UP, which was closest to the inflow structure, and highest at INT, which was the most downstream site (Table 2).
During the deployment of the acoustic devices, the discharges at UP had values up to 60 m3 s−1 (ADV, November) and 40 m3 s−1 (ADCP, September) (Figure 5a and Figure 6b). During the combined ADV and ADCP deployments at DS in June (DS-Jun), the maximum discharge was approximately 22 m3 s−1 (Figure 5b), while the discharges during the second ADCP deployment in December (DS-Dec) reached 40 m3 s−1 (Figure 6f). The DS-Dec deployment had higher discharges because the main inflow structure, S333, and a newer secondary structure, S333N (located north of S333), were both operational (released discharge into the canal) at that time, while during the DS-Jun deployment, only S333N was operational. At INT, the maximum discharge during the combined deployments did not exceed 14 m3 s−1 (Figure 5c).
In terms of water quality estimates, at UP, the TSS and TPP estimates (concentrations) did not exceed 18 mg L−1 and 17 μg L−1, respectively, near the canal bed (ADV, Figure 5a) and 15 mg L−1 and 13 μg L−1, respectively, in the water column (ADCP, Figure 6a and Figure 7a). At DS-Jun, concentrations were less than 14 mg L−1 for TSS and 16 μg L−1 for TPP near the canal bed (Figure 5b), while in the water column, the concentrations were less than 9 mg L−1 for TSS and 17 μg L−1 for TPP (Figure 6c and Figure 7b). At DS-Dec, the TSS concentration in the water column peaked at 11 mg L−1 (Figure 6e), and the higher concentrations compared to DS-Jun can be attributed to the higher discharges in December. No estimation for TPP in Dec-DS was performed because the number and sizes of measurement depths (ADCP bins) were different between the two deployments. DS-Jun had four 1-m bins, while DS-Dec had five 0.5-m bins, so they were not compatible beyond the calibration curve. Lastly, at INT, the near-bed concentrations did not exceed 5 mg L−1 for TSS and 12 μg L−1 for TPP (Figure 5c), while in the water column, the concentrations did not exceed 2.5 mg L−1 for TSS and 8 μg L−1 for TPP (Figure 6g and Figure 7c). These results generally showed higher TSS and TPP concentrations near the canal bed than in the water column.
An observation of the relationship between discharge and the estimates showed clear associations between higher discharges and higher concentrations of TSS and TPP at UP (Figure 5a and Figure 6a,b). The instantaneous and subsequent sustained discharge releases from the S333 structure caused a noticeable increase in TSS and TPP concentrations. Similarly, higher discharges matched higher TSS and TPP concentrations, while lower discharges matched lower TSS and TPP concentrations near the canal bed at DS-Jun (Figure 5b) and for TSS in the water column at DS-Dec (Figure 6e,f). This synchronization was not as strong in the water column at DS-Jun, probably due to fewer particles in the water column compared to the bed (Table 2), while at INT, there was no discernable relationship between discharge and TSS/TPP concentrations due to the relatively lower values of all three parameters.

4. Discussion

4.1. Justification for Not Correcting for Transmission Losses for ADCP Backscatter Processing

Studies that used ADCPs for TSS estimations and corrected for transmission losses had a wide range of concentrations, relatively high TSS concentrations, and/or greater water depths when compared to this study. For example, in a study performed in the Lembeh Strait, in Indonesia, the measured water depth was between 15 and 30 m, and the measured TSS ranged from 55 to 74 mg L−1 [14]. Similarly, in another study performed in the Tidung Island seawaters of Indonesia, the water depth was greater than 40 m, and the TSS ranged from <45 to 80 mg L−1 [48]. A study performed in the Hudson River in New York, USA, was conducted in a water depth of 18 m that had TSS concentrations that ranged from <10 and 100 mg L−1 [24,49]. Another study performed in San Francisco Bay, USA, had water depths between 7.3 and 16.1 m, and TSS concentrations that were greater than 100 mg L−1 [10]. Comparatively, the concentrations in our study were quite low (generally <10 mg L−1) and did not vary much between depths, and the water depth in the canal was approximately 5 m. Thus, the poor suitability of the relative backscatter (RB) for developing a useable calibration curve for TSS estimations in this study could be due to the relatively narrow range of acoustic backscatter and TSS in the water column in the L-29 Canal compared to other studies.
The correction introduced for transmission losses increases with increasing distance from the ADCP because the range (distance) from the device is considered. This means that more corrections are applied with increasing distance from the ADCP. Even though, in our study, the bin closest to the ADCP (Bin 1) had the highest RB values that corresponded with the generally higher TSS concentrations at the deeper water column (Table 2), the difference in RB between this bin and the farthest bin was not enough; when the corrections were introduced, the RB values became very similar among all the bins. This resulted in the lack of clear distinction between the RB values with depth, as shown in the results, while the backscatter, uncorrected for transmission losses, produced clearer distinctions with depth. Therefore, the application of such correction “overcorrected” the farthest bins. The use of the Downing near-field factor (Ψ) can overcorrect and bias the relative backscatter, and thus, it may be preferable to use uncorrected data [25], although in some cases, this overcorrection may be marginal [50]. Therefore, it was prudent to remove the two-way correction losses from the relative backscatter processing.

4.2. Vertical Profiles of Total Suspended Solids and Total Particulate Phosphorus in the Water Column of the L-29 Canal

The low values of TSS and TPP estimates from this study are not uncommon in the region of Central and South Florida [51,52]. For example, studies conducted in Lake Worth Lagoon (South Florida) and in the tributaries of the Indian River Lagoon (East Central Florida) generally had TSS concentration averages less than 12 mg L−1 [52,53] and an average TPP concentration in the water column of less than 70 μg L−1 [52]. Similarly, in a study performed along the L-40 Canal, north of the L-29 Canal (our study), the average TSS concentration did not exceed 18 mg L−1, and the average TPP concentration was 71 μg L−1 [51]. Such low TSS concentrations have been found in <5% of rivers globally [52]. However, the average TPP in TSS from our L-29 study, which was 2.05 μg mg−1, is higher than the range of global averages for rivers (1–2.01 μg mg−1) [54,55]. Therefore, despite the relatively low TSS in the L-29 Canal, the relatively higher TPP in TSS can threaten the ecological integrity of downstream waters.
The general increase in TSS and TPP estimates from the water surface to the canal bed was expected. We surmise that this is due to resuspension from the bed that caused higher concentrations of TSS [14,56] and the tendency for more suspended sediments to settle out in deeper waters [57]. Therefore, higher suspended sediment levels in deeper waters will result in higher TPP (μg L−1) in deeper waters as well. The opposite (decreasing) trend with depth observed for the TPP in TSS (μg mg−1) has been reported in other studies (e.g., [56]). A possible explanation for this trend is that the composition of suspended sediments varied from the water surface to the canal bed. The higher TPP in TSS suggests that the water surface had light, organic, and P-rich particles, while the deeper water had denser, more mineral, and P-poor particles. Light and flocculent P-rich organic material (aquatic plant detritus) has been found suspended on the water surface of farm canals [58,59].

4.3. Impact of Canal Discharges and Hydrology on Total Suspended Solids and Total Particulate Phosphorus

In southern Florida, canal discharge releases via inflow structures have been shown to increase the concentration and export of suspended sediments and particulate phosphorus in two steps. The first step, referred to as the first flush, marks the start of canal discharge release from the inflow structure after a greater than 24-h quiescent period. This first flush mobilizes and entrains materials deposited from the previous event as well as new material that may have accumulated or grown during the quiescent period. After the first flush is the second step, termed the cumulative flush, where sustained high discharges can mobilize more sediments and their associated P contents [59]. In this study, the effects of first flush and cumulative flush were also observed at the upstream site, closest to the S333 inflow structure, where the increases in TSS and TPP closely followed the discharges from the structure. Instantaneous discharges have been compared to storm events in their capacity to mobilize sediments and increase their concentrations and export [28]. In rivers, the first flush occurs at the start of a rain or storm event, where higher TSS can be recorded during the onset of very high flows that decrease over the duration of the event [52].
Even though the deployments and water quality estimates at the sites were performed at different times during the year, inferences can be made as to the origins and types of suspended sediments in the three canal sites by considering their proximities to the S333 inflow structure. The lower TPP in TSS at UP possibly occurred because high discharges from S333 exported all the fresh, low-density, and P-rich organic particles that accumulated or grew during the quiescent period and left behind the older and denser P-poor particles. Studies have shown that organic particles have a high phosphorus content and are easily transportable, while older particles are more mineral, P-poor, and heavier, which makes them harder to transport [59]. The higher TPP concentrations in the water column at DS suggest that discharges from the secondary inflow structure, S333N, as well as earlier discharges from the main S333 structure, transported P-rich organic particles from upstream which accumulated at the downstream site. Another reason is that the transporting power of the canal discharges was attenuated with increasing distance from the S333 inflow structure; therefore, only the lightest particles made it the farthest downstream. This decrease in transporting power downstream also implies the dominant effect of biogeochemical processes on the kind of particles present and the P availability of such particles. Biogeochemical processes such as adsorption–desorption, biologic uptake, and the decomposition of organic matter have been shown to dominate P availability during low-flow periods [60,61], and in this case, the downstream reaches of the canal can display these low flow conditions. Similarly, the highest TPP in TSS found in the canal interior, at the most downstream site from the S333 inflow structure, indicates that the particles therein may be of biological origin. A high phosphorus mass content has been associated with fresh biological organic matter [59]. These inferences are noteworthy because DS and INT are located close to the bridges that open into the Everglades National Park, and the potential export of P-rich particles into the park will be detrimental to Everglades restoration efforts.

5. Conclusions

This study used the acoustic backscatter (echo intensity), a biproduct of velocity measurement from two acoustic devices, an acoustic Doppler velocimeter (ADV) and an acoustic Doppler current profiler (ADCP), to estimate total suspended solids (TSS) and total particulate phosphorus (TPP) in the L-29 Canal in South Florida, USA. The echo intensity from the ADCP, uncorrected for beam spreading and attenuation due to water and sediments, proved sufficient to develop suitable calibration curves to estimate realistic TSS in canal locations that varied by distance from a managed inflow structure, S333. The suspended sediments used in the calibrations were then analyzed for total particulate phosphorus concentrations to obtain the TPP estimates in the water column (ADCP) and near the canal bed (ADV). Time series plots of discharges, TPP and TSS revealed that the site closest to S333, the upstream site, had the highest concentrations, driven by instantaneous and rapid discharge releases. Analyses of the TSS particles further revealed differences between the upstream and downstream sites. Suspended sediments further downstream had a higher phosphorus (P) content per suspended sediments (TPP in TSS μg mg−1) because they likely had higher organic matter (biological origin). Conversely, the particles upstream had a lower P content because of their proximity to the high discharges from S333, which transported the lighter P-rich organic particles downstream and left behind denser, P-poor mineral particles. Such important observations and inferences in this study were made possible using high-resolution data from acoustic devices. Therefore, the ability to evaluate how rapid changes in hydrology affect the potential export of sediments and phosphorus is critical for canal discharge management to enable the preservation of the integrity of downstream aquatic systems, such as the Florida Everglades.
This work is novel, because it extends the use of backscatter from acoustic devices, which has been traditionally used as a surrogate for TSS, to also using it as a surrogate for TPP. This is crucial for effective water quality and nutrient (phosphorus) monitoring because it reduces the human and material risks and costs of routine manual samplings and provides high spatiotemporal data. Additionally, this work practically tests the lower limits of using acoustic backscatter as a surrogate for TSS, because the concentrations were at the low end of detection (generally less than 10 mg L−1 of TSS), while other studies have focused on much higher concentrations. A general limitation of using acoustic devices as surrogates for TSS and TPP in low-particulate-level water is the possibility of weak signal amplitudes and low echo intensities, which may be indistinguishable from background noise. Another limitation is the potential difficulty in creating low TSS in the laboratory for the calibration of the low echo intensities that may occur under certain field conditions and in certain locations. In many locations, and specifically in the highly managed systems of South Florida, where flow conditions rapidly respond to pumping operations and other water management manipulations, short-term response measurements are needed, as they provide information not obtained by longer-term, less temporally robust data. Future studies should be conducted in several additional canal systems to better determine the range of applicability of these methods.

Author Contributions

Conceptualization, I.S.O., L.J.S. and D.C.F.; methodology, I.S.O., L.J.S. and D.C.F.; software, I.S.O.; validation, I.S.O., L.J.S. and D.C.F.; formal analysis, I.S.O.; investigation, I.S.O. and L.J.S.; resources, I.S.O. and L.J.S.; data curation, I.S.O.; writing—original draft preparation, I.S.O.; writing—review and editing, I.S.O., L.J.S. and D.C.F.; visualization, I.S.O. and D.C.F.; supervision, L.J.S. and D.C.F.; project administration, I.S.O. and L.J.S.; funding acquisition, I.S.O. and L.J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research is part of the first author’s dissertation research, which is funded by an FIU For Everglades Fellowship provided by the Everglades Foundation, NSF CREST Center for Aquatic Chemistry and Environment, Cristina Menendez Fellowship provided by the Miccosukee Tribe of Indians Endowment through the Institute of Environment, Graduate Student Research Support Program and Dissertation Year Fellowship provided by the University Graduate School, and the Department of Earth and Environment, all at Florida International University, and General Endowment Award from the Society for Freshwater Science.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Publicly available datasets were analyzed in this study. These data can be found here: https://www.sfwmd.gov/science-data/dbhydro and https://waterdata.usgs.gov/nwis (accessed on 1 November 2022).

Acknowledgments

The authors appreciate the assistance of the South Florida Water Management District (SFWMD) and the United States Geological Survey (USGS) for the published and publicly accessible data used in this study. The authors thank Assefa Melesse for reviewing this manuscript and for constructive comments. The authors also thank Luke Lamb-Wotton, Joshua Linenfelser, Jordan Massie, and Jonathan Rodemann for their help with fieldwork. This material is based upon work supported by the National Science Foundation under Grant No. HRD-1547798 and HRD-2111661. These NSF Grants were awarded to Florida International University as part of the Centers of Research Excellence in Science and Technology (CREST) Program. This is contribution #1532 from the Institute of Environment at Florida International University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wood, P.J. Biological Effects of Fine Sediment in the Lotic Environment. Environ. Manag. 1997, 21, 203–217. [Google Scholar] [CrossRef]
  2. Kronvang, B.; Laubel, A.; Grant, R. Suspended Sediment and Particulate Phosphorus Transport and Delivery Pathways in an Arable Catchment, Gelbæk Stream, Denmark. Hydrol. Process. 1997, 11, 627–642. [Google Scholar] [CrossRef]
  3. Gartner, J.W. Estimation Of Suspended Solids Concentrations Based On Acoustic Backscatter Intensity: Theoretical Background. In Proceedings of the Turbidity and Other Sediment Surrogates Workshop, Reno, NV, USA, 30 April–2 May 2002. [Google Scholar]
  4. Latosinski, F.G.; Szupiany, R.N.; García, C.M.; Guerrero, M.; Amsler, M.L. Estimation of Concentration and Load of Suspended Bed Sediment in a Large River by Means of Acoustic Doppler Technology. J. Hydraul. Eng. 2014, 140, 04014023. [Google Scholar] [CrossRef]
  5. Noe, G.B.; Harvey, J.W.; Saiers, J.E. Characterization of suspended particles in Everglades wetlands. Limnol. Oceanogr. 2007, 52, 1166–1178. [Google Scholar] [CrossRef]
  6. Noe, G.B.; Harvey, J.W.; Schaffranek, R.W.; Larsen, L.G. Controls of Suspended Sediment Concentration, Nutrient Content, and Transport in a Subtropical Wetland. Wetlands 2010, 30, 39–54. [Google Scholar] [CrossRef]
  7. Zhao, H.; Piccone, T.; Diaz, O. Supporting Information for Canal Evaluations; Technical Publication WR-2015-003; South Florida Water Management District: West Palm Beach, FL, USA, 2015. [Google Scholar]
  8. Dorioz, J.M.; Cassell, E.A.; Orand, A.; Eisenman, K.G. Phosphorus storage, transport and export dynamics in the Foron River watershed. Hydrol. Process. 1998, 12, 285–309. [Google Scholar] [CrossRef]
  9. Jordan-Meille, L.; Dorioz, J.M. Soluble phosphorus dynamics in an agricultural watershed. Agronomie 2004, 24, 237–248. [Google Scholar] [CrossRef] [Green Version]
  10. Gartner, J.W. Estimating suspended solids concentrations from backscatter intensity measured by acoustic Doppler current profiler in San Francisco Bay, California. Mar. Geol. 2004, 211, 169–187. [Google Scholar] [CrossRef]
  11. Rai, A.K.; Kumar, A. Continuous measurement of suspended sediment concentration: Technological advancement and future outlook. Measurement 2015, 76, 209–227. [Google Scholar] [CrossRef]
  12. Harrison, J.W.; Lucius, M.A.; Farrell, J.L.; Eichler, L.W.; Relyea, R.A. Prediction of stream nitrogen and phosphorus concentrations from high-frequency sensors using Random Forests Regression. Sci. Total. Environ. 2021, 763, 143005. [Google Scholar] [CrossRef]
  13. Öztürk, M. Sediment Size Effects in Acoustic Doppler Velocimeter-Derived Estimates of Suspended Sediment Concentration. Water 2017, 9, 529. [Google Scholar] [CrossRef] [Green Version]
  14. Dwinovantyo, A.; Manik, H.M.; Prartono, T.; Susilohadi, S. Quantification and Analysis of Suspended Sediments Concentration Using Mobile and Static Acoustic Doppler Current Profiler Instruments. Adv. Acoust. Vib. 2017, 2017, 1–14. [Google Scholar] [CrossRef] [Green Version]
  15. Sedimentology of Aqueous Systems; Cristiano, P.; Charlesworth, S.M. (Eds.) John Wiley & Sons: Hoboken, NJ, USA, 2010. [Google Scholar]
  16. Gray, J.R.; Gartner, J.W. Technological advances in suspended-sediment surrogate monitoring. Water Resour. Res. 2009, 45. Available online: https://onlinelibrary.wiley.com/doi/abs/10.1029/2008WR007063 (accessed on 12 December 2022). [CrossRef] [Green Version]
  17. Dominguez Ruben, L.G.; Szupiany, R.N.; Latosinski, F.G.; López Weibel, C.; Wood, M.; Boldt, J. Acoustic Sediment Estimation Toolbox (ASET): A software package for calibrating and processing TRDI ADCP data to compute suspended-sediment transport in sandy rivers. Comput. Geosci. 2020, 140, 104499. [Google Scholar] [CrossRef]
  18. David, C.F.; Friedrichs, C.T. Versatility of the Sontek ADV: Measurements of sediment fall velocity, sediment concentration, and TKE production from wave contaminated velocity data. In Proceedings of the Coastal Sediments, Clearwater Beach, FL, USA, 18–23 May 2003. [Google Scholar]
  19. Schoellhamer, D.H.; Wright, S.A. Continuous measurement of suspended-sediment. In Erosion and Sediment Transport Measurement in Rivers: Technological and Methodological Advances; IAHS: Wallingford, UK, 2003; p. 28. [Google Scholar]
  20. Hamilton, L.J.; Shi, Z.; Zhang, S.Y. Acoustic Backscatter Measurements of Estuarine Suspended Cohesive Sediment Concentration Profiles. J. Coast. Res. 1998, 14, 1213–1224. [Google Scholar]
  21. Kostaschuk, R.; Best, J.; Villard, P.; Peakall, J.; Franklin, M. Measuring flow velocity and sediment transport with an acoustic Doppler current profiler. Geomorphology 2005, 68, 25–37. [Google Scholar] [CrossRef]
  22. Thevenot, M.M.; Prickett, T.L.; Kraus, N.C. Tylers Beach, Virginia, Dredged Material Plume Monitoring Project, September 27 to October 4, 1991; Coastal Engineering Research Center: Vicksburg, MS, USA, 1992. [Google Scholar]
  23. Fugate, D.C.; Friedrichs, C.T. Determining concentration and fall velocity of estuarine particle populations using ADV, OBS and LISST. Cont. Shelf. Res. 2002, 22, 1867–1886. [Google Scholar] [CrossRef]
  24. Wall, G.R.; Nystrom, E.A.; Litten, S. Use of an ADCP to Compute Suspended-Sediment Discharge in the Tidal Hudson River, New York; No. 2006-5055; USGS: Troy, NY, USA, 2006. [Google Scholar]
  25. Landers, M.N.; Straub, T.D.; Wood, M.S.; Domanski, M.M. Sediment Acoustic Index Method for Computing Continuous Suspended-Sediment Concentrations (Internet); Techniques and Methods; vols. 3-C5; Report No.: 3-C5; U.S. Geological Survey: Reston, VA, USA, 2016; p. 83. Available online: http://pubs.er.usgs.gov/publication/tm3C5 (accessed on 13 June 2021).
  26. Chanson, H.; Takeuchi, M.; Trevethan, M. Using turbidity and acoustic backscatter intensity as surrogate measures of suspended sediment concentration in a small subtropical estuary. J. Environ. Manag. 2008, 88, 1406–1416. [Google Scholar] [CrossRef] [Green Version]
  27. Diamond, J.S.; Cohen, M.J. Complex patterns of catchment solute–discharge relationships for coastal plain rivers. Hydrol. Process. 2018, 32, 388–401. [Google Scholar] [CrossRef]
  28. Onwuka, I.S.; Scinto, L.J.; Mahdavi Mazdeh, A. Comparative Use of Hydrologic Indicators to Determine the Effects of Flow Regimes on Water Quality in Three Channels across Southern Florida, USA. Water 2021, 13, 2184. [Google Scholar] [CrossRef]
  29. Harvey, R.G.; Loftus, W.F.; Rehage, J.S.; Mazzotti, F.J. Effects of Canals and Levees in Everglades Ecosystems; EDIS: Gainesville, FL, USA, 2011. [Google Scholar]
  30. Carter, K.; Redfield, G.; Ansar, M.; Glenn, L.; Huebner, R.; Maxted, J.; Pettit, C.; VanArman, J. Canals in South Florida: A Technical Support Document; SFWMD: West Palm Beach, FL, USA, 2010. [Google Scholar]
  31. Velasco, D.W.; Huhta, C.A. Experimental Verification of Acoustic Doppler Velocimeter (ADV®) Performance in Fine-Grained, High Sediment Concentration Fluids; Sontek: San Diego, CA, USA, 2009. [Google Scholar]
  32. SonTek. ADVField Technical Manual; SonTek Technical Documentation; Sontek: San Diego, CA, USA, 2001. [Google Scholar]
  33. Current Profiling Solutions Brochure. Available online: https://www.ysi.com/File%20Library/Documents/Brochures%20and%20Catalogs/current-profiling-solutions.pdf. (accessed on 8 February 2023).
  34. Acoustic Doppler Velocimeters (ADVs®)—Xylem. Available online: https://www.xylem.com/siteassets/brand/sontek/resources/brochure/sontek-acoustic-doppler-velocimeters-2017.pdf (accessed on 8 February 2023).
  35. SonTek. Argonaut™-SL System Manual Firmware; Version 12.0; SonTek™/YSI Corporation: San Diego, CA, USA, 2009; p. 316. [Google Scholar]
  36. Argonaut-XR MULTI-CELL DOPPLER CURRENT PROFILER. Available online: https://www.ysi.com/File%20Library/Documents/Brochures%20and%20Catalogs/argonaut-xr-brochure.pdf (accessed on 8 February 2023).
  37. Downing, A.; Thorne, P.D.; Vincent, C.E. Backscattering from a suspension in the near field of a piston transducer. J. Acoust. Soc. Am. 1995, 97, 1614–1620. [Google Scholar] [CrossRef]
  38. Gartner, J.W.; Cheng, R.T. The promises and pitfalls of estimating total suspended solids based on backscatter intensity from acoustic Doppler current profilers. In Proceedings of the Seventh Federal Interagency Sedimentation Conference, Reno, NV, USA, 25–29 March 2001. [Google Scholar]
  39. McLean, A. Modified Water Deliveries: Improving Hydrologic Condition in Northeast Shark River Slough. 2015. Available online: https://www.nps.gov/ever/learn/nature/modwater.htm (accessed on 10 May 2020).
  40. Jacksonville District > Missions > Environmental > Ecosystem Restoration > Combined Operational Plan (Internet). Available online: https://www.saj.usace.army.mil/Missions/Environmental/Ecosystem-Restoration/G-3273-and-S-356-Pump-Station-Field-Test/ (accessed on 7 January 2023).
  41. Sarker, S.K.; Kominoski, J.S.; Gaiser, E.E.; Scinto, L.J.; Rudnick, D.T. Quantifying effects of increased hydroperiod on wetland nutrient concentrations during early phases of freshwater restoration of the Florida Everglades. Restor. Ecol. 2020, 28, 1561–1573. [Google Scholar] [CrossRef]
  42. Jenner, B. Beam Separation in ADPs and ADCPs. Available online: https://www.ysi.com/File%20Library/Documents/Technical%20Notes/beam-separation-in-adps-and-adcps.pdf (accessed on 5 July 2021).
  43. R: The R Project for Statistical Computing (Internet). Available online: https://www.r-project.org/ (accessed on 7 January 2023).
  44. MathWorks—Makers of MATLAB and Simulink (Internet). Available online: https://www.mathworks.com/ (accessed on 7 January 2023).
  45. Standard Operating Procedure for: Total Suspended Solids. 2007. Available online: https://beta-static.fishersci.com/content/dam/fishersci/en_US/documents/programs/scientific/technical-documents/white-papers/apha-total-suspended-solids-procedure-white-paper.pdf (accessed on 2 March 2021).
  46. United States Environmental Protection Agency. Methods for Chemical Analysis of Water and Wastes; Environmental Monitoring and Support Laboratory: Cincinnati, OH, USA, 1983; pp. 475–490. [Google Scholar]
  47. Christian D Lopez Roque (Person)—ScienceBase-Directory (Internet). Available online: https://www.sciencebase.gov/directory/person/6250 (accessed on 6 April 2022).
  48. Manik, H.M.; Firdaus, R. Quantifying Suspended Sediment using Acoustic Doppler Current Profiler in Tidung Island Seawaters. Pertanika. J. Sci. Technol. 2021, 29, 363–385. [Google Scholar] [CrossRef]
  49. Wall, G.R.; Nystrom, E.A.; Litten, S. Suspended sediment transport in the freshwater reach of the Hudson River estuary in eastern New York. Estuaries Coasts 2008, 31, 542–553. [Google Scholar] [CrossRef] [Green Version]
  50. Venditti, J.G.; Church, M.; Attard, M.E.; Haught, D. Use of ADCPs for suspended sediment transport monitoring: An empirical approach: ADCPS FOR SUSPENDED SEDIMENT MONITORING. Water Resour. Res. 2016, 52, 2715–2736. [Google Scholar] [CrossRef]
  51. Daroub, S.H.; Josan, M.S.; Lang, T.A.; Waldon, M.G. L40 Canal Water Quality Survey Study Progress Report for 2007; University of Florida: Gainesville, FL, USA, 2007. [Google Scholar]
  52. Trefry, J.H.; Fox, A.L. Extreme Runoff of Chemical Species of Nitrogen and Phosphorus Threatens a Florida Barrier Island Lagoon. Front. Mar. Sci. 2021, 8, 752945. [Google Scholar] [CrossRef]
  53. Trefry, J.H.; Trocine, R.P.; Bennett, H. Sediment sourcing study of Lake Worth Lagoon and C-51 basin, Palm Beach County. In Final Report to Palm Beach County and the South Florida Water Management District for Contract R2008-0985; South Florida Water Management District: West Palm Beach, FL, USA, 2009. [Google Scholar]
  54. Savenko, V.S. Chemical composition of sediment load carried by rivers. Geochem. Int. 2007, 45, 816–824. [Google Scholar] [CrossRef]
  55. Viers, J.; Dupré, B.; Gaillardet, J. Chemical composition of suspended sediments in World Rivers: New insights from a new database. Sci. Total. Environ. 2009, 407, 853–868. [Google Scholar] [CrossRef]
  56. Li, J.; Reardon, P.; McKinley, J.P.; Joshi, S.R.; Bai, Y.; Bear, K.; Jaisi, D.P. Water column particulate matter: A key contributor to phosphorus regeneration in a coastal eutrophic environment, the Chesapeake Bay. J. Geophys. Res. Biogeosciences. 2017, 122, 737–752. [Google Scholar] [CrossRef]
  57. Wang, Q.; Li, Y.; Ouyang, Y. Phosphorus fractionation and distribution in sediments from wetlands and canals of a water conservation area in the Florida Everglades. Water Resour. Res. 2011, 47. [Google Scholar] [CrossRef] [Green Version]
  58. Daroub, S.H.; Stuck, J.D.; Lang, T.A.; Diaz, O.A. Particulate Phosphorus in the Everglades Agricultural Area: I—Introduction and Sources; EDIS: Gainesville, FL, USA, 2002. [Google Scholar]
  59. Daroub, S.H.; Stuck, J.D.; Lang, T.A.; Diaz, O.A. Particulate Phosphorus in the Everglades Agricultural Area: II—Transport Mechanisms; EDIS: Gainesville, FL, USA, 2002. [Google Scholar]
  60. Withers, P.J.A.; Jarvie, H.P. Delivery and cycling of phosphorus in rivers: A review. Sci. Total Environ. 2008, 400, 379–395. [Google Scholar] [CrossRef] [PubMed]
  61. Moatar, F.; Abbott, B.W.; Minaudo, C.; Curie, F.; Pinay, G. Elemental properties, hydrology, and biology interact to shape concentration-discharge curves for carbon, nutrients, sediment, and major ions. Water Resour. Res. 2017, 53, 1270–1287. [Google Scholar] [CrossRef]
Figure 1. (a) Map of the Everglades region showing the Tamiami (L-29) Canal (in red) and (b) study sites in the Tamiami (L-29) Canal, where the acoustic devices were deployed. The sites were upstream (UP), downstream (DS), and interior (INT). Modified from Google Earth (2023).
Figure 1. (a) Map of the Everglades region showing the Tamiami (L-29) Canal (in red) and (b) study sites in the Tamiami (L-29) Canal, where the acoustic devices were deployed. The sites were upstream (UP), downstream (DS), and interior (INT). Modified from Google Earth (2023).
Sensors 23 02281 g001
Figure 2. An illustration of a typical deployment of the acoustic devices: acoustic Doppler velocimeter (ADV) and acoustic Doppler current profiler (ADCP). Across deployment events, the user determined the blanking distance and bin sizes of the ADCP ranging from 0.5 to 1 m, and the number of bins varied from 4 to 7. The height above the bed (hab) for the ADV deployments ranged from 0.14 to 0.5 m.
Figure 2. An illustration of a typical deployment of the acoustic devices: acoustic Doppler velocimeter (ADV) and acoustic Doppler current profiler (ADCP). Across deployment events, the user determined the blanking distance and bin sizes of the ADCP ranging from 0.5 to 1 m, and the number of bins varied from 4 to 7. The height above the bed (hab) for the ADV deployments ranged from 0.14 to 0.5 m.
Sensors 23 02281 g002
Figure 3. Water column time series of the vertical distribution of the relative acoustic backscatter (RB) (using the average of the beams of each preset bin), corrected for transmission losses, from the acoustic Doppler current profiler (ADCP) at the (a) upstream (UP) and (b) downstream (DS) sites. Time series of echo intensity (EI) at the UP (c) and DS (d) show greater differentiation between depths.
Figure 3. Water column time series of the vertical distribution of the relative acoustic backscatter (RB) (using the average of the beams of each preset bin), corrected for transmission losses, from the acoustic Doppler current profiler (ADCP) at the (a) upstream (UP) and (b) downstream (DS) sites. Time series of echo intensity (EI) at the UP (c) and DS (d) show greater differentiation between depths.
Sensors 23 02281 g003
Figure 4. Calibration curves developed using the echo intensity (EI) from the (a) ADV and (b) ADCP at the upstream (UP), downstream (DS), and interior (INT) sites in the L−29 Canal (for the ADCP, two data pairs/calibration points were removed at UP, while a pair each was removed at DS and IN to improve the fit of the regressions).
Figure 4. Calibration curves developed using the echo intensity (EI) from the (a) ADV and (b) ADCP at the upstream (UP), downstream (DS), and interior (INT) sites in the L−29 Canal (for the ADCP, two data pairs/calibration points were removed at UP, while a pair each was removed at DS and IN to improve the fit of the regressions).
Sensors 23 02281 g004
Figure 5. Time series of discharge and estimates of TSS and TPP from the ADV deployments at (a) upstream (UP), (b) downstream (DS), and (c) interior (INT) sites. The TSS and TPP estimates are those near the canal bed. During the field deployments at UP and INT, the measured backscatter obtained from the ADV was lower than the measured backscatter obtained during the lab calibrations because of the inability to artificially create very low TSS concentrations in the lab (<5 mg L−1).
Figure 5. Time series of discharge and estimates of TSS and TPP from the ADV deployments at (a) upstream (UP), (b) downstream (DS), and (c) interior (INT) sites. The TSS and TPP estimates are those near the canal bed. During the field deployments at UP and INT, the measured backscatter obtained from the ADV was lower than the measured backscatter obtained during the lab calibrations because of the inability to artificially create very low TSS concentrations in the lab (<5 mg L−1).
Sensors 23 02281 g005
Figure 6. Vertical distribution of estimates of total suspended sediments (TSS) in the water column and the corresponding discharges at the upstream, UP (a,b), downstream in June, DS−Jun (c,d), downstream in December, DS−Dec (e,f), and interior, INT (g,h), sites in the L−29 Canal in 2021. The depths are the sizes of the user-determined bins of the ADCP. The uncorrected echo intensity (EI) from the acoustic Doppler current profiler (ADCP) was used to generate the TSS estimates via calibration curves. The calibration curve developed for DS−Jun was also used to estimate TSS for DS−Dec. For clarity, these plots are portions of the deployment time series.
Figure 6. Vertical distribution of estimates of total suspended sediments (TSS) in the water column and the corresponding discharges at the upstream, UP (a,b), downstream in June, DS−Jun (c,d), downstream in December, DS−Dec (e,f), and interior, INT (g,h), sites in the L−29 Canal in 2021. The depths are the sizes of the user-determined bins of the ADCP. The uncorrected echo intensity (EI) from the acoustic Doppler current profiler (ADCP) was used to generate the TSS estimates via calibration curves. The calibration curve developed for DS−Jun was also used to estimate TSS for DS−Dec. For clarity, these plots are portions of the deployment time series.
Sensors 23 02281 g006
Figure 7. Vertical distribution of total particulate phosphorus estimates (TPP μg L−1) in the water column at the (a) upstream (UP), (b) downstream in June (DS−Jun), and (c) interior (INT) sites in the L−29 Canal. The TPP from the lab-analyzed TSS (TPP in TSS) that was used to create the calibration curves (unit of mg g−1) was multiplied by the TSS estimated from applying the calibration curves (units of mg L−1) to obtain the TPP (unit of μg L−1) in the water column.
Figure 7. Vertical distribution of total particulate phosphorus estimates (TPP μg L−1) in the water column at the (a) upstream (UP), (b) downstream in June (DS−Jun), and (c) interior (INT) sites in the L−29 Canal. The TPP from the lab-analyzed TSS (TPP in TSS) that was used to create the calibration curves (unit of mg g−1) was multiplied by the TSS estimated from applying the calibration curves (units of mg L−1) to obtain the TPP (unit of μg L−1) in the water column.
Sensors 23 02281 g007
Table 1. Characteristics of acoustic devices manufactured by SONTEK—a Xylem brand.
Table 1. Characteristics of acoustic devices manufactured by SONTEK—a Xylem brand.
ParametersAcoustic Doppler Velocimeters
(ADVs)
Argonaut Acoustic Doppler
Current Profilers (ADCPs)
References
Sampling volume distance/Max profiling range5–18 cm40 m[32,33]
Max depth range60–400 m10–200 m[33,34]
Min blanking distance-0.07 m[35]
Bin size-0.2–30 m[35]
Max number of bins-10 + 1[33,35]
Pressure sensor0.1% accuracy0.1% accuracy[34,36]
Temperature Resolution, Accuracy0.01 °C ± 0.1 °C0.01 °C ± 0.1 °C[34,36]
Velocity accuracy, resolution1% of measured velocity, 0.01 cm s−1±1% of measured velocity, 0.1 cm s−1[34,36]
Min and max velocity0.1 cm s−1–5 m s−1<0.01 m s−1–6.0 m s−1[32,34,35]
Min Signal to Noise Ratio (SNR) for reliable velocity measurements5 dB3 dB[32,35]
Min correlation coefficient for reliable velocity measurements≥70%-[32]
Table 2. Total suspended solids (TSS, gravimetrically determined) and total particulate phosphorus (TPP) content of the TSS (TPP in TSS μg mg−1) and in the water column (μg L−1), with increasing depth from the water surface of the L-29 Canal. The TSS values were used to calibrate the echo intensity (EI) from the acoustic Doppler velocimeter (ADV) and the acoustic Doppler current profiler (ADCP).
Table 2. Total suspended solids (TSS, gravimetrically determined) and total particulate phosphorus (TPP) content of the TSS (TPP in TSS μg mg−1) and in the water column (μg L−1), with increasing depth from the water surface of the L-29 Canal. The TSS values were used to calibrate the echo intensity (EI) from the acoustic Doppler velocimeter (ADV) and the acoustic Doppler current profiler (ADCP).
SiteDateDepth (m) 1TSS (mg L−1) *TPP
(μg mg−1)(μg L−1)
Upstream (UP)24 November 20210.252.101.192.50
0.751.751.272.23
1.253.350.792.66
1.752.500.922.31
2.253.150.912.85
2.753.350.892.97
3.253.250.882.87
4.252.750.952.62
Downstream (DS)22 June 20211.24.50 *2.5011.25
2.25.402.2312.02
3.25.102.0910.64
4.28.201.7814.60
5.026.01.1529.87
Interior (INT)19 August 20210.51.203.333.99
1.01.103.664.03
1.51.203.504.20
2.01.303.574.64
2.51.703.265.54
3.01.603.255.19
3.51.802.825.07
1 The lowest depth corresponded to the depth of the ADV and was used to convert the estimated TSS (from calibration) to TPP estimates near the canal bed. * The calibration curve generated using the TSS sampled in June at DS was applied to a second deployment in December.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Onwuka, I.S.; Scinto, L.J.; Fugate, D.C. High-Resolution Estimation of Suspended Solids and Particulate Phosphorus Using Acoustic Devices in a Hydrologically Managed Canal in South Florida, USA. Sensors 2023, 23, 2281. https://doi.org/10.3390/s23042281

AMA Style

Onwuka IS, Scinto LJ, Fugate DC. High-Resolution Estimation of Suspended Solids and Particulate Phosphorus Using Acoustic Devices in a Hydrologically Managed Canal in South Florida, USA. Sensors. 2023; 23(4):2281. https://doi.org/10.3390/s23042281

Chicago/Turabian Style

Onwuka, Ikechukwu S., Leonard J. Scinto, and David C. Fugate. 2023. "High-Resolution Estimation of Suspended Solids and Particulate Phosphorus Using Acoustic Devices in a Hydrologically Managed Canal in South Florida, USA" Sensors 23, no. 4: 2281. https://doi.org/10.3390/s23042281

APA Style

Onwuka, I. S., Scinto, L. J., & Fugate, D. C. (2023). High-Resolution Estimation of Suspended Solids and Particulate Phosphorus Using Acoustic Devices in a Hydrologically Managed Canal in South Florida, USA. Sensors, 23(4), 2281. https://doi.org/10.3390/s23042281

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop