Modelagem dinâmica espacial de mudanças no uso do solo urbano: contribuição metodológica

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Produtividade de vacas mestiças leiteiras em sistema semi-intensivo nos municípios de Irituia e Mãe do rio do nordeste paraense

Verificar a influência de fatores do meio ambiente sobre características de produção e reprodução em grupos genéticos de bovinos leiteiros explorados em algumas fazendas nos municípios de Irituia e Mãe do rio no Nordeste Paraense. Foi utilizada uma média de 454 animais em lactação mestiças de Holandês, Pardo-Suíça e Girolando pertencentes a duas fazendas com sistema de criação semi-intensivo na época menos chuvosa. As médias e desvio-padrão para produção de leite total foram iguais a 1097,36 ± 330,47 Kg, com coeficiente de variação igual a 25,64% sendo que o grupo genético Holandesa, a época de parto menos chuvosa e o ano de 2007 apresentaram maior produção. O período de lactação apresentou efeito linear e crescente sobre a produção de leite total. O período de lactação apresentou média e desvio-padrão iguais a 218,17 ± 43,17 dias, a maior média para o período de lactação foi observada na época mais chuvosa e a média do período de lactação diminuiu no decorrer dos anos 2007 e 2008. A média e desvio padrão de intervalo entre partos encontrado no rebanho foi igual a 398,975 ± 60,85 dias. A época menos chuvosa apresentou uma média de intervalo entre partos menor que época mais chuvosa e a média de intervalo entre partos foi reduzindo a partir de 2006. A média de idade ao primeiro parto encontrado foi de 38,57 ± 5,81 meses. Os coeficientes de determinação foram maiores que 0,85 (R²a > 0,85), sendo que para o grupo genético Holandesa, os ajustes foram melhores em relação aos ajustes dos outros dois grupos genéticos. Para o grupo genético Girolando, o formato da curva de lactação diferiu dos outros dois grupos genéticos dificultando o ajuste pelas funções Gama Incompleta, Linear Hiperbólica e Polinomial Inversa. Os ajustes promovidos pela função polinomial inversa apresentaram ligeiro desvio em relação às outras duas funções. O grupo genético Holandesa apresentou produção ao pico e persistência um pouco mais elevados em relação aos outros grupos genéticos.

Entre o "Zaca” e o "Madre”: ofensas raciais, processos identitários e discursos de mestiçagem em duas escolas de Belém

Esta dissertação consiste em uma etnografia do cotidiano escolar, orientada conceitualmente a partir dos marcadores de diferenciação cor/raça, gênero e classe social. A motivação da pesquisa que iniciou a produção deste trabalho partiu da seguinte indagação: processos ofensivos verbais são processos identitários que se constroem a partir de referencias pedagógicos e de discursos docentes em cada escola? Com base em uma trajetória de pesquisa iniciada em 2006, realizada inicialmente na Escola Alexandre Zacharias de Assumpção e posteriormente estendida ao Colégio Madre Zarife Sales, no bairro do Guamá, periferia de Belém do Pará, tento discutir o alcance destas instituições no cotidiano educacional vivido fora delas por seus discentes, ao identificar práticas dinâmicas de viver os espaços urbanos da cidade a partir do ir e voltar da escola. As duas escolas são definidas como escolas de referência no bairro, pois simbolizam, para as clientelas de estudantes por elas atendidas, a possibilidade de sucesso profissional pós-escola entre muitas famílias no bairro. Por outro lado, tento mostrar as contradições internas nas duas escolas para sugerir a ocorrência de acessos diferenciados aos serviços educacionais e ao próprio capital cultural docente conforme o turno em que se estuda. Discuto destacadamente a ocorrência de trocas de ofensas verbais raciais e não-raciais entre os estudantes conforme os marcadores de diferenciação assinalados, para tentar responder a problematização da pesquisa. Depois tento identificar as várias modalidades simbólicas do discurso de mestiçagem como um discurso abrangente que é por mim concebido como a face inversa das ofensas verbais raciais e do não-reconhecimento do racismo nelas existentes pelo corpo institucional, isto é, por docentes, funcionários e gestores.

Influência de estruturas geológicas bidimensionais no campo geoeletromagnético na presença do eletrojato equatorial

