The sources of ultra-high energy (UHE) cosmic rays, which can have energies up to 10^20 eV, remain a mystery. UHE neutrinos may provide important clues to understanding the nature of cosmic-ray sources. ARIANNA aims to detect UHE neutrinos via radio (Askaryan) emission from particle showers when a neutrino interacts with ice, which is an efficient method for neutrinos with energies between 10^16 eV and 10^20 eV. The ARIANNA radio detectors are located in Antarctic ice just beneath the surface. Neutrino observation requires that radio pulses propagate to the antennas at the surface with minimum distortion by the ice and firn medium. Using the residual hole from the South Pole Ice Core Project, radio pulses were emitted from a transmitter located up to 1.7 km below the snow surface. By measuring these signals with an ARIANNA surface station, the angular and polarization reconstruction abilities are quantified, which are required to measure the direction of the neutrino. After deconvolving the raw signals for the detector response and attenuation from propagation through the ice, the signal pulses show no significant distortion and agree with a reference measurement of the emitter made in an anechoic chamber. Furthermore, the signal pulses reveal no significant birefringence for our tested geometry of mostly vertical ice propagation. The origin of the transmitted radio pulse was measured with an angular resolution of 0.37 degrees indicating that the neutrino direction can be determined with good precision if the polarization of the radio-pulse can be well determined. In the present study we obtained a resolution of the polarization vector of 2.7 degrees. Neither measurement show a significant offset relative to expectation.

The proposed ARIANNA-200 neutrino detector, located at sea-level on the Ross Ice Shelf, Antarctica, consists of 200 autonomous and independent detector stations separated by 1 kilometer in a uniform triangular mesh, and serves as a pathfinder mission for the future IceCube-Gen2 project. The primary science mission of ARIANNA-200 is to search for sources of neutrinos with energies greater than 10^17 eV, complementing the reach of IceCube. An ARIANNA observation of a neutrino source would provide strong insight into the enigmatic sources of cosmic rays. ARIANNA observes the radio emission from high energy neutrino interactions in the Antarctic ice. Among radio based concepts under current investigation, ARIANNA-200 would uniquely survey the vast majority of the southern sky at any instant in time, and an important region of the northern sky, by virtue of its location on the surface of the Ross Ice Shelf in Antarctica. The broad sky coverage is specific to the Moore's Bay site, and makes ARIANNA-200 ideally suited to contribute to the multi-messenger thrust by the US National Science Foundation, Windows on the Universe - Multi-Messenger Astrophysics, providing capabilities to observe explosive sources from unknown directions. The ARIANNA architecture is designed to measure the angular direction to within 3 degrees for every neutrino candidate, which too plays an important role in the pursuit of multi-messenger observations of astrophysical sources.

The primary mission of the ARIANNA ultra-high energy neutrino telescope is to uncover astrophysical sources of neutrinos with energies greater than $10^{16}\mathrm{eV}$. A pilot array, consisting of seven ARIANNA stations located on the surface of the Ross Ice Shelf in Antarctica, was commissioned in November 2014. We report on the search for astrophysical neutrinos using data collected between November 2014 and February 2019. A straight-forward template matching analysis yielded no neutrino candidates, with a signal efficiency of 79%. We find a 90% confidence upper limit on the diffuse neutrino flux of $E^2\Phi=1.7\times 10^{-6}\mathrm{GeV cm^{-2}s^{-1}sr^{-1}}$ for a decade wide logarithmic bin centered at a neutrino energy of $10^{18}\mathrm{eV}$, which is an order of magnitude improvement compared to the previous limit reported by the ARIANNA collaboration. The ARIANNA stations, including purpose built cosmic-ray stations at the Moore's Bay site and demonstrator stations at the South Pole, have operated reliably. Sustained operation at two distinct sites confirms that the flexible and adaptable architecture can be deployed in any deep ice, radio quiet environment. We show that the scientific capabilities, technical innovations, and logistical requirements of ARIANNA are sufficiently well understood to serve as the basis for large area radio-based neutrino telescope with a wide field-of-view.

