Estrutura dinâmica das contracorrentes equatoriais do Oceano Atlântico ao longo de 44W

Estrutura dinâmica das contracorrentes equatoriais do Oceano Atlântico ao longo de 44W

Domingos Fernandes Urbano Neto^I; Ilson Carlos Almeida da Silveira^II

^ILaboratório de Dinâmica Oceânica – Departamento de Oceanografia Física, Química e Geológica – Instituto Oceanográfico - Universidade de São Paulo – Praça do Oceanográfico, 191– Cidade Universitária – CEP: 05508-900 – São Paulo - SP, Brasil, ou – Caixa Postal 9075 - 01065-970 – São Paulo - SP, Brasil – Tel.: [phone] – FAX: [phone]. E-mail: [e-mail]

^IILaboratório de Dinâmica Oceânica – Departamento de Oceanografia Física, Química e Geológica – Instituto Oceanográfico - Universidade de São Paulo – Praça do Oceanográfico, 191– Cidade Universitária – CEP: 05508-900 – São Paulo - SP, Brasil, ou – Caixa Postal 9075 - 01065-970 – São Paulo - SP, Brasil – Tel.: [phone] – FAX: [phone]. E-mail: [e-mail]

RESUMO

Estudos teóricos ou baseados em modelagem numérica indicam que a estrutura vertical do oceano tropical é dominada pelos dois primeiros modos baroclínicos. Entretanto, faltam resultados quantitativos que mostrem efetivamente a dominância destes modos através de comparações com observações de velocidade no Oceano Atlântico. Neste artigo, utilizamos perfilagens de Conductivity-Temperature-Depth (CTD) simultâneas às de velocidade (Pegasus), realizadas durante o Projeto Western Tropical Atlantic Experiment (WESTRAX), para avaliar a estrutura modal das contra-correntes equatoriais ao longo de 44W: as Contracorrentes Norte Equatorial (CCNE) e Subcorrente Norte Equatorial (SNE). Os resultados mostram que, de fato, os dois primeiros modos baroclínicos são mais importantes que o modo barotrópico e os demais modos baroclínicos. A soma dos três primeiros modos dinâmicos explica cerca de 50% da estrutura vertical das correntes. Talvez devido à proximidade da camada de limite oeste oceânica, e portanto sítio de origem tanto da CCNE quanto da SNE, os resultados apresentam nítida variação sazonal. De fato, os três primeiros modos baroclínicos respondem por cerca de 63% da estrutura vertical de velocidade observada, sendo que o primeiro modo é responsável por 35-45 % do sinal medido. Esses resultados corroboram aqueles encontrados para o cruzeiro de outono boreal em estudos anteriores. No entanto, para os cruzeiros de inverno boreal, onde uma estrutura de velocidade apresenta-se mais complexa, é necessária a introdução de maior número de modos dinâmicos. Tanto para o cruzeiro de outono quanto para o de inverno houve necessidade seis modos dinâmicos para reproduzir 63% da estrutura de velocidade observada. Os três primeiros modos dinâmicos são responsáveis por 30-40% do sinal observado. Estudos comparativos entre a estrutura dinâmica das correntes equatoriais em 44W e em porções mais centrais da bacia do Oceano Atlântico, devem ser conduzidos para estabelecer o papel dos modos de mais alta ordem e sua relação com a camada limite oeste.

Palavras-chave: Dinâmica Equatorial; Modos Dinâmicos; Meso-escala; Correntes Equatoriais.

ABSTRACT

It is commonly mentioned in the literature that the first and second baroclinic modes dominate the vertical structure in the tropical ocean. However, most of these studies refer to either theoretical or numerical modeling. There are no quantitative results that confirm the two internal mode dominance through comparisons with velocity observations in the Tropical Atlantic. In this article, we use hydrographic (CTD) and velocity (PEGASUS) profiles obtained during the "Western Tropical Atlantic Experiment-WESTRAX" to evaluate the modal structure of the equatorial currents along 44W, in particular, the North Equatorial Countercurrent (NECC) and the North Equatorial Undercurrent (NEUC). The results show that indeed the first two baroclinic modes explain about 50% of the vertical structure of the currents on average. We also find indication of a clear seasonal variation on the modal composition, a probable consequence of the proximity of the study region to the western boundary layer, and therefore, the source of both NECC and NEUC. As a mather of fact, the first three dynamical mode account for about 63% of the observed velocity vertical structure, while the first mode alone is responsible for 35-45 % of the measured velocity. These results reasonably agree with previous studies also using the fall WESTRAX cruise. However, the more complex boreal winter cruise velocity patterns yield an increased number of dynamical modes to reconstitute the velocity field. For the both fall and winter cruises, we need six dynamical modes to reproduce 63% of the observed velocity structure. The first three dynamical mode are responsible for 30-40% of the observed signal. A comparative study of the dynamical structure of the zonal equatorial currents at 44W and at more central portions of the Atlantic Ocean basin should be conducted in order to establish the role of the higher order modes and its relation with the western boundary layer.

