Search

Orientador: Prof. Ph.D. Carlos Roberto Sanquetta Coorientadores: Prof. Dra. Ana Paula Dalla Corte, Prof. Dr. Alexandre Behling Dissertação (mestrado) - Universidade Federal do Paraná, Setor de Ciências Agrárias, Programa de Pós-Graduação em Engenharia Florestal. Defesa : Curitiba, 02/03/2018 Inclui referências: p.127-140 Área de concentração: Manejo Florestal Resumo: O objetivo deste trabalho foi avaliar o potencial produtivo e definir Unidades de Manejo em plantios comerciais de acacia-negra, por meio da estruturacao de banco de dados do meio fisico, capaz de integrar as informacoes numericas e espaciais, mediante o uso de tecnicas estatisticas multivariadas e geoestatisticas, fornecendo subsidios para um melhor conhecimento potencial da especie e para o uso sustentado das terras. Para isso, o presente estudo foi conduzido em plantio comercial de acacia-negra situado no municipio de Cristal, localizado na regiao da Encosta do Sudeste, no estado do Rio Grande do Sul, e constituido por tres povoamentos seminais implantados nos anos de 2003, 2004 e 2005, perfazendo o total de 939,1 hectares de efetivo plantio. Os dados dendrometricos foram coletados em 91 unidades amostrais, nas quais foram tambem coletadas as variaveis do meio fisico. As espacializacoes do Indice de Sitio (IS_7) e das medidas de variabilidade foram utilizadas como diferentes estrategias para definicao das Unidades de Manejo. As medidas de variabilidade utilizadas foram os desvios padroes das variaveis dap (desvDAP), ht (desvHT), g (desvg), V (desvVOL) e o indice de uniformidade (PV50) do plantio. A dependencia espacial para a variavel Indice de Sitio, aos sete anos de idade, foi considerada fraca (8,6 %), portanto, nao justificando a modelagem geoestatistica para essa variavel. Entre as medidas de variabilidade, o melhor resultado foi encontrado utilizando o PV50, cuja dependencia espacial foi considerada forte (78,72 %). Os resultados revelam que o modelo geoestatistico exponencial permitiu visualizar o comportamento do PV50, possibilitando definir tres classes de uniformidade: I) 14,119 "¥ 23,119; II) 10,113 "¥ 14,119 e III) 1,112 "¥ 10,113. O estudo tambem mostrou que ha uma tendencia em locais onde a altitude e o teor de materia organica (MO) sao maiores apresentarem maiores valores para o PV50. Sendo assim, conclui-se que as Unidades de Manejo para acacia-negra podem ser aplicadas no estudo da relacao entre o indice de uniformidade (baseado no PV50) e atributos do meio fisico (altitude e materia organica do solo), possibilitando definir a espacializacao de tres Unidades de Manejo diferenciadas para acacia-negra, de acordo com as classes de PV50. Palavras-chave: Indice de Sitio, medidas de variabilidade, classes de uniformidade Abstract: Therefore, the objective of this work was to study the productive potential and to define management units in Black Wattle commercial plantations, by means of the database structure of the physical environment, capable of integrating the numerical and spatial information, through the use of multivariate and geostatistical techniques, providing subsidies for a better potential knowledge of this species and for sustained land use. For this, the present study was carried out in a commercial plantation of Black Wattle located in the region of Encosta do Sudeste, state of Rio Grande do Sul, and consisting of three seminal stands implanted in 2003, 2004 and 2005, making a total of 939,1 hectares of effective planting. The dendrometric data were collected in 91 sample units, in which the variables of the physical environment were also collected. Spatialisations of the Site Index (IS_7) and variability measures used as different strategies for the definition of Management Units. The variability measures used were the standard desviations of the variables dap (desvDAP), ht (desvHT), g (desvg), V (desvVOL) and the uniformity index (PV50) of the stand. The spatial dependence for the Site Index, at seven years, was considered weak (8,6 %), therefore, not justifying the geostatistical modeling for this variable. Among the measures pf variability, the best search result was found using the PV50, whose spatial dependence was considered strong (78,72 %). The results show that the exponential geostatistical model allowed to visualize the behavior of the PV50, allowing to define three classes of uniformity: I) 14,119 "¥ 23,119; II) 10,113 "¥ 14,119 e III) 1,112 "¥ 10,113. The study also showed that there is a tendency in places where altitude and organic matter content (MO) are higher present higher values for the PV50. Thus, it can be concluded that the management units for Black Wattle can be applied in the study of the relationship between the uniformity index (based on the PV50) and attributes of the physical environment (altitude and soil organic matter), making it possible to define spatialization o three different Management Units for Black Wattle, according to the classes of PV50. Keywords: Site Index, measures of variability, uniformity classes

