Coupling aerosol surface and bulk chemistry with a kinetic double layer model (K2-SUB): oxidation of oleic acid by ozone

We present a kinetic double layer model coupling aerosol surface and bulk chemistry (K2-SUB) based on the PRA framework of gas-particle interactions (Poschl-Rudich-Ammann, 2007). K2-SUB is applied to a popular model system of atmospheric heterogeneous chemistry: the interaction of ozone with oleic acid. We show that our modelling approach allows de-convoluting surface and bulk processes, which has been a controversial topic and remains an important challenge for the understanding and description of atmospheric aerosol transformation. In particular, we demonstrate how a detailed treatment of adsorption and reaction at the surface can be coupled to a description of bulk reaction and transport that is consistent with traditional resistor model formulations. From literature data we have derived a consistent set of kinetic parameters that characterise mass transport and chemical reaction of ozone at the surface and in the bulk of oleic acid droplets. Due to the wide range of rate coefficients reported from different experimental studies, the exact proportions between surface and bulk reaction rates remain uncertain. Nevertheless, the model results suggest an important role of chemical reaction in the bulk and an approximate upper limit of similar to 10(-11) cm(2) s(-1) for the surface reaction rate coefficient. Sensitivity studies show that the surface accommodation coefficient of the gas-phase reactant has a strong non-linear influence on both surface and bulk chemical reactions. We suggest that K2-SUB may be used to design, interpret and analyse future experiments for better discrimination between surface and bulk processes in the oleic acid-ozone system as well as in other heterogeneous reaction systems of atmospheric relevance.

Non-ergodicity of inviscid two-dimensional flow on a beta-plane and on the surface of a rotating sphere

It is shown that, for a sufficiently large value of β, two-dimensional flow on a doubly-periodic beta-plane cannot be ergodic (phase-space filling) on the phase-space surface of constant energy and enstrophy. A corresponding result holds for flow on the surface of a rotating sphere, for a sufficiently rapid rotation rate Ω. This implies that the higher-order, non-quadratic invariants are exerting a significant influence on the statistical evolution of the flow. The proof relies on the existence of a finite-amplitude Liapunov stability theorem for zonally symmetric basic states with a non-vanishing absolute-vorticity gradient. When the domain size is much larger than the size of a typical eddy, then a sufficient condition for non-ergodicity is that the wave steepness ε < 1, where ε = 2[surd radical]2Z/βU in the planar case and $\epsilon = 2^{\frac{1}{4}} a^{\frac{5}{2}}Z^{\frac{7}{4}}/\Omega U^{\frac{5}{2}}$ in the spherical case, and where Z is the enstrophy, U the r.m.s. velocity, and a the radius of the sphere. This result may help to explain why numerical simulations of unforced beta-plane turbulence (in which ε decreases in time) seem to evolve into a non-ergodic regime at large scales.

Retrieval characteristics of non-linear sea surface temperature from the Advanced Very High Resolution Radiometer

Criteria are proposed for evaluating sea surface temperature (SST) retrieved from satellite infra-red imagery: bias should be small on regional scales; sensitivity to atmospheric humidity should be small; and sensitivity of retrieved SST to surface temperature should be close to 1 K K−1. Their application is illustrated for non-linear sea surface temperature (NLSST) estimates. 233929 observations from the Advanced Very High Resolution Radiometer (AVHRR) on Metop-A are matched with in situ data and numerical weather prediction (NWP) fields. NLSST coefficients derived from these matches have regional biases from −0.5 to +0.3 K. Using radiative transfer modelling we find that a 10% increase in humidity alone can change the retrieved NLSST by between −0.5 K and +0.1 K. A 1 K increase in SST changes NLSST by <0.5 K in extreme cases. The validity of estimates of sensitivity by radiative transfer modelling is confirmed empirically.

Estimating the direct aerosol radiative perturbation: impact of ocean surface representation and aerosol non-sphericity

