Matching stages of heavy-ion collision models

Heavy-ion reactions and other collective dynamical processes are frequently described by different theoretical approaches for the different stages of the process, like initial equilibration stage, intermediate locally equilibrated fluid dynamical stage, and final freeze-out stage. For the last stage, the best known is the Cooper-Frye description used to generate the phase space distribution of emitted, noninteracting particles from a fluid dynamical expansion or explosion, assuming a final ideal gas distribution, or (less frequently) an out-of-equilibrium distribution. In this work we do not want to replace the Cooper-Frye description, but rather clarify the ways of using it and how to choose the parameters of the distribution and, eventually, how to choose the form of the phase space distribution used in the Cooper-Frye formula. Moreover, the Cooper-Frye formula is used in connection with the freeze-out problem, while the discussion of transition between different stages of the collision is applicable to other transitions also. More recently, hadronization and molecular dynamics models have been matched to the end of a fluid dynamical stage to describe hadronization and freeze-out. The stages of the model description can be matched to each other on space-time hypersurfaces (just like through the frequently used freeze-out hypersurface). This work presents a generalized description of how to match the stages of the description of a reaction to each other, extending the methodology used at freeze-out, in simple covariant form which is easily applicable in its simplest version for most applications.

Edge matching puzzles as hard SAT/CSP benchmarks

Recently, edge matching puzzles, an NP-complete problem, have rececived, thanks to money-prized contests, considerable attention from wide audiences. We consider these competitions not only a challenge for SAT/CSP solving techniques but also as an opportunity to showcase the advances in the SAT/CSP community to a general audience. This paper studies the NP-complete problem of edge matching puzzles focusing on providing generation models of problem instances of variable hardness and on its resolution through the application of SAT and CSP techniques. From the generation side, we also identify the phase transition phenomena for each model. As solving methods, we employ both; SAT solvers through the translation to a SAT formula, and two ad-hoc CSP solvers we have developed, with different levels of consistency, employing several generic and specialized heuristics. Finally, we conducted an extensive experimental investigation to identify the hardest generation models and the best performing solving techniques.

Inférence de source de traces d'essence retrouvées dans des débris d'incendie: Evaluation de la contribution de la chromatographie en phase gazeuse couplée à la spectrométrie de masse à rapport isotopique