Keywords: Equatorial Dynamics; Dynamical Modes; Mesoescale; Equatorial Currents.

INTRODUÇÃO

A região equatorial oeste do Oceano Atlântico tem sido fortemente investigada ao longo das últimas décadas. Correntes de Contorno Oeste desempenham um importante papel como elo de ligação na circulação termo-halina global. Processos como o de retroflexão do sistema Corrente Norte do Brasil (CNB), além de alimentarem em diferentes níveis picnoclínicos a Contracorrente Norte Equatorial (CCNE) e a Subcorrente Norte Equatorial (SNE), emitem anéis de núcleo quente em direção ao Mar do Caribe. Parte da Corrente Norte Equatorial (CNE), fluindo em direção à borda do continente sul-americano, reverte ciclonicamente e contribui de forma significante no transporte da CCNE. Torna-se evidente a complexidade na estrutura vertical das correntes desta região equatorial. Como conseqüência desta estrutura vertical, Philander (1989) resume que as regiões equatoriais, não só do Oceano Atlântico, se contrapõem às regiões subtropicais, apresentando os dois primeiros modos baroclínicos como dominantes, e o modo barotrópico possuindo menor importância. No entanto, estudos realizados a respeito da estrutura modal são baseados em resultados obtidos através de modelagem numérica aplicada a oceanos idealizados, baseados em modelos teóricos, ou em interpretação de simulações numéricas de modelagem processual.

Lighthill (1969) usou modelagem numérica aplicada a região tropical, fazendo uso de uma estratificação básica mais próxima de uma situação real, onde a termoclina se apresenta menos acentuada, mostrou que o primeiro modo baroclinico é o mais importante na estrutura modal de correntes equatoriais. Estudos analíticos realizados por Cane e Sarachik (1976); Cane e Sarachik (1977); Cane e Sarachik (1979) utilizam um modelo de águas rasas de primeiro modo baroclínico para examinar a geração de correntes lineares em um oceano equatorial. Cane (1979) estudou a resposta de um oceano equatorial sob influência do vento usando um modelo numérico de 1-1/2 camada. Esta estrutura de massa, obviamente implica num primeiro modo baroclínico importante.

No modelo de equações primitivas de Philander e Pacanowski (1980), uma estratificação básica é assumida e representada por um perfil de freqüência de Brunt-Väisälä idealizado. Segundo os autores, os modos naturais de oscilação desempenham um papel central no ajustamento do oceano às variações dos ventos de superfície. Afirmam ainda, que o segundo modo baroclínico é dominante neste processo por apresentar o maior valor de amplitude modal em regiões próximas à superfície, e que devido a isso, qualquer destúrbio nas camadas superiores do oceano resulta em uma grande projeção dentro do segundo modo baroclínico, que por esta razão, pode ser usado para representar a estrutura dinâmica do oceano tropical, não apenas no centro da bacia oceânica de baixas latitudes, mas em qualquer parte desta.

Mais recentemente, McClean e Klinck (1995), usando o modelo CME (WOCE Community Modeling Effort general circulation), mostraram que a região da CCNE é dominada pelo primeiro e segundo modos baroclínicos. Patti (2001) usou um modelo de três camadas – e portanto, de três modos dinâmicos – para reproduzir o processo de emissão de vórtices pela CNB, através da colisão de ondas de Rossby de primeiro e segundo modos na região do contorno oeste.

Em termos de análise da estrutura de velocidade observada no Atlântico Tropical, Bourlès e outros (1999) descreveram observações de perfilagem de velocidade ao longo de 44W a partir dos cruzeiros WESTRAX e um cruzeiro hidrográfico realizado pela NOAA (National Oceanic and Atmospheric Administration), em agosto de 1989 (WILSON et al., 1994). Estes autores examinaram o campo de velocidade na camada superior do oceano e apresentaram características inéditas da CCNE. O resultado mais importante deste estudo foi a descrição de um comportamento quase-sazonal do sistema zonal de correntes equatoriais. Durante os cruzeiros de outono/verão boreais, as CCNE e SNE estão fusionadas e apresentam um núcleo em subsuperfície. Durante os cruzeiros de inverno boreal, estas correntes estão separadas por pelo menos três graus de latitude. Estes resultados indicam a existência de uma estrutura modal baroclínica complexa, provavelmente devido à proximidade do contorno oeste.