The Askaryan Radio Array is an ultra-high energy (UHE, $>10^{17}$\,eV) neutrino detector designed to observe neutrinos by searching for the radio waves emitted by the relativistic byproducts of neutrino-nucleon interactions in Antarctic ice. In this paper, we present constraints on the diffuse flux of ultra-high energy neutrinos between $10^{16}-10^{21}$\,eV resulting from a search for neutrinos in two complementary analyses, both analyzing four years of data (2013-2016) from the two deep stations (A2, A3) operating at that time. We place a 90\% CL upper limit on the diffuse neutrino flux at $10^{18}$ eV of $EF(E)=3.9\times10^{-16}$~$\textrm{cm}^2$$\textrm{s}^{-1}$$\textrm{sr}^{-1}$. The analysis leverages more than quadruple the livetime of the previous ARA limit, and at this energy is the most sensitive reported by an in-ice radio neutrino detector by more than a factor of four. Looking forward, ARA has yet to analyze more than half of its archived data, and expects that dataset to double in size in the next three years as livetime accumulates in the full five-station array. With this additional livetime, we project ARA will have world-leading sensitivity to neutrinos above $10^{17}$ eV by 2022.

Ultra high energy neutrinos ($E_\nu > 10^{16.5}$eV$)$ are efficiently measured via radio signals following a neutrino interaction in ice. An antenna placed $\mathcal{O}$(15 m) below the ice surface will measure two signals for the vast majority of events (90% at $E_\nu$=$10^{18}$eV$)$: a direct pulse and a second delayed pulse from a reflection off the ice surface. This allows for a unique identification of neutrinos against backgrounds arriving from above. Furthermore, the time delay between the direct and reflected signal (D'n'R) correlates with the distance to the neutrino interaction vertex, a crucial quantity to determine the neutrino energy. In a simulation study, we derive the relation between time delay and distance and study the corresponding experimental uncertainties in estimating neutrino energies. We find that the resulting contribution to the energy resolution is well below the natural limit set by the unknown inelasticity in the initial neutrino interaction. We present an in-situ measurement that proves the experimental feasibility of this technique. Continuous monitoring of the local snow accumulation in the vicinity of the transmit and receive antennas using this technique provide a precision of $\mathcal{O}$(1 mm) in surface elevation, which is much better than that needed to apply the D'n'R technique to neutrinos.

The Askaryan Radio Array (ARA) experiment at the South Pole is designed to detect high-energy neutrinos which, via in-ice interactions, produce coherent radiation at frequencies up to 1000 MHz. In Dec. 2018, a custom high-amplitude radio-frequency transmitter was lowered into the 1700 m SPICE ice core to provide test sources for ARA receiver stations sensitive to vertical and horizontal polarizations. For these tests, signal geometries correspond to obliquely propagating radio waves from below. The ARA collaboration has recently measured the polarization-dependent time delay variation, and report more significant time delays for trajectories perpendicular to ice flow. Here we use fabric data from the SPICE ice core to construct a bounding model for the ice birefringence and the polarization time delays across ARA. The data-model comparison is consistent with the vertical girdle fabric at the South Pole having the prevailing horizontal crystallographic axis oriented near-perpendicular to ice flow. This study presents the possibility that ice birefringence can be used to constrain the range to a neutrino interaction, and hence aid in neutrino energy reconstruction, for in-ice experiments such as ARA.

We report on analysis of englacial radio-frequency (RF) pulser data received over horizontal baselines of 1--5 km, based on broadcasts from two sets of transmitters deployed to depths of up to 1500 meters at the South Pole. First, we analyze data collected usingtwo RF bicone transmitters 1400 meters below the ice surface, and frozen into boreholes drilled for the IceCube experiment in 2011. Additionally, in Dec., 2018, a fat-dipole antenna, fed by one of three high-voltage (~1 kV), fast (~(1-5 ns)) signal generators was lowered into the 1700-m deep icehole drilled for the South Pole Ice Core Experiment (SPICE), approximately 3 km from the geographic South Pole. Signals from transmitters were recorded on the five englacial multi-receiver ARA stations, with receiver depths between 60--200 m. We confirm the long, >1 km RF electric field attenuation length, test our observed signal arrival timing distributions against models, and measure birefringent asymmetries at the 0.15% level.

