Vulnerabilidade ao estresse térmico em assentamentos informais

Vulnerability to heat stress in informal settlements

Resumo

As mudanças climáticas globais e a rápida urbanização aumentam o risco de eventos climáticos extremos, tornando espaços urbanos focos de medidas de adaptação e mitigação. As desigualdades sociais inerentes a muitas cidades contribuem para acentuar fragilidades preexistentes, expondo a população a situações de risco ambiental e de saúde, por exemplo, como o estresse térmico. Em vista disso, o objetivo deste trabalho é verificar a vulnerabilidade ao estresse térmico no espaço aberto de um assentamento informal durante a ocorrência de onda de calor. Utilizou-se um índice multidimensional para avaliação da vulnerabilidade ao estresse térmico. Para tanto, realizou-se um levantamento de campo na favela de Paraisópolis, São Paulo, para coleta de dados meteorológicos e aplicação de questionários aos pedestres locais. Os resultados indicaram que o índice PET (temperatura fisiológica equivalente) subestimou a real avaliação dos pedestres em relação à sensação térmica. Considerando o índice multidimensional, aproximadamente um quarto da amostra encontra-se em alta ou muito alta situação de vulnerabilidade, sendo que os indicadores individuais que contribuíram para maior vulnerabilidade foram baixa escolaridade, idade avançada e presença de diabetes e/ou hipertensão. Dessa forma, há indicação de que populações residentes em favelas estão sujeitas ao estresse térmico e expostas a riscos à saúde.

Palavras-chave
Mudanças climáticas; Eventos climáticos extremos; Urbanização

Abstract

Global climate change and rapid urbanization increase the risk of extreme weather events and turn urban spaces into focal points for adaptation and mitigation measures. The social inequalities inherent in many cities exacerbate pre-existing vulnerabilities and expose the local population to environmental and health risk situations, for example, such as heat stress. Given this, the objective of this work is to assess the vulnerability to thermal stress in the open spaces of informal settlements during the occurrence of a heat wave. A multidimensional index was used to evaluate this vulnerability. To achieve this, a field survey was carried out in the Paraisópolis slum, São Paulo, to collect meteorological data and administer questionnaires to pedestrians. The results indicated that the PET (physiological equivalent temperature) index underestimated pedestrians’ actual thermal sensation. Considering the multidimensional index, approximately a quarter of the sample was found to be in a high or very high vulnerability situation, with the key individual factors contributing to the higher vulnerability being low education, advanced age, and the presence of diabetes and hypertension. Therefore, the findings suggest that populations residing in slums are susceptible to thermal stress and exposed to significant health risks.

Keywords
Climate change; Extreme weather events; Urbanization

Introdução

As mudanças climáticas globais e a rápida urbanização aumentam o risco de eventos climáticos extremos e desastres naturais nas cidades (IPCC, 2021), de modo que os espaços urbanos se tornaram foco de políticas de mitigação e adaptação (IPCC, 2022). As desigualdades sociais, inerentes a muitas regiões metropolitanas como São Paulo, tornam as cidades particularmente vulneráveis aos problemas de mudanças climáticas, sendo que fatores socioeconômicos, ambientais e geográficos convergem para acentuar as fragilidades preexistentes (Tarifa; Armani, 2001; Hardoy; Pandiella, 2009; Darela Filho et al., 2016; Di Giulio et al., 2017).

A diversidade das áreas urbanas propicia a formação de microclimas locais, que podem resultar em variações significativas de temperatura e ocasionar níveis distintos de exposição ao calor (Oke et al., 2017; GHHIN, 2020). Nesse sentido, o estresse térmico é entendido como a perturbação do mecanismo de termorregulação humana. No caso de estresse térmico por calor, o corpo produz mais calor do que dissipa e acarreta aumento da temperatura corporal e da frequência cardíaca (Cramer; Jay, 2016). Este desequilíbrio de energia interna no corpo humano pode resultar em desconforto, desmaios e até aumento do risco de morbidade e mortalidade (Di Napoli; Pappenberger; Cloke, 2019). No Brasil, a exposição a temperaturas extremamente altas tem sido associada ao aumento no risco de morte e hospitalização (Geirinhas et al., 2019; Libonati et al., 2022).

Os impactos do estresse térmico na saúde humana são, então, determinados não apenas por fatores climáticos, mas também por outras condições como local de residência e características pessoais, a saber, idade, gênero, peso corporal, condição física, doenças preexistentes, aclimatação, condições financeiras e condições de acesso a serviços de saúde (McGregor et al., 2015; WMO; WHO, 2022). A avaliação da vulnerabilidade a partir de um conjunto de indicadores permite a integração de dados quantitativos e qualitativos para compreender quais indivíduos são mais vulneráveis e de que maneira essa vulnerabilidade se manifesta. Assim, o conceito de vulnerabilidade é entendido como uma combinação entre exposição, sensibilidade e capacidade adaptativa (IPCC, 2014), sendo, portanto, uma noção fluida e dinâmica, que não será comum a todos os cidadãos e/ou residências de uma cidade.

O componente de exposição refere-se à intensidade, duração e frequência de exposição ao calor (Rathi et al., 2021; WMO; WHO, 2022). A sensibilidade, por sua vez, diz respeito a quão bem um indivíduo consegue lidar com o aumento de calor ou até que ponto a exposição o afetará fisicamente (Mallen; Stone; Lanza, 2019). Por fim, a capacidade adaptativa consiste no potencial de mudar para um estado mais desejável frente aos impactos ou riscos às mudanças climáticas (Eakin; Lemos; Nelson, 2014), evidenciando o nível de preparação para antecipar ou responder a eventos de perigo (IPCC, 2014; Sharma; Ravindranath, 2019). Desse modo, percebe-se que um aumento na exposição e sensibilidade acarreta aumento na vulnerabilidade, enquanto o acréscimo na capacidade adaptativa promove redução da vulnerabilidade (Wilhelmi; Hayden, 2010).

