Modelagem Bayesiana do risco de infecção tuberculosa para estudos com perdas de seguimento

OBJETIVO: Desenvolver um modelo estatístico baseado em métodos Bayesianos para estimar o risco de infecção tuberculosa em estudos com perdas de seguimento, comparando-o com um modelo clássico determinístico. MÉTODOS: O modelo estocástico proposto é baseado em um algoritmo de amostradores de Gibbs, utilizando as informações de perdas de seguimento ao final de um estudo longitudinal. Para simular o número desconhecido de indivíduos reatores ao final do estudo e perdas de seguimento, mas não reatores no tempo inicial, uma variável latente foi introduzida no novo modelo. Apresenta-se um exercício de aplicação de ambos os modelos para comparação das estimativas geradas. RESULTADOS: As estimativas pontuais fornecidas por ambos os modelos são próximas, mas o modelo Bayesiano apresentou a vantagem de trazer os intervalos de credibilidade como medidas da variabilidade amostral dos parâmetros estimados. CONCLUSÕES: O modelo Bayesiano pode ser útil em estudos longitudinais com baixa adesão ao seguimento.

Modelagem e simulação da dinâmica do arrastamento por ar

Os compostos orgânicos voláteis constituem uma fonte vulgar de contaminação da água subterrânea, a qual pode ser eliminada pela tecnologia do arrastamento por ar (air stripping) em colunas com enchimento desordenado e utilizando fluxos das fases em contra-corrente. Propõe-se neste trabalho uma nova metodologia de dimensionamento destas colunas, para qualquer tipo de enchimento e de contaminante, onde não há necessidade de se arbitrar nenhum diâmetro, onde se evita o recurso a ábacos experimentais e onde o regime hidráulico conveniente é seleccionado à partida. O procedimento proposto foi algoritmizado e convertido num programa em linguagem C++. Para verificar e testar não só o dimensionamento mas também o comportamento teórico estacionário e dinâmico construiu-se de raiz uma coluna experimental. Seleccionou-se como contaminante uma solução de clorofórmio em água destilada. A experimentação permite, ainda, corrigir o coeficiente de transferência de massa global teórico estimado pelas correlações de Onda e que depende de inúmeros parâmetros nem sempre controláveis experimentalmente. Apresenta-se, em seguida, um modelo original de simulação dinâmica do comportamento da coluna e que é constituído por um sistema de equações diferenciais não lineares (parâmetros distribuidos). No entanto, se os débitos forem arbitrados como constantes, o sistema passa a ser linear apesar de não possuir solução analítica evidente (p.e. por transformações integrais). A discretização por diferenças finitas permitiu superar estas dificuldades. Existe uma notável concordância entre os valores experimentais e os previstos no modelo.

Sete motivações teóricas para o uso da modelagem baseada no indivíduo em ecologia

A modelagem baseada no indivíduo tem sido crescentemente empregada para analisar processos ecológicos, desenvolver e avaliar teorias, bem como para fins de manejo da vida silvestre e conservação. Os modelos baseados no indivíduo (MBI) são bastante flexíveis, permitem o uso detalhado de parâmetros com maior significado biológico, sendo portanto mais realistas do que modelos populacionais clássicos, mais presos dentro de um rígido formalismo matemático. O presente artigo apresenta e discute sete razões para a adoção dos MBI em estudos de simulação na Ecologia: (1) a inerente complexidade de sistemas ecológicos, impassíveis de uma análise matemática formal; (2) processos populacionais são fenômenos emergentes, resultando das interações entre seus elementos constituintes (indivíduos) e destes com o meio; (3) poder de predição; (4) a adoção definitiva, por parte da Ecologia, de uma visão evolutiva; (5) indivíduos são entidades discretas; (6) interações são localizadas no espaço e (7) indivíduos diferem entre si.

Modelagem Sísmica via métodos das diferenças finitas: caso da bacia do Amazonas