A Monte Carlo simulation program for the radio detection of Ultra High Energy (UHE) neutrino interactions in the Antarctic ice as viewed by the Antarctic Impulsive Transient Antenna (ANITA) is described in this article. The program, icemc, provides an input spectrum of UHE neutrinos, the parametrization of the Askaryan radiation generated by their interaction in the ice, and the propagation of the radiation through ice and air to a simulated model of the third and fourth ANITA flights. This paper provides an overview of the icemc simulation, descriptions of the physics models used and of the ANITA electronics processing chain, data/simulation comparisons to validate the predicted performance, and a summary of the impact of published results.

The measurement of ultra-high energy (UHE) neutrinos (E $>$ \SI{e16}{eV}) opens a new field of astronomy with the potential to reveal the sources of ultra-high energy cosmic rays especially if combined with observations in the electromagnetic spectrum and gravitational waves. The ARIANNA pilot detector explores the detection of UHE neutrinos with a surface array of independent radio detector stations in Antarctica which allows for a cost-effective instrumentation of large volumes. Twelve stations are currently operating successfully at the Moore's Bay site (Ross Ice Shelf) in Antarctica and at the South Pole. We will review the current state of ARIANNA and its main results. We report on a newly developed wind generator that successfully operates in the harsh Antarctic conditions and powers the station for a substantial time during the dark winter months. The robust ARIANNA surface architecture, combined with environmentally friendly solar and wind power generators, can be installed at any deep ice location on the planet and operated autonomously. We report on the detector capabilities to determine the neutrino direction by reconstructing the signal arrival direction of a \SI{800}{m} deep calibration pulser, and the reconstruction of the signal polarization using the more abundant cosmic-ray air showers. Finally, we describe a large-scale design -- ARIA -- that capitalizes on the successful experience of the ARIANNA operation and is designed sensitive enough to discover the first UHE neutrino.

The ANtarctic Impulsive Transient Antenna (ANITA) NASA long-duration balloon payload completed its fourth flight in December 2016, after 28 days of flight time. ANITA is sensitive to impulsive broadband radio emission from interactions of ultra-high-energy neutrinos in polar ice (Askaryan emission). We present the results of two separate blind analyses searching for signals from Askaryan emission in the data from the fourth flight of ANITA. The more sensitive analysis, with a better expected limit, has a background estimate of $0.64^{+0.69}_{-0.45}$ and an analysis efficiency of $82\pm2\%$. The second analysis has a background estimate of $0.34^{+0.66}_{-0.16}$ and an analysis efficiency of $71\pm6\%$. Each analysis found one event in the signal region, consistent with the background estimate for each analysis. The resulting limit further tightens the constraints on the diffuse flux of ultra-high-energy neutrinos at energies above $10^{19.5}$ eV.

Recently, the ANITA collaboration reported on two upward-going extensive air shower events consistent with a primary particle that emerges from the surface of the ice. These events may be of $\nu_\tau$ origin, in which the neutrino interacts within the Earth to produce a $\tau$ lepton that emerges from the Earth, decays in the atmosphere, and initiates an extensive air shower. In this paper we estimate an upper bound on the ANITA acceptance to a diffuse $\nu_\tau$ flux detected via $\tau$-lepton-induced air showers within the bounds of Standard Model (SM) uncertainties. By comparing this estimate with the acceptance of Pierre Auger Observatory and IceCube and assuming SM interactions, we conclude that a $\nu_\tau$ origin of these events would imply a neutrino flux at least two orders of magnitude above current bounds.

Using data collected by the Askaryan Radio Array (ARA) experiment at the South Pole, we have used long-baseline propagation of radio-frequency signals to extract information on the radio-frequency index-of-refraction in South Polar ice. Owing to the increasing ice density over the upper 150--200 meters, rays are observed along two, nearly parallel paths, one of which is direct and a second which refracts through an inflection point, with differences in both arrival time and arrival angle that can be used to constrain the neutrino properties. We also observe indications, for the first time, of radio-frequency ice birefringence for signals propagating along predominantly horizontal trajectories, corresponding to an asymmetry of order 0.1% between the ordinary and extra-ordinary paths, numerically compatible with previous measurements of birefringent asymmetries for vertically-propagating radio-frequency signals at South Pole. Taken together, these effects offer the possibility of redundantly measuring the range from receiver to a neutrino interaction in Antarctic ice, if receiver antennas are deployed at shallow (25 m<z<100 m) depths. Such range information is essential in determining both the neutrino energy, as well as the incident neutrino direction.

