Detecting recombination in 4-taxa DNA sequence alignments with Bayesian hidden Markov models and Markov chain Monte Carlo

This article presents a statistical method for detecting recombination in DNA sequence alignments, which is based on combining two probabilistic graphical models: (1) a taxon graph (phylogenetic tree) representing the relationship between the taxa, and (2) a site graph (hidden Markov model) representing interactions between different sites in the DNA sequence alignments. We adopt a Bayesian approach and sample the parameters of the model from the posterior distribution with Markov chain Monte Carlo, using a Metropolis-Hastings and Gibbs-within-Gibbs scheme. The proposed method is tested on various synthetic and real-world DNA sequence alignments, and we compare its performance with the established detection methods RECPARS, PLATO, and TOPAL, as well as with two alternative parameter estimation schemes.

In situ analysis and structural elucidation of sainfoin (Onobrychis viciifolia) tannins for high-throughput germplasm screening

A rapid thiolytic degradation and cleanup procedure was developed for analyzing tannins directly in chlorophyll-containing sainfoin (Onobrychis viciifolia) plants. The technique proved suitable for complex tannin mixtures containing catechin, epicatechin, gallocatechin, and epigallocatechin flavan-3-ol units. The reaction time was standardized at 60 min to minimize the loss of structural information as a result of epimerization and degradation of terminal flavan-3-ol units. The results were evaluated by separate analysis of extractable and unextractable tannins, which accounted for 63.6−113.7% of the in situ plant tannins. It is of note that 70% aqueous acetone extracted tannins with a lower mean degree of polymerization (mDP) than was found for tannins analyzed in situ. Extractable tannins had between 4 and 29 lower mDP values. The method was validated by comparing results from individual and mixed sample sets. The tannin composition of different sainfoin accessions covered a range of mDP values from 16 to 83, procyanidin/prodelphinidin (PC/PD) ratios from 19.2/80.8 to 45.6/54.4, and cis/trans ratios from 74.1/25.9 to 88.0/12.0. This is the first high-throughput screening method that is suitable for analyzing condensed tannin contents and structural composition directly in green plant tissue.

An Exponentially Convergent Nonpolynomial Finite Element Method for Time-Harmonic Scattering from Polygons

In recent years nonpolynomial finite element methods have received increasing attention for the efficient solution of wave problems. As with their close cousin the method of particular solutions, high efficiency comes from using solutions to the Helmholtz equation as basis functions. We present and analyze such a method for the scattering of two-dimensional scalar waves from a polygonal domain that achieves exponential convergence purely by increasing the number of basis functions in each element. Key ingredients are the use of basis functions that capture the singularities at corners and the representation of the scattered field towards infinity by a combination of fundamental solutions. The solution is obtained by minimizing a least-squares functional, which we discretize in such a way that a matrix least-squares problem is obtained. We give computable exponential bounds on the rate of convergence of the least-squares functional that are in very good agreement with the observed numerical convergence. Challenging numerical examples, including a nonconvex polygon with several corner singularities, and a cavity domain, are solved to around 10 digits of accuracy with a few seconds of CPU time. The examples are implemented concisely with MPSpack, a MATLAB toolbox for wave computations with nonpolynomial basis functions, developed by the authors. A code example is included.

Deconvolution method for removing the distortion from time jitter in equivalent-time waveform samplers

A sampling oscilloscope is one of the main units in automatic pulse measurement system (APMS). The time jitter in waveform samplers is an important error source that affect the precision of data acquisition. In this paper, this kind of error is greatly reduced by using the deconvolution method. First, the probability density function (PDF) of time jitter distribution is determined by the statistical approach, then, this PDF is used as convolution kern to deconvolve with the acquired waveform data with additional averaging, and the result is the waveform data in which the effect of time jitter has been removed, and the measurement precision of APMS is greatly improved. In addition, some computer simulations are given which prove the success of the method given in this paper.

Comprehensive ground-based and in situ observations of substorm expansion phase onset

