1. Introdução
Existe uma vasta literatura científica acerca da ocorrência de precipitação no Nordeste do Brasil (NEB) (Yamazaki e Rao, 1977; Zhou e Lau, 2001; Carvalho et al., 2002, 2004; Ferreira e Mello, 2005; Souza e Cavalcanti, 2009; Reboita et al., 2010; Oliveira et al., 2013; 2016; Rodrigues et al., 2019; 2020), que é caracterizado por apresentar longos períodos de seca (Marengo et al., 2011; Costa et al., 2015; Dos Santos et al., 2019) e ocorrencias de eventos de precipitação intensa (Costa et al., 2015; Oliveira et al., 2017; Da Silva et al., 2019; Rodrigues et al., 2020). O NEB engloba os Estados do Maranhão, Piauí, Ceará, Rio Grande do Norte, Paraíba, Pernambuco, Alagoas, Sergipe e Bahia, totalizando uma área de aproximadamente 1.558.196 km². A região apresenta alta variabilidade espacial e temporal da distribuição de chuvas (Moura e Shukla, 1981; Zhou e Lau, 2001; Oliveira et al., 2017; Rodrigues et al., 2019), com acumulados anuais inferiores a 500 mm em regiões semiaridas e superiores a 1500 mm em regiões costeiras (Oliveira et al., 2017; Rodrigues et al., 2019). A principal estação chuvosa do NEB ocorre entre março e junho, enquanto o período seco, em grande parte do NEB, ocorre entre agosto e outubro (Rao et al. 1993; Rodrigues et al., 2019).
Dentre os Estados que englobam o NEB, Pernambuco e Rio Grande do Norte (RN) são os que possuem maior extensão de território com clima semiárido (Bsh), respectivamente, cerca de 61,4, e 61,2% (Alvares et al., 2013). Os Relatórios Científicos do IPCC (AR4 e AR5) expõem o semiárido do Brasil como uma das regiões do globo que são mais afetadas e sensíveis às mudanças climáticas atuais e futuras (Marengo, 2014).
O clima semiárido é predominante na região cobrindo desde áreas costeira do RN em cerca de 150 km da região litorânea, maior parte do Planalto da Borborema no estado da Paraíba, onde é considerado um dos lugares mais seco do Brasil (o vale do rio Paraíba) semelhante a Paulo Afonso na Bahia que também foi considerado com essa característica (Alvares et al., 2013). Além dos estados do Pernambuco, Ceará com regiões próximas a costa, Alagoas e Sergipe com 14,9 e 13,5% do território de Bsh. Em contrapartida, apenas o estado do Maranhão não apresenta clima Bsh (Alvares et al., 2013).
Apesar da ocorrência de secas ser uma característica marcante do clima do NEB; diversos sistemas meteorológicos favorecem a precipitação na região ocasionando máximos de precipitação anuais que variam de 300 a 2000 mm (Rodrigues et al., 2019). Entre esses principais sistemas pode-se citar a Zona de Convergência Intertropical (ZCIT) (Zhou e Lau, 2001; Souza e Cavalcanti, 2009), os Vórtices Ciclônicos de Altos Níveis (VCAN) (Kousky e Gan, 1981; Coutinho et al., 2010), os Distúrbios Ondulatórios de Leste (DOL) (Yamazaki e Rao, 1977; Torres e Ferreira, 2011; Gomes et al., 2015), Linhas de Instabilidade (LI) (Ferreira e Mello, 2005), Complexos Convectivos de Mesoescala (CCM) (Souza et al., 1998; Fedorova et al., 2009; Lyra et al., 2015), Sistemas Frontais (SF) (Kousky, 1979; Alves e Kayano, 1991) e a Zona de Convergência do Atlântico Sul (ZCAS) (Chaves e Cavalcanti, 2001; Carvalho et al., 2002; Carvalho et al., 2004).4 32
A ZCIT é o principal sistema que favorece a precipitação na região NEB. Atua entre as estações de verão e outono austral (Zhou e Lau, 2001), principalmente, no setor Norte do NEB (NNEB), em que em anos chuvosos pode atingir até 5° S (Melo et al., 2009). Os VCANs atuam no verão austral com máxima predominância no mês de janeiro (Kousky e Gan, 1981; Coutinho et al., 2010). Os DOL ocorrem durante todo ano (Diedhiou et al., 2010) ao longo da costa NNEB e Leste do NEB (LNEB); contudo são responsáveis por elevados índices pluviométricos nos meses de inverno em regiões do LNEB (Ramos, 1975; Amorim et al., 2019; Rodrigues et al., 2020). Os SF e a ZCAS atuam sobre o setor Sul do NEB (SNEB), em especial, nos meses de novembro a janeiro (Ferreira e Melo, 2005).
A precipitação do NEB é afetada por alguns modos de variabilidade em diferentes escalas de tempo, dos oceanos Pacífico e Atlântico. A influência dos modos de alta frequência no clima do NNEB é bem fundamentada na literatura (Souza e Nobre, 1998; Misra e Zhang, 2007). Basicamente, as alterações na circulação de grande escala leste-oeste (célula de Walker) e norte-sul (célula de Hadley) configuram um ambiente favorável ou desfavorável a atuação dos principais sistemas meteorológicos no NNEB, durante a estação de outono austral (Souza e Ambrizzi, 2002). No entanto, vale lembrar que o RN apresenta máximos de precipitação também no inverno.
Quanto aos modos de variabilidade baixa frequência no Atlântico, um dos existentes no setor Norte (Nobrega et al., 2014) é semelhante à Oscilação Decadal do Pacifico (Newman et al., 2016) com fases quente e fria, chamado de OMA (Oscilação Multidecadal do Atlântico) (Enfield et al., 2001). Durante o período de 1925 a 1960, a OMA apresentava-se na fase positiva e negativa entre 1964 a 1995 (32 anos) (Santos et al., 2016), e após 1995 iniciou-se a fase positiva (Trenberth e Zhand, 2019).
O segundo modo de variabilidade do Atlântico é caracterizado por centros de ação com fases opostas, um na região da Islândia e o outro na região do Atlântico Norte Central (Açores), Rogers (1984). No inverno do hemisfério norte (NH) é mais evidente, mas é possível a ocorrência nas demais estações do ano. Existe um índice para verificar esta oscilação (North Atlantic Oscilation Index), quando o valor é positivo, os centros de alta subtropical do Atlântico Norte e a baixa da Islândia apresentam-se anomalamente mais intensos. Justamente na variabilidade da alta subtropical do Atlântico Norte (ASAN) que se deve à influência em relação à precipitação da região NEB. Namais (1972) foi o trabalho pioneiro entre a precipitação no NEB e o NAO mostrando que anos chuvosos advêm durante a fase positiva do NAO. A fase negativa do NAO prevaleceu da década de 1950 até meados da década de 70. A partir do final da década de 1970 a 1990, NAO positivo predominou.
Especificamente, o RN encontra-se na área de atuação da ZCIT. No entanto, o clima do RN apresenta na faixa leste do RN, máximos de precipitação durante as estação do outono (março, abril e maio – MAM) referente a esse sistema, além da estação do inverno (junho, julho e agosto – JJA) austral. Estudos destacam fenômenos no inverno austral da região do NEB, que se aplica ao RN: a ocorrência de DOL, que se propagam para oeste sobre o Oceano Atlântico (Reboita et al., 2010; Torres e Ferreira, 2011; Santos et al., 2012; Kouadio et al., 2012; Gomes et al., 2015; Hounsou et al., 2015; Amorim et al., 2019). Além do DOL, a influência de brisa, especialmente, no período noturno, referente à formação de uma zona de convergência entre a brisa terrestre e o escoamento em baixos níveis (Kousky, 1979) e a variabilidade da Alta Subtropical do Atlântico Sul (ASAS) (Rao et al., 1993). Rao et al. (1993) sugerem que esse sistema é o principal mecanismo de transporte de vapor de água nos primeiros níveis da atmosfera.
O estado do RN apresenta faixa litorânea ao norte e a leste, e a maior população (877,662 habitantes) habita na capital Natal, que se encontra no litoral oriental do estado e ocupa uma área de 167,264 km² (IBGE, 2014). Recentemente, resultados mostraram evidências que eventos extremos de precipitação (EEP) diários apresentou tendência positiva nos últimos anos (Amorim et al., 2019), entre 1986 a 2016, na cidade de Natal. Alguns desses EEP contribuíram para seguidos desastres naturais, em decorrência a existência de áreas de alto risco natural (Moreira et al., 2014; Macedo et al., 2016; Rodrigues et al., 2020) na capital. Além disso, a população da capital também apresenta-se exposta a frequentes alagamentos devido à obras paradas e/ou incompletas do sistema de drenagem da cidade, que são potencializados na ocorrência dos EEP.
Amorim et al. (2019) destacam para a cidade de Natal, em vários anos de El Niños, a variabilidade anual na frequência de EEP diários apresentou-se diferente. Percebe-se que o RN apresenta variabilidadade da precipitação tanto no espaço quanto no tempo (Rodrigues et al., 2020). Sendo assim, estudos que venham a contribuir com informações fundamentais a ações de adaptações às alterações climáticas são de grande valia. Neste contexto, o objetivo principal deste estudo é investigar tendências climáticas em áreas homogêneas de precipitação do estado do RN, em especial, durante a estação chuvosa, além de correlaciona-las a modos de variabilidade dos oceanos Pacífico e Atlântico tropicais. As hipóteses deste trabalho são descritas pelas diferenças espaciais do comportamento da distribuição de chuvas sobre o RN, assim como a influência da variabilidade dos oceanos Pacífico e Atlântico tropicais em relação ao quantitativo da precipitação no RN.
2. Materiais e Métodos
2.1. Dados
Utilizou-se a quarta versão do produto precipitação mensal disponibilizado pela Universidade de Delaware para analisar o clima do RN. Foram obtidos vinte pontos de grade de espaçamento de 0,5° × 0,5° de latitude por longitude (Fig. 1). O produto apresenta dados interpolados para uma grade regular (Matsuura e Willmott, 2009) e cobre o período de 1900 a 2014. A construção dos dados partiu de múltiplas fontes (Global Historical Climatology Network – GHCN2; Serviço Atmosférico e Ambiental do Canadá; Instituto hidrometeorológico de São Petersburgo-Rússia cedido por Nikolay Shiklomanov; dados do projeto Greenland Climate Network – GC-Net, Steffen et al. (1996); registros das estações meteorológicas automáticas do Projeto Groenlândia disponibilizados por Charles R. Stearns na Universidade de Wisconsin-Madison; dados diários da índia do National Center for Atmospheric Research – NCAR; dados de precipitação do continente Africano, Nicholson (2001); dados mensais de precipitação do continente sul-americano, Webber e Willmott (1998); e do Global Surface Summary of Day – GSOD); no entanto, o GHCN2 é a fonte de dados principal do produto.
Utilizaram-se dados de temperatura da superfície do mar (TSM) mensal do National Oceanic & Atmospheric Administration – Extended Reconstructed Sea Surface Temperature V3b (NOAA_ERSST_V3b) dispostos em grade horizontal regular de 2,0° latitude × 2,0° longitude (Smith et al., 2007) para as áreas dos Oceanos Atlântico e Pacífico, no período de 1900 a 2012.
No Pacífico, as áreas tradicionais do Niño 3 (150° W – 90° W/ 5° S – 5° N) e Niño 3.4 (170° W – 120° W/ 5° S – 5° N) durante os meses de novembro a fevereiro (NDJF) (Trenberth, 1997) e para o Atlântico Tropical Norte (5 – 20°N, 60-30°W) e Sul (0-20° S, 30° W-10°) no trimestre março a maio (MAM).
Também se utilizaram séries temporais dos modos de variabilidade Atlantic Multidecadal Oscillation (AMO) para o período de 1900 a 2014 e North Atlantic Oscillation (NAO) de 1950 a 2014 disponibilizados pelos endereços eletrônicos https://climatedataguide.ucar.edu/climate-data/atlantic-multi-decadal-oscillation-amo e ftp://ftp.cpc.ncep.noaa.gov/wd52dg/data/indices/, respectivamente.
2.2 Métodos
Utilizou-se a análise de agrupamento para identificar áreas homogêneas durante 1900 a 2014. A avaliação do agrupamento realizou-se através do gráfico de Silhueta, desenvolvida por Rousseeuw (1987) e utilizada por diversos autores, Bernard et al. (2013); Santos et al. (2015); Wilcke e Bärring (2016); Mooney et al. (2017); Rahman et al. (2018); Salehnia et al. (2019), em estudos na área de ciências climáticas. Em seguida, investigaram-se as tendências climáticas na quadra chuvosa e possíveis quebras estruturais na série de dados, além de correlações com os índices oceânicos dos Oceanos Pacífico e Atlântico. A correlação cruzada foi utilizada para identificar a defasagem entre os modos AMO e NAO com a precipitação do RN.
2.2.1. Análise de agrupamento (Cluster)
A análise de agrupamento é um método multivariado que tem como objetivo separar grupos de acordo com uma medida de dissimilaridade entre eles (Wilks, 2006). Utilizou-se essa técnica para identificar regiões cujo comportamento da precipitação do RN possa ser considerado homogêneo durante 1900 a 2014. Aplicou-se a técnica aos dados de precipitação mensal (janeiro a dezembro) para vinte pontos de grade do domínio do estado exposto na Fig. 1.
Utilizou-se a distância euclidiana (Eq. 1) como medida de dissimilaridade comumente usadas em aplicações a dados meteorológicos.
em que Xi,k e Xj,k são as séries temporais para cada ponto de grade comparados dentro do grupo k; Wk =1, …, k é um fator de ponderação e os índices i e j representam os pontos de grade do grupo.
Além da medida de dissimilaridade, utilizou-se o método de agrupamento hierárquico de variância mínima de Ward (Eq. 2) definida segundo a Eq. 2:
sendo que G e g representam os grupos de pontos de grade; Xi é cada vetor de pontos de grade da variável K;
Em seguida utilizou-se o gráfico de Silhueta (Rousseeuw, 1987) para averiguar a construção de cada agrupamento (Santos et al., 2015). A Silhueta é a média do valor do contorno de cada observação. Observações agrupadas adequadamente possuem valores próximos de 1 e observações mal agrupadas, valores próximos de −1. O valor da Silhueta do objeto Xi (cada observação) é calculado pela Eq. (3):
em que a(Xi) é o comprimento médio de todas as distâncias dentro do mesmo grupo e b(Xi) é comprimento médio de todas distâncias das observações do grupo mais próximo.
A qualidade total do aglomerado pode ser calculada pelo coeficiente de Silhueta médio pela Eq. (4):
em que N é o número total de pontos de grade.
2.2.2. Análise da tendência
Utilizou-se para análise da tendência linear o teste não-paramétrico Mann-Kendall (Mann, 1945; Kendall, 1975). Aplicou-se o teste para o acumulado de precipitação do período chuvoso. O uso do teste Mann-Kendall (MK) tem sido comum para mostrar tendências hidrológicas (Santos e Santos e Silva, 2013; Costa et al., 2015; Oliveira et al., 2017), além de ser um método adequado no uso de análises de mudanças climáticas (Goossens e Berger, 1986). Utilizou-se o teste Pettitt (Pettitt, 1979) para identificar possível ponto único de mudança no acumulado de precipitação da estação chuvosa durante o período de estudo. Os níveis significância foram adotados de 1% e 5%.
3. Resultados e Discussão
Os resultados partiram da análise de agrupamento sugerindo três áreas homogêneas (Fig. 2) para o estado do RN: Oeste, Seridó (central) e Litoral. Em seguida, aplicaram-se os testes para série temporal da estação chuvosa de cada grupo, além da investigação aos modos de variabilidade dos Oceanos Pacífico e Atlântico.
3.1. Regiões homogêneas
Os resultados da análise de agrupamento são exibidos para duas (Fig. 2a) e três regiões (Fig. 2d) homogêneas, a preferência de três regiões (Fig. 2d) são mais consistentes com a distribuição espacial da climatologia de precipitação, concordando com a análise prévia de Santos e Silva et al. (2012). Além disso, diferenças significativas (valor-p < 0,01) foram obervadas no acumulado de precipitação nos meses da estação chuvosa, entre os grupos encontrados do Semiárido (Fig. 2f).