Atmospheric aerosols are now actively studied, in particular because of their radiative and climate impacts. Estimations of the direct aerosol radiative perturbation, caused by extinction of incident solar radiation, usually rely on radiative transfer codes and involve simplifying hypotheses. This paper addresses two approximations which are widely used for the sake of simplicity and limiting the computational cost of the calculations. Firstly, it is shown that using a Lambertian albedo instead of the more rigorous bidirectional reflectance distribution function (BRDF) to model the ocean surface radiative properties leads to large relative errors in the instantaneous aerosol radiative perturbation. When averaging over the day, these errors cancel out to acceptable levels of less than 3% (except in the northern hemisphere winter). The other scope of this study is to address aerosol non-sphericity effects. Comparing an experimental phase function with an equivalent Mie-calculated phase function, we found acceptable relative errors if the aerosol radiative perturbation calculated for a given optical thickness is daily averaged. However, retrieval of the optical thickness of non-spherical aerosols assuming spherical particles can lead to significant errors. This is due to significant differences between the spherical and non-spherical phase functions. Discrepancies in aerosol radiative perturbation between the spherical and non-spherical cases are sometimes reduced and sometimes enhanced if the aerosol optical thickness for the spherical case is adjusted to fit the simulated radiance of the non-spherical case.

Intercomparison of methods of coupling between convection and large-scale circulation. 2: Comparison over non-uniform surface conditions

As part of an international intercomparison project, the weak temperature gradient (WTG) and damped gravity wave (DGW) methods are used to parameterize large-scale dynamics in a set of cloud-resolving models (CRMs) and single column models (SCMs). The WTG or DGW method is implemented using a configuration that couples a model to a reference state defined with profiles obtained from the same model in radiative-convective equilibrium. We investigated the sensitivity of each model to changes in SST, given a fixed reference state. We performed a systematic comparison of the WTG and DGW methods in different models, and a systematic comparison of the behavior of those models using the WTG method and the DGW method. The sensitivity to the SST depends on both the large-scale parameterization method and the choice of the cloud model. In general, SCMs display a wider range of behaviors than CRMs. All CRMs using either the WTG or DGW method show an increase of precipitation with SST, while SCMs show sensitivities which are not always monotonic. CRMs using either the WTG or DGW method show a similar relationship between mean precipitation rate and column-relative humidity, while SCMs exhibit a much wider range of behaviors. DGW simulations produce large-scale velocity profiles which are smoother and less top-heavy compared to those produced by the WTG simulations. These large-scale parameterization methods provide a useful tool to identify the impact of parameterization differences on model behavior in the presence of two-way feedback between convection and the large-scale circulation.

Does reactive surface area depend on grain size? Results from pH 3, 25 circle C far-from-equilibrium flow-through dissolution experiments on anorthite and biotite

