Dynamical control of the mesosphere by orographic and non-orographic gravity wave drag during the extended northern winters of 2006 and 2009

A version of the Canadian Middle Atmosphere Model (CMAM) that is nudged toward reanalysis data up to 1 hPa is used to examine the impacts of parameterized orographic and non-orographic gravity wave drag (OGWD and NGWD) on the zonal-mean circulation of the mesosphere during the extended northern winters of 2006 and 2009 when there were two large stratospheric sudden warmings. The simulations are compared to Aura Microwave Limb Sounder (MLS) observations of mesospheric temperature, carbon monoxide (CO) and derived zonal winds. The control simulation, which uses both OGWD and NGWD, is shown to be in good agreement with MLS. The impacts of OGWD and NGWD are assessed using simulations in which those sources of wave drag are removed. In the absence of OGWD the mesospheric zonal winds in the months preceding the warmings are too strong, causing increased mesospheric NGWD, which drives excessive downwelling, resulting in overly large lower mesospheric values of CO prior to the warming. NGWD is found to be most important following the warmings when the underlying westerlies are too weak to allow much vertical propagation of the orographic gravity waves to the mesosphere. NGWD is primarily responsible for driving the circulation that results in the descent of CO from the thermosphere following the warmings. Zonal mean mesospheric winds and temperatures in all simulations are shown to be strongly constrained by (i.e. slaved to) the stratosphere. Finally, it is demonstrated that the responses to OGWD and NGWD are non-additive due to their dependence and influence on the background winds and temperatures.

Effects of spatial resolution in the simulation of daily and subdaily precipitation in the southwestern US

We evaluate the effects of spatial resolution on the ability of a regional climate model to reproduce observed extreme precipitation for a region in the Southwestern United States. A total of 73 National Climate Data Center observational sites spread throughout Arizona and New Mexico are compared with regional climate simulations at the spatial resolutions of 50 km and 10 km for a 31 year period from 1980 to 2010. We analyze mean, 3-hourly and 24-hourly extreme precipitation events using WRF regional model simulations driven by NCEP-2 reanalysis. The mean climatological spatial structure of precipitation in the Southwest is well represented by the 10 km resolution but missing in the coarse (50 km resolution) simulation. However, the fine grid has a larger positive bias in mean summer precipitation than the coarse-resolution grid. The large overestimation in the simulation is in part due to scale-dependent deficiencies in the Kain-Fritsch convective parameterization scheme that generate excessive precipitation and induce a slow eastward propagation of the moist convective summer systems in the high-resolution simulation. Despite this overestimation in the mean, the 10 km simulation captures individual extreme summer precipitation events better than the 50 km simulation. In winter, however, the two simulations appear to perform equally in simulating extremes.

On the longitudinal extent of magnetopause reconnection pulses

The open magnetosphere model of cusp ion injection, acceleration and precipitation is used to predict the dispersion characteristics for fully pulsed magnetic reconnection at a low-latitude magnetopause X-line. The resulting steps, as would be seen by a satellite moving meridionally and normal to the ionospheric projection of the X-line, are compared with those seen by satellites moving longitudinally, along the open/closed boundary. It is shown that two observed cases can be explained by similar magnetosheath and reconnection characteristics, and that the major differences between them are well explained by the different satellite paths through the events. Both cases were observed in association with poleward-moving transient events seen by ground-based radar, as also predicted by the theory. The results show that the reconnection is pulsed but strongly imply it cannot also be spatially patchy, in the sense of isolated X-lines which independently are intermittently active. Furthermore they show that the reconnection pulses responsible for the poleward-moving events and the cusp ion steps, must cover at least 3 h of magnetic local time, although propagation of the active reconnection region may mean that it does not extend this far at any one instant of time.

Vertical structure and physical processes of the Madden-Julian oscillation: exploring key model physics in climate simulations

