A simple column model to explore anticipated problems in variational assimilation of satellite observations

We investigate a simplified form of variational data assimilation in a fully nonlinear framework with the aim of extracting dynamical development information from a sequence of observations over time. Information on the vertical wind profile, w(z ), and profiles of temperature, T (z , t), and total water content, qt (z , t), as functions of height, z , and time, t, are converted to brightness temperatures at a single horizontal location by defining a two-dimensional (vertical and time) variational assimilation testbed. The profiles of T and qt are updated using a vertical advection scheme. A basic cloud scheme is used to obtain the fractional cloud amount and, when combined with the temperature field, this information is converted into a brightness temperature, using a simple radiative transfer scheme. It is shown that our model exhibits realistic behaviour with regard to the prediction of cloud, but the effects of nonlinearity become non-negligible in the variational data assimilation algorithm. A careful analysis of the application of the data assimilation scheme to this nonlinear problem is presented, the salient difficulties are highlighted, and suggestions for further developments are discussed.

The evaluation of CloudSat and CALIPSO ice microphysical products using ground-based cloud radar and lidar observations

In this paper, the statistical properties of tropical ice clouds (ice water content, visible extinction, effective radius, and total number concentration) derived from 3 yr of ground-based radar–lidar retrievals from the U.S. Department of Energy Atmospheric Radiation Measurement Climate Research Facility in Darwin, Australia, are compared with the same properties derived using the official CloudSat microphysical retrieval methods and from a simpler statistical method using radar reflectivity and air temperature. It is shown that the two official CloudSat microphysical products (2B-CWC-RO and 2B-CWC-RVOD) are statistically virtually identical. The comparison with the ground-based radar–lidar retrievals shows that all satellite methods produce ice water contents and extinctions in a much narrower range than the ground-based method and overestimate the mean vertical profiles of microphysical parameters below 10-km height by over a factor of 2. Better agreements are obtained above 10-km height. Ways to improve these estimates are suggested in this study. Effective radii retrievals from the standard CloudSat algorithms are characterized by a large positive bias of 8–12 μm. A sensitivity test shows that in response to such a bias the cloud longwave forcing is increased from 44.6 to 46.9 W m−2 (implying an error of about 5%), whereas the negative cloud shortwave forcing is increased from −81.6 to −82.8 W m−2. Further analysis reveals that these modest effects (although not insignificant) can be much larger for optically thick clouds. The statistical method using CloudSat reflectivities and air temperature was found to produce inaccurate mean vertical profiles and probability distribution functions of effective radius. This study also shows that the retrieval of the total number concentration needs to be improved in the official CloudSat microphysical methods prior to a quantitative use for the characterization of tropical ice clouds. Finally, the statistical relationship used to produce ice water content from extinction and air temperature obtained by the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) satellite is evaluated for tropical ice clouds. It is suggested that the CALIPSO ice water content retrieval is robust for tropical ice clouds, but that the temperature dependence of the statistical relationship used should be slightly refined to better reproduce the radar–lidar retrievals.

The statistical properties of tropical ice clouds generated by the West-African and Australian monsoons from ground-based radar-lidar observations