Este trabalho tem por objetivo apresentar os resultados da modelagem sísmica em meios com fortes descontinuidades de propriedades físicas, com ênfase na existência de difrações e múltiplas reflexões, tendo a Bacia do Amazonas como referência à modelagem. As condições de estabilidade e de fronteiras utilizadas no cálculo do campo de ondas sísmicas foram analisadas numericamente pelo método das diferenças finitas, visando melhor compreensão e controle da interpretação de dados sísmicos. A geologia da Bacia do Amazonas é constituída por rochas sedimentares depositadas desde o Ordoviciano até o Recente que atingem espessuras da ordem de 5 km. Os corpos de diabásio, presentes entre os sedimentos paleozóicos, estão dispostos na forma de soleiras, alcançam espessuras de centenas de metros e perfazem um volume total de aproximadamente 90000 Km³. A ocorrência de tais estruturas é responsável pela existência de reflexões múltiplas durante a propagação da onda sísmica o que impossibilita melhor interpretação dos horizontes refletores que se encontram abaixo destas soleiras. Para representar situações geológicas desse tipo foram usados um modelo (sintético) acústico de velocidades e um código computacional elaborado via método das diferenças finitas com aproximação de quarta ordem no espaço e no tempo da equação da onda. A aplicação dos métodos de diferenças finitas para o estudo de propagação de ondas sísmicas melhorou a compreensão sobre a propagação em meios onde existem heterogeneidades significativas, tendo como resultado boa resolução na interpretação dos eventos de reflexão sísmica em áreas de interesse. Como resultado dos experimentos numéricos realizados em meio de geologia complexa, foi observada a influência significativa das reflexões múltiplas devido à camada de alta velocidade, isto provocou maior perda de energia e dificultou a interpretação dos alvos. Por esta razão recomenda-se a integração de dados de superfície com os de poço, com o objetivo de obter melhor imagem dos alvos abaixo das soleiras de diabásio.

Uma proposta de protocolo para a obtenção de variáveis visando estudos de modelagem ecológica em sistemas pesqueiros fluviais da Amazônia

Pesquisas sobre as pescarias amazônicas apresentam problemas nas coletas de dados que impedem ou restringem as respostas resultantes da modelagem. Este trabalho tem como objetivo propor um protocolo de pesquisa, discriminando as principais variáveis necessárias à construção de modelos para analisar a pesca na região.

Observações hidrológicas e resultados de modelagem no espalhamento sazonal e espacial da pluma de água Amazônica

Neste estudo são apresentadas as análises espaciais e temporais dos dados sinóticos coletados durante quatro campanhas oceanográficas do programa “REVIZEE/SCORE-N” realizadas em março-abril de 1995, outubro-novembro de 1997, abril-maio de 1999 e agosto-setembro de 2001. Em complemento às análises dos dados hidrológicos, foram utilizados resultados do modelo numérico CLIPPER - Experimento ATL6 (1/6 deg.). A ênfase deste artigo é avaliar o deslocamento sazonal e espacial da pluma de água Amazônica nas regiões localizadas à sudoeste, noroeste e em frente à foz do rio Amazonas. Dados hidrológicos e resultados do modelo demonstram que durante os meses de março-maio a pluma do rio Amazonas move-se para sudoeste, ao longo da costa norte brasileiro, tendendo a diminuir o deslocamento da pluma para sudoeste durante os meses de agosto-setembro. De outubro a novembro, não foi observado ou registrado o deslocamento da pluma para sudoeste. Resultados do modelo revelam anéis gerados pela retroflexão da Corrente Norte do Brasil (CNB) em torno de 7º N, e também a influência de água proveniente do rio Amazonas em volta desses anéis durante o mês de agosto.

Complexação de íons de metais por matéria orgânica dissolvida: modelagem e aplicação em sistemas reais

