On the Potential Optical Signature of Convective Turbulence over the West Florida Shelf

: Atmospheric cold front propagation across the northern Gulf of Mexico is characterized by elevated surface wind velocities and a ~10–15 °C drop in surface air temperatures. These meteorological conditions result in significant heat energy losses from the surface ocean to the overlying atmosphere. These seasonally recurring cold ‐ air outbreak events may penetrate the southern portion of the West Florida continental shelf and initiate turbulent and convective overturn of the water column. Examination of true color images derived from ocean ‐ viewing, satellite ‐ based radiometer data reveals coincident and substantial surface water discolorations that are optically similar to smaller ‐ scale “whiting events,” despite the regional ‐ scale extent of the observed phenomenon (>25,000 km 2 ). Coupled air–sea numerical simulations suggest the surface water discoloration occurs and is sustained where the entire water column is dynamically unstable. The simulation results indicate significant density ( σ t ) inversions between the surface and bottom waters. Thus, the combined numerical model and remote sensing analysis suggest that convective turbulence may be contributing to the sustained ventilation of bottom waters containing a high concentration of suspended particulates. High ‐ temporal resolution true color images rendered from the GOES ‐ R Advanced Baseline Imager (ABI) data appear to support the surface water discoloration’s turbulent ‐ driven nature.


Introduction
Cold-air outbreak (CAO) events in the northern Gulf of Mexico (GOM) are known to result in substantial modulation of air-sea interactions and coastal ocean circulation as well as alteration of surface water masses [1][2][3][4]. CAO events are persistent in boreal autumn-winter but may occur during other seasons [5]. The southerly penetration of the continental, dry, and high-pressure air masses with accelerating winds out of the north has also been referred to as "northers" [6]. Each CAO event has varying degrees of southerly penetration into the GOM [7]. CAO events may only pervade the extreme northern GOM or they may penetrate as far south as Veracruz in southern Mexico [8].
The combination of elevated wind velocities and significantly reduced lower air temperatures during and following the initial passage of the atmospheric cold front drives substantial heat losses from the surface ocean [9]. The surface ocean buoyancy loss may then cause density-driven vertical mixing processes [10,11]. In a more general sense, such intense buoyancy losses result in turbulent fluid motions, that is, convective turbulence [12,13]. Turbulent convection in the ocean is often characterized by descending parcels of cold water that result from either diurnal cooling or seasonal cold-air outbreak events [12][13][14].
In regions of the continental shelf where this convective vertical ocean mixing penetrates down to the ocean bottom, the entire water column is then vertically wellmixed, i.e., homogenous with respect to density. Conceivably, significant and continuing heat losses at the surface under these well-mixed conditions may result in continuing density-driven instabilities and lead to sustained turbulent motions. Indeed, such a scenario for turbulent convection has been examined via numerical models intended to represent shallow waters off of the Florida coast during a CAO event [15].
Similar convective processes within shelf waters may resuspend fine-grained carbonates from the West Florida Shelf (WFS) seafloor during periods of intense buoyancy losses and vertical mixing. Accordingly, the term "whiting event" or "whitings" is often used to describe the elevated concentration of suspended calcium carbonate (CaCO3) particulates in discrete areas of the surface ocean. Whiting events are common on the Bahama platform [16,17], but the exact mechanism for their formation remains unknown. One common hypothesis is that some form of bottom sediment resuspension mechanism is involved and fish schools, tidal flows, small-scale turbulence, and even movement of water through subterranean networks have been proposed [16]. Convective turbulence is perhaps another mechanism that may resuspend carbonate particles, where such materials have accumulated on the seafloor. This paper combines advanced numerical air-sea simulations and coincident satellite images to elucidate the mechanistic functioning and the potential optical signatures of these CAO-induced turbulent mixing events. In the particular case of CAO events across the WFS, there is an evident and coincident discoloration of surface waters when examined using true color (TC) images that appears qualitatively similar to smaller-scale whiting events. Construction of TC images from the underlying radiometer data is an efficient way to broadly discern changes in surface ocean optical conditions that may not be as apparent when examining other products. For example, the phenomenon examined herein is not nearly as distinctive when analyzing satellite radiometer data in the context of standard, level-3, inversion products, such as the surface chlorophyll-a concentration. Precise determination of true color requires the application of ocean colorimetry, which is the quantitative determination of the ocean's apparent color based on the underlying spectral distribution of the emergent water-leaving radiance [18][19][20]. We employ colorimetry because it is a useful tool to rapidly delineate changes in surface ocean optics and potential delineations between water mass boundaries without any a priori assumption regarding the potential optical constituents and their respective impact on surface reflectance signals [18,21]. True color may also provide information that is lost in the standard inversion product fields but is nonetheless indicative of physical ocean processes and water mass boundaries [22]. Colorimetry also provides a useful framework for visible-band sensor fusion techniques, such as convolving GOES-R Advanced Baseline Imager (ABI) data with data from traditional ocean color sensors [23].

