This PhD thesis deals with phenomena which are closely related to the unique thermal structure of the polar summer mesosphere, namely Polar Mesosphere Summer Echoes (PMSE). PMSE are strong radar echoes commonly observed by VHF MST radars from thin layers in the 80-90 km altitude interval at high latitudes during summer. They follow a seasonal pattern of abrupt appearance in late May and a gradual disappearance in mid-August. This period corresponds roughly to the time between the completion of the summer time cooling of the polar mesopause to the time of reversal of the mesospheric circulation to autumn condition. In this connection, PMSE are associated with the extremely low temperatures, i.e. below 140 K, which are unique to the polar summer mesopause. Traditional theories of radar (partial) reflection and scattering have been unable to explain the PMSE and the exact mechanism for their occurrence remains unclear despite the steadily increasing interest in them over the past 20 years. Currently accepted theories regarding the mechanism giving rise to PMSE agree that one of the conditions needed for enhanced radar echoes is the presence of low-mobility charge carries such as large cluster ions and ice aerosols which capture the ambient electrons. It has been established that the PMSE are in some way associated with noctilucent clouds (NLC), layers of ice crystals, which constitute the highest observed clouds in the earth’s atmosphere. PMSE occurrence and dynamics are also found to be closely connected with the planetary and gravity waves. Observations of PMSE presented in this thesis have been carried out by the Esrange MST radar (ESRAD) located at Esrange (67°56’N, 21°04’E) just outside Kiruna in northernmost Sweden. The radar operates at 52 MHz with 72 kW peak power and a maximum duty cycle of 5%. The antenna consists of 12x12 array of 5-element Yagis with a 0.7l spacing. During the PMSE measurements the radar used a 16-bit complementary code having a baud length of 1mS. This corresponds to height resolution of 150 m. The sampling frequency was set at 1450 Hz. The covered height range was 80-90 km. The presence of PMSE was determined on the basis of the radar SNR (signal-to-noise ratio). The PMSE measurements have been made during May-August each year since 1997. PMSE seasonal and diurnal occurrence rates as well as dynamics have been studied in connection with tidal winds, planetary waves, temperature and water vapor content in the mesosphere (Papers I, IV and VI). Simultaneous and common-volume observations of PMSE and noctilucent clouds have been performed by radar, lidar and CCD camera (Paper V). Correlation between variations in PMSE and variations in extra ionization added by precipitating energetic electrons or high-energy particles from the Sun has been examined (Papers II and III). Possible influence of transport effects due to the electric field on PMSE appearance has been studied during a solar proton event (Paper III).

This thesis examines the utility of mathematical complexity measures for the analysis of stratospheric dynamics. Through theoretical considerations and tests with artificial data sets, e.g., the iteration of the logistic map, suitable parameters are determined for the application of the statistical entropy measures sample entropy (SE) and Rényi entropy (RE) to methane (a long-lived stratospheric tracer) data from simulations of the SOCOL chemistry-climate model. The SE is shown to be useful for quantifying the variability of recurring patterns in a time series and is able to identify tropical patterns similar to those reported by previous studies of the ``tropical pipe'' region. However, the SE is found to be unsuitable for use in polar regions, due to the non-stationarity of the methane data at extra-tropical latitudes. It is concluded that the SE cannot be used to analyse climate complexity on a global scale. The focus is turned to the RE, which is a complexity measure of probability distribution functions (PDFs). Using the second order RE and a normalisation factor, zonal PDFs of ten consecutive days of methane data are created with a Bayesian optimal binning technique. From these, the RE is calculated for every day (moving 10-day window). The results indicate that the RE is a promising tool for identifying stratospheric mixing barriers. In Southern Hemisphere winter and early spring, RE produces patterns similar to those found in other studies of stratospheric mixing. High values of RE are found to be indicative of the strong fluctuations in tracer distributions associated with relatively unmixed air in general, and with gradients in the vicinity of mixing barriers, in particular. Lower values suggest more thoroughly mixed air masses. The analysis is extended to eleven years of model data. Realistic inter-annual variability of some of the RE structures is observed, particularly in the Southern Hemisphere. By calculating a climatological mean of the RE for this period, additional mixing patterns are identified in the Northern Hemisphere. The validity of the RE analysis and its interpretation is underlined by showing that qualitatively similar patterns can be seen when using observational satellite data of a different tracer. Compared to previous techniques, the RE has the advantage that it requires significantly less computational effort, as it can be used to derive dynamical information from model or measurement tracer data without relying on any additional input such as wind fields. The results presented in this thesis strongly suggest that the RE is a useful new metric for analysing stratospheric mixing and its variability from climate model data. Furthermore, it is shown that the RE measure is very robust with respect to data gaps, which makes it ideal for application to observations. Hence, using the RE for comparing observations of tracer distributions with those from model simulations potentially presents a novel approach for analysing mixing in the stratosphere.