Lorsque de l'essence est employée pour allumer et/ou propager un incendie, l'inférence de la source de l'essence peut permettre d'établir un lien entre le sinistre et une source potentielle. Cette inférence de la source constitue une alternative intéressante pour fournir des éléments de preuve dans ce type d'événements où les preuves matérielles laissées par l'auteur sont rares. Le but principal de cette recherche était le développement d'une méthode d'analyse de spécimens d'essence par GC-IRMS, méthode pas routinière et peu étudiée en science forensique, puis l'évaluation de son potentiel à inférer la source de traces d'essence en comparaison aux performances de la GC-MS. Un appareillage permettant d'analyser simultanément les échantillons par MS et par IRMS a été utilisé dans cette recherche. Une méthode d'analyse a été développée, optimisée et validée pour cet appareillage. Par la suite, des prélèvements d'essence provenant d'un échantillonnage conséquent et représentatif du marché de la région lausannoise ont été analysés. Finalement, les données obtenues ont été traitées et interprétées à l'aide de méthodes chimiométriques. Les analyses effectuées ont permis de montrer que la méthodologie mise en place, aussi bien pour la composante MS que pour l'IRMS, permet de différencier des échantillons d'essence non altérée provenant de différentes stations-service. Il a également pu être démontré qu'à chaque nouveau remplissage des cuves d'une station-service, la composition de l'essence distribuée par cette station est quasi unique. La GC-MS permet une meilleure différenciation d'échantillons prélevés dans différentes stations, alors que la GC-IRMS est plus performante lorsqu'il s'agit de comparer des échantillons collectés après chacun des remplissages d'une cuve. Ainsi, ces résultats indiquent que les deux composantes de la méthode peuvent être complémentaires pour l'analyse d'échantillons d'essence non altérée. Les résultats obtenus ont également permis de montrer que l'évaporation des échantillons d'essence ne compromet pas la possibilité de grouper des échantillons de même source par GC-MS. Il est toutefois nécessaire d'effectuer une sélection des variables afin d'éliminer celles qui sont influencées par le phénomène d'évaporation. Par contre, les analyses effectuées ont montré que l'évaporation des échantillons d'essence a une forte influence sur la composition isotopique des échantillons. Cette influence est telle qu'il n'est pas possible, même en effectuant une sélection des variables, de grouper correctement des échantillons évaporés par GC-IRMS. Par conséquent, seule la composante MS de la méthodologie mise en place permet d'inférer la source d'échantillons d'essence évaporée. _________________________________________________________________________________________________ When gasoline is used to start and / or propagate an arson, source inference of gasoline can allow to establish a link between the fire and a potential source. This source inference is an interesting alternative to provide evidence in this type of events where physical evidence left by the author are rare. The main purpose of this research was to develop a GC-IRMS method for the analysis of gasoline samples, a non-routine method and little investigated in forensic science, and to evaluate its potential to infer the source of gasoline traces compared to the GC-MS performances. An instrument allowing to analyze simultaneously samples by MS and IRMS was used in this research. An analytical method was developed, optimized and validated for this instrument. Thereafter, gasoline samples from a large sampling and representative of the Lausanne area market were analyzed. Finally, the obtained data were processed and interpreted using chemometric methods. The analyses have shown that the methodology, both for MS and for IRMS, allow to differentiate unweathered gasoline samples from different service stations. It has also been demonstrated that each new filling of the tanks of a station generates an almost unique composition of gasoline. GC-MS achieves a better differentiation of samples coming from different stations, while GC-IRMS is more efficient to distinguish samples collected after each filling of a tank. Thus, these results indicate that the two components of the method can be complementary to the analysis of unweathered gasoline samples. The results have also shown that the evaporation of gasoline samples does not compromise the possibility to group samples coming from the same source by GC-MS. It is however necessary to make a selection of variables in order to eliminate those which are influenced by the evaporation. On the other hand, the carried out analyses have shown that the evaporation of gasoline samples has such a strong influence on the isotopic composition of the samples that it is not possible, even by performing a selection of variables, to properly group evaporated samples by GC-IRMS. Therefore, only the MS allows to infer the source of evaporated gasoline samples.

Monte Carlo study of the XY-model on quasi-periodic lattices

Monte Carlo Simulations were carried out using a nearest neighbour ferromagnetic XYmodel, on both 2-D and 3-D quasi-periodic lattices. In the case of 2-D, both the unfrustrated and frustrated XV-model were studied. For the unfrustrated 2-D XV-model, we have examined the magnetization, specific heat, linear susceptibility, helicity modulus and the derivative of the helicity modulus with respect to inverse temperature. The behaviour of all these quatities point to a Kosterlitz-Thouless transition occuring in temperature range Te == (1.0 -1.05) JlkB and with critical exponents that are consistent with previous results (obtained for crystalline lattices) . However, in the frustrated case, analysis of the spin glass susceptibility and EdwardsAnderson order parameter, in addition to the magnetization, specific heat and linear susceptibility, support a spin glass transition. In the case where the 'thin' rhombus is fully frustrated, a freezing transition occurs at Tf == 0.137 JlkB , which contradicts previous work suggesting the critical dimension of spin glasses to be de > 2 . In the 3-D systems, examination of the magnetization, specific heat and linear susceptibility reveal a conventional second order phase transition. Through a cumulant analysis and finite size scaling, a critical temperature of Te == (2.292 ± 0.003) JI kB and critical exponents of 0:' == 0.03 ± 0.03, f3 == 0.30 ± 0.01 and I == 1.31 ± 0.02 have been obtained.

Investigation of Frustrated Quasi-One-Dimensional Quantum Spin-Chain Materials

