By Shepherd, et all, posted on November 29th, 2013 in Articles, Climate

Earth is increasingly an “urbanized ” planet. The “World Population Clock ” registered a Population of 7,175,309,538 at 8:30 pm (LST) on Oct. 6, 2013. Current and future trends suggest that this population will increasingly reside in cities. Currently, 52 percent of the world population is urban, which means we are a majority “urbanized ” society. Figure 1 indicates this trend will continue, wit

Parallelization of Kinetic Theory Simulations

Numerical studies of shock waves in large scale systems via kinetic simulations with millions of particles are too computationally demanding to be processed in serial. In this work we focus on optimizing the parallel performance of a kinetic Monte Carlo code for astrophysical simulations such as core-collapse supernovae. Our goal is to attain a flexible program that scales well with the architecture of modern supercomputers. This approach requires a hybrid model of programming that combines a message passing interface (MPI) with a multithreading model (OpenMP) in C++. We report on our approach to implement the hybrid design into the kinetic code and show first results which demonstrate a significant gain in performance when many processors are applied.Comment: 10 pages, 3 figures, conference proceeding

Differential gene expression in the nucleus accumbens with ethanol self-administration in inbred alcohol-preferring rats

The current study examined the effects of operant ethanol (EtOH) self-administration on gene expression in the nucleus accumbens (ACB) and amygdala (AMYG) of inbred alcohol-preferring (iP) rats. Rats self-trained on a standard two-lever operant paradigm to administer either water-water, EtOH (15% v/v)-water, or saccharin (SAC; 0.0125% g/v)-water. Animals were killed 24 hr after the last operant session, and the ACB and AMYG dissected; RNA was extracted and purified for microarray analysis. For the ACB, there were 513 significant differences at the p < 0.01 level in named genes: 55 between SAC and water; 215 between EtOH and water, and 243 between EtOH and SAC. In the case of the AMYG (p < 0.01), there were 48 between SAC and water, 23 between EtOH and water, and 63 between EtOH and SAC group. Gene Ontology (GO) analysis indicated that differences in the ACB between the EtOH and SAC groups could be grouped into 15 significant (p < 0.05) categories, which included major categories such as synaptic transmission, cell and ion homeostasis, and neurogenesis, whereas differences between the EtOH and water groups had only 4 categories, which also included homeostasis and synaptic transmission. Several genes were in common between the EtOH and both the SAC and water groups in the synaptic transmission (e.g., Cav2, Nrxn, Gabrb2, Gad1, Homer1) and homeostasis (S100b, Prkca, Ftl1) categories. Overall, the results suggest that changes in gene expression in the ACB of iP rats are associated with the reinforcing effects of EtOH

Changes in Holocene climate and the intensity of Southern Hemisphere Westerly Winds based on a high-resolution palynological record from sub-Antarctic South Georgia

Sub-Antarctic South Georgia is a key region for studying climate variability in the Southern Hemisphere, because of its position at the core of the Southern Hemisphere Westerly Wind belt and between the Antarctic Circumpolar Current and the Polar Frontal Zone. Here, we present a 5.8-m long high-resolution pollen record from Fan Lake on Annenkov Island dominated by local sub-polar vegetation, with Acaena and Poaceae being present throughout the last 7000 years. Palynological and sedimentological analyses revealed a warm late Holocene ‘climate optimum’ between 3790 and 2750 cal. yr BP, which was followed by a gradual transition to cool and wet conditions. This cooling was interrupted by slightly warmer environmental conditions between 1670 and 710 cal. yr BP that partly overlap with the Northern Hemisphere ‘Medieval Climate Anomaly’. Increases in non-native and long-distance pollen grains transported from South America (e.g. Nothofagus, Podocarpus) indicate that stronger Southern Hemisphere Westerly Winds over South Georgia possibly occurred during some ‘colder’ phases of the late Holocene, most notably between c. 2210 and 1670 cal. yr BP and after 710 cal. yr BP

High potential for weathering and climate effects of non-vascular vegetation in the Late Ordovician