Silveira e outros (2000), realizando uma comparação entre os modos de variabilidade vertical estatísticos (Depth Empirical Orthogonal Functions - DEOFs) e os modos dinâmicos, calculados a partir de observações WESTRAX de outono boreal de 1990, da região de retroflexão da CNB, mostram que os dois primeiros modos estatísticos dessa função explicam 90% da variância na estrutura vertical. Um ajuste ótimo para os modos dinâmicos nos modos DEOF sugere que o primeiro e segundo modos dessa função são dominados (dinamicamente) pelos primeiro e segundo modos baroclínicos, respectivamente. Silveira et al (2000) mostraram ainda que a soma dos primeiros três modos dinâmicos é responsável por 75% da estrutura vertical dos dois primeiros modos da DEOF. Dessa forma, formularam um modelo de três camadas para representar o sistema de retroflexão da CNB.

Os resultados evidenciados anteriormente conduzem às seguintes questões:

• São os primeiro e segundo modos baroclínicos, como dito por Philander (1989), mais relevantes que o modo barotrópico na região equatorial do Oceano Atlântico Tropical Oeste?

• É a afirmação de Philander e Pacanowski (1980) – de que o segundo modo baroclínico é dominante – realmente válida? Pode este modo ser usado como representante da estrutura dinâmica para qualquer porção do Oceano Atlântico Tropical?

• Os três primeiros modos dinâmicos (barotrópico e dois primeiros baroclínicos) foram suficientes para reproduzir de forma satisfatória o processo de retroflexão do sistema CNB (SILVEIRA et al., 2000). O mesmo pode ser aplicado para o sistema de correntes zonais ao longo de 44ºW?

• De que forma o núcleo de subsuperfície da CCNE, relatado por Bourlès e outros (1999), está relacionado à estrutura modal da região?

Buscando responder essas questões, uma análise de modos normais foi realizada para um dos únicos conjuntos de dados quasesinóticos que envolvem medições diretas de correntes, existente na região tropical oeste do Oceano Atlântico, o conjunto WESTRAX. Observações hidrográficas e de velocidade obtidas de forma simultânea, tornam possível uma investigação detalhada e precisa a respeito da estrutura dinâmica ao longo da seção de 44ºW, ou seja, uma porção do oceano equatorial localizada entre dois importantes processos dinâmicos: a retroflexão do sistema CNB e o sistema de correntes zonais equatoriais.

CONJUNTO DE OBSERVAÇÕES

O conjunto de dados consiste de componentes de velocidade obtidas por perfilador de correntes PEGASUS (SPAIN et al., 1981) e amostragens simultâneas, em toda coluna de água, de temperatura e salinidade medidas através de CTD (Conductivity-Temperature-Depth). Estes dados foram obtidos durante o Western Tropical Atlantic Experiment – WESTRAX (BROWN et al., 1992). O programa experimental completo consiste de cinco cruzeiros realizados de 1990 a 1991 (Tabela 1). Neste trabalho, limitamo-nos à análise dos quatro primeiros cruzeiros.

O Projeto WESTRAX amostrou a região do Oceano Atlântico compreendida pela região à oeste de 44ºW e ao norte da costa sulamericana, até a latitude de 15N (Figura 1). As perfilagens foram dispostas em um conjunto de radiais, formando uma malha hidrográfica. Para este trabalho, foram analisadas as estações ao longo de 44ºW, uma vez que registraram o sinal do sistema zonal de correntes equatoriais próximo à sua região de origem, ou seja, a região entre a borda oeste e o centro da bacia.

Entre os mais recentes trabalhos utilizando o conjunto de dados provenientes do Projeto WESTRAX, encontra-se o de Bourlès e outros (1999). Os autores examinaram o campo de velocidade da camada superior do oceano, utilizando dados de ADCP de bordo (Shipboard Acoustic Doppler Current Profiler) e de CTDO2 (Conductivity-Temperature- Depth-Oxygen). O oceano superior foi dividido em porções, onde as fronteiras verticais foram, além da interface oceano-atmosfera, as superfícies das isopicnais de 24,50 e 26,75 kg m. Seções verticais de velocidade ao longo de 44W foram apresentadas e descritas em detalhe. Aqui, apresentamos seções da componente zonal de velocidade análogas àquelas apresentadas por Bourlès e outros (1999), porém utilizando os dados do perfilador PEGASUS (Figura 2). Ambas isopicnais foram sobrepostas aos contornos de velocidade e são indicadas pela linha espessa. As isotacas apresentam intervalo de 0.1 m s, com valores positivos indicando fluxo para leste e negativos para oeste. A localização das estações é sinalizada pelos círculos desenhados na parte superior do campo.

O painel da Figura 2 mostra claramente duas fases distintas no comportamento das correntes CCNE e SNE. Durante o período de inverno boreal (cruzeiros de jan90 e jan91; Tabela 1), o núcleo da CCNE está centrado em 3,2N. A SNE encontra-se deslocada para norte, com seu núcleo posicionado além da fronteira norte do campo, em 1990, e centrado em 6,8N, em 1991. Estão notoriamente separados por uma distância mínima de quatro graus de latitude. Por outro lado, nos cruzeiros de primavera e outono (set90 e Jun91), as duas correntes estão alinhadas verticalmente com um forte núcleo posicionado acima da isopicnal de 24,5 kg m (Figura 2 linha espessa). Bourlès e outros (1999) descrevem a presença de uma nítida separação vertical entre os núcleos dessas correntes durante o cruzeiro Jun91. Uma descontinuidade na isopicnal de 26,75 (Figura 2, linha espessa), em 3,2N, indica a presença da SNE durante esses cruzeiros.

