Modelagem e simulação do campo de velocidades do vento em tormentas elétricas

Esta tese propõe um modelo para descrever o campo de velocidade do vento originado numa tormenta elétrica. O escoamento descendente que ocorre durante uma tormenta elétrica madura é modelado como um campo de fluxo axial simétrico. A evolução da velocidade do vento com o tempo é considerada por meio de uma função de modulação empírica. As flutuações da velocidade ao redor da velocidade média em evolução são representadas por um processo aleatório do tipo ruído branco em faixa. Todos os parâmetros utilizados no modelo estão relacionados a variáveis meteorológicas, tais como a altura e diâmetro da nuvem cumulonimbus causadora do evento, a diferença de pressão e a distância entre o solo e a base da nuvem. Adicionalmente considera-se que a nuvem é transportada numa direção pré escolhida pelo vento de fundo, que pode ser devido à circulação geral na região ou a um ciclone extra-tropical. O campo de velocidades do vento resulta da superposição de ambos campos. O modelo é validado comparando as velocidades do vento determinadas em eventos simulados, com registros de estações meteorológicas. Adicionalmente, perfis verticais previstos pelo modelo são comparados com vários perfis para tormentas elétricas propostos na literatura. Finalmente, séries de velocidade máximas anuais a altura padrão de 10m foram simuladas para três estações meteorológicas no sul do Brasil, para as quais dados experimentais são conhecidos, em termos do número médio de tormentas elétricas por ano e das funções distribuição de probabilidade que caracterizam o vento de fundo. A comparação com séries de registros máximos existentes é satisfatória, dando credibilidade adicional ao modelo proposto.

Fatores intervenientes na capacidade de atendimento de praças de pedágio

Este trabalho tem como objetivo o levantamento e análise de fatores intervenientes na capacidade de processamento de veículos em cabines de praças de pedágio com o recolhimento manual de tarifas. Buscando o entendimento de como estes fatores interferem nos tempos de atendimento nas cabines foi realizada uma análise estatística e posterior modelagem, que utilizou redes neurais artificiais. Redes neurais artificiais são úteis no entendimento de problemas com alto grau de complexidade, que agregam diversas variáveis de entrada com relações não-lineares entre si. As variáveis de entrada escolhidas para a modelagem foram forma de pagamento, intensidade de fluxo, valor das tarifas e classes de veículos. A variável de saída foi o tempo de atendimento nas cabines de cobrança de pedágios. Foram obtidos três modelos que buscaram refletir a variação dos tempos de atendimento para um mesmo conjunto de dados de entrada: Modelo de Tempos Mínimos de Atendimento; Modelo de 85° Percentil de Tempos de Atendimento, e; Modelo de Tempos Máximos de Atendimento. As análises de sensibilidade dos modelos indicaram que tempos de atendimento são fortemente influenciados pelo fluxo de veículos nas praças. Quanto mais intenso o fluxo de veículos, tempos mínimos de atendimento tendem a sofrer leve aumento, indicando pequena perda de rendimento do processo. Perda de rendimento pode ser resultado de (i) necessidade de digitação das placas de licença dos veículos no sistema operacional das praças-dificuldade de visualização das mesmas em situação de filas, e (ii) desgaste físico dos arrecadadores. O desgaste físico dos arrecadadores também se apresenta como provável causa para o aumento de tempos mínimos de atendimento para fluxos altos. Quanto mais intenso o fluxo de veículos, menores são os tempos máximos de atendimento. Quanto maior o fluxo de veículos nas praças, as modelagens indicam uma maior estabilidade do sistema com relação ao processamento de veículos.

Utilização de imagens de satélite e modelagem numérica para o estudo da dispersão de poluentes nas usinas termoelétricas de Charqueadas e São Jerônimo

