Search

We describe an extension of the most recent version of the Planck Catalogue of Compact Sources (PCCS2), produced using a new multi-band Bayesian Extraction and Estimation Package (BeeP). BeeP assumes that the compact sources present in PCCS2 at 857 GHz have a dust-like spectral energy distribution (SED), which leads to emission at both lower and higher frequencies, and adjusts the parameters of the source and its SED to fit the emission observed in Planck's three highest frequency channels at 353, 545, and 857 GHz, as well as the IRIS map at 3000 GHz. In order to reduce confusion regarding diffuse cirrus emission, BeeP's data model includes a description of the background emission surrounding each source, and it adjusts the confidence in the source parameter extraction based on the statistical properties of the spatial distribution of the background emission. BeeP produces the following three new sets of parameters for each source: (a) fits to a modified blackbody (MBB) thermal emission model of the source; (b) SED-independent source flux densities at each frequency considered; and (c) fits to an MBB model of the background in which the source is embedded. BeeP also calculates, for each source, a reliability parameter, which takes into account confusion due to the surrounding cirrus. This parameter can be used to extract sub-samples of high-frequency sources with statistically well-understood properties. We define a high-reliability subset (BeeP/base), containing 26 083 sources (54.1% of the total PCCS2 catalogue), the majority of which have no information on reliability in the PCCS2. We describe the characteristics of this specific high-quality subset of PCCS2 and its validation against other data sets, specifically for: the sub-sample of PCCS2 located in low-cirrus areas; the Planck Catalogue of Galactic Cold Clumps; the Herschel GAMA15-field catalogue; and the temperature-and spectral-index-reconstructed dust maps obtained with Planck's Generalized Needlet Internal Linear Combination method. The results of the BeeP extension of PCCS2, which are made publicly available via the Planck Legacy Archive, will enable the study of the thermal properties of well-defined samples of compact Galactic and extragalactic dusty sources.

The largest temperature anisotropy in the cosmic microwave background (CMB) is the dipole, which has been measured with increasing accuracy for more than three decades, particularly with the Planck satellite. The simplest interpretation of the dipole is that it is due to our motion with respect to the rest frame of the CMB. Since current CMB experiments infer temperature anisotropies from angular intensity variations, the dipole modulates the temperature anisotropies with the same frequency dependence as the thermal Sunyaev-Zeldovich (tSZ) effect. We present the first, and significant, detection of this signal in the tSZ maps and find that it is consistent with direct measurements of the CMB dipole, as expected. The signal contributes power in the tSZ maps, which is modulated in a quadrupolar pattern, and we estimate its contribution to the tSZ bispectrum, noting that it contributes negligible noise to the bispectrum at relevant scales.

We describe an extension of the most recent version of the Planck Catalogue of Compact Sources (PCCS2), produced using a new multi-band Bayesian Extraction and Estimation Package (BeeP). BeeP assumes that the compact sources present in PCCS2 at 857 GHz have a dust-like spectral energy distribution (SED), which leads to emission at both lower and higher frequencies, and adjusts the parameters of the source and its SED to fit the emission observed in Planck's three highest frequency channels at 353, 545, and 857 GHz, as well as the IRIS map at 3000 GHz. In order to reduce confusion regarding diffuse cirrus emission, BeeP's data model includes a description of the background emission surrounding each source, and it adjusts the confidence in the source parameter extraction based on the statistical properties of the spatial distribution of the background emission. BeeP produces the following three new sets of parameters for each source: (a) fits to a modified blackbody (MBB) thermal emission model of the source; (b) SED-independent source flux densities at each frequency considered; and (c) fits to an MBB model of the background in which the source is embedded. BeeP also calculates, for each source, a reliability parameter, which takes into account confusion due to the surrounding cirrus. This parameter can be used to extract sub-samples of high-frequency sources with statistically well-understood properties. We define a high-reliability subset (BeeP/base), containing 26 083 sources (54.1% of the total PCCS2 catalogue), the majority of which have no information on reliability in the PCCS2. We describe the characteristics of this specific high-quality subset of PCCS2 and its validation against other data sets, specifically for: the sub-sample of PCCS2 located in low-cirrus areas; the Planck Catalogue of Galactic Cold Clumps; the Herschel GAMA15-field catalogue; and the temperature-and spectral-index-reconstructed dust maps obtained with Planck's Generalized Needlet Internal Linear Combination method. The results of the BeeP extension of PCCS2, which are made publicly available via the Planck Legacy Archive, will enable the study of the thermal properties of well-defined samples of compact Galactic and extragalactic dusty sources.