Orientador : Prof. Dr. Carlos Roberto Sanquetta Coorientador : Drª. Patrícia Póvoa de Mattos e Dr. Evaldo Muñoz Braz Tese (doutorado) - Universidade Federal do Paraná, Setor de Ciências Agrárias, Programa de Pós-Graduação em Engenharia Florestal. Defesa: Curitiba, 2//03/2016 Inclui referências ao final de cada capítulo: f: 20-22;53-59;86-89;110-113;138-143 Área de concentração: Manejo Florestal Resumo: A situação de fragilidade a que está sujeita a Floresta Ombrófila Mista (FOM) e a araucária demandam o aprofundamento do conhecimento sobre seus aspectos demográficos e de crescimento para subsidiar estratégias de manejo sustentável de tais recursos, que pode consistir em ferramenta imprescindível à conservação. Esta tese teve como objetivo geral determinar os padrões de recrutamento, mortalidade e crescimento em remanescentes de FOM no extremo sul paranaense, com destaque para as populações de araucária, e apresentar estratégia para o manejo sustentável e otimizado da espécie. O alvo deste estudo consiste em remanescentes de FOM localizados no município de General Carneiro, extremo sul do estado do Paraná. No local, foram instaladas oito unidades permanentes de 1 ha cada em 1999, em que todos os indivíduos arbóreos com diâmetro à altura do peito (dap) ? 10 cm são mensurados anualmente, contabilizando-se também mortalidade e recrutamento. O trabalho está dividido em quatro capítulos, que tratam respectivamente dos seguintes temas: dinâmica demográfica, com enfoque sobre a mortalidade e recrutamento da comunidade arbórea e araucária; dinâmica do crescimento e projeções da estrutura diamétrica da floresta por meio de modelos otimizados; análise de padrões de crescimento em araucárias centenárias, com base na dendrocronologia; e simulações de manejo da araucária, visando à obtenção de estrutura diamétrica otimizada que maximize a produção ao fim de um ciclo de corte por meio de desbastes iniciais. De modo geral, as taxas de recrutamento, seja da comunidade arbórea ou da araucária, foram superiores às de mortalidade. Observou-se tendência de efeito direto da densidade da floresta sobre as taxas de mortalidade e efeito inverso sobre o recrutamento. A precipitação apresentou relação inversa ao recrutamento, indicando maior ocorrência desse processo em períodos de menor umidade. A mortalidade foi inversamente relacionada à Razão de Umidade, indicando maior sensibilidade dos indivíduos a períodos secos. Não foi observada relação entre variáveis meteorológicas e os processos dinâmicos da araucária. Em relação ao incremento diamétrico da comunidade e da araucária observou-se que esse tende a ser maior em locais menos estocados, com padrão assimétrico, em que árvores de maior porte crescem a uma taxa maior. Os modelos gerados foram capazes de representar apropriadamente os processos dinâmicos de incremento e sobrevivência da comunidade e da araucária, podendo ser utilizados satisfatoriamente para a projeção da estrutura da floresta. A análise dos padrões de crescimento em araucárias, com base em 16 discos da base de indivíduos da espécie, revelou que o incremento apresenta relação inversa com o tamanho dos indivíduos até as classes intermediárias, estabilizando-se nas maiores classes e aumentando na última classe, possivelmente em razão dos indivíduos avaliados não terem atingido ainda o valor assintótico. A araucária apresenta autocorrelação temporal no crescimento, ou seja, indivíduos superiores tendem a apresentar maior crescimento desde as classes iniciais e esse padrão de crescimento tende a ser mantido ao longo do tempo. Em relação à otimização das estruturas diamétricas de araucária, as simulações indicaram que a aplicação de desbastes iniciais otimizados foi o método que proporcionou o maior incremento volumétrico ao longo dos dois ciclos de corte, maximizando o Valor Presente Líquido (VPL) da produção volumétrica. Estratégias de manejo florestal sustentável são necessárias para favorecer a conservação dos remanescentes florestais e da araucária. Palavras-chave: araucária, dinâmica florestal, incremento, modelagem, otimização. Abstract: The situation of fragility which Araucaria Forest (FOM) and araucaria populations have been subject requires a deeper understanding of their demographics and growth patterns to support sustainable management strategies of these resources, which may consist of an essential conservation tool. This thesis aimed to determine the patterns of recruitment, mortality and growth in FOM remnants, especially in populations of araucaria, in the Southern Paraná, and present a strategy for sustainable and optimized management of the species. The study was performed in FOM remnants in the municipality of General Carneiro, Paraná state, where eight permanent sample plots of 1 ha each were installed in 1999, where all trees with diameter at breast height (dbh) ? 10 cm are measured annually, also counting mortality and recruitment. This work is divided into four chapters, dealing respectively with the following topics: demographic dynamics, focusing on mortality and recruitment of tree community and araucaria; growth dynamics and projections of diameter distribution of forest through optimized individual models; growth patterns analysis of centenary araucarias, based on dendrochronology; and management simulations for araucaria in order to obtain an optimal diameter structure that maximizes the production at the end of a cutting cycle through initial thinning. Overall, recruitment rates, whether of the tree community or araucaria, were higher than mortality. There was direct effect trend of forest density on mortality rates and reverse effect on recruitment. Precipitation showed an inverse relation to recruitment, indicating higher rates of this process in lower humidity periods. Mortality was inversely related to humidity ratio, indicating higher sensitivity of individuals to dry periods. No relationship was found between meteorological variables and the dynamic processes of araucaria. Regarding the diametric increment of the community and the araucaria, was noted that this tends to be higher in less stocked places, with asymmetric pattern, in which larger trees grow at a higher rate. Individual models generated were able to represent properly the dynamic processes of increment and survival of the community and the araucaria, which can be used satisfactorily for the projection of the forest structure. Analysis of growth patterns in centenary araucarias revealed that the increment has an inverse relation to the size of individuals up to the intermediate classes, stabilizing in larger classes and increasing in the last class, possibly because the individuals have not reached yet their asymptotic value. Araucaria presents temporal autocorrelation in growth, i.e., superior individuals tend to have higher growth from the initial classes and this growth pattern tends to be maintained over time. Regarding the optimization of the diameter structure of araucaria, simulations showed that the application of optimized initial thinning was the method that provided the greatest volume increase over the two cutting cycles, maximizing net present value (NPV) of the volumetric production. Sustainable forest management strategies are required for the conservancy of forest remnants and araucaria populations. Keywords: araucaria, forest dynamics, increment, modeling, optimization.