O objetivo geral deste trabalho é desenvolver estudos relacionados à dispersão de poluentes, considerando a queima de carvão para geração de energia nas Usinas Termoelétricas de Charqueadas e São Jerônimo. O período de estudo foi do dia 17 a 23 de junho de 2003. Neste período houve a passagem de um sistema frontal. Sendo possível avaliar a dispersão dos poluentes em condições pré-frontal, frontal e pós-frontal. Para simular o comportamento dos poluentes neste período, foi utilizada uma subrotina de dispersão acoplada ao modelo RAMS (Regional Atmospheric Modeling System). Os resultados mostraram que nos dias classificados como pré-frontal e pós-frontal as concentrações do material particulado, dióxido de enxofre e do óxido de nitrogênio, atingiram seus valores máximos, pelo fato da umidade relativa do ar estar bastante baixa em torno de 60%, pressão atmosférica da ordem de 1021 hPa e a intensidade dos ventos fraca. No dia classificado como frontal, as concentrações estavam praticamente nulas, devido à passagem do sistema frontal que causou a queda na pressão atmosférica, aumento da umidade relativa e também pelo fato da ocorrência de precipitação atmosférica neste dia. As comparações dos resultados simulados, com os dados observados na estação de qualidade do ar mostraram-se satisfatórios. Com exceção do dia 22 de junho, que apresentou uma diferença da ordem de 90%.

Ensaios sobre emissões de corporate bonds (Debêntures) no mercado brasileiro

O objetivo da presente tese é analisar as três questões não respondidas ou tratadas na literatura nacional - os efeitos de rating, as medidas e proxies de liquidez e os desenhos de contratos nas emissões de diferentes ratings -, que ajudam a compreender a formação de taxa de juros de emissão de debêntures. Para alcançar esse objetivo, esta tese está estruturada em três ensaios empíricos. No ensaio de rating, encontraram-se evidências de que: i) o rating afeta o spread independentemente do indexador da emissão; ii) a origem - nacional ou internacional - dos ratings não é relevante; iii) a causalidade rating-spread não é clara nas emissões que ocorrem em ambiente favorável, em que se verifica maior incidência de classificações discrepantes de risco; iv) a expectativa do mercado internacional com relação ao ambiente econômico brasileiro, o tipo de setor e o volume da emissão são importantes variáveis de controle na determinação do spread. No ensaio de liquidez, foram obtidos os seguintes resultados: i) o tamanho de emissão e determinados tipos de emissores, como, por exemplo, setores de energia e petróleo, são proxies de liquidez; ii) controlando determinados tipos de emissores, as debêntures com maior tamanho de emissão são as mais líquidas; iii) a relação entre idade e liquidez não é clara; iv) a diferença entre preços máximos e mínimos das transações não é uma medida de liquidez apropriada. Por fim, no ensaio de contrato, foram detectadas diferenças entre as cláusulas de emissões de diferentes ratings; conseqüentemente, a padronização afeta a taxa de juros das emissoras, uma vez que o rating é um dos principais determinantes da taxa de juros de emissão. Além disso, a padronização pode não ser ótima nas emissões de baixo rating, pois os contratos-padrão não incluem cláusulas restritivas, tais como repactuação programada e garantia, que reduzem o custo de agência entre acionistas e debenturistas.

Análise do atraso fim-a-fim observado por diferentes classes de tráfego em redes Ethernet micro-segmentadas

A tecnologia de rede dominante em ambientes locais é ethernet. Sua ampla utilização e o aumento crescente nos requisitos impostos pelas aplicações são as principais razões para a evolução constante presenciada desde o surgimento deste tipo de rede simples e barata. Na busca por maior desempenho, configurações baseadas na disputa por um meio de transmissão compartilhado foram substituídas por topologias organizadas ao redor de switches e micro-segmentadas, ou seja, com canais de transmissão individuais para cada estação. Neste mesmo sentido, destacam-se os dispositivos de qualidade de serviço padronizados a partir de 1998, que permitem a diferenciação de tráfego prioritário. O uso desta tecnologia em ambientes de automação sempre foi refreado devido à imprevisibilidade do protocolo CSMA/CD. Entretanto, o uso de switched ethernet em conjunto com a priorização de tráfego representa um cenário bastante promissor para estas aplicações, pois ataca as duas fontes principais de indeterminismo: colisões e enfileiramento. Este trabalho procura avaliar esta estrutura em ambientes de controle distribuído, através de uma análise temporal determinística. Como resultado, propõe-se um modelo analítico de cálculo capaz de prever os valores máximos, médios e mínimos do atraso fim-a-fim observado nas comunicações de diferentes classes de tráfego. Durante a análise, outras métricas relevantes são investigadas como a variação no atraso (jitter) e a interferência entre classes de prioridades distintas.