Copper arsenite CuAs2O4 and Copper antimonite CuSb2O4 are S=1/2 (Cu2+ 3d9 electronic configuration) quasi-one-dimensional quantum spin-chain compounds. Both compounds crystallize with tetragonal structures containing edge sharing CuO6 octahedra chains which experience Jahn-Teller distortions. The basal planes of the octahedra link together to form CuO2 ribbon-chains which harbor Cu2+ spin-chains. These compounds are magnetically frustrated with competing nearest-neighbour and next-nearest-neighbour intrachain spin-exchange interactions. Despite the similarities between CuAs2O4 and CuSb2O4, they exhibit very different magnetic properties. In this thesis work, the physical properties of CuAs2O4 and CuSb2O4 are investigated using a variety of experimental techniques which include x-ray diffraction, magnetic susceptibility measurements, heat capacity measurements, Raman spectroscopy, electron paramagnetic resonance, neutron diffraction, and dielectric capacitance measurements. CuAs2O4 exhibits dominant ferromagnetic nearest-neighbour and weaker antiferromagnetic next-nearest-neighbour intrachain spin-exchange interactions. The ratio of the intrachain interactions amounts to Jnn/Jnnn = -4.1. CuAs2O4 was found to order with a ferromagnetic groundstate below TC = 7.4 K. An extensive physical characterization of the magnetic and structural properties of CuAs2O4 was carried out. Under the effect of hydrostatic pressure, CuAs2O4 was found to undergo a structural phase transition at 9 GPa to a new spin-chain structure. The structural phase transition is accompanied by a severe alteration of the magnetic properties. The high-pressure phase exhibits dominant ferromagnetic next-nearest-neighbour spin-exchange interactions and weaker ferromagnetic nearest-neighbour interactions. The ratio of the intrachain interactions in the high-pressure phase was found to be Jnn/Jnnn = 0.3. Structural and magnetic characterizations under hydrostatic pressure are reported and a relationship between the structural and magnetic properties was established. CuSb2O4 orders antiferromagnetically below TN = 1.8 K with an incommensurate helicoidal magnetic structure. CuSb2O4 is characterized by ferromagnetic nearest-neighbour and antiferromagnetic next-nearest-neighbour spin-exchange interactions with Jnn/Jnnn = -1.8. A (H, T) magnetic phase diagram was constructed using low-temperature magnetization and heat capacity measurements. The resulting phase diagram contains multiple phases as a consequence of the strong intrachain magnetic frustration. Indications of ferroelectricity were observed in the incommensurate antiferromagnetic phase.

Caractérisation cytogénétique et moléculaire des translocations chromosomiques dans la phase blastique de la leucémie myéloïde chronique

