Spatiotemporal heterogeneity of water flowpaths controls dissolved organic carbon sourcing in a snow-dominated, headwater catchment

The non-uniform distribution of water in snowdrift-driven systems can lead to spatial heterogeneity in vegetative communities and soil development, as snowdrifts may locally increase weathering. The focus of this study is to understand the coupled hydrological and biogeochemical dynamics in a heterogeneous, snowdrift-dominated headwater catchment (Reynolds Mountain East, Reynolds Creek Critical Zone Observatory, Idaho, USA). We determine the sources and fluxes of stream water and dissolved organic carbon (DOC) at this site, deducing likely flowpaths from hydrometric and hydrochemical signals of soil water, saprolite water, and groundwater measured through the snowmelt period and summer recession. We then interpret flowpaths using end-member mixing analysis in light of inferred subsurface structure derived from electrical resistivity and seismic velocity transects. Streamwater is sourced primarily from groundwater (averaging 25% of annual streamflow), snowmelt (50%), and water traveling along the saprolite/bedrock boundary (25%). The latter is comprised of the prior year\u27s soil water, which accumulates DOC in the soil matrix through the summer before flushing to the saprolite during snowmelt. DOC indices suggest that it is sourced from terrestrial carbon, and derives originally from soil organic carbon (SOC) before flushing to the saprolite/bedrock boundary. Multiple subsurface regions in the catchment appear to contribute differentially to streamflow as the season progresses; sources shift from the saprolite/bedrock interface to deeper bedrock aquifers from the snowmelt period into summer. Unlike most studied catchments, lateral flow of soil water during the study year is not a primary source of streamflow. Instead, saprolite and groundwater act as integrators of soil water that flows vertically in this system. Our results do not support the flushing hypothesis as observed in similar systems and instead indicate that temporal variation in connectivity may cause the unexpected dilution behavior displayed by DOC in this catchment

An improved formulation of the relativistic hydrodynamics equations in 2D Cartesian coordinates

A number of astrophysical scenarios possess and preserve an overall cylindrical symmetry also when undergoing a catastrophic and nonlinear evolution. Exploiting such a symmetry, these processes can be studied through numerical-relativity simulations at smaller computational costs and at considerably larger spatial resolutions. We here present a new flux-conservative formulation of the relativistic hydrodynamics equations in cylindrical coordinates. By rearranging those terms in the equations which are the sources of the largest numerical errors, the new formulation yields a global truncation error which is one or more orders of magnitude smaller than those of alternative and commonly used formulations. We illustrate this through a series of numerical tests involving the evolution of oscillating spherical and rotating stars, as well as shock-tube tests.Comment: 19 pages, 9 figure

Software "Proasu"

The paper describes software for automation of ASCTP project documentation development ("ProASU"). This software is implemented for development of text documentation in automate mode way that for every of document type (TS, MS, PS, IS, SS, OS) there are predefined individual scenarios of creation

Characteristic extraction in numerical relativity: binary black hole merger waveforms at null infinity

The accurate modeling of gravitational radiation is a key issue for gravitational wave astronomy. As simulation codes reach higher accuracy, systematic errors inherent in current numerical relativity wave-extraction methods become evident, and may lead to a wrong astrophysical interpretation of the data. In this paper, we give a detailed description of the Cauchy-characteristic extraction technique applied to binary black hole inspiral and merger evolutions to obtain gravitational waveforms that are defined unambiguously, that is, at future null infinity. By this method we remove finite-radius approximations and the need to extrapolate data from the near zone. Further, we demonstrate that the method is free of gauge effects and thus is affected only by numerical error. Various consistency checks reveal that energy and angular momentum are conserved to high precision and agree very well with extrapolated data. In addition, we revisit the computation of the gravitational recoil and find that finite radius extrapolation very well approximates the result at \scri. However, the (non-convergent) systematic differences to extrapolated data are of the same order of magnitude as the (convergent) discretisation error of the Cauchy evolution hence highlighting the need for correct wave-extraction.Comment: 41 pages, 8 figures, 2 tables, added references, fixed typos. Version matches published version

Spacelike distance from discrete causal order

Any discrete approach to quantum gravity must provide some prescription as to how to deduce continuum properties from the discrete substructure. In the causal set approach it is straightforward to deduce timelike distances, but surprisingly difficult to extract spacelike distances, because of the unique combination of discreteness with local Lorentz invariance in that approach. We propose a number of methods to overcome this difficulty, one of which reproduces the spatial distance between two points in a finite region of Minkowski space. We provide numerical evidence that this definition can be used to define a `spatial nearest neighbor' relation on a causal set, and conjecture that this can be exploited to define the length of `continuous curves' in causal sets which are approximated by curved spacetime. This provides evidence in support of the ``Hauptvermutung'' of causal sets.Comment: 32 pages, 16 figures, revtex4; journal versio

Selection at a single locus leads to widespread expansion of toxoplasma gondii lineages that are virulent in mice

The determinants of virulence are rarely defined for eukaryotic parasites such as T. gondii, a widespread parasite of mammals that also infects humans, sometimes with serious consequences. Recent laboratory studies have established that variation in a single secreted protein, a serine/threonine kinase known as ROPO18, controls whether or not mice survive infection. Here, we establish the extent and nature of variation in ROP18among a collection of parasite strains from geographically diverse regions. Compared to other genes, ROP18 showed extremely high levels of diversification and changes in expression level, which correlated with severity of infection in mice. Comparison with an out-group demonstrated that changes in the upstream region that regulates expression of ROP18 led to an historical increase in the expression and exposed the protein to diversifying selective pressure. Surprisingly, only three atypically distinct protein variants exist despite marked genetic divergence elsewhere in the genome. These three forms of ROP18 are likely adaptations for different niches in nature, and they confer markedly different virulence to mice. The widespread distribution of a single mouse-virulent allele among geographically and genetically disparate parasites may have consequences for transmission and disease in other hosts, including humans

The Current Status of Binary Black Hole Simulations in Numerical Relativity

Since the breakthroughs in 2005 which have led to long term stable solutions of the binary black hole problem in numerical relativity, much progress has been made. I present here a short summary of the state of the field, including the capabilities of numerical relativity codes, recent physical results obtained from simulations, and improvements to the methods used to evolve and analyse binary black hole spacetimes.Comment: 14 pages; minor changes and corrections in response to referee

Testing gravitational-wave searches with numerical relativity waveforms: Results from the first Numerical INJection Analysis (NINJA) project

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses.Comment: 56 pages, 25 figures; various clarifications; accepted to CQ

Sensitivity to Gravitational Waves from Compact Binary Coalescences Achieved during LIGO's Fifth and Virgo's First Science Run

We summarize the sensitivity achieved by the LIGO and Virgo gravitational wave detectors for compact binary coalescence (CBC) searches during LIGO's fifth science run and Virgo's first science run. We present noise spectral density curves for each of the four detectors that operated during these science runs which are representative of the typical performance achieved by the detectors for CBC searches. These spectra are intended for release to the public as a summary of detector performance for CBC searches during these science runs.Comment: 12 pages, 5 figure