The Use of SF 6 and GIS to Study Farfield Modeling of Ocean Outfall Plumes in Florida

In June, 2004 and February, 2007, in field tracer studies were conducted on the Hollywood and South Central outfalls, using sulfur hexafluoride (SF6) as a tracer. The objective of these studies was to determine if the tracer could be detected in the farfield at significant distance, and if so, could this data be used to construct a model of the farfield plume. Prior models for farfield plume movement do not appear to comport well with the conditions in southeast Florida. Extensive research was conducted in southeast Florida on 4 outfalls, which led to the development of nearfield dilution equations for same. However farfield modeling of outfall plumes was difficult to accomplish because the tracers used are not detectable for significant distances. The SF6 resolved that problem and as a result the Hollywood outfall was used to construct a model. Two methods were investigated for modeling the plume, 1) the Eureqa formulation method and 2) the Gamma-Curve method. The concentrations in the x-y plane were first found by using the Eureqa formulation to calculate the concentration at each grid point given its depth and the concentration of the centerline at the same latitude. The plume models were generated using MATLAB that matched with the results actually seen in the field.


Introduction
Six utilities in Palm Beach, Broward and Miami-Dade counties were discharging 396 million gallons per day (MGD) of wastewater in 2005.In 2009, the Florida Legislature passed a bill to cease outfall operations by 2025.In 2011 this deadline was modified slightly.The proposed solution was to reuse the wastewater.However how over 300 MGD of brackish water was to be reused was not answered.Instead more questions have been raised about the actual impact of the outfalls as the impact of policies implemented for the coastal ocean may have profound impacts on other activities as such decisions are closely linked to actions to be undertaken with respect to water and energy supplies to the rapidly growing population of southeast Florida.These costs are expected to excel $20 billion dollars and require over 2 Giga Watts of new power (proposed as nuclear reactors on Biscayne Bay [1]).
The State of Florida is presently faced with critical decisions regarding the future of its coastal environment as a result of terrestrial activities.Coastal ecosystems, including both inshore estuaries and offshore reefs, sustain the State of Florida's natural resource driven economy [2].Both must be protected to continue providing economic opportunities critical to the State.Increased nutrient loads could promote the growth of algae on the reefs, harming their habitat, their dependent wildlife and affecting tourist and commercial activities.Implication of anthropogenic causation have been suggested in the overgrowth of macroalgae and/or benthic cyanobacteria blooms which cause loss of coral and sponge cover and changes in reef fish demographics [3].In 1994-1995, a blooms of codium ishmocladum were reported in reefs off of Broward and Palm Beach Counties [4], and Caulerpa brachypus C. brachypus var.parvifolia in 2001 [5] [6].Lyngbya spp blooms were observed off of Broward County in 2003 [7].Blooms of Cladophora liniformis, Enteromorpha prolifera, and Centroceras clavulatum were observed in the spring of 2007 [8].
While it is asserted that a significant fraction of the nutrients impacting the coastal ecosystem are anthropogenic [9] [10], the current state of understanding concerning point and non-point sources of pollution, as well as natural process such as upwelling, is incomplete [11].Prior efforts by the authors has indicated that the canals and upwelling are significant courses of nutrient deposition [12]- [19], but other sources could include such far-field processes as coastal pollutants from the Florida Keys, the Caribbean, the Mississippi River, the southwest Florida shelf and atmospheric deposition, depending on the prevailing oceanographic conditions and location.Differentiating which, if any of these discharges may be the source(s) of nutrients is nearly impossible.Despite these numerous sources of nutrients, the outfalls have long been a target of divers as the cause of the algal blooms, despite the lack of data to implicate the outfalls [18] [20].
South Florida has six outfalls discharge (see Figure 1).Each of these outfalls discharges secondary treated wastewater at a depth of 100 feet below the surface (see Table 1).Prior researchers have developed equations to describe the nearfield dilution characteristics through the efforts from the Southeast Florida Ocean Outfall Experiments (SEFLOE) conducted from 1984-1994 [20]- [22].As a result, the southern four outfalls are among the most studied outfalls in the world.
Among the findings of SEFLOE were that the secondary treated wastewater is less dense than seawater and hence rises as a buoyant jet above the outfall (see Figure 2-acoustic photograph of the Broward outfall).The region of the rising plume where mixing is caused predominantly by buoyancy effects, is known as the nearfield region.Dilution exceeds 20:1 for all four outfalls before reaching the surface during the 10 th percentile current condition [20].Once the plume reaches the surface, it spreads out at the sea surface (top 1/3) as a less salty plume above the denser seawater as predicted by Richardson [23].The surface region is termed the farfield, which is the result of further mixing that occurs due to ambient ocean currents [24].Vertical transport is usually considered to be small, as diffusion effects are orders of magnitude smaller than advective horizontal transport [25].
Tracking nutrients from the ocean outfalls in the farfield is difficult as the nutrient concentrations are very low.For example, the total nitrogen concentration is typically less than 20 mg/L in the discharge.If it receives a 20:1 or more dilution in the rise of the plume to the surface, the total nitrogen concentration will be below 1 mg/L which is hard to track for any distance on the surface.During the SEFLOE effort, only entorocuccus was detectable at 800 m, while no other parameter was detectable at 400 m [20].Given that the plume rise to the surface, and tracking nutrients is difficult at the surface due to dilution, impacts to the reefs on the bottom are not anticipated in the nearfield vicinity of the outfalls.
A variety of investigators have made suggestions on modeling outfall plumes.Koh and Brooks [26] equated   outfalls to underwater chimneys with an emphasis on particulate deposition.They note that "it is practically never feasible to provide such complete treatment that an outfall is not necessary, nor is it common to depend entirely on an outfall with no treatment."As a result, modeling the fate and transport is assumed to follow a Gaussian dispersion which continually deforms under applied shear stress to generate a time-averaged pollutant concentration relationship between that can be modeled using a normal, Gaussian, distribution curve along its longitudinal axis [26] [27].
There has been much work on nearfield dilution, but far less on the far field.Initial plume models were developed for the nearfield rise based on thermal principles by Trent and Welty (1973), and subsequently improved on by others using buoyancy.Eulerian integral models simulate buoyant jets based on a set of control volumes fixed in space [28]- [33].Lagrangian models have also been used by treating buoyant jets as a series of non-interfering moving elements [34]- [36].Other efforts have been undertaken by Grace [37] and Carvalho et al. [38].Economopoulou and Economopoulou [39] [40] developed a graphical solution for outfalls that are useful in the nearfield Others developed or used software: VISJET [41], Delft3D and the Cornell mixing zone expert system (CORMIX [42]) and Visual Plumes [43] [44], Huang et al. [21] [22] presented a probabilistic approach to develop initial dilution criteria or standards for ocean outfall design in the nearfield as a part of SEFLOE, while.Lye and Sharp [45] used a Monte Carlo Simulation to predict the ability of outfall to meet compliance requirements in the nearfield.
With respect to farfield outfall studies, a series of investigators attempted to model the farfield without tracers [46] [47].Several found tracers to be not viable in Florida [48]- [52].Gray [53] investigated the impact on Australia's outfalls on fish populations, but detected no numerical impacts on fish larvae, which is not inconsistent with the finding of Bloetscher et al. [18] for Florida's outfalls.Outfall plumes that do not follow the flow regime of the Florida outfalls or their currents include plumes modeled by many others [54]- [66].As a result, there are no prior studies that occur in conditions similar to SE Florida, and due to cost to conduct tracer experiments, a limited data sample was available for this effort, fully acknowledging that studying the dispersion of effluent from coastal ocean sewage outfalls has continued to be a difficult problem to study in situ because the conditions are constantly changing as Jones, et al. [67] noted.