Os íons de metais pesados são conhecidos tanto pela sua importância fisiológica e industrial, bem como pelo risco ambiental e à saúde humana. Para elucidar o comportamento dessas espécies nos corpos hídricos, os quais recebem grande parte da descarga de metais, seja de origem antrópica ou por fontes naturais, é necessário entender as interações que elas apresentam com o meio, principalmente a especiação química. Um dos mais importantes processos pelos quais passam as espécies metálicas em corpos aquáticos naturais é a interação com a matéria orgânica dissolvida (MOD), que pode ser por adsorção, reações de troca iônica ou por complexação. Neste trabalho foram realizados vários experimentos com o objetivo de descrever o comportamento da complexação de três importantes cátions Cu(II), Cd(II) e Pb(II) com a matéria orgânica (ácido húmico comercial), sob condições diversas de força iônica em meio tamponado. Os resultados foram avaliados de acordo com o modelo de van den Berg e Kramer para a complexação de metais. O modelo foi aplicado na determinação da capacidade de complexação dos íons em amostras reais, oriundas de três rios maranhenses que integram a Amazônia legal: Itapecuru, Bacanga e Pericumã. Nas águas dos rios utilizou-se o parâmetro carbono orgânico dissolvido (COD) para expressar a MOD. Os resultados confirmaram forte interação entre a MOD e íons de metais pesados e que o modelo de van den Berg e Kramer é satisfatório para se estimar a constante de complexação (K) e a concentração de sítios de complexação (Lt). Nas amostras simuladas em laboratório a ordem de complexação dos metais foi Cu(II) > Cd(II) > Pb(II) e a capacidade de complexação mostrou ser linear em função da concentração de ácido húmico comercial. Acredita-se que por ter menor raio iônico, o íon Cu(II) possui maior afinidade com os sítios de complexação. Nas amostras retiradas dos corpos aquáticos, observou-se que o rio com maior concentração de COD (rio Bacanga) apresentou maior capacidade de complexação; entretanto, a ordem foi Pb(II) > Cu(II) > Cd(II) provavelmente devido à presença de íons Cu(II) em maior quantidade nas águas dos rios.

Distribuição geográfica potencial de espécies americanas do caranguejo "violinista" (Uca spp.) (Crustacea, Decapoda) com base em modelagem de nicho ecológico

Caranguejos do gênero Uca Leach, 1814 (caranguejo "violinista") são um grupo bem conhecido e caracterizado por um acentuado dimorfismo sexual e assimetria da quela do macho. Atualmente, estão descritas 97 espécies no mundo. Objetivou-se estimar a distribuição geográfica potencial de 4 espécies do gênero Uca que ocorrem na costa do continente Americano: Uca maracoani Latreille, 1802-1803, U. uruguayensis Nobili, 1901, U. panacea Novak & Salmon, 1974 e U. monilifera Rathbun, 1914. Para modelar a distribuição dessas espécies nas Américas foram utilizados pontos de ocorrência compilados da literatura. Para a modelagem foram utilizados os programas Maxent e GARP a partir de 10 variáveis climáticas e três variáveis topográficas. Todas as variáveis foram convertidas para uma malha com resolução de 0,0417 graus. Nos dois modelos (Maxent e GARP) as espécies apresentaram distribuição geográfica maior do que sugerido por outros trabalhos de registro de ocorrência, com exceção de U. monilifera. Segundo o critério de área sob a curva (AUC), os modelos gerados pelo GARP apresentaram melhores resultados do que os modelos do Maxent. Entretanto, avaliando em conjunto os resultados dos dois modelos é possível melhor estabelecer planos de conservação para espécies com habitat restrito (U. panaceae e U. monilifera), além de recomendar um aumento na amostragem de U. maracoani no nordeste brasileiro e U. uruguyaensis no sudeste brasileiro, a fim de detectar possíveis aumentos na sua distribuição geográfica com base nas predições dos modelos de nicho.

Modelagem de ocorrência e densidade de juvenis de Callinectes sapidus (Decapoda, Portunidae) em dois estuários do Rio Grande do Sul, com dados zero-inflacionados

Dados de contagem de juvenis de siri-azul (Callinectes sapidus Rathbun, 1896) coletados em dois estuários do Rio Grande do Sul são objeto do presente estudo. Por se encontrarem zero-inflacionados, esses dados motivaram a formulação de modelos hierárquicos, que quantificam o efeito das covariáveis categóricas mês e local sobre a probabilidade de ocorrência e densidade dessas populações, levando em conta a detecção imperfeita. Foram também desenvolvidos modelos não-hierárquicos para comparação. Uma abordagem Bayesiana foi adotada para a estimação dos parâmetros dos modelos por simulação Monte Carlo com Cadeias de Markov (MCMC). A comparação entre modelos foi feita com o Critério de Informação da Deviância (DIC). Os modelos hierárquicos apresentaram ajustes melhores que os modelos convencionais, mitigaram o problema do excesso de zeros e permitiram analisar simultaneamente as probabilidades de ocorrência e a densidade de juvenis de siri-azul. No estuário da Lagoa dos Patos, a probabilidade de ocorrência de juvenis na Classe 2 aumenta com a distância da desembocadura, enquanto em Tramandaí os pontos intermediários apresentam as maiores probabilidades. Em ambos os estuários a ocorrência é mais provável nos meses de verão e de inverno. A densidade de juvenis da Classe 2 apresenta marcada variação em relação aos meses do ano sendo, em geral, maior no estuário de Tramandaí.

