Numerical modeling of inhaled charged aerosol deposition in human airways

A new numerical modeling of inhaled charge aerosol has been developed based on a modified Weibel's model. Both the velocity profiles (slug and parabolic flows) and the particle distributions (uniform and parabolic distributions) have been considered. Inhaled particles are modeled as a dilute dispersed phase flow in which the particle motion is controlled by fluid force and external forces acting on particles. This numerical study extends the previous numerical studies by considering both space- and image-charge forces. Because of the complex computation of interacting forces due to space-charge effect, the particle-mesh (PM) method is selected to calculate these forces. In the PM technique, the charges of all particles are assigned to the space-charge field mesh, for calculating charge density. The Poisson's equation of the electrostatic potential is then solved, and the electrostatic force acting on individual particle is interpolated. It is assumed that there is no effect of humidity on charged particles. The results show that many significant factors also affect the deposition, such as the volume of particle cloud, the velocity profile and the particle distribution. This study allows a better understanding of electrostatic mechanism of aerosol transport and deposition in human airways.

Chain extension in electrospun polystyrene fibres: a SANS study

Small angle neutron scattering techniques were used to quantify the size and shape of the chain conformation in electrospun fibres of atactic polystyrene prepared from solutions in methyl ethyl ketone. Aligned arrays of fibres were collected onto a rotating collector with tangential velocity varying between 0 ms-1 and approximately 15 ms-1. The measured radii of gyration of the polystyrene chains were found to be slightly higher than that expected for samples prepared from solutions in the concentrated regime. The ratio of the radius of gyration parallel and perpendicular to the chain axis was found to be approximately 1.05 in contrast to the substantial macroscopic shape transformation intrinsic to electrospinning. When the tangential velocity of the rotating collector was greater than the flight velocity of the fibres (ca. 4 ms-1), a further extension of the polymer chains was observed with a ratio of the radii of gyration increasing to 1.20 at the highest collector speeds. It is proposed that the heterogeneous processes involved, particularly solvent evaporation and the formation of a polymer skin during electrospinning play a significant role in determining the level of molecular anisotropy in the fibres.

Surface activity of surfactin recovered and purified from fermentation broth using a two-step ultrafiltration (UF) process

B. subtilis under certain types of media and fermentation conditions can produce surfactin, a biosurfactant which belongs to the lipopeptide class. Surfactin has exceptional surfactant activity, and exhibits some interesting biological characteristics such as antibacterial activity, antitumoral activity against ascites carcinoma cells, and a hypocholesterolemic activity that inhibits cAMP phosphodiesterase, as well as having anti-HIV properties. A cost effective recovery and purification of surfactin from fermentation broth using a two-step ultrafiltration (UF) process has been developed in order to reduce the cost of surfactin production. In this study, competitive adsorption of surfactin and proteins at the air-water interface was studied using surface pressure measurements. Small volumes of bovine serum albumin (BSA) and β-casein solutions were added to the air-water interface on a Langmuir trough and allowed to stabilise before the addition of surfactin to the subphase. Contrasting interfacial behaviour of proteins was observed with β-casein showing faster initial adsorption compared to BSA. On introduction of surfactin both proteins were displaced but a longer time were taken to displace β-casein. Overall the results showed surfactin were highly surface-active by forming a β-sheet structure at the air-water interface after reaching its critical micelle concentration (CMC) and were effective in removing both protein films, which can be explained following the orogenic mechanism. Results showed that the two-step UF process was effective to achieve high purity and fully functional surfactin.

Development of orientation during electrospinning of fibres of poly(ε-caprolactone)

We explore the influence of a rotating collector on the internal structure of poly(ε-caprolactone) fibres electrospun from a solution in dichloroethane. We find that above a threshold collector speed, the mean fibre diameter reduces as the speed increases and the fibres are further extended. Small-angle and wide-angle X-ray scattering techniques show a preferred orientation of the lamellar crystals normal to the fibre axis which increases with collector speed to a maximum and then reduces. We have separated out the processes of fibre alignment on the collector and the orientation of crystals within the fibres. There are several stages to this behaviour which correspond to the situations (a) where the collector speed is slower than the fibre spinning rate, (b) the fibre is mechanically extended by the rotating collector and (c) where the deformation leads to fibre fracture. The mechanical deformation leads to a development of preferred orientation with extension which is similar to the prediction of the pseudo-affine deformation model and suggests that the deformation takes place during the spinning process after the crystals have formed.