In this paper, we present comprehensive ground-based and space-based in situ geosynchronous observations of a substorm expansion phase onset on 1 October 2005. The Double Star TC-2 and GOES-12 spacecraft were both located within the substorm current wedge during the substorm expansion phase onset, which occurred over the Canadian sector. We find that an onset of ULF waves in space was observed after onset on the ground by extending the AWESOME timing algorithm into space. Furthermore, a population of low-energy field-aligned electrons was detected by the TC-2 PEACE instrument contemporaneous with the ULF waves in space. These electrons appear to be associated with an enhancement of field-aligned Poynting flux into the ionosphere which is large enough to power visible auroral displays. The observations are most consistent with a near-Earth initiation of substorm expansion phase onset, such as the Near-Geosynchronous Onset (NGO) substorm scenario. A lack of data from further downtail, however, means other mechanisms cannot be ruled out.

Seasonality and extent of extratropical TST derived from in-situ CO measurements during SPURT

We present airborne in-situ trace gas measurements which were performed on eight campaigns between November 2001 and July 2003 during the SPURT-project (SPURenstofftransport in der Tropopausenregion, trace gas transport in the tropopause region). The measurements on a quasi regular basis allowed an overview of the seasonal variations of the trace gas distribution in the tropopause region over Europe from 35°-75°N to investigate the influence of transport and mixing across the extratropical tropopause on the lowermost stratosphere. From the correlation of CO and O3 irreversible mixing of tropospheric air into the lowermost stratosphere is identified. The CO distribution indicates that transport and subsequent mixing of tropospheric air across the extratropical tropopause predominantly affects a layer, which closely follows the shape of the local tropopause. In addition, the seasonal cycle of CO2 illustrates the strong coupling of that layer to the extratropical troposphere. Both, horizontal gradients of CO on isentropes as well as the CO-O3-distribution in the lowermost stratosphere reveal that the influence of quasi-horizontal transport and subsequent mixing weakens with distance from the local tropopause. The mixing layer extends to about 25 K in potential temperature above the local tropopause exhibiting only a weak seasonality. However, at large distances from the tropopause a significant influence of tropospheric air is still evident. The relation between N2O and CO2 indicates that a significant contribution of air originating from the tropical tropopause contributes to the background air in the extratropical lowermost stratosphere.

The accuracy of SST retrievals from AATSR: An initial assessment through geophysical validation against in situ radiometers, buoys and other SST data sets

The Advanced Along-Track Scanning Radiometer (AATSR) was launched on Envisat in March 2002. The AATSR instrument is designed to retrieve precise and accurate global sea surface temperature (SST) that, combined with the large data set collected from its predecessors, ATSR and ATSR-2, will provide a long term record of SST data that is greater than 15 years. This record can be used for independent monitoring and detection of climate change. The AATSR validation programme has successfully completed its initial phase. The programme involves validation of the AATSR derived SST values using in situ radiometers, in situ buoys and global SST fields from other data sets. The results of the initial programme presented here will demonstrate that the AATSR instrument is currently close to meeting its scientific objectives of determining global SST to an accuracy of 0.3 K (one sigma). For night time data, the analysis gives a warm bias of between +0.04 K (0.28 K) for buoys to +0.06 K (0.20 K) for radiometers, with slightly higher errors observed for day time data, showing warm biases of between +0.02 (0.39 K) for buoys to +0.11 K (0.33 K) for radiometers. They show that the ATSR series of instruments continues to be the world leader in delivering accurate space-based observations of SST, which is a key climate parameter.

Do early father–infant interactions predict the onset of externalising behaviours in young children? Findings from a longitudinal cohort study

Background Factors related to parents and parenting capacities are important predictors of the development of behavioural problems in children. Recently, there has been an increasing research focus in this field on the earliest years of life, however, relatively few studies have addressed the role of fathers, despite this appearing to be particularly pertinent to child behavioural development. This study aimed to examine whether father–infant interactions at age 3 months independently predicted child behavioural problems at 1 year of age. Method A sample of 192 families was recruited from two maternity units in the United Kingdom. Father–infant interactions were assessed in the family home and coded using the Global Rating Scales. Child behaviour problems were assessed by maternal report. Hierarchical and logistic regression analyses were used to examine associations between father–infant interaction and the development of behavioural problems. Results Disengaged and remote interactions between fathers and their infants were found to predict externalising behavioural problems at the age of 1 year. The children of the most disengaged fathers had an increased risk of developing early externalising behavioural problems [disengaged (nonintrusive) interactions – adjusted Odds Ratio 5.33 (95% Confidence Interval; 1.39, 20.40): remote interactions adj. OR 3.32 (0.92, 12.05)] Conclusions Disengaged interactions of fathers with their infants, as early as the third month of life, predict early behavioural problems in children. These interactions may be critical factors to address, from a very early age in the child’s life, and offer a potential opportunity for preventive intervention.