Estudo do transporte de poluentes na região do Polo Petroquímico, localizado em Triunfo-RS

Este trabalho estuda o transporte e a dispersão de SO2 (dióxido de Enxofre) e Material Particulado(PTS), emitidos pelas várias fontes da COPESUL, no Polo Petroquímico localizado em Triunfo, região metropolitana de Porto Alegre - RS, utilizando o modelo ISCST - Industrial Source Complex Term - Short Term, recomendado pela EPA para tratamento de dispersão. Neste estudo são apresentadas simulações para o período de 5 dias do ano de 2001, dias que apresentaram maiores concentrações ao nível do solo coletados no ponto receptor, sendo possível observar a dispersão da pluma de contaminantes e a localização dos máximos de concentração. Através do MNT(Modelo Numérico de Terreno), foi observado que a região analisada é relativamente plana, fazendo com que o modelo mostre resultados bem próximos a realidade, estando as fontes de emissão estudadas localizadas no centro do terreno. Os resultados indicam que os valores de NOx para o período estudado são os principais responsáveis pelo problema de poluição do ar na região, estes ultrapassaram em 2/5 dos casos o padrão secundário de 190 μg/m3. Os valores verificados foram de 198,65 μg/m3 no horário de 11 horas do dia 26/09/2001, e 218,59 μg/m3 no horário de 10 horas do dia 13/08/2001, estando os outros dentro dos limites exigidos pelo controle de qualidade do ar. Já para os poluentes SO2 e MP os padrões primário e secundário não foram ultrapassados em nenhum dos dias selecionados. Verificou-se que na posição em que a estação de recepção foi instalada, esta estava , na maioria dos dias fora da área modelada pela pluma de dispersão. Portanto os valores medidos não se referem as fontes da COPESUL e sim são provenientes de outras fontes de emissão da região do Polo Petroquímico do Sul como as indústrias de 2ª geração.

Estudo do comportamento mecânico de um solo laterítico estabilizado com cal, aplicado à pavimentação

Esta dissertação apresenta resultados de um estudo laboratorial sobre o comportamento mecânico de um solo laterítico estabilizado com cal cálcica e com cal dolomítica. Este tipo de mistura foi empregado em um pavimento experimental construído na rodovia BR-377/RS, próximo à cidade de Cruz Alta, no final dos anos 60. O estudo incluiu a caracterização do solo (Latossolo Vermelho), além de ensaios de compactação, difratometria de raios-X, compressão simples, tração por compressão diametral, tração na flexão, módulo de resiliência e durabilidade, do solo e das misturas solo-cal. Determinaram-se as influências do teor e tipo de cal e do tempo de cura, nas resistências à compressão simples e à tração e no módulo de resiliência. Estudou-se também a durabilidade de algumas misturas a ciclos de molhagem e secagem, e o efeito da demora na compactação na resistência à compressão simples. Adotaram-se teores de 3%, 4% e 5% de cal, e tempos de cura de até 168 dias. Os resultados da dosagem mostraram que o solo é reativo à cal, sendo necessário um teor de 3% de cal para sua estabilização. A mistura de solo + 4% de cal cálcica, (a mesma empregada no pavimento experimental), após 28 dias de cura, apresentou elevadas resistências à compressão simples (1.519 kPa), à tração na compressão diametral (216 kPa), bem como elevados módulos de resiliência (10.772 MPa) Contudo, observou-se que após terem atingido valores máximos, a resistência e a rigidez das misturas solo-cal diminuíram com tempos de cura mais longos, possivelmente devido à formação de CaCO3 e/ou produtos expansivos (etringita e taumasita). Complementarmente, realizou-se uma análise paramétrica de pavimentos com camadas de base e/ou sub-base de solo-cal, determinando-se curvas que relacionam a vida de fadiga da camada cimentada com a sua espessura e módulo de resiliência. Observou-se que a adoção de bases e sub-bases estabilizadas faz com que os revestimentos asfálticos trabalhem apenas à compressão, não sofrendo ruptura por fadiga. No global, as misturas solo-cal estudadas apresentaram comportamento aceitável para utilização em camadas de base e/ou sub-base de pavimentos, constituindo-se em alternativa para pavimentação, especialmente em regiões carentes de agregados.