Figura 2 Resultados da análise de agrupamento para o RN no período de 1900 a 2014 para dois e três grupos (ilustrações superiores e inferiores, respectivamente): (a) – (d) distribuição espacial das áreas homogêneas; (b) – (e) coeficientes de silhueta (qualidade dos agrupamentos); (c) – (f) ciclo anual de precipitação em milímetros, as letras O, S e L representam as regiões Oeste, Seridó e Litoral.
Os resultados da qualidade dos agrupamentos exibem os coeficientes de Silhueta (S(Xi)) de cada grupo diminuíndo ao aumentar a quantidade de grupos. Por exemplo, para dois grupos (Fig. 2b) os valores apresentaram-se acima de 0,60 (S(Xi) = 0,66 no Litoral e S(Xi) = 0,62 no Semiárido). Para três grupos (Fig. 2e) os valores S(Xi)[U1] diminuem para 0,45 e 0,32, respectivamente. O grupo do Litoral apresentou valor S(Xi) acima de 0,6 nas análises de agrupamento (Fig. 2b e 2e).
Valor da qualidade total do
Apesar da análise com dois grupos indicar maiores valores da qualidade total,
O ciclo anual médio dos grupos é exibido na Fig. 3, nos quais se observam a estação chuvosa dos grupos do Semiárido durante os meses de fevereiro a maio apresentando um acumulado médio de 561 mm para o Oeste (Fig. 3a) e 387 mm para o Seridó (Fig. 3b), concordando com Marengo et al. (2011). O grupo Litoral apresentou a estação chuvosa nos meses de abril a julho com acumulado médio de 772 mm (Fig. 3c). A estação seca compreende períodos semelhantes entre os grupos, no entanto, o grupo do Oeste apresenta o período de estiagem entre os meses de agosto a novembro com precipitação média de 20 mm, e os grupos do Seridó e Litoral a partir de setembro até dezembro com acumulados médios de 31 e 105 mm, respectivamente.

