Integrando a Camada Limite Planetária e os Fluxos de Superfície em Modelos Paramétricos de Ciclones Tropicais: Um Passo em Direção ao Realismo Termodinâmico

Descrição

Nas últimas quatro décadas, a modelagem numérica da atmosfera consolidou-se como um pilar fundamental para a previsão e compreensão de fenômenos meteorológicos complexos, com destaque para os ciclones tropicais. A necessidade de previsões precisas é impulsionada pelos devastadores impactos socioeconômicos associados a esses sistemas, o que exige ferramentas de diagnóstico cada vez mais acuradas para mitigar desastres. Historicamente, os modelos paramétricos, como os propostos por Schloemer (1954) e, particularmente, por Holland (1980), foram essenciais nesse desenvolvimento. Baseados em um número limitado de parâmetros-chave — como a pressão central e o raio de vento máximo —, esses modelos ofereceram representações matemáticas abstratas, porém coerentes e computacionalmente eficientes, dos campos de pressão e vento, estabelecendo uma base sólida para a avaliação de riscos. No entanto, a própria simplicidade que os tornou úteis representa também sua principal limitação, especialmente no que tange à sub-representação das características assimétricas dos ciclones. A crescente demanda por modelos mais complexos revela uma lacuna significativa na literatura: a sub-representação das características assimétricas. Frequentemente simplificados em favor de uma estrutura idealizada e axissimétrica, os componentes assimétricos são, na realidade, cruciais, pois resultam da interação com o cisalhamento vertical do vento, os gradientes de temperatura da superfície do mar e processos dinâmicos internos, como os ciclos de substituição da parede do olho. Essas assimetrias determinam a distribuição espacial dos perigos — como a localização do quadrante com ventos mais intensos, a ocorrência de chuvas torrenciais e o potencial de formação de tornados — influenciando diretamente a trajetória e a intensidade dos ciclones tropicais. O desafio central reside, portanto, em desenvolver modelos que equilibrem a eficiência computacional dos modelos paramétricos com o realismo físico dos modelos numéricos completos (como o WRF), que, apesar de seu detalhamento, apresentam alto custo computacional. Nesse contexto de busca por maior realismo físico, avanços teóricos seminais, como o de Emanuel (1986), revolucionaram a compreensão da termodinâmica dos ciclones tropicais. Emanuel reformulou o ciclone tropical não apenas como um vórtice em equilíbrio de gradiente, mas como uma máquina térmica de Carnot que extrai energia do oceano quente e a converte em energia cinética dos ventos. Esse processo é governado pelo mecanismo de Wind-Induced Surface Heat Exchange (WISHE), um ciclo de retroalimentação positiva no qual ventos mais fortes aumentam a evaporação da superfície do mar, fornecendo mais umidade para a convecção profunda, que, ao liberar calor latente, intensifica a circulação e fortalece ainda mais os ventos. Essa teoria destacou o papel central do equilíbrio do vento térmico e da conservação do momento angular, estabelecendo um limite superior teórico para a intensidade de um ciclone tropical — a chamada Intensidade Potencial (IP) — baseada nas condições ambientais. Inspirada nesses fundamentos, a pesquisa recente tem buscado incorporar mais dessa dinâmica nos modelos paramétricos. Um exemplo notável é o trabalho de Mendes et al. (2022), que aprimorou o modelo de Holland (1980) ao introduzir formulações baseadas na conservação do momento angular absoluto. Essa abordagem captura de forma mais elegante o processo de spin-up do vórtice, descrevendo como as parcelas de ar aceleram tangencialmente à medida que convergem para o centro do ciclone, resultando em perfis de vento mais realistas e superando algumas das limitações puramente empíricas do modelo original. O presente trabalho busca ir além, construindo sobre a base fisicamente mais robusta proposta por Mendes et al. (2022). Propõe-se a inclusão explícita de variáveis físicas cruciais que governam a interface oceano-atmosfera e a estrutura vertical inferior do sistema, elementos que são o verdadeiro motor da máquina térmica ciclônica. Especificamente, incorporam-se: (1) a altura da Camada Limite Planetária (CLP), a camada turbulenta na base da atmosfera através da qual o ar flui para o centro do ciclone; e (2) os fluxos de calor latente (FL) e sensível (FS) na superfície, que quantificam a taxa de extração de energia do oceano. O modelo proposto estrutura-se em torno de vários componentes físicos interligados. Primeiramente, na dinâmica do vórtice e no spin-up inicial, quantifica-se o momento angular absoluto por unidade de massa como a soma do momento angular relativo do furacão e do momento angular planetário. Sob condições ideais sem atrito, o momento angular absoluto conserva-se, permitindo que uma parcela de ar ao convergir para o centro adquira alta velocidade tangencial. Contudo, reconhece-se que o atrito com a superfície oceânica real dissipa energia, reduzindo os ventos observados. À medida que os ventos se intensificam, estabelece-se o equilíbrio do vento gradiente, que incorpora as forças de gradiente de pressão, Coriolis e centrífuga. Em regimes de ventos extremos (como na parede do olho), a força centrífuga domina sobre a de Coriolis, simplificando o equilíbrio para o vento ciclostrófico. Em relação à estrutura termodinâmica de núcleo quente, o modelo captura a anomalia térmica central gerada pela liberação de calor latente e pelo aquecimento adiabático por compressão. Esse gradiente radial de temperatura causa uma inversão do gradiente de pressão com a altitude, explicada pela relação do vento térmico. Utilizando a equação hipsométrica, relaciona-se a diferença de temperatura média entre o núcleo e o ambiente com a diferença de pressão entre a base e o topo da tempestade. Para os perfis de vento e pressão, o modelo adota funções analíticas que representam a queda pronunciada de pressão perto do centro, sendo o raio de vento máximo um parâmetro crítico que define a escala do núcleo intenso. Uma relação empírica, baseada na equação de Bernoulli, conecta a queda de pressão central ao aumento dos ventos máximos tangenciais. A partir da aproximação ciclostrófica e do perfil de pressão, deriva-se a distribuição da velocidade tangencial radial. O modelo também incorpora a circulação secundária. A velocidade radial descreve o influxo em espiral do ar na camada limite, causado pela quebra do equilíbrio de gradiente devido ao atrito com a superfície do mar. Esse influxo transporta o ar quente e úmido da periferia para a parede do olho. Pelo princípio da continuidade, o influxo intensifica-se à medida que o raio diminui, com formulação específica representando esse comportamento. A velocidade vertical, por sua vez, ocorre quando o ar convergente na camada limite atinge a parede do olho e é forçado a ascender violentamente, formando as correntes ascendentes. Dentro do olho predomina a subsidência, enquanto na parede do olho a velocidade vertical é fortemente positiva. Por fim, para o fechamento termodinâmico do modelo, garante-se que a intensidade da anomalia do núcleo quente esteja diretamente ligada à queda de pressão na superfície. Assumindo que a anomalia de pressão no topo do furacão é igual e oposta à da base (simplificação do equilíbrio hidrostático), e combinando com a relação termodinâmica derivada, obtém-se a diferença de temperatura média necessária para sustentar uma dada intensidade. Isso implica que uma queda de pressão maior exige um núcleo correspondentemente mais quente, garantindo consistência entre dinâmica e termodinâmica. Os objetivos da pesquisa são claramente definidos. O objetivo geral consiste em avaliar o impacto quantitativo da incorporação da altura da CLP e dos fluxos de calor latente e sensível no modelo de Mendes et al. (2022), transformando-o em um modelo paramétrico de múltiplas camadas, visando uma representação mais realista da dinâmica de ciclones tropicais, especialmente durante fases de rápida intensificação. Os objetivos específicos incluem: incorporar explicitamente essas variáveis; analisar se melhoram a representação da evolução da intensidade (taxa de queda da pressão central); avaliar se proporcionam simulação mais realista do campo de vento (evolução do raio de vento máximo); comparar os resultados com o modelo original e com dados observacionais de alta qualidade; e verificar se permitem capturar com maior fidelidade a estrutura termodinâmica e assimétrica, considerando modulações por cisalhamento vertical do vento. A expectativa é que o modelo proposto ofereça múltiplas vantagens e aplicações. No diagnóstico melhorado da estrutura e intensidade do ciclone, a abordagem baseada em balanços de força fundamentais (vento gradiente e ciclostrófico) e na conservação do momento angular permite uma representação mais robusta e fisicamente coerente da relação entre pressão e vento, superando abordagens puramente empíricas. Como ferramenta para a pesquisa da intensificação rápida — um dos maiores desafios na previsão de furacões — o modelo articula explicitamente a conexão entre a termodinâmica do núcleo quente e a dinâmica dos ventos, servindo como um "laboratório numérico" de baixo custo computacional para estudos de sensibilidade. Em termos de agilidade operacional, diferentemente dos modelos numéricos completos que exigem supercomputadores, o modelo analítico proposto pode ser resolvido quase instantaneamente em computadores convencionais, permitindo gerar rapidamente estimativas iniciais ou conjuntos de cenários. Na modelização de impactos costeiros, os perfis de vento e pressão gerados servem como forçantes para modelos de maré de tempestade e ondulação, com a base física garantindo forçamentos dinamicamente consistentes. Adicionalmente, o modelo constitui uma excelente ferramenta didática e de treinamento avançado, pois sua estrutura transparente permite visualizar de forma imediata como parâmetros como latitude ou temperatura do núcleo impactam a intensidade e estrutura do furacão. Conclui-se que a inclusão explícita da altura da Camada Limite Planetária e dos fluxos de calor latente e sensível representa um avanço significativo na modelagem paramétrica de ciclones tropicais. Acredita-se que essa parametrização permitirá uma representação mais fiel não apenas da intensidade, mas também da própria estrutura termodinâmica e assimétrica, uma vez que a distribuição espacial dos fluxos e a topografia da CLP são fortemente moduladas por condições assimétricas do ambiente, como o cisalhamento do vento. A hipótese central é que a inclusão explícita do motor termodinâmico conduzirá a uma simulação mais realista da evolução da intensidade e da estrutura do campo de vento, em comparação com o modelo original e com dados observacionais de alta qualidade, preenchendo assim uma lacuna importante na literatura e fornecendo uma ferramenta valiosa tanto para a pesquisa científica quanto para aplicações operacionais na bacia do Atlântico Tropical Norte.