Uma vez que é esperado o aumento da frequência e intensidade de riscos associados ao calor (Singh et al., 2019), é fundamental o mapeamento de perigo a estresse térmico para identificar e abordar questões de desigualdade social (GHHIN, 2020), assim como para aprimorar a tomada de decisões e apoiar iniciativas locais mais eficazes (WMO; WHO, 2022). No contexto da urbanização, os assentamentos informais caracterizam-se como espaços que carecem de serviços básicos e infraestrutura urbana. As habitações geralmente não se adequam a regulamentos de planejamento e construção e situam-se em áreas geografica e ambientalmente perigosas. Além disso, também há insegurança de posse em relação à terra e/ou moradias. Nesse sentido, as favelas se enquadram como a forma mais carente e precária de assentamentos informais. As populações residentes nesses espaços estão mais suscetíveis à exclusão social, espacial e econômica, com menor acesso aos benefícios e oportunidades disponíveis nas demais regiões urbanas (ONU, 2015a, 2015b).

Em vista disso, o estudo da percepção de desconforto ambiental em assentamentos informais permite uma análise mais abrangente da vulnerabilidade local. Ao investigar a resposta dos residentes ao microclima da comunidade, pode-se estabelecer uma melhor compreensão sobre como espaços que carecem de infraestrutura adequada, serviços básicos e planejamento urbano resiliente estão suscetíveis aos desafios climáticos iminentes (Sobral, 2005; Nobre et al., 2011). Nesse sentido, o objetivo deste trabalho é verificar a vulnerabilidade ao estresse térmico nos espaços abertos de assentamentos informais durante a ocorrência de onda de calor.

Materiais e métodos

A fim de atingir os objetivos deste estudo, adotou-se o método empírico dedutivo, partindo do pressuposto de que o observado durante a coleta de dados ocorra em escala representativa da favela de Paraisópolis. A metodologia empregada consistiu na coleta de dados em campo por meio de medições meteorológicas e aplicação de questionários, além da obtenção de dados secundários. Em seguida, utilizaram-se essas informações no cálculo de um índice de conforto térmico e no emprego de um índice de vulnerabilidade ao estresse térmico. Essas etapas serão detalhadas nas seções subsequentes.

Índice de vulnerabilidade ao estresse térmico

Diversas pesquisas dedicam-se à agregação de indicadores para criação de índices multidimensionais que avaliam a vulnerabilidade ao calor de modo geral (Wolf; McGregor, 2013; Lapola et al., 2019; Maragno; Fontana; Musco, 2020; Rathi et al., 2021). A partir de uma metodologia que considera características da população e do ambiente (Aubrecht; Özceylan, 2013), utilizaram-se informações demográficas (indicador 1: idade, indicador 2: doenças), socioeconômicas (indicador 3: escolaridade, indicador 4: índice de desenvolvimento humano (IDH)), comportamentais (indicador 5: tempo gasto sob radiação solar direta diariamente), de conforto térmico (indicador 6: de temperatura fisiológica equivalente – PET, do termo em inglês “physiological equivalent temperature” e percepção individual de sensação e conforto) e do ambiente (indicador 7: cobertura vegetal) para composição do índice de vulnerabilidade ao estresse térmico.

Estes indicadores podem ainda ser classificados nas dimensões exposição, sensibilidade e capacidade adaptativa (Quadro 1). É importante ressaltar que a metodologia de referência considera ainda outros indicadores de bases secundárias, isto é, dados já disponíveis em outras fontes, como no recenseamento demográfico. No âmbito desta pesquisa, entretanto, houve a predominância de dados primários, ou seja, aqueles coletados diretamente em campo. Todavia, visto que não há ponderação dos indicadores no cálculo do índice, o método pôde ser empregado em caráter exploratório e experimental. Além disso, originalmente, a metodologia é empregada em escala regional, visando comparar a vulnerabilidade de locais distintos, mas neste estudo aplicou-se um ensaio em escala local.

Tendo em vista as variáveis selecionadas, é fundamental compreender a contribuição de cada uma na vulnerabilidade ao estresse térmico. A idade é um fator importante de sensibilidade, uma vez que crianças e idosos são fisiologicamente mais suscetíveis aos riscos relacionados ao calor (McGregor et al., 2015; Singh et al., 2019; WMO; WHO, 2022). A presença de doenças, como diabetes e hipertensão, também faz com que os indivíduos sejam mais afetados pela exposição prolongada ao calor extremo (Patz et al., 2005; Guo et al., 2017; Dong et al., 2020; WMO; WHO, 2022).

Quadro 1
Indicadores que compõem o índice de vulnerabilidade ao estresse térmico

Dimensão	Indicador	Medida
Sensibilidade	1. Idade	Número, em anos, da idade
	2. Doenças (diabete e hipertensão)	Variável dummy, assume 1 no caso de presença de alguma das doenças e 0 caso contrário
	3. Escolaridade	Categorização nos grupos: 0=Analfabeto; 1=Ensino fundamental; 2=Ensino médio; 3=Ensino superior; 4=Pós-graduação
	4. IDH	Assume valor único para a amostra de acordo com o valor para Paraisópolis (0,361, Atlas Brasil, 2023)
Exposição	5. Tempo gasto sob radiação solar direta diariamente	Número, em horas, do tempo gasto
	6. Conforto térmico (PET e percepção individual de sensação e conforto)	Média entre PET e percepção individual de sensação e conforto. O PET assume valor único para a amostra de acordo com a maior frequência. A percepção individual de sensação é categorizada nos grupos: 0=Neutralidade; 1=Pouco calor; 2=Calor; 3=Muito calor. A percepção individual de conforto é categorizada nos grupos: 0=Confortável; 1=Pouco desconfortável; 2=Desconfortável; 3=Muito desconfortável
Capacidade adaptativa	7. Cobertura vegetal	Variável dummy, assume 1 no caso de ausência de vegetação na favela e 0 caso contrário. Assume valor único para a amostra de acordo com a situação de Paraisópolis (1)

Já o indicador de escolaridade aponta que pessoas com baixo nível educacional são propensas a maior vulnerabilidade ao estresse térmico no que diz respeito à compreensão dos riscos e alertas, bem como às relações do ambiente físico e a exposição ao calor (Harlan et al., 2006; Dong et al., 2020). A educação também está fortemente relacionada a outros fatores socioeconômicos, como o IDH. Este, por sua vez, indica a qualidade de vida das pessoas e o acesso a recursos materiais e informações. Desse modo, locais com menores valores tendem a apresentar realidades nas quais a exposição ao calor afetará significativamente os moradores, como condições habitacionais inadequadas, por exemplo (Patz et al., 2005; Guo et al., 2017; GHHIN, 2020).

