Estimation of optimal gravity wave parameters for climate models using data assimilation

There is a current need to constrain the parameters of gravity wave drag (GWD) schemes in climate models using observational information instead of tuning them subjectively. In this work, an inverse technique is developed using data assimilation principles to estimate gravity wave parameters. Because mostGWDschemes assume instantaneous vertical propagation of gravity waves within a column, observations in a single column can be used to formulate a one-dimensional assimilation problem to estimate the unknown parameters. We define a cost function that measures the differences between the unresolved drag inferred from observations (referred to here as the ‘observed’ GWD) and the GWD calculated with a parametrisation scheme. The geometry of the cost function presents some difficulties, including multiple minima and ill-conditioning because of the non-independence of the gravity wave parameters. To overcome these difficulties we propose a genetic algorithm to minimize the cost function, which provides a robust parameter estimation over a broad range of prescribed ‘true’ parameters. When real experiments using an independent estimate of the ‘observed’ GWD are performed, physically unrealistic values of the parameters can result due to the non-independence of the parameters. However, by constraining one of the parameters to lie within a physically realistic range, this degeneracy is broken and the other parameters are also found to lie within physically realistic ranges. This argues for the essential physical self-consistency of the gravity wave scheme. A much better fit to the observed GWD at high latitudes is obtained when the parameters are allowed to vary with latitude. However, a close fit can be obtained either in the upper or the lower part of the profiles, but not in both at the same time. This result is a consequence of assuming an isotropic launch spectrum. The changes of sign in theGWDfound in the tropical lower stratosphere, which are associated with part of the quasi-biennial oscillation forcing, cannot be captured by the parametrisation with optimal parameters.

A refined sampling strategy for intra-tooth stable isotope analysis of mammalian enamel

Serial sampling and stable isotope analysis performed along the growth axis of vertebrate tooth enamel records differences attributed to seasonal variation in diet, climate or animal movement. Because several months are required to obtain mature enamel in large mammals, modifications in the isotopic composition of environmental parameters are not instantaneously recorded, and stable isotope analysis of tooth enamel returns a time-averaged signal attenuated in its amplitude relative to the input signal. For convenience, stable isotope profiles are usually determined on the side of the tooth where enamel is thickest. Here we investigate the possibility of improving the time resolution by targeting the side of the tooth where enamel is thinnest. Observation of developing third molars (M3) in sheep shows that the tooth growth rate is not constant but decreases exponentially, while the angle between the first layer of enamel deposited and the enamel–dentine junction increases as a tooth approaches its maximal length. We also noted differences in thickness and geometry of enamel growth between the mesial side (i.e., the side facing the M2) and the buccal side (i.e., the side facing the cheek) of the M3. Carbon and oxygen isotope variations were measured along the M3 teeth from eight sheep raised under controlled conditions. Intra-tooth variability was systematically larger along the mesial side and the difference in amplitude between the two sides was proportional to the time of exposure to the input signal. Although attenuated, the mesial side records variations in the environmental signal more faithfully than the buccal side. This approach can be adapted to other mammals whose teeth show lateral variation in enamel thickness and could potentially be used as an internal check for diagenesis.

Heliospheric Observations of STEREO-Directed Coronal Mass Ejections in 2008 – 2010: Lessons for Future Observations of Earth-Directed CMEs