Tempestades severas ocorrem na Região Metropolitana de São Paulo (RMSP) principalmente nos meses quentes e úmidos do ano. Nesta tese, os mecanismos de disparo da convecção profunda são estudados por meio de análises observacionais e simulações numéricas com o Advanced Regional Prediction System (ARPS). A metodologia proposta compreende o uso da parametrização microfísica fria na simulação dos processos físicos que levam à formação de nuvens cumulonimbus, sem o uso da parametrização de cúmulos nas grades de altíssima resolução espacial. Nos eventos estudados, as primeiras células de precipitação observadas e simuladas surgiram em razão da interação entre o escoamento atmosférico na camada limite planetária e a topografia local. As células secundárias foram geralmente mais intensas, uma vez que elas surgiram após o aquecimento diabático adicional. O mecanismo de disparo das células secundárias foi a corrente ascendente induzida pela propagação horizontal das frentes de rajada em baixos níveis da atmosfera das correntes descendentes das células primárias. As frentes de rajada tiveram velocidade de propagação horizontal típica de 6 m s-1. No evento de 02 de fevereiro de 2004, células convectivas profundas foram simuladas com alto grau de realismo no domínio da grade de 3 km de resolução espacial. Observou-se que, neste caso, a frente de brisa marítima pôde atuar como guia de ondas para a colisão entre duas frentes de rajada. A propagação da frente de brisa marítima para o interior do continente ocorreu em conjunção a um forte gradiente de vapor dágua nos níveis mais baixos da troposfera. As células convectivas profundas secundárias surgiram e se desenvolveram exatamente nesta zona de interface, a qual representa o contraste entre as diferentes massas de ar marítima e continental. No evento de 04 de fevereiro de 2004, na grade de 1 km de resolução, a análise objetiva com as medidas das estações de superfície na RMSP correspondente às 1800 UTC indicou a presença de uma ilha de calor urbana com até 4 oC de aquecimento diferencial entre a Capital e vizinhanças. O principal efeito da assimilação destas medidas foi a redução do NCL em até 80 hPa, o que favoreceu o disparo da convecção naquela área. / Severe thunderstorms occur in the Metropolitan Area of São Paulo (MASP) mainly in the warm and wet months of the year. In this work, the triggering mechanisms of deep convection are studied through observed data and numerical simulations with the Advanced Regional Prediction System (ARPS). The proposed methodology focuses in the use of microphysics parameterization of cold clouds to simulate physical process linked to the life cycle of thunderstorms. The cumulus cloud parameterization isnt used in high resolution numerical grids. In the real case studies, both observed and simulated, early convective cells developed as a consequence of the interaction between the planetary boundary layer atmospheric flow and the local topography. The secondary convective cells were generally strongest, once they developed after additional surface diabatic heating. The triggering mechanism of these secondary cells was the updraft induced by gust fronts generated by downdrafts of primary cells. The gust fronts had a typical horizontal propagation velocity of 6 m s-1. In the February 02 2004 event, deep convective cells were simulated with high degree of realism with a 3 km resolution grid. It was observed that, in this case, the sea-breeze front could act as a wave guide to the collision between two different gust fronts. In addition, the sea breeze front propagated to the continental area together with a strong low level water vapor gradient. The secondary deep convective cells arose and developed exactly on this interface zone, which represents the contrast between the oceanic and continental air masses. The interface zone was marked by a water vapor mixing rate of 14 g kg-1. In the February 04 2004 event, the objective analysis, made with some MASP´s surface stations measurements at 1800 UTC in the 1 km resolution grid, indicates the presence of an urban heat island with up to 4 oC of differential heating between São Paulo city and its neighboring area. The main effect in assimilating these surface measurements was the lowering of the lift condensation level up to 80 hPa, which favored the triggering of convection in that area.

