A Hybrid Swarm Based Approach in Addressing University Timetabling Problems

This paper is concerned with the application of an automated hybrid approach in addressing the university timetabling problem. The approach described is based on the nature-inspired artificial bee colony (ABC) algorithm. An ABC algorithm is a biologically-inspired optimization approach, which has been widely implemented in solving a range of optimization problems in recent years such as job shop scheduling and machine timetabling problems. Although the approach has proven to be robust across a range of problems, it is acknowledged within the literature that there currently exist a number of inefficiencies regarding the exploration and exploitation abilities. These inefficiencies can often lead to a slow convergence speed within the search process. Hence, this paper introduces a variant of the algorithm which utilizes a global best model inspired from particle swarm optimization to enhance the global exploration ability while hybridizing with the great deluge (GD) algorithm in order to improve the local exploitation ability. Using this approach, an effective balance between exploration and exploitation is attained. In addition, a traditional local search approach is incorporated within the GD algorithm with the aim of further enhancing the performance of the overall hybrid method. To evaluate the performance of the proposed approach, two diverse university timetabling datasets are investigated, i.e., Carter's examination timetabling and Socha course timetabling datasets. It should be noted that both problems have differing complexity and different solution landscapes. Experimental results demonstrate that the proposed method is capable of producing high quality solutions across both these benchmark problems, showing a good degree of generality in the approach. Moreover, the proposed method produces best results on some instances as compared with other approaches presented in the literature.

Processo de cálculo automático de colunas mistas

Este relatório é elaborado no âmbito do estágio curricular no regime de parceria entre o Instituto Superior de Engenharia do Porto e o gabinete de projetos SE2P – Sociedade de Engenharia, Projetos e Planeamento, Lda. O tema do estágio é o estudo de colunas mistas aço-betão na vertente de projeto, mais concretamente através da criação dum processo de cálculo automático integrado com o modelo estrutural global. Aplicam-se os métodos de dimensionamento dos eurocódigos estruturais na avaliação da segurança dos estados limites. Avaliam-se as propriedades geométricas, determina-se a curva de interação, calcula-se a posição do eixo neutro donde se inferem os estados de tensão em flexão desviada ao nível da secção, necessários às verificações de segurança. O processo de cálculo tem por base uma folha de cálculo desenvolvida em Excel que, com ajuda do Visual Basic for Applications, comunica com o software de modelação e cálculo estrutural Autodesk Robot Professional Analisys. Cria-se fluxos de informação que possibilitam o cálculo iterativo de colunas mistas permitindo ajustar e por esta via otimizar as soluções de dimensionamento. A exportação de relatórios de cálculo detalhados para cada coluna constitui uma mais-valia para gabinetes de engenharia quando comparada com outras soluções de dimensionamento existentes no mercado. A aplicação do processo de cálculo desenvolvido a múltiplas combinações de soluções de colunas mistas permitiu analisar comparativamente os resultados obtidos e a criação de tabelas e ábacos que poderão ser úteis em cenários de pré-dimensionamento e de anteprojeto.

Avaliação económica da vacinação contra o vírus do papiloma humano 16 e 18 (CERVARIX®) : análise de custo-utilidade e custo-efetividade em mulheres dos 26 aos 55 anos em Portugal

RESUMO - O cancro do colo útero representa um importante problema de saúde pública em Portugal: é o terceiro cancro mais frequente nas mulheres entre os 15 e os 44 anos, originando a morte de 346 mulheres anualmente. Contudo, esta patologia é altamente evitável, nomeadamente, através da imunização contra a infeção HPV, que é a causa necessária para o desenvolvimento do cancro. A elevada prevalência da infeção em mulheres mais velhas sugere que a vacinação poderá ser uma estratégia custo-efetiva mesmo numa faixa etária superior. Para que seja racionalmente ponderada a comparticipação da vacina nestas mulheres é necessária a realização de um estudo fármaco-económico que comprove o custo-efetividade desta intervenção, já que o seu financiamento atual prevê apenas as mulheres não abrangidas pelo Programa Nacional de Vacinação, dos 18 aos 25 anos. Os objetivos do trabalho são realizar uma revisão da literatura sobre estudos de avaliação económica relativos à prevenção do CCU e avaliar a relação de custo-efetividade de vacinar mulheres contra o HPV entre os 26 e os 55 anos em comparação com a prática clínica corrente, em Portugal. É utilizado o Modelo Global Cervarix® e realiza-se uma análise de custo-utilidade e de custo-efetividade. Os resultados demonstraram que a vacinação em mulheres dos 26 aos 45 anos poderá ser uma opção custo-efetiva, permitindo um aumento de anos de vida, uma diminuição dos casos e mortes por CCU e um incremento de QALYs. O RCEI variou entre 7.914€/QALY e 29.049€/QALY com a vacinação aos 26 e aos 45 anos, respetivamente, para a alternativa de vacinação mais rastreio versus a situação atual de rastreio organizado e oportunístico, em Portugal.