Laboratory determined mineral weathering rates need to be normalised to allow their extrapolation to natural systems. The principle normalisation terms used in the literature are mass, and geometric- and BET specific surface area (SSA). The purpose of this study was to determine how dissolution rates normalised to these terms vary with grain size. Different size fractions of anorthite and biotite ranging from 180-150 to 20-10 mu m were dissolved in pH 3, HCl at 25 degrees C in flow through reactors under far from equilibrium conditions. Steady state dissolution rates after 5376 h (anorthite) and 4992 h (biotite) were calculated from Si concentrations and were normalised to initial- and final- mass and geometric-, geometric edge- (biotite), and BET SSA. For anorthite, rates normalised to initial- and final-BET SSA ranged from 0.33 to 2.77 X 10(-10) mol(feldspar) m(-2) s(-1), rates normalised to initial- and final-geometric SSA ranged from 5.74 to 8.88 X 10(-10) mol(feldspar) m(-2) s(-1) and rates normalised to initial- and final-mass ranged from 0.11 to 1.65 mol(feldspar) g(-1) s(-1). For biotite, rates normalised to initial- and final-BET SSA ranged from 1.02 to 2.03 X 10(-12) mol(biotite) m(-2) s(-1), rates normalised to initial- and final-geometric SSA ranged from 3.26 to 16.21 X 10(-12) mol(biotite) m(-2) s(-1), rates normalised to initial- and final-geometric edge SSA ranged from 59.46 to 111.32 x 10(-12) mol(biotite) m(-2) s(-1) and rates normalised to initial- and final-mass ranged from 0.81 to 6.93 X 10(-12) mol(biotite) g(-1) s(-1). For all normalising terms rates varied significantly (p <= 0.05) with grain size. The normalising terms which gave least variation in dissolution rate between grain sizes for anorthite were initial BET SSA and initial- and final-geometric SSA. This is consistent with: (1) dissolution being dominated by the slower dissolving but area dominant non-etched surfaces of the grains and, (2) the walls of etch pits and other dissolution features being relatively unreactive. These steady state normalised dissolution rates are likely to be constant with time. Normalisation to final BET SSA did not give constant ratios across grain size due to a non-uniform distribution of dissolution features. After dissolution coarser grains had a greater density of dissolution features with BET-measurable but unreactive wall surface area than the finer grains. The normalising term which gave the least variation in dissolution rates between grain sizes for biotite was initial BET SSA. Initial- and final-geometric edge SSA and final BET SSA gave the next least varied rates. The basal surfaces dissolved sufficiently rapidly to influence bulk dissolution rate and prevent geometric edge SSA normalised dissolution rates showing the least variation. Simple modelling indicated that biotite grain edges dissolved 71-132 times faster than basal surfaces. In this experiment, initial BET SSA best integrated the different areas and reactivities of the edge and basal surfaces of biotite. Steady state dissolution rates are likely to vary with time as dissolution alters the ratio of edge to basal surface area. Therefore they would be more properly termed pseudo-steady state rates, only appearing constant because the time period over which they were measured (1512 h) was less than the time period over wich they would change significantly. (c) 2006 Elsevier Inc. All rights reserved.

The composition and crystallinity of the near-surface regions of weathered alkali feldspars

Our ability to identify thin non-stoichiometric and amorphous layers beneath mineral surfaces has been tested by undertaking X-ray photoelectron spectroscopy (XPS) and transmission electron microscopy (TEM) work on alkali feldspars from pH 1 dissolution experiments. The outcomes of this work were used to help interpret XPS and TEM results from alkali feldspars weathered for <10,000 years in soils overlying the Shap Granite (north-west England). The chemistry of effluent solutions indicates that silica-rich layers a few nanometers in thickness formed during the pH I experiments. These layers can be successfully identified by XPS and have lower Al/Si, Na/Si, K/Si and Ca/Si values than the outermost similar to 9 nm of unweathered controls. Development of Al-Si non-stoichiometry is coupled with loss of crystal structure to produce amorphous layers that are identifiable by TEM where >similar to 2.5 nm thick, whereas the crystallinity of albite is retained despite leaching of Na to depths of tens to hundreds on nanometers. Integration of XPS data over the outermost 6-9 nm of naturally weathered Shap feldspars shows that they have stoichiometric Al/Si and K/Si ratios, which is consistent with findings of previous TEM work on the same material that they lack amorphous layers. There is some XPS evidence for loss of K from the outermost couple of nanometers of Shap orthoclase, and the possibility of leaching of Na from albite to greater depths cannot be excluded using the XPS or TEM results. This study demonstrates that the leached layer model, as formulated from laboratory experiments, is inapplicable to the weathering of alkali feldspars within acidic soils, which is an essentially stoichiometric reaction. (C) 2008 Elsevier Ltd. All rights reserved.

On the initiation of surface waves by turbulent shear flow