Orientadora : Profª. Drª. Ana Paula Dalla Corte Coorientador : Prof. Dr. Carlos Roberto Sanquetta Dissertação (mestrado) - Universidade Federal do Paraná, Setor de Ciências Agrárias, Programa de Pós-Graduação em Engenharia Florestal. Defesa: Curitiba, 19/02/2016 Inclui referências : f. 102-122 Área de concentração : Manejo florestal Resumo: O Índice de Área Foliar (IAF) é considerado como parte da estrutura mais sensível da floresta, pois as folhas do dossel regulam alguns processos fundamentais da produtividade florestal. A modelagem do IAF é considerada crítica e pouco acurada, principalmente no que se refere a florestas nativas. Diante disso, o presente estudo teve como objetivo testar diferentes métodos para estimativa do IAF por meio de detecção remota, em um fragmento de Floresta Ombrófila Mista, em São João do Triunfo, no Estado do Paraná. A área de estudo possui 4 parcelas permanentes, de modo que, foram distribuídas 81 unidades amostrais por parcela e 33 abrangendo a área total. A estimativa do IAF pelo método terrestre foi gerada automaticamente pelo sensor óptico CI-110, considerada como verdade de campo. O IAF estimado pelo método orbital foi gerado a partir de imagens do sensor Pléiades, baseado na relação empírica entre os índices de vegetação. Para o IAF obtido pelo método orbital foram testados os modelos propostos por SEBAL (2002) e Duchemin et al. (2006), e foram analisados os índices de vegetação NDVI, SR, SAVI e EVI. O IAF obtido a campo apresentou uma grande variabilidade, com valores entre 6,01 e 8,01, fator que está relacionado ao grande número de espécies existentes no local. Os índices de vegetação NDVI e SR não foram capazes de apresentar a variabilidade das espécies florestais para a estimativa do IAF, e os índices SAVI e EVI foram os que geraram melhores resultados, possibilitando a estimativa do IAF. O EVI apresentou correlação linear com o IAF obtido a campo variando de 0,54 a 0,77, e o SAVI entre 0,55 e 0,78. O valor médio do IAF gerado pelo modelo de Duchemin et al. (2006) foi bastante homogêneo, onde as parcelas variaram de 4,11 a 4,34, porém, não obteve correlação linear com os dados obtidos a campo. O modelo proposto por SEBAL (2002) apresenta uma maior variação entre as parcelas, exibindo valores entre 2,65 a 3,14, com r variando de 0,60 a 0,72, tendo forte correlação com os dados de IAF obtidos a campo. Logo, foram gerados modelos de regressão para estimativa do IAF baseados nos índices de vegetação SAVI e EVI, de modo que, todas as equações ajustadas foram consideradas satisfatórias, com R² maior do que 0,6 e Syx(%) de 4%. Os modelos gerados apresentaram resultados satisfatórios, o que demonstra a aplicabilidade desta metodologia. Palavras-chave: Floresta com Araucária; Fotografias Hemisféricas; Sensoriamento Remoto. Abstract: The Leaf Area Index (LAI) is considered as part of the most sensitive structure of the forest, because the leaves of the canopy regulate some fundamental processes of the forest productivity. The modeling of LAI is considered critical and little accurate, mainly with regard to native forests. Therefore, the present study aimed to test different methods for estimate of LAI by means of remote sending, in a Subtropical Ombrophilous Forest fragment, in São João do Triunfo, in Paraná state. The area of study has 4 permanent plots, so, 81 sampling units were distributed by plots and 33 covering the total area. The estimate of LAI by terrestrial method was generated automatically by the optical sensor CI-110, considered as truth field. The LAI estimated by orbital method was generated from images of Pléiades sensor, based on the empirical relation among the vegetation indexes. For the LAI obtained by orbital method it was tested the models proposed by SEBAL (2002) and Duchemin et al. (2006), and it was analyzed the indexes of vegetation NDVI, SR, SAVI and EVI. The LAI obtained by field presented a great variability, with values between 6.01 and 8.01, factor that is towed to the great number of existing species in the local. The indexes of vegetation NDVI and SR were not capable of presenting the variability of the forest species for estimate of LAI, and the indexes SAVI and EVI were the ones that generated better results, allowing the estimate of LAI. The EVI presented linear correlation with the LAI obtained by field varying from 0.54 to 0.77, and the SAVI between 0.55 and 0.78. The average value of LAI generated by the model of Duchemin et al. (2006) was fairly homogeneous, where the plots vary from 4.11 to 4.34; however, it was not obtained linear correlation with the data obtained by field. The model proposed by SEBAL (2002) presents a greater variation between the plots, displaying values between 2.65 and 3.14, with r varying from 0.60 to 0.72, and it has strong correlation with the LAI data obtained by field. Thus, it was generated models of regression for estimate of LAI based on the indexes of vegetation SAVI and EVI, in such a way that all the adjusted equations were considered satisfying, with R² greater than 0.6 and Syx of 4%. The generated models presented satisfying results, which demonstrates the applicability of this methodology. Keywords: Araucaria forest; Hemispherical Photographs; Remote Sensing.