La leucémie myéloïde chronique (LMC) est un modèle d’évolution tumorale dans les cancers humains. Le processus d’évolution de la LMC de la phase chronique (PC) à la phase blastique (PB) est caractérisé par un arrêt de différenciation et l’acquisition de la capacité d’autorenouvellement incontrôlé d’une cellule souche ou d’un progéniteur hématopoïétique. La LMC en PB est associée à la présence d’anomalies génétiques additionnelles à la fusion BCR-ABL1 qui résulte de la translocation chromosomique t(9;22). Contrairement aux patients en PC, les patients en PB de la LMC n’obtiennent pas une réponse moléculaire complète à long terme avec 1’Imatinib mesylate, un inhibiteur de la tyrosine kinase (ITK) BCR-ABL1. De plus, les ITKs de deuxième et troisième générations sont moins efficaces en PB de la LMC lorsque les cellules leucémiques ont acquis une résistance au traitement indépendante des mutations de BCR-ABL1. Les mécanismes moléculaires des voies de signalisation impliquées dans la progression de la LMC en PB ne sont pas entièrement élucidés. Le but de notre travail est de caractériser de nouvelles anomalies génétiques dans la PB de la LMC. Nous avons identifié en cytogénétique, quatre nouvelles translocations chromosomiques : t(1;21)(p36;q22), t(7;17)(p15;q22), t(8;17)(q11;q22) et t(2;12)(q31;p13) dans les cellules leucémiques de patients en PB de la LMC résistants au traitement. En utilisant des techniques d'hybridation in situ en fluorescence, de RT-PCR et de séquençage, nous avons délimité les régions à investiguer au niveau des points de cassure et identifié un réarrangement de plusieurs gènes codant pour des facteurs de transcription importants lors de l’hématopoïèse tels que RUNX1, ETV6, PRDM16 et HOXA. L’altération de ces gènes pourrait expliquer l’arrêt de différenciation et/ou l’acquisition de la capacité d’autorenouvellement caractéristiques de la LMC en PB. Nous avons identifié les fusions RUNX1-PRDM16, MSI2-HOXA, MSI2-SOX17 et ETV6-HOXD11, respectivement associées aux translocations chromosomiques t(1;21), t(7;17), t(8;17) et t(2;12). Ces fusions génèrent différents transcrits alternatifs qui maintiennent et altèrent le cadre ouvert de lecture. L’analyse des séquences des transcrits chimériques identifiés dans ce projet, incluant RUNX1-PRDM16, MSI2-HOXA9, MSI2-HOXA10, MSI2-HOXA11 et ETV6-HOXD11, nous a permis de prédire les domaines fonctionnels potentiellement présents au niveau des protéines chimériques prédites. Les transcrits de fusion qui respectent le cadre ouvert de lecture peuvent générer des domaines fonctionnels des deux partenaires. C’est le cas des deux transcrits identifiés pour la fusion RUNX1-PRDM16 où le domaine de liaison à l’ADN RHD (Runt homology domain) de RUNX1 est fusionné avec la quasi-totalité des domaines de PRDM16. Les transcrits de fusion qui ne respectent pas le cadre ouvert de lecture donnent des formes tronquées des transcrits RUNX1, MSI2 et ETV6. La juxtaposition des régions promotrices de ces derniers en 5’ de leurs partenaires entraîne l’activation de la forme courte oncogénique de PRDM16 dans la t(1;21) ou de différents gènes HOXA/D dans les t(7;17) et t(2;12), ainsi que l’expression aberrante d’un nouveau transcrit alternatif de SOX17 dans la t(8;17). Notre étude nous a permis d’identifier de nouveaux gènes de fusion et/ou une activation de gènes qui pourraient coopérer avec la fusion BCR-ABL1 dans la progression de la LMC et être impliqués dans la résistance au traitement de la LMC en phase avancée. La caractérisation des événements génétiques associés à la transformation blastique de la LMC est essentielle pour l’investigation des voies moléculaires impliquées dans cette phase de la maladie. Investiguer la résistance au traitement de ces patients pourrait aussi contribuer à identifier de nouvelles cibles thérapeutiques dans cette leucémie.

Étude du contrôle postural quasi-statique et dynamique en position debout des personnes ayant une lésion médullaire incomplète.

L’atteinte sensorimotrice découlant d’une lésion médullaire traumatique affecte la capacité à se tenir debout de façon sécuritaire. Chez les individus ayant une lésion médullaire incomplète, les chutes lors des tâches locomotrices sont fréquentes après la réadaptation, entraînant des blessures qui affectent la participation sociale et la qualité de vie. Une meilleure compréhension du contrôle postural en clinique aiderait à cibler des interventions efficaces à ce niveau. L’objectif général de cette thèse était donc d’étudier le contrôle postural debout lors de tâches variées chez les personnes avec lésion médullaire traumatique en utilisant une approche biomécanique. Les objectifs spécifiques étaient d’étudier divers aspects du contrôle postural en lien avec chacune des tâches et d’identifier les variables prédictives de la stabilité. Vingt-cinq (25) personnes ayant une lésion médullaire traumatique incomplète ont été recrutées. Elles ont été évaluées au laboratoire d’analyse de mouvement lors du maintien de la station debout quasi-statique yeux ouverts et fermés, de l’exécution d’un test des limites de stabilité multidirectionnelles, de la marche naturelle et de l’initiation et l’arrêt de la marche. Des mesures biomécaniques caractérisant le déplacement du centre de pression (COP) ainsi que le modèle du contrôle postural dynamique des forces stabilisantes et déstabilisantes ont été utilisés pour comparer le contrôle postural des patients à celui d’un groupe formé de 33 personnes en santé. Les résultats ont montré une diminution du contrôle postural quasi-statique et dynamique chez les personnes ayant une lésion médullaire comparativement aux personnes en santé. Cette diminution s’accompagnait d’une contribution élevée des informations visuelles à la station debout quasi-statique qui était associée au score du mini BESTest. Le déplacement du COP lors du test des limites de stabilité multidirectionnelle se caractérisait par une difficulté à suivre la direction indiquée. Ce manque de précision causait une augmentation du trajet nécessaire pour atteindre la distance maximale dans chacune des directions. Les résultats au maintien de la station debout quasi-statique et au test des limites de stabilité multidirectionnelle n’étaient pas corrélés. La phase unipodale de la marche des personnes ayant une lésion médullaire différait de celle des personnes en santé par une force stabilisante maximale moindre et une force déstabilisante plus grande alors que l’arrêt de la marche se révélait plus instable que l’initiation de la marche chez les personnes ayant une lésion médullaire. Pour l’ensemble des tâches de marche, la vitesse du COM et la distance entre le COP et la base de support étaient les facteurs explicatifs principaux des forces stabilisantes et déstabilisantes. En somme, les résultats confirment l’impression clinique d’une atteinte générale du contrôle postural debout des personnes ayant une lésion médullaire; ils en précisent les caractéristiques et ciblent les paramètres à considérer dans la rééducation.