An analytical model is developed for the initial stage of surface wave generation at an air-water interface by a turbulent shear flow in either the air or in the water. The model treats the problem of wave growth departing from a flat interface and is relevant for small waves whose forcing is dominated by turbulent pressure fluctuations. The wave growth is predicted using the linearised and inviscid equations of motion, essentially following Phillips [Phillips, O.M., 1957. On the generation of waves by turbulent wind. J. Fluid Mech. 2, 417-445], but the pressure fluctuations that generate the waves are treated as unsteady and related to the turbulent velocity field using the rapid-distortion treatment of Durbin [Durbin, P.A., 1978. Rapid distortion theory of turbulent flows. PhD thesis, University of Cambridge]. This model, which assumes a constant mean shear rate F, can be viewed as the simplest representation of an oceanic or atmospheric boundary layer. For turbulent flows in the air and in the water producing pressure fluctuations of similar magnitude, the waves generated by turbulence in the water are found to be considerably steeper than those generated by turbulence in the air. For resonant waves, this is shown to be due to the shorter decorrelation time of turbulent pressure in the air (estimated as proportional to 1/Gamma), because of the higher shear rate existing in the air flow, and due to the smaller length scale of the turbulence in the water. Non-resonant waves generated by turbulence in the water, although being somewhat gentler, are still steeper than resonant waves generated by turbulence in the air. Hence, it is suggested that turbulence in the water may have a more important role than previously thought in the initiation of the surface waves that are subsequently amplified by feedback instability mechanisms.

Multidecadal modulation of El Nino-Southern Oscillation (ENSO) variance by Atlantic Ocean sea surface temperatures

Observations suggest a possible link between the Atlantic Multidecadal Oscillation (AMO) and El Nino Southern Oscillation (ENSO) variability, with the warm AMO phase being related to weaker ENSO variability. A coupled ocean-atmosphere model is used to investigate this relationship and to elucidate mechanisms responsible for it. Anomalous sea surface temperatures (SSTs) associated with the positive AMO lead to change in the basic state in the tropical Pacific Ocean. This basic state change is associated with a deepened thermocline and reduced vertical stratification of the equatorial Pacific ocean, which in turn leads to weakened ENSO variability. We suggest a role for an atmospheric bridge that rapidly conveys the influence of the Atlantic Ocean to the tropical Pacific. The results suggest a non-local mechanism for changes in ENSO statistics and imply that anomalous Atlantic ocean SSTs can modulate both mean climate and climate variability over the Pacific.

Surface and interstitial Ti diffusion at the rutile TiO2(110) surface

Diffusion of Ti through the TiO2 (110) rutile surface plays a key role in the growth and reactivity of TiO2. To understand the fundamental aspects of this important process, we present an analysis of the diffusion of Ti adspecies at the stoichiometric TiO2(110) surface using complementary computational methodologies of density functional theory corrected for on-site Coulomb interactions (DFT+U) and a charge equilibration (QEq) atomistic potential to identify minimum energy pathways. We find that diffusion of Ti from the surface to subsurface (and vice versa) follows an intersticialcy exchange mechanism, involving exchange of surface Ti with the 6-fold coordinated Ti below the bridging oxygen rows. Diffusion in the subsurface between layers also follows an interstitialcy mechanism. The diffusion of Ti is discussed in light of continued attempts to understand the re-oxidation of non-stoichiometric TiO2(110) surfaces.

Effects of fluctuating daily surface fluxes on the time-mean oceanic circulation

The effect of fluctuating daily surface fluxes on the time-mean oceanic circulation is studied using an empirical flux model. The model produces fluctuating fluxes resulting from atmospheric variability and includes oceanic feedbacks on the fluxes. Numerical experiments were carried out by driving an ocean general circulation model with three different versions of the empirical model. It is found that fluctuating daily fluxes lead to an increase in the meridional overturning circulation (MOC) of the Atlantic of about 1 Sv and a decrease in the Antarctic circumpolar current (ACC) of about 32 Sv. The changes are approximately 7% of the MOC and 16% of the ACC obtained without fluctuating daily fluxes. The fluctuating fluxes change the intensity and the depth of vertical mixing. This, in turn, changes the density field and thus the circulation. Fluctuating buoyancy fluxes change the vertical mixing in a non-linear way: they tend to increase the convective mixing in mostly stable regions and to decrease the convective mixing in mostly unstable regions. The ACC changes are related to the enhanced mixing in the subtropical and the mid-latitude Southern Ocean and reduced mixing in the high-latitude Southern Ocean. The enhanced mixing is related to an increase in the frequency and the depth of convective events. As these events bring more dense water downward, the mixing changes lead to a reduction in meridional gradient of the depth-integrated density in the Southern Ocean and hence the strength of the ACC. The MOC changes are related to more subtle density changes. It is found that the vertical mixing in a latitudinal strip in the northern North Atlantic is more strongly enhanced due to fluctuating fluxes than the mixing in a latitudinal strip in the South Atlantic. This leads to an increase in the density difference between the two strips, which can be responsible for the increase in the Atlantic MOC.

