A Markov Chain Monte Carlo Algorithm for analysis of low signal-to-noise CMB data

We present a new Monte Carlo Markov Chain algorithm for CMB analysis in the low signal-to-noise regime. This method builds on and complements the previously described CMB Gibbs sampler, and effectively solves the low signal-to-noise inefficiency problem of the direct Gibbs sampler. The new algorithm is a simple Metropolis-Hastings sampler with a general proposal rule for the power spectrum, C_l, followed by a particular deterministic rescaling operation of the sky signal. The acceptance probability for this joint move depends on the sky map only through the difference of chi-squared between the original and proposed sky sample, which is close to unity in the low signal-to-noise regime. The algorithm is completed by alternating this move with a standard Gibbs move. Together, these two proposals constitute a computationally efficient algorithm for mapping out the full joint CMB posterior, both in the high and low signal-to-noise regimes.Comment: Submitted to Ap

The two-and three-point correlation functions of the polarized five-year WMAP sky maps

We present the two- and three-point real space correlation functions of the five-year WMAP sky maps, and compare the observed functions to simulated LCDM concordance model ensembles. In agreement with previously published results, we find that the temperature correlation functions are consistent with expectations. However, the pure polarization correlation functions are acceptable only for the 33GHz band map; the 41, 61, and 94 GHz band correlation functions all exhibit significant large-scale excess structures. Further, these excess structures very closely match the correlation functions of the two (synchrotron and dust) foreground templates used to correct the WMAP data for galactic contamination, with a cross-correlation statistically significant at the 2sigma-3sigma confidence level. The correlation is slightly stronger with respect to the thermal dust template than with the synchrotron template.Comment: 10 pages, 5 figures, published in ApJ. v2: New title, minor changes to appendix, and fixed some typos. v3: Matches version published in Ap

In field N transfer, build-up, and leaching in ryegrass-clover mixtures

Two field experiments investigating dynamics in grass-clover mixtures were conducted, using 15N- and 14C-labelling to trace carbon (C) and nitrogen (N) from grass (Lolium perenne L.) and clover (Trifolium repens L. and Trifolium pratense L.). The leaching of dissolved inorganic nitrogen (DIN), as measured in pore water sampled by suction cups, increased during the autumn and winter, whereas the leaching of dissolved organic nitrogen (DON) was fairly constant during this period. Leaching of 15N from the sward indicated that ryegrass was the direct source of less than 1-2 percent of the total N leaching measured, whereas N dynamics pointed to clover as an important contributor to N leaching. Sampling of roots indicates that the dynamics in smaller roots were responsible for N and C build-up in the sward, and that N became available for transfer among species and leaching from the root zone. The bi-directional transfer of N between ryegrass and clover could however not be explained only by root turnover. Other processes like direct uptake of organic N compounds, may have contributed

CMB likelihood approximation by a Gaussianized Blackwell-Rao estimator

We introduce a new CMB temperature likelihood approximation called the Gaussianized Blackwell-Rao (GBR) estimator. This estimator is derived by transforming the observed marginal power spectrum distributions obtained by the CMB Gibbs sampler into standard univariate Gaussians, and then approximate their joint transformed distribution by a multivariate Gaussian. The method is exact for full-sky coverage and uniform noise, and an excellent approximation for sky cuts and scanning patterns relevant for modern satellite experiments such as WMAP and Planck. A single evaluation of this estimator between l=2 and 200 takes ~0.2 CPU milliseconds, while for comparison, a single pixel space likelihood evaluation between l=2 and 30 for a map with ~2500 pixels requires ~20 seconds. We apply this tool to the 5-year WMAP temperature data, and re-estimate the angular temperature power spectrum, $C_{\ell}$, and likelihood, L(C_l), for l<=200, and derive new cosmological parameters for the standard six-parameter LambdaCDM model. Our spectrum is in excellent agreement with the official WMAP spectrum, but we find slight differences in the derived cosmological parameters. Most importantly, the spectral index of scalar perturbations is n_s=0.973 +/- 0.014, 1.9 sigma away from unity and 0.6 sigma higher than the official WMAP result, n_s = 0.965 +/- 0.014. This suggests that an exact likelihood treatment is required to higher l's than previously believed, reinforcing and extending our conclusions from the 3-year WMAP analysis. In that case, we found that the sub-optimal likelihood approximation adopted between l=12 and 30 by the WMAP team biased n_s low by 0.4 sigma, while here we find that the same approximation between l=30 and 200 introduces a bias of 0.6 sigma in n_s.Comment: 10 pages, 7 figures, submitted to Ap