Aimed at reducing deficiencies in representing the Madden-Julian oscillation (MJO) in general circulation models (GCMs), a global model evaluation project on vertical structure and physical processes of the MJO was coordinated. In this paper, results from the climate simulation component of this project are reported. It is shown that the MJO remains a great challenge in these latest generation GCMs. The systematic eastward propagation of the MJO is only well simulated in about one-fourth of the total participating models. The observed vertical westward tilt with altitude of the MJO is well simulated in good MJO models, but not in the poor ones. Damped Kelvin wave responses to the east of convection in the lower troposphere could be responsible for the missing MJO preconditioning process in these poor MJO models. Several process-oriented diagnostics were conducted to discriminate key processes for realistic MJO simulations. While large-scale rainfall partition and low-level mean zonal winds over the Indo-Pacific in a model are not found to be closely associated with its MJO skill, two metrics, including the low-level relative humidity difference between high and low rain events and seasonal mean gross moist stability, exhibit statistically significant correlations with the MJO performance. It is further indicated that increased cloud-radiative feedback tends to be associated with reduced amplitude of intraseasonal variability, which is incompatible with the radiative instability theory previously proposed for the MJO. Results in this study confirm that inclusion of air-sea interaction can lead to significant improvement in simulating the MJO.

Equation for the microwave backscatter cross section of aggregate snowflakes using the self-similar Rayleigh–Gans approximation

In this paper an equation is derived for the mean backscatter cross section of an ensemble of snowflakes at centimeter and millimeter wavelengths. It uses the Rayleigh–Gans approximation, which has previously been found to be applicable at these wavelengths due to the low density of snow aggregates. Although the internal structure of an individual snowflake is random and unpredictable, the authors find from simulations of the aggregation process that their structure is “self-similar” and can be described by a power law. This enables an analytic expression to be derived for the backscatter cross section of an ensemble of particles as a function of their maximum dimension in the direction of propagation of the radiation, the volume of ice they contain, a variable describing their mean shape, and two variables describing the shape of the power spectrum. The exponent of the power law is found to be −. In the case of 1-cm snowflakes observed by a 3.2-mm-wavelength radar, the backscatter is 40–100 times larger than that of a homogeneous ice–air spheroid with the same mass, size, and aspect ratio.

Rossby wave propagation on potential vorticity fronts with finite width

The horizontal gradient of potential vorticity (PV) across the tropopause typically declines with lead time in global numerical weather forecasts and tends towards a steady value dependent on model resolution. This paper examines how spreading the tropopause PV contrast over a broader frontal zone affects the propagation of Rossby waves. The approach taken is to analyse Rossby waves on a PV front of finite width in a simple single-layer model. The dispersion relation for linear Rossby waves on a PV front of infinitesimal width is well known; here an approximate correction is derived for the case of a finite width front, valid in the limit that the front is narrow compared to the zonal wavelength. Broadening the front causes a decrease in both the jet speed and the ability of waves to propagate upstream. The contribution of these changes to Rossby wave phase speeds cancel at leading order. At second order the decrease in jet speed dominates, meaning phase speeds are slower on broader PV fronts. This asymptotic phase speed result is shown to hold for a wide class of single-layer dynamics with a varying range of PV inversion operators. The phase speed dependence on frontal width is verified by numerical simulations and also shown to be robust at finite wave amplitude, and estimates are made for the error in Rossby wave propagation speeds due to the PV gradient error present in numerical weather forecast models.

ENSO related variation of equatorial MRG and Rossby waves and forcing from higher latitudes

The contrasting behaviour of westward-moving mixed Rossby-gravity (WMRG) and the first Rossby (R1) waves in El Niño (EN) and La Niña (LN) seasons is documented with a focus on the Northern Hemisphere winter. The eastward-moving variance in the upper troposphere is dominated by WMRG and R1 structures that appear to be Doppler-shifted by the flow and are referred to as WMRG-E and R1-E. In the East Pacific and Atlantic the years with stronger equatorial westerly winds have the stronger WMRG and WMRG- E. In the East Pacific, R1 is also a maximum in LN. However, R1-E exhibits an eastward-shift between LN and EN. The changes with ENSO phase provide a test-bed for the understanding of these waves. In the East Pacific and Atlantic, the stronger WMRG-E and WMRG with stronger westerlies are in accord with the dispersion relation with simple Doppler-shifting by the zonal flow. The possible existence of free waves can also explain stronger R1 in EN in the Eastern Hemisphere. 1-D free wave propagation theory based on wave activity conservation is also important for R1. However, this theory is unable to explain the amplitude maxima for other waves observed in the strong equatorial westerly regions in the Western Hemisphere, and certainly not their ENSO-related variation. The forcing of equatorial waves by higher latitude wave activity and its variation with ENSO phase is therefore examined. Propagation of extratropical eastward-moving Rossby wave activity through the westerly ducts into the equatorial region where it triggers WMRG-E is favoured in the stronger westerlies, in LN in the East Pacific and EN in the Atlantic. It is also found that WMRG is forced by Southern Hemisphere westward-moving wavetrains arching into the equatorial region where they are reflected. The most significant mechanism for both R1 and R1-E appear to be lateral forcing by subtropical wavetrains.

