Estimating subseasonal variability and trends in global atmosphere using reanalysis data

A new measure of subseasonal variability is introduced that provides a scale-dependent estimation of vertically and meridionally integrated atmospheric variability in terms of the normal modes of linearized primitive equations. Applied to the ERA-Interim data, the new measure shows that subseasonal variability decreases for larger zonal wave numbers. Most of variability is due to balanced (Rossby mode) dynamics but the portion associated with the inertio-gravity (IG) modes increases as the scale reduces. Time series of globally integrated variability anomalies in ERA-Interim show an increase in variability after year 2000. In recent years the anomalies have been about 2% above the 1981–2010 average. The relative increase in variability projecting on the IG modes is larger and more persistent than for the Rossby modes. Although the IG part is a small component of the subseasonal variability, it is an important effect likely reflecting the observed increase in the tropical precipitation variability. ©2018. The Authors

Galactic interstellar filaments as probed by LOFAR and Planck

Recent Low Frequency Array (LOFAR) observations at 115-175 MHz of a field at medium Galactic latitudes (centered at the bright quasar 3C196) have shown striking filamentary structures in polarization that extend over more than 4 degrees across the sky. In addition, the Planck satellite has released full sky maps of the dust emission in polarization at 353GHz. The LOFAR data resolve Faraday structures along the line of sight, whereas the Planck dust polarization maps probe the orientation of the sky projected magnetic field component. Hence, no apparent correlation between the two is expected. Here we report a surprising, yet clear, correlation between the filamentary structures, detected with LOFAR, and the magnetic field orientation, probed by the Planck satellite. This finding points to a common, yet unclear, physical origin of the two measurements in this specific area in the sky. A number of follow-up multi- frequency studies are proposed to shed light on this unexpected finding.Comment: 6 pages, 4 figures, accepted for publication in MNRAS Letter

Foregrounds for observations of the cosmological 21 cm line: II. Westerbork observations of the fields around 3C196 and the North Celestial Pole

In the coming years a new insight into galaxy formation and the thermal history of the Universe is expected to come from the detection of the highly redshifted cosmological 21 cm line. The cosmological 21 cm line signal is buried under Galactic and extragalactic foregrounds which are likely to be a few orders of magnitude brighter. Strategies and techniques for effective subtraction of these foreground sources require a detailed knowledge of their structure in both intensity and polarization on the relevant angular scales of 1-30 arcmin. We present results from observations conducted with the Westerbork telescope in the 140-160 MHz range with 2 arcmin resolution in two fields located at intermediate Galactic latitude, centred around the bright quasar 3C196 and the North Celestial Pole. They were observed with the purpose of characterizing the foreground properties in sky areas where actual observations of the cosmological 21 cm line could be carried out. The polarization data were analysed through the rotation measure synthesis technique. We have computed total intensity and polarization angular power spectra. Total intensity maps were carefully calibrated, reaching a high dynamic range, 150000:1 in the case of the 3C196 field. [abridged]Comment: 20 pages, 22 figures, accepted for publication in A&A. A version with full resolution figures is available at http://www.astro.rug.nl/~bernardi/NCP_3C196/bernardi.pd

Fast Large-Scale Reionization Simulations

We present an efficient method to generate large simulations of the Epoch of Reionization (EoR) without the need for a full 3-dimensional radiative transfer code. Large dark-matter-only simulations are post-processed to produce maps of the redshifted 21cm emission from neutral hydrogen. Dark matter haloes are embedded with sources of radiation whose properties are either based on semi-analytical prescriptions or derived from hydrodynamical simulations. These sources could either be stars or power-law sources with varying spectral indices. Assuming spherical symmetry, ionized bubbles are created around these sources, whose radial ionized fraction and temperature profiles are derived from a catalogue of 1-D radiative transfer experiments. In case of overlap of these spheres, photons are conserved by redistributing them around the connected ionized regions corresponding to the spheres. The efficiency with which these maps are created allows us to span the large parameter space typically encountered in reionization simulations. We compare our results with other, more accurate, 3-D radiative transfer simulations and find excellent agreement for the redshifts and the spatial scales of interest to upcoming 21cm experiments. We generate a contiguous observational cube spanning redshift 6 to 12 and use these simulations to study the differences in the reionization histories between stars and quasars. Finally, the signal is convolved with the LOFAR beam response and its effects are analyzed and quantified. Statistics performed on this mock data set shed light on possible observational strategies for LOFAR.Comment: 18 pages, 21 figures, submitted to MNRAS For high-resolution images follow "http://www.astro.rug.nl/~thomas/eormap.pdf