Tropical-extratropical interactions over Southern Africa: three cases of heavy summer season rainfall

The synoptic evolution of three tropical–extratropical (TE) interactions, each responsible for extreme rainfall events over southern Africa, is discussed in detail. Along with the consideration of previously studied events, common features of these heavy rainfall producing tropical temperate troughs (TTTs) over southern Africa are discussed. It is found that 2 days prior to an event, northeasterly moisture transports across Botswana, set up by the Angola low, are diverted farther south into the semiarid region of subtropical southern Africa. The TTTs reach full maturity as a TE cloud band, rooted in the central subcontinent, which is triggered by upper-level divergence along the leading edge of an upper-tropospheric westerly wave trough. Convection and rainfall within the cloud band is supported by poleward moisture transports with subtropical air rising as it leaves the continent and joins the midlatitude westerly flow. It is shown that these systems fit within a theoretical framework describing similar TE interactions found globally. Uplift forcing for the extreme rainfall of each event is investigated. Unsurprisingly, quasigeostrophic uplift is found to dominate in the midlatitudes with convective processes strongest in the subtropics. Rainfall in the semiarid interior of South Africa appears to be a result of quasigeostrophically triggered convection. Investigation of TTT formation in the context of planetary waves shows that early development is sometimes associated with previous anticyclonic wave breaking south of the subcontinent, with full maturity of TTTs occurring as a potential vorticity trough approaches the continent from the west. Sensitivity to upstream wave perturbations and effects on anticyclonic wave breaking in the South Indian Ocean are also observed.

The Valanginian δ13C excursion may not be an expression of a global oceanic anoxic event

Marine and terrestrial sediments of the Valanginian age display a distinct positive δ13C excursion, which has recently been interpreted as the expression of an oceanic anoxic episode (OAE) of global importance. Here we evaluate the extent of anaerobic conditions in marine bottom waters and explore the mechanisms involved in changing carbon storage on a global scale during this time interval. We asses redox-sensitive trace-element distributions (RSTE; uranium, vanadium, cobalt, arsenic and molybdenum) and the quality and quantity of preserved organic matter (OM) in representative sections along a shelf-basin transect in the western Tethys, in the Polish Basin and on Shatsky Rise. OM-rich layers corresponding in time to the δ13C shift are generally rare in the Tethyan sections and if present, the layers are not thicker than several centimetres and total organic carbon (TOC) contents do not surpass 1.68 wt..%. Palynological observations and geochemical properties of the preserved OM suggest a mixed marine and terrestrial origin and deposition in an oxic environment. In the Polish Basin, OM-rich layers show evidence for an important continental component. RSTE exhibit no major enrichments during the δ13C excursion in all studied Tethyan sections. RSTE enrichments are, however, observed in the pre-δ13C excursion OM-rich “Barrande” levels of the Vocontian Trough. In addition, all Tethyan sections record higher Mn contents during the δ13C shift, indicating rather well-oxygenated bottom waters in the western Tethys and the presence of anoxic basins elsewhere, such as the restricted basins of the North Atlantic and Weddell Sea. We propose that the Valanginian δ13C shift is the consequence of a combination of increased OM storage in marginal seas and on continents (as a sink of 12C-enriched organic carbon), coupled with the demise of shallow-water carbonate platforms (diminishing the storage capacity of 13C-enriched carbonate carbon). As such the Valanginian provides a more faithful natural analogue to present-day environmental change than most other Mesozoic OAEs, which are characterized by the development of ocean-wide dysaerobic to anaerobic conditions.

Modulation of visually evoked cortical fMRI responses by phase of ongoing occipital alpha oscillations