15Nitrogen uptake from shallow- versus deep-rooted plants in multi-species mixtures and monoculture grassland

Only few studies have explored the importance of functional diversity in temperate agricultural grasslands in relation to nitrogen (N) uptake. This study investigates the consequence of growing deep-rooted plants together with grass-clover mixtures in terms of N uptake efficiency from deep soil layers. The objective was to compare the N uptake of the shallow-rooted grassland species Lolium perenne and Trifolium repens; and the deep-rooted species Cichorium intybus and Medicago sativa in monocultures and mixtures. We hypothesized that growing deep-rooted plant species in mixture with shallow-rooted species increases the N uptake from deep soil layers partly through competition. A 15N tracer study was carried out with 15N enriched ammonium-sulphate placed at three different soil depths (40, 80 and 120 cm). To recover 15N, above-ground plant biomass was harvested after 10 days. We described the decline of 15N uptake with depth by using an exponential decay function. The studied plant communities showed the same relative decline in 15N uptake by increasing soil depths, but different capacities in total 15N uptake. Monoculture L. perenne foraged less 15N in all depths compared to the other four plant communities. The relative 15N uptake of individual plant species grown in mixture decreased stronger with depth than in monoculture. Thus, both findings rejected our hypothesis

The effects of the exhaust plume on the lightning triggering conditions for launch vehicles

Apollo 12 and Atlas Centaur 67 are two launch vehicles that have experienced triggered lightning strikes. Serious consequences resulted from the events; in the case of Atlas Centaur 67, the vehicle and the payload were lost. These events indicate that it is necessary to develop launch rules which would prevent such occurrences. In order to develop valid lightning related rules, it is necessary to understand the effects of the plume. Some have assumed that the plume can be treated as a perfect conductor, and have computed electric field enhancement factors on that basis. The authors have looked at the plume, and believe that these models are not correct, because they ignore the fluid motion of the conducting plates. The authors developed a model which includes this flow character. In this model, the external field is excluded from the plume as it would be for any good conductor, but, in addition, the charge must distribute so that the charge density is zero at some location in the exhaust. When this condition is included in the calculation of triggering enhancement factors, they can be two to three times larger than calculated by other methods which include a conductive plume but don't include the correct boundary conditions. Here, the authors review the relevant features of rocket exhausts for the triggered lightning problem, present an approach for including flowing conductive gases, and present preliminary calculations to demonstrate the effect that the plume has on enhancement factors

15Nitrogen transfer from legumes to neighbouring plants in multi-species grassland

This study investigates the N transfer from legumes to neighbouring plants, grasses, legumes and herbs in a temperate grassland. In a field experiment white clover (Trifolium repens), red clover (Trifolium pratense) and lucerne (Medicago sativa) were leaf-labelled with 15N enriched urea. The 15N tracer was measured in above-ground plant tissue of eight neighbouring plants in two subsequent harvests in 2008. The three legumes donated 15N to all neighbouring plants, of which grasses, white and red clover were strong receivers. Results show that N transfer increases with N application and from the 1st to the 2nd cut

Constraints on mode couplings and modulation of the CMB with WMAP data

We investigate a possible asymmetry in the statistical properties of the cosmic microwave background temperature field and to do so we construct an estimator aiming at detecting a dipolar modulation. Such a modulation is found to induce correlations between multipoles with $\Delta\ell=1$. Applying this estimator, to the V and W bands of the WMAP data, we found a significant detection in the V band. We argue however that foregrounds and in particular point sources are the origin of this signal.Comment: 14 pages, 14 figure