We present a study of coronal mass ejections (CMEs) which impacted one of the STEREO spacecraft between January 2008 and early 2010. We focus our study on 20 CMEs which were observed remotely by the Heliospheric Imagers (HIs) onboard the other STEREO spacecraft up to large heliocentric distances. We compare the predictions of the Fixed-Φ and Harmonic Mean (HM) fitting methods, which only differ by the assumed geometry of the CME. It is possible to use these techniques to determine from remote-sensing observations the CME direction of propagation, arrival time and final speed which are compared to in-situ measurements. We find evidence that for large viewing angles, the HM fitting method predicts the CME direction better. However, this may be due to the fact that only wide CMEs can be successfully observed when the CME propagates more than 100∘ from the observing spacecraft. Overall eight CMEs, originating from behind the limb as seen by one of the STEREO spacecraft can be tracked and their arrival time at the other STEREO spacecraft can be successfully predicted. This includes CMEs, such as the events on 4 December 2009 and 9 April 2010, which were viewed 130∘ away from their direction of propagation. Therefore, we predict that some Earth-directed CMEs will be observed by the HIs until early 2013, when the separation between Earth and one of the STEREO spacecraft will be similar to the separation of the two STEREO spacecraft in 2009 – 2010.

Multiwavelength Observations of EXO 0748-676. I. Reprocessing of X-Ray Bursts

We present high time-resolution multiwavelength observations of X-ray bursts in the low-mass X-ray binary UY Vol. Strong reprocessed signals are present in the ultraviolet and optical, lagged and smeared with respect to the X-rays. The addition of far-ultraviolet coverage for one burst allows much tighter constraints on the temperature and geometry of the reprocessing region than previously possible. A blackbody reprocessing model for this burst suggests a rise in temperatures during the burst from 18,000 to 35,000 K and an emitting area comparable to that expected for the disk and/or irradiated companion star. The lags are consistent with those expected. The single-zone blackbody model cannot reproduce the ratio of optical to ultraviolet flux during the burst, however. The discrepancy seems too large to explain with deviations from a local blackbody spectrum and more likely indicates that a range of reprocessing temperatures are required. Comparable results are derived from other bursts, and in particular the lag and smearing both appear shorter when the companion star is on the near side of the disk as predicted. The burst observed by HST also yielded a spectrum of the reprocessed light. It is dominated by continuum, with a spectral shape consistent with the temperatures derived from lightcurve modeling. Taken as a whole, our observations confirm the standard paradigm of prompt reprocessing distributed across the disk and companion star, with the response dominated by a thermalized continuum rather than by emission lines.

Optimal planning of cocoa clonal selection programmes

Data from three cocoa (Theobroma cacao) clonal selection trials are used to investigate the genetic and environmental components of variation in yield and the percentage of total pods affected by black pod disease (Phytophtora pod rot). Simulations based on these estimated components of variation are then used to discuss the best choice in future of numbers of clones, replicates and years of harvest to maximise selection advances in the traits measured. The three main conclusions are the need to increase the number of clones at the expense of the number of replicates of each clone, the diminishing returns from additional years of harvesting and the importance of widening the genetic base of the clones chosen to be tested.

Simplified estimation of ray-path mirroring height for HF radiowaves reflected from the ionospheric F-region

From Milsom's equations, which describe the geometry of ray-path hops reflected from the ionospheric F-layer, algorithms for the simplified estimation of mirror-reflection height are developed. These allow for hop length and the effects of variations in underlying ionisation (via the ratio of the F2- and E-layer critical frequencies) and F2-layer peak height (via the M(3000)F2-factor). Separate algorithms are presented which are applicable to a range of signal frequencies about the FOT and to propagation at the MUF. The accuracies and complexities of the algorithms are compared with those inherent in the use of a procedure based on an equation developed by Shimazaki.

Portable smartphone quantitation of prostate specific antigen (PSA) in a fluoropolymer microfluidic device

We present a new, power-free and flexible detection system named MCFphone for portable colorimetric and fluorescence quantitative sandwich immunoassay detection of prostate specific antigen (PSA). The MCFphone is composed by a smartphone integrated with a magnifying lens, a simple light source and a miniaturised immunoassay platform, the Microcapillary Film (MCF). The excellent transparency and flat geometry of fluoropolymer MCF allowed quantitation of PSA in the range 0.9 to 60 ng/ml with < 7 % precision in 13 minutes using enzymatic amplification and a chromogenic substrate. The lower limit of detection was further improved from 0.4 to 0.08 ng/ml in whole blood samples with the use of a fluorescence substrate. The MCFphone has shown capable of performing rapid (13 to 22 minutes total assay time) colorimetric quantitative and highly sensitive fluorescence tests with good %Recovery, which represents a major step in the integration of a new generation of inexpensive and portable microfluidic devices with commercial immunoassay reagents and off-the-shelf smartphone technology.