Dynamo Magnétohydrodynamique en champ moyen

De nos jours, il est bien accepté que le cycle magnétique de 11 ans du Soleil est l'oeuvre d'une dynamo interne présente dans la zone convective. Bien qu'avec la puissance de calculs des ordinateurs actuels il soit possible, à l'aide de véritables simulations magnétohydrodynamiques, de résoudre le champ magnétique et la vitessse dans toutes les directions spatiales, il n'en reste pas moins que pour étudier l'évolution temporelle et spatiale de la dynamo solaire à grande échelle, il reste avantageux de travailler avec des modèles plus simples. Ainsi, nous avons utilisé un modèle simplifié de la dynamo solaire, nommé modèle de champ moyen, pour mieux comprendre les mécanismes importants à l'origine et au maintien de la dynamo solaire. L'insertion d'un tenseur-alpha complet dans un modèle dynamo de champ moyen, provenant d'un modèle global-MHD [Ghizaru et al., 2010] de la convection solaire, nous a permis d'approfondir le rôle que peut jouer la force électromotrice dans les cycles magnétiques produits par ce modèle global. De cette façon, nous avons pu reproduire certaines caractéristiques observées dans les cycles magnétiques provenant de la simulation de Ghizaru et al., 2010. Tout d'abord, le champ magnétique produit par le modèle de champ moyen présente deux modes dynamo distincts. Ces modes, de périodes similaires, sont présents et localisés sensiblement aux mêmes rayons et latitudes que ceux produits par le modèle global. Le fait que l'on puisse reproduire ces deux modes dynamo est dû à la complexité spatiale du tenseur-alpha. Par contre, le rapport entre les périodes des deux modes présents dans le modèle de champ moyen diffère significativement de celui trouvé dans le modèle global. Par ailleurs, on perd l'accumulation d'un fort champ magnétique sous la zone convective dans un modèle où la rotation différentielle n'est plus présente. Ceci suggère que la présence de rotation différentielle joue un rôle non négligeable dans l'accumulation du champ magnétique à cet endroit. Par ailleurs, le champ magnétique produit dans un modèle de champ moyen incluant un tenseur-alpha sans pompage turbulent global est très différent de celui produit par le tenseur original. Le pompage turbulent joue donc un rôle fondamental au sein de la distribution spatiale du champ magnétique. Il est important de souligner que les modèles dépourvus d'une rotation différentielle, utilisant le tenseur-alpha original ou n'utilisant pas de pompage turbulent, parviennent tous deux à produire une dynamo oscillatoire. Produire une telle dynamo à l'aide d'un modèle de ce type n'est pas évident, a priori. Finalement, l'intensité ainsi que le type de profil de circulation méridienne utilisés sont des facteurs affectant significativement la distribution spatiale de la dynamo produite.

Modeling the interplay between land use, bioenergy production and climate change