Immersive Projection for biochemical pattern matching

In all biological processes, protein molecules and other small molecules interact to function and form transient macromolecular complexes. This interaction of two or more molecules can be described by a docking event. Docking is an important phase for structure-based drug design strategies, as it can be used as a method to simulate protein-ligand interactions. Various docking programs exist that allow automated docking, but most of them have limited visualization and user interaction. It would be advantageous if scientists could visualize the molecules participating in the docking process, manipulate their structures and manually dock them before submitting the new conformations to an automated docking process in an immersive environment, which can help stimulate the design/docking process. This also could greatly reduce docking time and resources. To achieve this, we propose a new virtual modelling/docking program, whereby the advantages of virtual modelling programs and the efficiency of the algorithms in existing docking programs will be merged.

The quasi-nondispersive regimes of long extratropical baroclinic rossby waves over (slowly varying) topography

Actual energy paths of long, extratropical baroclinic Rossby waves in the ocean are difficult to describe simply because they depend on the meridional-wavenumber-to-zonal-wavenumber ratio tau, a quantity that is difficult to estimate both observationally and theoretically. This paper shows, however, that this dependence is actually weak over any interval in which the zonal phase speed varies approximately linearly with tau, in which case the propagation becomes quasi-nondispersive (QND) and describable at leading order in terms of environmental conditions (i.e., topography and stratification) alone. As an example, the purely topographic case is shown to possess three main kinds of QND ray paths. The first is a topographic regime in which the rays follow approximately the contours f/h(alpha c) = a constant (alpha(c) is a near constant fixed by the strength of the stratification, f is the Coriolis parameter, and h is the ocean depth). The second and third are, respectively, "fast" and "slow" westward regimes little affected by topography and associated with the first and second bottom-pressure-compensated normal modes studied in previous work by Tailleux and McWilliams. Idealized examples show that actual rays can often be reproduced with reasonable accuracy by replacing the actual dispersion relation by its QND approximation. The topographic regime provides an upper bound ( in general a large overestimate) of the maximum latitudinal excursions of actual rays. The method presented in this paper is interesting for enabling an optimal classification of purely azimuthally dispersive wave systems into simpler idealized QND wave regimes, which helps to rationalize previous empirical findings that the ray paths of long Rossby waves in the presence of mean flow and topography often seem to be independent of the wavenumber orientation. Two important side results are to establish that the baroclinic string function regime of Tyler and K se is only valid over a tiny range of the topographic parameter and that long baroclinic Rossby waves propagating over topography do not obey any two-dimensional potential vorticity conservation principle. Given the importance of the latter principle in geophysical fluid dynamics, the lack of it in this case makes the concept of the QND regimes all the more important, for they are probably the only alternative to provide a simple and economical description of general purely azimuthally dispersive wave systems.