Parametrização do transporte dispersivo de solutos em solos estruturados: heterogeneidade do meio, percurso de transporte e modelagem

O transporte dispersivo de solutos no solo é parametrizado por meio de equações difusivas e de modelos que partem da existência de uma relação de causalidade entre as características físicas do meio poroso, da sua condição de umidade e da manifestação macroscópica do processo nele ocorrente. A validade dessa proposição, como critério na elaboração de modelos dos processos que governam o transporte de solutos no solo, é discutida heuristicamente. Com base em índices obtidos a partir de uma análise de momentos de curvas de avanço (BTCs), é possível avaliar objetivamente o vínculo que as características do sistema estabelecem com o atendimento do pressuposto do equilíbrio local e, assim, com o sucesso ou não da adoção de uma abordagem de modelagem. Assim, BTCs assimétricas nada mais são do que a manifestação macroscópica de um processo incompleto, e o alcance de uma distribuição simétrica da concentração de solutos só pode ser garantido com a adoção de um "comprimento mínimo de mistura".

Distribución a mesoescala de comunidades bentónicas y clorofila-a en sedimentos superficiales de áreas potenciales para reconstrucción paleoceanográfica. Otoño 2003.

Se analizó la distribución espacial de la clorofila-a y otras propiedades de los sedimentos superficiales, así como del macrobentos (>0,5 mm) durante la Primera Etapa del Crucero de Evaluación de Recursos Demersales y Paleoceanografía BIC Olaya 0305-06, del 20 al 25 de mayo 2003. La materia orgánica “fresca” de origen fitoplanctónico se distribuyó de forma heterogénea en ambas áreas, alcanzando concentraciones muy elevadas sobre la plataforma exterior y talud superior, hasta el momento, las más altas determinadas para el margen continental peruano. Se encontró una asociación positiva entre contenido de materia orgánica “fresca” superficial y presencia de sedimentos laminados o bandeados en ambas áreas. En general, el área frente a Callao presentó un mayor enriquecimiento orgánico, menores biomasas de macrofauna y Thioploca spp. y mayor abundancia de nemátodes que el área de Pisco. Se discuten posibles mecanismos que expliquen la retención, distribución y acumulación de materia orgánica en los sedimentos superficiales de ambas áreas. La distribución batimétrica de la macrofauna sugiere un debilitamiento de las condiciones hipóxicas frente a Pisco, a partir de los 400 m.

Simulação de cenários alternativos de uso da terra em uma microbacia utilizando técnicas de modelagem e geoprocessamento

O modelo SWAT (Soil and Water Assessment Tool) foi aplicado na simulação de cenários alternativos de uso da terra na microbacia hidrográfica do Ribeirão dos Marins, em Piracicaba (SP), no biênio 1999/2000. Dois cenários foram simulados. No primeiro, o uso atual foi mantido numa faixa de mata ciliar de 30 m em toda a extensão dos cursos d'água e de 50 m ao redor das nascentes, de acordo com o Código Florestal. No segundo cenário, como as pastagens ocupavam 30,9 % da área da microbacia, nas encostas mais íngremes, com alto potencial erosivo, as áreas de pastagem foram substituídas por vegetação florestal. As simulações dos dois cenários foram comparadas com as condições do cenário atual em termos de produção de sedimentos. Os cenários geraram diferentes padrões espaciais da produção de sedimentos. Uma redução de 94,0 % na produção de sedimentos foi obtida com a substituição da pastagem por vegetação nativa (cenário 2). No cenário 1, a redução foi de 10,8 %. Esses resultados evidenciam a necessidade de tratar a paisagem em bacias hidrográficas de forma global, identificando as "áreas sensíveis ambientalmente", onde são necessárias práticas de controle dos processos erosivos e não somente práticas de proteção dispensada aos cursos d'água por meio da mata ciliar.