Detection and extraction of signals from the epoch of reionization using higher-order one-point statistics

Detecting redshifted 21-cm emission from neutral hydrogen in the early Universe promises to give direct constraints on the epoch of reionization (EoR). It will, though, be very challenging to extract the cosmological signal (CS) from foregrounds and noise which are orders of magnitude larger. Fortunately, the signal has some characteristics which differentiate it from the foregrounds and noise, and we suggest that using the correct statistics may tease out signatures of reionization. We generate mock data cubes simulating the output of the Low Frequency Array (LOFAR) EoR experiment. These cubes combine realistic models for Galactic and extragalactic foregrounds and the noise with three different simulations of the CS. We fit out the foregrounds, which are smooth in the frequency direction, to produce residual images in each frequency band. We denoise these images and study the skewness of the one-point distribution in the images as a function of frequency. We find that, under sufficiently optimistic assumptions, we can recover the main features of the redshift evolution of the skewness in the 21-cm signal. We argue that some of these features ¿ such as a dip at the onset of reionization, followed by a rise towards its later stages ¿ may be generic, and give us a promising route to a statistical detection of reionization

GMRT observation towards detecting the Post-reionization 21-cm signal

We have analyzed 610 MHz GMRT observations towards detecting the redshifted 21-cm signal from z=1.32. The multi-frequency angular power spectrum C_l(Delta nu) is used to characterize the statistical properties of the background radiation across angular scales ~20" to 10', and a frequency bandwidth of 7.5 MHz with resolution 125 kHz. The measured C_l(Delta nu) which ranges from 7 mK^2 to 18 mK^2 is dominated by foregrounds, the expected HI signal C_l^HI(Delta nu) ~10^{-6}- 10^{-7} mK^2 is several orders of magnitude smaller. The foregrounds, believed to originate from continuum sources, is expected to vary smoothly with Delta nu whereas the HI signal decorrelates within ~0.5 MHz and this holds the promise of separating the two. For each l, we use the interval 0.5 < Delta nu < 7.5 MHz to fit a fourth order polynomial which is subtracted from the measured C_l(Delta nu) to remove any smoothly varying component across the entire bandwidth Delta nu < 7.5 MHz. The residual C_l(Delta nu), we find, has an oscillatory pattern with amplitude and period respectively ~0.1 mK^2 and Delta nu = 3 MHz at the smallest l value of 1476, and the amplitude and period decreasing with increasing l. Applying a suitably chosen high pass filter, we are able to remove the residual oscillatory pattern for l=1476 where the residual C_l(Delta nu) is now consistent with zero at the 3-sigma noise level. We conclude that we have successfully removed the foregrounds at l=1476 and the residuals are consistent with noise. We use this to place an upper limit on the HI signal whose amplitude is determined by x_HI b where x_HI and b are the HI neutral fraction and the HI bias respectively. A value of x_HI b greater than 7.95 would have been detected in our observation, and is therefore ruled out at the 3-sigma level. (abridged)Comment: 29 pages, 13 figures, Accepted to MNRA

Constraining the epoch of reionization with the variance statistic: simulations of the LOFAR case

Several experiments are underway to detect the cosmic redshifted 21-cm signal from neutral hydrogen from the Epoch of Reionization (EoR). Due to their very low signal-to-noise ratio, these observations aim for a statistical detection of the signal by measuring its power spectrum. We investigate the extraction of the variance of the signal as a first step towards detecting and constraining the global history of the EoR. Signal variance is the integral of the signal's power spectrum, and it is expected to be measured with a high significance. We demonstrate this through results from a simulation and parameter estimation pipeline developed for the Low Frequency Array (LOFAR)-EoR experiment. We show that LOFAR should be able to detect the EoR in 600 hours of integration using the variance statistic. Additionally, the redshift ($z_r$) and duration ($\Delta z$) of reionization can be constrained assuming a parametrization. We use an EoR simulation of $z_r = 7.68$ and $\Delta z = 0.43$ to test the pipeline. We are able to detect the simulated signal with a significance of 4 standard deviations and extract the EoR parameters as $z_r = 7.72^{+0.37}_{-0.18}$ and $\Delta z = 0.53^{+0.12}_{-0.23}$ in 600 hours, assuming that systematic errors can be adequately controlled. We further show that the significance of detection and constraints on EoR parameters can be improved by measuring the cross-variance of the signal by cross-correlating consecutive redshift bins.Comment: 13 pages, 14 figures, Accepted for publication in MNRA

