Modeling the Paranagua Estuarine Complex, Brazil: tidal circulation and cotidal charts

RESEARCH ARTICLES

Modeling the Paranagua Estuarine Complex, Brazil: tidal circulation and cotidal charts

Ricardo de Camargo^I; Joseph Harari^II

^IInstituto de Astronomia, Geofísica e Ciências Atmosférica da Universidade de São Paulo, Departamento de Ciências Atmosféricas (Rua do Matão, 1226, 05508-900, São Paulo, SP, Brasil)

^IIInstituto Oceanográfico da Universidade de São Paulo (Praça do Oceanográfico, 191, 05508-900, São Paulo, SP, Brasil)

ABSTRACT

The tidal circulation in Paranagua Bay (Parana State, Southern Brazil) was studied based on the Princeton Ocean Model. The model domain covered the near shore region and the estuarine area, with about 1 km grid resolution in cross-shore and along-shore directions. Homogeneous and diagnostic distributions for temperature and salinity were used and 12 tidal constituents were considered to specify the elevations at the open boundaries. Tidal analysis of 29-days time series of elevations and currents for each grid point generated corange and cophase lines as well as the correspondent axes of the current ellipses for each constituent. These computed values reproduced well the observed amplifications and phase lags of surface elevations and currents. Residual flows show the formation of tidal eddies, related to coastal geometry and bottom topography.

Descriptors: Paranagua Bay, Numerical modeling, Cotidal charts, Tidal current ellipses.

RESUMO

A circulação de maré na Baía de Paranaguá (Estado do Paraná, sul do Brasil) foi estudada através do Princeton Ocean Model. O domino do modelo abrange a região costeira adjacente e a área estuarina, com resolução de aproximadamente I km nas direções perpendicular e paralela à costa. Distribuições homogêneas e diagnosticas para temperatura e salinidade foram usadas e 12 constituintes de maré especificaram as elevações de superfície nos contornos abertos. Análises de maré de séries temporais de 29 dias de elevações e correntes para cada ponto de grade geraram linhas cotidais de amplitude e de fase, assim como elipses de correntes, para cada constituinte. Os valores obtidos pelo modelo reproduziram satisfatoriamente as amplificações e defasagens observadas nas elevações e correntes de superfície. Fluxos residuais mostram a formação de vórtices de maré, relacionados com a geometria da costa e a topografia do fundo.

Descritores: Baía de Paranaguá, Modelagem numérica, Mapas cotidais, Elipses de maré.

Introduction

The Paranagua estuarine complex, hereafter called Paranaguá Bay, is located on the Southern Brazilian coast, around 25°30'S and 048°30'W (Fig. 1). Both coastal geometry and bottom topography yield to a considerable amplification of the tidal range, from the mouth to the interior points. The present work is part of a recent study of tidal and wind-driven circulation at the bay using a three-dimensional numerical model (Blumberg & Mellor, 1983, 1987). For the purposes of this paper, only tidal characteristics will be investigated, on a comprehensive study which also leads to a system for tidal predictions over the entire estuary.

Similar researches of tidal features in estuaries can be found in Flather & Heaps (1975) for the Morecambe Bay (UK), in Oey et al (1985a, b and c) for the Hudson-Raritan Estuary (USA) and in Galperin & Mellor (1990a and b) for the Delaware Bay (USA), among others. For Paranaguá Bay, some previous works on numerical modeling of estuarine circulation can be mentioned, such as Marone & Jamiyanaa (1997), Oliveira (1998) and Camargo (1998); the two formers considered only tidal aspects, while the latter included the wind-driven circulation.

Brief description of Paranagua Bay

The estuary has two harbors, one of them with great importance for grain exportation, and also an oil terminal; there is a rich mangrove ecosystem and a large number of small communities living with local resources.

The estuary is formed by two branches (Fig. 1); the east-west branch has an extension of more than 40 km and includes both mentioned harbors; the north-south branch is shorter (~30 km) but wider than the east-west branch (8 km maximum). The mean depth is about 3 meters, but the navigational channels reach 12 to 15 m, maintained by frequent dredging; there are many areas of low depths and the amount of sediment is considerable.

The estuary is partially mixed and laterally homogeneous, with 612 km- of liquid surface and a drainage area of more than 3000 km². In a cross section in front of Paranagua Harbor, the observed values of surface water temperature are between 23° to 30°C in summer and between 18° to 25°C in winter; the salinity extremes are 12 and 29 in summer, 20 and 34 in winter (Knoppers el al., 1987).

The range of the astronomical tide is about 1.7 m in the mouth, but in the most interior areas it may reach 2.7 m; this tidal influence is associated with a tidal prism of 1.3 km³. The ratio between the main diurnal and semi-diurnal constituents is 0.23, but it is necessary to mention the non-linear effects generating shallow water constituents with large amplitudes at the fourth-diurnal frequencies; they are bigger than the diurnal ones and represent an important effect in surface elevations and associated currents (Marone et al, 1995^*; Camargo, 1998).

