Diagnosis of variability and trends in a global precipitation dataset using a physically motivated statistical model

A physically motivated statistical model is used to diagnose variability and trends in wintertime ( October - March) Global Precipitation Climatology Project (GPCP) pentad (5-day mean) precipitation. Quasi-geostrophic theory suggests that extratropical precipitation amounts should depend multiplicatively on the pressure gradient, saturation specific humidity, and the meridional temperature gradient. This physical insight has been used to guide the development of a suitable statistical model for precipitation using a mixture of generalized linear models: a logistic model for the binary occurrence of precipitation and a Gamma distribution model for the wet day precipitation amount. The statistical model allows for the investigation of the role of each factor in determining variations and long-term trends. Saturation specific humidity q(s) has a generally negative effect on global precipitation occurrence and with the tropical wet pentad precipitation amount, but has a positive relationship with the pentad precipitation amount at mid- and high latitudes. The North Atlantic Oscillation, a proxy for the meridional temperature gradient, is also found to have a statistically significant positive effect on precipitation over much of the Atlantic region. Residual time trends in wet pentad precipitation are extremely sensitive to the choice of the wet pentad threshold because of increasing trends in low-amplitude precipitation pentads; too low a choice of threshold can lead to a spurious decreasing trend in wet pentad precipitation amounts. However, for not too small thresholds, it is found that the meridional temperature gradient is an important factor for explaining part of the long-term trend in Atlantic precipitation.

Oestrogenic activity of p-hydroxybenzoic acid (common metabolite of paraben esters) and methylparaben in human breast cancer cell lines

This paper addresses the question of whether p-hydroxybenzoic acid, the common metabolite of parabens, possesses oestrogenic activity in human breast cancer cell lines. The alkyl esters of p-hydroxybenzoic acid (parabens) are used widely as preservatives in consumer products to which the human population is exposed and have been shown previously to possess oestrogenic activity and to be present in human breast tumour tissue, which is an oestrogen-responsive tissue. Recent work has shown p-hydroxybenzoic acid to give an oestrogenic response in the rodent uterotrophic assay. We report here that p-hydroxybenzoic acid possesses oestrogenic activity in a panel of assays in human breast cancer cell lines. p-Hydroxybenzoic acid was able to displace [H-3]oestradiol from cytosolic oestrogen receptor of MCF7 human breast cancer cells by 54% at 5 x 10(6)-fold molar excess and by 99% at 10(7)-fold molar excess. It was able to increase the expression of a stably integrated oestrogen responsive reporter gene (ERE-CAT) at a concentration of 5 x 10(-4) M in MCF7 cells after 24 h and 7 days, which could be inhibited by the anti-oestrogen ICI 182 780 (Faslodex, fulvestrant). Proliferation of two human breast cancer cell lines (MCF7, ZR-75-1) could be increased by 10(-5) M p-hydroxybenzoic acid. Following on from previous studies showing a decrease in oestrogenic activity of parabens with shortening of the linear alkyl chain length, this study has compared the oestrogenic activity of p-hydroxybenzoic acid where the alkyl grouping is no longer present with methylparaben, which has the shortest alkyl group. Intrinsic oestrogenic activity of p-hydroxybenzoic acid was similar to that of methylparaben in terms of relative binding to the oestrogen receptor but its oestrogenic activity on gene expression and cell proliferation was lower than that of methylparaben. It can be concluded that removal of the ester group from parabens does not abrogate its oestrogenic activity and that p-hydroxybenzoic acid can give oestrogenic responses in human breast cancer cells. Copyright (C) 2005 John Wiley & Sons, Ltd.

Intra-individual variation of serum prostate specific antigen levels in men with benign prostate biopsies

To determine the intra-individual (physiological) variation of prostate-specific antigen (PSA) measurements in men after a benign prostatic biopsy. Sixty-four men were prospectively assessed, all of whom had a benign prostatic biopsy within the preceding 13 months. The degree of intra-individual variability was established by calculating the coefficient of variation on four PSA levels obtained from each patient weekly over a month. Six patients were subsequently diagnosed with prostate cancer and their data are presented separately. In the remaining 58 patients the median (range) individual mean PSA value was 6.3 (0.5-34.1) ng/mL. The median (range) coefficient of variation within the group was 9.5 (2.4-76.1)%. There was a clear linear relationship between mean PSA level and the standard deviation. In 48 of the 63 patients analysed, the coefficient of variation for serum PSA values in the group as a whole was greater than the variation claimed for the assay technique. The significance of the linear relationship between PSA and the standard deviation is discussed, with particular reference to those men who had a benign prostate biopsy.