Rossby wave propagation tracks in southern hemisphere mean basic flows associated to generalized frosts over southern South America

Based on previous observational studies on cold extreme events over southern South America, some recent studies suggest a possible relationship between Rossby wave propagation remotely triggered and the occurrence of frost. Using the concept of linear theory of Rossby wave propagation, this paper analyzes the propagation of such waves in two different basic states that correspond to austral winters with maximum and minimum generalized frost frequency of occurrence in the Wet Pampa (central-northwest Argentina). In order to determine the wave trajectories, the ray tracing technique is used in this study. Some theoretical discussion about this technique is also presented. The analysis of the basic state, from a theoretical point of view and based on the calculation of ray tracings, corroborates that remotely excited Rossby waves is the mechanism that favors the maximum occurrence of generalized frosts. The basic state in which the waves propagate is what conditions the places where they are excited. The Rossby waves are excited in determined places of the atmosphere, propagating towards South America along the jet streams that act as wave guides, favoring the generation of generalized frosts. In summary, this paper presents an overview of the ray tracing technique and how it can be used to investigate an important synoptic event, such as frost in a specific region, and its relationship with the propagation of large scale planetary waves.

Cell coupling in mouse pancreatic beta-cells measured in intact islets of Langerhans

The perforated whole-cell configuration of the patch-clamp technique was applied to functionally identified beta-cells in intact mouse pancreatic islets to study the extent of cell coupling between adjacent beta-cells. Using a combination of current- and voltage-clamp recordings, the total gap junctional conductance between beta-cells in an islet was estimated to be 1.22 nS. The analysis of the current waveforms in a voltage-clamped cell ( due to the. ring of an action potential in a neighbouring cell) suggested that the gap junctional conductance between a pair of beta-cells was 0.17 nS. Subthreshold voltage-clamp depolarization (to -55 mV) gave rise to a slow capacitive current indicative of coupling between beta-cells, but not in non-beta-cells, with a time constant of 13.5 ms and a total charge movement of 0.2 pC. Our data suggest that a superficial beta-cell in an islet is in electrical contact with six to seven other beta-cells. No evidence for dye coupling was obtained when cells were dialysed with Lucifer yellow even when electrical coupling was apparent. The correction of the measured resting conductance for the contribution of the gap junctional conductance indicated that the whole-cell K(ATP) channel conductance (G(K,ATP)) falls from approximately 2.5 nS in the absence of glucose to 0.1 nS at 15 mM glucose with an estimated IC(50) of approximately 4 mM. Theoretical considerations indicate that the coupling between beta-cells within the islet is sufficient to allow propagation of [Ca(2+)](i) waves to spread with a speed of approximately 80 mu m s(-1), similar to that observed experimentally in confocal [Ca(2+)](i) imaging.

Nonlinear refractive index of RECOB (RE = Gd and La) crystals