Land use has become a force of global importance, considering that 34% of the Earth’s ice-free surface was covered by croplands or pastures in 2000. The expected increase in global human population together with eminent climate change and associated search for energy sources other than fossil fuels can, through land-use and land-cover changes (LUCC), increase the pressure on nature’s resources, further degrade ecosystem services, and disrupt other planetary systems of key importance to humanity. This thesis presents four modeling studies on the interplay between LUCC, increased production of biofuels and climate change in four selected world regions. In the first study case two new crop types (sugarcane and jatropha) are parameterized in the LPJ for managed Lands dynamic global vegetation model for calculation of their potential productivity. Country-wide spatial variation in the yields of sugarcane and jatropha incurs into substantially different land requirements to meet the biofuel production targets for 2015 in Brazil and India, depending on the location of plantations. Particularly the average land requirements for jatropha in India are considerably higher than previously estimated. These findings indicate that crop zoning is important to avoid excessive LUCC. In the second study case the LandSHIFT model of land-use and land-cover changes is combined with life cycle assessments to investigate the occurrence and extent of biofuel-driven indirect land-use changes (ILUC) in Brazil by 2020. The results show that Brazilian biofuels can indeed cause considerable ILUC, especially by pushing the rangeland frontier into the Amazonian forests. The carbon debt caused by such ILUC would result in no carbon savings (from using plant-based ethanol and biodiesel instead of fossil fuels) before 44 years for sugarcane ethanol and 246 years for soybean biodiesel. The intensification of livestock grazing could avoid such ILUC. We argue that such an intensification of livestock should be supported by the Brazilian biofuel sector, based on the sector’s own interest in minimizing carbon emissions. In the third study there is the development of a new method for crop allocation in LandSHIFT, as influenced by the occurrence and capacity of specific infrastructure units. The method is exemplarily applied in a first assessment of the potential availability of land for biogas production in Germany. The results indicate that Germany has enough land to fulfill virtually all (90 to 98%) its current biogas plant capacity with only cultivated feedstocks. Biogas plants located in South and Southwestern (North and Northeastern) Germany might face more (less) difficulties to fulfill their capacities with cultivated feedstocks, considering that feedstock transport distance to plants is a crucial issue for biogas production. In the fourth study an adapted version of LandSHIFT is used to assess the impacts of contrasting scenarios of climate change and conservation targets on land use in the Brazilian Amazon. Model results show that severe climate change in some regions by 2050 can shift the deforestation frontier to areas that would experience low levels of human intervention under mild climate change (such as the western Amazon forests or parts of the Cerrado savannas). Halting deforestation of the Amazon and of the Brazilian Cerrado would require either a reduction in the production of meat or an intensification of livestock grazing in the region. Such findings point out the need for an integrated/multicisciplinary plan for adaptation to climate change in the Amazon. The overall conclusions of this thesis are that (i) biofuels must be analyzed and planned carefully in order to effectively reduce carbon emissions; (ii) climate change can have considerable impacts on the location and extent of LUCC; and (iii) intensification of grazing livestock represents a promising venue for minimizing the impacts of future land-use and land-cover changes in Brazil.

Verfahren zur echtzeitfähigen Simulation dynamischer Emissionsverläufe ausgewählter Schadstoffe von Ottomotoren

Im Rahmen dieser Arbeit werden Modellbildungsverfahren zur echtzeitfähigen Simulation wichtiger Schadstoffkomponenten im Abgasstrom von Verbrennungsmotoren vorgestellt. Es wird ein ganzheitlicher Entwicklungsablauf dargestellt, dessen einzelne Schritte, beginnend bei der Ver-suchsplanung über die Erstellung einer geeigneten Modellstruktur bis hin zur Modellvalidierung, detailliert beschrieben werden. Diese Methoden werden zur Nachbildung der dynamischen Emissi-onsverläufe relevanter Schadstoffe des Ottomotors angewendet. Die abgeleiteten Emissionsmodelle dienen zusammen mit einer Gesamtmotorsimulation zur Optimierung von Betriebstrategien in Hybridfahrzeugen. Im ersten Abschnitt der Arbeit wird eine systematische Vorgehensweise zur Planung und Erstellung von komplexen, dynamischen und echtzeitfähigen Modellstrukturen aufgezeigt. Es beginnt mit einer physikalisch motivierten Strukturierung, die eine geeignete Unterteilung eines Prozessmodells in einzelne überschaubare Elemente vorsieht. Diese Teilmodelle werden dann, jeweils ausgehend von einem möglichst einfachen nominalen Modellkern, schrittweise erweitert und ermöglichen zum Abschluss eine robuste Nachbildung auch komplexen, dynamischen Verhaltens bei hinreichender Genauigkeit. Da einige Teilmodelle als neuronale Netze realisiert werden, wurde eigens ein Verfah-ren zur sogenannten diskreten evidenten Interpolation (DEI) entwickelt, das beim Training einge-setzt, und bei minimaler Messdatenanzahl ein plausibles, also evidentes Verhalten experimenteller Modelle sicherstellen kann. Zum Abgleich der einzelnen Teilmodelle wurden statistische Versuchs-pläne erstellt, die sowohl mit klassischen DoE-Methoden als auch mittels einer iterativen Versuchs-planung (iDoE ) generiert wurden. Im zweiten Teil der Arbeit werden, nach Ermittlung der wichtigsten Einflussparameter, die Model-strukturen zur Nachbildung dynamischer Emissionsverläufe ausgewählter Abgaskomponenten vor-gestellt, wie unverbrannte Kohlenwasserstoffe (HC), Stickstoffmonoxid (NO) sowie Kohlenmono-xid (CO). Die vorgestellten Simulationsmodelle bilden die Schadstoffkonzentrationen eines Ver-brennungsmotors im Kaltstart sowie in der anschließenden Warmlaufphase in Echtzeit nach. Im Vergleich zur obligatorischen Nachbildung des stationären Verhaltens wird hier auch das dynami-sche Verhalten des Verbrennungsmotors in transienten Betriebsphasen ausreichend korrekt darge-stellt. Eine konsequente Anwendung der im ersten Teil der Arbeit vorgestellten Methodik erlaubt, trotz einer Vielzahl von Prozesseinflussgrößen, auch hier eine hohe Simulationsqualität und Ro-bustheit. Die Modelle der Schadstoffemissionen, eingebettet in das dynamische Gesamtmodell eines Ver-brennungsmotors, werden zur Ableitung einer optimalen Betriebsstrategie im Hybridfahrzeug ein-gesetzt. Zur Lösung solcher Optimierungsaufgaben bieten sich modellbasierte Verfahren in beson-derer Weise an, wobei insbesondere unter Verwendung dynamischer als auch kaltstartfähiger Mo-delle und der damit verbundenen Realitätsnähe eine hohe Ausgabequalität erreicht werden kann.