Satellite Data
Satellite datasets presented herein were obtained from three sensor types: (1) the Ocean and Land Colour Imager (OLCI) onboard the Sentinel-3A/B satellites; (2) the Visible Infrared Imaging Radiometer Suite (VIIRS) on the Suomi-National Polar-Orbiting Partnership (SNPP) and NOAA-20 satellites; and (3) the Advanced Baseline Imager (ABI) on the Geostationary Operational Environmental Satellite (GOES)-East. The OLCI and VIIRS level-1 sensor data were processed using the Naval Research Laboratory's (NRLs) Automated Optical Processing System (AOPS) [24]. The Sentinel-3A and B satellites have a staggered orbital coverage. The platform selected for the OLCI sensor analysis was whichever one provided the optimal spatial coverage of the WFS on that date.
For the TC scenes depicted in this paper, an intermediate processing radiance, ["rho_s"; ρs], was used to generate the color information. This product was selected because the cloud cover in some scenes made the full atmospheric (aerosol) correction difficult, and it was necessary to show the entire sequence of color images during the CAO event described herein. Using the intermediate product provides more data to view the full series of daily surface ocean discoloration images, particularly in the CAO sequence's initial phase when cloud cover is persistent. The ρs product is designated as quasi-surface reflectance in the NASA product documentation [25]. Standard level-3 remote sensing product data were rendered using the full atmospheric correction procedure in AOPS and are presented in Section 3.5.
TC fields shown in this paper were constructed following the colorimetric method [18][19][20]. Reflectance from the sensors at each visible band was interpolated to reconstruct the hyperspectral information (Δλ = 1 nm, 400-700 nm) for each pixel. The interpolated spectra were integrated with the CIE 1931 standard color matching functions [26] to produce the X, Y, and Z color primaries and then converted to the standard RGB color table for visual display as a three-channel jpeg image file. It is important to note that the conversion from X, Y, and Z primaries to standard RGB followed an established convention [27], and was not an arbitrary scaling of selected radiometer bands.
For the respective ocean color sensors (OLCI, VIIRS), TC fields were constructed from data obtained during each day corresponding to the time period of the numerical air-sea simulations (Section 2.2). Changes in the TC fields were detected first by visual inspection of the jpeg images. The mapped color gradients were then spatially compared to the physical simulation results. The initial analysis focused on heat flux terms and sea surface temperature (SST). A subsequent analysis then examined potential density and temperature inversions between the simulated surface and the bottom waters. Contour plots of the physical model results were overlaid upon mapped TC images. This procedure's motivation is that color gradients in coastal TC images often correspond to gradients in the physical oceanographic variables and may indicate the underlying physical oceanographic processes [22].
In this paper, we also examine visible band data from the GOES-ABI sensor. The geostationary GOES data have a very high temporal resolution (every 5 min) and can provide supplemental information for scenes where the surface ocean is detected. However, ocean color (that is, ocean surface water-leaving radiance detection in the visible) is not an intended application for this sensor. In order to repurpose GOES-ABI data to provide some continuity with dedicated ocean sensors, a method was developed to convolve GOES-ABI data with dedicated ocean color sensor information, and this procedure is described in Jolliff et al. [23]. Preliminary GOES-ABI level 1B CONUS data (~500 m horizontal resolution) were obtained from the NOAA Comprehensive Large-Array Data Stewardship System (CLASS), and processing details are described elsewhere [28]. GOES-ABI has two broad bands in the visible, one centered at 470 nm and one at 640 nm. ABI scene color information for the green wavelengths (centered on 555 nm) was estimated based on examining the coincident VIIRS and OLCI data. Under some circumstances, this method allows for GOES-ABI data to be used as a very high temporal visualization of surface ocean color changes.