The meteorological station of Paranagua registers the higher precipitation rates of the South Brazilian coast, with mean value of 1988 mm/year (Climanálise, 1986). There are notable dry and wet seasons, respectively from January to March and from June to August. In spite of the large amount of rain, the fresh water discharge in the east-west branch varies from 70 to 200 m³ (in the north-south branch there are few measurements); these values indicate a large water absorption by the soil in the whole region.

Finally, the winds arc mainly controlled by the South Atlantic Subtropical High and the passage of synoptic and sub-synoptic systems; the prevailing winds blow from ENE. E, ESE and SE, with mean intensity of 4 m/s (Camargo & Marone, 1995^**); the sea breeze cell has amplitudes of about 2 m/s from November to March (Camargo el al., 1996).

The model and the domain of interest

The POM96 version of the numerical model of oceanic circulation described in Blumberg & Mellor (1983, 1987) was employed. The model is 3-D, based on a bottom-following σ coordinate system, and includes: an external-internal splitting mode algorithm; free-surface; a second moment turbulence closure model, for providing vertical mixing coefficients; Smagorinsky paramelrizations for the horizontal diffusion; explicit leapfrog method for the numerical integration in time and horizontal; and implicit method for the vertical integration of the equations. Other useful information about the model can be found in Galperin & Mellor (1990).

The study area of this work was represented by a constant horizontal resolution mesh (926 m), oriented cross and along-shelf, respectively with 70 and 90 points. Although the thermohaline distribution has been considered constant along the time integration, the model run in a 3-D domain. The vertical levels were defined according to Oey el al. (1985a) and Galperin & Mellor (1990), being 11 a layers located at 0, 0.03125, 0.0625, 0.125, 0.25. 0.5, 0.75, 0.875, 0.9375, 0.96875 and 1 times each local depth. The external mode time step was 15s and the corresponding internal mode time step was 300s. A 31-day simulation with this configuration required about 8 h of CPU on a workstation IBM RiscSystem 6000 model 43 P-140.

Tidal forcing and boundary conditions

The specification of the tidal oscillations at the open boundaries was based on cotidal charts of the main constituents for the Southeastern Brazilian shelf, obtained by Harari & Camargo (1994). These authors considered Ql, O1, P1, Kl, N2, M2, S2, K2 and M3 as the most important constituents for tidal motions in the Southeastern Brazilian shelf and obtained the correspondent values of amplitudes and phases relative to Greenwich, with a 13.9 km resolution. This data set was linearly interpolated for the borders of the Paranagua Bay domain.

Besides these. other high frequency constituents - MN4, M4 and MS4 - had to be included, because their large amplification in the interior areas of the estuary. These fourth-diurnal constituents are comparable to the diurnal ones (or even larger) and represent an important contribution for the tidal circulation. Their amplitudes and phases at the borders were estimated through a few available harmonic analyses of tidal series measured in the estuary and one in the adjacent shelf.

The complete set of 12 tidal constituents was used to compute tidal oscillations at the borders, although no information on tidal currents was available. Due to the lack of current data, surface tidal oscillations were imposed at the open boundaries and first internal grid points, as done by O'Connor (1991). This solution for boundary conditions led to good results for tidal simulations in Santos (Brazil, 24°S), yielding the obtention of corange and cophase lines for each constituent (Harari & Camargo, 1998). In the present case, the tidal influence at the borders was performed through the use of a partially clamped boundary condition, as presented in Chapman (1985) and in Blumberg & Kantha (1985); in this solution, the tidal signal is imposed as a restoring term added to radiational conditions, e.g. a nudging term, as follows:

where c = (gh)^1/2 is the phase velocity (dependent on the gravity acceleration g and local depth h), n is the normal to the border, η is the free surface elevation at the border, η_tide is the correspondent harmonic tide and T_f is associated with the weight of the nudging term (note that T_f → 0 indicates clamped boundaries and T_f → ∞ leads to a pure radiational condition). The term on the right represents damping which tends to force the value of η at the boundary to the known value η_tide with a time scale of the order of T_f.

Sensitivity studies were conducted to determine an appropriate determination of T_f: where values of 900s, 1800s, 3600s and 7200s were employed. The value of 1800s was considered the best choice for this domain, which keeps the tidal characteristics at the boundaries and allows the inclusion of wind fields at the surface (although not used in the present work).

The boundary condition for the vertically intenrated currents (the external mode) was the pure radiational condition applied only in the last border lines (which expression is the same as the equation mentioned above, with T_f=0). For the internal modes velocities, no-gradient boundary conditions were adopted.

Results