River ecosystems at risk

Worldwide water managers are increasingly challenged to allocate sufficient and affordable water supplies to different water use sectors without further degrading river ecosystems and their valuable services to mankind. Since 1950 human population almost tripled, water abstractions increased by a factor of four, and the number of large dam constructions is about eight times higher today. From a hydrological perspective, the alteration of river flows (temporally and spatially) is one of the main consequences of global change and further impairments can be expected given growing population pressure and projected climate change. Implications have been addressed in numerous hydrological studies, but with a clear focus on human water demands. Ecological water requirements have often been neglected or addressed in a very simplistic manner, particularly from the large-scale perspective. With his PhD thesis, Christof Schneider took up the challenge to assess direct (dam operation and water abstraction) and indirect (climate change) impacts of human activities on river flow regimes and evaluate the consequences for river ecosystems by using a modeling approach. The global hydrology model WaterGAP3 (developed at CESR) was applied and further developed within this thesis to carry out several model experiments and assess anthropogenic river flow regime modifications and their effects on river ecosystems. To address the complexity of ecological water requirements the assessment is based on three main ideas: (i) the natural flow paradigm, (ii) the perception that different flows have different ecological functions, and (iii) the flood pulse concept. The thesis shows that WaterGAP3 performs well in representing ecologically relevant flow characteristics on a daily time step, and therefore justifies its application within this research field. For the first time a methodology was established to estimate bankfull flow on a 5 by 5 arc minute grid cell raster globally, which is a key parameter in eFlow assessments as it marks the point where rivers hydraulically connect to adjacent floodplains. Management of dams and water consumption pose a risk to floodplains and riparian wetlands as flood volumes are significantly reduced. The thesis highlights that almost one-third of 93 selected Ramsar sites are seriously affected by modified inundation patterns today, and in the future, inundation patterns are very likely to be further impaired as a result of new major dam initiatives and climate change. Global warming has been identified as a major threat to river flow regimes as rising temperatures, declining snow cover, changing precipitation patterns and increasing climate variability are expected to seriously modify river flow regimes in the future. Flow regimes in all climate zones will be affected, in particular the polar zone (Northern Scandinavia) with higher river flows during the year and higher flood peaks in spring. On the other side, river flows in the Mediterranean are likely to be even more intermittent in the future because of strong reductions in mean summer precipitation as well as a decrease in winter precipitation, leading to an increasing number of zero flow events creating isolated pools along the river and transitions from lotic to lentic waters. As a result, strong impacts on river ecosystem integrity can be expected. Already today, large amounts of water are withdrawn in this region for agricultural irrigation and climate change is likely to exacerbate the current situation of water shortages.

Oportunidades de negocios bilaterales entre Colombia e Indonesia