A Terra atua como um grande magneto esférico, cujo campo assemelha-se àquele gerado por um dipolo magnético. Este campo apresenta mudanças de intensidade que variam com a localização e a hora local. A parte principal do campo geomagnético se origina no interior da Terra através de processos eletromagnéticos. Extensivos estudos mostraram ainda que existem contribuições de origem externa ao planeta, principalmente de origem solar. Dentre estas fontes há anomalias do campo magnético que surgem a partir de um aumento diurno da corrente elétrica em uma estreita faixa da ionosfera, de direção leste-oeste, centrada no equador magnético e denominada Eletrojato Equatorial (EEJ). Ocasionalmente estas correntes podem apresentar reversões de fluxo, sendo denominadas Contra-Eletrojato (CEJ). Vários autores têm estudado os efeitos do EEJ e CEJ sobre as observações geoeletromagnéticas. Eles estão interessados no efeito combinado do EEJ e estruturas geológicas condutivas 1-D e 2-D. Nestes trabalhos a estrutura 2-D sempre se apresentava paralela ao eletrojato, o que é uma hipótese bastante restritiva ao se modelar ambientes geológicos mais realistas, em que corpos bidimensionais podem ter qualquer strike em relação ao EEJ. Neste trabalho apresentamos a solução deste problema sem esta restrição. Assim, mostramos os campos geoeletromagnéticos devidos a estruturas bidimensionais que possuam strike oblíquo em relação ao EEJ, através de perfis dos campos elétrico e magnéticos calculados na superfície e formando direção arbitrária à heterogeneidade condutiva 2-D. Com esta resposta avaliamos ainda qual a influência que estruturas bidimensionais exercem sobre a resposta magnetotelúrica, sob influência do Eletrojato Equatorial. Durante o desenvolvimento deste trabalho, utilizamos o método de elementos finitos, tendo por fonte eletromagnética o EEJ e o CEJ, que por sua vez foram representados por uma combinação de distribuições gaussianas de densidade de corrente. Estas fontes foram decompostas nas direções paralela e perpendicular à estrutura 2-D, resultando nos modos de propagação TE₁ e TE₂ e TM acoplados, respectivamente. Resolvemos o modo acoplado aplicando uma Transformada de Fourier nas equações de Maxwell e uma Transformada Inversa de Fourier na solução encontrada. De acordo com os experimentos numéricos realizados em um modelo interpretativo da Anomalia Condutiva da Bacia do Parnaíba, formado por uma enorme estrutura de 3000 ohm-m dentro de um corpo externo condutivo (1 ohm-m), concluímos que a presença do CEJ causa uma inversão na anomalia, se compararmos com o resultado do EEJ. Concluímos também que para as frequências mais altas as componentes do campo elétrico apresentam menor influência da parte interna do corpo 2-D do que da parte externa. Já para frequências mais baixas este comportamento se observa com as componentes do campo magnético. Com relação à frequência, vimos os efeitos do “skin-depth”, principalmente nas respostas magnéticas. Além disso, quando a estrutura 2-D está paralela ao eletrojato, o campo elétrico é insensível à estrutura interna do modelo para todos os valores de frequência utilizados. Com respeito ao ângulo θ_h entre a heterogeneidade e a fonte, vimos que o modo TM se manifesta naturalmente quando θ_h é diferente de 0°. Neste caso, o modo TE é composto por uma parte devido à componente da fonte paralela à heterogeneidade e a outra devido à componente da fonte perpendicular, que é acoplada ao modo TM. Assim, os campos calculados têm relação direta com o valor de θ_h. Analisando a influência do ângulo entre a direção do perfil dos campos e o strike da heterogeneidade verificamos que, à medida que θ_h se aproxima de 90°, os campos primários tornam-se variáveis para valores de θ_p diferentes de 90°. Estas variações causam uma assimetria na anomalia e dão uma idéia da inclinação da direção do perfil em relação aos corpos. Finalmente, concluímos que uma das influências que a distância entre o centro do EEJ e o centro da estrutura 2-D, causa sobre as componentes dos campos está relacionado às correntes reversas do EEJ e CEJ, pois a 500 km do centro da fonte estas correntes têm máxima intensidade. No entanto, com o aumento da distância, as anomalias diminuem de intensidade. Nas sondagens MT, nós também usamos o EEJ e o CEJ como fonte primária e comparamos nossos resultados com a resposta da onda plana. Deste modo observamos que as componentes do campo geoeletromagnético, usadas para calcular a impedância, têm influência do fator de acoplamento entre os modos TE₂ e TM. Além disso, esta influência se torna maior em meios resistivos e nas frequências mais baixas. No entanto, o fator de acoplamento não afeta os dados magnetotelúricos em frequências maiores de 10^-2 Hz. Para frequências da ordem de 10^-4 Hz os dados MT apresentam duas fontes de perturbação: a primeira e mais evidente é devido à presença fonte 2-D (EEJ e CEJ), que viola a hipótese da onda plana no método MT; e a segunda é causada pelo acoplamento entre os modos TE₂ e TM, pois quando a estrutura bidimensional está obliqua à fonte 2-D temos correntes elétricas adicionais ao longo da heterogeneidade. Concluimos assim, que o strike de uma grande estrutura condutiva bidimensional relativamente à direção do EEJ ou CEJ tem de fato influência sobre o campo geomagnético. Por outro lado, para estudos magnetotelúricos rasos (frequências maiores de 10^-3 Hz) o efeito do ângulo entre a estrutura geológica 2-D e a direção do EEJ não é tão importante. Contudo, em estudos de litosfera frequências menores de 10^-3 Hz) o acoplamento entre os modos TE₂ e TM não pode ser ignorado.

