Alguns aspectos da modelação ecológica na costa portuguesa (Aveiro) : efeitos físicos na distribuição de nutrientes e biomassa fitoplanctónica

O presente trabalho propõe estudar alguns aspectos da modelação ecológica na Costa Portuguesa, nomeadamente o efeito físico na distribuição de nutrientes e biomassa fitoplanctónica. O principal propósito foi implementar e validar um modelo acoplado tridimensional físico e ecológico, para a costa portuguesa, e aplicá-lo numa área limitada, a norte, por Vila do Conde e a sul pela Figueira da Foz centrada na região de Aveiro, para estudar a distribuição, vertical e horizontal, de temperatura, nutrientes e biomassa fitoplanctónica. A região em estudo está situada na costa oeste da Península Ibérica e faz parte da Região de Afloramento do Atlântico Norte. É caracterizada por condições meteorológicas onde a predominância e prevalência, durante uma grande parte do ano, de ventos de norte/noroeste, constitui um dos principais elementos forçadores do transporte para o largo das águas costeiras e consequente subida das águas mais frias e profundas, ricas em nutrientes. A estas condições juntam-se as boas condições de luminosidade e temperatura essenciais ao desenvolvimento fitoplanctónico, que servirá de alimento às espécies marinhas. Este facto, torna esta região, uma zona de elevada riqueza biológica favorável ao desenvolvimento de várias espécies marinhas, transformando este local num ecossistema de forte produtividade. O modelo foi calibrado e validado para a área em estudo e simulou com sucesso, a resposta do sistema à situação de ventos favoráveis ao afloramento costeiro para a região em estudo. Quando comparada com dados observados, os resultados mostram que o modelo é capaz de prever satisfatoriamente as distribuições superficiais e na coluna de água: da temperatura, dos nutrientes, do oxigénio e da clorofila-a. Os resultados evidenciam o crucial papel desempenhado pelos processos físicos no aumento de fitoplâncton que ocorre ao longo de uma estreita área da costa norte portuguesa, mostrando a estreita ligação entre a distribuição costeira de fitoplâncton e a distribuição costeira de temperatura, à superfície. A produtividade do fitoplâncton que ocorre nas águas costeiras não só é atribuída à disponibilidade de nutrientes mas também à intensidade luminosa. A luz e a intensidade dos ventos de norte/noroeste são os factores chave no controlo dos blooms de fitoplâncton observados nesta região de afloramento, sobretudo no Verão, nomeadamente na camada de mistura pouco profunda e junto à nutriclina. Deste modo os modelos numéricos em associação com dados in situ e imagens de satélite poderão ser considerados uma excelente ferramenta para a análise e previsão de cenários, presentes e futuros, de acções praticadas sobre o meio ambiente.

Modelos autoregressivos por limiares em séries de contagem

A modelação e análise de séries temporais de valores inteiros têm sido alvo de grande investigação e desenvolvimento nos últimos anos, com aplicações várias em diversas áreas da ciência. Nesta tese a atenção centrar-se-á no estudo na classe de modelos basedos no operador thinning binomial. Tendo como base o operador thinning binomial, esta tese focou-se na construção e estudo de modelos SETINAR(2; p(1); p(2)) e PSETINAR(2; 1; 1)T , modelos autorregressivos de valores inteiros com limiares autoinduzidos e dois regimes, admitindo que as inovações formam uma sucessão de variáveis independentes com distribuição de Poisson. Relativamente ao primeiro modelo analisado, o modelo SETINAR(2; p(1); p(2)), além do estudo das suas propriedades probabilísticas e de métodos, clássicos e bayesianos, para estimar os parâmetros, analisou-se a questão da seleção das ordens, no caso de elas serem desconhecidas. Com este objetivo consideraram-se algoritmos de Monte Carlo via cadeias de Markov, em particular o algoritmo Reversible Jump, abordando-se também o problema da seleção de modelos, usando metodologias clássica e bayesiana. Complementou-se a análise através de um estudo de simulação e uma aplicação a dois conjuntos de dados reais. O modelo PSETINAR(2; 1; 1)T proposto, é também um modelo autorregressivo com limiares autoinduzidos e dois regimes, de ordem unitária em cada um deles, mas apresentando uma estrutura periódica. Estudaram-se as suas propriedades probabilísticas, analisaram-se os problemas de inferência e predição de futuras observações e realizaram-se estudos de simulação.

Population genetics of C. Maenas : oceanography and larval dispersal