Methods
Given highly diluted environments in the nearfield, defining farfield effects requires the ability to measure very small quantities of highly diluted contributions on the bottom of the ocean floor.Therefore, to resolve the problem of detecting plume movement in the farfield, NOAA proposed a conservative tracer to prevent impacts on the reef community whenever the terrestrial plume touched the bottom.The indicator chemical that could be detected in-situ at very low levels was sulfur hexafluoride (SF 6 ).SF 6 was suggested as the tracer due to NOAA's success in using it in the open ocean.Its usefulness in the coastal ocean had not been tested until the Hollywood experiment discussed here in 2004 [68].The detection limit of SF 6 is less than 1 part per trillion by volume (pptv) using headspace analysis, and two orders of magnitude lower still using trap and purge techniques [69] which has made it possible to perform tracer studies over large space scales and long time scales.The continuous release was performed to create a quasi-stationary tracer field in the ocean.Wannikhof et al. [68] notes this application adds to the utility of the tracer, which had only previously been used as a tracer of gas exchange processes in the ocean [70]- [72], lakes [73] [74], and in rivers [75]; as a surface water dispersion tracer [76], as a tag for surface water in Lagrangian studies ([77]- [79]; and to determine mixing in the deep ocean [80] [81].Because Wannikhof et al. [68] noted success with the Hollywood tracer in June 2004, a second tracer study of the Boynton outfall was conducted in 2007.The intention was to use these data points to develop a model of the farfield plume to evaluate the potential effects of the outfalls on benthic communities. The Hollywood outfall is located approximately 2.5 km offshore and approximately 8.5 km south of Port Everglades.The outfall terminates at a depth of approximately 30 m and discharges roughly 46 MGD of treated waste water through a 60 inch diameter armored concrete pipe.A conservative tracer sulfur hexafluoride gas (SF 6 ) was injected via bubbling into the wastewater stream.SF 6 was used because it is a good tracer to study point source releases as it has low background levels in the environment.It is nontoxic, stable, and measurable at low concentrations [68].The injection of SF 6 began on June 4, 2004 (10 a.m.EDT) and ended on June 9, 2004 (3 p.m. EDT) as described in Wannikhof et al. [68].
The ambient current velocity of the ocean during the sampling, data was available from two different locations.One was a current profiler located 5 km downstream of the outfall at 11 m depth.Here, the current was predominantly northward with little vertical shear and with a flow of about 20 cm•s −1 .Another measurement of northward velocity was taken 12 km offshore just north (downstream) of the study region with a flow of 80 cm•s −1 .For the development of this model 20 cm•s −1 was used.In the offshore, there is no noticeable difference in salinity [82] [83].
The continuous injection of SF 6 into the outfall pipes was performed to create an environment of elevated SF 6 concentrations along the predominant flow path of the plume to facilitate tracking of the farfield SF 6 levels.The sampling survey of the ocean waters commenced on three days after initial dosing commenced on the Hollywood outfall in a continuous mode using a towed "fish" containing a Seabird Micro Cat C (onductivity)-T (emperature) recorder and submersible water pump and samples as described in Wannikhof et al. [68].The same basic procedures were used in February 2007 for the second SF 6 tracer experiment conducted by NOAA and Florida Atlantic University (FAU).
Garvine [84] suggests that the plumes may be modeled and characterized by using averaged conditions to develop a simpler model as opposed to partial or stochastic differential equations or may be specified only probabilistically.Using this theory, the constituent was assumed effluent be released in steady state to the ocean surface at a rate (M) and at a location (x 0 , y 0, z 0 ) as defined by Huang et al. (1996), where the effects of advection are assumed to be significantly larger than the amount of dispersion in the x direction.Pire-Schmidt et al. [13] used limited data to develop a useful predictive model of an inlet plume, based on the same tracer study methods using the 2007 data.The plume was successfully modeled with a Gaussian plume model that appears to mimic the response.It was noted that the tracer concentrated in a series of boluses that migrated north of the inlet [13], which is not expected because the inlet pulses with the tidal cycle.
The dispersion coefficients σ y and σ z are the standard deviations of the horizontal and vertical Gaussian distribution respectively [27].As such, the dispersion coefficients represent the area under a Gaussian curve that encompasses 68% of the variance present in the distribution.The most common procedure for estimating the dispersion coefficients was introduced by Pasquill [85] and modified by Gifford [86].In ocean plume modeling, dispersion coefficients are the driving force behind the advection-dispersion equation.This equation can be solved using a finite element method that is forward in time and central in space for dispersion.Using finite elements, the solution of a 2D diffusion equation for an instantaneous point source of mass per unit length M/L released at time t = t 0 and location (x, y) = (x o , y o ) in an unbounded domain is: ) The value of σ y and σ z are based on units of (m 2 •s −1 ).Equation ( 1) can be used to predict the concentration at a distance x, y, and z away from the point source given all parameters.Finding these parameters however, is an art of its own requiring professional expertise, mathematical formulation, and computer modeling.Shear and horizontal dispersion as discussed by Okubo [87]- [91] were reviewed, but the conditions did not appear to be replicated for the outfalls.This work will therefore only use the principles of ADE to model the plume equation in the x-z plane.The x-y plane will be modeled using data-fitting software to derive an empirical equation based on the Hollywood outfall data