Ao norte do sistema CCNE existe um escoamento para oeste associado à borda sul da CNE. Esta corrente aparece mais evidente durante os cruzeiros de primavera e outono (de set90 e jun91), quando a CCNE desloca-se para sul devido ao posicionamento da Zona de Confluência Intertropical (ZCIT).

Como as perfilagens de PEGASUS cobriram toda coluna de água, ao contrário dos dados de ADCP que atingem profundidades de no máximo 500 metros, apresentamos painéis complementares que mostram a circulação abaixo dos 1000 metros de profundidade. Atentamos que estes gráficos apresentam escala vertical distinta da porção superior, no entanto, fica evidente um padrão de velocidades baixas na maior parte da coluna de água e que as regiões de maior cisalhamento encontramse confinadas aos primeiros 1000 metros de profundidade.

DECOMPOSIÇÃO MODAL

O cenário aqui envolvido consiste de um oceano sob aproximação hidrostática, continuamente estratificado, com tampa rígida e fundo plano como condições de contorno para a superfície e fundo, respectivamente. A estrutura vertical de velocidade foi decomposta em modos dinâmicos, seguindo formalismo apresentado por Silveira e outros (2000). As componentes da equação do movimento linearizadas, sob aproximação de Boussinesq, podem ser desacopladas em estruturas vertical e horizontal, aplicando-se o método de separação de variáveis, descrito em LeBlond e Mysak (1978):

onde P_i, V_i e W_i são as amplitudes do i-ésimo modo vertical da pressão e das velocidades zonal e vertical, respectivamente; F_i é o i-ésimo modo vertical; e N é a freqüência de estratificação ou Freqüência de Brünt-Väisälä.

Os modos verticais satisfazem a equação de autovalores (li)

As condições de contorno para F₁ são fundo plano e tampa rígida em superfície:

A Equação (4) juntamente com as condições de contorno representadas por (5) constituem um problema de autovalores para l_i. Assim, existe um número infinito de soluções F_i(z), para i=0,1,2,..., cada uma associada a um autovalor l_i real e discreto. Portanto, os modos dinâmicos são funcionais da freqüência de Brünt-Väisälä.

As amplitudes U_i são obtidas através da projeção dos perfis de velocidades observadas por PEGASUS nos modos dinâmicos, através da integração:

onde u é a componente zonal da velocidade observada.

Nas seções a seguir descreve-se uma discussão sobre a estrutura de estratificação da região WESTRAX e cálculo dos modos dinâmicos associados a esta.

ESTRUTURA VERTICAL DA ESTRATIFICAÇÃO

Perfis médios sazonais de densidade potencial (s_q) foram calculados a partir dos dados de salinidade e temperatura obtidos por CTD. Um perfil médio dos quatro cruzeiros também foi calculado e denominado de média anual. Perfis de Freqüência de Brunt-Väisälä (N) foram obtidos através da equação:

onde g é a aceleração devido à gravidade, r₀ é a densidade potencial média da água do mar. Os perfis de N e da densidade convencional s_q estão dispostos na Figura 3. Os perfis sazonais das freqüências de empuxo estão representados (Figura 4), pela curva mais espessa. A curva mais delgada representa o perfil médio.

Em todos os cruzeiros, um pico acentuado confinado na parte mais superior da coluna de água e centrado em 125 m, indica uma alta estratificação em subsuperfície. A presença da picnoclina permanente, típica de baixas latitudes (EMERY et al., 1983; HOURY et al., 1987), está representada pelo forte gradiente entre 2 e 9 ciclos por hora (cph).

O perfil médio anual (Figura 4 – curva delgada) foi sobreposto aos perfis médios de cada cruzeiro como referência climatológica. Existem variações sazonais nos diferentes cruzeiros. O cruzeiro de set90 apresenta o pico próximo de 10 cph, mantendo altos valores de estratificação até a superfície. Situação similar é encontrada no perfil de Jun91, porém, próximo à superfície, a freqüência cai rapidamente de 10 para 6 cph. Em todos os cruzeiros, em regiões mais profundas que 250 metros, encontram-se valores de freqüência tão baixos quanto 2 cph, permanecendo até o fundo. As curvas sazonais se sobrepõem à curva média anual nestas profundidades.