Decifrar a complexa interacção entre os ciclos de vida de espécies marinhas e a oceanografia revela-se fundamental para a compreensão do fluxo genético e da conectividade no meio marinho. Nas espécies marinhas com desenvolvimento indirecto o fluxo de genes entre populações depende da distância que separa as populações, bem como da interacção entre a duração do desenvolvimento larvar, do comportamento das larvas e dos padrões de circulação oceânica. A conectividade larvar influencia uma variedade de processos como a dinâmica de stocks e de populações, a distribuição e limites geográficos das espécies, a estrutura genética das populações e a dispersão de espécies invasivas e reveste-se consequentemente de uma importância fundamental na identificação das unidades populacionais evolucionariamente relevantes e para a gestão e conservação marinhas. Os marcadores genéticos e os Modelos Individuais Acoplados a Modelos Físico-Biológicos (“ICPBMs”) são actualmente ferramentas fundamentais para o estudo dos padrões de dispersão larvar e para avaliar o nível de conectividade populacional. A presente tese respeita à avaliação das escalas espaciais de conectividade de populações de uma espécie costeira, o caranguejo Carcinus maenas, e utiliza conjuntamente informação de marcadores genéticos, análise de séries temporais de fornecimento de larvas e um modelo numérico de circulação oceânica. O primeiro capítulo introduz a temática da conectividade em espécies marinhas e inclui algumas referências aos métodos moleculares, analíticos e de modelação seguidos ao longo da tese. Através da utilização de múltiplas ferramentas – avaliação da estrutura genética geográfica de C. maenas na sua distribuição nativa com recurso a marcadores de DNA (microssatélites) (Capítulo 2), avaliação da estrutura genética temporal das larvas que formam os eventos de fornecimento larvar à Ria de Aveiro, NW Portugal (Capítulo 3), descrição da variabilidade inter-anual do fornecimento larvar à Ria de Aveiro, NW Portugal (Capítulo 4) e validação de um modelo ICPBM que descreve os padrões observados de fornecimento (Capítulo 5) – esta tese espera poder contribuir para uma melhor compreensão dos mecanismos que regulam o fluxo de genes e a conectividade entre populações de organismos marinhos. No Capítulo 6 são apresentadas as principais conclusões da investigação. A análise genética com recurso a microssatélites indicou que as populações de C. maenas são geneticamente homogéneas ao longo de várias centenas de km, dentro da distribuição nativa da espécie. Paralelamente, não foram encontrados indícios da existência de reprodução por “sweepstakes” em C. maenas de populações da costa oeste da Península Ibérica, visto que não se obtiveram diferenças genéticas significativas entre os eventos larvares. Também não se encontrou qualquer estrutura familiar entre as larvas que formam cada episódio de fornecimento, e não houve nenhuma redução significativa da variabilidade genética das larvas quando comparada com a de caranguejos adultos. A análise de séries temporais de suprimento de larvas na Ria de Aveiro em cinco anos estudados indica que este é um fenómeno episódico e variável, sendo os maiores episódios de fornecimento coincidentes com as marés vivas e acentuados por fortes ventos de sul. O modelo ICPBM foi validado com sucesso e parece fornecer uma estimativa realística das escalas espaciais e temporais de dispersão larvar, de acordo com as observações da estrutura genética e da ausência de reprodução por “sweepstake” em C. maenas da costa oeste da Península Ibérica

Studies on magnetocaloric and magnetic coupling effects

O presente trabalho apresenta novas metodologias desenvolvidas para a análise das propriedades magnéticas e magnetocalóricas de materiais, sustentadas em considerações teóricas a partir de modelos, nomeadamente a teoria de transições de fase de Landau, o modelo de campo médio molecular e a teoria de fenómeno crítico. São propostos novos métodos de escala, permitindo a interpretação de dados de magnetização de materiais numa perspectiva de campo médio molecular ou teoria de fenómeno crítico. É apresentado um método de estimar a magnetização espontânea de um material ferromagnético a partir de relações entropia/magnetização estabelecidas pelo modelo de campo médio molecular. A termodinâmica das transições de fase magnéticas de primeira ordem é estudada usando a teoria de Landau e de campo médio molecular (modelo de Bean-Rodbell), avaliando os efeitos de fenómenos fora de equilíbrio e de condições de mistura de fase em estimativas do efeito magnetocalórico a partir de medidas magnéticas. Efeitos de desordem, interpretados como uma distribuição na interacção magnética entre iões, estabelecem os efeitos de distribuições químicas/estruturais nas propriedades magnéticas e magnetocalóricas de materiais com transições de fase de segunda e de primeira ordem. O uso das metodologias apresentadas na interpretação das propriedades magnéticas de variados materiais ferromagnéticos permitiu obter: 1) uma análise quantitativa da variação de spin por ião Gadolínio devido à transição estrutural do composto Gd5Si2Ge2, 2) a descrição da configuração de cluster magnético de iões Mn na fase ferromagnética em manganites da família La-Sr e La-Ca, 3) a determinação dos expoentes críticos β e δ do Níquel por métodos de escala, 4) a descrição do efeito da pressão nas propriedades magnéticas e magnetocalóricas do composto LaFe11.5Si1.5 através do modelo de Bean-Rodbell, 5) uma estimativa da desordem em manganites ferromagnéticas com transições de segunda e primeira ordem, 6) uma descrição de campo médio das propriedades magnéticas da liga Fe23Cu77, 7) o estudo de efeitos de separação de fase na família de compostos La0.70-xErxSr0.30MnO3 e 8) a determinação realista da variação de entropia magnética na família de compostos de efeito magnetocalórico colossal Mn1-x-yFexCryAs.

High-resolution seismic reflection investigation of gas accumulation and seepage in the tidal channels of the Ria of Aveiro barrier lagoon (Portugal)