Aplicação de deconvolução homomórfica a dados sísmicos

Um registro sísmico é frequentemente representado como a convolução de um pulso-fonte com a resposta do meio ao impulso, relacionada ao caminho da propagação. O processo de separação destes dois componentes da convolução é denominado deconvolução. Existe uma variedade de aproximações para o desenvolvimento de uma deconvolução. Uma das mais comuns é o uso da filtragem linear inversa, ou seja, o processamento do sinal composto, através de um filtro linear, cuja resposta de frequência é a recíproca da transformada de Fourier de um dos componentes do sinal. Obviamente, a fim de usarmos a filtragem inversa, tais componentes devem ser conhecidas ou estimadas. Neste trabalho, tratamos da aplicação a sinais sísmicos, de uma técnica de deconvolução não linear, proposta por Oppenheim (1965), a qual utiliza a teoria de uma classe de sistemas não lineares, que satisfazem um princípio generalizado de superposição, denominados de sistemas homomórficos. Tais sistemas são particularmente úteis na separação de sinais que estão combinados através da operação de convolução. O algoritmo da deconvolução homomórfica transforma o processo de convolução em uma superposição aditiva de seus componentes, com o resultado de que partes simples podem ser separadas mais facilmente. Esta classe de técnicas de filtragem representa uma generalização dos problemas de filtragem linear. O presente método oferece a considerável vantagem de que não é necessário fazer qualquer suposição prévia sobre a natureza do pulso sísmico fonte, ou da resposta do meio ao impulso, não requerendo assim, as considerações usuais de que o pulso seja de fase-mínima e que a distribuição dos impulsos seja aleatória, embora a qualidade dos resultados obtidos pela análise homomórfica seja muito sensível à razão sinal/ruído, como demonstrado.

Aplicação do método de Kalman a dados geofísicos

O filtro de Kalman é aplicado para filtragem inversa ou problema de deconvolução. Nesta dissertação aplicamos o método de Kalman, considerado como uma outra visão de processamento no domínio do tempo, para separar sinal-ruído em perfil sônico admitido como uma realização de um processo estocástico não estacionário. Em um trabalho futuro estudaremos o problema da deconvolução. A dedução do filtro de Kalman destaca a relação entre o filtro de Kalman e o de Wiener. Estas deduções são baseadas na representação do sistema por variáveis de estado e modelos de processos aleatórios, com a entrada do sistema linear acrescentado com ruído branco. Os resultados ilustrados indicam a aplicabilidade dessa técnica para uma variedade de problemas de processamento de dados geofísicos, por exemplo, ideal para well log. O filtro de Kalman oferece aos geofísicos de exploração informações adicionais para o processamento, problemas de modelamento e a sua solução.

Efeitos do strike das estruturas geológicas de duas dimensões nas pseudos-seções de IP resistividade