El trabajo que se presenta a continuación muestra a grandes rasgos cómo se han desarrollado las relaciones comerciales entre Colombia e Indonesia desde sus inicios hasta la actualidad y qué nuevas oportunidades se pueden presentar para potencializar en mayor medida el comercio bilateral entre las dos naciones. De acuerdo a lo anterior, el desarrollo del trabajo parte en primer lugar de un análisis económico, social y del comercio tanto de Colombia como de Indonesia, en donde se estudian diferentes indicadores que muestran en qué nivel se encuentra cada país en estos aspectos y como es su comercio desde y hacia el mundo. En segundo lugar, se presenta un enfoque en el desarrollo del comercio bilateral, en donde se comienza a evidenciar las posibles nuevas oportunidades de intercambio entre los dos países, posteriormente se realiza un estudio del macro ambiente de cada país, el cual permite hacer un análisis de riesgo país, que arroja temas a tener en cuenta al momento de establecer una relación comercial y económica que permita evidenciar el riesgo que representa una nación para la otra. Por último, se realiza un DOFA con su matriz, de la cual surgen una serie de estrategias para aprovechar las nuevas oportunidades de intercambio comercial que se pudieron evidenciar durante todo el trabajo.

Paisatge i models urbans contemporanis. Les comarques gironines (1979-2006): del "desarrollismo" a la globalització

This doctoral thesis offers a quantitative and qualitative analysis of the changes in the urban shape and landscape of the Girona Counties between 1979 and 2006. The theoretical part of the research lies within the framework of the dispersed city phenomenon, and is based on the hypothesis of convergence towards a global urban model. The empirical part demonstrates this proposition with a study of 522 zone development plans in the Girona Counties. The results point to the consolidation of the dispersed city phenomenon, as shown by the sudden increase in built-up space, the spread of urban development throughout the territory, and the emergence of a new, increasingly generic landscape comprising three major morphological types: urban extensions, low density residential estates and industrial zones. This reveals shortcomings of planning for urban growth, weakening of the city as a public project, and a certain degradation of the Mediterranean city model.

The response of tropospheric circulation to perturbations in lower-stratospheric temperature

A multiple regression analysis of the NCEP-NCAR reanalysis dataset shows a response to increased solar activity of a weakening and poleward shift of the subtropical jets. This signal is separable from other influences, such as those of El Nino-Southern Oscillation (ENSO) and the North Atlantic Oscillation (NAO), and is very similar to that seen in previous studies using global circulation models (GCMs) of the effects of an increase in solar spectral irradiance. The response to increased stratospheric (volcanic) aerosol is found in the data to be a weakening and equatorward shift of the jets. The GCM studies of the solar influence also showed an impact on tropospheric mean meridional circulation with a weakening and expansion of the tropical Hadley cells and a poleward shift of the Ferrel cells. To understand the mechanisms whereby the changes in solar irradiance affect tropospheric winds and circulation, experiments have been carried out with a simplified global circulation model. The results show that generic heating of the lower stratosphere tends to weaken the subtropical jets and the tropospheric mean meridional circulations. The positions of the jets, and the extent of the Hadley cells, respond to the distribution of the stratospheric heating, with low-latitude heating forcing them to move poleward, and high-latitude or latitudinally uniform heating forcing them equatorward. The patterns of response are similar to those that are found to be a result of the solar or volcanic influences, respectively, in the data analysis. This demonstrates that perturbations to the heat balance of the lower stratosphere, such as those brought about by solar or volcanic activity, can produce changes in the mean tropospheric circulation, even without any direct forcing below the tropopause.

Variability of the Atlantic thermohaline circulation described by three-dimensional empirical orthogonal functions

We describe the use of bivariate 3d empirical orthogonal functions (EOFs) in characterising low frequency variability of the Atlantic thermohaline circulation (THC) in the Hadley Centre global climate model, HadCM3. We find that the leading two modes are well correlated with an index of the meridional overturning circulation (MOC) on decadal timescales, with the leading mode alone accounting for 54% of the decadal variance. Episodes of coherent oscillations in the sub-space of the leading EOFs are identified; these episodes are of great interest for the predictability of the THC, and could indicate the existence of different regimes of natural variability. The mechanism identified for the multi-decadal variability is an internal ocean mode, dominated by changes in convection in the Nordic Seas, which lead the changes in the MOC by a few years. Variations in salinity transports from the Arctic and from the North Atlantic are the main feedbacks which control the oscillation. This mode has a weak feedback onto the atmosphere and hence a surface climatic influence. Interestingly, some of these climate impacts lead the changes in the overturning. There are also similarities to observed multi-decadal climate variability.

Potential predictability of rapid changes in the Atlantic meridional overturning circulation