Using simultaneous electroencephalography as a measure of ongoing activity and functional magnetic resonance imaging (fMRI) as a measure of the stimulus-driven neural response, we examined whether the amplitude and phase of occipital alpha oscillations at the onset of a brief visual stimulus affects the amplitude of the visually evoked fMRI response. When accounting for intrinsic coupling of alpha amplitude and occipital fMRI signal by modeling and subtracting pseudo-trials, no significant effect of prestimulus alpha amplitude on the evoked fMRI response could be demonstrated. Regarding the effect of alpha phase, we found that stimuli arriving at the peak of the alpha cycle yielded a lower blood oxygenation level-dependent (BOLD) fMRI response in early visual cortex (V1/V2) than stimuli presented at the trough of the cycle. Our results therefore show that phase of occipital alpha oscillations impacts the overall strength of a visually evoked response, as indexed by the BOLD signal. This observation complements existing evidence that alpha oscillations reflect periodic variations in cortical excitability and suggests that the phase of oscillations in postsynaptic potentials can serve as a mechanism of gain control for incoming neural activity. Finally, our findings provide a putative neural basis for observations of alpha phase dependence of visual perceptual performance.

Impact of global warming on the Asian winter monsoon in a coupled GCM

The Asian winter monsoon (AWM) response to the global warming was investigated through a long-term integration of the transient greenhouse warming with the ECHAM4/OPYC3 CGCM. The physics of the response was studied through analyses of the impact of the global warming on the variations of the ocean and land contrast near the ground in the Asian and western Pacific region and the east Asian trough and jet stream in the middle and upper troposphere. Forcing of transient eddy activity on the zonal circulation over the Asian and western Pacific region was also analyzed. It is found that in the global warming scenario the winter northeasterlies along the Pacific coast of the Eurasian continent weaken systematically and significantly, and intensity of the AWM reduces evidently, but the AWM variances on the interannual and interdecadal scales are not affected much by the global warming. It is suggested that the global warming makes the climate over the most part of Asia to be milder with enhanced moisture in winter. In the global warming scenario the contrasts of the sea level pressure and the near-surface temperature between the Asian continent and the Pacific Ocean become significantly smaller, northward and eastward shifts and weakening of the east Asian trough and jet stream in the middle and upper troposphere are found. As a consequence, the cold air in the AWM originating from the east Asian trough and high latitudes is less powerful. In addition, feedback of the transient activity also makes a considerable contribution to the higher-latitude shift of the jet stream over the North Pacific in the global warming scenario.

Synoptic scale disturbances of the Indian Summer Monsoon as simulated in a high resolution climate model