It has been hypothesized that predecessors of today’s bryophytes significantly increased global chemical weathering in the Late Ordovician, thus reducing atmospheric CO2 concentration and contributing to climate cooling and an interval of glaciations. Studies that try to quantify the enhancement of weathering by non-vascular vegetation, however, are usually limited to small areas and low numbers of species, which hampers extrapolating to the global scale and to past climatic conditions. Here we present a spatially explicit modelling approach to simulate global weathering by non-vascular vegetation in the Late Ordovician. We estimate a potential global weathering flux of 2.8 (km3 rock) yr−1, defined here as volume of primary minerals affected by chemical transformation. This is around three times larger than today’s global chemical weathering flux. Moreover, we find that simulated weathering is highly sensitive to atmospheric CO2 concentration. This implies a strong negative feedback between weathering by non-vascular vegetation and Ordovician climate

The timescale of early land plant evolution

Establishing the timescale of early land plant evolution is essential for testing hypotheses on the coevolution of land plants and Earth's System. The sparseness of early land plant megafossils and stratigraphic controls on their distribution make the fossil record an unreliable guide, leaving only the molecular clock. However, the application of molecular clock methodology is challenged by the current impasse in attempts to resolve the evolutionary relationships among the living bryophytes and tracheophytes. Here, we establish a timescale for early land plant evolution that integrates over topological uncertainty by exploring the impact of competing hypotheses on bryophyte-tracheophyte relationships, among other variables, on divergence time estimation. We codify 37 fossil calibrations for Viridiplantae following best practice. We apply these calibrations in a Bayesian relaxed molecular clock analysis of a phylogenomic dataset encompassing the diversity of Embryophyta and their relatives within Viridiplantae. Topology and dataset sizes have little impact on age estimates, with greater differences among alternative clock models and calibration strategies. For all analyses, a Cambrian origin of Embryophyta is recovered with highest probability. The estimated ages for crown tracheophytes range from Late Ordovician to late Silurian. This timescale implies an early establishment of terrestrial ecosystems by land plants that is in close accord with recent estimates for the origin of terrestrial animal lineages. Biogeochemical models that are constrained by the fossil record of early land plants, or attempt to explain their impact, must consider the implications of a much earlier, middle Cambrian-Early Ordovician, origin

Validation of Coevolving Residue Algorithms via Pipeline Sensitivity Analysis: ELSC and OMES and ZNMI, Oh My!

Correlated amino acid substitution algorithms attempt to discover groups of residues that co-fluctuate due to either structural or functional constraints. Although these algorithms could inform both ab initio protein folding calculations and evolutionary studies, their utility for these purposes has been hindered by a lack of confidence in their predictions due to hard to control sources of error. To complicate matters further, naive users are confronted with a multitude of methods to choose from, in addition to the mechanics of assembling and pruning a dataset. We first introduce a new pair scoring method, called ZNMI (Z-scored-product Normalized Mutual Information), which drastically improves the performance of mutual information for co-fluctuating residue prediction. Second and more important, we recast the process of finding coevolving residues in proteins as a data-processing pipeline inspired by the medical imaging literature. We construct an ensemble of alignment partitions that can be used in a cross-validation scheme to assess the effects of choices made during the procedure on the resulting predictions. This pipeline sensitivity study gives a measure of reproducibility (how similar are the predictions given perturbations to the pipeline?) and accuracy (are residue pairs with large couplings on average close in tertiary structure?). We choose a handful of published methods, along with ZNMI, and compare their reproducibility and accuracy on three diverse protein families. We find that (i) of the algorithms tested, while none appear to be both highly reproducible and accurate, ZNMI is one of the most accurate by far and (ii) while users should be wary of predictions drawn from a single alignment, considering an ensemble of sub-alignments can help to determine both highly accurate and reproducible couplings. Our cross-validation approach should be of interest both to developers and end users of algorithms that try to detect correlated amino acid substitutions

A Study of Time-Dependent CP-Violating Asymmetries and Flavor Oscillations in Neutral B Decays at the Upsilon(4S)