We explore the potential predictability of rapid changes in the Atlantic meridional overturning circulation (MOC) using a coupled global climate model (HadCM3). Rapid changes in the temperature and salinity of surface water in the Nordic Seas, and the flow of dense water through Denmark Strait, are found to be precursors to rapid changes in the model's MOC, with a lead time of around 10 years. The mechanism proposed to explain this potential predictability relies on the development of density anomalies in the Nordic Seas which propagate through Denmark Strait and along the deep western boundary current, affecting the overturning. These rapid changes in the MOC have significant, and widespread, climate impacts which are potentially predictable a few years ahead. Whilst the flow through Denmark Strait is too strong in HadCM3, the presence of such potential predictability motivates the monitoring of water properties in the Nordic Seas and Denmark Strait.

Decadal predictability of the Atlantic Ocean in a coupled GCM: forecast skill and optimal perturbations using Linear Inverse Modelling

The decadal predictability of three-dimensional Atlantic Ocean anomalies is examined in a coupled global climate model (HadCM3) using a Linear Inverse Modelling (LIM) approach. It is found that the evolution of temperature and salinity in the Atlantic, and the strength of the meridional overturning circulation (MOC), can be effectively described by a linear dynamical system forced by white noise. The forecasts produced using this linear model are more skillful than other reference forecasts for several decades. Furthermore, significant non-normal amplification is found under several different norms. The regions from which this growth occurs are found to be fairly shallow and located in the far North Atlantic. Initially, anomalies in the Nordic Seas impact the MOC, and the anomalies then grow to fill the entire Atlantic basin, especially at depth, over one to three decades. It is found that the structure of the optimal initial condition for amplification is sensitive to the norm employed, but the initial growth seems to be dominated by MOC-related basin scale changes, irrespective of the choice of norm. The consistent identification of the far North Atlantic as the most sensitive region for small perturbations suggests that additional observations in this region would be optimal for constraining decadal climate predictions.

Will Extratropical Storms Intensify in a Warmer Climate?

Extratropical cyclones and how they may change in a warmer climate have been investigated in detail with a high-resolution version of the ECHAM5 global climate model. A spectral resolution of T213 (63 km) is used for two 32-yr periods at the end of the twentieth and twenty-first centuries and integrated for the Intergovernmental Panel on Climate Change (IPCC) A1B scenario. Extremes of pressure, vorticity, wind, and precipitation associated with the cyclones are investigated and compared with a lower-resolution simulation. Comparison with observations of extreme wind speeds indicates that the model reproduces realistic values. This study also investigates the ability of the model to simulate extratropical cyclones by computing composites of intense storms and contrasting them with the same composites from the 40-yr ECMWF Re-Analysis (ERA-40). Composites of the time evolution of intense cyclones are reproduced with great fidelity; in particular the evolution of central surface pressure is almost exactly replicated, but vorticity, maximum wind speed, and precipitation are higher in the model. Spatial composites also show that the distributions of pressure, winds, and precipitation at different stages of the cyclone life cycle compare well with those from ERA-40, as does the vertical structure. For the twenty-first century, changes in the distribution of storms are very similar to those of previous study. There is a small reduction in the number of cyclones but no significant changes in the extremes of wind and vorticity in both hemispheres. There are larger regional changes in agreement with previous studies. The largest changes are in the total precipitation, where a significant increase is seen. Cumulative precipitation along the tracks of the cyclones increases by some 11% per track, or about twice the increase in global precipitation, while the extreme precipitation is close to the globally averaged increase in column water vapor (some 27%). Regionally, changes in extreme precipitation are even higher because of changes in the storm tracks.

Annular variability and eddy-zonal flow interactions in a simplified atmospheric GCM. Part 1: Characterization of high and low frequency behaviour

Experiments have been performed using a simplified, Newtonian forced, global circulation model to investigate how variability of the tropospheric jet can be characterized by examining the combined fluctuations of the two leading modes of annular variability. Eddy forcing of this variability is analyzed in the phase space of the leading modes using the vertically integrated momentum budget. The nature of the annular variability and eddy forcing depends on the time scale. At low frequencies the zonal flow and baroclinic eddies are in quasi equilibrium and anomalies propagate poleward. The eddies are shown primarily to reinforce the anomalous state and are closely balanced by the linear damping, leaving slow evolution as a residual. At high frequencies the flow is strongly evolving and anomalies are initiated on the poleward side of the tropospheric jet and propagate equatorward. The eddies are shown to drive this evolution strongly: eddy location and amplitude reflect the past baroclinicity, while eddy feedback on the zonal flow may be interpreted in terms of wave breaking associated with baroclinic life cycles in lateral shear.