O metano é um gás de estufa potente e uma importante fonte de energia. A importância global e impacto em zonas costeiras de acumulações e escape de gás metano são ainda pouco conhecidas. Esta tese investiga acumulações e escape de gás em canais de maré da Ria de Aveiro com dados de cinco campanhas de reflexão sísmica de alta resolução realizadas em 1986, 1999, 2002 e 2003. Estas incluem três campanhas de Chirp (RIAV99, RIAV02 e RIAV02A) e duas campanhas de Boomer (VOUGA86 e RIAV03). O processamento dos dados de navegação incluíram filtros de erros, correcções de sincronização de relógios de sistemas de aquisição de dados, ajuste de “layback” e estimativa da posição de “midpoint”. O processamento do sinal sísmico consistiu na correcção das amplitudes, remoção de ruído do tipo “burst”, correcções estáticas, correcção do “normal move-out”, filtragem passabanda, desconvolução da assinatura e migração Stolt F-K. A análise da regularidade do trajecto de navegação, dos desfasamentos entre horizontes e dos modelos de superfícies foi utilizada para controlo de qualidade, e permitiu a revisão e melhoria dos parâmetros de processamento. A heterogeneidade da cobertura sísmica, da qualidade do sinal, da penetração e da resolução, no seu conjunto constrangeram o uso dos dados a interpretações detalhadas, mas locais, de objectos geológicos da Ria. É apresentado um procedimento para determinar a escolha de escalas adequadas para modelar os objectos geológicos, baseado na resolução sísmica, erros de posicionamento conhecidos e desfasamentos médios entre horizontes. As evidências de acumulação e escape de gás na Ria de Aveiro incluem turbidez acústica, reflexões reforçadas, cortinas acústicas, domas, “pockmarks” e alinhamentos de “pockmarks” enterradas, horizontes perturbados e plumas acústicas na coluna de água (flares). A estratigrafia e a estrutura geológica controlam a distribuição e extensão das acumulações e escape de gás. Ainda assim, nestes sistemas de baixa profundidade de água, as variações da altura de maré têm um impacto significativo na detecção de gás com métodos acústicos, através de alterações nas amplitudes originais de reflexões reforçadas, turbidez acústica e branqueamento acústico em zonas com gás. Os padrões encontrados confirmam que o escape de bolhas de gás é desencadeado pela descida da maré. Há acumulações de gás em sedimentos Holocénicos e no substrato de argilas e calcários do Mesozóico. Evidências directas de escape de gás em sondagens em zonas vizinhas, mostraram gás essencialmente biogénico. A maioria do gás na área deve ter sido gerado em sedimentos lagunares Holocénicos. No entanto, a localização e geometria de estruturas de escape de fluidos em alguns canais de maré, seguem o padrão de fracturas do substrato Mesozóico, indicando uma possível fonte mais profunda de gás e que estas fracturas funcionam como condutas preferenciais de migração dos fluidos e exercem um controlo estrutural na ocorrência de gás na Ria.

Numerical modelling of radionuclide migration in near-surface systems

Em todo o mundo são usados, hoje em dia, modelos numéricos hidrogeoquímicos para simular fenómenos naturais e fenómenos decorrentes de actividades antrópicas. Estes modelos ajudam-nos a compreender o ambiente envolvente, a sua variabilidade espacial e evolução temporal. No presente trabalho apresenta-se o desenvolvimento de modelos numéricos hidrogeoquímicos aplicados no contexto do repositório geológico profundo para resíduos nucleares de elevada actividade. A avaliação da performance de um repositório geológico profundo inclui o estudo da evolução geoquímica do repositório, bem como a análise dos cenários de mau funcionamento do repositório, e respectivas consequências ambientais. Se se escaparem acidentalmente radionuclídeos de um repositório, estes poderão atravessar as barreiras de engenharia e barreiras naturais que constituem o repositório, atingindo eventualmente, os ecosistemas superficiais. Neste caso, os sedimentos subsuperficiais constituem a última barreira natural antes dos ecosistemas superficiais. No presente trabalho foram desenvolvidos modelos numéricos que integram processos biogeoquímicos, geoquímicos, hidrodinâmicos e de transporte de solutos, para entender e quantificar a influência destes processos na mobilidade de radionuclídeos em sistemas subsuperficiais. Os resultados alcançados reflectem a robustez dos instrumentos numéricos utilizados para desenvolver simulações descritivas e predictivas de processos hidrogeoquímicos que influenciam a mobilidade de radionuclídeos. A simulação (descritiva) de uma experiência laboratorial revela que a actividade microbiana induz a diminuição do potencial redox da água subterrânea que, por sua vez, favorece a retenção de radionuclídeos sensíveis ao potencial redox, como o urânio. As simulações predictivas indicam que processos de co-precipitação com minerais de elementos maioritários, precipitação de fases puras, intercâmbio catiónico e adsorção à superfície de minerais favorecem a retenção de U, Cs, Sr e Ra na fase sólida de uma argila glaciar e uma moreia rica em calcite. A etiquetagem dos radionuclídeos nas simulações numéricas permitiu concluir que a diluição isotópica joga um papel importante no potencial impacte dos radionuclídeos nos sistemas subsuperficiais. A partir dos resultados das simulações numéricas é possivel calcular coeficientes de distribuição efectivos. Esta metodologia proporciona a simulação de ensaios de traçadores de longa duração que não seriam exequíveis à escala da vida humana. A partir destas simulações podem ser obtidos coeficientes de retardamento que são úteis no contexto da avaliação da performance de repositórios geológicos profundos.

Métodos robustos em geoestatística