The effects of a non-zero background state on atmospheric dynamics is explored through simple models and observations. Firstly, we examine the effects of moisture gradients on the stability and propagation of Rossby waves in a mid-latitude -plane. We begin by a consistent derivation of the forced quasi-geostrophic equations on a -plane to understand the constraints placed by geostrophy on the time scale of condensation. We see that the presence of meridional gradients of moisture results in a slowdown of the waves. On the introduction of zonal gradients of moisture, the waves become unstable, and for certain parameters which are representative of the real atmosphere, they propagate eastward and mature on an intra-seasonal timescale. The mechanism of the in hence of moisture on waves is understood by thinking of condensation as providing an \equivalent" potential vorticity (PV) gradient which opposes the dynamical PV gradient. Secondly, we look at the effects of a mean background ow on the Matsuno-Gill response in the spherical shallow water system. The mean ow is prescribed to resemble the climatological upper tropospheric zonal wind structure in the atmosphere. As the strength of the ow increases, the equatorially trapped Matsuno-Gill response rst transforms into a poleward propagating Rossby wavetrain. As the strength of the mean ow reaches values similar to that observed in the atmosphere, the stationary wave response becomes a zonally oriented quadrupole structure. This structure bears a striking resemblance to the observed upper level structure of the Madden-Julian oscillation (MJO). The time evolution of this quadrupole structure is quick enough to be relevant on MJO timescales, and the structure is quite robust across a range of values for the drag coefficient. Finally, we look at the role played by low frequency variability in the Pacific in the recent expansion of the Hadley cell. We find that the dominant effect of the low frequency variability is a stationary dispersive Rossby wavetrain extending from the tropical Paci. We further find that most of the observed expansion of the Hadley cell can be accounted for by this low frequency variability. We nd that large scale changes such as the changes in the equator-pole temperature gradient or midlatitude static stability need not be invoked to understand the observed expansion.

<p dir="ltr">It is known that tornadic storms favor environments characteristic of high values of thermal instability, adequate vertical wind shear, abundant near-surface moisture supply, and strong storm-relative helicity at the lowest 1-km boundary layer. These mesoscale environmental conditions and associated storm behaviors are strongly governed by large-scale synoptic patterns and sensitive to variabilities in near-surface characteristics, which are less known in the current research community. This study aims to advance the relatively underexplored area regarding the interaction between surface characteristics, mesoscale environmental conditions, and large-scale synoptic patterns driving tornadic storms in the U.S. </p><p dir="ltr">We first investigate the impact of surface drag on the structure and evolution of these boundaries, their associated distribution of near-surface vorticity, and tornadogenesis and maintenance. Comparisons between idealized simulations without and with drag introduced in the mature stage of the storm prior to tornadogenesis reveal that the inclusion of surface drag substantially alters the low-level structure, particularly with respect to the number, location, and intensity of surface convergence boundaries. Substantial drag-generated horizontal vorticity induces rotor structures near the surface associated with the convergence boundaries in both the forward and rear flanks of the storm. Stretching of horizontal vorticity and subsequent tilting into the vertical along the convergence boundaries lead to elongated positive vertical vorticity sheets on the ascending branch of the rotors and the opposite on the descending branch. The larger near-surface pressure deficit associated with the faster development of the near-surface cyclone when drag is active creates a downward dynamic vertical pressure gradient force that suppresses vertical growth, leading to a weaker and wider tornado detached from the surrounding convergence boundaries. A conceptual model of the low-level structure of the tornadic supercell is presented that focuses on the contribution of surface drag, with the aim of adding more insight and complexity to previous conceptual models.</p><p dir="ltr">We then examine the behaviors and dynamics of TLVs in response to a range of surface drag strengths in idealized simulations and explore their sensitivities to different storm environments. We find that the contribution of surface drag on TLV development is strongly governed by the interaction between surface rotation, surface convergence boundaries, and the low-level mesocyclone. Surface drag facilitates TLV formation by enhancing near-surface vortices and low-level lifting, mitigating the need for an intense updraft gradient developing close to the ground. As surface drag increases, a wider circulation near the surface blocks the inflow from directly reaching the rotating core, leading to a less tilted structure that allows the TLV position beneath the pressure minima aloft. Further increase in drag strength discourages TLV intensification by suppressing vertical stretching due to a negative vertical pressure perturbation gradient force, and it stops benefiting from the support of surrounding convergence boundaries and the overlying low-level updraft, instead becoming detached from them. We hence propose a favorable condition for TLV formation and duration where a TLV forms a less tilted structure directly beneath the low-level mesocyclone but also evolves near surrounding surface boundaries, which scenario strongly depends on underlying surface drag strength. </p><p dir="ltr">Beyond near-surface characteristics, we further explore how these storm-favorable environmental conditions may interact with the larger-scale synoptic patterns and how these interactions may affect the tornadic storm potential in the current warming climate. We employ hierarchical clustering analysis to classify the leading synoptic patterns driving tornadic storms across different geographic regions in the U.S. We find that the primary synoptic patterns are distinguishable across geographic regions and seasonalities. The intense upper-level jet streak described by the high values of eddy kinetic energy (EKE) associated with the dense distribution of Z500 contours dominates the tornado events in the southeast U.S. in the cold season (November-March). Late Spring and early Summer Tornado events in the central and south Great Plains are dominated by deep trough systems to the west axes of the tornado genesis position, while more summer events associated with weak synoptic forcing are positioned closer to the lee side of Rocky Mountain. Moreover, the increasing trend in tornado frequency in the southeastern U.S. is mainly driven by synoptic patterns with intense forcing, and the decreasing trends in portions of the Great Plains are associated with weak synoptic forcing. This finding indicates that the physical mechanisms driving the spatial trends of tornado occurrences differ across regions in the U.S.</p>