The model was initialized from rest at 00hGMT of March 1^st. 1997 for a 31-days simulation. The spin-up time related to the tidal forcing was about 30 hours; the subsequent 696 hours (29 days) time series of surface elevations and currents, for each grid point, were analyzed through the harmonic method (Franco, 1988). Comparisons between model results and observed harmonic constants at five points presented good agreement (Tab. 1 with positions signed in Fig. 1). With these observations, one can note the amplification rates of all constituents, except the diurnal ones; the model suggested satisfactory values for these amplification rates, except for MS4.

For surface currents, two points with measurements indicated a reasonable agreement (Tab. 2); in any case, it is a satisfactory representation, considering some limitations of the simulations, specially the grid size and uncertainties in the boundary data. In general, the model underestimated the current intensities, which are directly related to the horizontal resolution.

The composition of the harmonic constants of surface elevation and surface currents yield the charting of corange and cophase lines, as well as the axes of surface currents ellipses of each constituent. The major constituents M2, S2, M3 and M4 were selected for presentation (Figs 2 to 5). For these constituents, the amplification rates (from the coastal region to the interior areas) are very significant, being respectively 50%, 40%, over 100% and again over 100%. The region of the mouth of the estuary and the narrow interior channels presented the largest intensities of tidal currents, for all mentioned constituents.

Another feature of interest of the simulation is the residual current field associated to the tidal cycles. The configuration shown in Figure 6 was obtained through the computation of the mean currents over 10 tidal cycles, for each grid point, during the spring tide of the period of 6 to 10 March 1997. In this figure, one can note the presence of tidal eddies as the ones described by Zimmerman (1981), which can be generated by coastal geometry and bottom topography (since there are no wind and no thermohaline effects in the present simulation). In any case, the nonlinear inertia effects are fundamental to the formation of these eddies, as stated by Tee (1976) and Oey el al. (1985a). For example, the eddies formed inside the bay appear to be mainly induced by bottom topography, while the eddies formed in the region of the mouth appear to be induced by the combination of both bathymetry and coastal line geometry. For the present case, the residual currents can reach 10 to 15% of the total tidal flow, being very important for pollutant dispersion and sediment transport in the bay.

Discussions and conclusions

The use of partially clamped boundary conditions allowed the modeling of tidal propagation in the Paranagua Bay; the major advantage of this solution is the possibility of inclusion of surface wind fields with the same formulation.

In terms of surface elevation, the tidal ranges and amplification rates, from the mouth to the interior points, were well represented, as well as the correspondent phase lags. For surface tidal currents, the comparison of modeled and observed harmonic constants indicated a reasonable representation; in respect to the discrepancies in the comparisons of currents, it is necessary to mention that the adopted horizontal resolution is not fine enough to represent some cross-channel vertical sections (mainly bottom features and coastline details). In any case, the corange and cophase lines besides the axes of surface ellipses, for each tidal constituent, yield to tidal forecasts with considerable precision, for the whole domain.

Residual tidal currents were determined for the entire region, indicating the integrated effects of tidal cycles. The obtained tidal eddies are result of nonlinear interactions between the tidal flow, the coastal geometry and the bottom topography (since no winds and thermohaline effects were considered).

The next steps of this research will consider (i) the influence of surface wind fields in driving currents within the estuary and (ii) the inclusion of baroclinic effects through the use of temperature and salinity as prognostic variables (e.g., varying during tidal cycles and depending on wind conditions). These future simulations are planned to run on higher resolution domains, in order to improve the representation of effects due to bottom topography and coastal geometry.

Acknowledgements

The authors are very grateful to Dr. Eduardo Marone, of the Center for Marine Studies of the Federal University of Parana, for his cooperation and for supplying data field, and to Dr. Edmo J. D. Campos, of the Institute of Oceanography of the University of Sao Paulo, for his cooperation and computational support. The authors also would like to thank to FAPESP and CNPq, for providing fellowships to the authors.

(Manuscript received 10 December 2001; revised 04 November 2002; accepted 12 February 2003)

Contr. nº 850 do Inst. oceanogr. da Usp.

References