LOFAR Deep Fields: Probing faint Galactic polarised emission in ELAIS-N1

We present the first deep polarimetric study of Galactic synchrotron emission at low radio frequencies. Our study is based on 21 observations of the European Large Area Infrared Space Observatory Survey-North 1 (ELAIS-N1) field using the Low-Frequency Array (LOFAR) at frequencies from 114.9 to 177.4 MHz. These data are a part of the LOFAR Two-metre Sky Survey Deep Fields Data Release 1. We used very low-resolution ($4.3'$) Stokes QU data cubes of this release. We applied rotation measure (RM) synthesis to decompose the distribution of polarised structures in Faraday depth, and cross-correlation RM synthesis to align different observations in Faraday depth. We stacked images of about 150 hours of the ELAIS-N1 observations to produce the deepest Faraday cube at low radio frequencies to date, tailored to studies of Galactic synchrotron emission and the intervening magneto-ionic interstellar medium. This Faraday cube covers $\sim36~{\rm deg^{2}}$ of the sky and has a noise of $27~{\rm \mu Jy~PSF^{-1}~RMSF^{-1}}$ in polarised intensity. This is an improvement in noise by a factor of approximately the square root of the number of stacked data cubes ($\sim\sqrt{20}$), as expected, compared to the one in a single data cube based on five-to-eight-hour observations. We detect a faint component of diffuse polarised emission in the stacked cube, which was not detected previously. Additionally, we verify the reliability of the ionospheric Faraday rotation corrections estimated from the satellite-based total electron content measurements to be of $~\sim0.05~{\rm rad~m^{-2}}$. We also demonstrate that diffuse polarised emission itself can be used to account for the relative ionospheric Faraday rotation corrections with respect to a reference observation.Comment: 15 pages, 15 figures, accepted for publication in A&

Foregrounds for observations of the cosmological 21 cm line: I. First Westerbork measurements of Galactic emission at 150 MHz in a low latitude field

We present the first results from a series of observations conducted with the Westerbork telescope in the 140--160 MHz range with a 2 arcmin resolution aimed at characterizing the properties of the foregrounds for epoch of reionization experiments. For the first time we have detected fluctuations in the Galactic diffuse emission on scales greater than 13 arcmin at 150 MHz, in the low Galactic latitude area known as Fan region. Those fluctuations have an $rms$ of 14 K. The total intensity power spectrum shows a power--law behaviour down to $\ell \sim 900$ with slope $\beta^I_\ell = -2.2 \pm 0.3$. The detection of diffuse emission at smaller angular scales is limited by residual point sources. We measured an $rms$ confusion noise of $\sim$3 mJy beam$^{-1}$. Diffuse polarized emission was also detected for the first time at this frequency. The polarized signal shows complex structure both spatially and along the line of sight. The polarization power spectrum shows a power--law behaviour down to $\ell \sim 2700$ with slope $\beta^P_\ell = -1.65 \pm 0.15$. The $rms$ of polarization fluctuations is 7.2 K on 4 arcmin scales. By extrapolating the measured spectrum of total intensity emission, we find a contamination on the cosmological signal of $\delta T= \sqrt{\ell (\ell+1) C^I_\ell / 2\pi} \sim 5.7$ K on 5 arcmin scales and a corresponding $rms$ value of $\sim$18.3 K at the same angular scale. The level of the polarization power spectrum is $\delta T \sim 3.3$ K on 5 arcmin scales. Given its exceptionally bright polarized signal, the Fan region is likely to represent an upper limit on the sky brightness at moderate and high Galactic latitude.Comment: Minor corrections made to match the final version printed on A&A. A version with high resolution figures is available at http://www.astro.rug.nl/~bernardi/FAN/fan.pd