O objectivo principal da presente tese consiste no desenvolvimento de estimadores robustos do variograma com boas propriedades de eficiência. O variograma é um instrumento fundamental em Geoestatística, pois modela a estrutura de dependência do processo em estudo e influencia decisivamente a predição de novas observações. Os métodos tradicionais de estimação do variograma não são robustos, ou seja, são sensíveis a pequenos desvios das hipóteses do modelo. Essa questão é importante, pois as propriedades que motivam a aplicação de tais métodos, podem não ser válidas nas vizinhanças do modelo assumido. O presente trabalho começa por conter uma revisão dos principais conceitos em Geoestatística e da estimação tradicional do variograma. De seguida, resumem-se algumas noções fundamentais sobre robustez estatística. No seguimento, apresenta-se um novo método de estimação do variograma que se designou por estimador de múltiplos variogramas. O método consiste em quatro etapas, nas quais prevalecem, alternadamente, os critérios de robustez ou de eficiência. A partir da amostra inicial, são calculadas, de forma robusta, algumas estimativas pontuais do variograma; com base nessas estimativas pontuais, são estimados os parâmetros do modelo pelo método dos mínimos quadrados; as duas fases anteriores são repetidas, criando um conjunto de múltiplas estimativas da função variograma; por fim, a estimativa final do variograma é definida pela mediana das estimativas obtidas anteriormente. Assim, é possível obter um estimador que tem boas propriedades de robustez e boa eficiência em processos Gaussianos. A investigação desenvolvida revelou que, quando se usam estimativas discretas na primeira fase da estimação do variograma, existem situações onde a identificabilidade dos parâmetros não está assegurada. Para os modelos de variograma mais comuns, foi possível estabelecer condições, pouco restritivas, que garantem a unicidade de solução na estimação do variograma. A estimação do variograma supõe sempre a estacionaridade da média do processo. Como é importante que existam procedimentos objectivos para avaliar tal condição, neste trabalho sugere-se um teste para validar essa hipótese. A estatística do teste é um estimador-MM, cuja distribuição é desconhecida nas condições de dependência assumidas. Tendo em vista a sua aproximação, apresenta-se uma versão do método bootstrap adequada ao estudo de observações dependentes de processos espaciais. Finalmente, o estimador de múltiplos variogramas é avaliado em termos da sua aplicação prática. O trabalho contém um estudo de simulação que confirma as propriedades estabelecidas. Em todos os casos analisados, o estimador de múltiplos variogramas produziu melhores resultados do que as alternativas usuais, tanto para a distribuição assumida, como para distribuições contaminadas.

Estimação robusta em modelos lineares de equações simultâneas

Os Modelos de Equações Simultâneas (SEM) são modelos estatísticos com muita tradição em estudos de Econometria, uma vez que permitem representar e estudar uma vasta gama de processos económicos. Os estimadores mais usados em SEM resultam da aplicação do Método dos Mínimos Quadrados ou do Método da Máxima Verosimilhança, os quais não são robustos. Em Maronna e Yohai (1997), os autores propõem formas de “robustificar” esses estimadores. Um outro método de estimação com interesse nestes modelos é o Método dos Momentos Generalizado (GMM), o qual também conduz a estimadores não robustos. Estimadores que sofrem de falta de robustez são muito inconvenientes uma vez que podem conduzir a resultados enganadores quando são violadas as hipóteses subjacentes ao modelo assumido. Os estimadores robustos são de grande valor, em particular quando os modelos em estudo são complexos, como é o caso dos SEM. O principal objectivo desta investigação foi o de procurar tais estimadores tendo-se construído um estimador robusto a que se deu o nome de GMMOGK. Trata-se de uma versão robusta do estimador GMM. Para avaliar o desempenho do novo estimador foi feito um adequado estudo de simulação e foi também feita a aplicação do estimador a um conjunto de dados reais. O estimador robusto tem um bom desempenho nos modelos heterocedásticos considerados e, nessas condições, comporta-se melhor do que os estimadores não robustos usados no estudo. Contudo, quando a análise é feita em cada equação separadamente, a especificidade de cada equação individual e a estrutura de dependência do sistema são dois aspectos que influenciam o desempenho do estimador, tal como acontece com os estimadores usuais. Para enquadrar a investigação, o texto inclui uma revisão de aspectos essenciais dos SEM, o seu papel em Econometria, os principais métodos de estimação, com particular ênfase no GMM, e uma curta introdução à estimação robusta.

Contributos para a análise estatística de séries de contagem

A análise das séries temporais de valores inteiros tornou-se, nos últimos anos, uma área de investigação importante, não só devido à sua aplicação a dados de contagem provenientes de diversos campos da ciência, mas também pelo facto de ser uma área pouco explorada, em contraste com a análise séries temporais de valores contínuos. Uma classe que tem obtido especial relevo é a dos modelos baseados no operador binomial thinning, da qual se destaca o modelo auto-regressivo de valores inteiros de ordem p. Esta classe é muito vasta, pelo que este trabalho tem como objectivo dar um contributo para a análise estatística de processos de contagem que lhe pertencem. Esta análise é realizada do ponto de vista da predição de acontecimentos, aos quais estão associados mecanismos de alarme, e também da introdução de novos modelos que se baseiam no referido operador. Em muitos fenómenos descritos por processos estocásticos a implementação de um sistema de alarmes pode ser fundamental para prever a ocorrência de um acontecimento futuro. Neste trabalho abordam-se, nas perspectivas clássica e bayesiana, os sistemas de alarme óptimos para processos de contagem, cujos parâmetros dependem de covariáveis de interesse e que variam no tempo, mais concretamente para o modelo auto-regressivo de valores inteiros não negativos com coeficientes estocásticos, DSINAR(1). A introdução de novos modelos que pertencem à classe dos modelos baseados no operador binomial thinning é feita quando se propõem os modelos PINAR(1)T e o modelo SETINAR(2;1). O modelo PINAR(1)T tem estrutura periódica, cujas inovações são uma sucessão periódica de variáveis aleatórias independentes com distribuição de Poisson, o qual foi estudado com detalhe ao nível das suas propriedades probabilísticas, métodos de estimação e previsão. O modelo SETINAR(2;1) é um processo auto-regressivo de valores inteiros, definido por limiares auto-induzidos e cujas inovações formam uma sucessão de variáveis independentes e identicamente distribuídas com distribuição de Poisson. Para este modelo estudam-se as suas propriedades probabilísticas e métodos para estimar os seus parâmetros. Para cada modelo introduzido, foram realizados estudos de simulação para comparar os métodos de estimação que foram usados.