When a lightning flash is propagating in the atmosphere it is known that especially the negative leaders emit a large number of very high frequency (VHF) radio pulses. It is thought that this is due to streamer activity at the tip of the growing negative leader. In this work, we have investigated the dependence of the strength of this VHF emission on the altitude of such emission for two lightning flashes as observed by the Low Frequency ARray (LOFAR) radio telescope. We find for these two flashes that the extracted amplitude distributions are consistent with a power-law, and that the amplitude of the radio emissions decreases very strongly with source altitude, by more than a factor of 2 from 1 km altitude up to 5 km altitude. In addition, we do not find any dependence on the extracted power-law with altitude, and that the extracted power-law slope has an average around 3, for both flashes., (© 2021 The Authors.)

The common phenomenon of lightning still harbors many secrets such as what are the conditions for lightning initiation and what is driving the discharge to propagate over several tens of kilometers through the atmosphere forming conducting ionized channels called leaders. Since lightning is an electric discharge phenomenon, there are positively and negatively charged leaders. In this work we report on measurements made with the LOFAR radio telescope, an instrument primarily build for radio-astronomy observations. It is observed that a negative leader rather suddenly changes, for a few milliseconds, into a mode where it radiates 100 times more VHF power than typical negative leaders after which it spawns a large number of more typical negative leaders. This mode occurs during the initial stage, soon after initiation, of all lightning flashes we have mapped (about 25). For some flashes this mode occurs also well after initiation and we show one case where it is triggered twice, some 100 ms apart. We postulate that this is indicative of a small (order of 5 km[Formula: see text]) high charge pocket. Lightning thus appears to be initiated exclusively in the vicinity of such a small but dense charge pocket., (© 2021. The Author(s).)

Since their introduction 22 years ago, lightning mapping arrays (LMA) have played a central role in the investigation of lightning physics. Even in recent years with the proliferation of digital interferometers and the introduction of the LOw Frequency ARray (LOFAR) radio telescope, LMAs still play an important role in lightning science. LMA networks use a simple windowing technique that records the highest pulse in either 80 μs or 10 μs fixed windows in order to apply a time-of-arrival location technique. In this work, we develop an LMA-emulator that uses lightning data recorded by LOFAR to simulate an LMA, and we use it to test three new styles of pulse windowing. We show that they produce very similar results as the more traditional LMA windowing, implying that LMA lightning mapping results are relatively independent of windowing technique. In addition, each LMA station has its GPS-conditioned clock. While the timing accuracy of GPS receivers has improved significantly over the years, they still significantly limit the timing measurements of the LMA. Recently, new time-of-arrival techniques have been introduced that can be used to self-calibrate systematic offsets between different receiving stations. Applying this calibration technique to a set of data with 32 ns uncertainty, observed by the Colorado LMA, improves the timing uncertainty to 19 ns. This technique is not limited to LMAs and could be used to help calibrate future multi-station lightning interferometers., (© 2021. The Authors. Earth and Space Science published by Wiley Periodicals LLC on behalf of American Geophysical Union.)

An analysis is presented of electric fields in thunderclouds using a recently proposed method based on measuring radio emission from extensive air shower events during thunderstorm conditions. This method can be regarded as a tomography of thunderclouds using cosmic rays as probes. The data cover the period from December 2011 till August 2014. We have developed an improved fitting procedure to be able to analyze the data. Our measurements show evidence for the main negative-charge layer near the -10° isotherm. This we have seen for a winter as well as for a summer cloud where multiple events pass through the same cloud and also the vertical component of the electric field could be reconstructed. On the day of measurement of some cosmic-ray events showing evidence for strong fields, no lightning activity was detected within 100 km distance. For the winter events, the top heights were between 5 and 6 km, while in the summer, typical top heights of 9 km were seen. Large horizontal components in excess of 70 kV/m of the electric fields are observed in the middle and top layers., (©2020. The Authors.)

We use the Low Frequency Array (LOFAR) to probe the dynamics of the stepping process of negatively charged plasma channels (negative leaders) in a lightning discharge. We observe that at each step of a leader, multiple pulses of vhf (30-80 MHz) radiation are emitted in short-duration bursts (<10  μs). This is evidence for streamer formation during corona flashes that occur with each leader step, which has not been observed before in natural lightning and it could help explain x-ray emission from lightning leaders, as x rays from laboratory leaders tend to be associated with corona flashes. Surprisingly, we find that the stepping length is very similar to what was observed near the ground, however with a stepping time that is considerably larger, which as yet is not understood. These results will help to improve lightning propagation models, and eventually lightning protection models.