Segundo Emery e outros (1983) e Houry e outros (1987), regiões de médias e altas latitudes apresentam variabilidade sazonal de grande amplitude em contraposição às regiões de baixas latitudes. Desse modo, seria adequado, então, assumir uma estratificação básica para o conjunto de dados WESTRAX? Para verificar se as oscilações sazonais descritas anteriormente, assim como a proximidade do equador, afetam as amplitudes modais e os raios de deformação, impedindo a aplicação do perfil médio, apresentamos a próxima seção.

RAIOS DE DEFORMAÇÃO E MODOS VERTICAIS

Os autovalores l_i, os quais formam um conjunto infinito porém discreto, são relacionados aos raios de deformação internos Rd_i por:

Desta forma, para cada raio de deformação interno existe um autovalor e um automodo vertical associados. A Tabela 2 apresenta os valores dos raios de deformação internos para os primeiros sete modos baroclínicos obtidos para cada cruzeiro e também para o perfil médio anual. Os valores de Rd_i dos dois primeiros modos baroclínicos variam cerca de 3,5% da média ao longo do ano, enquanto os raios dos modos de maior ordem (3 – 7) praticamente inexistem. Os valores de raios obtidos a partir do perfil médio anual refletem esta similaridade, e sugerem pouca sensitividade ao cálculo relativamente às variações de N na porção superior da coluna de água.

Como a região de estudo apresenta-se dentro do cinturão equatorial (±1,5º de latitude), a seção hidrográfica foi separada em duas partes, afim de calcular os raios de deformação internos e automodos normais para duas faixas de latitudes separadas pelo valor de 5ºN. A parte mais próxima do equador utilizou o valor de q₀ centrado em 3ºN enquanto a segunda faixa utilizou um q₀ em 7ºN. A Tabela 3 mostra a variação dos valores de Rd_i nas diferentes faixas de latitude em comparação com os valores obtidos para a seção completa, compreendida entre 0º e 10ºN. Para este cálculo, foi utilizado o perfil médio anual da freqüência de estratificação. A escala dos movimentos, para o primeiro modo baroclínico, varia de 300km na região mais ao sul da seção (q₀ de 3ºN), vale aproximadamente 200km para um q₀ de 5ºN, e cai para um valor da ordem de 100km para a região mais ao norte da seção (q₀ de 7ºN). Isto é decorrente da variação de f₀, inversamente proporcional à Rd_i na Equação (8).

Assim como para cada raio de deformação existe um autovalor associado, também existe um autovetor ou modo (F_i(z)) correspondente. A Figura 5 apresenta um painel composto pelo perfil da freqüência de estratificação em toda coluna de água em conjunto com os sete primeiros automodos baroclínicos (linha cheia). As linhas tracejadas e pontilhadas são os modos calculados para as faixas de latitude mais próximas e mais distantes do equador, respectivamente. Fica evidente que a estrutura modal não se altera devido a proximidade do equador, apenas os raios baroclínicos crescem em função do decréscimo do parâmetro de Coriolis.

RECOMPOSIÇÃO MODAL DOS PERFIS DE VELOCIDADE

A recomposição dos perfis de velocidade, usando um número finito de modos, foi realizada através da equação (LEBLOND; MYSAK, 1978):

onde u* é a componente zonal de velocidade recomposta; n corresponde ao modo dinâmico em que a somatória é truncada. Seções de velocidade ao longo de 44ºW foram construídas através destes perfis recompostos utilizando até oito modos dinâmicos, ou seja, o modo barotrópico (n=0) e os sete primeiros modos baroclínicos (n=1,2, ... ,7).

Para identificar os modos dinâmicos que individualmente respondem pela estrutura de velocidade, bem como em termos da média ponderada destes modos na recomposição da estrutura do escoamento, foram calculadas as raízes do erro médio quadrático (remq) seguindo formalismo apresentado por Pinardi e Robinson (1987),

No entanto, valores de remq expressos em percentuais são mais intuitivos e foram calculados utilizando a relação:

A Tabela 4 apresenta os valores das raízes do erro-médio quadrático calculados para os perfis truncados com os oito modos dinâmicos.

Painéis com as seções da componente zonal de velocidade reconstruída para cada cruzeiro são apresentados nas figuras 6 a 9. Os campos exibem a recomposição modal (Equação (9)) de dois até oito modos dinâmicos (modo barotrópico e sete primeiros modos baroclínicos) em comparação com a seção de velocidade observada por PEGASUS (painel mais inferior à direita). As isotacas apresentam intervalo de 0.1m s^-1. Valores positivos indicam fluxo para leste e negativos para oeste. A posição das estações é indicada pelos círculos desenhados na parte superior dos campos.