Os levantamentos de IP-resistividade efetuados na Serra dos Carajás não foram executados ortogonalmente às estruturas geológicas, pois utilizaram linhas anteriormente abertas pelas equipes de geoquímica. Este fato motivou este estudo teórico da influência da direção das linhas de medidas de IP-resistividade em relação ao "strike" da estrutura. Foi usado o programa de elementos finitos de Rijo (1977), desenvolvido para levantamentos perpendiculares às estruturas com as adaptações necessárias. A modificação principal foi na rotina de transformação inversa de Fourier. Para o caso simples dos levantamentos perpendiculares, a transformada inversa é uma integral discreta com apenas sete pontos. No entanto, para as medidas obliquas, o integrando é oscilatório, e portanto, a integral a ser calculada é mais complexa. Foi adaptado um método apresentado por Ting e Luke (1981), usando dezoito pontos em cada integração. Foi constatado que o efeito da direção da linha em relação ao "strike" é desprezível para ângulos maiores que 60 graus. Para ângulos menores, o efeito consiste no alongamento da anomalia, com pequenas alterações em seu centro. Não há uma maneira simples de compensar este efeito com mudanças nos parâmetros do modelo.

Estudo da hidrogeologia no município de Palestina do Pará usando o método de eletroresistividade

Esse trabalho de geofísica aplicada à hidrogeologia, foi realizado numa área localizada na porção sudeste do estado do Pará, mas precisamente na cidade de Palestina à margem esquerda do rio Araguaia. Ironicamente existem problemas de abastecimento de água potável para a população da cidade. Esse problema deve-se ao fato que, no momento, é antieconômico o tratamento e transporte de água do rio para a cidade. Considerando que a extração de água subterrânea é economicamente mais viável para a solução desse problema, o Departamento de Geofísica e o curso de Pós-graduação de Geofísica da Universidade Federal do Pará (UFPa), por solicitação da prefeitura daquele município e da Fundação Nacional de Saúde (FNS), participaram no levantamento geofísico de eletroresistividade para investigar os possíveis locais de ocorrência de água subterrânea. Foram realizadas, na área, 21 Sondagens Elétricas Verticais (SEVs) medidas na superfície, utilizando o arranjo Schlumberger com abertura máxima dos eletrodos AB e MN de 500m e 50m, respectivamente. Para a sua interpretação foram usados os programas de modelagem direta (EGSLIB/SEV1D) e inversa (EGSLIB/SEV1DIN) de RIJO (1994), de modelos de n camadas horizontais, homogêneas e isotrópicas. Com base nos resultados da interpretação das 21 SEVs, e também das informações da geologia local foram confeccionadas três mapas de contorno e seis secções geoelétricas. Dos seis horizontes que aparecem nessas secções, o mais favorável para a exploração hidrogeológica é o quarto, com espessura variando de 10 a 79 m, estando a uma profundidade média de 50m e com valores de resistividade variando de 150 — 850 Ωm, sendo correlacionavel com a Formação Potí. A analise final dos mapas e das secções geoelétricas, indicou que o local mais apropriado para a perfuração de poços, é uma área próximo à atual área de captação.

Interação genótipo x ambiente para a produção de leite na espécie bubalina utilizando inferência Bayesiana por meio de Amostradores de Gibbs