Orientadora : Profª. Drª. Ana Paula Dalla Corte Coorientador : Prof. Dr. Carlos Roberto Sanquetta Dissertação (mestrado) - Universidade Federal do Paraná, Setor de Ciências Agrárias, Programa de Pós-Graduação em Engenharia Florestal. Defesa: Curitiba, 11/03/2016 Inclui referências : f. 59-71 Área de concentração : Tecnologia e utilização de produtos florestais Resumo: A acurácia e precisão na avaliação quantitativa de povoamentos para fins comerciais são fundamentais, especialmente a estimativa do volume individual. Nesse contexto, as equações de volume e funções de afilamento são importantes ferramentas para a estimativa do volume individual, ambos com métodos tradicionalmente aplicados e sólida teoria. Por outro lado, a inovação matemática é muito dinâmica, com diversos outros recursos disponíveis para o estudo de aproximação de função (regressão), métodos esses que ainda carecem de estudos de base, podendo potencialmente melhorar as estimativas no campo de modelagem florestal. Assim, o trabalho tem como objetivo principal analisar a acurácia de técnicas de aprendizado de máquina em relação a um modelo volumétrico e a uma função de afilamento, para a espécie Acacia mearnsii De Wild. Foram utilizados dados de cubagem Acacia mearnsii, variando de 1 a 10 anos. A base de dados foi dividida em 60% para ajuste e o restante para validação. Foram ajustadas equações com o modelo volumétrico de Schumacher e Hall e com a função de afilamento de Hradetzky, comparados com três algoritmos de aprendizado de máquina: k-vizinho mais próximo (k-NN), Random Forest (RF) e Redes Neurais Artificiais (RNA) para a estimativa do volume total e diâmetro referente à altura relativa. Os modelos foram ranqueados conforme estatísticas de erros, bem como observadas as distribuições destes. Para a estimativa do volume em função do dap e altura, a RNA e o modelo de Schumacher e Hall apresentaram melhores resultados no ranqueamento do que o k-NN e RF. Os métodos de aprendizado de máquina aplicados se mostraram mais acurados que o polinômio de Hradetzky para estimativas da forma da árvore, tais como o diâmetro ao longo do fuste e volume total. Os modelos de AM se mostraram adequados como alternativa na modelagem tradicionalmente aplicada na mensuração florestal, contudo a sua utilização deve ser cuidadosa devida a maior possibilidade de supertreinamento a base de ajuste. Palavras-chave: Modelagem. Inteligência Artificial. Mineração de dados. Abstract: Accuracy and precision are essential topics when it comes to the quantitative evaluation of a forest stand with commercial purposes, especially in regards to the estimation of individual volume. In one hand, due to its solid theory, volume equations and taper functions are important tools for estimating individual volume. On the other hand, the mathematical breakthrough is dynamic, having several resources for the study of approximation functions (regression). However, these new methods still lack baseline studies and may potentially improve estimates in the forest modeling field. Under this circumstances, this study aims to analyze the accuracy of machine learning techniques in regards to a volumetric model and a taper function for the species Acacia mearnsii De Wild. Acacia mearnsii scaling data, from 1 to10 years, were used in order to achieve the objective of this work. The database was divided into 60% for adjustment and the remainder for validation. Equations were adjusted with the volume model of Schumacher and Hall and the taper function Hradetzky and were compared to three machine learning algorithms: nearest neighbor models (k-NN), Random Forest (RF) and Artificial Neural Networks (ANN) to estimate of the total volume and diameter concerning relative height. The models were ranked according to its statistics errors and distributions. ANN and Schumacher and Hall model showed better results in regards to the estimation of volume as a function of dbh and height. The applied machine learning methods were more accurate than Hradetzky polynomial to estimates of tree shape, such as the diameter along the stem and total volume. The machine learning models were satisfactory as an alternative to the traditional methods in forest measurement. However, due to its higher possibility of overtraining the adjustment basis, it should be carefully used. Keywords: Modelling. Artificial Intelligence. Data mining.t