Tendo em vista a dimensão exposição, o indicador de conforto térmico, composto pelo índice PET, calculado em função dos dados meteorológicos medidos durante o levantamento de campo, e a percepção individual, obtida por meio das respostas ao questionário. A média entre estes dois indicadores é uma forma de ponderar o índice objetivo e a percepção subjetiva. A variável comportamental de tempo gasto sob radiação solar direta contribui para vulnerabilidade, visto que quanto maior este tempo, maiores os riscos relacionados a queimaduras, insolação e desidratação (Sobral, 2005; Nobre et al., 2011). Por fim, a presença significativa de vegetação no espaço de estudo é um dos elementos capazes de contribuir para a redução dos efeitos do calor, principalmente devido ao sombreamento e ao maior fluxo de calor latente (Mills, 2007; Oke et al., 2017).

Após definição e seleção dos indicadores, normalizaram-se os dados para os valores obtidos ficarem entre 0 e 1 (Equação 1 para relação funcional positiva, e Equação 2 para relação funcional negativa). Algumas abordagens para atribuição de pesos aos indicadores podem incluir julgamento de especialistas (Özceylan; Coskun, 2012), técnicas estatísticas de multivariáveis, como análise de componentes principais (Cutter; Boruff; Shirley, 2003; Conlon et al., 2020) ou abordagens aditivas, por exemplo. Nesse sentido, consideraram-se pesos iguais para os indicadores, de modo que estes tenham igual importância no índice de vulnerabilidade, seguindo a abordagem aditiva (Bao; Li; Yu, 2015). Esta alternativa é comumente utilizada na composição de indicadores, visto que contribui para evitar subjetividade adicional no desenvolvimento do índice (UNEP, 2002; Vincent, 2004; Wilhelmi; Hayden, 2010).

Agregaram-se, por média simples, os valores dos indicadores, de modo que o índice de vulnerabilidade resultou num valor no intervalo entre 0 e 1, em que 1 representa máxima vulnerabilidade ao estresse térmico. Finalmente, dividiram-se os resultados do índice em cinco classes de vulnerabilidade (muito baixa, baixa, média, alta, muito alta) a partir do método de quebras naturais de Jenks (Jenks, 1967).

$d_i=\frac{A_i-m_i}{M_i-m_i}$ (Eq. 1)

$d_i=\frac{M_i-A_i}{M_i-m_i}$ (Eq. 2)

Em que:

d_i é o valor normalizado da variável;

A_i é o valor mensurado da variável;

M_i é o valor máximo da variável; e

m_i é o valor mínimo da variável.

Levantamento de dados meteorológicos

Realizou-se a pesquisa com o método empírico dedutivo, a partir de um trabalho de campo executado no dia 22 de setembro de 2023, entre 09h00 e 16h00, em uma área de grande circulação de pessoas, para medição das seguintes variáveis meteorológicas: temperatura do ar, umidade relativa do ar, direção do vento, velocidade do vento e temperatura de globo cinza. A escolha do dia relaciona-se com a atuação de uma onda de calor em grande parte do Brasil, incluindo o município de São Paulo. Conforme esperado para esse tipo de sistema, registrou-se um dia com baixa cobertura de nuvens e elevada incidência de radiação solar direta.

Utilizou-se um termo higrômetro Hobo Pro v2 com precisão de 0,21 °C para temperatura do ar e de 2,5% para umidade relativa do ar, sendo ambos os sensores localizados no interior de um abrigo de coloração branca e programados para coleta dos dados a cada 5 minutos. Instalaram-se esses equipamentos a um metro do chão. Dadas as condições específicas da área de estudo, não foi possível proceder de modo exato conforme a norma da Agência Nacional de Águas (2016), devido a limitações dos equipamentos e à não existência de área com vegetação. Entretanto, as instalações e medições seguiram boas práticas, além das especificações possíveis. Ademais, mensurou-se, a cada hora, a velocidade do vento com um anemômetro de hélice e a direção do vento com o auxílio de uma bússola.

A partir dos dados de temperatura de globo cinza, calculou-se a temperatura radiante média (TRM) (Aparicio et al., 2016) (Equação 3).

$TRM={[{(T_g+273)}^4+\frac{0,25\ast {10}^8}{{\epsilon }_g}\ast {\left(\frac{|T_g-T_a|}{D}\right)}^{\frac{1}{4}}\ast (T_g-T_a)]}^{\frac{1}{4}}-273$ (Eq. 3)

Em que:

T_g é a temperatura de globo;

é a emissividade do termômetro de globo;

T_a é a temperatura do ar; e

D é o diâmetro do globo.

Em seguida, considerando a temperatura do ar, umidade relativa do ar, velocidade do vento e temperatura radiante média, aplicou-se o índice de temperatura fisiológica equivalente (PET) para avaliação do conforto térmico no espaço aberto. Para cálculo deste índice, que tem como base o balanço energético do ser humano, utilizou-se o software RayMan, considerando os seguintes parâmetros: sexo masculino, metabolismo 100 W, altura 1,65 m, peso 70 kg, resistência térmica da roupa 0,5 clo. Justifica-se o uso de tais parâmetros a partir dos dados preponderantes da amostra analisada. Além disso, o valor para o metabolismo representa a condição de atividade aeróbica leve, como a caminhada, situação em que a maioria das pessoas que responderam ao questionário encontrava-se (Pereira; Souza Junior, 2010).

Aplicação de questionários

Realizou-se uma pesquisa por meio de questionários com os pedestres do local de levantamento. Aplicou-se a entrevista entre 09h00 e 16h00 a partir de abordagens às pessoas que transitavam no local e que concordaram em responder às perguntas propostas. Com essa etapa, coletaram-se informações para caracterização dos residentes da favela de Paraisópolis, bem como para avaliar a percepção dos pedestres quanto à sensação térmica e conforto térmico do local. Além disso, adquiriram-se também dados necessários para o desenvolvimento do índice de vulnerabilidade. Nesse sentido, elaborou-se um questionário com 13 perguntas (Figura 1), orientado pela norma ISO 10551 (ISO, 2019) e submetido na Plataforma Brasil (Parecer número 7.474.894).

Figura 1
Questionário aplicado aos transeuntes do local de levantamento

Caracterização da área de estudo

Realizou-se o estudo em Paraisópolis (Figuras 2 e 3), a maior favela do município de São Paulo. A comunidade localiza-se no bairro nobre Morumbi, no distrito de Vila Andrade, evidenciando um cenário de intensa desigualdade social e segregação espacial. A população da favela é de aproximadamente 58 mil habitantes, sendo estes formados por maioria preta e parda e com aproximadamente 40% de pessoas pertencentes à classe D/E (Instituto Favela Diz, 2022; IBGE, 2022).