Modelagem numérica para o processo da evaporação da água do solo

Foram realizados vários ensaios laboratoriais para avaliar o desempenho de um modelo numérico em simular o processo unidimensional da evaporação da água do solo. Este modelo foi desenvolvido a partir da linearização da equação de Richards no espaço e do uso da técnica iterativa de Newton-Raphson, para resolver essa equação não-linear no tempo, e leva em conta o fluxo de vapor dentro do perfil de solo e a taxa de evaporação na sua superfície. Os resultados mostraram que o modelo foi capaz de simular satisfatoriamente o processo, tanto em meios porosos homogêneos quanto estratificados

Modelagem da capacidade de suporte de carga e quantificação dos efeitos das operações mecanizadas em um Latossolo Amarelo cultivado com cafeeiros

O Estado de Minas Gerais é o maior produtor de café do Brasil, com mais de 50 % da produção total, sendo grande parte da área cultivada originalmente ocupada por vegetação de cerrado. Essa região é destaque na cafeicultura em razão da alta tecnologia empregada na mecanização, que submete o solo ao tráfego de máquinas, tornando preocupante a possibilidade de disseminação da compactação do solo. Objetivou-se com este estudo: (a) propor modelos de capacidade de suporte de carga de um Latossolo Amarelo distrófico típico (LAd) cultivado com cafeeiro em função da pressão de preconsolidação e da umidade, (b) determinar, com o uso destes modelos, em três profundidades e em duas épocas de amostragem, a influência das operações de manejo na estrutura do Latossolo Amarelo nas estações: seca e chuvosa. As amostragens foram realizadas em três camadas (0-3, 10-13 e 25-28 cm), sendo 30 amostras coletadas no local onde não houve tráfego e 10 amostras na linha de tráfego para cada equipamento, coletadas nos meses de agosto de 2002 e março de 2003 para quantificar seu efeito na estrutura do solo nas duas épocas: seca e chuvosa, respectivamente. As amostras indeformadas foram utilizadas nos ensaios de compressão uniaxial. Foram determinados também granulometria, matéria orgânica, densidade de partículas, limite de contração e limite de plasticidade do solo. O modelo de capacidade de suporte de carga do LAd é expresso pela equação sigmap = 10(2,72-1,17U), em que sigmap significa a pressão de preconsolidação e U a umidade. Os efeitos das operações mecanizadas na lavoura cafeeira foram quantificados com o modelo de capacidade de suporte de carga e com o seu uso foi possível quantificar a influência das operações de manejo na estrutura do Latossolo Amarelo distrófico típico.

Modelagem de atributos físico-hídricos do solo numa bacia hidrográfica da região do Alto Rio Grande, MG

O presente trabalho objetivou desenvolver modelos para estimativa dos atributos físico-hídricos do solo porosidade drenável (PD) e capacidade total de retenção de água (CTA), com base em atributos físicos de fácil e rápida obtenção (textura, argila dispersa em água, volume total de poros, densidade do solo e densidade de partícula) e a determinação de matéria orgânica. Os atributos foram coletados na bacia hidrográfica do Ribeirão Marcela, na camada de 0 a 15 cm, obedecendo aos grids de 240 x 240 e 60 x 60 m, totalizando 165 pontos amostrados. Trabalhou-se com regressão múltipla linear, constituindo-se variáveis por meio da combinação dos diferentes atributos entre si e estimando os respectivos coeficientes pelo método dos mínimos quadrados, utilizando-se o programa SAS for Windows. Os modelos gerados para CTA e PD apresentaram boa qualidade estatística, com elevados coeficientes de determinação, baixos erros de estimativa e variáveis significativas, o que permitiu sua aplicação na estimativa desses atributos na bacia hidrográfica do Ribeirão Marcela, representativa do domínio dos Latossolos na região do Alto Rio Grande, fornecendo assim subsídios para implantação de projetos agrícolas e ambientais, especialmente no auxílio à parametrização de modelos de simulação hidrossedimentológica, os quais estão em desenvolvimento para aplicação na região.