Lightning is a dangerous yet poorly understood natural phenomenon. Lightning forms a network of plasma channels propagating away from the initiation point with both positively and negatively charged ends-called positive and negative leaders 1 . Negative leaders propagate in discrete steps, emitting copious radio pulses in the 30-300-megahertz frequency band 2-8 that can be remotely sensed and imaged with high spatial and temporal resolution 9-11 . Positive leaders propagate more continuously and thus emit very little high-frequency radiation 12 . Radio emission from positive leaders has nevertheless been mapped 13-15 , and exhibits a pattern that is different from that of negative leaders 11-13,16,17 . Furthermore, it has been inferred that positive leaders can become transiently disconnected from negative leaders 9,12,16,18-20 , which may lead to current pulses that both reconnect positive leaders to negative leaders 11,16,17,20-22 and cause multiple cloud-to-ground lightning events 1 . The disconnection process is thought to be due to negative differential resistance 18 , but this does not explain why the disconnections form primarily on positive leaders 22 , or why the current in cloud-to-ground lightning never goes to zero 23 . Indeed, it is still not understood how positive leaders emit radio-frequency radiation or why they behave differently from negative leaders. Here we report three-dimensional radio interferometric observations of lightning over the Netherlands with unprecedented spatiotemporal resolution. We find small plasma structures-which we call 'needles'-that are the dominant source of radio emission from the positive leaders. These structures appear to drain charge from the leader, and are probably the reason why positive leaders disconnect from negative ones, and why cloud-to-ground lightning connects to the ground multiple times.

Lightning mapping technology has proven instrumental in understanding lightning. In this work we present a pipeline that can use lightning observed by the LOw-Frequency ARray (LOFAR) radio telescope to construct a 3-D map of the flash. We show that LOFAR has unparalleled precision, on the order of meters, even for lightning flashes that are over 20 km outside the area enclosed by LOFAR antennas (∼3,200 km 2 ), and can potentially locate over 10,000 sources per lightning flash. We also show that LOFAR is the first lightning mapping system that is sensitive to the spatial structure of the electrical current during individual lightning leader steps.

Cosmic rays are the highest-energy particles found in nature. Measurements of the mass composition of cosmic rays with energies of 10(17)-10(18) electronvolts are essential to understanding whether they have galactic or extragalactic sources. It has also been proposed that the astrophysical neutrino signal comes from accelerators capable of producing cosmic rays of these energies. Cosmic rays initiate air showers--cascades of secondary particles in the atmosphere-and their masses can be inferred from measurements of the atmospheric depth of the shower maximum (Xmax; the depth of the air shower when it contains the most particles) or of the composition of shower particles reaching the ground. Current measurements have either high uncertainty, or a low duty cycle and a high energy threshold. Radio detection of cosmic rays is a rapidly developing technique for determining Xmax (refs 10, 11) with a duty cycle of, in principle, nearly 100 per cent. The radiation is generated by the separation of relativistic electrons and positrons in the geomagnetic field and a negative charge excess in the shower front. Here we report radio measurements of Xmax with a mean uncertainty of 16 grams per square centimetre for air showers initiated by cosmic rays with energies of 10(17)-10(17.5) electronvolts. This high resolution in Xmax enables us to determine the mass spectrum of the cosmic rays: we find a mixed composition, with a light-mass fraction (protons and helium nuclei) of about 80 per cent. Unless, contrary to current expectations, the extragalactic component of cosmic rays contributes substantially to the total flux below 10(17.5) electronvolts, our measurements indicate the existence of an additional galactic component, to account for the light composition that we measured in the 10(17)-10(17.5) electronvolt range.

We present measurements of radio emission from cosmic ray air showers that took place during thunderstorms. The intensity and polarization patterns of these air showers are radically different from those measured during fair-weather conditions. With the use of a simple two-layer model for the atmospheric electric field, these patterns can be well reproduced by state-of-the-art simulation codes. This in turn provides a novel way to study atmospheric electric fields.

We report the results of a joint analysis of data from BICEP2/Keck Array and Planck. BICEP2 and Keck Array have observed the same approximately 400  deg^{2} patch of sky centered on RA 0 h, Dec. -57.5°. The combined maps reach a depth of 57 nK deg in Stokes Q and U in a band centered at 150 GHz. Planck has observed the full sky in polarization at seven frequencies from 30 to 353 GHz, but much less deeply in any given region (1.2  μK deg in Q and U at 143 GHz). We detect 150×353 cross-correlation in B modes at high significance. We fit the single- and cross-frequency power spectra at frequencies ≥150  GHz to a lensed-ΛCDM model that includes dust and a possible contribution from inflationary gravitational waves (as parametrized by the tensor-to-scalar ratio r), using a prior on the frequency spectral behavior of polarized dust emission from previous Planck analysis of other regions of the sky. We find strong evidence for dust and no statistically significant evidence for tensor modes. We probe various model variations and extensions, including adding a synchrotron component in combination with lower frequency data, and find that these make little difference to the r constraint. Finally, we present an alternative analysis which is similar to a map-based cleaning of the dust contribution, and show that this gives similar constraints. The final result is expressed as a likelihood curve for r, and yields an upper limit r&#95;{0.05}<0.12 at 95% confidence. Marginalizing over dust and r, lensing B modes are detected at 7.0σ significance.