COAMPS and Buoy Data
The Coupled Ocean-Atmosphere Mesoscale Prediction System (COAMPS) is an integrated ocean-atmosphere numerical modeling system [29,30] that provides operational support for the U.S. Navy [31]. The nested atmospheric and ocean models exchange information on fluxes of thermal energy and other data via the Earth System Modeling Framework (ESMF). The ocean model component, the Navy Coastal Ocean Model (NCOM) [32], was configured at 3 km horizontal resolution for the northern Gulf of Mexico. The ocean grid consisted of 50 vertical levels, with 36 sigma layers between the surface and 190-m depth (which are bottom-following in water shallower than 190 m), and 14 fixed-depth layers between 190 m and the maximum depth of 5500 m. Bathymetric data were obtained from the Navy's two-minute-resolution Digital Bathymetric Data Base (DBDB2). Initial and boundary conditions for NCOM were obtained from a global version of NCOM [30]. The atmospheric model inner nest was configured at 6 km horizontal resolution and 60 vertical levels. Atmospheric boundary conditions for each of the forecasts were provided by the Navy Operational Global Atmospheric Prediction System (NOGAPS) [33]. Additional details on model coupling, information exchange, and oceanatmosphere data assimilation can be found elsewhere [34,35]. The COAMPS simulation was executed for December 2018 with an internal ocean model timestep of twenty seconds. Simulation results were output every simulated hour for the duration of the model execution.
Meteorological and oceanographic buoy data were obtained from NOAA's National Data Buoy Center (https://www.ndbc.noaa.gov). The location of each buoy used is indicated herein. Some buoys were operated and maintained by the University of South Florida Coastal Ocean Monitoring and Prediction System (USF COMPS).

Meteorological and Physical Oceanographic Setting
The CAO event on 21 December 2018 provides a prototypical example of a cold front penetration to the southern portion of the WFS. Observed wind velocities and simulated wind stress peak on 21 December, with observed 10-m winds approaching 20 m s −1 before maintaining a steady direction out of the north, were consistent with previous CAO events ( Figure 1A). COAMPS simulates the propagation of the cold air mass and the ocean's response. The results maintain fidelity to buoy observations ( Figure 1B; Buoy 42023 [C13] is approximately ~95 km from shore and over the 50-m isobath). On 21 December, the airsea temperature difference (defined as the 2-m air temperature minus sea surface temperature [SST]) approaches −8 °C in the buoy data. The simulated temperature differences do not become this large (−6 °C). However, both the COAMPS results and the buoy data show a rapid increase in the magnitude of air-sea temperature differences beginning on 21 December and corresponding to the penetration of cold air to the southern portion of the WFS.
The simulated heat losses from the surface ocean also peak on 21 December and exceed 800 W m −2 . The latent heat flux losses calculated from the buoy 42023 data (not shown) also peak in excess of 600 W m −2 on 21 December and appear to be the dominant heat flux term during the CAO's initial phase. After the initial peak in simulated heat flux, the following four days (22-25 December) still show simulated heat losses in the ~400 W m −2 range ( Figure 1B). These heat losses are significant and comparable in magnitude to higher latitude sites of deep open ocean convection [36].
These buoy observations verify the COAMPS simulation results and confirm the characteristic features of a CAO event over the WFS. The COAMPS simulations then provide a more synoptic view of the regional event. On 20 December, the COAMPS simulated near-surface air temperatures (2-m) reveal the southwest to the northeast orientation of the progressing cold front (Figure 2a). Simulated latent heat flux values at that time range from nominal to ~200 W m −2 across the shelf (Figure 2b). The cold air mass pushes south rapidly, and this is shown in subsequent hourly results from the COAMPS simulation. By 21 simulated hours later, latent heat fluxes across the WFS are >400 W m −2, and lower air temperatures have broadly declined by ~10 °C (Figure 2c,d).