These proceedings address a recent publication by the ANITA collaboration of four upward- pointing cosmic-ray-like events observed in the first flight of ANITA. Three of these events were consistent with stratospheric cosmic-ray air showers where the axis of propagation does not inter- sect the surface of the Earth. The fourth event was consistent with a primary particle that emerges from the surface of the ice suggesting a possible {\tau}-lepton decay as the origin of this event. These proceedings follow-up on the modeling and testing of the hypothesis that this event was of {\tau} neutrino origin.

The Antarctic Impulsive Transient Antenna (ANITA), a NASA long-duration balloon payload, searches for radio emission from interactions of ultra-high-energy neutrinos in polar ice. The third flight of ANITA (ANITA-III) was launched in December 2014 and completed a 22-day flight. We present the results of three analyses searching for Askaryan radio emission of neutrino origin. In the most sensitive of the analyses, we find one event in the signal region on an expected a priori background of $0.7^{+0.5}_{-0.3}$. Though consistent with the background estimate, the candidate event remains compatible with a neutrino hypothesis even after additional post-unblinding scrutiny.

We report on an upward traveling, radio-detected cosmic-ray-like impulsive event with characteristics closely matching an extensive air shower. This event, observed in the third flight of the Antarctic Impulsive Transient Antenna (ANITA), a NASA-sponsored long-duration balloon payload, is consistent with a similar event reported in a previous flight. These events may be produced by the atmospheric decay of an upward-propagating $\tau$-lepton produced by a $\nu_{\tau}$ interaction, although their relatively steep arrival angles create tension with the standard model (SM) neutrino cross section. Each of the two events have $a~posteriori$ background estimates of $\lesssim 10^{-2}$ events. If these are generated by $\tau$-lepton decay, then either the charged-current $\nu_{\tau}$ cross section is suppressed at EeV energies, or the events arise at moments when the peak flux of a transient neutrino source was much larger than the typical expected cosmogenic background neutrinos.

The balloon-borne HiCal radio-frequency (RF) transmitter, in concert with the ANITA radio-frequency receiver array, is designed to measure the Antarctic surface reflectivity in the RF wavelength regime. The amplitude of surface-reflected transmissions from HiCal, registered as triggered events by ANITA, can be compared with the direct transmissions preceding them by O(10) microseconds, to infer the surface power reflection coefficient $\cal{R}$. The first HiCal mission (HiCal-1, Jan. 2015) yielded a sample of 100 such pairs, resulting in estimates of $\cal{R}$ at highly-glancing angles (i.e., zenith angles approaching $90^\circ$), with measured reflectivity for those events which exceeded extant calculations. The HiCal-2 experiment, flying from Dec., 2016-Jan., 2017, provided an improvement by nearly two orders of magnitude in our event statistics, allowing a considerably more precise mapping of the reflectivity over a wider range of incidence angles. We find general agreement between the HiCal-2 reflectivity results and those obtained with the earlier HiCal-1 mission, as well as estimates from Solar reflections in the radio-frequency regime. In parallel, our calculations of expected reflectivity have matured; herein, we use a plane-wave expansion to estimate the reflectivity R from both a flat, smooth surface (and, in so doing, recover the Fresnel reflectivity equations) and also a curved surface. Multiplying our flat-smooth reflectivity by improved Earth curvature and surface roughness corrections now provides significantly better agreement between theory and the HiCal 2a/2b measurements.