Figura 2
Mapa de localização da favela de Paraisópolis e do levantamento de campo

Figura 3
Entorno da área de estudo

O espaço é caracterizado por um aglomerado de habitações em situação precária de construção, que, em geral, são de alvenaria, sem revestimento, sem isolamento térmico, mal ventiladas e com uma parcela de edificações localizada ainda em área de risco de alagamento. A topografia da região é irregular e o sistema viário interno é configurado por vielas estreitas, escadarias, becos sem saídas e com poucas vias adequadas ao tráfego de automóveis. Essa configuração urbana contribui para acentuar as fragilidades já existentes e expor a comunidade a riscos e vulnerabilidades de diversas ordens (Silva; Ribeiro, 2006; Katz, 2008; Alessi, 2009; Castilho, 2013; Pizarro, 2014; Maziviero; Silva, 2018).

Análise dos resultados

Caracterização da amostra

A partir da aplicação do questionário aos transeuntes do local de levantamento, coletou-se a entrevista de 62 pessoas. Nesse ponto, é importante destacar que esta amostra é pequena em relação à população da favela de Paraisópolis. Este total de entrevistas decorre de limitações de tempo, equipamento e pessoal para a realização da pesquisa. Todavia, ressalta-se que, apesar dos limites de representatividade, o estudo pôde prosseguir e ser empregado em caráter exploratório e experimental.

Observou-se distribuição bastante próxima da amostra em relação aos sexos feminino e masculino. Constatou-se que a maioria dos entrevistados é jovem, com idade entre 20 e 29 anos, sendo a menor idade registrada 15 anos e a maior 78 anos. Nesse sentido, notou-se a presença de 13% de idosos (acima de 60 anos) na amostra e a ausência do grupo de risco de crianças. Tendo em vista a altura dos entrevistados, verificou-se a concentração entre 1,55 m e 1,74 m, sendo a média de 1,65 m (±0,09 m). Já no peso da amostra, identificaram-se distribuições similares nas faixas de 55 a 64 kg, 65 a 74 kg e 75 a 84 kg, sendo que estes intervalos compreendem 78% do total das pessoas (Figura 4). Devido ao dia de elevadas temperaturas e alta incidência de radiação solar, característico de um cenário de onda de calor, a vestimenta predominante das pessoas entrevistadas foi de 0,5 clo.

Figura 4
Caracterização da amostra em relação a sexo, idade, altura e peso

Considerando a escolaridade da amostra, observou-se que a maioria possui ensino médio completo, seguido por ensino fundamental completo. Nesta variável, é importante destacar que cerca de 6% dos entrevistados são analfabetos. Em relação à ocorrência das doenças diabete e hipertensão, verificou-se que estas atingem aproximadamente um quarto das pessoas. Já no que diz respeito ao tempo gasto diariamente sob luz solar direta, grande parte das pessoas passa menos de 30 minutos ou entre 1 e 2 horas, enquanto 10% da amostra expõe-se por aproximadamente 8 horas. Além disso, a maioria dos entrevistados é residente de Paraisópolis há pelo menos 10 anos, sendo que 88% dos transeuntes estão na comunidade pelo tempo mínimo de 2 anos (Figura 5), caracterizando uma amostra bastante aclimatada ao local.

Figura 5
Caracterização da amostra em relação à escolaridade, presença de doenças, tempo gasto sob luz solar direta diariamente e tempo de residência na favela de Paraisópolis

Meteorologia e conforto térmico

A partir do levantamento realizado entre 09h00 e 16h00 do dia 22 de setembro de 2023, verificou-se variação da temperatura do ar entre 22,7 °C, valor obtido no início das medições, e 32,6 °C, valor registrado às 14h25min. A temperatura média do período foi de 31,2 °C (±2,1 °C), com relativa estabilidade ao entorno de 32,0 °C entre 12h00 e 16h00. Para a umidade relativa do ar, notou-se valor máximo de 72%, às 09h00, e valor mínimo de 34%, às 11h40min. A umidade relativa média do intervalo foi de 41% (±7%). Em relação à temperatura radiante média, registrou-se o menor valor no início do levantamento, 26,5 °C, e o maior valor às 11h25min, 52,2 °C (Figura 6).

Figura 6
Curvas de temperatura do ar, temperatura radiante média, temperatura fisiológica equivalente e umidade relativa do ar

Considerando a calibração para o município de São Paulo das faixas de temperatura do índice PET, tem-se o intervalo de 18,0 °C a 26,0 °C, correspondendo à neutralidade, de 26,0 °C a 34,0 °C a pouco calor, de 34,0 °C a 43,0 °C a calor e acima de 43 °C a muito calor (Monteiro, 2018). Tendo em vista os resultados obtidos para o PET no período do levantamento, verificaram-se 2% das ocorrências classificadas como neutralidade, 41% como pouco calor e 57% como calor (Figura 6). Nesse sentido, registrou-se o maior valor às 11h25min, 41,4 °C, e o menor valor às 09h00, 20,2 °C, em concordância ao comportamento da TRM. É importante salientar que, durante o período da manhã, o termômetro de globo cinza ficou exposto à radiação solar direta. Esse fato contribui para explicar os valores elevados de TRM e, consequentemente, do PET durante esse período.

Além dos dados medidos e calculados, tem-se também a avaliação dos pedestres entrevistados quanto à sensação térmica e ao conforto térmico (Figura 6). Nesse sentido, observou-se a distribuição das respostas relativas a calor, muito calor, desconfortável e muito desconfortável ao longo de todo o período de estudo. Notaram-se ainda respostas para sensação térmica de neutralidade concentradas no início da manhã e respostas pontuais para pouco calor no final da manhã.

Comparando os resultados dos questionários e dos dados calculados de PET (Figura 7), verificou-se que, enquanto nenhum valor de PET se enquadra na faixa de muito calor, a maioria dos entrevistados (66%), na verdade, classificou a sensação térmica como muito calor. Essas respostas foram registradas, sobretudo, entre 13h00 e 15h00. Apenas 37% da amostra residente em Paraisópolis caracterizou essa situação como muito desconfortável e 29% como desconfortável. Além disso, em relação à tolerância às condições do tempo, 18% dos entrevistados consideraram o ambiente intolerável. Já para a preferência, 61% da amostra gostaria que o tempo estivesse um pouco mais frio, apesar da sensação de muito calor.