Creating the Base Map of Measured Concentrations
Using the same assumptions as Pire-Schmidt et al. [13], and Pire-Schmidt [12], the process involved taking the data sets, and mining them for information.As noted during the 2007 SF 6 event [13] two sources of bias exist in the sample data: the SF 6 adhering to sampling hardware on eastward (outbound) sampling tracks which results in artificially increased SF 6 concentrations on eastward tracks and reflection of the SF 6 pollutant from the southeast Florida coast which lies to the west of the outfall location.Pire-Schmidt [12] outlined how this bias was addressed, but generally the following assumptions were used: • The SF 6 adherence issue biases maximum SF 6 concentrations towards the east.
• The coastal reflection issue biases maximum SF 6 concentrations towards the west.
• The biases are generally of equal magnitude.
• Rotating east-biased sample points westward and west-biased sample points eastward will sufficiently remove biases.• A coordinate of the centerline of the SF 6 tracer plume in the Hollywood outfall experiment is generally lo-cated at the point of highest concentration along a given latitude (after sample points are bias-corrected).
• The centerline of the SF 6 tracer plume in the Hollywood outfall experiment has its origin (passes through) the actual Hollywood outfall location (i.e. the origin of the Hollywood tracer centerline is not significantly affected by the hydrodynamic phenomena at Port Everglades).• Beyond 13 kilometers north of the Hollywood outfall origin, the tracer plume is completely mixed in the water column.