In-situ spatio-temporal measurements of the detailed azimuthal substructure of the substorm current wedge

The substorm current wedge (SCW) is a fundamental component of geomagnetic substorms. Models tend to describe the SCW as a simple line current flowing into the ionosphere towards dawn and out of the ionosphere towards dusk, linked by a westward electrojet. We use multi-spacecraft observations from perigee passes of the Cluster 1 and 4 spacecraft during a substorm on 15 Jan 2010, in conjunction with ground-based observations, to examine the spatial structuring and temporal variability of the SCW. At this time, the spacecraft travelled east-west azimuthally above the auroral region. We show that the SCW has significant azimuthal sub-structure on scales of 100~km at altitudes of 4,000-7,000~km. We identify 26 individual current sheets in the Cluster 4 data and 34 individual current sheets in the Cluster 1 data, with Cluster 1 passing through the SCW 120-240~s after Cluster 4 at 1,300-2,000~km higher altitude. Both spacecraft observed large-scale regions of net upward and downward field-aligned current, consistent with the large-scale characteristics of the SCW, although sheets of oppositely directed currents were observed within both regions. We show that the majority of these current sheets were closely aligned to a north-south direction, in contrast to the expected east-west orientation of the pre-onset aurora. Comparing our results with observations of the field-aligned current associated with bursty bulk flows (BBFs) we conclude that significant questions remain for the explanation of SCW structuring by BBF driven ``wedgelets". Our results therefore represent constraints on future modelling and theoretical frameworks on the generation of the SCW.

Folding of graphene slit like pore walls—a simple method of improving CO2 separation from mixtures with CH4 or N2

We report for the first time a detailed procedure for creating a simulation model of energetically stable, folded graphene-like pores and simulation results of CO2/CH4 and CO2/N2 separation using these structures. We show that folding of graphene structures is a very promising method to improve the separation of CO2 from mixtures with CH4 and N2. The separation properties of the analyzed materials are compared with carbon nanotubes having similar diameters or S/V ratio. The presented results have potential importance in the field of CO2 capture and sequestration.

Earth's ion upflow associated with polar cap patches: global and in-situ observations

We report simultaneous global monitoring of a patch of ionization and in situ observation of ion upflow at the center of the polar cap region during a geomagnetic storm. Our observations indicate strong fluxes of upwelling O+ ions originating from frictional heating produced by rapid antisunward flow of the plasma patch. The statistical results from the crossings of the central polar cap region by Defense Meteorological Satellite Program F16–F18 from 2010 to 2013 confirm that the field-aligned flow can turn upward when rapid antisunward flows appear, with consequent significant frictional heating of the ions, which overcomes the gravity effect. We suggest that such rapidly moving patches can provide an important source of upwelling ions in a region where downward flows are usually expected. These observations give new insight into the processes of ionosphere-magnetosphere coupling.

Sonic to ultrasonic Q of sandstones and limestones: Laboratory measurements at in situ pressures