Quasi-biennial oscillation and tracer distributions in a coupled chemistry-climate model

[ 1] We have used a fully coupled chemistry-climate model (CCM), which generates its own wind and temperature quasi-biennial oscillation (QBO), to study the effect of coupling on the QBO and to examine the QBO signals in stratospheric trace gases, particularly ozone. Radiative coupling of the interactive chemistry to the underlying general circulation model tends to prolong the QBO period and to increase the QBO amplitude in the equatorial zonal wind in the lower and middle stratosphere. The model ozone QBO agrees well with Stratospheric Aerosol and Gas Experiment II and Total Ozone Mapping Spectrometer satellite observations in terms of vertical and latitudinal structure. The model captures the ozone QBO phase change near 28 km over the equator and the column phase change near +/- 15 degrees latitude. Diagnosis of the model chemical terms shows that variations in NOx are the main chemical driver of the O-3 QBO around 35 km, i.e., above the O-3 phase change.

Precision of quasi-optical null-balanced bridge techniques for transmission and reflection coefficient measurements

The precision of quasioptical null-balanced bridge instruments for transmission and reflection coefficient measurements at millimeter and submillimeter wavelengths is analyzed. A Jones matrix analysis is used to describe the amount of power reaching the detector as a function of grid angle orientation, sample transmittance/reflectance and phase delay. An analysis is performed of the errors involved in determining the complex transmission and reflection coefficient after taking into account the quantization error in the grid angle and micrometer readings, the transmission or reflection coefficient of the sample, the noise equivalent power of the detector, the source power and the post-detection bandwidth. For a system fitted with a rotating grid with resolution of 0.017 rad and a micrometer quantization error of 1 μm, a 1 mW source, and a detector with a noise equivalent power 5×10−9 W Hz−1/2, the maximum errors at an amplitude transmission or reflection coefficient of 0.5 are below ±0.025.

Response of the northern stratospheric polar vortex to the seasonal alignment of QBO phase transitions

This study considers the strength of the Northern Hemisphere Holton-Tan effect (HTE) in terms of the phase alignment of the quasi-biennial oscillation (QBO) with respect to the annual cycle. Using the ERA-40 Reanalysis, it is found that the early winter (Nov–Dec) and late winter (Feb–Mar) relation between QBO phase and the strength of the stratospheric polar vortex is optimized for subsets of the 44-year record that are chosen on the basis of the seasonality of QBO phase transitions at the 30 hPa level. The timing of phase transitions serves as a proxy for changes in the vertical structure of the QBO over the whole depth of the tropical stratosphere. The statistical significance of the Nov–Dec (Feb–Mar) HTE is greatest when 30 hPa QBO phase transitions occur 9–14 (4–9) months prior to the January of the NH winter in question. This suggests that there exists for both early and late winter a vertical structure of tropical stratospheric winds that is most effective at influencing the interannual variability of the polar vortex, and that an early (late) winter HTE is associated with an early (late) progression of QBO phase towards that structure. It is also shown that the seasonality of QBO phase transitions at 30 hPa varies on a decadal timescale, with transitions during the first half of the calendar year being relatively more common during the first half of the tropical radiosonde wind record. Combining these two results suggests that decadal changes in HTE strength could result from the changing seasonality of QBO phase transitions. Citation: Anstey, J. A., and T. G. Shepherd (2008), Response of the northern stratospheric polar vortex to the seasonal alignment of QBO phase transitions, Geophys. Res. Lett., 35, L22810, doi:10.1029/2008GL035721.

Interferometric technique for measuring terahertz antenna phase patterns