Spatial Mapping of Experimental Data
The primary tool for generating spatial interpolation maps was ESRI Arc Map v.10.0 software using its Geostatistical Analyst Tool.Data sets originally provided in a spreadsheet form were imported as points into the Arc-Map software environment.Spatial data sets were also projected to the NAD 1983 HARN State Plane Florida East FIPS 0901 coordinate system.
To map the SF 6 concentrations, the Kernel Smoothing interpolation with barriers method under the ArcMap Geostatistical Analyst tool was used.The Kernel Smoothing method is a variant of a first-order local polynomial interpolation [92] [93] and allows for the incorporation of barrier features (i.e. the Florida coastline) into the interpolation.This interpolation method was selected as results appeared to capture the theoretical north-south anisotropy of the tracer plume better than other available interpolation methods.The biases noted above were corrected by separated the incoming (west biased) and outgoing (east biased) data sets and rotating them through a common origin.Initial interpolation algorithms were applied to the rotated concentration data using Inverse Distance Weighting (IDW) and Kriging, with Kriging being the most typically used algorithm for plume modeling as determined during the literature review.The predictions created by the IDW algorithms produced results that when viewed using a dozen classifications bins highly indicated the presence of the tent pole effect around the data points Kernel smoothing was the applied to achieve the best Root Square Mean Error (RSME).
To map the seafloor depths in the Hollywood outfall study area, the Ordinary Kriging interpolation under the ArcMap Geostatistical Analyst tool was used.Bathymetry is an important factor due to the reflection caused by the seabed.After the generation of interpolation maps for SF 6 concentrations and seafloor depths in the Hollywood region, results were exported to a rectangular grid of 2368 data points, covering an area of approximately 8.5 × 10 8 square meters that reflected the actual data gathered during the experiment.