Com o objetivo de verificar a existência da interação genótipo x ambiente, sob a forma de heterogeneidade de variâncias para a produção de leite na espécie bubalina e o seu impacto na avaliação genética dos animais, utilizando a inferência Bayesiana por meio de Amostrador de Gibbs, foram utilizados 5.484 registros de produção de leite referentes à produções de 2.994 búfalas predominantemente Murrah, filhas de 150 reprodutores, acasalados com 1130 matrizes, cujos partos ocorreram entre os anos de 1974 e 2004. Os registros foram provenientes do Programa de Melhoramento Genético dos Bubalinos (PROMEBUL) com a adição de registros provenientes do rebanho da EMBRAPA Amazônia Oriental -EAO, localizada em Belém, Pará. Foram estabelecidas classes de rebanho-ano de parto e de acordo com o desvio padrão de cada classe, os registros de produção de leite foram classificados em classes de alto e baixo desvio-padrão fenotípico. Posteriormente, os dados foram analisados desconsiderando e considerando as classes de desvio-padrão. O modelo utilizado empregou os efeitos fixos referentes às classes de rebanho-ano, mês de parto e covariáveis idade da fêmea ao parto e duração da lactação, além do efeito aleatório de animal, ambiente permanente e ambiente temporário. Para os efeitos fixos, foi assumido distribuição à priori uniforme e para os componentes de (co)variâncias foram assumidas distribuições priori qui-quadrado inversa e Wishart invertida. As médias observadas e desvio-padrão para produção de leite nas classes de alto e baixo desvio-padrão e em análise geral, foram iguais a 1870,21±758,78, 1900,50±587,76 e 1885,48±677,98, respectivamente. As médias posteriores para os componentes de variâncias foram maiores na classe de alto desvio-padrão. A herdabilidade obtida na classe de alto desvio-padrão foi próxima do valor observado na análise geral e inferior ao valor encontrado na classe de baixo desvio-padrão fenotípico. A correlação genética para produção de leite entre as classes de desvio-padrão foi igual a 0,58. As correlações de Spearman entre os valores genéticos para a produção de leite obtidos em análise geral com os valores obtidos nas classes de alto e baixo desvio padrão foram iguais a 0,94 e 0,93, respectivamente, para todos os reprodutores. Para uma amostra dos 10 melhores reprodutores, as mesmas correlações foram iguais a 0,94 e 0,47, respectivamente. Tais resultados revelam presença de heterogeneidade de variâncias entre rebanhos e esta heterogeneidade de variâncias é resultante de fatores ambientais, que podem levar a uma classificação errônea dos melhores reprodutores geneticamente para a produção leite.

Formulação e implementação de um pacote para interpretar mapas magnéticos utilizando os polinômios de Walsh

O presente trabalho consiste na formulação de uma metodologia para interpretação automática de dados de campo magnético. Desta forma, a sua utilização tornará possível a determinação das fronteiras e magnetização de cada corpo. Na base desta metodologia foram utilizadas as características de variações abruptas de magnetização dos corpos. Estas variações laterais abruptas serão representadas por polinômios descontínuos conhecidos como polinômios de Walsh. Neste trabalho, muitos conceitos novos foram desenvolvidos na aplicação dos polinômios de Walsh para resolver problemas de inversão de dados aeromagnéticos. Dentre os novos aspectos considerados, podemos citar. (i) - O desenvolvimento de um algoritmo ótimo para gerar um jôgo dos polinômios "quase-ortogonais" baseados na distribuição de magnetização de Walsh. (ii) - O uso da metodologia damped least squares para estabilizar a solução inversa. (iii) - Uma investigação dos problemas da não-invariância, inerentes quando se usa os polinômios de Walsh. (iv) - Uma investigação da escolha da ordem dos polinômios, tomando-se em conta as limitações de resolução e o comportamento dos autovalores. Utilizando estas características dos corpos magnetizados é possível formular o problema direto, ou seja, a magnetização dos corpos obedece a distribuição de Walsh. É também possível formular o problema inverso, na qual a magnetização geradora do campo observado obedece a série de Walsh. Antes da utilização do método é necessária uma primeira estimativa da localização das fontes magnéticas. Foi escolhida uma metodologia desenvolvida por LOURES (1991), que tem como base a equação homogênea de Euler e cujas exigências necessárias à sua utilização é o conhecimento do campo magnético e suas derivadas. Para testar a metodologia com dados reais foi escolhida uma região localizada na bacia sedimentar do Alto Amazonas. Os dados foram obtidos a partir do levantamento aeromagnético realizado pela PETROBRÁS.

Avaliação do polimorfismo INDEL no gene TYMS em associação a resposta quanto ao uso de fluoropirimidinas em pacientes portadores de neoplasias do trato gastrointestinal