The Antarctic Impulsive Transient Antenna (ANITA) is a NASA long-duration balloon experiment with the primary goal of detecting ultra-high-energy ($>10^{18}\,\mbox{eV}$) neutrinos via the Askaryan Effect. The fourth ANITA mission, ANITA-IV, recently flew from Dec 2 to Dec 29, 2016. For the first time, the Tunable Universal Filter Frontend (TUFF) boards were deployed for mitigation of narrow-band, anthropogenic noise with tunable, switchable notch filters. The TUFF boards also performed second-stage amplification by approximately 45 dB to boost the $\sim\,\mu\mbox{V-level}$ radio frequency (RF) signals to $\sim$ mV-level for digitization, and supplied power via bias tees to the first-stage, antenna-mounted amplifiers. The other major change in signal processing in ANITA-IV is the resurrection of the $90^{\circ}$ hybrids deployed previously in ANITA-I, in the trigger system, although in this paper we focus on the TUFF boards. During the ANITA-IV mission, the TUFF boards were successfully operated throughout the flight. They contributed to a factor of 2.8 higher total instrument livetime on average in ANITA-IV compared to ANITA-III due to reduction of narrow-band, anthropogenic noise before a trigger decision is made.

The primary science goal of the NASA-sponsored ANITA project is measurement of ultra-high energy neutrinos and cosmic rays, observed via radio-frequency signals resulting from a neutrino- or cosmic ray- interaction with terrestrial matter (atmospheric or ice molecules, e.g.). Accurate inference of the energies of these cosmic rays requires understanding the transmission/reflection of radio wave signals across the ice-air boundary. Satellite-based measurements of Antarctic surface reflectivity, using a co-located transmitter and receiver, have been performed more-or-less continuously for the last few decades. Satellite-based reflectivity surveys, at frequencies ranging from 2--45 GHz and at near-normal incidence, yield generally consistent reflectivity maps across Antarctica. Using the Sun as an RF source, and the ANITA-3 balloon borne radio-frequency antenna array as the RF receiver, we have also measured the surface reflectivity over the interval 200-1000 MHz, at elevation angles of 12-30 degrees, finding agreement with the Fresnel equations within systematic errors. To probe low incidence angles, inaccessible to the Antarctic Solar technique and not probed by previous satellite surveys, a novel experimental approach ("HiCal-1") was devised. Unlike previous measurements, HiCal-ANITA constitute a bi-static transmitter-receiver pair separated by hundreds of kilometers. Data taken with HiCal, between 200--600 MHz shows a significant departure from the Fresnel equations, constant with frequency over that band, with the deficit increasing with obliquity of incidence, which we attribute to the combined effects of possible surface roughness, surface grain effects, radar clutter and/or shadowing of the reflection zone due to Earth curvature effects.

Ultra-high energy neutrinos are interesting messenger particles since, if detected, they can transmit exclusive information about ultra-high energy processes in the Universe. These particles, with energies above $10^{16}\mathrm{eV}$, interact very rarely. Therefore, detectors that instrument several gigatons of matter are needed to discover them. The ARA detector is currently being constructed at South Pole. It is designed to use the Askaryan effect, the emission of radio waves from neutrino-induced cascades in the South Pole ice, to detect neutrino interactions at very high energies. With antennas distributed among 37 widely-separated stations in the ice, such interactions can be observed in a volume of several hundred cubic kilometers. Currently 3 deep ARA stations are deployed in the ice of which two have been taking data since the beginning of the year 2013. In this publication, the ARA detector "as-built" and calibrations are described. Furthermore, the data reduction methods used to distinguish the rare radio signals from overwhelming backgrounds of thermal and anthropogenic origin are presented. Using data from only two stations over a short exposure time of 10 months, a neutrino flux limit of $3 \cdot 10^{-6} \mathrm{GeV} / (\mathrm{cm^2 \ s \ sr})$ is calculated for a particle energy of 10^{18}eV, which offers promise for the full ARA detector.

Employing a 40-kW radio-frequency transmitter just west of Delta, UT, and operating at 54.1 MHz, the TARA (Telescope Array RAdar) experiment seeks radar detection of extensive air showers (EAS) initiated by ultra-high energy cosmic rays (UHECR). For UHECR with energies in excess of $10^{19}$ eV, the Doppler-shifted "chirps" resulting from EAS shower core radar reflections should be observable above background (dominantly galactic) at distances of tens of km from the TARA transmitter. In order to stereoscopically reconstruct cosmic ray chirps, two remote, autonomous self-powered receiver stations have been deployed. Each remote station (RS) combines both low power consumption as well as low cost. Triggering logic, the powering and communication systems, and some specific details of hardware components are discussed.