Laboratory measurements of the attenuation and velocity dispersion of compressional and shear waves at appropriate frequencies, pressures, and temperatures can aid interpretation of seismic and well-log surveys as well as indicate absorption mechanisms in rocks. Construction and calibration of resonant-bar equipment was used to measure velocities and attenuations of standing shear and extensional waves in copper-jacketed right cylinders of rocks (30 cm in length, 2.54 cm in diameter) in the sonic frequency range and at differential pressures up to 65 MPa. We also measured ultrasonic velocities and attenuations of compressional and shear waves in 50-mm-diameter samples of the rocks at identical pressures. Extensional-mode velocities determined from the resonant bar are systematically too low, yielding unreliable Poisson's ratios. Poisson's ratios determined from the ultrasonic data are frequency corrected and used to calculate the sonic-frequency compressional-wave velocities and attenuations from the shear- and extensional-mode data. We calculate the bulk-modulus loss. The accuracies of attenuation data (expressed as 1000/Q, where Q is the quality factor) are +/- 1 for compressional and shear waves at ultrasonic frequency, +/- 1 for shear waves, and +/- 3 for compressional waves at sonic frequency. Example sonic-frequency data show that the energy absorption in a limestone is small (Q(P) greater than 200 and stress independent) and is primarily due to poroelasticity, whereas that in the two sandstones is variable in magnitude (Q(P) ranges from less than 50 to greater than 300, at reservoir pressures) and arises from a combination of poroelasticity and viscoelasticity. A graph of compressional-wave attenuation versus compressional-wave velocity at reservoir pressures differentiates high-permeability (> 100 mD, 9.87 X 10(-14) m(2)) brine-saturated sandstones from low-permeability (< 100 mD, 9.87 X 10 (14) m(2)) sandstones and shales.

Assessing the sources and magnitude of diurnal nitrate variability in the San Joaquin River (California) with an in-situ optical nitrate sensor and dual nitrate isotopes

We investigated diurnal nitrate (NO3-) concentration variability in the San Joaquin River using an in situ optical NO3- sensor and discrete sampling during a 5-day summer period characterized by high algal productivity. Dual NO3- isotopes (delta N-15(NO3) and delta O-18(NO3)) and dissolved oxygen isotopes (delta O-18(DO)) were measured over 2 days to assess NO3- sources and biogeochemical controls over diurnal time-scales. Concerted temporal patterns of dissolved oxygen (DO) concentrations and delta O-18(DO) were consistent with photosynthesis, respiration and atmospheric O-2 exchange, providing evidence of diurnal biological processes independent of river discharge. Surface water NO3- concentrations varied by up to 22% over a single diurnal cycle and up to 31% over the 5-day study, but did not reveal concerted diurnal patterns at a frequency comparable to DO concentrations. The decoupling of delta N-15(NO3) and delta O-18(NO3) isotopes suggests that algal assimilation and denitrification are not major processes controlling diurnal NO3- variability in the San Joaquin River during the study. The lack of a clear explanation for NO3- variability likely reflects a combination of riverine biological processes and time-varying physical transport of NO3- from upstream agricultural drains to the mainstem San Joaquin River. The application of an in situ optical NO3- sensor along with discrete samples provides a view into the fine temporal structure of hydrochemical data and may allow for greater accuracy in pollution assessment.

Suprathermal electron evolution in a Parker spiral magnetic field

Suprathermal electrons (>70 eV) form a small fraction of the total solar wind electron density but serve as valuable tracers of heliospheric magnetic field topology. Their usefulness as tracers of magnetic loops with both feet rooted on the Sun, however, most likely fades as the loops expand beyond some distance owing to scattering. As a first step toward quantifying that distance, we construct an observationally constrained model for the evolution of the suprathermal electron pitch-angle distributions on open field lines. We begin with a near-Sun isotropic distribution moving antisunward along a Parker spiral magnetic field while conserving magnetic moment, resulting in a field-aligned strahl within a few solar radii. Past this point, the distribution undergoes little evolution with heliocentric distance. We then add constant (with heliocentric distance, energy, and pitch angle) ad-hoc pitch-angle scattering. Close to the Sun, pitch-angle focusing still dominates, again resulting in a narrow strahl. Farther from the Sun, however, pitch-angle scattering dominates because focusing is effectively weakened by the increasing angle between the magnetic field direction and intensity gradient, a result of the spiral field. We determine the amount of scattering required to match Ulysses observations of strahl width in the fast solar wind, providing an important tool for inferring the large-scale properties and topologies of field lines in the interplanetary medium. Although the pitch-angle scattering term is independent of energy, time-of-flight effects in the spiral geometry result in an energy dependence of the strahl width that is in the observed sense although weaker in magnitude.