Comportamento mecânico de britas empregadas em pavimentação

Camadas granulares têm um importante papel no comportamento global de estruturas de pavimentos, especialmente quando possuem revestimentos delgados ou não possuem revestimento. Para estabelecer métodos de projeto mais eficientes e critérios construtivos, é necessário que a resposta das camadas granulares sob a ação do tráfego seja bem compreendida e levada em consideração. Assim, esta dissertação relata resultados de uma pesquisa sobre o comportamento quanto à resistência ao cisalhamento e deformações permanentes de três britas (uma graduada e duas uniformes com diferentes tamanhos máximos de partícula). Foram realizados ensaios triaxiais estáticos e dinâmicos obtendo-se os parâmetros de resistência de Mohr-Coulomb, módulos de Young, módulos de resiliência e o comportamento quanto a deformações permanentes. Para a brita graduada, desenvolveu-se um estudo complementar sobre o efeito do grau de compactação (90% e 100%) no seu comportamento mecânico. Uma observação conjunta dos resultados estáticos e dinâmicos associada a uma análise mecanística proporciona uma visão global dos resultados obtidos. A brita graduada, com grau de compactação de 100%, mostrou-se um agregado de excelente qualidade, apresentando bom comportamento quanto a deformações permanentes e elevada resistência ao cisalhamento, porém os seu módulos de resiliência foram apenas razoáveis Surpreendentemente, a granulometria uniforme com menor “tamanho máximo de agregado” apresentou um bom comportamento mecânico, porém devido aos seus baixos módulos de resiliência, não pode ser utilizada como camada de base. O seu uso fica restrito a pavimentos permeáveis sob tráfego de veículos leves, conforme demonstrado pelos resultados do monitoramento de um estacionamento experimental analisado nesta dissertação. As britas: uniforme com maior “tamanho máximo de agregado”(GU2) e graduada (GG1) com grau de compactação de 90%, apresentaram baixa resistência ao cisalhamento e alta deformabilidade plástica, comprometendo o seu uso como camada de base. De forma geral, foi confirmada a influência do grau de compactação e da granulometria no comportamento mecânico de materiais granulares. Nem o Índice de Suporte Califórnia e nem a resistência ao cisalhamento mostraram-se bons indicativos do comportamento de materiais granulares quanto a deformações permanentes.

Ecologia e aplicabilidade de métodos para avaliação da produção primária de Zizaniopsis bonariensis : uma macrófita aquática emergente, no sistema hidrológico do Taim, RS

Este estudo analisou variações mensais na dinâmica populacional, biomassa aérea, taxas de crescimento e produção primária de Zizaniopsis bonariensis, uma macrófita aquática emergente, no Sistema Hidrológico do Taim, RS, durante o ano de 2004. Indivíduos fotossinteticamente ativos tiveram todas suas folhas etiquetadas, quantificações de incrementos em biomassa foram determinadas de modo não destrutivo através de quadros experimentais “in situ”. Variáveis alométricas foliares como comprimento de tecido clorofilado e senescente e diâmetro basal de centenas de indivíduos foram mensuradas em campo e convertidas em biomassa através da regressão da área da superfície foliar contra biomassa. Devido à emergência contínua de indivíduos de Z.bonariensis durante todo o ano, coortes mensais foram recrutadas e monitoradas. Variações na biomassa desta espécie também foram avaliadas de modo destrutivo através de blocos de colheita. Os dados obtidos foram utilizados para calcular a Produção Primária Aérea Anual (PPAA) com base em nove métodos de estimativa: Pico Anual de Biomassa; Pico Anual de Biomassa modificado; Milner & Hughes (IBP); Smalley; Valiela et al; Curva de Allen; Soma dos Máximos de Biomassa; Soma de Incrementos e Crescimento Instantâneo. A taxa media mensal de emergência de brotos foi de 11,7ind m-2, sem evidências de pulsos de emergência em função de mudanças sazonais. A densidade média foi de 48 ind m-2. A longevidade das coortes foi estimada em 10,5 meses, e taxa de mortalidade aumentou após o período reprodutivo, sugerindo preliminarmente influência da sazonalidade sobre estas variáveis A biomassa média das coortes aumentou durante o período reprodutivo, demonstrando alto investimento reprodutivo, sendo o maior valor registrado em outubro, 39,96 gPS m-2, pela coorte de janeiro. A taxa de crescimento relativo se mostrou maior durante o primeiro mês após emergência entre todas as coortes, 0,048 g PSm-2d-1. A biomassa clorofilada dos quadros experimentais decresceu no inverno, voltando a aumentar durante a primavera. A PPAA foi estimada em 2870 gPS m-2. Dos métodos utilizados, apenas o método de Milner & Hughes (não-destrutivo) se mostrou adequado, enquanto os baseados em dados de colheita (destrutivo) subestimaram a produção real de Z. bonariensis.