Orientador : Prof. Dr. Sylvio Péllico Netto Coorientador : Prof. Ph.D. Carlos Roberto Sanquetta Coorientador : Profª. Dª. Ana Paula Dalla Corte Tese (doutorado) - Universidade Federal do Paraná, Setor de Ciências Agrárias, Programa de Pós-Graduação em Engenharia Florestal. Defesa: Curitiba, 12/09/2016 Inclui referências : f. 92-114 Área de concentração : Manejo florestal Resumo: Neste trabalho é discutido o tema aditividade dos componentes da biomassa de árvores. Foram apresentados procedimentos para estabelecer sistemas de equações de biomassa aditivos, propriedade indispensável e que determina ser as estimativas de biomassa dos componentes equivalente à biomassa total. Para avaliar e discutir tal contexto várias indagações foram feitas, que estão apresentadas no texto, e que foram devidamente respondidas por meio de informações da literatura e pelas avaliações experimentais discutidas e avaliadas no conteúdo deste trabalho. Um experimento em florestas de Acacia mearnsii De Wild. (acácia negra) foi conduzido nos locais que concentram os plantios da espécie no estado do Rio Grande do Sul. Para tais condições foram ajustadas equações de biomassa total e para os componentes copa e fuste por meio de weighted-nonlinear seemingly unrelated regressions (WNSUR, estimador de sistemas de equações) e por meio de generalized nonlinear least squares (GNLS, estimador independente). Na análise do desempenho dos estimadores, especialmente no que se refere às estimativas de biomassa e das estatísticas de ajuste das equações, pode-se concluir a respeito da indiferença do emprego de equações ajustadas por meio dos estimadores independentes e de sistemas de equações. Por outro lado, essa conclusão diverge quando avaliada a consistência biológica e a eficiência dos estimadores. Estimar as equações de biomassa para os componentes e para a biomassa total independentemente não é realista, pois do ponto de vista biológico as estimativas de biomassa são inconsistentes, ou seja, não são aditivas. As estimativas de biomassa dos componentes e da biomassa total resultantes de equações ajustadas por meio de sistemas de equações viabilizaram intervalos de confiança menores em relação às equações ajustadas independentemente, sendo mais eficientes. O estimador WNSUR apresenta melhores propriedades biológicas e estatísticas para estimar equações alométricas de biomassa dos componentes e biomassa total quando comparado com a estimação independente, sendo o método que deve ser utilizado. Por fim, os exemplos apresentados na revisão da literatura são úteis para a compreensão dos princípios envolvidos e também para propósitos de difusão das técnicas. Palavras-chave: Consistência biológica. Estimador eficiente. Derivação de pesos. Intervalo de confiança. Abstract:This study discusses the additivity of the components of tree biomass. Procedures have been presented to establish additive biomass equation systems, essential property to determine the estimated biomass of components equivalent to total biomass. To evaluate and discuss this context, several inquiries were made, which are presented in the text, which were duly answered by the literature information and by the experimental evaluations discussed in the content of this work. An experiment in forests of Acacia mearnsii De Wild. (black wattle) was conducted on sites that concentrate the plantations of the species in the state of Rio Grande do Sul. For such conditions, equations for total biomass, canopy and stem components were fitted by weighted-nonlinear seemingly unrelated regressions (WNSUR, estimator of equation systems) and by generalized nonlinear least squares (GNLS, independent estimator). In the performance of the estimators' analysis, especially with regard to biomass estimates and fitted statistics of the equations, it can be concluded about the equivalence of using equations adjusted by means of independent estimators or by system of equations. Moreover, this conclusion diverges when evaluating the biological consistency and efficiency of estimators. Estimating the biomass equations of the components and the total biomass independently is not realistic, because biologically these biomass estimates are inconsistent, i.e. they are not additive. The biomass estimates resultant from components and from total biomass, when adjusting system of equations, turns possible getting smaller confidence intervals than using equations adjusted independently, what make them more efficient. The estimator WNSUR provides better biological properties and statistics to estimate biomass allometric equation of components and total biomass when compared to independent estimators, thus this should be the method to be used. Finally, the examples presented in the literature review are useful for understanding the principals involved and dissemination purposes of these techniques. Keywords: Biological consistency. Efficient estimator. Derivation of weights. Confidence interval.