A quasi-optical interferometric technique capable of measuring antenna phase patterns without the need for a heterodyne receiver is presented. It is particularly suited to the characterization of terahertz antennas feeding power detectors or mixers employing quasi-optical local oscillator injection. Examples of recorded antenna phase patterns at frequencies of 1.4 and 2.5 THz using homodyne detectors are presented. To our knowledge, these are the highest frequency antenna phase patterns ever recovered. Knowledge of both the amplitude and phase patterns in the far field enable a Gauss-Hermite or Gauss-Laguerre beam-mode analysis to be carried out for the antenna, of importance in performance optimization calculations, such as antenna gain and beam efficiency parameters at the design and prototype stage of antenna development. A full description of the beam would also be required if the antenna is to be used to feed a quasi-optical system in the near-field to far-field transition region. This situation could often arise when the device is fitted directly at the back of telescopes in flying observatories. A further benefit of the proposed technique is simplicity for characterizing systems in situ, an advantage of considerable importance as in many situations, the components may not be removable for further characterization once assembled. The proposed methodology is generic and should be useful across the wider sensing community, e.g., in single detector acoustic imaging or in adaptive imaging array applications. Furthermore, it is applicable across other frequencies of the EM spectrum, provided adequate spatial and temporal phase stability of the source can be maintained throughout the measurement process. Phase information retrieval is also of importance to emergent research areas, such as band-gap structure characterization, meta-materials research, electromagnetic cloaking, slow light, super-lens design as well as near-field and virtual imaging applications.

Influence of the quasi-biennial oscillation on the extratropical winter stratosphere in an atmospheric general circulation model and in reanalysis data

The interannual variability of the stratospheric polar vortex during winter in both hemispheres is observed to correlate strongly with the phase of the quasi-biennial oscillation (QBO) in tropical stratospheric winds. It follows that the lack of a spontaneously generated QBO in most atmospheric general circulation models (AGCMs) adversely affects the nature of polar variability in such models. This study examines QBO–vortex coupling in an AGCM in which a QBO is spontaneously induced by resolved and parameterized waves. The QBO–vortex coupling in the AGCM compares favorably to that seen in reanalysis data [from the 40-yr ECMWF Re-Analysis (ERA-40)], provided that careful attention is given to the definition of QBO phase. A phase angle representation of the QBO is employed that is based on the two leading empirical orthogonal functions of equatorial zonal wind vertical profiles. This yields a QBO phase that serves as a proxy for the vertical structure of equatorial winds over the whole depth of the stratosphere and thus provides a means of subsampling the data to select QBO phases with similar vertical profiles of equatorial zonal wind. Using this subsampling, it is found that the QBO phase that induces the strongest polar vortex response in early winter differs from that which induces the strongest late-winter vortex response. This is true in both hemispheres and for both the AGCM and ERA-40. It follows that the strength and timing of QBO influence on the vortex may be affected by the partial seasonal synchronization of QBO phase transitions that occurs both in observations and in the model. This provides a mechanism by which changes in the strength of QBO–vortex correlations may exhibit variability on decadal time scales. In the model, such behavior occurs in the absence of external forcings or interannual variations in sea surface temperatures.

Climatology and Forcing of the Quasi-Biennial Oscillation in the MAECHAM5 Model

The quasi-biennial oscillation (QBO) in the equatorial zonal wind is an outstanding phenomenon of the atmosphere. The QBO is driven by a broad spectrum of waves excited in the tropical troposphere and modulates transport and mixing of chemical compounds in the whole middle atmosphere. Therefore, the simulation of the QBO in general circulation models and chemistry climate models is an important issue. Here, aspects of the climatology and forcing of a spontaneously occurring QBO in a middle-atmosphere model are evaluated, and its influence on the climate and variability of the tropical middle atmosphere is investigated. Westerly and easterly phases are considered separately, and 40-yr ECMWF Re-Analysis (ERA-40) data are used as a reference where appropriate. It is found that the simulated QBO is realistic in many details. Resolved large-scale waves are particularly important for the westerly phase, while parameterized gravity wave drag is more important for the easterly phase. Advective zonal wind tendencies are important for asymmetries between westerly and easterly phases, as found for the suppression of the easterly phase downward propagation. The simulation of the QBO improves the tropical upwelling and the atmospheric tape recorder compared to a model without a QBO. The semiannual oscillation is simulated realistically only if the QBO is represented. In sensitivity tests, it is found that the simulated QBO is strongly sensitive to changes in the gravity wave sources. The sensitivity to the tested range of horizontal resolutions is small. The stratospheric vertical resolution must be better than 1 km to simulate a realistic QBO.