• Blumberg, A. F. & Kantha, L. H. 1985. Open boundary conditions for circulation models. J. Hydr. Engin., 111(2):237-255.
• Blumberg, A. F. & Mellor, G. L. 1983. Diagnostic and prognostic numerical circulation studies of the South Atlantic Bight. J. Geophys. Res., 88(C8):4579-4582.
• Blumberg, A. F. & Mellor, G. L. 1987. A description of a three-dimensional coastal ocean circulation model. In: Heaps, N. ed. Three-dimensional coastal ocean models, Washington, D.C., American Geophysical Union, 4:1-16.
• Camargo, R. 1998. Estudo numérico das circulações atmosférica e oceânica na região da Baía de Paranaguá. PhD Thesis. Universidade de São Paulo, Instituto Astronômico e Geofísico. 181 p.
• Camargo, R.; Marone, E. & Silva Dias, P. L. 1996. Detecção do sinal de brisa no registro de vento de Pontal do Sul (PR). In: CONGRESSO BRASILEIRO DE METEOROLOGIA, 9. Rio de Janeiro, 1996. Anais. Rio de Janeiro, Sbmet. p.1036-1040.
• Chapman, D. C. 1985. Numerical treatment of cross-shelf open boundaries in a barotropic coastal ocean model. J. phys. Oceanogr., 15(8):1060-1075.
• Climanálise. 1986. Boletim de Monitoramento c Análise Climática. Número Especial. São Paulo, MCT/INPE.
• Flather, R. A. & Heaps, N. S. 1975. Tidal computations for Morecambe Bay. Geophys. .J1. R. astr. Soc. 41:489-5 17.
• Franco, A. S. 1988. Tides: fundamentals, analysis and prediction. 2 ed. São Paulo, FETH. 249p.
• Galperin, B. & Mellor, G. L. 1990a.    A time-dependent, three-dimensional model      of the Delaware Bay and river system.      Part 1: Description of the model and tidal    analysis. Estuar. coast, shelf Sci. 31(3):23 1-253.
• Galperin, B. & Mellor, G. L. 1990b. A time-dependent, three-dimensional model of the Delaware Bay and river system. Part 2: Three-dimensional flow fields and residual circulation. Estuar., coast, shelf Sci., 31 (3):255-281.
• Harari, J. & Camargo, R. 1994. Simulação da propagação das nove principais componentes de maré na plataforma sudeste brasileira através de modelo numérico hidrodinâmico. Bolm Inst. oceanogr., S Paulo. 42( 1):35-54.
• Harari, J. & Camargo, R. 1998 Modelagem numérica da região costeira de Santos (SP): Circulação de maré. Rev. bras, oceanogr., 46(2):135-156.
• Knoppers, B. A.; Brandini, F. P. & Thamm, C. A. 1987. Ecological studies in the Bay of Paranaguá II. Some physical and chemical characteristics. Nerítica, Curitiba, 2(1):1-36.
• Marone, E. & Jamiyanaa, D. 1997. Tides and tidal modelling in the Paranaguá Bay, PR. Neritica, Curitiba, 11(1):95-108.
• O'Connor, W. P. 1991. A numerical model of tides and storm surges in the Rio de la Plata Estuary. Continent. Shelf Res., 1 1(12):1491-1508.
• Oey, L.-Y.; Mellor, G. L. & Hires, R. I. 1985a. Tidal Modelling of the Hudson-Raritan Estuary. Estuar. Coast, shelf Sci.,20:511-527.
• Oey, L.-Y.; Mellor, G. L. & Hires, R. l., 1985b. A Three-dimensional simulation of the Hudson-Raritan Estuary. Part I: Description of the model and model simulations. J. Phys. Oceanogr., 15(12):1676-1692.
• Oey, L.-Y.; Mellor, G. L. & Hires, R. I. 1985c. A Three-Dimensional Simulation of the Hudson-Raritan Estuary. Part II: Comparison with observation. J. Phys. Oceanogr., 15(12):1693-1709.
• Oliveira, F. C. T. 1998. Estudo de correntes de maré no complexo estuarino da Baía de Paranaguá. MSc. Dissertation. Universidade Federal do Rio de Janeiro. 92p.
• Tee, K. T. 1976. Tide-induced residual current, a 2-D nonlinear numerical tidal model. J. Mar. Res., 34(4):603-628.
• Zimmerman, J. T. F. 1981. Dynamics, diffusion and geomorphological significance of tidal residual eddies. Nature, 290:549-555.

(*

) Marone. E.: Guimarães. M. R. F.; Camargo. R., Prata Júnior. V. P. & Klingenfuss. M. S. 1995. Caracterização física das condições oceanográlicas, meteorológicas e costeiras das zonas estuarinas da Baía de Paranaguá - PR. Brasil. In: CONGRESO LAT1NO AMER1CANO DE CIÊNCIAS DEL, MARCOLACMAR. 6. Mar del Plata. 1995. Mar del Plata. Universidad Nacional de Mar del Plata, p. 129.

(**

) Camargo. R. & Marone, E. 1995. Comportamento do vento de superfície em Pontal do Sul (PR. Brasil) -Análise preliminar. In: CONGRESO LATINOAMERICANO DE CIENCIAS DEL MAR- COLACMAR. 6. Mar del Plata. 1995. Resumos. Mar del Plata. Universidad Nacional de Mar del Plata, p. 63

Publication Dates

• Publication in this collection
  23 Apr 2013
• Date of issue
  2003

History

• Accepted
  12 Feb 2003
• Received
  10 Dec 2001
• Reviewed
  04 Nov 2002