Figura 7
Frequência das respostas dos pedestres para sensação térmica, conforto térmico, tolerância e preferência, além da temperatura fisiológica equivalente segundo classificação de Monteiro (2018)

Índice de vulnerabilidade ao estresse térmico

A partir do método de quebras naturais de Jenks, que determina os valores que melhor dividem os dados em grupos homogêneos, agregaram-se os resultados do índice de vulnerabilidade para verificação da vulnerabilidade ao estresse térmico em Paraisópolis em 5 classes. Nesse sentido, constatou-se predominância de pessoas em vulnerabilidade muito baixa e baixa (Figura 8). Entretanto, notou-se que aproximadamente um quarto da amostra encontra-se em situação de vulnerabilidade alta ou muito alta. Para estas pessoas, observou-se que os indicadores que contribuíram para maior vulnerabilidade foram baixa escolaridade, idade elevada e presença de diabetes e/ou hipertensão. A média do índice para a amostra é de 0,44 (±0,11), sendo o menor valor 0,30 e o maior 0,74.

Figura 8
Frequência dos intervalos do índice de vulnerabilidade ao estresse térmico

Discussão

Durante o dia da coleta de dados, 22 de setembro de 2023, grande parte do Brasil, incluindo o município de São Paulo, estava sob o efeito de uma onda de calor (INMET, 2023a). Este evento extremo ocorre predominantemente nas estações de primavera e verão e tem se tornado mais duradouro e intenso nos últimos anos (Bitencourt et al., 2016). Nesse contexto, tendências positivas no número de ondas de calor em cidades brasileiras, como São Paulo, por exemplo, foram associadas ao efeito da urbanização (Geirinhas et al., 2017).

Somado a este cenário, o mês de setembro na capital paulista foi o setembro mais quente considerando toda a série histórica da estação convencional de referência climatológica (INMET, 2023b). Como consequência desse aquecimento urbano, a Secretaria de Estado da Saúde de São Paulo (SES-SP) registrou, nos primeiros sete meses de 2023, um aumento de 102% nos atendimentos relacionados à exposição ao calor com relação ao mesmo período em 2022. Os grupos mais afetados compreendem pessoas com mais de 60 anos, crianças com menos de 4 anos e indivíduos com deficiências cognitivas (SP, 2023).

Tendo em vista as medições realizadas na favela de Paraisópolis, as elevadas temperaturas registradas são coerentes a este cenário generalizado de aumento do calor, representando risco à saúde da população. Nesse sentido, a partir da comparação entre os resultados obtidos para o PET e as respostas dos pedestres para a sensação térmica, constatou-se que os transeuntes apresentaram uma percepção de maior calor do que o índice conseguiu demonstrar. Este fato é observado também em pesquisas semelhantes, em que índices tendem a subestimar a real percepção das pessoas e não demonstrar uma alta taxa de acerto em relação à avaliação da população (Gobo, 2017). Possíveis medidas mitigadoras de estresse térmico em favelas geralmente estão relacionadas aos materiais de construção das residências, especialmente telhados (Nutkiewicz et al., 2022).

Além disso, apesar de aproximadamente dois terços dos entrevistados avaliarem a sensação térmica como de muito calor, apenas cerca de um terço classificou essa sensação como muito desconfortável. Estes resultados evidenciam que as condições para o conforto térmico dependem de uma série de fatores além do clima local, como adaptação cultural e expectativas pessoais, por exemplo (Cramer; Jay, 2016). Dessa forma, ressalta-se a complexidade de estudos de conforto térmico em ambientes externos, visto que essa análise engloba diversos aspectos subjetivos.

Em relação à análise de vulnerabilidade ao estresse térmico, os resultados apontam para que Paraisópolis, já caracterizado como um espaço frágil em termos socioeconômicos, também se encontra suscetível aos impactos decorrentes de alterações climáticas. Os indicadores considerados ressaltam que componentes individuais como baixa escolaridade, idade avançada e presença de certas doenças contribuem para expor os indivíduos a maior vulnerabilidade ao estresse térmico. Além disso, indicadores comuns a toda a amostra acentuam esse cenário de vulnerabilidade, como as condições meteorológicas do levantamento de campo, a ausência de vegetação na comunidade (Shinzato; Duarte, 2018) e o valor médio de IDH (Atlas Brasil, 2023).

No contexto do índice de vulnerabilidade utilizado, é importante ainda ressaltar que a investigação de outros indicadores é relevante para a melhor compreensão da realidade local. A análise de mais doenças preexistentes, principalmente cardiovasculares, respiratórias e psiquiátricas, permite rastrear indivíduos mais suscetíveis a riscos relacionados ao calor, visto que se constatou este indicador como contribuinte para maior vulnerabilidade. A inclusão de mais características do espaço urbano no exame da comunidade também é fundamental para ampliar a pesquisa. As construções precárias, por exemplo, com materiais e técnicas que dificultam o isolamento térmico, resultam em menor inércia térmica e tornam os residentes expostos a maiores variações de temperatura e a condições climáticas extremas. A resiliência do ambiente construído a partir do planejamento urbano é um ponto fundamental no combate à vulnerabilidade das grandes cidades (Bortoli; Vila, 2020; Carvalho et al., 2020).

Ademais, é válido enfatizar a existência de diversas metodologias para análise de vulnerabilidade ao calor de modo geral na literatura. Esta pluralidade permite a aplicação dos métodos em diversos cenários, considerando as especificidades locais e a possibilidade de acesso e/ou levantamento dos dados demandados. Entretanto, a variedade de metodologias dificulta a comparação dos resultados entre pesquisas distintas. Diante disso, percebe-se a necessidade de utilização do mesmo índice de vulnerabilidade em locais com características diferentes, a fim de permitir comparações e análises em escalas regionais.

Conclusão

Este estudo destaca a urgência de abordar a vulnerabilidade ao estresse térmico nos espaços urbanos. A diversidade das cidades resulta em microclimas locais que precisam ser investigados com atenção a todas as suas especificidades. As desigualdades sociais deste cenário precisam ser consideradas no contexto das mudanças climáticas globais, uma vez que eventos climáticos extremos tendem a acentuar as fragilidades preexistentes em assentamentos informais, como a favela de Paraisópolis.