Orientador: Celso Carnieri Co-orientador: Carlos Roberto Sanquetta Co-orientador: Julio Eduardo Arce Tese (doutorado) - Universidade Federal do Paraná, Setor de Ciencias Agrárias, Programa de Pós-Graduaçao em Engenharia Florestal. Defesa: Curitiba, 2004 Inclui bibliografia Área de concentraçao: Manejo florestal Este estudo foi realizado nas áreas de reflorestamento da Indústria Pedro N. Pizzatto, localizada na região de General Carneiro, Centro-Sul do Paraná, com o objetivo geral de modelar o planejamento florestal otimizado de longo prazo visando maximizar o lucro da empresa florestal em questão, levando em consideração o estoque de carbono existente nas florestas com o decorrer dos tratamentos silviculturais. O estudo apresentou os seguintes objetivos específicos: estimar o volume de sortimento retirado nos desbastes e no corte raso, em 11 projetos, utilizando 4 simulações de regimes de manejo além do regime de manejo atual da empresa; estimar a quantidade de carbono (C), existente acima do solo, que é retirada nos desbastes simulados para cada regime .de manejo e para os projetos, assim como a quantidade remanescente dessa variável; apresentar a quantidade total de carbono que seria fixada aos 20 anos para os projetos e regimes de manejo estudados, verificando se existe diferença significativa entre eles; avaliar economicamente os projetos submetidos aos diferentes regimes de manejo; identificar a melhor opção de manejo florestal para cada projeto da empresa, maximizando o lucro total e mantendo um nível de carbono na área. De posse dos dados advindos de 11 Projetos da referida empresa, utilizou-se o Software SISPINUS para prognosticar o volume, a altura e o diâmetro desses projetos, submetidos a 5 diferentes regimes de manejo, para um horizonte de planejamento de 21 anos. Através dos resultados obtidos de diâmetro e altura estimou-se a quantidade de carbono remanescente e retirado nos anos do horizonte de planejamento, assim como a . quantidade total de carbono aos 20 anos de cada projeto com o objetivo de verificar se existe diferença significativa entre as médias encontradas para os projetos e para os regimes de manejo. Na programação linear adotou-se 36 cenários, sendo que estes variaram em função da demanda de madeira (exigida pela empresa, 20% menor e 30% menor), da taxa de desconto anual (8, 10, 12 e 16,5% a.a.) e da estratégia comercial (Compra e Venda maximizando o VPL, DownGrade minimizando a compra e DownGrade maximizando o VPL). Adotou-se também um modelo de programação em metas, onde foram estabelecidas metas de carbono (baseada no resultado obtido com o melhor cenário), metas econômicas (baseada no faturamento anual da empresa) e metas de demanda para cada sortimento. Após a análise dos resultados encontrados as seguintes conclusões podem ser tiradas: o maior volume dos sortimentos estudados (serraria, tomo pequeno e tomo grande) concentrou-se nos primeiros anos do horizonte de planejamento, em todos os regimes de manejo adotados, devido à idade avançada dos projetos no ano considerado como base para início do estudo; a quantidade remanescente total de carbono, considerando todos os projetos, apresentou-se instável durante o horizonte de planejamento, sendo a maior quantidade encontrada nos primeiros e nos últimos anos do horizonte de planejamento, nos 5 regimes de manejo , adotados no estudo; comparando-se estatisticamente os regimes de manejo estudados verificou-se que estes apresentam diferenças estatísticas significantes ao nível de 5% de probabilidade, sendo que o regime da empresa e o regime 4 apresentaram as menores médias e os regimes I, 2 e 3 as maiores médias para quantidade de carbono total aos 20 anos de idade; com relação às médias de carbono total encontradas aos 20 anos para os projetos estudados pode-se dizer que apenas o Projeto 12 apresentou resultado significativo, sendo a maior média de todos os projetos, o que pode ser justificado pelo seu alto índice de sítio; o cenário que apresentou o melhor resultado para a função objetivo no modelo de programação linear foi o de número 5. Este cenário considera uma taxa de desconto de 8% ao ano no cálculo do Valor Presente Líquido, a compra e venda da madeira e uma demanda 30% do que a utilizada pela empresa; os projetos apresentaram áreas muito quebradas, sendo que algumas muito pequenas. Isto pode ser causado pelo número grande de restrições impostas pelo modelo e a quebra dos projetos grandes pode estar sugerindo a criação de novas unidades de manejo. This study was carried out in lands of the Industrias Pedro N. Pizzatto, located in General Carneiro county, Southern-Central Parana State, Brazil. The general objective of the study was develop a long-term optimization forest planning model which maximizes the revenue of the forest company, taking into consideration the carbon stocked in the forests under different management regimes. The specific objectives of the study were: to estimate the timber volume harvested by the occasion of thinnings and final cut in 11 stands, using 4 simulated management regimes besides the standard regime adopted by the company; to estimate the aboveground carbon stock (C) removed by the simulated thinnings for each management regime and stand as well as the remaining stock of this variable; to estimate the total amount of carbon captured by the forest stands in a 20-year period for each management regime and comparing the statistic of them; to evaluate the economic result of the different management regimes; to identify the best management option for each stand in the study area which maximizes the total revenue keeping a certain level of carbon in the forest. The data come from 11 stands of the company, and the SISPINUS growth and yield simulator was used to predict volume, mean tree height and dbh of these stands managed under five different options for a planning horizon of 21 years. Removed and remaining carbon stocks of each stand were estimated through dbh and height for each management event and for the total period of time aiming at identifying significant differences among the management regimes. Linear programming techniques were applied to 36 scenarios defined as a function the log demand (required by the company, 20% less and 30% less), the annual interest rate (8, 10, 12, and 16.5% per annum), as well as the commercial strategy (purchase and sale maximizing NPV -Net Present Value, DownGrade minimizing purchase and DownGrade maximizing NPV). A goal programming approach was also used to indicate some carbon stock goals for the best previously selected scenario, economical goals (based upon the annual income ofthe company) and demand goals of each log use (assortment). The results indicated the greatest volume of the assortments (sawnlog, small veneer log, long veneer log) were concentrated in the first years of the planning horizon for all management regimes due to the age of the stands used for the beginning of the analysis (the current existing stands are only mature ones). It was also noticed that the total carbon stock remaining in the stands under study was instable during the planning horizon and concentrated in the last years of the planning horizon for the five regimes under analysis. The statistical comparison among the management regimes showed significant differences at 0.05 probability level. The current management regime used by the company and regime 4 gave the lowest means, whereas regimes I, 2, and 3 gave the greatest means of total carbon stock in 20 years. The best condition was reached by stand 12, the greatest mean for the most favorable site index. The best scenario indicated by the objective function in the linear programming model was number 5. This option takes into consideration an 8% interest rate for calculating NPV, purchase and sale oflogs in the market and a demand of 30% of the woody raw material consumed by the factory. The stands originated after the application of the model showed irregular areas, some very small. This might be caused by the great number of restrictions imposed by the model and the conversion of large stands into small ones might be suggesting the necessity of creating new and variable management units.