Ecotoxicological evaluation of pesticide mixtures

Neste trabalho o risco da aplicação de moluscicidas (metaldeído e metiocarbe) para o isópode terrestre Porcellionides pruinosus foi avaliado usando como parâmetros a mortalidade e biomarcadores de exposição. O tempo até à morte dos isópodes (após contacto com os moluscicidas) foi muito curto, especialmente no caso da exposição ao metiocarbe. Os vários biomarcadores revelaram-se úteis para a compreensão do modo de acção dos dois moluscicidas neste isópode, particularmente o efeito do carbamato metiocarbe na inibição da enzima acetilcolinesterase (AChE). Os efeitos de combinações binárias de três produtos de protecção das plantas (PPP) dimetoato, glifosato e espirodiclofeno foram avaliados testando o comportamento de evitamento de P. pruinosus, o sucesso reprodutivo do colêmbolo Folsomia candida e o crescimento das plantas Brassica rapa e Triticum aestivum, usando os dois modelos de referência de concentração de adição (CA) e acção independente (IA). O modelo MIXTOX foi usado para avaliar possíveis desvios (devido a interacções entre os pesticidas) dos dois modelos de referência. Os resultados obtidos permitem constatar que estes PPP quando aplicados segundo a dose recomendada não acarretam efeitos perniciosos para os organismos testados. Foi detectado sinergismo na mistura feita com glifosato e espirodiclofeno no isópode P. pruinosus e na mistura com glifosato e dimetoato na planta T. aestivum. Um ecossistema terrestre em pequena escala (“STEM”) foi desenvolvido, contendo um solo agrícola mediterrânico. Nestes STEM, minhocas (Eisenia andrei), P. pruinosus, B. rapa e “bait-lamina” foram incorporados no sentido de avaliar os efeitos da aplicação de dimetoato com espirodiclofeno e glifosato com dimetoato. A dose recomendada de aplicação dos PPP, quer na exposição individual quer nas misturas binárias não teve quaisquer efeitos nas espécies testadas. As minhocas foram sensíveis à aplicação conjunta de dimetoato com espirodiclofeno (10 vezes a dose recomendada) na sua distribuição vertical ao longo da coluna do STEM, e foi detectado sinergismo (i.e. mais minhocas escaparam do que a predição feita pelo modelo IA). Em todas as misturas binárias feitas com glifosato e dimetoato registou-se um decréscimo no consumo de “bait-lamina”, indicando sinergismo (menos “bait-lamina” consumidos que o esperado). Dos quatro biomarcadores (Catalase, AChE, GST e LPO) avaliados nos isópodes, verificaram-se diferenças significativas na actividade da enzima AChE (quando dimetoato foi aplicado no solo) e LPO (aumento da actividade devido à aplicação de glifosato e dimetoato).

Hydraulics of stepped chutes: experimental-numerical-theoretical study