Numbers in motion

We explore the debates surrounding the constructive and discursive capabilities of accounting information focusing in particular on the reception volatility of numbers once they are produced and ‘exposed’ to various communities of minds. Drawing on Goffman’s (1974) frame analysis and Vollmer’s (2007) work on the three-dimensional character of numerical signs, we explore how numbers can go through gradual or instantaneous transformations, get caught up in public debates and become ‘agents’ or ‘captives’ in creating social order and in some cases social drama. In our analysis we also relate to the work of Durkheim (1993, 2002) on the sociology of morality to illustrate how numbers can become indicators of moral transgression. The study explores both historical and contemporary examples of controversies and recent accounting scandals to demonstrate how preparers (of financial information) can lose control over numbers which then acquire new meanings through social context and collective (re)framing. The main contribution of the study is to illustrate how the narratives attached to numbers are malleable and fluid across both time and space.

Low rank representation on Riemannian manifold of symmetric positive definite matrices

Sparse coding aims to find a more compact representation based on a set of dictionary atoms. A well-known technique looking at 2D sparsity is the low rank representation (LRR). However, in many computer vision applications, data often originate from a manifold, which is equipped with some Riemannian geometry. In this case, the existing LRR becomes inappropriate for modeling and incorporating the intrinsic geometry of the manifold that is potentially important and critical to applications. In this paper, we generalize the LRR over the Euclidean space to the LRR model over a specific Rimannian manifold—the manifold of symmetric positive matrices (SPD). Experiments on several computer vision datasets showcase its noise robustness and superior performance on classification and segmentation compared with state-of-the-art approaches.

Monitoring one-electron photo-oxidation of guanine in DNA crystals using ultrafast infrared spectroscopy

To understand the molecular origins of diseases caused by ultraviolet and visible light, and also to develop photodynamic therapy, it is important to resolve the mechanism of photoinduced DNA damage. Damage to DNA bound to a photosensitizer molecule frequently proceeds by one-electron photo-oxidation of guanine, but the precise dynamics of this process are sensitive to the location and the orientation of the photosensitizer, which are very difficult to define in solution. To overcome this, ultrafast time-resolved infrared (TRIR) spectroscopy was performed on photoexcited ruthenium polypyridyl–DNA crystals, the atomic structure of which was determined by X-ray crystallography. By combining the X-ray and TRIR data we are able to define both the geometry of the reaction site and the rates of individual steps in a reversible photoinduced electron-transfer process. This allows us to propose an individual guanine as the reaction site and, intriguingly, reveals that the dynamics in the crystal state are quite similar to those observed in the solvent medium.

Environmental control for layers