Evidenciou-se a importância de considerar a percepção das pessoas na análise de conforto térmico, uma vez que o PET subestimou a real avaliação dos pedestres acerca do dia de análise. Todavia, como cada índice possui limitações, recomenda-se, em pesquisas futuras, a comparação do PET com outros índices para ambiente externo, como o UTCI (“universal thermal climate index”). Além disso, a partir da metodologia empregada, constatou-se que um quarto da amostra se encontra em alta ou muito alta vulnerabilidade ao estresse térmico, sendo os principais indicadores individuais agravantes a idade avançada, baixa escolaridade e presença de diabetes e/ou hipertensão.

Nesse sentido, aponta-se a amostra reduzida obtida durante o levantamento de campo como uma limitação da pesquisa, sendo necessários estudos futuros com amostras mais abrangentes. Além disso, como não se investigou o tempo que os pedestres estavam no ambiente externo no momento da entrevista, levanta-se a hipótese da influência desse período de exposição na tolerância às condições meteorológicas. Assim, é importante a inclusão de outras perguntas no questionário, que permitam aprofundar a compreensão da relação entre sensação térmica, percepção de conforto térmico e expectativa. Ademais, também recomenda-se a adoção de condições similares durante todo o período de levantamento de campo, evitando situações de sombreamento, por exemplo, a fim de minimizar a influência das variações diurnas nas medições.

Em linhas gerais, o estudo contribui para a ampliação do conhecimento acerca da favela de Paraisópolis e consequentemente contribui para amparar o desenvolvimento de medidas de mitigação e adaptação mais eficazes. Considera-se que os resultados obtidos são uma indicação de que a população residente em comunidades está sujeita ao estresse térmico e exposta a riscos à saúde. Todavia, é necessária a condução de mais pesquisas para obtenção de evidências desta realidade e dos efeitos das alterações climáticas à saúde humana.

Referências

• AGÊNCIA NACIONAL DE ÁGUAS. Manual de procedimento para instalação, operação e manutenção de estações pluviométricas convencionais 2016. Disponível em: https://www.gov.br/ana/pt-br/assuntos/monitoramento-e-eventos-criticos/monitoramento-hidrologico/monitoramento-hidrologico-do-setor-eletrico/resolucao-conjunta-ana-aneel-127-2022/ManualInstalaoEstaesPluviomtricas.pdf Acesso em: 10 ago. 2023.
  » https://www.gov.br/ana/pt-br/assuntos/monitoramento-e-eventos-criticos/monitoramento-hidrologico/monitoramento-hidrologico-do-setor-eletrico/resolucao-conjunta-ana-aneel-127-2022/ManualInstalaoEstaesPluviomtricas.pdf
• ALESSI, N. A. Formam-se favelas e ganham importância no cenário urbano de São Paulo: Heliópolis e Paraisópolis. São Paulo, 2009. 169 f. Dissertação (Mestrado em Geografia humana) - Faculdade de Filosofia, Letras e Ciências Humanas, Universidade de São Paulo, São Paulo, 2009.
• APARICIO, P. et al The globe thermometer in comfort and environmental studies in buildings. Revista de La Construcción, v. 15, n. 3, p. 57-66, dez. 2016.
• ATLAS BRASIL. [Consulta]. Disponível em: http://www.atlasbrasil.org.br/consulta Acesso em: 19 out. 2023.
  » http://www.atlasbrasil.org.br/consulta
• AUBRECHT, C.; ÖZCEYLAN, D. Identification of heat risk patterns in the U.S. National Capital Region by integrating heat stress and related vulnerability. Environment International, v. 56, p. 65-77, jun. 2013.
• BAO, J.; LI, X.; YU, C. The construction and validation of the heat vulnerability index: a review. International Journal of Environmental Research and Public Health, v. 12, n. 7, p. 7220-7234, jun. 2015.
• BITENCOURT, D. P. et al Frequência, duração, abrangência espacial e intensidade das ondas de calor no Brasil. Revista Brasileira de Meteorologia, v. 31, n. 4, p. 506-517, dez. 2016.
• BORTOLI, K. C. R. de; VILLA, S. B. Adequação ambiental como atributo facilitador da resiliência no ambiente construído em Habitações de Interesse Social. Ambiente Construído, Porto Alegre, v. 20, n. 1, p. 391-422, jan./mar. 2020.
• CARVALHO, W. K. M. et al. Mudanças climáticas na metrópole paulista: uma análise de planos diretores e leis urbanísticas. Ambiente Construído, Porto Alegre, v. 20, n. 4, p. 143-156, out./dez. 2020.
• CASTILHO, J. V. A favelização do espaço urbano em São Paulo: estudo de caso - Heliópolis e Paraisópolis. São Paulo, 2013. 258 f. Dissertação (Mestrado em Habitat) - Faculdade de Arquitetura e Urbanismo, Universidade de São Paulo, São Paulo, 2013.
• CONLON, K. C. et al. Mapping human vulnerability to extreme heat: a critical assessment of heat vulnerability indices created using principal components analysis. Environmental Health Perspectives, v. 128, n. 9, p. 097001, set. 2020.
• CRAMER, M. N.; JAY, O. Biophysical aspects of human thermoregulation during heat stress. Autonomic Neuroscience, v. 196, p. 3-13, abr. 2016.
• CUTTER, S. L.; BORUFF, B. J.; SHIRLEY, W. L. Social vulnerability to environmental hazards. Social Science Quarterly, v. 84, n. 2, p. 242-261, maio 2003.
• DARELA FILHO, J. P. et al Socio-climatic hotspots in Brazil: how do changes driven by the new set of IPCC climatic projections affect their relevance for policy? Climatic Change, v. 136, n. 3, p. 413-425, mar. 2016.
• DI GIULIO, G. M. et al Mudanças climáticas, riscos e adaptação na megacidade de São Paulo, Brasil. Sustentabilidade em Debate, Brasília, v. 8, n. 2, p. 75-87, ago. 2017.
• DI NAPOLI, C.; PAPPENBERGER, F.; CLOKE, H. L. Verification of heat stress thresholds for a health-based heat-wave definition. Journal of Applied Meteorology and Climatology, v. 58, n. 6, p. 1177-1194, jun. 2019.
• DONG, J. et al Heatwave-induced human health risk assessment in megacities based on heat stress-social vulnerability-human exposure framework. Landscape and Urban Planning, v. 203, p. 103907-103915, nov. 2020.
• EAKIN, H.; LEMOS, M. C.; NELSON, D. Differentiating capacities as a means to sustainable climate change adaptation. Global Environmental Change, v. 27, p. 1-8, jul. 2014.
• GEIRINHAS, J. L. et al. Characterizing the atmospheric conditions during the 2010 heatwave in Rio de Janeiro marked by excessive mortality rates. Science of the Total Environment, v. 650, p. 796-808, fev. 2019.
• GEIRINHAS, J. L. et al Climatic and synoptic characterization of heat waves in Brazil. International Journal of Climatology, v. 38, n. 4, p. 1760-1776, set. 2017.
• GLOBAL HEAT HEALTH INFORMATION NETWORK (org.). Heat in the city: dialogue outcome brief. Dialogue Outcome Brief. 2020. Disponível em: https://ghhin.org/resources/heat-in-the-city-dialogue-outcome-brief/ Acesso em: 30 ago. 2023.
  » https://ghhin.org/resources/heat-in-the-city-dialogue-outcome-brief/