Introdução: O câncer consiste em um problema de saúde publica mundial, com estimativa de 27 milhões de casos novos e 17 milhões de mortes por câncer no ano de 2030. No Brasil, as estimativas para o câncer no ano de 2014, apontam a ocorrência de aproximadamente 580 mil casos novos. As fluoropirimidinas são usadas nos principais esquemas quimioterápicos direcionados a tumores do trato gastrointestinal. Nos últimos anos muito se tem investigado sobre causas de respostas individuais diferenciadas em relação ao tratamento quimioterápico; dessa forma têm-se buscado uma terapia individualizada que possa maximizar a eficácia dos medicamentos e minimizar os efeitos adversos associados aos fármacos. Objetivamos buscar a associação do INDEL (rs16430) do gene TYMS com o padrão de resposta ao tratamento oncológico de fármacos com base em fluoropirimidinas, de maneira a contribuir para o desenvolvimento da medicina personalizada. Material: Estudadas 151 amostras de sangue periférico de pacientes oncológicos tratados com fluoropirimidinas, da população da região amazônica brasileira com elevado grau de mistura interétinica. Foi genotipado um polimorfismo INDEL (rs16430) no gene TYMS envolvido na resposta ao tratamento com uso de fluoropirimidinas. Resultados: A investigação relatou que a maioria dos pacientes tinha doença avançada no momento do diagnóstico; 32,7% receberam tratamento com intenção paliativa; e 22,8% tratamento neoadjuvante. Nossos resultados evidenciam que o polimorfismo INDEL no gene TYMS demonstrou ter um efeito de proteção à progressão tumoral (p=0,033). Pacientes tratados com fluoropirimidinas que eram homozigotos selvagens (INS/INS) apresentaram uma proteção à progressão tumoral de 24% comparado com outros genótipos desse polimorfismo. As estimativas de ancestralidade genômica global da amostra investigada foram: 62,4% europeia; 25,2% nativo americana e 12,4% africana. Foi possível estabelecer uma correlação inversa entre o aumento da ancestralidade ameríndia e a presença de metástase (p=0,024). Conclusão: Estudos farmacogenéticos podem proporcionar uma terapia personalizada reduzindo mortalidade por toxicidades e aumentando a eficácia terapêutica, desta forma proporcionando um tratamento oncológico com melhores resultados clínicos.

Influência da temperatura nos parâmetros de modelos bi-paramétricos que predizem isotermas de adsorção de umidade do guaraná (Paullinia cupana) em pó

Oito modelos matemáticos bi-paramétricos, existentes na literatura e com larga aplicação na predição de isotermas de adsorção foram submetidos à análise. O guaraná (Paullinia cupana) em pó objeto deste estudo foi obtido em "spray dryer", a partir de um extrato hidroalcoólico. Ajustaram-se os pontos experimentais das isotermas de adsorção de umidade do produto à 15°C, 25°C e 35°C, por análise de regressão não-linear. Para estudar o efeito da temperatura nos parâmetros dos modelos utilizaram-se regressões dos tipos: linear, exponencial, logarítmica e inversa. Utilizou-se para fazer os ajustes o aplicativo STATGRAPHICS 5.1. Entre os modelos testados os que apresentam melhores resultados foram as equações de Handerson, Oswin e Mizrahi.

Matéria orgânica, capacidade de troca catiônica e acidez potencial no solo com dezoito cultivares de cana-de-açúcar

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Métodos de avaliação do estado nutricional para o algodoeiro no centro-oeste do Brasil

Pós-graduação em Agronomia (Produção Vegetal) - FCAV

Geologia, mineralogia e geoquímica dos fosfatos de Sapucaia (Bonito - PA)