Estudo de caracterização e avaliação da influência do tempo na redução carbotérmica do pó de aciaria elétrica

O pó de aciaria elétrica (PAE) é um resíduo sólido gerado na coleta do material particulado no processo de fabricação do aço em forno elétrico a arco. A classificação segundo norma NBR 10004 para o mesmo é classe I – Perigoso, porque apresenta metais pesados (chumbo e cádmio) em sua constituição, que lixiviam em água excedendo os limites máximos estabelecidos por esta norma. Os elementos ferro e zinco são os constituintes majoritários presentes no PAE. O conhecimento da forma como estes elementos estão combinados é de suma importância antes do estudo de novas tecnologias para processamento do resíduo. Este trabalho teve como objetivos principais caracterizar química, física, estrutural e morfologicamente o PAE e avaliar a influência do tempo na redução carbotérmica deste resíduo. Na etapa de caracterização do PAE as seguintes técnicas foram utilizadas: análise granulométrica, análise química, microscopia eletrônica de varredura (MEV), análise por dispersão de energia (EDS), análise digitalizada de mapeamento de raios-x, difração de raios-x (DRX) e espectroscopia Mössbauer. Os experimentos de redução foram realizados em forno mufla na temperatura de 1050°C, em três tempos de redução: 15, 30 e 60 minutos. O agente redutor utilizado foi um carvão gaúcho, proveniente da mina de Faxinal, beneficiado e contendo 17% de cinza. A quantidade de carbono utilizada nos experimentos foi de 50% em excesso em relação a quantidade estequiométrica. O grau de metalização do ferro das amostras reduzidas foi avaliado via técnicas de análise química via-úmida e espectroscopia Mössbauer. O grau de remoção do zinco nas amostras reduzidas foi avaliado via técnica de análise química por Espectrometria de Emissão Atômica por Plasma – ICP. Através da etapa de caracterização verificou-se que o ferro e o zinco são os elementos majoritários presentes no PAE, e estão combinados na forma dos seguintes óxidos: Fe3O4, ZnFe2O4, FeCr2O4, Ca0,15Fe2,85O4 e ZnO. Nos ensaios de redução verificou-se que inicialmente ocorre a redução dos óxidos de ferro. O máximo grau de metalização obtido no tempo de 60 minutos foi em torno de 50%. A análise estatística mostrou que objetivando maximizar a variável de resposta (grau de metalização) e garantir economia o tempo de 15 minutos se mostrou o mais adequado. Para as fases com zinco, verificou-se que após a decomposição da ZnFe2O4 e tempos maiores de redução teve início a redução da ZnO. Com 60 minutos na temperatura de 1050°C o grau de remoção do zinco foi de 95%.

Derivações de Shamsuddin simples de K[X1,...,Xn] e ideais maximais cíclicos à esquerda da álgebra de Weyl An(k)

Seja K um corpo de característica zero e seja An(K) a n-ésima Álgebra de Weyl sobre K. Neste trabalho, discutimos a existência de ideais maximais à esquerda de An (K) gerados por operadores de ordem 1. Primeiramente, estabelecemos uma relação entre derivações simples de K[X1: ..., Xn] e ideais principais maximais à esquerda de Ân(K). Para n> 2, caracterizamos as derivações de Shamsuddin de K[X1, ..., Xn] que são simples. Depois, mostrámos que se d é uma derivação de Shamsuddin simples de K[X1, ..., Xn], então existe 9 E K[X1, ..., Xn] tal que Ân.(d+g) é um ideal maximal principal à esquerda.