The Z-scan technique is employed to obtain the nonlinear refractive index (n (2)) of the Ca(4)REO(BO(3))(3) (RECOB, where RE = Gd and La) single crystals using 30 fs laser pulses centered at 780 nm for the two orthogonal orientations determined by the optical axes (X and Z) relative to the direction of propagation of the laser beam (k//Y// crystallographic b-axis). The large values of n (2) indicate that both GdCOB and LaCOB are potential hosts for Yb:RECOB lasers operating in the Kerr-lens mode locking (KLM) regime.

Protein lethality investigated in terms of long range dynamical interactions

The relationship between network structure/dynamics and biological function constitutes a fundamental issue in systems biology. However, despite many related investigations, the correspondence between structure and biological functions is not yet fully understood. A related subject that has deserved particular attention recently concerns how essentiality is related to the structure and dynamics of protein interactions. In the current work, protein essentiality is investigated in terms of long range influences in protein-protein interaction networks by considering simulated dynamical aspects. This analysis is performed with respect to outward activations, an approach which models the propagation of interactions between proteins by considering self-avoiding random walks. The obtained results are compared to protein local connectivity. Both the connectivity and the outward activations were found to be strongly related to protein essentiality.

Multiple Pathways Analysis of Brain Functional Networks from EEG Signals: An Application to Real Data

In the present study, we propose a theoretical graph procedure to investigate multiple pathways in brain functional networks. By taking into account all the possible paths consisting of h links between the nodes pairs of the network, we measured the global network redundancy R (h) as the number of parallel paths and the global network permeability P (h) as the probability to get connected. We used this procedure to investigate the structural and dynamical changes in the cortical networks estimated from a dataset of high-resolution EEG signals in a group of spinal cord injured (SCI) patients during the attempt of foot movement. In the light of a statistical contrast with a healthy population, the permeability index P (h) of the SCI networks increased significantly (P < 0.01) in the Theta frequency band (3-6 Hz) for distances h ranging from 2 to 4. On the contrary, no significant differences were found between the two populations for the redundancy index R (h) . The most significant changes in the brain functional network of SCI patients occurred mainly in the lower spectral contents. These changes were related to an improved propagation of communication between the closest cortical areas rather than to a different level of redundancy. This evidence strengthens the hypothesis of the need for a higher functional interaction among the closest ROIs as a mechanism to compensate the lack of feedback from the peripheral nerves to the sensomotor areas.

Time-resolved study electronic and thermal contributions to the nonlinear refractive index of Nd3+: SBN laser crystals

The propagation of an optical beam through dielectric media induces changes in the refractive index, An, which causes self-focusing or self-defocusing. In the particular case of ion-doped solids, there are thermal and non-thermal lens effects, where the latter is due to the polarizability difference, Delta alpha, between the excited and ground states, the so-called population lens (PL) effect. PL is a pure electronic contribution to the nonlinearity, while the thermal lens (TL) effect is caused by the conversion of part of the absorbed energy into heat. In time-resolved measurements such as Z-scan and TL transient experiments, it is not easy to separate these two contributions to nonlinear refractive index because they usually have similar response times. In this work, we performed time-resolved measurements using both Z-scan and mode mismatched TL in order to discriminate thermal and electronic contributions to the laser-induced refractive index change of the Nd3+-doped Strontium Barium Niobate (SrxBa1-xNb2O6) laser crystal. Combining numerical simulations with experimental results we could successfully distinguish between the two contributions to An. (C) 2007 Elsevier B.V. All rights reserved.

Adiabatic propagation in potential structures

In this work the adiabatic approximation is applied to the propagation of matter waves in confined geometries like those experimentally realized in recent atom optical experiments. Adiabatic propagation along a channel is assumed not to mix the various transverse modes. Nonadiabatic corrections arise from the potential squeezing and bending. Here we investigate the effect of the former. Detailed calculations of two-dimensional propagation are carried out both exactly and in an adiabatic approximation. This offers the possibility to analyze the validity of adiabaticity criteria. A semiclassical (sc) approach, based on the sc Massey parameter is shown to be inadequate, and the diffraction due to wave effects must be included separately. This brings in the Fresnel parameter well known from optical systems. Using these two parameters, we have an adequate understanding of adiabaticity on the system analyzed. Thus quantum adiabaticity must also take cognizance of the intrinsic diffraction of matter waves.

Contagion in financial networks: a network theory and agent-based approaches to modeling the spread of risk in financial systems

Esta dissertação estuda a propagação de crises sobre o sistema financeiro. Mais especi- ficamente, busca-se desenvolver modelos que permitam simular como um determinado choque econômico atinge determinados agentes do sistema financeiro e apartir dele se propagam, transformando-se em um problema sistêmico. A dissertação é dividida em dois capítulos,além da introdução. O primeiro capítulo desenvolve um modelo de propa- gação de crises em fundos de investimento baseado em ciência das redes.Combinando dois modelos de propagação em redes financeiras, um simulando a propagação de perdas em redes bipartites de ativos e agentes financeiros e o outro simulando a propagação de perdas em uma rede de investimentos diretos em quotas de outros agentes, desenvolve-se um algoritmo para simular a propagação de perdas através de ambos os mecanismos e utiliza-se este algoritmo para simular uma crise no mercado brasileiro de fundos de investimento. No capítulo 2,desenvolve-se um modelo de simulação baseado em agentes, com agentes financeiros, para simular propagação de um choque que afeta o mercado de operações compromissadas.Criamos também um mercado artificial composto por bancos, hedge funds e fundos de curto prazo e simulamos a propagação de um choque de liquidez sobre um ativo de risco securitizando utilizado para colateralizar operações compromissadas dos bancos.