O estudo de escoamentos turbulentos em descarregadores em degraus tem sido um desafio para os investigadores. A macro-rugosidade do leito, a ondulação da superfície livre, a intermitência da localização da secção inicial de entrada de ar e o escoamento bi-fásico a jusante da secção inicial de entrada de ar fazem com que a caracterização do escoamento deslizante sobre turbilhões em descarregadores em degraus não seja simples. Actualmente, é possível combinar técnicas de medição fiáveis com simulações numéricas e análise teórica. Nesta dissertação, o estudo experimental baseia-se em resultados experimentais obtidos em duas instalações experimentais: a instalação A, do Laboratório Nacional de Engenharia Civil (LNEC), dotada de um descarregador em degraus com declive de 1V:0.75H e a instalação B, do Instituto Superior Técnico (IST), dotada de um descarregador em degraus com declive de 1V:2H. Uma sonda de ar, um tubo de Pitot modificado e vários hidrómetros permitiram o estudo do escoamento deslizante sobre turbilhões mono- e bifásico ao longo do descarregador em degraus da instalação A. Na bacia de dissipação de energia, a observação visual do escoamento foi auxiliada pelas leituras efectuadas em tomadas de pressão localizadas na soleira da bacia e pelas leituras da altura do escoamento efectuadas em réguas graduadas localizadas nas paredes da bacia. Na instalação B foram utilizados um tubo de Pitot e vários hidrómetros para estudar a região não arejada do escoamento deslizante sobre turbilhões no descarregador em degraus. O estudo numérico é baseado em simulações numéricas da região não arejada do escoamento deslizante sobre turbilhões sobre descarregadores em degraus com o código comercial de CFD FLOW-3D® de modo a reproduzir as condições ensaiadas experimentalmente. As simulações beneficiaram da técnica de blocos múltiplos (multi-block) num sistema de coordenadas cartesianas, da determinação da superfície livre pelo método TruVOF e da utilização de dois modelos de turbulência: os modelos k- e RNG k-. Por último, o estudo teórico consistiu em desenvolver um modelo simplificado 1D para determinar as características hidráulicas principais do trecho não arejado do escoamento deslizante sobre turbilhões em descarregadores em degraus. O modelo foi desenvolvido a partir das equações de Navier-Stokes, conjuntamente com resultados experimentais e numéricos. Os resultados apresentados nesta dissertação contribuem para o conhecimento do escoamento deslizante sobre turbilhões em descarregadores em degraus, nomeadamente na região não arejada, na secção inicial de entrada de ar e na região arejada. A hidráulica dos dissipadores de energia a jusante de descarregadores em degraus, em particular das bacias tipo III do USBR, é também objecto de estudo. Em relação à região não arejada do escoamento em descarregadores com declive acentuado, e com base em resultados experimentais e numéricos, são propostas expressões para estimar o desenvolvimento da altura equivalente de água, da espessura da camada limite, da concentração média de ar, do coeficiente de energia cinética, da dissipação de energia, do factor de resistência e do coeficiente n da fórmula de Manning. São ainda propostas expressões adimensionais para a energia cinética turbulenta e sua dissipação. Para declives moderados, são propostas expressões para estimar o desenvolvimento da altura equivalente de água, do coeficiente de energia cinética e da energia específica residual. São ainda apresentados valores do expoente 1/N da expressão adimensional da distribuição de velocidades, quer para descarregadores com declive acentuado quer com declive moderado. Em conformidade com outros estudos centrados em escoamentos de parede e com derivações teóricas, para a região não arejada do escoamento em descarregadores em degraus de acentuado declive, observa-se que o factor de resistência depende da macro-rugosidade criada pelos degraus e da geometria da secção transversal e que o coeficiente n da fórmula de Manning aumenta com a rugosidade. A descrição estatística da turbulência do escoamento é igualmente explorada, contribuindo para o conhecimento da estrutura do escoamento. Observou-se que para números de Reynolds rugoso não superiores a 6.8x104 a energia cinética turbulência e a sua dissipação cumprem leis de semelhança. Estas expressões adimensionais estão de acordo com os resultados obtidos por outros autores para escoamentos completamente desenvolvidos em canais abertos e no escoamento em rios com leito de gravilha. Em acréscimo, a taxa de dissipação de energia, quer para descarregadores de declive acentuado quer de moderado declive, é baixa. Por último, observa-se que os valores da média temporal da concentração de ar entre 0 e 1 medidos na região não arejada do escoamento dizem respeito não só ao ar capturado entre ondas de água, na zona de ondulação da superfície livre, mas também ao ar emulsionado no escoamento, i.e., sob a forma de bolhas de ar, quando perto da secção média inicial de entrada de ar, devido à diferença entre localizações instantânea e média temporal. Foram revistas metodologias e fórmulas para estimar a localização da secção inicial de entrada de ar e apresentadas expressões para estimar a concentração média de ar e a altura equivalente de água nessa secção. Relativamente à região de escoamento arejado em descarregadores em degraus com declive acentuado, os resultados experimentais apresentados nesta dissertação permitiram estimar a influência da definição da superfície livre nos parâmetros hidráulicos da região do escoamento arejado e estimar a máxima elevação do escoamento nesta região do escoamento. Com base nos resultados experimentais obtidos na bacia de dissipação de energia do tipo III do USBR localizada a jusante do descarregador em degraus da instalação A, observou-se que os perfis da altura piezométrica e da altura do escoamento tendem a seguir o perfil recomendado pelo USBR para bacias tipo III. A excepção ocorre à entrada da bacia, onde as alturas piezométricas apresentadas nesta dissertação excedem largamente as apresentadas pelo USBR. É ainda observado que, tal como entre as bacias tipo I e tipo III do USBR, o ressalto hidráulico estabiliza muito mais rapidamente numa bacia tipo III a jusante de um descarregador em degraus do que uma bacia tipo I a jusante do mesmo descarregador em degraus. Finalmente, observa-se que os blocos de amortecimento a colocar no descarregador não têm influência visível nos resultados da altura piezométrica nem da altura do escoamento ao longo da bacia. Relativamente às simulações numéricas do escoamento não arejado, a proximidade entre resultados experimentais e numéricos permite validar o modelo teórico e a integração numérica usados no FLOW-3D®. As simulações desenvolvidas também mostraram que o modelo de turbulência k- permite representar as características do escoamento não arejado em descarregadores em degraus, uma vez que não foram observadas diferenças significativas entre as simulações com este modelo e com o modelo RNG k-. Finalmente, observou-se que o modelo de entrada de ar usado no FLOW-3D® é válido para estimar a localização da secção inicial de entrada de ar. Por último, a proximidade entre os resultados obtidos da aplicação do modelo teórico desenvolvido no âmbito desta dissertação e os resultados experimentais indica que as hipóteses e simplificações consideradas no desenvolvimento do modelo são adequadas.

Location-routing problems of semi-obnoxious facilities: approaches and decision support