Figura 3 Série temporal do acumulado de precipitação da quadra chuvosa das áreas homogêneas de precipitação do RN no período de 1900 a 2014: (a) Oeste, (b) Seridó e (c) Litoral. A linha tracejada azul representa o valor médio da precipitação.
A presença da ZCIT na região NEB é o principal sistema meteorológico responsável pela ocorrência de precipitação da estação chuvosa na região do RN, principalmente, nos grupos do Semiárido (Oeste e Seridó) e nos meses de março e abril do grupo do Litoral. A partir de maio e inverno autral, destaca-se a influência dos distúrbios ondulatórios de leste (Gomes et al., 2015), principalmente, no mês junho no Litoral que apresenta o máximo de precipitação média de 210 mm. Além disso, eventos de precipitação intensa durante o inverno austral foram associados a uma circulação anticiclônica anômala que predomina na região em 850 hPa (Oliveira et al., 2013), onde os ventos dessa circulação convergem com vento do Oceano Atlântico, criando uma área de confluência transportando umidade para o grupo do Litoral.
Durante o período de estiagem da costa norte do LNEB, que envolve o grupo do Litoral foi associado ao predomínio de VCANs (Oliveira et al., 2013). Esses sistemas meteorológicos foram associados eventos de precipitação intensas durante a estação seca da costa atuando a partir de novembro e durante o verão austral.
3.2. Análise de tendência linear
Aplicou-se o teste da tendência para o acumulado de precipitação do período chuvoso de 1900 a 2014 em cada grupo e os resultados são apresentados na Fig. 3. A análise da tendência não mostrou resultados significativos nas regiões do Semiárido (Oeste e Seridó). Em contrapartida, o Litoral sofreu um acréscimo de 130 mm na estação chuvosa com nível de significância de 1%. O teste de Pettitt (P) identificou uma mudança significativa (valor-p < 0,05) na série a partir da década de 60. O acumulado de precipitação da estação chuvosa até meados dessa década não apresentava mudança no valor médio (linha tracejada azul, Fig. 3c), partir de 1964 a estação chuvosa passou a apesentar valores mais elevados refletindo no valor médio de 714 para 844 mm (Fig. 3c).
A precipitação do NEB é afetada por alguns modos de variabilidade em diferentes escalas de tempo, dos oceanos Pacífico e Atlântico. A mudança no acumulado de precipitação da estação chuvosa do Litoral do RN pode ser atribuída algum desses modos de variabilidade. Nóbrega et al. (2014) observaram que a precipitação extrema de regiões que se aproximam da área litorânea no estado da Paraíba, como as regiões da Borborema/Agreste e Litoral confrontadas ao Sertão, destacam-se apenas com os modos do Atlântico.
Este comportamento que Nóbrega et al. (2014) observaram de regiões que se aproximam da área litorânea foi evidenciado nos dados da estação chuvosa do RN (Tabela 1) as tradicionais áreas do Pacífico (Niño3 e Niño3.4, durante os meses de novembro a fevereiro) e Atlântico Tropical Norte e Sul (trimestre março a maio).
Tabela 1 Coeficiente de correlação com significância estatística de 1% entre o período de 1900 a 2012 entre a precipitação da estação chuvosa das áreas do Oeste (FMAM), Seridó (FMAM) e Litoral (AMJJ) e as anomalias de TSM (ATSM) das tradicionais áreas dos Nino3, Nino3.4, Atlântico tropical Norte e Sul.
Anomalias de TSM | Oeste | Seridó | Litoral |
---|---|---|---|
Nino 3 | −0,35 (<0,01) | −0,29 (<0,01) | 0,03 (0,7718) |
Nino 3.4 | −0,35 (<0,01) | −0,28 (<0,01) | 0,01 (0,8803) |
Atlântico Norte | −0,34 (<0,01) | −0,26 (<0,01) | 0,01 (0,9228) |
Atlântico Sul | 0,41 (<0,01) | 0,45 (<0,01) | 0,49 (<0,01) |
Os resultados com modos de alta frequência mostraram correlações com significância estatística (valor-p < 0,01) entre as áreas do Semiárido (Oeste e Seridó) e as tradicionais áreas dos oceanos (Tabela 1). Em contrapartida, na área do Litoral a correlação apresentou significância estatística (valor-p < 0,01) apenas com Atlântico Tropical Sul – ATLS (Tabela 1), resultado semelhante foi observado por Amorim et al., 2014.
Após o grupo do Litoral apresentar correlação significativa (valor-p < 0,01) com ATLS, e estudos mostrarem que esse setor do Atlântico é reconhecido por apresentar uma variabilidade independente de outras regiões oceânicas (Chang et al., 1998; Czaja et al., 2002; Barreiro et al. 2004, 2005; Trzaska et al. 2007), diferentemente do ATLN que é fortemente influenciado pelo Pacífico Equatorial sendo positivamente correlacionado com o índice Niño 3 (Enfield e Mayer 1997; Alexander e Scott 2002; Huang et al., 2005a), investigou-se a tendência linear durante o trimestre MAM para o ATLS. Aplicou-se o teste da tendência para as anomalias de TSM (ATSM) do ATLS durante o período MAM de 1900 a 2012, e os resultados são apresentados na Fig. 4.