• GOBO, J. P. A. Bioclimatologia subtropical e modelização do conforto humano: da escala local à regional São Paulo, 2017. Tese (Doutorado em Geografia Física) - Faculdade de Filosofia, Letras e Ciências Humanas, Universidade de São Paulo, São Paulo, 2017.
• GOOGLE EARTH. Imagens street view de Paraisópolis, SP, Brasil. 2023. Disponível em: https://earth.google.com/web/ Acesso em: Acesso em: 22 out. 2023.
  » https://earth.google.com/web/
• GUO, Y. et al. Heat wave and mortality: a multicountry, multicommunity study. Environmental Health Perspectives, v. 125, n. 8, p. 087006, ago. 2017.
• HARDOY, D. J.; PANDIELLA, G. Urban poverty and vulnerability to climate change in Latin America. Environment and Urbanization v. 21, n. 1, p. 203-224, abr. 2009.
• HARLAN, S. L. et al. Neighborhood microclimates and vulnerability to heat stress. Social Science & Medicine, v. 63, n. 11, p. 2847-2863, dez. 2006.
• INSTITUTO BRASILEIRO DE GEOGRAFIA E ESTATÍSTICA. Panorama do Censo de 2022 2022. Disponível em: https://censo2022.ibge.gov.br/panorama/ Acesso em: 17 out. 2023.
  » https://censo2022.ibge.gov.br/panorama/
• INSTITUTO FAVELA DIZ. Pesquisa traça perfil de moradores de Heliópolis e Paraisópolis 2022. Disponível em: https://agenciabrasil.ebc.com.br/geral/noticia/2022-11/pesquisa-traca-perfil-de-moradores-de-heliopolis-e-paraisopolis Acesso em: 17 out. 2023.
  » https://agenciabrasil.ebc.com.br/geral/noticia/2022-11/pesquisa-traca-perfil-de-moradores-de-heliopolis-e-paraisopolis
• INSTITUTO NACIONAL DE METEOROLOGIA. Aviso de: onda de calor. 2023a. Disponível em: https://alertas2.inmet.gov.br/44758 Acesso em: 18 set. 2023.
  » https://alertas2.inmet.gov.br/44758
• INSTITUTO NACIONAL DE METEOROLOGIA. Balanço de São Paulo (SP) em setembro de 2023 2023b. Disponível em: https://portal.inmet.gov.br/uploads/notastecnicas/Balan%C3%A7o-SP-setembro2023_rvsd.pdf Acesso em: 24 out. 2023.
  » https://portal.inmet.gov.br/uploads/notastecnicas/Balan%C3%A7o-SP-setembro2023_rvsd.pdf
• INTERGOVERNMENTAL PANEL ON CLIMATE CHANGE. Climate change 2014: Impacts, Adaptation, and Vulnerability. Cambridge: Cambridge University Press, 2014. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change.
• INTERGOVERNMENTAL PANEL ON CLIMATE CHANGE. Climate Change 2021: The Physical Science Basis. Cambridge: Cambridge University Press, 2021. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change.
• INTERGOVERNMENTAL PANEL ON CLIMATE CHANGE. Climate Change 2022: Impacts, Adaptation, and Vulnerability. Cambridge: Cambridge University Press, 2022. Contribution of Working Group II to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change.
• INTERNATIONAL ORGANIZATION FOR STANDARDIZATION. 10551: ergonomics of the physical environment — subjective judgement scales for assessing physical environments. 2. ed. [S.l.]: ISO, 2019.
• JENKS, G. F. The data model concept in statistical mapping. International Yearbook of Cartography, v. 7, p. 186-190, 1967.
• KATZ, D. Configurações espaciais da interface entre os habitantes e a natureza da cidade: o caso da favela de Paraisópolis. São Paulo, 2008. 156 f. Dissertação (Mestrado em Arquitetura e Urbanismo) - Faculdade de Arquitetura e Urbanismo, Universidade São Judas Tadeu, São Paulo, 2008.
• LAPOLA, D. M. et al. Heat stress vulnerability and risk at the (super) local scale in six Brazilian capitals. Climatic Change, v. 154, n. 3-4, p. 477-492, maio 2019.
• LIBONATI, R. et al. Drought–heatwave nexus in Brazil and related impacts on health and fires: a comprehensive review. Annals of the New York Academy of Sciences, v. 1517, n. 1, p. 44-62, set. 2022.
• MALLEN, E.; STONE, B.; LANZA, K. A methodological assessment of extreme heat mortality modeling and heat vulnerability mapping in Dallas, Texas. Urban Climate, v. 30, p. 100528, dez. 2019.
• MARAGNO, D.; FONTANA, M. D.; MUSCO, F. Mapping heat stress vulnerability and risk assessment at the neighborhood scale to drive urban adaptation planning. Sustainability, v. 12, n. 3, p. 1056, fev. 2020.
• MAZIVIERO, M. C.; SILVA, A. S. O caso do Complexo Paraisópolis em gestões: diferenças conceituais em programas de intervenção em favelas em São Paulo. Revista Brasileira de Gestão Urbana, v. 10, n. 3, p. 500-520, jun. 2018.
• MCGREGOR, G. R. et al (ed.). Heatwaves and health: guidance on warning-system development. Geneva: World Meteorological Organization, 2015.
• MILLS, G. Cities as agentes of global change. International Journal of Climatology v. 27, n. 14, p. 1849-1857, ago. 2007.
• MONTEIRO, L. M.. Conforto térmico em espaços urbanos abertos: verificações modelares como aportes à exploração de abordagens. São Paulo, 2018. 517 f. Tese de livre docência - Faculdade de Arquitetura e Urbanismo, Universidade de São Paulo, São Paulo, 2018.
• NOBRE, C. A. et al Vulnerabilidades das megacidades brasileiras às mudanças climáticas: região metropolitana de São Paulo. São José dos Campos: INPE, 2011.
• NUTKIEWICZ, A. et al. Cool roofs can mitigate cooling energy demand for informal settlement dwellers. Renewable and Sustainable Energy Reviews, v. 159, p. 112183, maio 2022.
• OKE, T. R. et al. Urban climate Cambridge: Cambridge University Press, 2017.
• ORGANIZAÇÃO DAS NAÇÕES UNIDAS. Documentos temáticos da Habitat III: 20 - Habitação. ONU - HABITAT III. Maio 2015a. Disponível em: https://habitat3.org/wp-content/uploads/20-Habita%C3%A7%C3%A3o_final.pdf Acesso em: 19 nov. 2024.
  » https://habitat3.org/wp-content/uploads/20-Habita%C3%A7%C3%A3o_final.pdf