A case of four-coordinate silver(I) adopting a high energy planar geometry. Experiment and theory

The ligands PhL and MeL are obtained by condensing 2-formylpyridine with benzil dihydrazone and diacetyl dihydrazone, respectively, in 2: 1 molar proportion. With silver( I), PhL yields a double-stranded dinuclear cationic helicate 1 in which the metal is tetrahedral but MeL gives a cationic one-dimensional polymeric complex 2 where silver( I) is distorted square planar and the ligand backbone is nearly planar. In both complexes, metal: ligand ratio is 1: 1. Ab initio calculations on the ligands at the HF/6-31+G* level reveal that while PhL strongly prefers a helical conformation, MeL has a natural inclination to remain in a planar conformation. Density functional theory calculations on model silver( I) complexes show that formation of the linear polymer in the case of MeL is also an important factor in imposing the planar geometry of Ag(I) in 2.

The threading of [(CuCN)(2)(mu-4,4 '-bpy)] sheets by CuCN chains

Yellow (CuCN)(2)[(CuCN)(2)(mu-4,4'-bpy)], formed in the hydrothermal reaction of CuCN with 4,4'-bipyridine at 453 K, contains two types of infinite CuCN chains. One set of CuCN chains is linked by 4,4'-bpy ligands to form almost flat sheets of composition [(CuCN)(2)(mu-4,4'-bpy)]. Holes in these sheets are aligned to allow pairs of approximately linear, infinite -(CuCN)- chains to thread through them. The closest interatomic approach between copper atoms in the threading chains and host sheets (similar to2.74 Angstrom) does not appear to represent a significant covalent bond as it leads to only a small distortion of the -(CuCN)- chains from linearity The relationship of this material to the previously determined structures of the host [(CuCN)(2)(mu-4,4'-bpy)] sheets and (CuCN)(3)[(CuCN)(2)(mu-4,4'-bPY)](2), in which these sheets are threaded by single -(CuCN)- chains, is discussed.

Simulating the detailed chemical composition of secondary organic aerosol formed on a regional scale during the TORCH 2003 campaign in the southern UK

Following on from the companion study (Johnson et al., 2006), a photochemical trajectory model (PTM) has been used to simulate the chemical composition of organic aerosol for selected events during the 2003 TORCH (Tropospheric Organic Chemistry Experiment) field campaign. The PTM incorporates the speciated emissions of 124 nonmethane anthropogenic volatile organic compounds (VOC) and three representative biogenic VOC, a highly-detailed representation of the atmospheric degradation of these VOC, the emission of primary organic aerosol (POA) material and the formation of secondary organic aerosol (SOA) material. SOA formation was represented by the transfer of semi and non-volatile oxidation products from the gas-phase to a condensed organic aerosol-phase, according to estimated thermodynamic equilibrium phase-partitioning characteristics for around 2000 reaction products. After significantly scaling all phase-partitioning coefficients, and assuming a persistent background organic aerosol (both required in order to match the observed organic aerosol loadings), the detailed chemical composition of the simulated SOA has been investigated in terms of intermediate oxygenated species in the Master Chemical Mechanism, version 3.1 ( MCM v3.1). For the various case studies considered, 90% of the simulated SOA mass comprises between ca. 70 and 100 multifunctional oxygenated species derived, in varying amounts, from the photooxidation of VOC of anthropogenic and biogenic origin. The anthropogenic contribution is dominated by aromatic hydrocarbons and the biogenic contribution by alpha-and beta-pinene (which also constitute surrogates for other emitted monoterpene species). Sensitivity in the simulated mass of SOA to changes in the emission rates of anthropogenic and biogenic VOC has also been investigated for 11 case study events, and the results have been compared to the detailed chemical composition data. The role of accretion chemistry in SOA formation, and its implications for the results of the present investigation, is discussed.

Mathematical modelling of competitive LDL/VLDL binding and uptake by hepatocytes

Elevated levels of low-density-lipoprotein cholesterol (LDL-C) in the plasma are a well-established risk factor for the development of coronary heart disease. Plasma LDL-C levels are in part determined by the rate at which LDL particles are removed from the bloodstream by hepatic uptake. The uptake of LDL by mammalian liver cells occurs mainly via receptor-mediated endocytosis, a process which entails the binding of these particles to specific receptors in specialised areas of the cell surface, the subsequent internalization of the receptor-lipoprotein complex, and ultimately the degradation and release of the ingested lipoproteins' constituent parts. We formulate a mathematical model to study the binding and internalization (endocytosis) of LDL and VLDL particles by hepatocytes in culture. The system of ordinary differential equations, which includes a cholesterol-dependent pit production term representing feedback regulation of surface receptors in response to intracellular cholesterol levels, is analysed using numerical simulations and steady-state analysis. Our numerical results show good agreement with in vitro experimental data describing LDL uptake by cultured hepatocytes following delivery of a single bolus of lipoprotein. Our model is adapted in order to reflect the in vivo situation, in which lipoproteins are continuously delivered to the hepatocyte. In this case, our model suggests that the competition between the LDL and VLDL particles for binding to the pits on the cell surface affects the intracellular cholesterol concentration. In particular, we predict that when there is continuous delivery of low levels of lipoproteins to the cell surface, more VLDL than LDL occupies the pit, since VLDL are better competitors for receptor binding. VLDL have a cholesterol content comparable to LDL particles; however, due to the larger size of VLDL, one pit-bound VLDL particle blocks binding of several LDLs, and there is a resultant drop in the intracellular cholesterol level. When there is continuous delivery of lipoprotein at high levels to the hepatocytes, VLDL particles still out-compete LDL particles for receptor binding, and consequently more VLDL than LDL particles occupy the pit. Although the maximum intracellular cholesterol level is similar for high and low levels of lipoprotein delivery, the maximum is reached more rapidly when the lipoprotein delivery rates are high. The implications of these results for the design of in vitro experiments is discussed.

Parallel Hybrid Monte Carlo Algorithms for Matrix Computations

In this paper we consider hybrid (fast stochastic approximation and deterministic refinement) algorithms for Matrix Inversion (MI) and Solving Systems of Linear Equations (SLAE). Monte Carlo methods are used for the stochastic approximation, since it is known that they are very efficient in finding a quick rough approximation of the element or a row of the inverse matrix or finding a component of the solution vector. We show how the stochastic approximation of the MI can be combined with a deterministic refinement procedure to obtain MI with the required precision and further solve the SLAE using MI. We employ a splitting A = D – C of a given non-singular matrix A, where D is a diagonal dominant matrix and matrix C is a diagonal matrix. In our algorithm for solving SLAE and MI different choices of D can be considered in order to control the norm of matrix T = D –1C, of the resulting SLAE and to minimize the number of the Markov Chains required to reach given precision. Further we run the algorithms on a mini-Grid and investigate their efficiency depending on the granularity. Corresponding experimental results are presented.

Model selection approaches for non-linear system identification: a review

The identification of non-linear systems using only observed finite datasets has become a mature research area over the last two decades. A class of linear-in-the-parameter models with universal approximation capabilities have been intensively studied and widely used due to the availability of many linear-learning algorithms and their inherent convergence conditions. This article presents a systematic overview of basic research on model selection approaches for linear-in-the-parameter models. One of the fundamental problems in non-linear system identification is to find the minimal model with the best model generalisation performance from observational data only. The important concepts in achieving good model generalisation used in various non-linear system-identification algorithms are first reviewed, including Bayesian parameter regularisation and models selective criteria based on the cross validation and experimental design. A significant advance in machine learning has been the development of the support vector machine as a means for identifying kernel models based on the structural risk minimisation principle. The developments on the convex optimisation-based model construction algorithms including the support vector regression algorithms are outlined. Input selection algorithms and on-line system identification algorithms are also included in this review. Finally, some industrial applications of non-linear models are discussed.

An efficient flux difference splitting algorithm for unsteady duct flows of a real gas

A numerical scheme is presented for the solution of the Euler equations of compressible flow of a real gas in a single spatial coordinate. This includes flow in a duct of variable cross-section, as well as flow with slab, cylindrical or spherical symmetry, as well as the case of an ideal gas, and can be useful when testing codes for the two-dimensional equations governing compressible flow of a real gas. The resulting scheme requires an average of the flow variables across the interface between cells, and this average is chosen to be the arithmetic mean for computational efficiency, which is in contrast to the usual “square root” averages found in this type of scheme. The scheme is applied with success to five problems with either slab or cylindrical symmetry and for a number of equations of state. The results compare favourably with the results from other schemes.

An efficient numerical method for compressible flows of a real gas using arithmetic averaging

An efficient numerical method is presented for the solution of the Euler equations governing the compressible flow of a real gas. The scheme is based on the approximate solution of a specially constructed set of linearised Riemann problems. An average of the flow variables across the interface between cells is required, and this is chosen to be the arithmetic mean for computational efficiency, which is in contrast to the usual square root averaging. The scheme is applied to a test problem for five different equations of state.

An efficient numerical scheme for the Euler equations

A finite difference scheme based on flux difference splitting is presented for the solution of the Euler equations for the compressible flow of an ideal gas. A linearised Riemann problem is defined, and a scheme based on numerical characteristic decomposition is presented for obtaining approximate solutions to the linearised problem. An average of the flow variables across the interface between cells is required, and this average is chosen to be the arithmetic mean for computational efficiency, leading to arithmetic averaging. This is in contrast to the usual ‘square root’ averages found in this type of Riemann solver, where the computational expense can be prohibitive. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids nonphysical, spurious oscillations. The scheme is applied to a shock tube problem and a blast wave problem. Each approximate solution compares well with those given by other schemes, and for the shock tube problem is in agreement with the exact solution.

An efficient shock capturing algorithm for compressible flows in a duct of variable cross-section

A numerical scheme is presented for the solution of the Euler equations of compressible flow of a gas in a single spatial co-ordinate. This includes flow in a duct of variable cross-section as well as flow with slab, cylindrical or spherical symmetry and can prove useful when testing codes for the two-dimensional equations governing compressible flow of a gas. The resulting scheme requires an average of the flow variables across the interface between cells and for computational efficiency this average is chosen to be the arithmetic mean, which is in contrast to the usual ‘square root’ averages found in this type of scheme. The scheme is applied with success to five problems with either slab or cylindrical symmetry and a comparison is made in the cylindrical case with results from a two-dimensional problem with no sources.

Uniformization of affine dependence programs for parallel embedded system design

This paper is concerned with the uniformization of a system of afine recurrence equations. This transformation is used in the design (or compilation) of highly parallel embedded systems (VLSI systolic arrays, signal processing filters, etc.). In this paper, we present and implement an automatic system to achieve uniformization of systems of afine recurrence equations. We unify the results from many earlier papers, develop some theoretical extensions, and then propose effective uniformization algorithms. Our results can be used in any high level synthesis tool based on polyhedral representation of nested loop computations.

Nocturnal ranging behaviour of urban hedgehogs, Erinaceus europaeus, in relation to risk and reward

Urban areas have both positive and negative influences on wildlife. For terrestrial mammals, one of the principle problems is the risk associated with moving through the environment whilst foraging. In this study, we examined nocturnal patterns of movement of urban-dwelling hedgehogs (Erinaceus europaeus) in relation to (i) the risks posed by predators and motor vehicles and (ii) nightly weather patterns. Hedgehogs preferentially utilised the gardens of semi-detached and terraced houses. However, females, but not males, avoided the larger back gardens of detached houses, which contain more of the habitat features selected by badgers. This difference in the avoidance of predation risk is probably associated with sex differences in breeding behaviour. Differences in nightly movement patterns were consistent with strategies associated with mating behaviour and the accumulation of fat reserves for hibernation. Hedgehogs also exhibited differences in behaviour associated with the risks posed by humans; they avoided actively foraging near roads and road verges, but did not avoid crossing roads per se. They were, however, significantly more active after midnight when there was a marked reduction in vehicle and foot traffic. In particular, responses to increased temperature, which is associated with increased abundance of invertebrate prey, were only observed after midnight. This variation in the timing of bouts of activity would reduce the risks associated with human activities. There were also profound differences in both area ranged and activity with chronological year which warrant further investigation.