Para o cruzeiro de jan90, o campo composto por dois modos dinâmicos (Figura 6 - 2MODOS) apresenta feições pobres e com valores de velocidade baixos. O modo barotrópico (BT) responde individualmente por 13,1% do sinal medido por PEGASUS e o primeiro modo baroclínico (1BC) por 35,8% (Tabela 4). Os dois primeiros modos dinâmicos (2MD) respondem juntos por apenas 37,9%. Tanto o campo truncado em três modos (Figura 6 - 3MODOS), quanto aquele em quatro modos (Figura 6 - 4MODOS), ainda não conseguem reproduzir o núcleo de subsuperfície da SNE. Isto está relacionado à pequena importância dinâmica do segundo e do terceiro modos baroclínicos, que respondem, respectivamente por 2,6% e 3,3%. Forte melhora da reprodução do sinal acontece quando a somatória modal inclui o quarto modo baroclínico (4BC), o qual explica 6,2% do sinal medido. O campo de cinco modos (Figura 6 - 5MODOS), explica 61,0% da variabilidade da estrutura vertical. Essa diferença é evidenciada tanto pela comparação qualitativa entre o campos truncados em quatro e em cinco modos (Figura 6), quanto quantitativamente pelo salto de 46,1% para 61,0% na Tabela 4. Os modos baroclínicos de maior ordem contribuem com valores inexpressivos inferiores a 2%, o que pode ser constatado pelo fato de não haver grandes diferenças entre os campos truncados em 5, 6, 7 e 8 modos.

Para o cruzeiro de set90 (Figura 7), a contribuição dos modos manteve o mesmo comportamento. O 1BC domina com 46,2% seguido pelo 2BC com 12,9%. O 4BC é mais alto que o 3BC. Porém, o modo BT mostra o menor valor entre todos os cruzeiros analisados. Os três primeiros modos dinâmicos (3MD), respondem, então, por 66,9% da estrutura vertical. Com a significante colaboração do 4BC, salta para o expressivo valor de 75,5%. Chega a 84,5% do sinal medido utilizando os oito primeiros modos dinâmicos, que é o mais alto valor percentual da tabela. Mais uma vez, o núcleo de subsuperfície do sistema CCNE-SNE somente é recomposto após a inclusão do quarto modo baroclínico.

Entretanto, para o cruzeiro de jan91 (Figura 8), a estrutura dinâmica mostra-se mais complexa e diferente daquela do cruzeiro jan90. Este fato foi inesperado, uma vez que ambos ocorreram exatamente na mesma estação: o inverno boreal. Nesse cruzeiro, é o segundo modo baroclínico aquele com maior amplitude, respondendo individualmente por 22,4%. O BT e 1BC responderam por inexpressivos 6,7% e 5,0%, respectivamente. O 5BC é responsável por 7,8% do sinal perfilado. O núcleo da SNE só atinge valores superiores a 0,5 m s no campo recomposto por seis modos dinâmicos. Todos os outros modos entre o BT e o 5BC respondem com o mesmo grau de importância, dando a sensação de que a energia do sistema foi espalhada nos primeiros modos dinâmicos, deixando o segundo modo dominar. O campo composto por oito modos dinâmicos explicou 67,7% da variabilidade do sinal amostrado. Foi o menor valor percentual entre os campos recompostos por oito modos dinâmicos. Isto indica que mais energia está espalhada em modos de maior ordem que não foram levados em consideração neste estudo.

Novamente, para o cruzeiro de jun91 (Figura 9), o 1BC é dominante, respondendo por 34,1%. De modo similar ao cruzeiro de set90, o segundo e o quarto modos baroclínicos se mostram importantes, respondendo por 21,3% e 4,3%, respectivamente. O modo BT aparece com seu menor valor entre todos os cruzeiros, 3,1%. Todos os outros modos colaboram com valores muito baixos, incapazes de ultrapassar valores de 1,1%. Da análise qualitativa podemos perceber que o núcleo de subsuperfície da CCNE é dissolvido após a inclusão do 7BC.

CONCLUSÕES

Utilizando perfilagens simultâneas de velocidade, salinidade e temperatura, o Método de Análise dos Modos Normais foi aplicado na região da CCNE/SNE. Os perfis da componente zonal de velocidade foram recompostos e comparados com o sinal medido. Análises qualitativas e quantitativas com reconstituição por truncamento a um número finito e crescente de modos dinâmicos, baseadas em figuras de mérito estatísticas, forneceram informações necessárias para responder as seguintes questões: a) São os primeiro e segundo modos baroclínicos, como dito por Philander (1989), mais relevantes que o modo barotrópico na região equatorial do Oceâno Atlântico Tropical Oeste? b) É a afirmação de Philander e Pacanowski (1980) de que o segundo modo baroclínico é dominante realmente válida? Pode este modo ser usado como representante da estrutura dinâmica para qualquer porção do Oceano Atlântico Tropical? c) Os três primeiros modos dinâmicos (barotrópico e dois primeiros baroclínicos) foram suficientes para reproduzir de forma satisfatória o processo de retroflexão do sistema CNB (SILVEIRA et al., 2000). O mesmo pode ser aplicado para o sistema de correntes zonais ao longo de 44ºW? d) De que forma o núcleo de subsuperfície da CCNE, relatado por Bourlès e outros (1999), está relacionado à estrutura modal da região?