In the present paper we characterize the statistical properties of non-precipitating tropical ice clouds (deep ice anvils resulting from deep convection and cirrus clouds) over Niamey, Niger, West Africa, and Darwin, northern Australia, using ground-based radar–lidar observations from the Atmospheric Radiation Measurement (ARM) programme. The ice cloud properties analysed in this paper are the frequency of ice cloud occurrence, cloud fraction, the morphological properties (cloud-top height, base height, and thickness), the microphysical and radiative properties (ice water content, visible extinction, effective radius, terminal fall speed, and concentration), and the internal cloud dynamics (in-cloud vertical air velocity). The main highlight of the paper is that it characterizes for the first time the probability density functions of the tropical ice cloud properties, their vertical variability and their diurnal variability at the same time. This is particularly important over West Africa, since the ARM deployment in Niamey provides the first vertically resolved observations of non-precipitating ice clouds in this crucial area in terms of redistribution of water and energy in the troposphere. The comparison between the two sites also provides an additional observational basis for the evaluation of the parametrization of clouds in large-scale models, which should be able to reproduce both the statistical properties at each site and the differences between the two sites. The frequency of ice cloud occurrence is found to be much larger over Darwin when compared to Niamey, and with a much larger diurnal variability, which is well correlated with the diurnal cycle of deep convective activity. The diurnal cycle of the ice cloud occurrence over Niamey is also much less correlated with that of deep convective activity than over Darwin, probably owing to the fact that Niamey is further away from the deep convective sources of the region. The frequency distributions of cloud fraction are strongly bimodal and broadly similar over the two sites, with a predominance of clouds characterized either by a very small cloud fraction (less than 0.3) or a very large cloud fraction (larger than 0.9). The ice clouds over Darwin are also much thicker (by 1 km or more statistically) and are characterized by a much larger diurnal variability than ice clouds over Niamey. Ice clouds over Niamey are also characterized by smaller particle sizes and fall speeds but in much larger concentrations, thereby carrying more ice water and producing more visible extinction than the ice clouds over Darwin. It is also found that there is a much larger occurrence of downward in-cloud air motions less than 1 m s−1 over Darwin, which together with the larger fall speeds retrieved over Darwin indicates that the life cycle of ice clouds is probably shorter over Darwin than over Niamey.

Two-photon excitation with pico-second fluorescence lifetime imaging to detect nuclear association of flavanols

Two-photon excitation enabled for the first time the observation and measurement of excited state fluorescence lifetimes from three flavanols in solution, which were ∼1.0 ns for catechin and epicatechin, but <45 ps for epigallocatechin gallate (EGCG). The shorter lifetime for EGCG is in line with a lower fluorescence quantum yield of 0.003 compared to catechin (0.015) and epicatechin (0.018). In vivo experiments with onion cells demonstrated that tryptophan and quercetin, which tend to be major contributors of background fluorescence in plant cells, have sufficiently low cross sections for two-photon excitation at 630 nm and therefore do not interfere with detection of externally added or endogenous flavanols in Allium cepa or Taxus baccata cells. Applying two-photon excitation to flavanols enabled 3-D fluorescence lifetime imaging microscopy and showed that added EGCG penetrated the whole nucleus of onion cells. Interestingly, EGCG and catechin showed different lifetime behaviour when bound to the nucleus: EGCG lifetime increased from <45 to 200 ps, whilst catechin lifetime decreased from 1.0 ns to 500 ps. Semi-quantitative measurements revealed that the relative ratios of EGCG concentrations in nucleoli associated vesicles: nucleus: cytoplasm were ca. 100:10:1. Solution experiments with catechin, epicatechin and histone proteins provided preliminary evidence, via the appearance of a second lifetime (τ2 = 1.9–3.1 ns), that both flavanols may be interacting with histone proteins. We conclude that there is significant nuclear absorption of flavanols. This advanced imaging using two-photon excitation and biophysical techniques described here will prove valuable for probing the intracellular trafficking and functions of flavanols, such as EGCG, which is the major flavanol of green tea.

Prediction and functional analysis of native disorder in proteins from the Three Kingdoms of Life

An automatic method for recognizing natively disordered regions from amino acid sequence is described and benchmarked against predictors that were assessed at the latest critical assessment of techniques for protein structure prediction (CASP) experiment. The method attains a Wilcoxon score of 90.0, which represents a statistically significant improvement on the methods evaluated on the same targets at CASP. The classifier, DISOPRED2, was used to estimate the frequency of native disorder in several representative genomes from the three kingdoms of life. Putative, long (>30 residue) disordered segments are found to occur in 2.0% of archaean, 4.2% of eubacterial and 33.0% of eukaryotic proteins. The function of proteins with long predicted regions of disorder was investigated using the gene ontology annotations supplied with the Saccharomyces genome database. The analysis of the yeast proteome suggests that proteins containing disorder are often located in the cell nucleus and are involved in the regulation of transcription and cell signalling. The results also indicate that native disorder is associated with the molecular functions of kinase activity and nucleic acid binding.

The effect of nonlinearity on the variational assimilation of satellite observations using a simple column model