Creating a Model
Data on current velocity was gathered from bottom instruments.Coastal inlets also were investigated to see if tidal periods had an effect on determining the centerline equation from interpolation.The reason why inlets are of concern is that during tidal periods large quantities of water are channeled through the inlet as a gradient change.Specifically the velocity component of the tide is of concern to this project.The underlying assumption is that advection which is driven by velocity will cause the spread of the tracer in the direction of the velocity.When considering South Florida, the initial velocity travels parallel to the coast in a northern direction.When peak tide periods are present the possibility of an additional velocity component acting transversely to primary advection current can significantly alter the direction of the centerline plume.To determine if tidal influence could be an influencing factor on the centerline and concentrations an evaluation of the location of test sites to proximity of inlets were made.

Model in the X-Plane
As noted above, it was suggested that under perfect conditions, the behavior of the plume is such that it can be characterized as having a Gaussian distribution in respect to its longitudinal axis with the highest peak at the centerline, while exponentially decaying along that same axis.The first step was to develop the centerline, taking into account the reflection is caused by the plume being trapped at the surface, hitting the shallow ground and/or the westbound coast and reflecting off those surfaces creating a higher concentration than what is actually present.The highest concentrations were mapped by the latitude and longitude.Ultimately, the coastline reflection effect was handled by deleting the westbound samples, and assuming that the concentration west of the centerline (towards the shoreline) would have a concentration equal to or greater than that of the centerline for the same reasons as outlined in Pire-Schmidt [12] and Pire-Schmidt et al. [13].This implies that the effects of the reflection sources are proportional to depth; hence, depth is a parameter in the model.To determine the magnitude of the effects due to image sources, one can simply reflect the left side of the Gaussian distribution to generate a graph that would exponentially increase the concentration at shallower depths.

X-Y Plane
Conceptually, the concentrations in the x-y plane were found by modifying the "ideal" plume to account for the different factors that alter the dispersion/diffusion of it.In practice, concentrations were found by calculating the concentration, using the data-fitting modeling software Eureqa v.0.95 at each grid point given the depth and the centerline concentration associated with the grid point.The Eureqa v.0.95 application was developed by the Cornell Creative Machines Lab [94] as a tool that searches mathematical relationships between input and output variables.Assumptions used for modeling purposes were as follows: • The initial dilution S, is the dilution (S o ) found at the end of the near-field regime (L b ), which is assumed to be 40:1 based prior modeling using and average current speed of 0.2 m/s.• The average current velocity u a , was 20 cm/s or 0.2 m/s (based on estimated, average, current meter data).
• Vertical diffusion has minimal affect on the plume until the plume and ambient waters are nearly identical, and therefore can be ignored in the initial phases of the farfield.For purposes on the analysis this will be assumed to occur when S >100:1 and beyond 13 km from the outfall discharge point.• Initial plume thickness remains a thickness of 0.3 of the total depth for some undefined distance (<13 km), but becomes a complete mix system after that point.• Results of SEFLOE II are used where applicable.
The actual data points were used in Eureqa.Identification of a mathematical relationship between these variables assisted with the characterization of pollutant plume dispersion in the lateral (east-west) direction without the need for determining horizontal dispersion coefficients directly.The relationship search was simplified by allowing seafloor depth to serve as a surrogate for east-west distance of a point from the Florida coast, based on the high positive correlation between the shallowness of seafloor depth and nearness to the southeast Florida coast.The ideal equation minimizes complexity while still exhibiting optimal performance in various statistical error and correlation measures, which obfuscates the need to determine the exact relationship between a variable and the system.The concentration at a point (x,y) on the grid changed by C x,y -C xo,yo regardless of the mechanism.