Nossos resultados confirmam a afirmação de Philander (1989) de que os dois primeiros modos baroclínicos são dominantes em regiões equatoriais. Para o cruzeiro de jun91, o modo barotrópico explica apenas 3,1% da estrutura vertical de velocidade observada contra os 34,1% e 21,3% dos primeiro e segundo modos baroclínicos, respectivamente. Apenas para o cruzeiro de jan90, o modo barotrópico apresentou maior importância (13,1%).

Os resultados obtidos por Philander e Pacanowski (1980), através de modelagem numérica, indicaram que o segundo modo baroclínico é dominante devido ao fato de possuir maior amplitude em comparação com o primeiro e o terceiro modos baroclínicos, e conseqüentemente, pode representar a dinâmica do oceano equatorial em qualquer região da bacia atlântica. Nossos resultados mostram claramente que, ao longo de 44ºW, o segundo modo baroclínico não é dominante para três dos quatro cruzeiros analisados, e que para o cruzeiro de jan90, este modo foi responsável por apenas 2% da estrutura vertical do sinal de velocidade medido. Apenas para o cruzeiro de jan91, o segundo modo dinâmico apresentou maior importância que os demais, explicando 22,4% do sinal de velocidade em 44ºW. Comparado com o primeiro modo baroclínico, que respondeu por aproximadamente 35% nos cruzeiros de jan90 e jun91, e por 46% em jan91, o segundo modo baroclínico não domina. Este resultado corrobora a clássica afirmação de Lighthill (1969) de que o primeiro modo baroclínico domina a estrutura nos primeiros 100 metros de coluna de água.

A análise de modos normais revelou a existência de uma variação sazonal relacionada com a estrutura de velocidade. Os cruzeiros de set90 e jun91 (outono e verão boreais) apresentaram uma dominância clara dos dois primeiros modos baroclínicos, dando suporte a estudos prévios baseados em modelagem e teoria. De fato, os três primeiros modos dinâmicos foram responsáveis por 63% da estrutura de velocidade observada. Estes resultados concordam com os que Silveira et al. (2000) encontraram quando analisaram o cruzeiro de outono boreal (set90) para a região de retroflexão da CNB. Por outro lado, os cruzeiros de inverno boreal (jan90 e jan91) apresentaram um padrão de velocidade mais complexo. Para estes dois cruzeiros, são necessários seis (primeiros) modos dinâmicos para reproduzir 63% da estrutura de velocidade observada. Os três primeiros modos dinâmicos responderam por 30%-40% do sinal observado. Isto se deve ao fato de que, durante esses períodos, o processo de retroflexão da CNB não ocorreu, deixando de alimentar o sistema zonal de correntes equatoriais. Os valores de velocidade zonal são fracos e as CCNE e SNE apresentam-se separadas. Obviamente, quanto mais complexo o padrão de circulação, maior é o número de modos dinâmicos que detêm a energia da estrutura vertical.

A questão que se impõe: se a estrutura modal observada aqui, com participação significativa de modos de ordem mais alta, é resultado do fato de o meridiano de 44ºW estar próximo à camada limite oeste do Oceano Atlântico Tropical, e portanto, sítio de origem das CCNE e SNE. Comparação com a estrutura de velocidade observada em porções mais centrais da bacia oceânica, assim como a análise modal dessas, deve ser conduzida no intuito de confirmar os resultados aqui obtidos, que denotam maior complexidade dinâmica do que a teoria e a modelagem numérica aplicadas a região lhes conferiram até o momento.

Agradecimentos

Os autores expressam agradecimento ao Dr. Wendell Brown, por gentilmente ceder os dados WESTRAX, e à Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), pelos recursos concedidos através dos processos 2002/07292-2 e 2000/06985-9.

Submetido em 15, julho, 2003/Aceito em 13, abril, 2004

Submited July 15, 2003/Accepted April 13, 2004

NOTAS SOBRE OS AUTORES

Domingos Fernandes Urbano Neto é bacharel em Oceanologia na área de Gerenciamento Ambiental pela FURG, vinculado ao Laboratório de Oceanografia Física (de 1996 a 1999). Doutorando Direto no programa de Pós-Graduação em Oceanografia Física, desde 2000, pelo Instituto Oceanográfico da USP. Membro do Laboratório de Dinâmica Oceânica atua em estudos de circulação de meso-escala através de análise de dados e modelagem numérica. Bolsista-sanduíche CAPES, desde julho de 2003, no Department of Earth, Atmospheric, and Planetary Sciences do Massachusetts Institute of Technology (MIT), nos Estados Unidos, sob orientação da Profa. Dra. Paola Malanotte-Rizzoli e Dr. Markus Jochum.