Cloud imagery is not currently used in numerical weather prediction (NWP) to extract the type of dynamical information that experienced forecasters have extracted subjectively for many years. For example, rapidly developing mid-latitude cyclones have characteristic signatures in the cloud imagery that are most fully appreciated from a sequence of images rather than from a single image. The Met Office is currently developing a technique to extract dynamical development information from satellite imagery using their full incremental 4D-Var (four-dimensional variational data assimilation) system. We investigate a simplified form of this technique in a fully nonlinear framework. We convert information on the vertical wind field, w(z), and profiles of temperature, T(z, t), and total water content, qt (z, t), as functions of height, z, and time, t, to a single brightness temperature by defining a 2D (vertical and time) variational assimilation testbed. The profiles of w, T and qt are updated using a simple vertical advection scheme. We define a basic cloud scheme to obtain the fractional cloud amount and, when combined with the temperature field, we convert this information into a brightness temperature, having developed a simple radiative transfer scheme. With the exception of some matrix inversion routines, all our code is developed from scratch. Throughout the development process we test all aspects of our 2D assimilation system, and then run identical twin experiments to try and recover information on the vertical velocity, from a sequence of observations of brightness temperature. This thesis contains a comprehensive description of our nonlinear models and assimilation system, and the first experimental results.

On spectral invariance of single scattering albedo for water droplets and ice crystals at weakly absorbing wavelengths

The single scattering albedo w_0l in atmospheric radiative transfer is the ratio of the scattering coefficient to the extinction coefficient. For cloud water droplets both the scattering and absorption coefficients, thus the single scattering albedo, are functions of wavelength l and droplet size r. This note shows that for water droplets at weakly absorbing wavelengths, the ratio w_0l(r)/w_0l(r0) of two single scattering albedo spectra is a linear function of w_0l(r). The slope and intercept of the linear function are wavelength independent and sum to unity. This relationship allows for a representation of any single scattering albedo spectrum w_0l(r) via one known spectrum w_0l(r0). We provide a simple physical explanation of the discovered relationship. Similar linear relationships were found for the single scattering albedo spectra of non-spherical ice crystals.

Rhinocetin, a venom-derived integrin-specific antagonist inhibits collagen-induced platelet and endothelial cell functions

Snaclecs are small non-enzymatic proteins present in viper venoms reported to modulate haemostasis of victims through effects on platelets, vascular endothelial and smooth muscle cells. In this study, we have isolated and functionally characterised a snaclec which we named rhinocetin from the venom of West African gaboon viper, Bitis gabonica rhinoceros. Rhinocetin was shown to comprise α and β chains with the molecular masses of 13.5 and 13kDa respectively. Sequence and immunoblot analysis of rhinocetin confirmed this to be a novel snaclec. Rhinocetin inhibited collagen-stimulated activation of human platelets in dose dependent manner, but displayed no inhibitory effects on glycoprotein VI (collagen receptor) selective agonist, CRP-XL-, ADP- or thrombin-induced platelet activation. Rhinocetin antagonised the binding of monoclonal antibodies against the α2 subunit of integrin α2β1 to platelets and coimmunoprecipitation analysis confirmed integrin α2β1 as a target for this venom protein. Rhinocetin inhibited a range of collagen induced platelet functions such as fibrinogen binding, calcium mobilisation, granule secretion, aggregation and thrombus formation. It also inhibited integrin α2β1 dependent functions of human endothelial cells. Together, our data suggest rhinocetin to be a modulator of integrin α2β1 function and thus may provide valuable insights into the role of this integrin in physiological and pathophysiological scenarios including haemostasis, thrombosis and envenomation.

A neurofuzzy classifier for two class problems

A neurofuzzy classifier identification algorithm is introduced for two class problems. The initial fuzzy base construction is based on fuzzy clustering utilizing a Gaussian mixture model (GMM) and the analysis of covariance (ANOVA) decomposition. The expectation maximization (EM) algorithm is applied to determine the parameters of the fuzzy membership functions. Then neurofuzzy model is identified via the supervised subspace orthogonal least square (OLS) algorithm. Finally a logistic regression model is applied to produce the class probability. The effectiveness of the proposed neurofuzzy classifier has been demonstrated using a real data set.