Hamburg atmospheric general circulation model ECHAM3 at T106 resolution (1.125' lat.Aon.) has considerable skill in reproducing the observed seasonal reversal of mean sea level pressure, the location of the summer heat low as well as the position of the monsoon trough over the Indian subcontinent. The present-day climate and its seasonal cycle are realistically simulated by the model over this region. The model simulates the structure, intensity, frequency, movement and lifetime of monsoon depressions remarkably well. The number of monsoon depressions/storms simulated by the model in a year ranged from 5 to 12 with an average frequency of 8.4 yr-', not significantly different from the observed climatology. The model also simulates the interannual variability in the formation of depressions over the north Bay of Bengal during the summer monsoon season. In the warmer atmosphere under doubled CO2 conditions, the number of monsoon depressions/cyclonic storms forming in Indian seas in a year ranged from 5 to 11 with an average frequency of 7.6 yr-', not significantly different from those inferred in the control run of the model. However, under doubled CO2 conditions, fewer depressions formed in the month of June. Neither the lowest central pressure nor the maximum wind speed changes appreciably in monsoon depressions identified under simulated enhanced greenhouse conditions. The analysis suggests there will be no significant changes in the number and intensity of monsoon depressions in a warmer atmosphere.

Quantifying the climatological relationship between extratropical cyclone intensity and atmospheric precursors

We introduce a novel technique in which linear regression analysis is applied to clusters of tracked cyclones to statistically assess the factors controlling cyclone development. We illustrate this technique by evaluating the differences between cyclones forming in the west and east North Atlantic (herein termed west and east Atlantic cyclones). Enhanced cyclone intensity 2 days after genesis is found to be associated with deeper upper-level troughs upstream of the cyclone center at the genesis time in both west and east Atlantic cyclones. However, whilst west Atlantic cyclones are also enhanced by the presence of strong fronts, east Atlantic cyclones are not. Instead, east Atlantic cyclones exhibit an enhancement when diabatically generated midlevel potential vorticity is present (with the enhancement being of approximately equal magnitude to that associated with the potential vorticity in the upper-level trough). This is consistent with the paradigm of latent heat release in the warm conveyor belt region playing an important role in the development of east Atlantic cyclones.

Runge-Kutta IMEX schemes for the Horizontally Explicit/Vertically Implicit (HEVI) solution of wave equations

Many operational weather forecasting centres use semi-implicit time-stepping schemes because of their good efficiency. However, as computers become ever more parallel, horizontally explicit solutions of the equations of atmospheric motion might become an attractive alternative due to the additional inter-processor communication of implicit methods. Implicit and explicit (IMEX) time-stepping schemes have long been combined in models of the atmosphere using semi-implicit, split-explicit or HEVI splitting. However, most studies of the accuracy and stability of IMEX schemes have been limited to the parabolic case of advection–diffusion equations. We demonstrate how a number of Runge–Kutta IMEX schemes can be used to solve hyperbolic wave equations either semi-implicitly or HEVI. A new form of HEVI splitting is proposed, UfPreb, which dramatically improves accuracy and stability of simulations of gravity waves in stratified flow. As a consequence it is found that there are HEVI schemes that do not lose accuracy in comparison to semi-implicit ones. The stability limits of a number of variations of trapezoidal implicit and some Runge–Kutta IMEX schemes are found and the schemes are tested on two vertical slice cases using the compressible Boussinesq equations split into various combinations of implicit and explicit terms. Some of the Runge–Kutta schemes are found to be beneficial over trapezoidal, especially since they damp high frequencies without dropping to first-order accuracy. We test schemes that are not formally accurate for stiff systems but in stiff limits (nearly incompressible) and find that they can perform well. The scheme ARK2(2,3,2) performs the best in the tests.

20th century intraseasonal Asian monsoon dynamics viewed from Isomap

The Asian summer monsoon is a high dimensional and highly nonlinear phenomenon involving considerable moisture transport towards land from the ocean, and is critical for the whole region. We have used daily ECMWF reanalysis (ERA-40) sea-level pressure (SLP) anomalies to the seasonal cycle, over the region 50-145°E, 20°S-35°N to study the nonlinearity of the Asian monsoon using Isomap. We have focused on the two-dimensional embedding of the SLP anomalies for ease of interpretation. Unlike the unimodality obtained from tests performed in empirical orthogonal function space, the probability density function, within the two-dimensional Isomap space, turns out to be bimodal. But a clustering procedure applied to the SLP data reveals support for three clusters, which are identified using a three-component bivariate Gaussian mixture model. The modes are found to appear similar to active and break phases of the monsoon over South Asia in addition to a third phase, which shows active conditions over the Western North Pacific. Using the low-level wind field anomalies the active phase over South Asia is found to be characterised by a strengthening and an eastward extension of the Somali jet whereas during the break phase the Somali jet is weakened near southern India, while the monsoon trough in northern India also weakens. Interpretation is aided using the APHRODITE gridded land precipitation product for monsoon Asia. The effect of large-scale seasonal mean monsoon and lower boundary forcing, in the form of ENSO, is also investigated and discussed. The outcome here is that ENSO is shown to perturb the intraseasonal regimes, in agreement with conceptual ideas.

Modelling the effects of gravity waves on stratocumulus clouds observed during VOCALS-UK

During the VOCALS campaign spaceborne satellite observations showed that travelling gravity wave packets, generated by geostrophic adjustment, resulted in perturbations to marine boundary layer (MBL) clouds over the south-east Pacific Ocean (SEP). Often, these perturbations were reversible in that passage of the wave resulted in the clouds becoming brighter (in the wave crest), then darker (in the wave trough) and subsequently recovering their properties after the passage of the wave. However, occasionally the wave packets triggered irreversible changes to the clouds, which transformed from closed mesoscale cellular convection to open form. In this paper we use large eddy simulation (LES) to examine the physical mechanisms that cause this transition. Specifically, we examine whether the clearing of the cloud is due to (i) the wave causing additional cloud-top entrainment of warm, dry air or (ii) whether the additional condensation of liquid water onto the existing drops and the subsequent formation of drizzle are the important mechanisms. We find that, although the wave does cause additional drizzle formation, this is not the reason for the persistent clearing of the cloud; rather it is the additional entrainment of warm, dry air into the cloud followed by a reduction in longwave cooling, although this only has a significant effect when the cloud is starting to decouple from the boundary layer. The result in this case is a change from a stratocumulus to a more patchy cloud regime. For the simulations presented here, cloud condensation nuclei (CCN) scavenging did not play an important role in the clearing of the cloud. The results have implications for understanding transitions between the different cellular regimes in marine boundary layer (MBL) clouds.

Using synoptic type analysis to understand New Zealand climate during the Mid-Holocene

Diagnosing the climate of New Zealand from low-resolution General Circulation Models (GCMs) is notoriously difficult due to the interaction of the complex topography and the Southern Hemisphere (SH) mid-latitude westerly winds. Therefore, methods of downscaling synoptic scale model data for New Zealand are useful to help understand past climate. New Zealand also has a wealth of palaeoclimate-proxy data to which the downscaled model output can be compared, and to provide a qualitative method of assessing the capability of GCMs to represent, in this case, the climate 6000 yr ago in the Mid-Holocene. In this paper, a synoptic weather and climate regime classification system using Empirical Orthogonal Function (EOF) analysis of GCM and reanalysis data was used. The climate regimes are associated with surface air temperature and precipitation anomalies over New Zealand. From the analysis in this study, we find at 6000 BP that increased trough activity in summer and autumn led to increased precipitation, with an increased north-south pressure gradient ("zonal events") in winter and spring leading to drier conditions. Opposing effects of increased (decreased) temperature are also seen in spring (autumn) in the South Island, which are associated with the increased zonal (trough) events; however, the circulation induced changes in temperature are likely to have been of secondary importance to the insolation induced changes. Evidence from the palaeoclimate-proxy data suggests that the Mid-Holocene was characterized by increased westerly wind events in New Zealand, which agrees with the preference for trough and zonal regimes in the models.

Development of convection along the SPCZ within a Madden-Julian oscillation

A subtropical Rossby-wave propagation mechanism is proposed to account for the poleward and eastward progression of intraseasonal convective anomalies along the South Pacific convergence zone (SPCZ) that is observed in a significant proportion of Madden–Julian oscillations (MJOs). Large-scale convection, associated with an MJO, is assumed to be already established over the Indonesian region. The latent heating associated with this convection forces an equatorial Rossby-wave response with an upper-tropospheric anticyclone centred over, or slightly to the west of, the convection. Large potential-vorticity (PV) gradients, associated with the subtropical jet and the tropopause, lie just poleward of the anticyclone, and large magnitude PV air is advected equatorwards on the eastern side of the anticyclone. This ‘high’ PV air, or upper-tropospheric trough, is far enough off the equator that it has associated strong horizontal temperature gradients, and it induces deep ascent on its eastern side, at a latitude of about 15–30°. If this deep ascent is over a region susceptible to deep convection, such as the SPCZ, then convection may be forced or triggered. Hence convection develops along the SPCZ as a forced response to convection over Indonesia. The response mechanism is essentially one of subtropical Rossby-wave propagation. This hypothesis is based on a case study of a particularly strong MJO in early 1988, and is tested by idealized modelling studies. The mechanism may also be relevant to the existence of the mean SPCZ, as a forced response to mean Indonesian convection.