Effect of growth time on the surface and adhesion properties of Lactobacillus rhamnosus GG

Aims: To investigate the changes in the surface properties of Lactobacillus rhamnosus GG during growth, and relate them with the ability of the Lactobacillus cells to adhere to Caco-2 cells. Methods and Results: Lactobacillus rhamnosus GG was grown in complex medium, and cell samples taken at four time points and freeze dried. Untreated and trypsin treated freeze dried samples were analysed for their composition using SDS-PAGE analysis and Fourier transform infrared spectroscopy (FTIR), hydrophobicity and zeta potential, and for their ability to adhere to Caco-2 cells. The results suggested that in the case of early exponential phase samples (4 and 8 h), the net surface properties, i.e. hydrophobicity and charge, were determined to a large extent by anionic hydrophilic components, whereas in the case of stationary phase samples (13 and 26 h), hydrophobic proteins seemed to play the biggest role. Considerable differences were also observed between the ability of the different samples to adhere to Caco-2 cells; maximum adhesion was observed for the early stationary phase sample (13 h). The results suggested that the adhesion to Caco-2 cells was influenced by both proteins and non-proteinaceous compounds present on the surface of the Lactobacillus cells. Conclusion: The surface properties of Lact. rhamnosus GG changed during growth, which in return affected the ability of the Lactobacillus cells to adhere to Caco-2 cells. Significance and Impact of the Study: The levels of adhesion of Lactobacillus cells to Caco-2 cells were influenced by the growth time and reflected changes on the bacterial surface. This study provides critical information on the physicochemical factors that influence bacterial adhesion to intestinal cells.

Non-hemagglutinating flaviviruses: molecular mechanisms for the emergence of new strains via adaptation to European ticks

Tick-borne encephalitis virus (TBEV) causes human epidemics across Eurasia. Clinical manifestations range from inapparent infections and fevers to fatal encephalitis but the factors that determine disease severity are currently undefined. TBEV is characteristically a hemagglutinating (HA) virus; the ability to agglutinate erythrocytes tentatively reflects virion receptor/fusion activity. However, for the past few years many atypical HA-deficient strains have been isolated from patients and also from the natural European host tick, Ixodes persulcatus. By analysing the sequences of HA-deficient strains we have identified 3 unique amino acid substitutions (D67G, E122G or D277A) in the envelope protein, each of which increases the net charge and hydrophobicity of the virion surface. Therefore, we genetically engineered virus mutants each containing one of these 3 substitutions; they all exhibited HA-deficiency. Unexpectedly, each genetically modified non-HA virus demonstrated increased TBEV reproduction in feeding Ixodes ricinus, not the recognised tick host for these strains. Moreover, virus transmission efficiency between infected and uninfected ticks co-feeding on mice was also intensified by each substitution. Retrospectively, the mutation D67G was identified in viruses isolated from patients with encephalitis. We propose that the emergence of atypical Siberian HA-deficient TBEV strains in Europe is linked to their molecular adaptation to local ticks. This process appears to be driven by the selection of single mutations that change the virion surface thus enhancing receptor/fusion function essential for TBEV entry into the unfamiliar tick species. As the consequence of this adaptive mutagenesis, some of these mutations also appear to enhance the ability of TBEV to cross the human blood-brain barrier, a likely explanation for fatal encephalitis. Future research will reveal if these emerging Siberian TBEV strains continue to disperse westwards across Europe by adaptation to the indigenous tick species and if they are associated with severe forms of TBE.