VLBI observations of the Crab pulsar with the 64-m radio telescope at Kalyazin (Russia) and the 46-m radio telescope of the Algonquin Radio Observatory (Canada) at 2.2 GHz and single-dish observations of the millisecond pulsar B1937+21 with the GBT (USA) at 2.1 GHz were conducted to probe the interstellar medium and study the properties of giant pulses. The VLBI data were processed with a dedicated software correlator, which allowed us to obtain the visibility of single giant pulses. Two frequency scales of 50 and 450 kHz were found in the diffraction spectra of giant pulses from the Crab pulsar. The location of the scattering region was estimated to be close to the outer edge of the nebula. No correlation was found between the power spectra of giant pulses at left- and right-hand circular polarization. We explain this lack of correlation through the influence of the strong magnetic field on circularly polarized emission in the region close to the Crab pulsar. Combining the measurement of the decorrelation bandwidth with that of the scattering time of giant pulses for B1937+21, we found three frequency scales of 1.7, 3.8, and 16.5 MHz. The scattering time of giant pulses of B1937+21 at 2.1 GHz was found to be 40+-4 ns. We obtained an upper limit of the intrinsic width of giant pulses from B1937+21 of less than 8 ns. The frequency dependences of the scattering times for the Crab pulsar and PSR B1937+21 were found to be different. They are characterized by exponents of -3.5 and -4.2, respectively. We attribute the difference to the large influence of scattering in the Crab nebula.

Giant radio pulses of the millisecond pulsar B1937+21 were recorded with the S2 VLBI system at 1.65 GHz with NASA/JPL's 70-m radio telescope at Tidbinbilla, Australia. These pulses have been observed as strong as 65000 Jy with widths <= 15 ns, corresponding to a brightness temperature T_b >= 5 x 10^39 K, the highest observed in the universe. The vast majority of these pulses occur in a 5.8 mcs and 8.2 mcs window at the very trailing edges of the regular main pulse and interpulse profiles, respectively. Giant pulses occur in general with a single spike. Only in one case out of 309 was the structure clearly more complex. The cumulative distribution is fit by a power law with index -1.40 +/- 0.01 with a low-energy but no high-energy cutoff. We estimate that giant pulses occur frequently but are only rarely detected. When corrected for the directivity factor, 25 giant pulses are estimated to be generated in one neutron star revolution alone. The intensities of the giant pulses of the main pulses and interpulses are not correlated with each other nor with the intensities or energies of the main pulses and interpulses themselves. Their radiation energy density can exceed 300 times the plasma energy density at the surface of the neutron star and can even exceed the magnetic field energy density at that surface. We therefore do not think that the generation of giant pulses is linked to the plasma mechanisms in the magnetosphere. Instead we suggest that it is directly related to discharges in the polar cap region of the pulsar.

We report on a study of microstructure and its quasi-periodicities of three pulsars at 1.65 GHz with the S2 VLBI system at a resolution of 62.5 ns, by far the highest for any such statistical study yet. For PSR B1929+10 we found in the average cross-correlation function (CCF) broad microstructure with a characteristic timescale of 95+-10 mcs and confirmed microstructure with characteristic timescales between 100 and 450 mcs for PSRs B0950+08 and B1133+16. On a finer scale PSRs B0950+08, B1133+16 (component II) and B1929+10 show narrow microstructure with a characteristic timescale in the CCFs of ~10 mcs, the shortest found in the average CCF or autocorrelation function (ACF) for any pulsar, apart perhaps for the Crab pulsar. Histograms of microstructure widths are skewed heavily toward shorter timescales but display a sharp cutoff. The shortest micropulses have widths between 2 and 7 mcs. No nanopulses or unresolved pulse spikes were detected. Cross-power spectra of single pulses show a large range of complexity with single spectral features representing classic quasi-periodicities and broad and overlapping features with essentially no periodicities at all. Significant differences were found for the two components of PSR B1133+16 in every aspect of our statistical analysis of micropulses and their quasi-periodicities. Asymmetries in the magnetosphere and the hollow cone of emission above the polar cap of the neutron star may be responsible for these differences.