West Florida shelfcirculation and temperature budget for the 1999 spring transition
Robert H. Weisberg*
College of Marine Science, University of South Florida
140 7th Ave. South, St. Petersburg, Florida, 33701
Revised: August 2001
*Corresponding author. Tel: 1-727-553-1568; Fax: 1-727-553-1189
E-mail address: email@example.com
Mid-latitude continental shelves undergo a spring transition as the net surface heat flux changes from cooling to warming. Using in-situ data and a numerical circulation model we investigate the circulation and temperature budget on the West Florida Continental Shelf (WFS) for the spring transition of 1999. The model is a regional adaptation of the primitive equation, Princeton Ocean Model forced by NCEP reanalysis wind and heat flux fields and by river inflows. Based on agreements between the modeled and observed fields we use the model to draw inferences on how the surface momentum and heat fluxes affect the seasonal and synoptic scale variability. We account for a strong southeastward current at mid-shelf by the baroclinic response to combined wind and buoyancy forcing, and we show how this local forcing leads to annually occurring cold and low salinity tongues. Through term-by-term analyses of the temperature budget we describe the WFS temperature evolution in spring. Heat flux largely controls the seasonal transition, whereas ocean circulation largely controls the synoptic scale variability. These two processes, however, are closely linked. Bottom topography and coastline geometry are important in generating regions of convergence and divergence. Rivers contribute to the local hydrography and are important ecologically. Along with upwelling, river inflows facilitate frontal aggregation of nutrients and the spring formation of a high concentration chlorophyll plume near the shelf break (the so-called "Green River") coinciding with the cold, low salinity tongues. These features originate by local, shelf-wide forcing; the Loop Current is not an essential ingredient.
Keywords: Circulation, Temperature budget, West Florida Shelf, Modeling and observations
Located at the eastern edge of the Gulf of Mexico, the West Florida Continental Shelf (WFS) is one of the broadest continental shelves in North America. Between its southern and northern ends, bounded by the Florida Keys and the Florida Big Bend, respectively, the WFS isobaths vary smoothly, and they generally parallel the coastline. This geometry changes along the Florida Panhandle in the north where the coastline undergoes a right angle bend, and the shelf width decreases to a minimum at the DeSoto Canyon.
Long-term observations (Weisberg et al., 1996) show that the WFS circulation, forced by tides, winds, buoyancy, and possible interactions with the Gulf of Mexico Loop Current, varies on time scales from semi-diurnal to inter-annual. Monthly mean currents at mid-shelf suggest a seasonal cycle with along-shore flows either to the southeast in spring, or to the northwest in late summer to early autumn. Weisberg et al. (1996) hypothesized that these seasonal currents are of a baroclinic nature based on an observed thermal wind shear and the seasonal reversal of the across-shelf density gradient. As a consequence of the spring transition in surface heat flux from cooling to warming, they argued that spatial differences in heating (from the coast to offshore by increasing depth and from the south to north by solar declination) forms a mid-shelf cold tongue and a seasonally maximum across-shelf density gradient that supports a southeastward current. Here we examine this locally forced, seasonal circulation hypothesis by focusing on the spring transition for 1999, a year when the Loop Current, as evidenced in relatively flat isopycnal topography at the shelf break, did not have a strong direct influence the WFS. Our objective is to describe the circulation and temperature budget for spring 1999 with respect to the shelf-wide winds, surface heat fluxes, and river inflows.
An argument used to describe the transition from wintertime horizontal stratification to summertime vertical stratification for mid-latitude shelves is that decreased winds and increased solar heating conspire to form a thermocline. The details of this process are not well understood, however. Chapman and Gawarkiewicz (1993) reason that nonlinearity in the equation of state can account for the elimination of horizontal stratification by spatially uniform heating, but Morey and O’Brien (2001) point out that this argument is valid only for certain salinity and temperature ranges. Other processes must also be important. Morey and O’Brien (2001), using a two-dimensional model with a sloping bottom, argue that the surface heat flux divided by the water depth is the critical factor in the seasonal transition, essentially the differential heating argument advanced earlier. The degree to which their argument is valid in the fully three-dimensional sense and the regional partition between ocean dynamical and local heating affects are also topics of our paper.
The observational record [e.g. Niiler, 1976; Koblinsky and Niiler, 1980; Mitchum and Sturges, 1982; Cragg et al, 1983; Marmorino, 1983; Mitchum and Clark, 1986] shows that the WFS circulation and sea level variations are highly correlated with the synoptic scale wind stress variations. The passage of cold fronts also affects the local temperature balance (e.g., Price, 1976). Along with these local synoptic scale variations are baroclinic effects that originate with the Loop Current at the shelf break [e.g., Paluszkiewicz et al, 1983]. What remains unclear are the relative importance between the momentum and buoyancy that are input either locally, or at the shelf break.
Such questions are of multi-disciplinary interest since, despite its oligotrophic description, the WFS supports highly productive ecosystems. These include episodic toxic dinoflagellate blooms (red tides) near the coast (Steidinger, 1983; Vargo et al, 1987), a seasonal chlorophyll plume near the shelf break (Gilbes et al, 1995), and important commercial and recreational fisheries throughout the WFS. Parallel programs of in situ measurements and numerical model experiments are presently in place for an improved understanding of the circulation and how it affects seasonally varying water properties and influences organism growth and distribution.
This paper focuses on local wind and buoyancy forcing during the spring transition of 1999, independent of the Loop Current. We use the primitive equation, Princeton Ocean Model (POM) described by Blumberg and Mellor (1987) forced by National Center for Environmental Predication (NCEP) reanalysis winds and net surface heat flux and by river inflows. The only role of the adjacent Gulf of Mexico is to set the vertical distribution of temperature and salinity for initializing the model density field. Once begun, the integration proceeds solely on the basis of local forcing. By running twin experiments, one with heat flux and the other without, we explore the relative importance of wind and buoyancy in effecting the seasonal and synoptic scale variability.
Section 2 describes the model and forcing fields. Section 3 compares model results with in-situ observations. Based upon these comparisons the model is used in section 4 to describe the seasonal mean circulation on the WFS for spring 1999 and the evolution of the corresponding temperature and salinity fields. Section 5 presents a term-by-term analysis of the three-dimensional temperature budget. The results are summarized and discussed in section 6.
2. Model and Forcing fields
We use the POM for the following reasons. First, it has an embedded turbulence closure submodel (Mellor and Yamada, 1974, 1982; Galperin et al., 1988) for parameterizing vertical turbulence mixing. Second, it employs a sigma coordinate in the vertical which, with the turbulence closure sub-model, is well suited to study the nonlinear dynamics over a shallow, gently sloping continental shelf. Third, its orthogonal curvilinear coordinates in the horizontal are convenient for resolving the near-shore regions.
Previous WFS POM applications include Yang et al. (1999a,b), Li and Weisberg (1999a,b), Weisberg et al. (2000), and Weisberg et al. (2001). Yang et al. (1999a) studied the WFS response to climatological monthly mean wind forcing. Qualitative agreements were found between the model and observations to some extent, but monthly mean wind stress alone could not account for the southeastward current observed at mid-shelf in spring (Weisberg et al., 1996). Li and Weisberg (1999a,b) focused on synoptic scale winds, respectively describing the kinematics and dynamics of WFS responses to idealized upwelling favorable wind forcing under constant density. The inner-shelf length scale was found to be a frictional one, consistent with the analytical work of Mitchum and Clarke (1986). The same model was also applied to a specific upwelling case study with both constant density and stratified conditions (Weisberg et al. 2000). A comparison of in-situ data with model results confirmed a simple Ekman-geostrophic route to spin-up and identified regional upwelling centers promoted by coastline and isobath geometries. The utility of this model in replicating the longer-term synoptic scale variability and the sensitivity of the WFS response to its background density state was then shown by Weisberg et al. (2001). An asymmetry in the inner-shelf responses to upwelling and downwelling favorable winds was found that helps to clarify the scale of the inner-shelf response. Consistent with Ekman dynamics, the inner shelf is the region with divergent bottom boundary layer. Since thermal wind effects can either enhance or decrease the bottom boundary layer development, the inner-shelf under stratified conditions can respond asymmetrically to upwelling and downwelling favorable winds.
Many WSF circulation questions remain. For instance, the relative effects of buoyancy and wind forcing have not been considered. While previous drifter and current meter measurements hint at seasonally varying circulation patterns (Tolbert and Salsman,1964; Gaul, 1967; Williams et al., 1967; Weisberg et al., 1996; Yang et al., 1999b), there are no definitive measures of these. The communication of water between the deep ocean and the shelf and the communication of water between different regions of the shelf remain poorly understood. To address these questions it is necessary to explore a broader region and consider more complete forcing functions than in previous studies.
The model domain (Fig 1) extends from the Florida Keys in the southeast to west of the Mississippi River in the northwest, and it has one open boundary arcing between these two locations for which a radiation boundary condition (Orlanski, 1976) is used. The model domain includes the major rivers that impact the WSF and the Desoto Canyon region where the shelf is narrowest, and its orthogonal curvilinear grid has horizontal resolution that varies from less than 2km near the coast to 6 km near the open boundary. Vertically, the sigma coordinate has 21 layers with higher resolution near the surface and bottom to better resolve the frictional boundary dynamics. In total, the model has 121´ 81´ 21 grid points. Horizontal diffusivities are parameterized using the Smagorinsky (1963) formulation with a coefficient of 0.2. Bottom stress, t b, is calculated by a quadratic law with variable drag coefficient having a minimum value of 0.0025. A mode splitting technique is used for computational efficiency (Blumberg and Mellor, 1987). Here we use external and internal time steps of 6 seconds and 360 seconds, rspectively.
The model is initialized at rest with horizontally uniform stratification. Stratification above 200m is based on temperature and salinity observations taken during a March 1999 trans-shelf hydrographic survey [from the Ecology of Harmful Algal Blooms (ECOHAB) Program]. Stratification below 200m is based on climatology. From this initial zero-baroclinicity state, the model spins up rapidly, generating baroclinicity in balance with the wind and buoyancy forcing. An alternative is to begin with a baroclinic field and allow the model currents to come into balance diagnostically with this field before proceeding with the spring simulation. The hydrographic data are not sufficient for this, however, and spurious currents due to incorrect density would corrupt the experiment. Consistent with our objective of determining the WFS responses to local, shelf-wide forcing only, our initial baroclinicity-free state is a sensible choice.
Tidal forcing is excluded in the present application since we are not considering high frequency variability. It is recognized that tidal mixing can affect the synoptic and seasonal scales when the tidal currents are large, but here the tidal currents are only a few cm s-1. Modeled and observed tidal current analyses will be reported on separately.2.2 Atmospheric forcing
Different from previous WFS model studies that considered wind forcing only, here we include both wind and thermohaline forcing. The wind and heat flux fields are from the NCEP daily reanalysis product for the period February 28, 1999 to June 1, 1999. These values, with a grid resolution of 2.5o´ 2.5o, are interpolated onto the model grid. The NCEP winds agree well with in-situ buoy winds for the spring 1999 season. Unlike the winds, however, coarse resolution renders the NCEP heat flux unrealistic because of smaller scale WFS temperature structures. We correct for this using a relaxation method (e.g., Ezer et al, 1992; Chu et al, 1999). Thus, the surface heat flux forcing is given by
where QH is the net heat flux, q obs is an interpolation of the monthly-mean satellite observed sea surface temperature, and Cp is the specific heat. The salinity flux in this study is set to be zero. The relaxation coefficient, C, or the reciprocal of the restoring time per unit area, is set at 1m/day. Such relaxation prevents deviations from observed monthly-mean SST in an attempt to force realistic baroclinic flow structures. These structures are facilitated by turbulence mixing through the coefficients KM and KH computed with the Mellor and Yamada (1982) 2.5 level turbulence closure sub-model.2.3 Lateral boundary forcing.
Gulf of Mexico Loop Current forcing is excluded in this study for two reasons. First, previous observations and model studies concluded that persistent forcing of the middle and inner-shelf by the Loop Current is minimal (Marmorino, 1982, 1983 a,b). Second, modeling the effects on the WFS of an aperiodically varying Loop Current and its associated eddies remains a great challenge (Marmorino, 1982; Cooper, 1987), presupposing that the Loop Current itself is being described properly. To better assess the role of the Loop Current as a WFS boundary condition it will be necessary to nest a regional model with a larger Gulf of Mexico/Caribbean/Atlantic Ocean model. This is beyond the scope of the present paper that focuses on local, shelf-wide forcing only. We find, however, that local forcing is capable of driving much of the observed synoptic and seasonal scale variability.
Seven major rivers are introduced into the model domain for land derived buoyancy forcing. These are the Mississippi, Mobile, Apalachicola, Suwannee, Hillsborough, Peace and Shark rivers. We use the technique of Kourafalou et al (1996)(also see Pullen 2000), whereby interpolated monthly mean mass flux data for these rivers are input to the top sigma level at the grid cells closest to the rivers’ locations.
We define the spring season here as March 1 to May 31, and we focus on this period for 1999. As an initial value problem we begin from a state of rest on February 28. With no initial baroclinicity the spin-up phase proceeds rapidly over the course of a few pendulum days, consistent with the barotropic response arguments for a gently sloping shelf of Clarke and Brink (1985). Under the conditions of surface cooling that occurs prior to the spring warming transition in mid-March, convective mixing very efficiently adjusts the initial density field on the shallow shelf. In other words, the "memory" of initial density field for this spring transition experiment is short, and sensitivity experiments that we performed using longer spin-up times showed very little difference from the present model results.
3. Model and Data Comparisons
3.1 Sea level
Since the model is forced without tides, all of the model and data comparisons are shown after low-pass filtering to exclude tidal and inertial period oscillations. Sea surface height comparisons are given in Fig. 2 at four different tide-gauge stations from Pensacola in the northwest to Naples in the southeast. Agreement is good at all of these with squared correlation coefficients exceeding 0.80. We conclude that coastal sea level for this three-month period responds primarily to local, shelf-wide forcing.
Comparisons are made between the modeled and observed velocity vector time series at the 50m, 30m, 25m, 20m, and 10m isobaths (moorings CM2, EC3, NA2, EC4, EC5, and EC6 in Fig. 1). As examples, we show the modeled and observed vector time series at the 50m, 25m, and 10m locations in Figs. 3-5, respectively. The observations are from acoustic Doppler current profilers (ADCP), and for each location we show comparisons at three different depths: near-surface, mid-water column, and near-bottom. These comparisons are quantified by a complex correlation analysis (e.g., Kundu 1976). Defining the modeled and observed velocity vectors in the Argand plane as and , respectively, the complex squared correlation coefficient is
where the overbar denotes a time average. The complex correlation has an amplitude and a phase, the amplitude being the correlation coefficient and the phase being the angle (measured counterclockwise) between the modeled and the observed currents. Like sea level, the modeled and observed currents also compare well as measured by the two sets of numbers provided with each plot. The left-hand sets are the seasonal mean east and north velocity components. The right-hand sets are the squared correlation coefficient, phase angle, and regression coefficient. At all stations and depths the squared correlation coefficients range between 0.62-0.82, and the orientations agree to within –10 to +20 degrees. More importantly, as seen directly from the time series, the model gets the sense of the velocity rotation correct in both the surface and bottom Ekman layers. A deficiency in the modeled currents, however, lies in the amplitudes. The regression coefficients show that the model underestimates the observed velocity fluctuations by between 20 to 50 percent.
Notwithstanding the amplitude disparity, the model reproduces the patterns of current variability reasonably well. The systematic underestimate of the currents may be the result of the low resolution NCEP forcing fields, and hence too much smoothing when interpolating these fields onto the model grid. Model performance also degrades between the shallowest and deepest comparison sites, i.e., the 10m isobath currents agree better than the 50m isobath currents. This is expected based on the frictional scale of the inner-shelf response to wind forcing (e.g., Weisberg et al, 2001). The 50m isobath is at the outer half of the inner-shelf so we anticipate a decreased correlation there.
The modeled and observed velocity comparisons are summarized in Fig. 6 where we show the mid-depth seasonal mean vectors and hodograph ellipses at all of the mooring locations. The mean vectors compare reasonably well, and while the ellipse semi-major axes are off by between 20-50%, the orientations and eccentricities tend to agree. On the basis of these agreements we now use the model to discuss the WFS circulation in spring 1999.
4. Mean Circulation
4.1 Flow fields
The seasonal mean circulation, obtained by averaging the model flow fields from March 1 to May 31, is presented in Fig 7. Three-dimensional flow features arise because of the WFS geometry that includes partial blocking by the Florida Keys, the coastline changes of the Big Bend, and the intrusion of Desoto Canyon. These features are depicted in horizontal velocity field maps given for sigma layers 2 (near-surface), 10 (mid-level), 20 (near-bottom), and for the depth-average.
The mid-level and depth-averaged fields are similar (and with the same scale), and they show the general nature of the 1999 spring season currents. A jet exists with axis situated between the mid-shelf and the shelf break. This jet originates along the northern coast east of the Mississippi River, and it flows along the relatively narrow Florida Panhandle shelf as a closely confined coastal feature. The coastal jet bifurcates at Cape San Blas into a mid-shelf part that heads along-isobath toward the southeast and a coastal part that hugs the Big Bend coastline. The mid-shelf part is consistent with the spring season southeastward current described by Weisberg et al. (1996). This mid-shelf current again bifurcates upon approaching the Florida Keys with a portion turning toward Florida Bay and another portion turning farther offshore. Why this occurs and how it affects the Florida Bay/Florida Keys region are evident in the near-bottom and near-surface flow fields (note the scale changes). Near the bottom we see a convergence of vectors on Florida Bay and the Florida Keys. In contrast to this, near the surface we see flow paralleling the Florida Keys. The near bottom flow upwells, feeding the near surface flow, and the second bifurcation is due to the recirculation associated with the upwelling. Li and Weisberg (1999a) discussed a similar recirculation due to the Florida Keys for the case of upwelling winds under a constant density setting. From the surface current map we can appreciate why surface drifters originating in the north do not penetrate the southeast portion of the WFS; the so-called 'Forbidden Zone' described by Yang et al. (1999b). We also note that the currents within the Big Bend and those that flow southward near-shore between Cedar Key and Sarasota are much weaker than the currents at mid-shelf. As will be shown in the next section, this is a consequence of a surface heat flux-induced cyclonic circulation that adds destructively (constructively) to the wind driven flow near-shore (offshore).
Three-dimensionality in the flow field is further evident in the vertical component of velocity shown for near-bottom, middle, and near-surface sigma levels in Fig. 8. Upwelling is prominent north of the Florida Keys with a maximum near the Dry Tortugas, as is often evident in satellite SST imagery (e.g., Weisberg et al., 2000). Upwelling is also prominent west of DeSoto Canyon, around Cape San Blas, and along the west Florida coastline. Other regions show downwelling. This vertical circulation is a function of both the bottom boundary layer and the interior. The near bottom w can be reconciled based on the near bottom horizontal velocity field (Fig. 7) and the bottom kinematic boundary condition. Mid-water column divergence then accounts for the transition of w from its near-bottom to near-surface values. Upwelling occurs along the entire near-shore domain near the bottom, whereas upwelling is more localized near the surface.
4.2 Temperature and salinity field
The modeled surface temperature and salinity fields sampled at the end of the model run on May 31 are presented in Fig. 9. Since the initial model temperature and salinity fields are horizontally uniform, this figure shows the combined effects of the momentum and buoyancy fluxes in changing the surface temperature and salinity. The two principal features are the mid-shelf cold tongue that extends southeastward from Cape San Blas and the low salinity tongue that also extends southeastward, but displaced seaward of the cold tongue. Both of these features occur annually on the WFS. The cold tongue is imposed to some extent through the surface heat flux relaxation, whereas the low salinity tongue is a fully prognostic result of the model.
The low salinity tongue derives as a river plume, accumulating fresh water primarily from the Mississippi River with additions from the Mobile and Apalachicola Rivers. Modulated interannually, the low salinity tongue extends southward each year, and in some years (1993, for instance – see Dowgiallo et al., 1995) it can be traced around the Florida peninsular to the Carolinas. Using Coastal Zone Color Scanner data between 1979 and 1986, Gilbes et al. (1995) reported a spring chlorophyll plume at mid-shelf (termed the ‘Green River’) that also extends southeast from Cape San Blas. Their explanations for the plume included: i. nutrient fluxes from the Apalachicola River; ii. nutrient fluxes from the Mississippi and Mobile Rivers; iii. seasonal changes in steric height between the shelf and deeper Gulf of Mexico; and iv. circulation of water from the Loop Current. Our results help to clarify these speculations. Consistent with Weisberg et al. (1996) the ‘Green River' is associated with both the cold tongue and the low salinity tongue. The low salinity tongue, at least in 1999, appears to originate at the Mississippi River. It is advected eastward, on average, by the spring coastal jet and then southeastward where the jet bifurcates at Cape San Blas. This bifurcation, in part, is related baroclinically to the cold tongue; hence the offset between the low salinity and cold tongues. Since the Loop Current is not included in our model, and since it did not extend to the northeast Gulf of Mexico in spring 1999, we can rule out the Loop Current as a primary conveyance of these salinity and temperature features. What remains unclear is the relative importance of land-derived versus upwelling-derived nutrients in fueling the chlorophyll plume. Aggregation of nutrients and phytoplankton by the frontal regions of the cold tongue/low salinity tongue complex can nevertheless account for the 'Green River'.
To explore the roles of momentum and heat flux in the formation of these surface temperature and salinity features we ran a model twin experiment forced by NCEP winds only. Fig 10 shows the seasonal mean circulation with wind forcing only for comparison with Fig.7 with both wind and heat flux forcing. Differences occur in all fields. The most obvious is the location of the jet, as seen in either the depth-averaged or the middle level fields. With wind forcing only the jet has a stronger coastal expression over the entire WFS, as contrasted with a stronger mid-shelf jet under wind and heat flux forcing. Baroclinicity explains the difference. With wind and heat flux forcing a dynamical feedback occurs with the cold tongue causing a cyclonic baroclinic circulation that adds constructively (destructively) with the wind forced circulation at mid-shelf (near-shore). A baroclinic circulation also arises without heat flux, but this adds constructively with the wind forced circulation near-shore. This latter effect for wind forcing only is seen in the surface temperature field of Fig. 11. Surface temperature is colder everywhere since the ocean circulation merely rearranges it. By virtue of mean upwelling the near-shore region is colder than the offshore region which causes a southward baroclinic circulation that adds constructively with the southward wind-driven circulation.
Given the importance of surface heat flux on the seasonal circulation it is instructive to see how the baroclinic contribution evolves. This is shown in Fig. 12 for the March, April, and May monthly mean, depth-averaged flows relative to the seasonal mean. March shows an anticyclonic circulation with southeastward flow near-shore. During April, once the effects of positive heat flux set in, we see the formation of the mid-shelf jet and a cessation or reversal of the near-shore flow. This further develops into a strong mid-shelf/shelf-break current in May as the cyclonic baroclinic flow becomes fully developed. This is consistent with the Weisberg et al (1996) hypothesis on the origin of the springtime southeastward current.
Temperature cross-sections at various positions along the WFS provide further information on the cold tongue evolution. These are shown for wind and heat flux forcing and for wind forcing only on March 15, April 15 and May 15, in Figs. 13, 14, and 15, respectively, at four transects offshore of DeSoto Canyon, Cape San Blas, the Big Bend, and Sarasota. Since the heat flux is initially out of the ocean the March 15 transects, with or without heat flux, are similar at depth. They differ on the shelf where surface cooling, coupled with efficient convective mixing, produces typical wintertime horizontal stratification. By April 15, with a reversal in the sign of the heat flux, the two cases depart everywhere. The Big Bend transect shows a dome of cold water at mid-shelf in the heat flux case. This originates primarily by the burial (by water being heated from above) of the cold water previously formed in the Big Bend. Water is warmer in the shallows near-shore for the same reason that they were colder in the previous season (a similar heat flux over a shallower layer will result in a larger internal energy change per unit area). There are also secondary contributions from upwelling at the shelf break. Here the shelf break upwelling is from the Big Bend as opposed to the DeSoto Canyon, and this is in part a consequence of bottom topography. The shelf break in the Big Bend occurs about 20m deeper than in the DeSoto Canyon. Hence, less upwelling is required for water to rise onto the shelf in the Big Bend than in the DeSoto Canyon. This finding is clear in continuous animations where we also see that the advective time scale between the DeSoto Canyon and the Big Bend is longer than the synoptic weather time scale. For spring 1999, at least, we do not see a direct connection between the DeSoto Canyon and Big Bend regions as would be evidenced by continuous isotherms. Note that while 210C waters occur on the shelf offshore of both the DeSoto Canyon and Big Bend regions these waters are not seen on the shelf at the intervening Cape San Blas cross section. Local upwelling in the Big Bend region (due to a relatively deep shelf break) is therefore important, and it may help to explain the relatively high productivity of the Florida Middle Grounds. A similar subsurface cold-water core is seen farther south off Sarasota and for the same reasons. It is smaller in magnitude because of differential cooling and heating from north to south, i.e., when the heat flux reversed from cooling to warming the temperatures off Sarasota were already warmer than those in the Big Bend. Dynamically, the doming of cold water at mid shelf (and hence the cold tongue) induces the cyclonic baroclinic circulation. This is further facilitated by the no heat flux bottom boundary condition that bends the isotherms perpendicular to the bottom. These combined effects of differential heating (near-shore to offshore by bottom slope and southward to northward by solar declination) and dynamical feedbacks (baroclinic pressure gradients associated with the temperature gradients) conspire to cause the spring season temperature transition from well mixed to stratified with a cold tongue at mid shelf.
The May 15 transects show the continued evolution of these processes. By now the mid-shelf isotherm doming has relaxed due to continued heating and so has the mid-shelf circulation due to decreased winds and baroclinicity. The location of the jet is now farther offshore at the shelf break, and the modeled temperature agrees with observed ECOHAB Program hydrography.
4.3 Lagrangian drifters
Model-simulated Lagrangian drifter tracks provide additional information on the evolution of the cold and low salinity tongues and the regions for possible communication between the deep ocean and the shelf. Dynamically passive particles released in the model on April 1stare tracked through May 31st. These are shown in Figs. 16 and 17 for particles originating either near the surface or the bottom, respectively. Most evident for the particles released near the surface is advection by the mid-shelf jet and the bifurcations at Cape San Blas and north of the Florida Keys. The particle trajectories agree with the Eulerian maps previously shown and with the cold and low salinity tongue features. With the exception of those in the mid-shelf jet the particle displacements are not very large. The near-surface drifters also show a propensity to flow either along-isobath or offshore. The near-bottom released drifters show subtle, but important differences, particularly in the Florida Big Bend/Middle Ground region. There we see onshore flow between the shelf break and the mid-shelf. This is consistent with the changes in hydrography discussed in the previous section, supporting our speculation that the Big Bend may be an important region for deep ocean/shelf interactions.
5. Temperature Budget
5.1 The temperature equation
To analyze the temperature budget, the temperature equation is recast from its modeled flux divergence, sigma level form to an advective, z-level form. Thus, we diagnose
a b c d vf hfx hfy
which equates the local rate of change of temperature (a) to a combination of the flow field advective rate of change (b+c+d), and the rates of change by vertical diffusion (vf), and horizontal diffusion (hfx+hfy). The temperature balance is explored using: i. time series of vertical averages, ii. time series of vertical profiles at four different locations, and iii. term by term spatial maps of the sea surface temperature rates of change.
The four analysis locations (see Fig 1) are chosen with respect to the cold tongue and the upwelling region north of the Florida Keys. Point A is on the 50m isobath at the seaward side of the Big Bend shelf subsurface cold dome. Point B is on the 40m isobath west of Tampa Bay where the cold tongue begins to taper off. Point C is on the 15m isobath offshore of Sarasota on the inshore side of the cold tongue. Point D is north of the Florida Keys. A term-by-term analysis of (3) quantifies the contributions by each physical process in changing the temperature.
5.2 Depth-averaged balances.
A depth-averaged temperature equation is obtained by vertically integrating equation (3). Since horizontal temperature diffusion is generally at least an order of magnitude less than the other terms, the depth-averaged diffusion term is essentially the depth-averaged vertical diffusion, , where Q is the net surface heat flux, r and CP are the seawater density and specific heat, and H is the water depth. The temperature variations depend on both ocean advection and diffusion. With a two-dimensional model, Morey and O’Brien (2001) proposed that temperature, in a depth-averaged sense, could be well represented without advection, i.e.,
The validity of this assertion for various locations on the WFS can be evaluated by a fully three-dimensional analysis.
Time series of the depth-averaged ocean advection and diffusion terms and their sum for locations A-D are shown in Figs. 18a-18d, respectively. While not shown in these figures, we note that the sum of these two terms is nearly exactly equal to the negative of the local rate of temperature change, which is essentially a check on our budget analysis. At station A, the depth-averaged rates of change are relatively small, and a clear distinction exists between the balances occurring on seasonal and synoptic time scales. The seasonal change (as given by the three-month mean values) is primarily by surface heat flux, whereas the synoptic variability (as given by the standard deviations) is primarily by ocean advection. Although a relatively small contribution to the seasonal change, the ocean advection does provide a cooling influence that offsets the surface heating by about 4% at this location. Thus, of the total change in vertically averaged temperature (1.660C), the contributions by ocean circulation and surface heat flux are a cooling of -0.070C and a warming of 1.730C, respectively. In contrast with the seasonal change, for synoptic scale variability we see that the ocean circulation is twice as effective as the surface heat flux in changing the vertically averaged temperature at station A.
Somewhat different results are found at station B. While surface heat flux and ocean circulation again control the seasonal and synoptic scale variability, respectively, the relative magnitudes and the signs change. Of the total change in vertically averaged temperature (3.050C), both the ocean circulation and the surface heat flux at the southern terminus of the cold tongue contribute to warming at 0.680C and 2.370C, respectively. The ocean circulation also plays a proportionately larger role in the synoptic scale variability at station B than at station A.
The relative importance of the ocean circulation versus surface heat flux continues to change approaching shallower water. At the 15m isobath, station C, of the total change in vertically averaged temperature (4.760C), the contributions by ocean circulation and surface heat flux are cooling of -1.240C and a warming of 6.000C, respectively. The shallower the water, the larger the total spring transition change. However, the role of the ocean circulation in this change varies in both magnitude and sign with location. For instance, at station D in the upwelling region north of the Florida Keys, the ocean circulation provides a cooling influence of 2.790C that partially offsets the surface heat flux warming influence of 8.110C. Without the upwelling influence by ocean circulation the water temperature north of the Florida Keys would be much warmer. Along with these increased effects by the ocean circulation on the seasonal mean temperature change, the magnitude of the controlling ocean circulation effect on the synoptic scale variability also increases with decreasing water depth. The temperature budget on all time scales is a complex one requiring a fully three-dimensional description.
5.3 Vertical profiles of the term by term balances
The temperature budget three-dimensionality is further explored through time series of the depth profiles of the individual terms that comprise the temperature balance at these four stations (Figs.19-22). In each of these figures the left hand panels show the horizontal and vertical components of the ocean advection and their sum, and the right hand panels show the diffusion, the diffusion plus the advection (which is nearly exactly equal to the negative of the local rate of change of temperature), and the temperature. With the exception of the initial portion of the record when the surface heat flux is out of the ocean and convective mixing is evident, the immediate impact of the surface heat flux is primarily in the near surface region (Fig. 19). Without convective mixing, it is then only through turbulence mixing, brought about by the ocean circulation dynamics, that the surface heat flux effect accumulates downward to make its contribution to the depth-averaged spring transition. From the vertical distributions of the seasonal means shown to the right of each panel we see that turbulence mixing distributes the surface heat flux over about 20m depth. Along with its role in turbulence mixing, the direct role of the ocean dynamics through advection on the synoptic scale variability is seen in both horizontal and vertical directions. Omission of any coordinate direction would compromise the model’s ability to describe the temperature evolution. For instance, the strong stratification at the end May at station A is due largely to the combined effects of horizontal and vertical advection. Similar conclusions follow from station B (Fig. 20). Station C (Fig. 21), in shallow water, shows the additional effect of turbulence mixing in the bottom Ekman layer. Here, temperature is elevated near the bottom by ocean circulation-induced mixing. The role of the ocean circulation in promoting mixing is seen in the covariability between the advection and diffusion terms. These interactions are even more evident at station D (Fig. 22) where warming by near-bottom turbulence mixing is required to partially offset the cooling influences by horizontal advection and upwelling. Another interesting feature of the shallow water regions is the tendency for large changes in stratification by the ocean circulation. Two downwelling events culminating in destratification and warming by advection are seen at station C around April 5th and May 10th. Similarly, but in an opposite sense, stratification caused by horizontal advection is seen at this station at the end of May.
5.4 Term by term contributions to the seasonal change in SST
The contributions to the seasonal mean rate of change of sea surface temperature (SST) by each of the advection and diffusion terms are discussed with respect to Fig. 23. The left panels show the advection terms and their sum, and the right panels show the horizontal and vertical diffusion terms, the sum of all terms. Of particular interest here is the spatial distribution of the processes controlling SST over the model domain. Horizontal diffusion is minimal everywhere except for north of the Florida Keys where it is still the smallest of the terms. Vertical diffusion controlled by the surface heat flux, warms everywhere with maximum warming tendencies stretching along mid-shelf from the Big Bend to the Florida Keys. This warming is the seasonal signal associated with spring, which displays the spatial inhomogeneities mentioned earlier. The difference between this map and the local rate of change of SST is the ocean circulation. By advecting relatively cool water southward, advection provides a cooling influence over the WFS with largest values north of the Florida Keys. With regard to total advection, Cape San Blas divides the entire domain into two sections. To the west, and extending along the north Florida coastline to the Mississippi River, advection tends to provide a warming influence, whereas to the east, along the west Florida coastline, advection tends to provide a cooling influence. Without advection, the center of the cold tongue on the WFS would be some 3-50C warmer.
Mid-latitude continental shelves undergo a spring transition as the net surface heat flux changes from cooling to warming. Using in-situ data and a numerical circulation model we investigate the circulation and temperature budget on the WFS, including the northeast Gulf of Mexico shelf from the Mississippi River to the Florida Keys, for the spring transition of 1999. The data consist of sea level from coastal stations, velocity profiles from instruments moored across the shelf between the 50m and 10m isobaths, and hydrography from ship surveys. The model is a regional adaptation of the primitive equation, POM forced by NCEP reanalysis wind stress and heat flux fields and by river inflows. Based on agreements between the modeled and observed fields we use the model to draw inferences on how the surface momentum and heat fluxes affect the seasonal and synoptic scale variability.
Spring season features of the WFS include a mid-shelf southeastward current, cold and low salinity tongues, and a high chlorophyll plume. We account for the southeastward current in 1999 by the combined responses to local, shelf-wide wind and buoyancy forcing. Wind stress drives a circulation that tends to be strongest near-shore. Heat flux provides a cyclonic contribution that adds constructively (destructively) at mid-shelf (near-shore), thus forming the observed mid-shelf jet. This heat flux-induced baroclinic circulation is related to the spring season cold tongue. By advecting Mississippi (and other) River water it forms the low salinity tongue that is displaced seaward of the cold tongue. Convergence of nutrients and associated phytoplankton growth then accounts for the high chlorophyll concentrations ('Green River') that are co-located with these surface features. These findings support the hypothesis advanced by Weisberg et al. (1996) on the origin of the southeastward current and cold tongue through differential heating from the coast to offshore (by shoaling topography) and from south to north (by solar declination). Since we arrive at these features with a model experiment that explicitly omits the Gulf of Mexico Loop Current we argue that the Loop Current is not an essential element of these spring transition features.
Through term-by-term analyses of the temperature budget we describe the evolution of WFS temperature in spring. Surface heat flux largely controls the seasonal transition, whereas ocean circulation largely controls the synoptic scale variability. These two processes are closely linked, however. Since the ocean circulation controls the turbulence mixing, the effects of the ocean circulation are of increasing importance with decreasing water depth. For instance, warming by turbulence mixing near the bottom only occurs when the bottom Ekman layer is well developed by strong currents. Thus, the water column can warm even under upwelling influence if the mixing is large enough. The water column is also found to either stratify or destratify in response to ocean circulation changes. Examples of these effects are shown. While our temperature analyses support the Morey and O'Brien (2001) one dimensional temperature balance argument advanced on the basis of a two dimensional model we demonstrate that the temperature balance is more complicated in time and space and requires fully three-dimensional thermodynamics.
Bottom topography and coastline geometry are important in generating regions of convergence and divergence and hence upwelling centers. In particular, we show that the region north of the Florida Keys has strong upwelling in spring and we speculate on the importance of the Florida Big Bend as a region for communication between the deeper Gulf of Mexico and the WFS. The shelf break there is about 20m deeper than at the DeSoto Canyon thereby requiring less upwelling for deeper waters to broach the shelf. This may be one reason why the Florida Middle Grounds are productive. The northeast Gulf of Mexico shelf from Cape San Blas to the Mississippi River also shows mean upwelling during spring 1999.
This work was supported by grants from the Office of Naval Research, grant # N00014-98-1-0158 and the National Oceanic and Atmospheric Administration, grant #NA76RG0463. We also benefited from an evolving program with previously supported from the United States Geological Survey and the Minerals Management Service. Ocean Circulation Group staff (Messrs. R. Cole, J. Donovan, C. Merz, and P. Smith) contributed to the field work and analyses, as did students W. Hemme and J. Virmani. We thank F. Muller-Karger for providing assistance with the monthly mean SST fields and Z. Li for helpful discussions on the model diagnostics. NCEP reanalysis data provided by NOAA-CIRES Climate Diagnostics Center, Boulder, CO was obtained from their web site at http://www.cdc.noaa.gov/
Blumberg, A. F., and G. L. Mellor, A description of a three-dimensional coastal ocean
circulation model, Three-Dimensional Coastal Ocean Models, Vol. 4, N. Heaps
(ed.), 208-233, AGU, Washington, D. C., 1987.
Chew, F., On the offshore circulation and a convergence mechanism in the red tide
region off the west coast of Florida. EOS Trans. AGU 36, 963-974, 1955a.
Chu, P. Edmons, N. and Fan, C. Dynamical mechanisms for the south china sea seasonal
circulation and thermohaline varibilities. J. Phys. Oceanogr., 29, 2971-2989, 1999
Chapman, D.C and G. Gawarkiewicz, On the establishment of the seasonal pycnocline in the
Middle Atlantic Bight, J. Phys. Oceanogr., 23, 2487-2492, 1993
Clark, A. J., and K. H. Brink, The response of stratified, frictional flow of shelf and slope
waters to fluctuating large-scale, low-frequency wind forcing, J. Phys. Oceanogr.,
15, 439-453, 1985.
Clarke, A. J., and S. Van Gorder, A method for estimating wind-driven frictional,
time-dependent, stratified shelf and slope water flow, J. Phys. Oceanogr., 16,
Cooper, C., A numerical modeling study of low-frequency circulation on the west
Florida shelf, Coastal Engineering, 11, 29-56, 1987.
Cragg, J., G. Mitchum, and W. Sturges, Wind-induced sea-surface slopes on the West
Florida shelf, J. Phys. Oceanogr., 13, 2201-2212, 1983.
Dowgiallo, M.J., ed., Coastal oceanographic effects of the summer 1993 Mississippi River
flooding, Special NOAA Report, March 1994, 77pp.
Ezer, T., and G. Mellor, A numerical study of the variability and separation of the Gulf
Stream, induced by surface atmopheric focing and lateral boundary flows. J.
Phys. Oceanogr. 22, 660-682, 1992
Galperin, B., L. H. Kantha, S. Hassid, and A. Rosati, A quasi-equilibrium turbulent
energy model for geophysical flows, J. Atmos. Sci., 45, 55-62,1988.
Gilbes, F., C. Thomas, J. J. Walsh, and F. E. Muller-Karger, An episodic chlorophyll
plume on the West Florida Shelf, Continental Shelf Research, 16, 1201-1224,
Kourafalou, V. H., L.-Y. Oey, J.D. Wang, and T.L. Lee, The fate of river discharge on the
continental shelf. 1: Modeling the river plume and the inner shelf coastal current. J.Geophys. Res., 101, 3415-3434, 1996
Kundu, P.K. An analysis of inertial oscillations observed near the Oregon coast. J. Phys.
Oceanogr. 6, 879-893, 1976
Li, Z., and R. H. Weisberg, West Florida continental shelf response to upwelling
Favorable wind forcing, part II: kinematic description. J.Geophys. Res., 104, 13507-13527, 1999
Li, Z., and R. H. Weisberg, West Florida continental shelf response to upwelling
favorable wind forcing, part II: dynamical analyses, J. Geophys. Res., 104, 23427-23442,1999.
Marmorino, G. O, Variability of current, temperature, and bottom pressure across the
West Florida continental shelf, winter, 1981-1982, J. Geophys. Res., 88, c7,
Marmorino, G. O, Summertime coastal currents in the Northeastern Gulf of Mexico,
J. Phys. Oceanogr., 13, 65-77, 1983b.
Mitchum, T. G., and A. J. Clarke, The frictional nearshore response to forcing by
synoptic scale winds, J. Phys. Oceanogr., 16, 934-946, 1986a.
Mitchum, T. G., and A. J. Clarke, Evaluation of frictional, wind forced long wave theory
on the West Florida shelf, J. Phys. Oceanogr., 16, 1029-1037, 1986b.
Mitchum, T. G., and W. Sturges, Wind-driven currents on the West Florida shelf,
J. Phys. Oceanogr., 12, 1310-1317, 1982.
Morey, S. L. and J. J. O’Brien, The spring transition of thermal stratification on a mid-latitude
continental shelf. Submitted manuscript, 2001
Niiler, P. P., Observations of low-frequency currents on the West Florida continental
shelf, Memoires Societe Royale des Sciences de Liege, 6, X, 331-358, 1976.
Orlanski, I., A simple boundary condition for unbounded hyperbolic flows, J. Comput.
Phys., 21, 251-269, 1976.
Pullen, J.D., Modeling studies of the coastal circulation off northern California, Ph.D.
Dissertation, Oregon State University, 2000
Sturges, W., and J. C. Evans, On the variability of the Loop Current in the Gulf of
Mexico, J. Mar. Res., 41, 639-653, 1983.
Weisberg, R. H., B. D. Black, H. Yang, Seasonal modulation of the West Florida
continental shelf circulation, Geophys. Res. Lett., 23,2247-2250, 1996a.
Weisberg, R. H., B. D. Black, J. C. Donovan and R. D. Cole, The west-central Florida
shelf hydrography and circulation study a report on data collected using a surface
moored acoustic Doppler current profiler, October 1993-January 1995,
Department of Marine Science, University of South Florida, St. Petersburg,
Florida, Data Report, January 1996, 1996b.
Weisberg, R. H., B. D. Black, Z. Li, A upwelling case study on the Florida’s west coast.
J. Geophys. Res. 105, 11459-11469, 2000
Weisberg, R. H., Z. Li, F. Muller-Karger, West Florida shelf response to local wind
forcing, April 1998. J.Geophy. Res., in press, 2001
Yang, H., R. H. Weisberg, Response of the west-Florida continental shelf to climatological
monthly mean wind forcing, J. Geophy. Res., 104, 5301-5320, 1999
Yang, H, R. H. Weisberg, P. P. Niiler, W. Sturges, W. Johnson, Lagrangian circulation
and forbidden zone on the West Florida Shelf. Continental shelf Res. 19. 1221-
List of Figures
Figure 1 The regional model grid (upper panel) and bathymetry and station locations (lower panel). Sea level comparisons are with Florida tide gauges at Pensacola, Apalachicola, St. Petersburg, and Naples. Velocity comparisons are with acoustic Doppler current profiles from instruments moored at the 50m, 30m, 25m, 20m, and 10m isobaths (1-6). Temperature is described along transects I-IV, and the temperature budget is diagnosed at Stations A, B, C, and D.
Figure 2 Comparisons between modeled (bold) and observed (thin) sea level at Pensacola, Apalachicola, St. Petersburg, and Naples as quantified by a squared correlation coefficient, along with the NCEP wind velocity sampled at station A.
Figure 3 Comparisons between modeled and observed currents at the 50m isobath (mooring CM2) sampled at depths of 3m, 25m and 40m, along with the NCEP wind velocity sampled at station A. Each vector current time series is accompanied by its seasonal mean east and north velocity components (left hand couplet), and each model/data comparison is quantified by its squared complex correlation coefficient, phase angle (or angular deviation of the model vector from the data vector measured counterclockwise), and regression coefficient (right hand triplet).
Figure 4 Comparisons between modeled and observed currents at the 25m isobath (mooring NA2) sampled at depths of 3m, 12m and 20m, along with the NCEP wind velocity sampled at station A. Quantitative comparisons are as in Figure 3.
Figure 5 Comparisons between modeled and observed currents at the 10m isobath (mooring EC5) sampled at depths of 2m, 5m and 8m, along with the NCEP wind velocity sampled at station A. Quantitative comparisons are as in Figure 3.
Figure 6 Comparisons between modeled (bold) and observed (thin) seasonal mean velocity vectors and hodograph ellipses at mid-depth for all six mooring locations on the WFS between the 50m and 10m isobaths.
Figure 7 Modeled seasonal mean velocity vectors for the depth averaged and near-surface, mid-water column, and near-bottom sigma levels, k=2, 10, and 20, respectively.
Figure 8 Modeled seasonal mean vertical velocity component fields (converted to the z-plane) sampled at the near-surface, mid-water column, and near-bottom sigma layers, k=6, 12, and 16, respectively. Bold lines indicate the upwelling. Thin lines indicate the downwelling.
Figure 9 Sea surface temperature and sea surface salinity fields at the end of spring 1999 model simulation (May 31, 1999).
Figure 10 Modeled seasonal mean velocity vectors for the depth averaged and near-surface, mid-water column, and near-bottom sigma levels, k=2, 10, and 20, respectively, from the model twin experiment forced with wind stress only.
Figure 11 Sea surface temperature at the end of spring 1999 simulation (May 31, 1999) for the model twin experiment forced by wind stress only.
Figure 12 Evolution of the monthly mean, depth averaged velocity vectors for March, April, and May relative to the spring 1999 seasonal mean.
Figure 13 Modeled temperature sections sampled on March 15 across transects originating at DeSoto Canyon, Cape San Blas, Florida Big Bend, and Sarasota. The contour interval is 10C.
Figure 14 Same as Figure 13 except sampled on April 15.
Figure 15 Same as Figure 13 except sampled on May 15.
Figure 16 Modeled near-surface Lagrangian drifter trajectories for the period April 1, 1999 to May 31, 1999. Drifters were started along nine different transects with initial positions given by solid dots.
Figure 17 Same as Figure 16 except for near-bottom Lagrangian drifter trajectories.
Figure 18a The relative contributions to the depth-averaged temperature balance by ocean circulation and diffusion at station A. Three time series are shown: the advection, the diffusion, and their sum (which equals the local rate of change of depth-averaged temperatuere). Accompanying each time series are their seasonal means and standard deviations in units of 0C day-1 as measures of the seasonal and synoptic scale variability.
Figure 18b Same as Figure 18a except for station B.
Figure 18c Same as Figure 18a except for station C.
Figure 18d Same as Figure 18a except for station D.
Figure 19 Time series of the depth profiles of the individual terms that comprise the temperature balance at station A. The left hand panels show the horizontal and vertical components of the ocean advection and their sum, and the right hand panels show the diffusion, the diffusion plus the advection, and the temperature. To the right of each panel is the seasonal mean profile. The contour interval for each of the budget terms is 0.05OC day-1, and the contour interval for temperature is 0.5OC. Shading indicates warming and clear indicates cooling.
Figure 20 Same as Figure 19 except for station B.
Figure 21 Same as Figure 19 except for station C.
Figure 22 Same as Figure 19 except for station D.
Figure 23 The contributions made to the seasonal mean rate of change of SST by each of the advection and diffusion terms. The left panels show the advection terms and their sum, and the right panels show the horizontal and vertical diffusion terms, the sum of all terms. The contour interval is 0.020C day-1. Bold lines indicate warming and thin lines indicate cooling.