Computer modelling and estimation of recruitment patterns of non-branching coral colonies at three sites in Wakatobi Marine Park, S.E. Sulawesi, Indonesia; implications for coral reef conservation

We have studied growth and estimated recruitment of massive coral colonies at three sites, Kaledupa, Hoga and Sampela, separated by about 1.5 km in the Wakatobi Marine National Park, S.E. Sulawesi, Indonesia. There was significantly higher species richness (P<0.05), coral cover (P<0.05) and rugosity (P<0.01) at Kaledupa than at Sampela. A model for coral reef growth has been developed based on a rational polynomial function, where dx/dt is an index of coral growth with time; W is the variable (for example, coral weight, coral length or coral area), up to the power of n in the numerator and m in the denominator; a1……an and b1…bm are constants. The values for n and m represent the degree of the polynomial, and can relate to the morphology of the coral. The model was used to simulate typical coral growth curves, and tested using published data obtained by weighing coral colonies underwater in reefs on the south-west coast of Curaçao [‘Neth. J. Sea Res. 10 (1976) 285’]. The model proved an accurate fit to the data, and parameters were obtained for a number of coral species. Surface area data was obtained on over 1200 massive corals at three different sites in the Wakatobi Marine National Park, S.E. Sulawesi, Indonesia. The year of an individual's recruitment was calculated from knowledge of the growth rate modified by application of the rational polynomial model. The estimated pattern of recruitment was variable, with little numbers of massive corals settling and growing before 1950 at the heavily used site, Sampela, relative to the reef site with little or no human use, Kaledupa, and the intermediate site, Hoga. There was a significantly greater sedimentation rate at Sampela than at either Kaledupa (P<0.0001) or Hoga (P<0.0005). The relative mean abundance of fish families present at the reef crests at the three sites, determined using digital video photography, did not correlate with sedimentation rates, underwater visibility or lack of large non-branching coral colonies. Radial growth rates of three genera of non-branching corals were significantly lower at Sampela than at Kaledupa or at Hoga, and there was a high correlation (r=0.89) between radial growth rates and underwater visibility. Porites spp. was the most abundant coral over all the sites and at all depths followed by Favites (P<0.04) and Favia spp. (P<0.03). Colony ages of Porites corals were significantly lower at the 5 m reef flat on the Sampela reef than at the same depth on both other reefs (P<0.005). At Sampela, only 2.8% of corals on the 5 m reef crest are of a size to have survived from before 1950. The Scleractinian coral community of Sampela is severely impacted by depositing sediments which can lead to the suffocation of corals, whilst also decreasing light penetration resulting in decreased growth and calcification rates. The net loss of material from Sampela, if not checked, could result in the loss of this protective barrier which would be to the detriment of the sublittoral sand flats and hence the Sampela village.

Development of a method fon non-destructive testing of fruits using scanning laser vibrometry (SLV)

Quality control on fruits requires reliable methods, able to assess with reasonable accuracy and possibly in a non-destructive way their physical and chemical characteristics. More specifically, a decreased firmness indicates the presence of damage or defects in the fruit or else that the fruit has exceeded its “best before date”, becoming unsuitable for consumption. In high-value exotic fruits, such as mangoes, where firmness cannot be easily measured from a simple observation of texture, colour changes and unevenness of fruits surface, the use of non-destructive techniques is highly recommendable. In particular, the application of Laser vibrometry, based on the Doppler effect, a non-contact technique sensitive to differences in displacements inferior to the nanometre, appears ideal for a possible on-line control on food. Previous results indicated that a phase shift can be in a repeatable way associated with the presence of damage on the fruit, whilst a decreased firmness results in significant differences in the displacement of the fruits under the same excitation signal. In this work, frequency ranges for quality control via the application of a sound chirp are suggested, based on the measurement of the signal coherence. The variations of the average vibration spectrum of a grid of points, or of point-by-point signal velocity allows the go-no go recognition of “firm” and “over-ripe” fruits, with notable success in the particular case of mangoes. The future exploitation of this work will include the application of this method to allow on-line control during conveyor belt distribution of fruits.