Representações reduzidas por redes neurais com aprendizado local não-supervisionado sobre campos receptivos restritos

Algoritmos ótimos na extração de componentes principais com aprendizado não-supervisionado em redes neurais de múltiplos neurônios de saída são não-locais, ou seja, as modificações em uma dada sinapse entre dois neurônios dependem também da atividade de outros neurônios. Esta rede ótima extrairá as principais componentes dos dados e submetidos à sua primeira camada. As principais componentes são as projeções destes vetores nos autovalores máximos da matriz de correlação Gij = (eiej), onde a média (-) é sobre a distribuição de e. Existem fortes evidências indicando que sinapses biológicas só se modificam via regras locais, como por exemplo a regra de Hebb. Mas se aplicarmos regras locais numa rede com múltiplas saídas, todos os neurônios da saída serão equivalentes e darão respostas redundantes. A rede será bastante ineficiente. Um modo de contornar este problema é através da restrição dos campos receptivos dos neurônios de saída. Se cada neurônio acessar diferentes partes dos estímulos de entrada, a redundância diminui significativamente. Em contrapartida, ao mesmo tempo que a redundância diminui, também diminui a informação contida em cada neurônio; assim, devemos balancear os dois efeitos otimizando o campo receptivo. O valor ótimo, em geral, depende da natureza dos estímulos, sua estatística, e também do ruído intrínseco à rede. Objetivamos com este trabalho determinar a estrutura ótima de campos receptivos com aprendizado não-supervisionado para uma rede neural de uma camada em diversas condições medindo seu desempenho a partir de técnicas de reconstrução.

Resistência ao cisalhamento de solos não saturados : análise experimental e teórica

Um equipamento triaxial convencional originalmente utilizado para ensaios em amostras de solos saturados saturada foi modificado para que fosse possível a realização de ensaios triaxiais em solos não saturados. As principais modificações foram a instalação de uma pedra cerâmica de alto valor de entrada de ar (300 kPa) na nova base da câmara triaxial e uma nova linha de aplicação de pressão para aplicação de pressão de ar na amostra. A técnica de translação de eixos é utilizada para evitar cavitação no sistema de aplicação de pressão. O solo estudado nesta pesquisa consiste de um colúvio de arenito localizado nas encostas da Formação Serra Geral, entre os municípios de Timbé do Sul (SC) e São José dos Ausentes (RS). Ensaios triaxiais convencionais foram realizados nas condições drenado e não drenado em amostras indeformadas para tensões de confinamento iniciais variando de 50 até 500 kPa; ensaios triaxiais convencionais foram realizados na condição drenada em amostras remoldadas para as tensões de confinamento iniciais de 50, 100 e 200kPa. Ensaios triaxiais com sucção controlada foram realizados na condição drenado em amostras indeformadas para as tensões normais líquidas de 50, 100, 150 e 200 kPa em níveis de sucção variando de 25 até 150 kPa. A partir dos resultados dos ensaios triaxiais determinou-se os parâmetros de resistência ao cisalhamento: c’ – intercepto coesivo efetivo; f´ - ângulo de atrito interno efetivo; fb – ângulo de atrito interno que quantifica a contribuição da sucção na resistência ao cisalhamento. O comportamento encontrado foi de uma envoltória de resistência ao cisalhamento não saturada bilinear com acréscimo na resistência ao cisalhamento para baixos níveis de sucção até aproximadamente 75 kPa, e após um decréscimo na resistência ao cisalhamento. A tensão normal líquida não influencia no valor de fb e a sucção não provoca alteração no valor de f´. O valor de fb apresenta valores superiores ao valor de f´, quando se utiliza o critério de ruptura clássico de resistência ao cisalhamento. Alguns critérios de determinação dos valores máximos de ruptura foram aplicados para um melhor entendimento do comportamento resistente deste material.