Modulation of potassium ion channel proteins utilising antibodies

The application of antibodies to living cells has the potential to modulate the function of specific proteins by virtue of their high specificity. This specificity has proven effective in determining the involvement of many proteins in neuronal function where specific agonists and antagonists do not exist, e.g. ion channel subunits. We discuss a way to utilise subunit specific antibodies to target individual channel subunits in electrophysiological experiments to determine functional roles within native neurones. Utilising this approach, we have investigated the role of the voltage-gated potassium channel Kv3.1b subunit within a region of the brainstem important in the regulation of autonomic function. We provide some useful control experiments in order to help validate this method. We conclude that antibodies can be extremely valuable in determining the functions of specific proteins in living neurones in neuroscience research.

Immunopharmacology - antibodies for specific modulation of proteins involved in neuronal function

The application of antibodies to living neurones has the potential to modulate function of specific proteins by virtue of their high specificity. This specificity has proven effective in determining the involvement of many proteins in neuronal function where specific agonists and antagonists do not exist, e.g. ion channel subunits. We discuss studies where antibodies modulate functions of voltage gated sodium, voltage gated potassium, voltage gated calcium hyperpolarisation activated cyclic nucleotide (HCN gated) and transient receptor potential (TRP) channels. Ligand gated channels studied in this way include nicotinic acetylcholine receptors, purinoceptors and GABA receptors. Antibodies have also helped reveal the involvement of different intracellular proteins in neuronal functions including G-proteins as well as other proteins involved in trafficking, phosphoinositide signalling and neurotransmitter release. Some suggestions for control experiments are made to help validate the method. We conclude that antibodies can be extremely valuable in determining the functions of specific proteins in living neurones in neuroscience research.

Complex formation in solutions of oppositely charged polyelectrolytes at different polyion compositions and salt content

The formation of complexes in solutions containing positively charged polyions (polycations) and a variable amount of negatively charged polyions (polyanions) has been investigated by Monte Carlo simulations. The polyions were described as flexible chains of charged hard spheres interacting through a screened Coulomb potential. The systems were analyzed in terms of cluster compositions, structure factors, and radial distribution functions. At 50% charge equivalence or less, complexes involving two polycations and one polyanion were frequent, while closer to charge equivalence, larger clusters were formed. Small and neutral complexes dominated the solution at charge equivalence in a monodisperse system, while larger clusters again dominated the solution when the polyions were made polydisperse. The cluster composition and solution structure were also examined as functions of added salt by varying the electrostatic screening length. The observed formation of clusters could be rationalized by a few simple rules.

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

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

A unified neurofuzzy model for classification

This work proposes a unified neurofuzzy modelling scheme. To begin with, the initial fuzzy base construction method is based on fuzzy clustering utilising a Gaussian mixture model (GMM) combined with the analysis of covariance (ANOVA) decomposition in order to obtain more compact univariate and bivariate membership functions over the subspaces of the input features. The mean and covariance of the Gaussian membership functions are found by the expectation maximisation (EM) algorithm with the merit of revealing the underlying density distribution of system inputs. The resultant set of membership functions forms the basis of the generalised fuzzy model (GFM) inference engine. The model structure and parameters of this neurofuzzy model are identified via the supervised subspace orthogonal least square (OLS) learning. Finally, instead of providing deterministic class label as model output by convention, a logistic regression model is applied to present the classifier’s output, in which the sigmoid type of logistic transfer function scales the outputs of the neurofuzzy model to the class probability. Experimental validation results are presented to demonstrate the effectiveness of the proposed neurofuzzy modelling scheme.

Padé approximants for the acoustical characteristics of rigid frame porous media

In this paper it is shown that a number of theoretical models of the acoustical properties of rigid frame porous media, especially those involving ratios of Bessel functions of complex argument, can be accurately approximated and greatly simplified by the use of Padé approximation techniques. In the case of the model of Attenborough [J. Acoust. Soc. Am. 81, 93–102 (1987)] rational approximations are produced for the characteristic impedance, propagation constant, dynamic compressibility, and dynamic density, as a function of frequency and the material parameters. The model proposed by Stinson and Champoux