Figura 4 Série temporal das ATSM do ATLS para o período de MAM no período de 1900 a 2012. A linha tracejada azul representa o valor médio da anomalia: (a) Período de 1900 a 2012 e (d) Período de 1960 a 2012.
A análise da tendência mostrou um aumento nos valores das ATSM com significância estatística (valor-p < 0,01) para região do ATLS. A região que apresentou correlações significativas com o litoral do RN sofreu um acréscimo de 0,6 °C entre o período de 1900 a 2012 (Fig. 4a). O teste de Pettitt (P) identificou uma mudança significativa (valor-p < 0,01) na série a partir da década de 60 (Fig. 3c). As ATSM do período MAM apresentava valor médio de −0,73 °C (linha tracejada azul, Fig. 4a), partir de 1960 ATSM passou a apresentar valores mais elevados refletindo no valor médio de 0,11 °C.
A mesma análise se repetiu para o período a partir de 1960 a 2012 (Fig. 4b) para o trimestre MAM do ATLS. Os resultados mostram um aumento significativo (valor-p < 0,01) de 0,2 °C das ATSM a partir de meados da década de 80 (Fig. 4b). As ATSM apresentada anteriormente exibem um valor médio de −0,3 °C, a partir de 1985 ATSM passa a apesentar valor médio de 0,1 °C (Fig. 4b).
A influência que os modos de variabilidade dos oceanos apresentam no clima, alguns não ocorrem de maneira simultânea. O resultado da correlação sem defasagem de tempo entre os dados OMA e NAO com a precipitação na escala mensal confirma esse comportamento, apesar de apresentarem significância estatística com valores-p < 0,01, os valores dos coeficientes de correlação apresentam-se abaixo de 0,2 para os três grupos do RN.
A correlação cruzada entre o modo OMA e a precipitação dos grupos do RN exibiram defasagens a partir de um mês até doze meses anteriores. O segundo modo de variabilidade do Atlântico (NAO), a correlação cruzada não mostrou resultados satisfatórios.
O modo OMA apresentou correlações significativas com grupos do Semiárido (Oeste e Seridó) apenas nas defasagens entre dois e seis meses que antecedem a estação chuvosa, ou seja, durante julho a dezembro do ano anterior, as correlações foram −0,38 com o grupo Oeste e −0,44 com Seridó). Apesar do resultado mostrar correlações com valores baixos, a relação negativa exibida com dados da precipitação dos grupos do RN apresentaram significância estatística (valor-p < 0,01). No entanto, o Litoral mostrou correlações de −0,23 entre quatro e seis meses que antecederam o período de AMJJ com nível de significância de 10%. A Fig. 5 ilustra a série temporal mensal do OMA e NAO.