X-Z Plane
Evaluating the x-z plane of the outfall plume requires information about the general behavior of an effluent under buoyancy-dominated conditions.At the outfall pipe, the jet stream will bend upwards within a short distance from the outlet due to the difference in buoyancy between the effluent and the surrounding environment.Once it hits the surface, the plume is advected with the ambient velocity field at the speed u of the surface plume layer [21].The subsequent dilution/dispersion of outfall plumes is essentially governed by buoyant spreading and oceanic turbulent diffusion.Results from the SEFLOE field tests indicated that the dilution/dispersion of the Hollywood outfall plume was dominated by buoyant spreading over a range of several hundred meters from the outfall [20] [21].It is also known that the plume will continue to spread in 2D in the x-z plane until it is fully mixed along its water column.Based on cast data during the 2007 experiment, 13 km downstream of the outfall the water column appears to be completely mixed.Since the concentration at 13 km is known at the centerline (y = 0), the z-dispersion coefficient can be back calculated.

Hollywood Model
The centerline model was developed as follows based on the data shown in Figure 3: Figure 3 shows the plot of this equation versus the data points.
Figure 4 shows the Eureqa plot for the equation to describe the horizontal diffusion:  ( ) 2.51 , 0, 0 0.451 0.00845 Both predicted and observed values of SF 6 concentrations were plotted along a 45 degree line to provide a graphical illustration of how the model is performing.A review of the graph indicates that the plotted points generally follow the path of the 45 degree line.This suggests that the model is performing relatively well.More points appear on the left of the 45 degree line which is an indication that the model is slightly overestimating SF 6 concentrations.
A line plot was also done for the analysis (see Figure 5).The graph shows how the predicted values (red), matches those of the observed values (blue), indicating a good fit.
The software program MATLAB was used to generate a boxplot of the residual errors (see Figure 6).The plot showed four (4) points which could be considered as outliers.Removing these 4 points reduces the RSME value from 27.77 to 13.13 and the MAE value from 17.08 to 10.70.The statistical performance measures asso-  ciated with this formulation are presented in Table 2.

X-Z Plane
Using the fact that based on the Hollywood cast samples, at 13 kilometers north of the Hollywood outfall the ocean plume the SF 6 had completely mixed along the water column, the dispersion coefficient in the z direction was back calculated.Additionally, the concentration values at the centerline were plotted against distance and an exponential regression line of best fit was drawn.A trial and error method was employed to find the dispersion coefficient in the z direction because it was known that at 13 km, the ocean plume was completely mixed in the water column.These centerline concentrations were the inputs to obtain the following expression for the x-z plane: ) Figure 7 shows a comparison of the SF 6 concentrations from the Hollywood outfall plume obtained using two  different methods: the Kernel Smoothing method and the mathematical model.Both tend to follow the same pattern as far as surface concentration distribution is concerned.Additionally, the magnitudes of the concentrations seem to be consistent with one another.The mathematical model however displays series of boluses that appear along the longitudinal axis with decreasing magnitude.It is uncertain what the cause to these boluses but they may be caused by the different depths present at the ocean bed since the Kriging method does not take depth into account when formulating its graph.The stars represent the location of the ocean outfall.An error analysis was performed for the calculation of the concentration plume for the Hollywood outfall in the x and y directions.No error analysis was done for the z direction due to limited sampled data points.A grid system was used for comparison.Table 3 summarizes how the model correlates to the raw data points.

