Print version ISSN 1679-8759
Braz. j. oceanogr. vol.60 no.1 São Paulo Jan./Mar. 2012
Circulation and salt intrusion in the Piaçaguera Channel, Santos (SP)
Luiz Bruner de MirandaI; Eugenio Dalle Olle; Alessandro Luvizon BérgamoII; Lourval dos Santos SilvaIII; Fernando Pinheiro AnduttaIV
IInstituto Oceanográfico da Universidade de São Paulo (Praça do Oceanográfico, 191 05508-120 São Paulo, SP, Brasil) E-mail: firstname.lastname@example.org
IIUniversidade Federal do Pará (Rua Augusto Corrêa, 01, 66075-110 Belém, PA, Brasil)
IIIFaculdade de Tecnologia "Rubens Lara" (Av. Bartolomeu de Gusmão, 110, 11045-908 Santos, SP, Brasil)
IVSchool of Engineering and Physical Sciences, James Cook University (Townsville QLD 4811, Australia)
Analysis of thermohaline properties and currents sampled at an anchor station in the Piaçaguera Channel (Santos Estuary) in the austral winter was made in terms of tidal (neap and spring tidal cycles) and non-tidal conditions, with the objective to characterize the stratification, circulation and salt transport due to the fortnightly tidal modulation. Classical methods of observational data analysis of hourly and nearly synoptic observations and analytical simulations of nearly steady-state salinity and longitudinal velocity profiles were used. During the neap tidal cycle the flood (v<0) and ebb (v>0) velocities varied in the range of -0.20 m/s to 0.30 m/s associated with a small salinity variation from surface to bottom (26.4 psu to 30.7 psu). In the spring tidal cycle the velocities increased and varied in the range of -0.40 m/s to 0.45 m/s, but the salinity stratification remained almost unaltered. The steady-state salinity and velocity profiles simulated with an analytical model presented good agreement (Skill near 1.0), in comparison with the observational profiles. During the transitional fortnightly tidal modulation period there was no changes in the channel classification (type 2a - partially mixed and weakly stratified), because the potential energy rate was to low to enhance the halocline erosion. These results, associated with the high water column vertical stability (RiL >20) and the low estuarine Richardson number (RiE=1.6), lead to the conclusions: i) the driving mechanism for the estuary circulation and mixing was mainly balanced by the fresh water discharge and the tidal forcing associated with the baroclinic component of the gradient pressure force; ii) there was no changes in the thermohaline and circulation characteristics due to the forthnigtly tidal modulation; and iii) the nearly steady-state of the vertical salinity and velocity profiles were well simulated with a theoretical classical analytical model.
Descriptors: Thermohaline properties, Circulation, Theoretical-observational profiles, Salt transport, Stratification.
A análise de dados termohalinos e correntes medidos em uma estação fixa no Canal de Piaçaguera (Estuário de Santos) no inverno foi feita em termos de condições cíclicas da maré (quadratura e sizígia) e quase-estacionária, com o objetivo de caracterizar a estratificação da massa de água estuarina, sua circulação e transporte de sal forçados pela modulação quinzenal da maré. Foram utilizados métodos clássicos de análise de dados observacionais horários e quase sinóticos e de simulações analíticas de perfis estacionários de salinidade e do componente longitudinal da velocidade. Durante o ciclo de maré de quadratura as velocidades de enchente (v<0) e vazante (v>0) variaram de -0.20 m/s a 0.30 m/s, associadas à pequena variação de salinidade entre a superfície e o fundo (26.4 psu a 30.7 psu). No ciclo de sizígia a velocidade aumentou de -0.40 m/s a 0.45 m/s, mas a estratificação de salinidade permaneceu praticamente a mesma. Os perfis estacionários teóricos de salinidade e de velocidade apresentaram boa concordância (Skill próximo a 1,0) quando comparados aos perfis observacionais. Durante a modulação quinzenal da maré não houve alteração na classificação do canal estuarino (tipo 2a-parcialmente misturado e fracamente estratificado), pois a taxa de aumento da energia potencial não foi suficiente para ocasionar a erosão da haloclina. Esses resultados, associados à alta estabilidade vertical (RiL >20) e ao número de Richardson estuarino (1,6), permitem as seguintes conclusões: i) o mecanismo que forçou a circulação e os processos de mistura foi principalmente o balanço da descarga fluvial com a maré, associado ao componente baroclínico da força de gradiente de pressão; ii) não houve variações nas principais características termohalinas e da circulação devido à modulação quinzenal da maré; e iii) os perfis quase estacionários de salinidade e da velocidade foram adequadamente simulados com um modelo analítico clássico.
Descritores: Propriedades termohalinas, Circulação, Perfis teóricos-observacionais, Transporte de sal, Estratificação.
Since earlier times maritime countries used the estuarine environments to enhance urban communities and industrial development. The Santos estuary is an outstanding Brazilian example; it is bordered by the cities of Santos, São Vicente, Guarujá, Bertioga and Cubatão which is one the most industrialized of South America.
The Piaçaguera Channel is located in the upper reaches of the Santos estuary (SE) and its northern region is oriented in the N-S direction, and at this end industrial plantations have maritime terminals facilities (Fig. 1). In view of preceding studies in the Santos estuary, a need to broaden the knowledge of this channel emerged.
The hydrodynamics and mixing processes of the SE play a key role in the transport of properties concentration, pollutants and in the erosion, transport and deposition of river sediments. These processes are driven mainly by the estuarine circulation. However, due to water-bottom interactions there is a constant feedback between the hydrodynamic, mixing processes, estuarine substrate and morphology, at both the micro and macro-scale levels. It is a complex estuarine channel whose geometry and river discharge has been drastically altered during the last century by urban and industry development, land reclamation, dredging and effluent receptor from several industries.
In this paper the knowledge of the estuarine hydrodynamics related with the non-tidal and tidal longitudinal circulation, stratification and mixing characteristics will be applied. Further, with the steady-state analytical solutions of the longitudinal velocity and salinity profiles and the Stratification-circulation diagram (HANSEN; RATTRAY, 1965, 1966), it will be possible to estimate numerical values of the kinematic eddy viscosity and diffusion coefficients and the relative proportion of the advection and diffusion processes to the upward salt transport.
The objective of the study covered in this article is the analysis of specific series of observational data on the thermohaline properties and velocity to improve the knowledge on the dynamics at the upper reaches of this channel, during a winter fortnightly tidal modulation. In order to achieve such goal the study aims specifically at: i) characteristics on the tidal and non-tidal thermohaline properties, currents, stability of the water column and changes in the potential energy anomaly ii) channel classification; iii) advective salt transport and: iv) analytical simulation of velocity and salinity steady-state profiles, enabling to estimate some of its free parameters (kinematic viscosity and diffusion coefficients).
Finally, it is worth to mentioning that channels water column and bottom sediments have a high contamination level (CETESB, 1981, 2001). Hence, according to earlier investigations (OFFICER, 1977; HAMILTON; WILSON, 1980), the knowledge of estuarine circulation, transport processes and mixing are essential to the actions necessary to investigate the biological and chemical aspects of water quality, sediment erosion, transport and sedimentation as a tool for pollutant control.
In 1950's and 1960's pioneering work on the measurements of hydrographic properties, currents and solids in suspension was made by the Hydraulic Laboratory of the Polytechnic School of Engineering of the University of São Paulo (SP) reported in CTH (1959, 1966). In the 1970's measurements were carried out by consulting companies. The stratification and salt balance was used to classify the main channel of the system and to estimate the flushing time which varied from 2.1 to 5.0 days (CETESB, 1977). At the end of this decade measurements were made with the technical collaboration of the Instituto Oceanográfico da Universidade de São Paulo (IOUSP) and the Superintendência de Saneamento Básico do Estado de São Paulo (SABESP), with the objective to obtain complementary data for the outfall submarine project planned for the Santos Bay (Fig. 1). Studies focusing on sediments, ecological and chemical aspects of the Santos estuary were published by Fúlfaro and Ponçano (1976). Heavy metal concentration in the Santos Bay was initially investigated in 1983 (FÚLFARO et al., 1983) and further by Figueira et al. (2004). Updated depositional history studies and isotopic composition of the organic matter in the SE and the Santos Bay may be found in recent publications.
A pioneering study on the numerical modeling of the SE and the inner continental shelf forced by the main tidal components was published by Harari and Camargo (1998). Although taking into account only barotropic conditions it is possible to find out in this work some interesting features: stronger ebbing circulation east side of the bay and convergent circulation due to flooding currents in the channels of São Vicente and Bertioga; this work was further improved by Harari et al. (2001) and Harari et al. (2002) with the theoretical simulations of the barotropic circulation in the Piaçaguera Channel validated with observational results. A description of theoretical and applied studies related to the numerical modeling of the effluent dispersion using the SE as the region to be investigated was presented by Diniz and Meyer (2004).
Variations on tidal and non-tidal circulation, thermohaline characteristics, salt transport, classification and tidal effects on the nutrients, phytoplankton concentration, suspended matter and chlorophyll-a, due to the tidal modulation during consecutive neap-spring tidal cycles in the SE channels (São Vicente, Santos and Bertioga) were investigated by Miranda et al. (1998), Gianesella et al. (2000), Bernardes and Miranda (2001) and Moser et al. (2005).
Data sampling and methodology
Time series of hourly profiles of hydrographic properties (S,T,p) and current velocity (speed and direction) were recorded at an anchor station in the Piaçaguera Channel (Station: lat. 23'º53,987'S; long. 046'º22,640'W), at a depth closely to 11 m (Fig. 1). This station is located south of the maritime terminals and near the Cubatão River outflow (Fig. 1). According to the field experiment report the oceanographic sampling was made with a CTD/Current meter Falmouth, model 2D-ACM. The station was monitored during 19 days (460 h), from June 23 to July 11, 2001, covering 36 semidiurnal tidal cycles. However, due to weather conditions, maritime traffic and equipment failure the time series presented discontinuities which ranged from one hour up to several hours. Despite such discontinuities time series of salinity, density and longitudinal velocity component were analysed. From this time series a set of Eulerian hourly profiles of hydrographic properties and velocity of three complete semidiurnal tidal cycles (37 hours), for the neap (June 28-29, 2001) and spring (July 04-06, 2001), were used to compute the non-tidal nearly steady-state salinity and the v-velocity component profiles.
The temperature (T in oC) and salinity (S in psu) profiles were analysed taking into account its tidal and non-tidal characteristics. The density was calculated with the International Equation of State of Sea Water. The velocity vector was decomposed in longitudinal (N-S or v-component) and secondary (E-W or u-component) components according to the local referential system Oxyz (Ox, Oy and Oz with the origin on the free surface and oriented positively toward East, South and upward, respectively); ebbing and flooding currents are indicated by v>0 and v<0, respectively. All raw data were submitted to quality control and then interpolated in terms of the non-dimensional depth (Z=z/|h(t)|) at intervals -0.1 from the surface (Z=0) down to the bottom (Z=-1.0) and at hourly time intervals.
The description of the steady-state dynamic response and salt stratification due to the combination of the river discharge, density-driven forcing and wind was approximated by the two-dimensional Hansen and Rattray (1965) steady-state analitical solution, further extended by Officer (1976, 1977), Prandle (1985, 2004) and others, to include bottom frictional effects. The Hansen and Rattray's classical model consider the balance of the barotropic and baroclinic modes, with wind stress forcing on the surface and no-slip at the bottom as upper and lower boundary conditions, respectively. An idealized rectilinear coastal-plain estuary (B and h being the width and depth, respectively) and laterally homogeneous is assumed in the model. Bottom friction due to the tidal currents is the predominant cause of turbulent mixing but have no other influence upon the net circulation on the estuary.
The central regime solutions, satisfying boundary conditions and assuming constant kinematic eddy viscosity and diffusion coefficients, for steady-state v-velocity (vc) and salinity (Sc) profiles in terms of the non-dimensional depth (Z), take the following forms obtained from the derivation and integration of the current function (HANSEN; RATTRAY ,1965);
The v-component profile of the velocity vc=vc(y,Z) (eq. 1) is expressed as the sum of three circulation modes: i) the gravitational circulation in which the pressure gradient due to the longitudinal salinity (density) gradient (temporally and vertically constant) is balanced by frictional forces, which is associated with an analog of the Rayleigh number (Ra) times the n parameter defined as the ratio of the tidal diffusion term to the total upstream net salt transport due to the advection and tidal diffusion (SIMPSON et al., 2001); ii) the velocity generated by the river discharge (vf=Qf/Bh); and iii) the contribution of the surface wind-stress forcing, tW. The quantity Nz is the vertical kinematic coefficient of eddy viscosity treated as a constant. In the theory the density r is approximated by the linear equation of state of sea water. It should be pointed out that the analytical expression of the equation (1), is equivalent to the one derived latter by Officer (1976).
In the salinity profile Sc=Sc(y,Z) (eq. 2) the quantity So is the time-mean salinity averaged along the water depth. The quantities x and M are the dimensionless longitudinal coordinate and the ratio of the tidal mixing to the river flow given by x=vfy/Kho and M=KzKho/(vfh)2, where Kz and Kho are the vertical and horizontal kinematic coefficient of eddy salt diffusion considered as constants, respectively. Taking into account the hypothesis of the Hansen and Rattray's theory the eddy coefficients (Nz, Kz and Kho) and the wind stress will be treated as free parameters, to be used to validate the theoretical profiles against near-steady state profiles derived from experimental data. In order to quantify the agreement between the theoretical results of the velocity and salinity profiles and the in situ data, the method suggested by Wilmott (1981) based in the Skill parameter defined by Warner et al. (2005) was used,
The quantities Xobs and Xmodel are the observational properties (v-velocity component or salinity) and the corresponding values calculated by the model (vc and Sc) and (Xm)obs its depth mean observational values. According to the equation (3) the Skill parameter varies from 1 to zero indicating the best fit and a complete disagreement between observation and the theoretical result, respectively. The validation skill parameter as defined (eq. 3) was used to evaluated numerical simulations of the Curimataú River estuary (RN-Brasil) with a three dimensional model (Delft 3D- Flow) in comparison with set of time series measurements (ANDUTTA et al., 2006).
The classical Stratification-circulation diagram theoretically derived by (HANSEN; RATTRAY, 1966) was used in the estuary classification and to calculate the relative contribution of the advective and diffusion processes to the upestuary salt transport through its key parameter (n). The nearly-steady salinity [S=S(x,z)=S(x,Z)] and the longitudinal velocity component [v=v(x,z)=v(x,Z)] used to calculate the stratification (pe)circulation parameters (pc) were obtained from the hourly timemean values of vertical profiles based on the algorithms of Miranda et al. (2002) and Bérgamo et al. (2002). These parameters were also calculated from the theoretical profiles simulations of the v-velocity component and salinity. Taking into account the residual values of velocity and salinity the advective salt transport was calculated with the assumption that the channel was laterally homogeneous.
Stability characteristics during the tidal cycle and the main forcing mechanisms for the circulation and mixing processes were analyzed by the classical Richardson layer and estuarine numbers, RiL and RiE, respectively (FISCHER, 1976; BOWDEN, 1978):
in the first expression g, h(t), ΔρV(t), rm(t) and (vm)2 are the gravity acceleration, the depth, the bottom minus the surface densities, the mean-depth density, and the square of the mean-depth velocity. In the RiE definition ΔρH, hm, vf=va, rm and vrms, are the differences of the densities at the mouth and the estuary head, va is the residual velocity and (vrms)3 is the root mean square cubed of the v-velocity component which controls the intensity of tidal stirring; due to the lack of the fresh water discharges the velocity (vf) was approximate by the residual velocity va.
The RiL=20 is the upper limit for which the turbulent mixing occurs near the halocline in partially mixed estuaries. Bellow this critical number (RiL<20) the bottom turbulence became effective to the vertical mixing process in the water column. Bellow RiL=2 the turbulent mixing is isotropic fully developed and the water column is unstable and for RiL >20 the water column is stable with low vertical mixing (DYER; NEW, 1986).
The energy of the tidal Straining Process due to the Periodic Stratification (SIPS) due to the fortnightly tidal modulation and its influence on water column stability was calculated in terms with the potential energy anomaly f (NUNES et al., 1989; SIMPSON et al. (1990)
where ρ=r(z) and rm are the density profile and its mean-depth, respectively. f is the work per volume unit required to mix of the water column (J/m3, in the SI system); f=0 for a well mixed estuary and increases for a high stratified estuary. Taking into account that for a non-compressible horizontal flow the local density variation is related to the advection term by , and it is readily shown that the time rate of is given by (SIMPSON et al., op. cit.).
where v=v(z) and vm are the v-velocity component profile and its mean depth, respectively. In terms of the non-dimensional depth (Z) this equation may be written as
and is to be expected to have low and high values for a well mixed and a highly stratified estuary, respectively. The time rate of f will be used to compare the tidal stirring intensity changes over the fortnightly tidal modulation of the Piaçaguera Channel based on the observational results of v=v(y,z)=v(y,Z).
Understanding the development and breakdown of stratification is a key to shallow sea and estuarine oceanography. The stratification intensity in the water column is crucial in controlling the intensity of vertical mixing and hence the vertical fluxes of water properties such as heat, salt, momentum, and nutrient concentrations. This intensity may be of critical importance in controlling the biological production because inhibiting vertical displacement might influence the degree of light exposure experienced by marine organisms (SIMPSON et al., 1990).
The forcing mechanisms
The fresh water input
Fresh water is discharged into the SE from a complex drainage basin with several rivers with springs in the highlands maritime border (Serra do Mar). Fresh water feeding the Piaçaguera Channel are mainly from the following rivers: Piaçaguera, Cubatão, Quilombo and Mogi. The processes of erosion and transport carry high concentration of suspended matter which are deposited in the mangrove forest and the in the estuary bed.
The complexity of this drainage basin and the lack of the gauging stations of river discharges make the calculation of fresh water input to the channel uncertain. For the whole estuary the historical data published by CETESB (1977) gives an annual mean estimate of 42.8 m3/s which is close agreement with the one based on the climatological values of air temperature, precipitation and evapotranspiration rates (41.0 m3/s) unpublished data calculated by Dr. Mario Pereira da Silva (personal communication). For the Cubatão River there was a gauging station (3E-077 - Ponte Preta - lat. 23o53'39"S; long. 046o27'26"W) with only a few years of measurements which gives an annual climatological mean of 7.5 m3/s.
Hourly tidal heights at the station of Conceiçãozinha (lat. 23o58,1'S; long. 046o17,5'W) were used for the tidal analysis and prediction of tidal spectrum obtained according to Pawlowicz et al. (2002). These results agree very well with the previous ones (FRANCO, 2000) indicating that the estuary is microtidal (ranges less than 2 m) and semidiurnal mixed (form number 0.32). Prominent peaks occur at diurnal (K1, O1), semi-diurnal (M2, S2) and quarter-diurnal (M4). Smaller peaks also appear at the M6 tidal harmonic.
Tidal oscillations from June 1st to July 31, 2001 corresponding in the relative time scale from 0 h to 1.400 h, and the experimental period neap-spring tidal cycle is shown by the horizontal bar (Fig. 2).
During the neap tide observation period northerly winds were strong (speed up to 60 km/h) and almost damped the tidal oscillations and a quick decrease in the tidal oscillations of 0.6 m was observed (Fig. 2). By this time the measurements had to be interrupted for several hours for safety reasons. In the beginning of the subsequent spring tide it was observed a small increase in the amplitude of the tidal oscillations and the mean sea level remained low and superimposed to the tidal oscillations with amplitudes of almost 1.5 m.
Experimental and theoretical results
Neap and spring tidal experiments
The Eulerinan vertical salinity profiles variations [S=S(Z,t)] generated by the influence of tidal diffusion, tidal currents and the river discharge in the neap-spring tidal experiments (28-29 and 04-06 of June and July, 2001, respectively) covering three semidiurnal tidal cycles are presented in Fig. 3. These profiles show the advective influence of the ebbing and flooding tidal currents, the halocline intensification and erosion due to vertical mixing intensification during the tidal cycle and smaller differences between neap and spring tidal cycles. The surface and bottom salinities (in psu) varies in the intervals: 25.6<S<29.7 and 30.2 <S<30.7 (Fig. 3-A) for the neap and are slightly smaller for the spring tides (25.9 <S<29.2 and 29.5 <S<30.4) (Fig. 3-B).
Due to the typically estuary small temperature variations the thermohaline characteristics of the neap-spring tidal cycles are shown by the scatter T-S Diagram (Fig. 4). The images of the T-S plots show the estuarine water mass characterized by a small temperature variation (20.3oC to 22.0oC) associated with salinities ranging from 25.6 psu to 30.7 psu. There is a distinctive T-S scatter for the neap (o) and spring (x) tidal cycles with slightly higher temperatures in the last one.
The time-depth mean profiles of salinity [S=<S(Z)>] at the neap and spring tidal cycles reflect the small differences in the Eulerian profiles despite the fortnightly tidal modulation and the low vertical salinity stratification (bottom minus surface salinities equal to 2.5) with very close time-depth mean values (29.8 psu and 29.2 psu), respectively. These profiles will be shown later in comparison to the steady-state simulations.
The Eulerian profiles of the v-velocity component [v=v(Z,t)] associated with the salinity profiles (Fig. 3) showed distinctive patterns with asymmetric variations from -0.18 m/s to 0.30 m/s (Fig. 5-A) and -0.39 m/s to 0.46 m/s (Fig. 5-B) for the neap-spring tidal conditions, respectively (Fig. 5-A,B).
During the neap-spring tidal experiment there are profiles showing ebbing (v>0) and flooding (v<0) unidirectional motions and some exhibiting the classical bidirectional circulation with depths of no-motion which characterizes the compensating flooding currents due to the up-estuary baroclinic forcing (Fig. 5). The non-tidal velocity profiles indicate the two-layer estuarine circulation down (v>0) and up-estuary (v<0) in the upper and lower layers, respectively, with no-motion depths at Z=-0.4 and Z=-0.5 which is the characteristic bidirectional motion of the partially mixed estuaries due to the fresh water discharge, baroclinic forcing and tidal stirring. The nearly steady-state residual current (va) are down-estuary with speeds of 0.009 m/s and 0.003 m/s for the neap and spring tidal cycles, respectively. These profiles will be shown comparatively to the nearly steady-state analytical simulations.
The down (>0) and up-estuary (<0) salt transport intensities due to the residual velocity component and the baroclinic forcing, under the assumption that the channel is laterally homogeneous, decreased from 3.0 kg/s.m to 0.85 kg/s.m and -0.32 kg/s.m to -0.11 kg/s.m for the neap and spring tidal cycles, respectively.
The water column stability variation and the ratio of the buoyancy input due to the fresh water discharge and the tidal stirring, during the neap and spring tidal cycle was investigate by the Richardson layer and estuarine numbers, RiL, and RiE (eq. 4). In the neap and spring tidal cycles the water column tends to be stable (RiL>20, but there was some values with RiL close or lower than 20, which indicated a possible occurrence of turbulent mixing across the halocline. For spring tidal cycle the calculated Richardson Estuarine number is 1.6 and the fresh water discharge and tidal currents are the main forcing mechanisms.
Time variations of the longitudinal velocity component (v) during the neap and spring tidal conditions on the surface (Z=0.0), middle depth (Z=-0.5) and near the bottom (Z=-0.9) are shown in Fig. 6. At Z=0.0 and in the neap tide, the time variation of the current has an asymmetric behavior of the flooding (v<0) and ebbing (v>0) currents (with speeds from -0.06 m/s to 0.27 m/s) not well correlated with the time variations in deeper layers, and at the end of the sampling period the circulation in this layer was towards north (v<0) in opposition with the near surface circulation (Fig. 6A). However, the currents at the spring tide are in phase at all depths and well correlated with the tidal oscillation and we may observe an increase in the velocity shear towards the bottom (Fig. 6-B).
The current roses (not shown) were draw for the neap and spring tidal cycles and its analysis for the depths Z=0.0, Z=-0.5 and Z=-0.9 indicated that the secondary currents (u-velocity component) were almost absent in the surface layer and presented a small increase up to 0,01 m/s towards the bottom.
Steady-state analytical simulations
In the Hansen and Rattray's analytical steady-state solutions of the v-velocity component (vc) and salinity (Sc) (eqs. 1 and 2), in addition to the free parameters (Nz, Kz, Kho and ), we took into account the following physical quantities obtained from the experimental results: i) the mean time-depth or residual velocity (va), replacing the velocity driven by the river discharge vf; ii) the mean water depth (h); iii) the n parameter, obtained from the Stratificationcirculation Diagram (results presented in item bellow), and; iv) the distance of the anchor station from the mouth (x) measured on the nautical chart. The remaining parameters: the time-depth mean salinity Sa and the longitudinal density gradient were adjusted to obtain the best theoretical versus observational profiles fit as indicate by the Skill parameter (eq. 3). All numeric values used in the theoretical simulation are presented in Table 1 and we can notice that in the spring tide simulations the kinematics eddy coefficients Nz (Kz) are two times higher (ten times lower) than the one for the neap tide simulations, however, for the horizontal eddy diffusion coefficient Kho its numerical values were the same for both tidal cycles.
The best fits of the theoretical v-velocity component (vc) and salinity profiles (Sc) in comparison with the corresponding observational profiles are shown in Figures 7 and 8 for neap and spring tidal cycles, respectively.
Although some adjustments had to be made in the free parameters, small discrepancies were observed in the salinity profiles (Fig. 7-A) and the Skill parameter was near one (0.96). Theoretical salinity values minus the observational ones for the neap tidal cycle varies 0.58 psu at the surface, less than 0.05 psu along the halocline and 0.30 psu at the bottom. For the v-velocity profile a good correspondence (less than 0.008 m/s) between the theoretical results and the observations was found (Fig. 7-B), and the Skill=1.0 indicates a very good fit. Inspection of the experimental and theoretical vertical profiles for the spring tidal cycle (Fig. 8-AB) shows the model also gives a good agreement for the steady-state v-velocity component (Skill=0.98). However, although the salinity profiles have similar configuration, for this tidal regime the theoretical salinity values minus the observational ones varies from 0.07 psu and 0.13 psu on the surface and bottom, respectively, with higher deviation along the halocline: from 0.20 up to 0.40 psu, and the validation parameter is close to 0.97.
According to the parameters used in the simulations (Table 1), it is worth to mentioning that: i) it was assumed a weak and a negligible wind stress at the neap and spring tidal cycles, respectively, and ii) as indicated above, some adjustments were necessary in the free parameters (Nz, Kz and Kho).
Channel Classification and tidal straining
The images of the stratification [pe=(Sf-Ss)/So] and circulation [pc=(vs/va)] parameters of the Hansen & Rattray's classification diagram obtained from data of Figs. 7 and 8 for the neap (pe= 0.073 and pc=11.42) and spring (pe=0.066 and pc=12.49) tidal experiments with n=0.85, classify the Piaçaguera Channel as type 2a (partially mixed and weakly stratified) for both neap-spring tidal conditions (Fig. 9). Hence, 85% and 15% of up-estuary salt transport are due to tidal diffusion and the baroclinic shear, respectively.
The Stratification-circulation diagram was also used to classify the channel using the results of the steady-state analytical simulations of the vertical profiles of salinity and the v-velocity component, respectively, and the results were very close (Fig. 9), with the n parameter almost the same (0.92). This result is independent on the tidal fortnightly modulation and may allow the conclusion that the tidal energy reaching the channel was insensitive to fortnightly time scale tidal modulation. These results are very different when compared to the ones from the lower reaches of the estuarine Bertioga Channel (the NE entrance of the Santos estuary), which was strongly influenced by the fortnightly tidal modulation and its classification changed from type 2b (n=0.7) to type 2a n=0.9) during the transition from the neap to spring tidal cycles based on winter time observations by Miranda et al. (1998).
The energy of the tidal straining process due to the fortnightly tidal modulation and its influence on water column stability was calculated in terms of the time rate of the potential energy anomaly , according to Simpson et al. (1990). Such calculations used the steady-state results of the velocity profiles v=v(Z) and its depth-mean value (vm). The N-S density gradient used in the equation (7) was the same of the analytical simulations (1.2x10-3 kg/m4). Our results indicate that there was a decrease in the time rate of the potential energy anomaly from 1.50x10-3 J/m3.s to 0.56x10-3 J/m3.s, during the neap to spring tidal cycle transitional period, respectively. This decrease explain why the occurrence of the halocline break-down wasn't observed and the estuarine channel stratification remained partially mixed (Type 2a) in the spring tidal cycle.
Almost the same result may also be obtained using the correlation of ub (near bottom velocity) and the estuary depth (h) with the longitudinal density gradient used as a parameter, according to the Nunes et al. (1989). This correlation represents the balance of the potential energy due to the estuarine circulation and tidal stirring which is taken as the boundary between mixed and stratified regimes. By plotting the spring tide results (near bottom velocity at Z=0.9 as ub=0.1 m/s from Fig. 6) and h=11 m (Table 1) in the Figure 3 of Nunes et al. (op. cit), the image of these values in the ub versus h correlation lies well above the density gradient line (1.2x10-3 kg/m4) of the analytical simulations, which confirms that the production of turbulent kinetic energy was't enough to promote the vertical mixing of the water column in the spring tide.
During the fortnightly tidal modulation observation period the Piaçaguera channel remained as a partially mixed estuary and weakly stratified (type 2a) with almost de same mixing parameter calculated with experimental (n=0.85) and theoretical data n=0.92). This independence on the tidal forcing may be justified by the decrease of the potential energy rate which ranged from 1.5x10-3 J/m3.s to 0.60x10-3 J/m3.s over the neap to spring tidal cycle, respectively.
Although the channel was type 2a, the time variation of the layer Richardson number during the neap tide was much higher than the upper limit (RiL>20) for which the turbulent mixing occurs in the halocline. However, in the spring tidal cycle were calculated values lower than the critical number (RiL<20) and the bottom turbulence became effective. Due to the low value of the potential energy rate the vertical mixing process was inhibited in the spring tidal cycle.
Time variations of the longitudinal velocity component (v) during the neap tidal conditions at the surface (Z=0) were not well correlated with the underneath layers (Z=-0.5 and Z=-0.9) variations due to a cold front passage generating winds with speeds up to 60 km/h. However, in the spring-tide a well correlation of the v-velocity component in the water column and the tidal oscillation was observed and they are almost in phase.
The classical momentum balance was appropriate to describe the essential steady-state dynamics of the estuarine channel which was mainly due to the fresh water discharge, the baroclinic gradient pressure and the dissipative frictional forces. During those experiments the baroclinic forcing towards the estuary-head reaches the northern part of the estuarine channel located at almost 17 km from its mouth. However, due to its low dynamics the Piaçaguera Channel was almost undisturbed by the fortnightly tidal modulation during this winter time observation period and may contribute to the sediment trapping in the inner estuary region.
The decrease in the salt transport components down and up-estuary due to the residual currents and the baroclinic forcing in the transitional period form the neap and spring tidal conditions indicates the tendency towards a balance in the spring tide.
We are grateful to the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) for the Research Productivity Fellowship (Proc.302069/2004-6) and Ph.D. and M.Sc. fellowships concession to the authors and the Fundação de Estudos e Pesquisa Aquáticas (FUNDESPA). The collaboration of the graduate students of the Department of Oceanography, Chemistry and Geology of the University of São Paulo for their assistance with the field work is also very much appreciated, and to José Roberto Bairão Leite for processing the spectral tidal analysis. Minor editing was made by Patricia Leite de Miranda. We also thank an anonymous reviewer whose comments helped the article improvement.
ANDUTTA, F. P; MIRANDA, L. B.; CASTRO, B. M.; FONTES, R. F. C. Numerical simulation of the hydrodynamic in the Curimataú Estuary, RN Brazil.In: SIMPÓSIO BRASILEIRO DE OCEANOGRAFIA, 3. Oceanography and Global Changes, São Paulo, 2006. p. 545-558. [ Links ]
BÉRGAMO, A. L.; MIRANDA, L. B. DE; CORRÊA, M. A. Estuário: Programas para processamento e análise de dados hidrográficos e correntográficos. Rel. Téc. Instituto Oceanográfico, São Paulo , n. 49, p. 1-16, 2002. [ Links ]
BERNARDES, M. E. C.; MIRANDA, L. B. Circulação estacionária e estratificação de Sal em canais estuarinos: Simulação com modelos analíticos.
Rev. bras. Oceanogr, v. 49, n. 1/2, p. 115-132, 2001.
BOWDEN, K. F. Mixing Processes in Estuaries. In: KJERFVE, B., (Ed.). Estuarine Transport Processes, Belle W. Columbia, N.C.: Baruch Library In Marine Science, University of South Carolina Press, 1978. p. 11-16 (n. 7). [ Links ]
CETESB. Estudo da Circulação das Águas no Estuário de Santos. São Paulo: CETESB, 1977. 69 p. (Relatório Técnico). [ Links ]
CETESB. Metais pesados na Baía de Santos e Estuários de Santos e São Vicente. São Paulo: CETESB, 1981. 231 p. (Relatório Técnico). [ Links ]
CETESB. Contaminação do Sistema Estuário de Santos/São Vicente (SP) - Brasil. São Paulo: CETESB, 2001. 81 p. + anexos e mapas. (Relatório Técnico). [ Links ]
CTH. Baía de Santos. Levantamento de correntes na barra do porto para a Companhia Docas de Santos. São Paulo: USP-Escola Politécnica, Laboratório de Hidráulica-DAEE, 1959. 12 p. +mapas. [ Links ]
CTH. Estuário Santista. Pesquisas sobre assoreamento verificado na faixa portuária para a Companhia Docas de Santos. São Paulo: USP-Escola Politécnica, Laboratório de Hidráulica-DAEE. 1966. v.1, 112 p. [ Links ]
DINIZ, G. L.; MEYER, J. F. C. A. Estudo e modelagem matemática da dispersão de poluentes em sistemas estuarinos. Biomatemática, v. 14, p. 103-118, 2004. [ Links ]
DYER, K. R.; A. L. NEW, Intermittency in estuarine mixing. In: WOLFE, D. A. (Ed.). Estuarine Variability. Orlando, Fla.: Academic Press, 1986. p. 321-339. [ Links ]
FIGUEIRA, R. C.; TESSLER, M. G.; MAHIQUES, M. M. DE; FUKUMOTO, M. M.; CAZZOLI, S.; GOYA, Y. Sedimentação recente e níveis de metais pesados em testemunhos do Estuário Santista, SP - Brasil. CONGRESSO BRASILEIRO DE PESQUISAS AMBIENTAIS E DA SAÚDE, CBPAS, 4. Santos, SP, 2004. Anais...p. 73-76. (CD Rom). [ Links ]
FISCHER, H. B. Mixing and dispersion in Estuaries. Ann. Rev. Fluid Mech., v. 8, p. 107-133, 1976. [ Links ]
FRANCO, A. S. MARÉS: Programa para Previsão e Análise. In: São Paulo: Manual, BSP, 2000. 36 p. [ Links ]
FÚLFARO, V. J.; PONÇANO, W. L. Sedimentação atual do estuário e Baía de Santos: um modelo geológico aplicado a projetos de expansão da zona portuária. In: CONGRESSO BRASILEIRO DE GEOLOGIA E ENGENHARIA, 1, 1976. Rio de Janeiro. Anais... Rio de Janeiro: ABGE, v. 2, p. 67-90, 1976. [ Links ]
FÚLFARO, V. J.; REQUEJO, C. S.; LANDIM, P. M. B.; FÚLFARO, R. Distribuição de elementos metálicos nos sedimentos da baía de Santos. In: SIMPÓSIO REGIONAL DE GEOLOGIA, SBG, 4. São Paulo, 1983. p. 275-289. [ Links ]
GIANESELLA, S. M. F.; SALDANHA-CORRÊA, F. M. P.; TEIXEIRA, C. Tidal effects on nutrients and phytoplankton distribution in the Bertioga Channel, São Paulo, Brazil. Aquat. Ecosys. Health Manag, v. 3, p. 533-544, 2000. [ Links ]
HAMILTON, A .D.; WILSON, R. E. Nontidal circulation and mixing processes in the lower Potomac Estuary. Estuaries, v. 3, n. 1, p.11-19, 1980. [ Links ]
HANSEN, D. V.; RATTRAY, M., JR. Gravitational circulation in straits and estuaries. J. mar. Res., v. 23, n. 1, p. 102-122, 1965. [ Links ]
HANSEN, D. V.; RATTRAY M., JR. New dimensions in estuary classification. Liminol. Oceanogr., v. 11, n. p .319-325, 1966. [ Links ]
HARARI, J.; CAMARGO, R. Modelagem numérica da região costeira de Santos (SP): Circulação de Maré. Rev. bras. Oceanogr., S Paulo, v. 46, n. 2, p. 135-156, 1998. [ Links ]
HARARI, J.; CAMARGO; R.; MIRANDA, L. B. DE. Modelagem numérica hidrodinâmica tridimensional da região costeira e estuarina de São Vicente e Santos (SP). Pesquisa Naval, RJ, Suplemento, p. 79-97, 2001. [ Links ]
HARARI, J; MIRANDA, L. B.; BÉRGAMO, A. L.; CORRÊA, M. A. Hidrografia, correntometria e modelagem numérica em apoio a obras de engenharia costeira no canal de Piaçaguera (Santos-SP). In: CONGRESSO BRASILEIRO DE PESQUISAS AMBIENTAIS E DA SAÚDE, CBPAS, 2002, Santos, SP. p. 1-6. (CD Rom). [ Links ]
MIRANDA, L. B. DE. ; CASTRO, B. M.; KJERFVE, B. Circulation and mixing in the Bertioga channel (SP, Brazil) due to tidal forcing. Estuaries, v. 21, n. 2, p. 204-214, 1998. [ Links ]
MIRANDA, L. B. DE; CASTRO, B. M.; KJERFVE, B. Princípios de Oceanografia Física de Estuários. S. Paulo: EDUSP, 2002. 424 p. (Acadêmica 42). [ Links ]
MOSER, G. A. O.; GIANESELLA, S. M. F.; BARRERA-ALBA, J. J.; BÉRGAMO, A. L.; SALDANHA-CORRÊA, F. M. P.; MIRANDA, L. B. DE; HARARI. J. Instantaneous transport of salt, nutrients, suspended matter and chlorophyll-a in the tropical estuarine system of Santos. Braz. J. Oceanogr., v. 53, p. 115-127, 2005. [ Links ]
NUNES, R. A.; LENNON, G. W.; SAMARASINGHE, J. R. S. The negative role of turbulence in Estuarine mass Transport. Estuar. coast. Shelf Sci., v. 28, p. 361-377, 1989. [ Links ]
OFFICER, C. B. Physical Oceanography of Estuaries and Associated Coastal Waters. New York: John Wiley & Sons, 1976. 465 p. [ Links ]
OFFICER, C. B. Longitudinal Circulation and Mixing Relations in Estuaries. Estuaries, Geophysics, and the Environment. Studies in Geophysics. The National Research Council, 1977. p. 13-21. [ Links ]
PAWLOWICZ, R.; BEARDSLEY, B.; LENTZ, S. Classical tidal harmonic analysis including error estimates in MATLAB using T_TIDE. Computers, Geosciences, v. 28, p. 929-937, 2002. [ Links ]
PRANDLE, D. On salinity regimes and the vertical structure of residual flows in narrow tidal estuaries. Estuar. coast. Shelf Sci., v. 20, p. 615-635, 1985. [ Links ]
PRANDLE, D. Saline intrusion in partially mixed estuaries. Estuar. coast. Shelf Sci., v. 59, p. 385-397, 2004. [ Links ]
SIMPSON, J. H.; VENNELL, R.; SOUZA, A. J. The salt fluxes in a tidally-energetic estuary. Estuar.coast. Shelf Sci., v. v. 52, p. 131-142, 2001. [ Links ]
SIMPSON, J. H.; BROWN, J.; MATTEWS, J.; ALLEN, G. Tidal straining, density currents, and stirring in the control of Eeuarine stratification. Estuaries, v. 13, n. 2, p.125-132, 1990. [ Links ]
WARNER, J. C.; GEYER, W. R.; LERCZAK, J. A. Numerical modeling of an estuary: A comprehensive skill assessment. J. Geophys. Res., v. 110(CO5001), p. 1-13. 2005. [ Links ]
WILMOTT, C. J. On the validation models. Phys. Geogr., n. 2, p. 184-194, 1981. [ Links ]
Manuscript received 16 June 2010;
revised 23 October 2011;
accepted 20 November 2011