• ORGANIZAÇÃO DAS NAÇÕES UNIDAS. Documentos temáticos da Habitat III: 22 - Assentamentos informais. ONU - HABITAT III. Maio 2015b. Disponível em: https://habitat3.org/wp-content/uploads/22-Assentamentos-Informais_final.pdf Acesso em: 19 nov. 2024.
  » https://habitat3.org/wp-content/uploads/22-Assentamentos-Informais_final.pdf
• ÖZCEYLAN, D.; COSKUN, E. The relationship between turkey›s provinces› development levels and social and economic vulnerability to disasters. Journal of Homeland Security and Emergency Management, v. 9, n. 1, jan. 2012.
• PATZ, J. A. et al. Impact of regional climate change on human health. Nature, v. 438, n. 7066, p. 310-317, nov. 2005.
• PEREIRA, B.; SOUZA JUNIOR, T. P. de. Metabolismo celular e exercício físico: aspectos bioquímicos e nutricionais. São Paulo: Phorte Editora, 2010.
• PIZARRO, E. P. Interstícios e interfaces urbanos como oportunidades latentes: o caso da favela de Paraisópolis, São Paulo. São Paulo, 2014. 188 f. Dissertação (Mestrado em Tecnologia da Arquitetura) - Faculdade de Arquitetura e Urbanismo, Universidade de São Paulo, São Paulo, 2014.
• RATHI, S. K. et al A heat vulnerability index: spatial patterns of exposure, sensitivity and adaptive capacity for urbanites of four cities of india. International Journal of Environmental Research and Public Health, v. 19, n. 1, p. 283, dez. 2021.
• SÃO PAULO. Secretaria da Saúde registra aumento de 102% nos atendimentos por efeitos do calor Disponível em: https://www.saopaulo.sp.gov.br/ultimas-noticias/secretaria-da-saude-registra-aumento-de-102-nos-atendimentos-por-efeitos-do-calor/ Acesso em: 24 out. 2023.
  » https://www.saopaulo.sp.gov.br/ultimas-noticias/secretaria-da-saude-registra-aumento-de-102-nos-atendimentos-por-efeitos-do-calor/
• SHARMA, J.; RAVINDRANATH, N. H. Applying IPCC 2014 framework for hazard-specific vulnerability assessment under climate change. Environmental Research Communications, v. 1, n. 5, p. 051004, jun. 2019.
• SHINZATO, P.; DUARTE, D. H. S. Impacto da vegetação nos microclimas urbanos e no conforto térmico em espaços abertos em função das interações solo-vegetação-atmosfera. Ambiente Construído, Porto Alegre, v. 18, n. 2, p. 197-215, abr. 2018.
• SILVA, E. N. da; RIBEIRO, H. Alterações da temperatura em ambientes externos de favela e desconforto térmico. Revista de Saúde Pública, v. 40, n. 4, p. 663-670, abr. 2006.
• SINGH, R. et al Heatwave Guide for Cities The Hague: Red Cross Red Crescent Climate Centre, 2019.
• SOBRAL, H. R. Heat island in São Paulo, Brazil: effects on health. Critical Public Health, v. 15, n. 2, p. 147-156, jun. 2005.
• TARIFA, J. R.; ARMANI, G. As unidades climáticas urbanas da cidade de São Paulo: os climas e a (re)produção do espaço nas metrópoles. In: ENCUENTRO DE GEÓGRAFOS DE AMERICA LATINA, 8., Santiago, 2001. Anais [...] Santiago, 2001.
• UNITED NATIONS ENVIRONMENT PROGRAMME. Assessing human vulnerability to environmental change: concepts, issues, methods and case studies. Nairobi: UNEP, 2002.
• VINCENT, K. Creating an index of social vulnerability to climate change for Africa Norwich: Tyndall Centre for Climate Change Research, Working Paper, 2004. (Working Paper, n. 56).
• WILHELMI, O. V.; HAYDEN, M. H. Connecting people and place: a new framework for reducing urban vulnerability to extreme heat. Environmental Research Letters, v. 5, n. 1, p. 014021, mar. 2010.
• WOLF, T.; MCGREGOR, G. The development of a heat wave vulnerability index for London, United Kingdom. Weather and Climate Extremes, v. 1, n. 2, p. 59-68, set. 2013.
• WORLD METEOROLOGICAL ORGANIZATION; WORLD HEALTH ORGANIZATION. Health and climate: heat preparedness through early warning systems. heat preparedness through early warning systems. 2022. Disponível em: https://ghhin.org/resources/final-report-from-the-g7-health-communique-to-action-health-and-climate-heat-preparedness-through-early-warning-systems/ Acesso em: 30 ago. 2023.
  » https://ghhin.org/resources/final-report-from-the-g7-health-communique-to-action-health-and-climate-heat-preparedness-through-early-warning-systems/

Editado por

• Editor:
  Enedir Ghisi

Datas de Publicação

• Publicação nesta coleção
  16 Jun 2025
• Data do Fascículo
  2025

Histórico

• Recebido
  04 Out 2024
• Aceito
  21 Dez 2024