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.

NuRadioMC is a Monte Carlo framework designed to simulate ultra-high energy neutrino detectors that rely on the radio detection method. This method exploits the radio emission generated in the electromagnetic component of a particle shower following a neutrino interaction. NuRadioMC simulates everything from the neutrino interaction in a medium, the subsequent Askaryan radio emission, the propagation of the radio signal to the detector and finally the detector response. NuRadioMC is designed as a modern, modular Python-based framework, combining flexibility in detector design with user-friendliness. It includes a state-of-the-art event generator, an improved modelling of the radio emission, a revisited approach to signal propagation and increased flexibility and precision in the detector simulation. This paper focuses on the implemented physics processes and their implications for detector design. A variety of models and parameterizations for the radio emission of neutrino-induced showers are compared and reviewed. Comprehensive examples are used to discuss the capabilities of the code and different aspects of instrumental design decisions.

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 past decade has welcomed the emergence of cosmic neutrinos as a new messenger to explore the most extreme environments of the universe. The discovery measurement of cosmic neutrinos, announced by IceCube in 2013, has opened a new window of observation that has already resulted in new fundamental information that holds the potential to answer key questions associated with the high-energy universe, including: what are the sources in the PeV sky and how do they drive particle acceleration; where are cosmic rays of extreme energies produced, and on which paths do they propagate through the universe; and are there signatures of new physics at TeV-PeV energies and above? The planned advancements in neutrino telescope arrays in the next decade, in conjunction with continued progress in broad multimessenger astrophysics, promise to elevate the cosmic neutrino field from the discovery to the precision era and to a survey of the sources in the neutrino sky. The planned detector upgrades to the IceCube Neutrino Observatory, culminating in IceCube-Gen2 (an envisaged $400M facility with anticipated operation in the next decade, described in this white paper) are the cornerstone that will drive the evolution of neutrino astrophysics measurements.

We report the observation of radar echoes from the ionization trails of high-energy particle cascades. These data were taken at the SLAC National Accelerator Laboratory, where the full electron beam ($\sim$10$^9$ e$^-$ at $\sim$10 GeV/e$^-$) was directed into a plastic target to simulate an ultra high-energy neutrino interaction. This target was interrogated with radio waves, and coherent radio reflections from the cascades were detected, with properties consistent with theoretical expectations. This is the first definitive observation of radar echoes from high-energy particle cascades, which may lead to a viable neutrino detection technology for energies $\gtrsim 10^{16}$ eV.

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.

RNO is the mid-scale discovery instrument designed to make the first observation of neutrinos from the cosmos at extreme energies, with sensitivity well beyond current instrument capabilities. This new observatory will be the largest ground-based neutrino telescope to date, enabling the measurement of neutrinos above $10^{16}$ eV, determining the nature of the astrophysical neutrino flux that has been measured by IceCube at higher energies, similarly extending the reach of multi-messenger astrophysics to the highest energies, and enabling investigations of fundamental physics at energies unreachable by particle accelerators on Earth.

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.

The Askaryan Radio Array (ARA) is an ultra-high energy (UHE) neutrino telescope at the South Pole consisting of an array of radio antennas aimed at detecting the Askaryan radiation produced by neutrino interactions in the ice. Currently, the experiment has five stations in operation that have been deployed in stages since 2012. This contribution focuses on the development of a search for a diffuse flux of neutrinos in two ARA stations (A2 and A3) from 2013-2016. A background of $\sim 0.01-0.02$ events is expected in one station in each of two search channels in horizontal- and vertical-polarizations. The expected new constraints on the flux of ultra-high energy neutrinos based on four years of analysis with two stations improve on the previous limits set by ARA by a factor of about two. The projected sensitivity of ARA's five-station dataset is beginning to be competitive with other neutrino telescopes at high energies near $10^{10.5}\,$GeV.

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.