Conclusions
The objective of the plume modeling study was to 1) determine if a farfield tracer could used for data gathering and 2 ) if that data would permit modeling and calibration of a farfield, coastal ocean plume that adequately describes the advection and dispersion of a near-coastal conservative pollutant plume in the southeast Florida region.The following are the conclusions that can be gleaned from the experiment: 1) SF 6 can successfully be used for monitoring terrestrial discharges to the coastal ocean in the farfield.The tracer captured data ≈ 10 to 60 km downstream.
2) From 20 to 60 km northward of the outfall maximum concentrations decreased less than twofold.
3) To track the far, farfield, the SF 6 that may be required may overwhelm detection instruments in the nearfield.
4) Complete mix of terrestrial discharges occur less than 13 km downstream of the discharge terminus.5) Terrestrial activities such as the port is a complicating factor for coastal ocean discharges, especially if they are less than 13 km from the terminus as their impact on complete mixing is unclear.
The development of the plume model in this study used data gathered from the experiment conducted by Wanninkh of et al. [68] on the Hollywood outfall as well as data other NOAA operations [83] [95].
The results of the modeling showed that the data could be translated into a model.The magnitudes of the concentrations seem to be consistent with one another.The error analysis performed on the Hollywood model (Eureqa method) indicated that the model is performing exceptionally well for that region.Four performance measures were used to evaluate the model; Root Mean Squared Error (RMSE), Mean Absolute Error (MAE), Absolute Error (AE) and Correlation Coefficient (ρ).An analysis was also done of the gamma-curve method using the same procedures.The models showed a slight inclination towards overestimating the SF 6 concentrations.
There were several issues that were discovered through the experiment that need to be dealt with in the future.A complicating factor was Port Everglades.The tidal considerations create the potential for the change from a continuous flow stream to become a pulsed stream as a result of SF 6 moving to the port during incoming tides.Hence the boluses.No samples were taken to identify the variation in SF 6 in the port, but a repeated experiment should be undertaken to study this issue.
The impact of trapping as a result of the Florida Current being east of the outfall discharge, differential current velocities below the plume and the complete mix assumptions.A predictive Bayesian exercise would appear to be worthwhile.
Finally wind is an issue not addressed during the experiment.The wind was shoreward.As a result, the off-gassing SF 6 could have been reincorporated into the surface waves, creating a false footprint of the plume, which is one reason the west side of the data was ignored.Hence the SF 6 should not be used as the only tracer in future studies intended to trace a specific point source.
Varying physical conditions in the ocean mean that the experiments should be repeated to demonstrate veracity and robustness of the model.A future experiment, under different current speeds is suggested.However if such an experiment is conducted, the following should be addressed: 1) For quantitative interpretation, current (ADCP) measurements and wind data must be collected in conjunction with the SF 6 injection to determine the coastal current conditions.
2) A second, conservative, non-gaseous tracer should be included along with the SF 6 to get an additional confirmation of the nutrient sources.
3) Reconnaissance of discharge pipelines should be undertaken to insure that there are no leaks and that appurtenances do not leak SF 6 gas.
4) Background information on significant current regime disruptors, like Port Everglades, should be gathered and characterized along with other discharges.
One major finding was that dilution in the farfield (60 km) was extensive (greater than 10,000:1), and within

Figure 1 .
Figure 1.Location of ocean outfalls in SE Florida.

Figure 3 .
Figure 3. Exponential regression line of CL values for Hollywood outfall.

Figure 5 .
Figure 5. Line plot of predicted and observed values for Hollywood outfall.

Table 2 .
Error measures for selected horizontal dispersion formulation.

Table 3 .
Summary statistics for model: Hollywood model (Eureqa) error analysis.