A presente tese resulta de um trabalho de investigação cujo objectivo se centrou no problema de localização-distribuição (PLD) que pretende abordar, de forma integrada, duas actividades logísticas intimamente relacionadas: a localização de equipamentos e a distribuição de produtos. O PLD, nomeadamente a sua modelação matemática, tem sido estudado na literatura, dando origem a diversas aproximações que resultam de diferentes cenários reais. Importa portanto agrupar as diferentes variantes por forma a facilitar e potenciar a sua investigação. Após fazer uma revisão e propor uma taxonomia dos modelos de localização-distribuição, este trabalho foca-se na resolução de alguns modelos considerados como mais representativos. É feita assim a análise de dois dos PLDs mais básicos (os problema capacitados com procura nos nós e nos arcos), sendo apresentadas, para ambos, propostas de resolução. Posteriormente, é abordada a localização-distribuição de serviços semiobnóxios. Este tipo de serviços, ainda que seja necessário e indispensável para o público em geral, dada a sua natureza, exerce um efeito desagradável sobre as comunidades contíguas. Assim, aos critérios tipicamente utilizados na tomada de decisão sobre a localização destes serviços (habitualmente a minimização de custo) é necessário adicionar preocupações que reflectem a manutenção da qualidade de vida das regiões que sofrem o impacto do resultado da referida decisão. A abordagem da localização-distribuição de serviços semiobnóxios requer portanto uma análise multi-objectivo. Esta análise pode ser feita com recurso a dois métodos distintos: não interactivos e interactivos. Ambos são abordados nesta tese, com novas propostas, sendo o método interactivo proposto aplicável a outros problemas de programação inteira mista multi-objectivo. Por último, é desenvolvida uma ferramenta de apoio à decisão para os problemas abordados nesta tese, sendo apresentada a metodologia adoptada e as suas principais funcionalidades. A ferramenta desenvolvida tem grandes preocupações com a interface de utilizador, visto ser direccionada para decisores que tipicamente não têm conhecimentos sobre os modelos matemáticos subjacentes a este tipo de problemas.

Estudo óptico de cristais e heterostruturas de CdTe

Neste trabalho privilegiam-se as técnicas ópticas, nomeadamente a fotoluminescência a transmissão, a reflexão e a difusão de Raman na caracterização de defeitos e impurezas em cristais e heterostruturas de CdTe/GaAs. No primeiro capítulo efectua-se uma revisão sucinta das propriedades dos compostos II-VI e são colocados os problemas a investigar. O segundo capítulo é dedicado a alguns aspectos teóricos relevantes para a análise dos resultados obtidos nos capítulos quatro a nove. O equipamento utilizado e as experiências realizadas são descritas no terceiro capítulo. A caracterização por espectroscopia de luminescência do tipo de transições e dos defeitos envolvidos nas amostras de fábrica é efectuada com detalhe no quarto capítulo. Neste capítulo são analisadas ainda as emissões devidas a desvios estequiométricos causados pelo recozimento com excesso e defeito de Cd. No capítulo cinco são estudadas por fotoluminescência amostras dopadas intencionalmente com oxigénio por difusão e mostra-se que este elemento se comporta como trapa isoelectrónica no CdTe tal como acontece no ZnTe. Neste capítulo são estudadas também amostras dopadas com ferro pelo mesmo método e são apresentadas as dificuldades em colocar este dopante em sítios substitucionais, nomeadamente no sítio do Cd. No sexto capítulo é estudada a região 1.4 eV, evidenciando o seu comportamento em função da temperatura e da potência de excitação, o seu perfil e a interacção electrãorede. No sétimo capítulo mostra-se que a técnica micro-Raman com luz visível coerente não permite extrair conclusões fiáveis acerca das inclusões de Te na superfície do CdTe, uma vez que a radiação ao ser focada nas amostras, induz a formação de aglomerados nos quais o Te é o elemento dominante. Neste capítulo calcula-se também a concentração de portadores livres através da interacção do plasmão com o fonão óptico longitudinal. O oitavo capítulo é dedicado ao estudo de camadas de CdTe/GaAs com diferentes espessuras nomeadamente na análise da distribuição das deslocações e da deformação na superfície em função da espessura. Os resultados obtidos são comparados através da largura a meia altura das curvas de DCXRD (“Double Crystal X Ray Diffraction”) e dos espectros de reflectância. Por fotoluminescência, são caracterizados os defeitos introduzidos durante o crescimento, são utilizados e desenvolvidos modelos complementares na distinção do tipo de transições ópticas obtidas. No nono capítulo, por espectroscopia de absorção e de reflexão em cristais de CdTe e em camadas de CdTe/GaAs na região reststrahlen são determinadas as frequências dos modos ópticos longitudinal e transversal com bastante precisão. Os resultados obtidos — pelas relações de dispersão de Kramers- Kronig — são simulados pelo modelo do oscilador harmónico classico, mostrando que ambos os métodos descrevem de forma semelhante o comportamento do CdTe nessa região sendo possível determinar as frequências ópticas transversal e longitudinal, as constantes dieléctrica óptica e estática e os coeficientes de amortecimento e de Szigeti. Nas heterostruturas e nas camadas mais espessas determina-se também a concentração de portadores de carga n. No décimo capítulo resumem-se as conclusões do trabalho e são abordados aspectos relacionados com desenvolvimentos possíveis de futuros trabalhos.

Aplicações das redes de Bragg na biomecânica