Testbeam experiment 576 (T576) at the SLAC National Accelerator Laboratory sought to make the first measurement of coherent radio reflections from the ionization produced in the wake of a high-energy particle shower. The >10 GeV electron beam at SLAC End Station A was directed into a large high-density polyethylene target to produce a shower analogous to that produced by an EeV neutrino interaction in ice. Continuous wave (CW) radio was transmitted into the target, and receiving antennas monitored for reflection of the transmitted signal from the ionization left in the wake of the shower. We detail the experiment and report on preliminary hints of a signal consistent with a radio reflection at a significance of 2.36$\sigma$. We recommend another test-beam measurement in order to verify the signal.

Ultra-high energy neutrinos are detectable through impulsive radio signals generated through interactions in dense media, such as ice. Subsurface in-ice radio arrays are a promising way to advance the observation and measurement of astrophysical high-energy neutrinos with energies above those discovered by the IceCube detector ($\geq$1 PeV) as well as cosmogenic neutrinos created in the GZK process ($\geq$100 PeV). Here we describe the $\textit{NuPhase}$ detector, which is a compact receiving array of low-gain antennas deployed 185 m deep in glacial ice near the South Pole. Signals from the antennas are digitized and coherently summed into multiple beams to form a low-threshold interferometric phased array trigger for radio impulses. The NuPhase detector was installed at an Askaryan Radio Array (ARA) station during the 2017/18 Austral summer season. $\textit{In situ}$ measurements with an impulsive, point-source calibration instrument show a 50% trigger efficiency on impulses with voltage signal-to-noise ratios (SNR) of $\le$2.0, a factor of $\sim$1.8 improvement in SNR over the standard ARA combinatoric trigger. Hardware-level simulations, validated with $\textit{in situ}$ measurements, predict a trigger threshold of an SNR as low as 1.6 for neutrino interactions that are in the far field of the array. With the already-achieved NuPhase trigger performance included in ARASim, a detector simulation for the ARA experiment, we find the trigger-level effective detector volume is increased by a factor of 1.8 at neutrino energies between 10 and 100 PeV compared to the currently used ARA combinatoric trigger. We also discuss an achievable near term path toward lowering the trigger threshold further to an SNR of 1.0, which would increase the effective single-station volume by more than a factor of 3 in the same range of neutrino energies.

Ongoing experimental efforts in Antarctica seek to detect ultra-high energy neutrinos by measurement of radio-frequency (RF) Askaryan radiation generated by the collision of a neutrino with an ice molecule. An array of RF antennas, deployed either in-ice or in-air, is used to infer the properties of the neutrino. To evaluate their experimental sensitivity, such experiments require a refractive index model for ray tracing radio-wave trajectories from a putative in-ice neutrino interaction point to the receiving antennas; this gives the degree of signal absorption or ray bending from source to receiver. The gradient in the density profile over the upper 200 meters of Antarctic ice, coupled with Fermat's least-time principle, implies ray "bending" and the existence of "forbidden" zones for predominantly horizontal signal propagation at shallow depths. After re-deriving the formulas describing such shadowing, we report on experimental results that, somewhat unexpectedly, demonstrate the existence of electromagnetic wave transport modes from nominally shadowed regions. The fact that this shadow-signal propagation is observed both at South Pole and the Ross Ice Shelf in Antarctica suggests that the effect may be a generic property of polar ice, with potentially important implications for experiments seeking to detect neutrinos.

The Askaryan Radio Array (ARA) reports an observation of radio emission coincident with the "Valentine's Day" solar flare on Feb. 15$^{\rm{th}}$, 2011 in the prototype "Testbed" station. We find $\sim2000$ events that passed our neutrino search criteria during the 70 minute period of the flare, all of which reconstruct to the location of the sun. A signal analysis of the events reveals them to be consistent with that of bright thermal noise correlated across antennas. This is the first natural source of radio emission reported by ARA that is tightly reconstructable on an event-by-event basis. The observation is also the first for ARA to point radio from individual events to an extraterrestrial source on the sky. We comment on how the solar flares, coupled with improved systematic uncertainties in reconstruction algorithms, could aid in a mapping of any above-ice radio emission, such as that from cosmic-ray air showers, to astronomical locations on the sky.