L’un des objectifs de cette thèse est d’interpréter les observations des principaux hydrocarbures(C2H2 et C2H6) effectuées par Cassini (NASA/ESA) sur Jupiter et Saturne. Les modèles photochimiques à une dimension sont insuffisants pour interpréter ces observations spatialement résolues. J’ai développé le premier modèle photochimique saisonnier à deux dimensions (altitude-latitude) des planètes géantes qui calcule leur composition chimique.En l’absence de transport méridional, la composition chimique de Saturne suit les variations d’ensoleillement. Les abondances de C2H2 et C2H6 mesurées par Cassini (Guerletet al., 2009) sont reproduites jusqu’aux latitudes moyennes, à des pressions supérieures à0,1mbar. Les écarts notés dans l’hémisphère sud suggèrent la présence de dynamique ou d’une chimie entre les ions et les espèces neutres. J’ai couplé, pour la première fois, mon modèle photochimique avec le modèle radiatif de Greathouse et al. (2008). Nous prédisons un décalage du pic saisonnier de température, par rapport aux précédents modèles, d’une demi-saison à haute altitude et aux hautes latitudes.Jupiter présente de faibles variations saisonnières de composition chimique, uniquement contrôlées par son excentricité. Les distributions méridionales observées de C2H2 etC2H6 présentent des tendances opposées (Nixon et al., 2010). Mon modèle est en accord avec les observations de C2H6 lorsque j’invoque une combinaison de diffusion méridionale et de circulation stratosphérique, tout en provoquant un plus grand désaccord avec les observations de C2H2. La chimie ionique pourrait principalement affecter C2H2 et jouer un rôle important dans l’atmosphère de Jupiter. / One of the goals of this thesis is to interpret the observations of the main hydrocarbons(C2H2 and C2H6) from Cassini (NASA/ESA) on Jupiter and Saturn. The one-dimensional photochemical models are insufficient to explain these spatially resolved observations. I have developed the first two-dimensional (altitude-latitude) seasonal photochemical model for the giant planets, which predicts their chemical composition.Without meridional transport, Saturn’s chemical composition follows the insolation variations. The C2H2 and C2H6 abundances measured by Cassini (Guerlet et al., 2009)are reproduced from the equator up to mid-latitudes, at pressures higher than 0.1mbar.At higher latitudes, the disagreements suggest either a stratospheric circulation cell orthe signature of ion-neutral chemistry. For the first time, I have coupled our seasonal photochemical model with the seasonal radiative model of Greathouse et al. (2008). I predict that the seasonal temperature peak is shifted half a season earlier, with respect to previous models, at high latitudes in the higher stratosphere.Jupiter shows weak seasonal variations of chemical composition, only controlled by its orbital eccentricity. The observed meridional distributions of C2H2 and C2H6 show opposition trends (Nixon et al., 2010). C2H6 observed distribution is reproduced when Isuppose a combination of meridional diffusion and stratospheric circulation, while causingat the same time a stronger agreement with the C2H2 observations. Accounting for theion-neutral chemistry might preferentially affect C2H2 and potentially play a key role on hydrocarbon abundances in Jupiter’s stratosphere.