O presente trabalho tem como objectivo o estudo, desenvolvimento e aplicações na área da biomecânica de sensores intrínsecos baseados em redes de Bragg em fibras ópticas (FBG). As aplicações são feitas em modelos biomecânicos in vitro tais como: implantes de anca, prótese de joelho, placas de osteossíntese e implantes dentários. A optimização do desenvolvimento de próteses e respectivos elementos de fixação é actualmente dependente da geração e validação experimental de seus modelos computacionais. A validação destes modelos é normalmente feita utilizando-se dados de ensaios não invasivos e invasivos em modelos sintéticos. Em ensaios in vitro os sensores convencionais têm um princípio de funcionamento eléctrico e apresentam por vezes dimensões inadequadas. Existem situações exploradas no presente trabalho, tais como sensoriamento de superfícies irregulares e junções ou ainda análises de deformações internas, onde é recomendável a utilização de sensores FBG, pois apresentam dimensões reduzidas e flexibilidade o que permite efectuar medidas localizadas. O desenvolvimento de um protocolo de utilização de FBG e a sua aplicação no contexto apresentado demonstrou-se mais adequado, pela precisão e segurança futura oferecidas. Foi desenvolvida uma metodologia experimental para medidas de deformações utilizando FBG ao longo de uma placa de osteossíntese metálica aparafusada a um fémur sintético fracturado. Foi efectuada a monitorização da cura do cimento ósseo utilizado como fixador do prato tibial na artroplastia total do joelho através da medida da sua contracção e temperatura. Foi também desenvolvido um sistema refrigerador com resposta às leituras de temperatura com vista a evitar a necrose do osso. Foram efectuados estudos de deformação nesse cimento após a sua cura, como resultado da aplicação de cargas mecânicas estáticas. Foram efectuados estudos da cura de cimento ósseo aplicado a próteses de anca e também de deformações nestas próteses. Foi ainda efectuado o estudo comparativo de vários implantes dentários através da medida da distribuição de deformações como resposta a excitações mecânicas impulsivas. Para a desmodulação das FBG foram inicialmente utilizados sistemas comerciais. Entretanto algumas aplicações não puderam ser implementadas com estes sistemas comerciais devido à baixa reflectividade das FBG utilizadas, mas fundamentalmente devido à necessidade de executar testes com uma taxa de aquisição maior do que os 5 Hz disponíveis (cerca de 15 kHz). Por estes motivos foi desenvolvido um sistema optoelectrónico completo de desmodulação de FBG baseado num filtro sintonizável e que tem como característica principal a alta taxa de aquisição (até 1,2 MHz) mas também se destaca pela facilidade na reconfiguração dos parâmetros de leitura, pela apresentação duma interface de utilizador amigável e pela capacidade de operar com até 5 FBG na mesma fibra óptica.

Embryonic and larval ecology and biochemistry of offshore decapods

Understanding the biology of offshore species is hardened by the difficulties of sampling in the deep-sea environment. Additionally, due to the vastness of the open ocean, knowledge of early life histories of pelagic larvae is still relatively scarce. In decapod species with bentho-pelagic lifestyle, the transition from life in the seafloor to the water column not only is associated with drastic morphological metamorphosis, but also with changes in behavior and feeding ecology. The purpose of the present thesis was to investigate physiological, biochemical and behavioral adaptation occurring during early development of such species. The Norway lobster, Nephrops norvegicus, and the crab Monodaeus couchi were used as a model as these two species are encountered off the NE Atlantic shelf at depth greater than 300 m. Chapter 1 introduces the challenges faced by both adult and larvae inhabiting such remote habitats, including the effect of food availability on development and oceanographic processes on dispersal and recruitment. The thesis follows early life histories, starting with within-brood variability in the fatty acid (FA) profile displayed by developing N. norvegicus embryos. There were no differences in the FA composition of embryos sampled from both sides of the brooding chamber in most females. However, all females exhibited significant differences in the FA profiles of embryos sampled from different pleopods. Potential causes for the variations recorded may be differential female investment during oocyte production or shifts in FA catabolism during the incubation period promoted by embryo’s location within the brooding chamber. Next, feeding rates and digestive enzymes activity of the early stage larvae was investigated in N. norvegicus. Both stages were able to maximize food intake when larvae were scarce and showed increased feeding rate following periods of starvation. Amylase activity indicated that carbohydrates are not the primary energy reserve and that feeding may be required soon after hatching to trigger amylase activity. Protease activity indicated that protein reserves are catabolized under starvation. These results indicate that larvae may maximize prey ingestion in the presence of plankton patches with higher food abundance and minimize the deleterious effects induced by previous periods of intermittent starvation or unsuitable prey densities/types. Additionally, changes in enzymatic activity may allow newly hatched N. norvegicus larvae to metabolize protein reserves to overcome short-term starvation. Vertical migration behavior and the influence of oceanographic properties were studied next. All zoeal stages of M. couchi displayed reverse diel vertical migration. Abundance of early stages was correlated with chlorophyll a levels. An ontogenic shift in vertical distribution explained the results; earlier zoeal stages remain in the food-rich upper water column while later stages migrate to the bottom for settlement. This vertical migration behavior is likely to affect horizontal distribution of larvae. Indeed, global current patterns will result in low inter-annual variations in decapod larvae recruitment, but short term variations such as upwelling events will cause deviation from the expected dispersal pattern. Throughout development, from the embryo to metamorphosis into benthic juvenile, offshore decapods face many challenges. For the developing individual survivorship will depend heavily on food availability but also on the reserves passed on by the mother. Even though vertical migration behavior can allow the larvae to take advantage of depth varying currents for transport, the effect of general circulation pattern will superimpose local current and influence feeding conditions and affect dispersal and recruitment.