Coincident Satellite True Color
These air-sea interactions are simultaneous to observed synoptic changes in the TC fields across much of the WFS. Before the cold front arrives on 18 December, the OLCI sensor-derived TC image shows homogeneity of color between the shelf area and the open GOM, save for the locations very near the Florida coast (<25 km; Figure 3A). The following two days of OLCI-based TC imagery (21 and 22 December; Figure 3A,B) are more obscured by clouds but nonetheless reveal a sharp color discontinuity at mid-shelf, and much of the shelf appears turbid. The following three days of OLCI images ( Figure 4) are less impacted by clouds and show the persistent discoloration across much of the shelf area. The TC signal discontinuity then begins to dissipate on 26 December (image is not shown). However, the shelf discoloration's persistence appears to be synchronous with the period during which the COAMPS simulations suggest large heat flux losses (~400 W m −2 ) and sustained negative air-sea temperature differences ( Figure 1B).  Figure 1). The satellite platform (Sentinel A or B) was selected based on the best spatial coverage for the WFS on that date.
The VIIRS sensors (NOAA-20 and SNPP satellites) also capture the TC gradients across the shelf and the shelf-wide discoloration 21-24 December ( Figure 5). The VIIRS images emphasize both the temporal persistence and spatial homogeneity of the surface water discoloration. The VIIRS TC images on 24 and 25 December (Figure 4 B, C depict apparent parcels of turbid water exiting the WFS to the south and entering the Florida Straits. This apparent "streamer" of turbid water appears to be consistent with our hypothesis that the anomalous color images indicate surface turbidity, likely due to suspended particulates. The TC image was overlaid upon the DBDB2 bathymetry ( Figure 6). The anomalous water color appears to extend out to approximately the ~60-70 m isobathymetric contours and extend farther north over much of the WFS. The discontinuity in the color field corresponds to a sharp increase in the gradient of the bathymetry.

Merger of COAMPS Results with Satellite Images
COAMPS simulations provide the physical context for the shelf water discoloration. As the cold front penetrates across the shelf, the intense heat loss at the air-sea interface sets up a dynamically unstable regime in the water column. For example, differences between the simulated surface and bottom temperatures and densities (σt) on 20 December 1700 UTC show some temperature inversions (where the surface water temperature is colder than the bottom temperature) on the WFS developing shoreward of the 100-m isobathymetric contour (Figure 7a). These areas are consistent with developing density inversions (Figure 7b).  The temperature and density inversions shown in the COAMPS simulations were generally consistent with the patterns observed in the TC satellite images. Shown here, the COAMPS temperature inversions on 21 December (Figure 9a) were spatially consistent with the WFS water discoloration that was most clearly visible on 23 December in the VIIRS TC image (Figure 9b). The overlay of the temperature inversion contours upon the TC image shows that the discoloration was spatially consistent with the simulated instability (Figure 9c). This consistency suggests that the discoloration may result from convective turbulence occurring over the water column that results in either the direct resuspension of shelf particulates from the seafloor or the surface ventilation of bottom waters that contain a large concentration of resuspended particles.

High-Temporal Resolution GOES-R ABI Data
The discoloration of the WFS surface waters was not, in the initial stages, spatially homogenous. For example, an oblong region of the surface ocean, an approximately 147 km 2 area, remained discordant from the surrounding waters (OLCI-1527 UTC; Figure  10A). This area of the surface ocean had a color signal of dark blue that more closely resembled the open GOM's deeper waters, and it did not match the brightness color pattern seen in the surrounding surface waters. This area of the surface ocean was observed again by the VIIRS sensor a few hours later (VIIRS-1757 UTC; Figure 10B). In this subsequent TC image, there was no significant difference in the color pattern from the surrounding waters: the oblong color feature seen in the OLCI image a few hours previously was no longer discernable. This change suggests that the combined convective and wind-wave turbulence may result in some initial spatial inhomogeneity in the observed color pattern due to spatial differences in the resuspension of near-bottom particulate materials. As the turbulence persists, however, these inhomogeneities were removed over a short interval of time (~hours). The color-enhanced ABI image sequence (GOES-East, [23]) clearly shows the temporal progress of the presumed vertical mixing (Supplementary Video S1; Figure 11). The darker ocean area was discernable from 1400-1600 UTC ( Figure 11A,B) and became less distinct through 1800 UTC ( Figure 11C,D). The rapid hourly change in the surface water color properties over a large ocean area (~147 km 2 ) suggests vertical mixing was the dominant physical process, as opposed to horizontal advection and diffusion. This change is consistent with the hypothesis that the discoloration is due primarily to the vertical resuspension of particulates.

Examination of Level-3 Inversion Products
The optical distinctiveness of this surface water discoloration is not as apparent in some of the standard inversion products.   The standard CHL algorithm is based on remote sensing reflectance (Rrs) ratios [37]. The magnitude of Rrs in each waveband, neglecting wavelength and angular complexities, is proportional to the ratio of Inherent Optical Properties (IOPs) following: (1) where g is an empirical term to represent the wavelength and geometric complexities [38], bbt is the total backscattering coefficient, and at is the total absorption coefficient [39]. Under the hypothesis that the physical instability results in the turbulent resuspension of very fine-grained shelf carbonates, this would increase the backscattering coefficient [40]. Under the additional simplifying assumption that at >> bbt and Rrs  bbt/(at), then the substantial and proportionate increase in bbt in the spectral bands used for the CHL algorithms may largely cancel out when the band ratios are examined. Thus, empirical band-ratio algorithms may register a modest change, although the optical conditions have changed significantly. This concept is consistent with the observation that resuspended particulates that act as very efficient backscattering agents of downwelling light also tend to increase the amplitude of the hyperspectral visible signal of the emerging, waterleaving radiance, as opposed to significantly changing the spectral shape [41].
Examination of the hyperspectral signals interpolated from the OLCI multispectral data appear consistent with the hypothesis that increased particle backscattering was amplifying the Rrs signal; however, there is some degree of change in the spectral shape as well. Extracted Rrs spectra reconstructed from OLCI visible bands were selected for 9 pixels around the location of buoy 42023 and spatially averaged for 18 December and 24 December ( Figure 14). These spectra are compared to a representative turbid, coastal water spectrum obtained from the IOCCG reference spectra [42]. The 24 December spectrum is stronger than either of the other two spectra ( Figure 14A). Increased reflectance values would result in a very high brightness value when the signal is converted to a true color reference frame. Differences in spectral shape are more apparent when the spectra are normalized by their respective maximum values ( Figure 14B). The 18 December and 24 December spectra peak are near the 490 nm wavelength, however, the 24 December spectrum has more relative signal strength towards the green and red, whereas the 18 December has more relative blue (shorter wavelength) signal. Comparing the 24 December turbidity signal to the reference turbidity signal, the reference turbidity peak is shifted towards the red due to increased amounts of absorption from phytoplankton pigments and Chromophoric Dissolved Organic Matter (CDOM) in the shorter wavelengths, as is typical of coastal turbidity [43]. Consistent with the earlier examination of these different turbidity signatures [44], the 24 December turbidity peak is comparatively shifted towards the blue. This consistency suggests the 24 December "blueshifted" turbidity signal is due to a large increase in particle backscattering that is not matched by a commensurate increase in blue light absorption, as would be expected of coastal turbidity plumes from rivers and estuaries that are rich in phytoplankton pigments and CDOM. The band ratio chlorophyll algorithm [37,45] was applied to the spectra, and the resulting CHL values are indicated in Figure 14B. The CAO-induced turbidity results in a modest doubling of the chlorophyll product (chlorophyll tends to vary on a log-scale [46]), whereas the reference turbidity results in a satellite chlorophyll value of ~35 mg m −3 . Thus, using the band-ratio chlorophyll product as a general indicator of turbidity or water clarity would be misleading in this particular case.
In contrast to CHL, the IOP inversion products are more informative and support the general inference of increased backscattering across the WFS. The Quasi-Analytic Algorithm (QAA) [47] was used to examine the total backscattering and absorption coefficients in the 490 nm band. To accommodate an increase in 490 nm Rrs from 0.005 to 0.038 Sr −1 , one would require at least an order-of-magnitude increase in the backscattering coefficient if changes in absorption are marginal. Indeed, the QAA products support this inference ( Figure 15). Backscattering coefficients, as estimated by the QAA product, increased from ~0.005 to ~0.05 m −1 in the area of buoy 42023. The corresponding QAA total absorption coefficient products do not register nearly as dramatic of a change (Figure 16), although values are elevated on 24 December. The increase is not commensurate with the order-of-magnitude increase in the backscattering signal; thus, the reflectance values are accordingly elevated as well. This result also suggests the backscattering signal is primarily due to non-pigmented particles, that is, inorganic particulates. The likeliest source for such particles over such a vast area of the shelf is the underlying sediments.

Discussion
The surface ocean optical phenomenon examined in this paper is similar to the whiting events studied on the southern portion of the WFS [48]. The WFS whiting events examined by Long et al. [48] consist of <50 km 2 patches of discolored surface water. Some characteristics of smaller-scale whitings reported by Long et al. are very similar to what is reported here. Specifically, there appears to be no large and substantial difference in surface CHL, as estimated by satellite, before or during the CAO event for large areas of the shelf. Long et al. also found no significant elevation of observed surface CHL within smaller whiting areas as compared to the surrounding surface seawater. Long et al. also found that remotely-sensed reflectance values (Rrs) were substantially increased within the discolored surface waters. This increase in Rrs was matched by a corresponding increase in particulate backscattering, while the overall shape of the Rrs spectra was largely unchanged.
Nonetheless, the space-time scales of the phenomenon reported herein and those of Long et al. are substantially different. The discoloration described here is orders of magnitude larger (~25,000 km 2 ) and appears to be a regional shelf response to physical water column conditions. The discoloration case shown here lasts five days (21)(22)(23)(24)(25) and begins to dissipate thereafter; Long et al. report whiting patches on the WFS that last for weeks to months. The results shown here are more consistent with those of Conmy et al. [49]. These authors examined RGB (color) images constructed from MODIS data before and after a "winter storm event" in December 2004. They found that the post-storm event color brightness increased nearshore and attributed this color pattern to resuspended particulates resulting from elevated wind-driven turbulence.
Given these differences, it is less likely that the discoloration reported here is the result of either biologically mediated carbonate precipitation, as suggested by Long et al. for smaller patches of white water, or a synchronous and massive bloom of planktic cocclithophorids, which would also result in a similar surface water discoloration [50]. The regional-scale event described in this paper is more likely the result of a synoptic change in the water column's physical stability. The WFS is overlaid by largely unconsolidated sediments consisting of biogenic carbonates and quartz [51]. Biogenic calcium carbonate constitutes >75% of the total sediment weight for large portions of the WFS [52], an area that corresponds to the surface water discoloration shown herein (Figures 4-6). Biogenic carbonate particles, such as detached coccoliths from planktic coccolithophorids, have strong optical scattering properties [53]. The sediments also contain carbonate exoskeleton fragments that are eroded into progressively finer grains of carbonate [51]. Fine-scale carbonate particulates may remain suspended in the water column for an extended period of time and contribute to the persistent appearance of water column turbidity [48].
Thus, we hypothesize that the convective turbulence contributes to and sustains the persistent resuspension of fine-grained shelf carbonates, and this results in the water color discoloration seen in the satellite true color images as well as the increased backscattering coefficients and remotely-sensed reflectance. However, convective turbulence is not likely the sole mechanism of resuspension. Elevated winds and an excited surface wave field during the initial phase of the CAO event undoubtedly cause some resuspension in coastal areas through the combined effects of wind-driven currents, wave action, and Langmuir circulation [54]. Also, a downcoast or northerly wind stress will result in an upwelling mode of coastal circulation along the WFS [55]. The combination of upwelling-favorable wind stress and the breakdown of any shelf stratification may result in the elimination of discrete frictional boundary layers between the surface and bottom and then result in turbulent mixing.
Nevertheless, simple scaling arguments do lend support for the dominance of the convective turbulent mechanism over the extended duration of the CAO event. Ocean turbulence generated by surface winds decays away from the ocean surface as a function of depth [56]. Whereas this may be a sufficient source of bottom turbulence in nearshore areas (<~25 m depth) to resuspend sediments, it is a less satisfying explanation for sustained turbulence over deeper regions of the shelf (~60 m). The burst of intense wind stress corresponds to the initial propagation of the front (21 December). During the remaining period of the CAO (22)(23)(24)(25), winds are less active ( Figure 1A). In contrast, for areas where the water column density profile is unstable, turbulent energy is proportional to the surface buoyancy flux all across that unstable layer with little attenuation with depth [14].
Thus, we hypothesize that the surface water discoloration during the CAO event may be sustained by convective turbulence where such turbulence penetrates to the ocean bottom. This hypothesis would further suggest that the TC fields or the IOP backscattering coefficient products may indicate where this turbulence is occurring on the shelf during CAO events. As demonstrated herein, the process of combining TC images derived from ocean color sensors with numerical air-sea simulations can illuminate potential mechanisms that lead to observed surface ocean color patterns as well as identify phenomena that may be overlooked if one is solely examining standard level-3 ocean remote sensing products such as the chlorophyll-a concentration.
A more detailed spatiotemporal budgeting of turbulent energy is being analyzed from the COAMPS results during the CAO event for a follow-up paper; however, further testing of this hypothesis will require fieldwork and mooring deployments that occur during a CAO event. The focus within this paper has been the December 2018 CAO case study. However, CAO events are common and appear to be associated with similar episodes of shelf water discoloration. Under the presumption that the discoloration results from sediment resuspension, this suggests that CAO events may have critical biogeochemical consequences for shelf waters. In addition to particulates, resuspension events may release nutrients trapped in pore waters that may lead to subsequent biological productivity [57].

Conclusions
The correspondence between the simulated shelf-wide density inversions (σt) and the sustained surface ocean water discoloration strongly suggests that convective turbulence over the WFS areas shoreward of approximately the ~60-m isobathymetric contour contributes to the substantial water column turbidity through either the resuspension of bottom sediments or the ventilation of the bottom waters rich in suspended particulates. We hypothesize that where significant heat losses drive convective turbulence across the entirety of the water column, resuspended particulate materials consisting largely of carbonates create a distinct discoloration of the surface waters. Thus, areas of full-water column convective overturn and mixing may be detected with visible remote sensing methods during a CAO event, and true color images are particularly informative. Confirmation of this hypothesis will require in situ observations of surface water optical properties, above water radiometry, and analysis of the constituent particulate materials coincident with the satellite observations during an episodic CAO event. Further numerical modeling studies with sediment models integrated into the air-sea modeling system (COAMPS) would also elucidate the potential mechanisms leading to the observed phenomenon.