O depósito mineral de Sapucaia, situado no município de Bonito, região nordeste do Estado do Pará, é parte de um conjunto de ocorrências de fosfatos de alumínio lateríticos localizados predominantemente ao longo da zona costeira dos estados do Pará e Maranhão. Estes depósitos foram alvos de estudo desde o início do século passado, quando as primeiras descrições de “bauxitas fosforosas” foram mencionadas na região NW do Maranhão. Nas últimas décadas, com o crescimento acentuado da demanda por produtos fertilizantes pelo mercado agrícola mundial, diversos projetos de exploração mineral foram iniciados ou tiveram seus recursos ampliados no território brasileiro, dentre estes destaca-se a viabilização econômica de depósitos de fosfatos aluminosos, como o de Sapucaia, que vem a ser o primeiro projeto econômico mineral de produção e comercialização de termofosfatos do Brasil. Este trabalho teve como principal objetivo caracterizar a geologia, a constituição mineralógica e a geoquímica do perfil laterítico alumino-fosfático do morro Sapucaia. A macrorregião abrange terrenos dominados em sua maioria por rochas pré-cambrianas a paleozóicas, localmente definidas pela Formação Pirabas, Formação Barreiras, Latossolos e sedimentos recentes. A morfologia do depósito é caracterizada por um discreto morrote alongado que apresenta suaves e contínuos declives em suas bordas, e que tornam raras as exposições naturais dos horizontes do perfil laterítico. Desta forma, a metodologia aplicada para a caracterização do depósito tomou como base o programa de pesquisa geológica executada pela Fosfatar Mineração, até então detentora dos respectivos direitos minerais, onde foram disponibilizadas duas trincheiras e amostras de 8 testemunhos de sondagem. A amostragem limitou-se à extensão litológica do perfil laterítico, com a seleção de 44 amostras em intervalos médios de 1m, e que foram submetidas a uma rota de preparação e análise em laboratório. Em consonância com as demais ocorrências da região do Gurupi, os fosfatos de Sapucaia constituem um horizonte individualizado, de geometria predominantemente tabular, denominado simplesmente de horizonte de fosfatos de alumínio ou crosta aluminofosfática, que varia texturalmente de maciça a cavernosa, porosa a microporosa, que para o topo grada para uma crosta ferroalumino fosfática, tipo pele-de-onça, compacta a cavernosa, composta por nódulos de hematita e/ou goethita cimentados por fosfatos de alumínio, com características similares aos do horizonte de fosfatos subjacente. A crosta aluminofosfática, para a base do perfil, grada para um espesso horizonte argiloso caulinítico com níveis arenosos, que repousa sobre sedimentos heterolíticos intemperizados de granulação fina, aspecto argiloso, por vezes sericítico, intercalados por horizontes arenosos, e que não possuem correlação aparente com as demais rochas aflorantes da geologia na região. Aproximadamente 40% da superfície do morro é encoberta por colúvio composto por fragmentos mineralizados da crosta e por sedimentos arenosos da Formação Barreiras. Na crosta, os fosfatos de alumínio estão representados predominantemente pelo subgrupo da crandallita: i) série crandallita-goyazita (média de 57,3%); ii) woodhouseíta-svanbergita (média de 15,8%); e pela iii) wardita-millisita (média de 5,1%). Associados aos fosfatos encontram-se hematita, goethita, quartzo, caulinita, muscovita e anatásio, com volumes que variam segundo o horizonte laterítico correspondente. Como os minerais pesados em nível acessório a raro estão zircão, estaurolita, turmalina, anatásio, andalusita e silimanita. O horizonte de fosfatos, bem como a crosta ferroalumínio-fosfática, mostra-se claramente rica em P₂O₅, além de Fe₂O₃, CaO, Na₂O, SrO, SO₃, Th, Ta e em terras-raras leves como La e Ce em relação ao horizonte saprolítico. Os teores de SiO₂ são consideravelmente elevados, porém muito inferiores aqueles identificados no horizonte argiloso sotoposto. No perfil como um todo, observa-se uma correlação inversa entre SiO₂ e Al₂O₃; entre Al₂O₃ e Fe₂O₃, e positiva entre SiO₂ e Fe₂O₃, que ratificam a natureza laterítica do perfil. Diferente do que é esperado para lateritos bauxíticos, os teores de P₂O₅, CaO, Na₂O, SrO e SO₃ são fortemente elevados, concentrações consideradas típicas de depósitos de fosfatos de alumínio ricos em crandallita-goyazita e woodhouseítasvanbergita. A sucessão dos horizontes, sua composição mineralógica, e os padrões geoquímicos permitem correlacionar o presente depósito com os demais fosfatos de alumínio da região, mais especificamente Jandiá (Pará) e Trauíra (Maranhão), bem como outros situados além do território brasileiro, indicando portanto, que os fosfatos de alumínio de Sapucaia são produtos da gênese de um perfil laterítico maturo e completo, cuja rocha fonte pode estar relacionada a rochas mineralizadas em fósforo, tais como as observadas na Formação Pimenteiras, parcialmente aflorante na borda da Bacia do Parnaíba. Possivelmente, o atual corpo de minério integrou a paleocosta do mar de Pirabas, uma vez que furos de sondagem às proximidades do corpo deixaram claro a relação de contato lateral entre estas unidades.