Light patterns have less effect on numbers of eggs laid by current stocks than on those of forty years ago, but the principles have not changed. Ovarian activity is stimulated by increasing photoperiods and suppressed by decreasing photoperiods. The light pattern used during rearing can still have large effects on age at 50% lay, even for modern stocks. Early sexual maturity maximises egg numbers but gives smaller eggs. Late maturity maximises egg size at the expense of numbers. The relationship between egg output (g/hen d) and age at first egg is curvilinear, with maximum yield occurring in flocks maturing in about the centre of their potential range. Fancy patterns of increasing daylength after maturity are probably not justified. A flock held on a constant 14h day will lay as many eggs as one given step up lighting. Intermittent lighting saves about 5% of feed consumption with no loss of output, provided that the feed has adequate amino acid content to allow for the reduced feed intake. Producers with light-proof laying houses should be taking advantage of intermittent lighting. The recommended light intensity for laying houses is still 10 lx, although the physiological threshold for response to changes in photoperiod is closer to 2 lx. Very dim (0.05 lx) light filtering into blacked out houses will not stimulate the hypothalamic receptors responsible for photo-sexual responses, but may affect the bird's biological clock, which can alter its response to a constant short photoperiod. Feed intake shows a curvilinear dependence on environmental temperature. At temperatures below the panting threshold, performance can be maintained by adjusting the feed so as to maintain an adequate intake of critical amino acids. Above the panting threshold, the hen is unable to take in enough energy to maintain normal output. There is no dietary modification which can effectively offset this problem. Diurnally cycling temperatures result in feed intake and egg production equivalent to that observed under a constant temperature equal to the mean of the cycle. When the poultry house is cooler at night than by day, it helps to provide light so that the birds can feed during the cooler part of the cycle.

Nanospheres on a silver plate

A Viewpoint on: 'Surface Geometry of C60 on Ag(111)' H. I. Li, K. Pussi, K. J. Hanna, L.-L. Wang, D. D. Johnson, H.-P. Cheng, H. Shin, S. Curtarolo, W. Moritz, J. A. Smerdon, R. McGrath, and R. D. Diehl. . Published in Physical Review Letters 103, 056101 (2009) on July 27, 2009.

Integrating Hamiltonian systems defined on the Lie groups SO(4) and SO(1,3)

This paper examines optimal solutions of control systems with drift defined on the orthonormal frame bundle of particular Riemannian manifolds of constant curvature. The manifolds considered here are the space forms Euclidean space E³, the spheres S³ and the hyperboloids H³ with the corresponding frame bundles equal to the Euclidean group of motions SE(3), the rotation group SO(4) and the Lorentz group SO(1,3). The optimal controls of these systems are solved explicitly in terms of elliptic functions. In this paper, a geometric interpretation of the extremal solutions is given with particular emphasis to a singularity in the explicit solutions. Using a reduced form of the Casimir functions the geometry of these solutions is illustrated.

A helmet mounted display system with active gaze control for visual telepresence

The objective of a Visual Telepresence System is to provide the operator with a high fidelity image from a remote stereo camera pair linked to a pan/tilt device such that the operator may reorient the camera position by use of head movement. Systems such as these which utilise virtual reality style helmet mounted displays have a number of limitations. The geometry of the camera positions and of the displays is generally fixed and is most suitable only for viewing elements of a scene at a particular distance. To address such limitations, a prototype system has been developed where the geometry of the displays and cameras is dynamically controlled by the eye movement of the operator. This paper explores why it is necessary to actively adjust the display system as well as the cameras and justifies the use of mechanical adjustment of the displays as an alternative to adjustment by electronic or image processing methods. The electronic and mechanical design is described including optical arrangements and control algorithms. The performance and accuracy of the system is assessed with respect to eye movement.

Dynamic vergence telepresence system incorporating active gaze control

Visual Telepresence system which utilize virtual reality style helmet mounted displays have a number of limitations. The geometry of the camera positions and of the display is fixed and is most suitable only for viewing elements of a scene at a particular distance. In such a system, the operator's ability to gaze around without use of head movement is severely limited. A trade off must be made between a poor viewing resolution or a narrow width of viewing field. To address these limitations a prototype system where the geometry of the displays and cameras is dynamically controlled by the eye movement of the operator has been developed. This paper explores the reasons why is necessary to actively adjust both the display system and the cameras and furthermore justifies the use of mechanical adjustment of the displays as an alternative to adjustment by electronic or image processing methods. The electronic and mechanical design is described including optical arrangements and control algorithms, An assessment of the performance of the system against a fixed camera/display system when operators are assigned basic tasks involving depth and distance/size perception. The sensitivity to variations in transient performance of the display and camera vergence is also assessed.