Determinantes da liquidez de corporate bonds no mercado brasileiro

O objetivo desta dissertação é investigar o mercado secundário de debêntures do Brasil, para responder quais as características dos títulos afetam sua liquidez e quais as características de liquidez podem ser observadas nas debêntures brasileiras. Cinco medidas de liquidez foram utilizadas: número de dias que ocorreram transações, número de transações, volume relativo de transações em relação ao montante emitido, diferença entre os preços máximos e mínimos transacionados e a volatilidade do rendimento. Para cada medida de liquidez, verificou-se a influência de oito características das debêntures: rating, volume emitido, prazo de vencimento, segmento do emissor, listagem em bolsa, idade da emissão e tipo de emissão (incentivada e sob instrução de esforços restritos). Foram coletadas 998 emissões públicas de debêntures e suas respectivas transações até 18 meses após a emissão, no período de janeiro de 2007 a agosto de 2015. A base de dados, que somou 53.085 observações, fundamentou-se nas cotações de mercado fornecidas diariamente pelo Sistema Nacional de Debêntures. Como resultado, verificou-se que o volume da emissão, tipo de emissão (incentivada ou restrita) e determinados segmentos são variáveis de liquidez. Adicionalmente constatou-se que, controlando os segmentos dos emissores, debêntures com maior volume emitido são mais líquidas. E mais, a relação entre idade e liquidez não é clara e a diferença entre preços máximos e mínimos das transações não é uma medida de liquidez apropriada. Por fim, verificou-se que a grande concentração de títulos emitidos sob esforços restritos reduziu a liquidez do mercado em comparação com o estudo de Sheng e Saito (2008), apesar do aumento do volume emitido no período. Em contrapartida, a emissão de títulos incentivados elevou o nível de transações no mercado secundário.

Impactos das mudanças climáticas na produção de açúcar e seus efeitos sobre a economia do Brasil

As mudanças climáticas irão impactar as produções agrícolas de maneiras distintas em cada região do mundo. Os impactos sobre culturas açucareiras (cana-de-açúcar e beterraba açucareira) provocarão mudanças na oferta de açúcar com consequentes mudanças econômicas nos países produtores. Para compreender efeitos dos choques de produtividade no setor sucroenergético e na economia brasileira (PIB e bem-estar) realizou-se uma revisão da literatura sobre os impactos previstos nestas culturas para cenários de mudanças climáticas. Levantaram-se 21 trabalhos, divididos em 10 regiões do mundo de acordo com seus padrões de produção e consumo, que permitiram estabelecer três cenários de choques em produtividade, mínimo, médio e máximo, para cada região. Levantaram-se também choques de produtividade mínimos e máximos para as grandes culturas (trigo, milho, arroz e soja) com base em trabalhos divulgados pelo IPCC (Intergovernmental Panel on Climate Change). Para traduzir os choques de produtividade em retornos econômicos, foram simulados em um modelo econômico os três cenários de choques em produtividade sobre as culturas açucareiras, e outros seis cenários combinando tais choques aos choques em produtividade sobre as grandes culturas. O modelo econômico de equilíbrio geral computável GTAP (Global Trade Analysis Project) foi utilizado, considerando o ano base de 2011, para analisar os efeitos sobre o setor sucroenergético e sobre a economia brasileira. Observou-se que as mudanças climáticas tendem a promover ganhos de produtividade nas culturas açucareiras e na produção de açúcar de diversas regiões do globo. No caso do Brasil, os cenários de choques de produtividade sobre as grandes culturas e culturas açucareiras em conjunto promoveram ganhos de PIB e bem-estar pouco superiores aos cenários de choques somente nas culturas açucareiras. Os resultados sugerem que as mudanças climáticas pouco influenciarão os setores de cana-de-açúcar e açúcar do Brasil. A participação destes setores no PIB é pequena de forma que as variações que promoverão para o PIB e bem-estar brasileiro serão modestas, embora positivas. Para futuros trabalhos sugere-se a incorporação de novas culturas e da pecuária aos cenários, estudos mais regionalizados sobre os impactos do clima futuro nas produtividades agrícolas, e a adoção de modelos econômicos dinâmicos nas análises de cenários de choques de produtividade.