Figura 5 índices: (a) da OMA de janeiro de 1900 a dezembro de 2014 e (b) NAO, de janeiro de 1950 a dezembro de 2014. Fonte: NOAA.
O comportamento do Atlântico, em específico o ATLS, pode ter contribuído para mudança no acumulado médio da precipitação do Litoral do RN. Estudos mais detalhados estãos em andamento, além de outros bancos de dados estão sendo explorados para complementar os resultados discutidos neste estudo.
A fase negativa (1964 a 1995) da OMA, além da tendência de aumento das ATSM do ATLS são indícios para mudança no acumulado de precipitação discutida para o Litoral do RN. Paralelo a esses resultados, o crescimento populacional e a urbanização são fatores que precisam ser analisados, pois suas influências são poucas discutidas na literatura. E a população do RN foi estimada em 3.474,998 habitantes (IBGE, 2014) e a maior população (877,662 habitantes) reside no Litoral, onde se localiza a capital Natal.
4. Considerações Finais
Apresentou-se uma avaliação da precipitação do produto de precipitação da Universidade de Delaware. Essa fonte de dados capturou a variabilidade espacial e temporal da precipitação mensal, destacando-se áreas mais secas e mais chuvosas presentes no Estado do RN.
Os resultados da análise de agrupamento foram consistentes com a literatura sugerindo três grupos para o RN: Oeste, Seridó e Litoral. A diferença entre os grupos do Semiárido (Oeste e Seridó) apresentou-se na intensidade da precipitação durante a estação chuvosa.
Os resultados da análise da tendência mostraram evidências que o acumulado de precipitação na região onde se localiza a capital do RN apresentou mudança. Paralelo a esse resultado, o Atlântico Tropical Sul apresentou um aumento nas anomalias de temperatura da superfície na década em que a mudança no acumulado médio da precipitação do grupo do Litoral foi observada. A importante influência exercida pelo Atlântico Tropical Sul à região litorânea discutida na literatura foi mais uma vez ilustrada nos resultados deste estudo.
Além dos modos de variabilidade presentes no Atlântico, o crescimento populacional também precisa ser explorado para verificar possíveis contribuições às alterações climáticas nas regiões litorânes, regiões estas que apresentam cada vez mais ocupação de áreas de risco, refletindo em mais desastres naturais nunca antes registrados.
Em geral, os resultados obtidos neste estudo fornecerão o conhecimento prévio do comportamento da precipitação sobre o RN. Espera-se que a subdivisão espacial definida nesse estudo, seja considerada por pesquisadores que desejam estudar a precipitação sobre o RN, uma vez que as sub-regiões aqui definidas caracterizam a variabilidade espacial presente no Estado de acordo com a precipitação. Além disso, os resultados deste estudo foram capazes de identificar tendência positiva, estatisticamente significativa (valor-p < 0,01), dos valores anuais de precipitação acumulada sobre o litoral do RN. Este resultado foi associado com a tendência positiva observada sobre as anomalias de TSM no ATLS, que foi estatísticamente significante (valor-p < 0,01) para o período de MAM.
Espera-se que os resultados deste estudo sejam benéficos para amparar a formulação de políticas públicas para minimizar a vulnerabilidade da população do RN, sobretudo aquela que vive em áreas de risco, por meio da prevenção e mitigação de eventos catastróficos que podem ocorrer devido ao excesso ou escassez de chuvas.