We present a measurement of time-dependent CP-violating asymmetries in neutral B meson decays collected with the BABAR detector at the PEP-II asymmetric-energy B Factory at the Stanford Linear Accelerator Center. The data sample consists of 29.7 ${\rm fb}^{-1}$ recorded at the $\Upsilon(4S)$ resonance and 3.9 ${\rm fb}^{-1}$ off-resonance. One of the neutral B mesons, which are produced in pairs at the $\Upsilon(4S)$, is fully reconstructed in the CP decay modes $J/\psi K^0_S$, $\psi(2S) K^0_S$, $\chi_{c1} K^0_S$, $J/\psi K^{*0}$ ($K^{*0}\to K^0_S\pi^0$) and $J/\psi K^0_L$, or in flavor-eigenstate modes involving $D^{(*)}\pi/\rho/a_1$ and $J/\psi K^{*0}$ ($K^{*0}\to K^+\pi^-$). The flavor of the other neutral B meson is tagged at the time of its decay, mainly with the charge of identified leptons and kaons. The proper time elapsed between the decays is determined by measuring the distance between the decay vertices. A maximum-likelihood fit to this flavor eigenstate sample finds $\Delta m_d = 0.516\pm 0.016 {\rm (stat)} \pm 0.010 {\rm (syst)} {\rm ps}^{-1}$. The value of the asymmetry amplitude $\sin2\beta$ is determined from a simultaneous maximum-likelihood fit to the time-difference distribution of the flavor-eigenstate sample and about 642 tagged $B^0$ decays in the CP-eigenstate modes. We find $\sin2\beta=0.59\pm 0.14 {\rm (stat)} \pm 0.05 {\rm (syst)}$, demonstrating that CP violation exists in the neutral B meson system. (abridged)Comment: 58 pages, 35 figures, submitted to Physical Review

Evidence for the Rare Decay B -> K*ll and Measurement of the B -> Kll Branching Fraction

We present evidence for the flavor-changing neutral current decay $B\to K^*\ell^+\ell^-$ and a measurement of the branching fraction for the related process $B\to K\ell^+\ell^-$, where $\ell^+\ell^-$ is either an $e^+e^-$ or $\mu^+\mu^-$ pair. These decays are highly suppressed in the Standard Model, and they are sensitive to contributions from new particles in the intermediate state. The data sample comprises $123\times 10^6$ $\Upsilon(4S)\to B\bar{B}$ decays collected with the Babar detector at the PEP-II $e^+e^-$ storage ring. Averaging over $K^{(*)}$ isospin and lepton flavor, we obtain the branching fractions ${\mathcal B}(B\to K\ell^+\ell^-)=(0.65^{+0.14}_{-0.13}\pm 0.04)\times 10^{-6}$ and ${\mathcal B}(B\to K^*\ell^+\ell^-)=(0.88^{+0.33}_{-0.29}\pm 0.10)\times 10^{-6}$, where the uncertainties are statistical and systematic, respectively. The significance of the $B\to K\ell^+\ell^-$ signal is over $8\sigma$, while for $B\to K^*\ell^+\ell^-$ it is $3.3\sigma$.Comment: 7 pages, 2 postscript figues, submitted to Phys. Rev. Let

Sea Surface Temperature of the mid-Piacenzian Ocean:A Data-Model Comparison

The mid-Piacenzian climate represents the most geologically recent interval of long-term average warmth relative to the last million years, and shares similarities with the climate projected for the end of the 21st century. As such, it represents a natural experiment from which we can gain insight into potential climate change impacts, enabling more informed policy decisions for mitigation and adaptation. Here, we present the first systematic comparison of Pliocene sea surface temperature (SST) between an ensemble of eight climate model simulations produced as part of PlioMIP (Pliocene Model Intercomparison Project) with the PRISM (Pliocene Research, Interpretation and Synoptic Mapping) Project mean annual SST field. Our results highlight key regional and dynamic situations where there is discord between the palaeoenvironmental reconstruction and the climate model simulations. These differences have led to improved strategies for both experimental design and temporal refinement of the palaeoenvironmental reconstruction