Ilson Carlos Almeida da Silveira é professor-Doutor do Departamento de Oceanografia Física, Química e Geológica do Instituto Oceanográfico da Universidade de São Paulo, desde 1997. Graduado em Oceanografia pela UERJ, em 1984. Mestre em Oceanografia Física pelo IOUSP, em 1990. Nos EUA, Doutor pela University of New Hampshire, em 1996 e pós-doutoramento em Rhode Island, em 1997. Sua linha de pesquisa compreende a circulação de meso-escala do oceano, voltada para teoria de instabilidade hidrodinâmica e para o estudo da dinâmica das correntes de contorno oeste, principalmente da Corrente do Brasil, através do desenvolvimento e utilização de modelos analíticos, semi-analíticos e numéricos.

• BOURLÈS, B et al. Upper layer currents in the western tropical North Atlantic (1989-1991). Journal of Physical Oceanography, [S.l.], v. 104, n. C1, 13611375, 1999.
• BROWN, W.S. et al. A western tropical atlantic experiment (WESTRAX). Oceanography, [S.l.], v. 5, n. 1, p. 7377, 1992.
• CANE, M. The response of an equatorial ocean to simple wind-stress patterns, Parts I and II. Journal of Marine Research, [S.l.], n. 37, p. 233-299, 1979.
• ______; SARACHIK, E. Forced baroclinic ocean motions, I. Journal of Marine Research, [S.l.], v. 34, n. 4, p. 629665, 1976.
• _____; ______. Forced baroclinic ocean motions, II. Journal of Marine Research, [S.l.], v. 35, n. 2, p. 395432, 1977.
• _____; ______. Forced baroclinic ocean motions, III. Journal of Marine Research, [S.l.], v. 37, p. 355398, 1979.
• EMERY, W. J.; LEE ,W. G.; MAGAARD, L. Geographic and seasonal distributions of Brunt-Väisälä frequency and Rossby Radii in the North Pacific and North Atlantic. Journal of Physical Oceanography, [S.l.], v. 14, p. 294317, 1983.
• HOURY, S. et al. Brunt-Väisälä Frequency and Rossby Radii in the South Atlantic. Journal of Physical Oceanography, [S.l.], v. 17, p. 16191626, 1987.
• LEBLOND, P. H.; MYSAK, L. A. Waves in the ocean [S.l.]: Elsevier, 1978.
• LIGHTHILL, M. J. Dynamics response of the Indian Ocean to the onset of the southwest monsoon. Philosophical Transactions of the Royal Society of London, [S.l.], v. 265, p. 45-92, 1969.
• MCCLEAN, J. L.; KLINCK, J. M. Description and vorticity analysis of the 50-day oscillations in the western tropical Atlantic region of the CME model. Journal of Phisical Oceanography, [S.l.], v. 25, p. 24982517, 1995.
• MCCREARY, J. P. A linear stratified ocean model of the Equatorial Undercurrent. Philosophical Transactions of the Royal Society of London, [S.l.], v. 298, p. 603635, 1980. Series A, Mathematical and Physical Sciences.
• PATTI, E. A geração dos anéis da Corrente Norte do Brasil devido ao forçamento pelo vento em um modelo isopicnal. 2001. 115 f. Tese (Doutorado)-Departamento de Oceanografia Física, Química e Geológica, Instituto Oceanográfico, Universidade de São Paulo, São Paulo, 2001.
• PHILANDER, S. G. H. El Niño La Niña, and the Southern oscillation San Diego: Academic Press, 1989.
• ______; PACANOWSKI, R. C. The Generation of Equatorial Currents. Journal of Geophysical Research, [S.l.], v. 85, n. C2, p. 11231136, 1980.
• PINARDI, N.; ROBINSON, A. R. Dynamics of deep thermocline jets in the POLYMODE region. Journal of Physical Oceanography, [S.l.], v. 17, p. 11631188, 1987.
• SILVEIRA, I. C. A., BROWN, W. S.; FLIERL, G. R.. Dynamics of the North Brazil Current Retroflection from the WESTRAX observations. Journal of Geophysics Research, [S.l.], v. 105, n. C12, p. 28559-28583, 2000.
• SPAIN, P. F.; DORSON, D. L.; ROSSBY, H. T. PEGASUS: a simple, acoustically tracked, velocity profiler. Deep-Sea Research, [S.l.], v. 28A, n. 12, p. 15531567, 1981.
• WILSON, W. D.; JOHNS, E.; MOLINARI R.L. Upper layer circulation in the western tropical North Atlantic Ocean during August 1989. Journal of Geophysics Research, [S.l.], v. 99, n. C11, p. 22, p. 22513-22523, 1994.

Datas de Publicação

• Publicação nesta coleção
  24 Jan 2007
• Data do Fascículo
  Jun 2003

Histórico

• Aceito
  13 Abr 2004
• Recebido
  15 Jul 2003