Hostname: page-component-788cddb947-wgjn4 Total loading time: 0 Render date: 2024-10-12T01:33:54.075Z Has data issue: false hasContentIssue false

In situ, broadband measurement of the radio frequency attenuation length at Summit Station, Greenland

Published online by Cambridge University Press:  30 May 2022

J. A. Aguilar
Affiliation:
Université Libre de Bruxelles, Science Faculty CP230, B-1050 Brussels, Belgium
P. Allison
Affiliation:
Department of Physics, Center for Cosmology and AstroParticle Physics, Ohio State University, Columbus, OH 43210, USA
J. J. Beatty
Affiliation:
Department of Physics, Center for Cosmology and AstroParticle Physics, Ohio State University, Columbus, OH 43210, USA
D. Besson
Affiliation:
Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045, USA National Nuclear Research University MEPhI, Kashirskoe Shosse 31, 115409 Moscow, Russia
A. Bishop
Affiliation:
Wisconsin IceCube Particle Astrophysics Center (WIPAC) and Department of Physics, University of Wisconsin-Madison, Madison, WI 53703, USA
O. Botner
Affiliation:
Department of Physics and Astronomy, Uppsala University, Uppsala, SE-752 37 Sweden
S. Bouma
Affiliation:
Erlangen Center for Astroparticle Physics (ECAP), Friedrich-Alexander-University Erlangen-Nuremberg, 91058 Erlangen, Germany
S. Buitink
Affiliation:
Astrophysical Institute, Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium
M. Cataldo
Affiliation:
Erlangen Center for Astroparticle Physics (ECAP), Friedrich-Alexander-University Erlangen-Nuremberg, 91058 Erlangen, Germany
B. A. Clark
Affiliation:
Department of Physics and Astronomy, Michigan State University, East Lansing MI 48824, USA
Z. Curtis-Ginsberg
Affiliation:
Department of Physics, Enrico Fermi Institute, Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA
A. Connolly
Affiliation:
Department of Physics, Center for Cosmology and AstroParticle Physics, Ohio State University, Columbus, OH 43210, USA
P. Dasgupta
Affiliation:
Université Libre de Bruxelles, Science Faculty CP230, B-1050 Brussels, Belgium
S. de Kockere
Affiliation:
Dienst ELEM, Vrije Universiteit Brussel, B-1050 Brussels, Belgium
K. D. de Vries
Affiliation:
Dienst ELEM, Vrije Universiteit Brussel, B-1050 Brussels, Belgium
C. Deaconu
Affiliation:
Department of Physics, Enrico Fermi Institute, Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA
M. A. DuVernois
Affiliation:
Wisconsin IceCube Particle Astrophysics Center (WIPAC) and Department of Physics, University of Wisconsin-Madison, Madison, WI 53703, USA
C. Glaser
Affiliation:
Department of Physics and Astronomy, Uppsala University, Uppsala, SE-752 37 Sweden
A. Hallgren
Affiliation:
Department of Physics and Astronomy, Uppsala University, Uppsala, SE-752 37 Sweden
S. Hallmann
Affiliation:
DESY, Platanenallee 6, 15738 Zeuthen, Germany
J. C. Hanson
Affiliation:
Whittier College, Whittier, CA 90602, USA
B. Hendricks
Affiliation:
Department of Physics, Department of Astronomy & Astrophysics, Penn State University, University Park, PA 16801, USA
C. Hornhuber
Affiliation:
Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045, USA
K. Hughes
Affiliation:
Department of Physics, Enrico Fermi Institute, Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA
A. Karle
Affiliation:
Wisconsin IceCube Particle Astrophysics Center (WIPAC) and Department of Physics, University of Wisconsin-Madison, Madison, WI 53703, USA
J. L. Kelley
Affiliation:
Wisconsin IceCube Particle Astrophysics Center (WIPAC) and Department of Physics, University of Wisconsin-Madison, Madison, WI 53703, USA
I. Kravchenko
Affiliation:
Department of Physics and Astronomy, University of Nebraska-Lincoln, NE 68588, USA
R. Krebs
Affiliation:
Department of Physics, Department of Astronomy & Astrophysics, Penn State University, University Park, PA 16801, USA
R. Lahmann
Affiliation:
Erlangen Center for Astroparticle Physics (ECAP), Friedrich-Alexander-University Erlangen-Nuremberg, 91058 Erlangen, Germany
U. Latif
Affiliation:
Dienst ELEM, Vrije Universiteit Brussel, B-1050 Brussels, Belgium
J. Mammo
Affiliation:
Department of Physics and Astronomy, University of Nebraska-Lincoln, NE 68588, USA
Z. S. Meyers
Affiliation:
Erlangen Center for Astroparticle Physics (ECAP), Friedrich-Alexander-University Erlangen-Nuremberg, 91058 Erlangen, Germany DESY, Platanenallee 6, 15738 Zeuthen, Germany
K. Michaels
Affiliation:
Department of Physics, Enrico Fermi Institute, Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA
K. Mulrey
Affiliation:
Department of Astrophysics/IMAPP, Radboud University, PO Box 9010, 6500 GL, The Netherlands
A. Nelles
Affiliation:
Erlangen Center for Astroparticle Physics (ECAP), Friedrich-Alexander-University Erlangen-Nuremberg, 91058 Erlangen, Germany DESY, Platanenallee 6, 15738 Zeuthen, Germany
A. Novikov
Affiliation:
Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045, USA
A. Nozdrina
Affiliation:
Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045, USA
E. Oberla
Affiliation:
Department of Physics, Enrico Fermi Institute, Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA
B. Oeyen
Affiliation:
Department of Physics and Astronomy, Ghent University, B-9000 Gent, Belgium
Y. Pan
Affiliation:
Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, USA
H. Pandya
Affiliation:
Astrophysical Institute, Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium
I. Plaisier
Affiliation:
Erlangen Center for Astroparticle Physics (ECAP), Friedrich-Alexander-University Erlangen-Nuremberg, 91058 Erlangen, Germany DESY, Platanenallee 6, 15738 Zeuthen, Germany
N. Punsuebsay
Affiliation:
Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, USA
L. Pyras
Affiliation:
Erlangen Center for Astroparticle Physics (ECAP), Friedrich-Alexander-University Erlangen-Nuremberg, 91058 Erlangen, Germany DESY, Platanenallee 6, 15738 Zeuthen, Germany
D. Ryckbosch
Affiliation:
Department of Physics and Astronomy, Ghent University, B-9000 Gent, Belgium
O. Scholten
Affiliation:
Dienst ELEM, Vrije Universiteit Brussel, B-1050 Brussels, Belgium Kapteijn Institute, University of Groningen, Groningen, The Netherlands
D. Seckel
Affiliation:
Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, USA
M. F. H. Seikh
Affiliation:
Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045, USA
D. Smith*
Affiliation:
Department of Physics, Enrico Fermi Institute, Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA
D. Southall
Affiliation:
Department of Physics, Enrico Fermi Institute, Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA
J. Torres
Affiliation:
Department of Physics, Center for Cosmology and AstroParticle Physics, Ohio State University, Columbus, OH 43210, USA
S. Toscano
Affiliation:
Université Libre de Bruxelles, Science Faculty CP230, B-1050 Brussels, Belgium
D. Tosi
Affiliation:
Wisconsin IceCube Particle Astrophysics Center (WIPAC) and Department of Physics, University of Wisconsin-Madison, Madison, WI 53703, USA
D. J. Van Den Broeck
Affiliation:
Astrophysical Institute, Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium Dienst ELEM, Vrije Universiteit Brussel, B-1050 Brussels, Belgium
N. van Eijndhoven
Affiliation:
Dienst ELEM, Vrije Universiteit Brussel, B-1050 Brussels, Belgium
A. G. Vieregg
Affiliation:
Department of Physics, Enrico Fermi Institute, Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA
C. Welling
Affiliation:
Erlangen Center for Astroparticle Physics (ECAP), Friedrich-Alexander-University Erlangen-Nuremberg, 91058 Erlangen, Germany DESY, Platanenallee 6, 15738 Zeuthen, Germany
D. R. Williams
Affiliation:
Department of Physics and Astronomy, University of Alabama, Tuscaloosa, AL 35487, USA
S. Wissel
Affiliation:
Department of Physics, Department of Astronomy & Astrophysics, Penn State University, University Park, PA 16801, USA Physics Department California Polytechnic State University, San Luis Obispo CA 93407, USA
R. Young
Affiliation:
Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045, USA
A. Zink
Affiliation:
Erlangen Center for Astroparticle Physics (ECAP), Friedrich-Alexander-University Erlangen-Nuremberg, 91058 Erlangen, Germany
*
Author for correspondence: D. Smith, E-mail: danielsmith@uchicago.edu
Rights & Permissions [Opens in a new window]

Abstract

Over the last 25 years, radiowave detection of neutrino-generated signals, using cold polar ice as the neutrino target, has emerged as perhaps the most promising technique for detection of extragalactic ultra-high energy neutrinos (corresponding to neutrino energies in excess of 0.01 Joules, or 1017 electron volts). During the summer of 2021 and in tandem with the initial deployment of the Radio Neutrino Observatory in Greenland (RNO-G), we conducted radioglaciological measurements at Summit Station, Greenland to refine our understanding of the ice target. We report the result of one such measurement, the radio-frequency electric field attenuation length $L_\alpha$. We find an approximately linear dependence of $L_\alpha$ on frequency with the best fit of the average field attenuation for the upper 1500 m of ice: $\langle L_\alpha \rangle = ( ( 1154 \pm 121) - ( 0.81 \pm 0.14) \, ( \nu /{\rm MHz}) ) \,{\rm m}$ for frequencies ν ∈ [145 − 350] MHz.

Type
Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © The Author(s), 2022. Published by Cambridge University Press

1. Introduction

We report a measurement of the radio frequency electric field attenuation length of deep glacial ice at the US National Science Foundation's Summit Station in Greenland. This measurement is of interest to the ultra-high energy neutrino (UHEN) community due to the development of the Radio Neutrino Observatory in Greenland (RNO-G), a particle astrophysics experiment that uses the ice as a target material in the search for astrophysical and cosmogenic neutrinos (Aguilar and others, Reference Aguilar2021).

The IceCube experiment has placed a flux upper limit for astrophysical neutrinos of $E^2 \phi \leq 2\times 10^{-8}\, {\rm GeV}/( {\rm cm}^2$ s sr) at $E_\nu = 1$ EeV (Aartsen and others, Reference Aartsen2018). At such fluxes, a particle detector requires an active volume of O (10  km3) or larger for a discovery-level detection within a detector's lifetime. A sparsely instrumented array of radio antennas, deployed in and on an extensive dielectric medium, can satisfy this volume requirement. An UHEN interaction creates an extensive electromagnetic shower that produces impulsive radio emission via the Askaryan effect (Askaryan, Reference Askaryan1961). If the interaction occurs in an environment of low radio attenuation, a relatively small number of radio antennas can probe target volumes at the scale needed for UHEN observations.

Glacial ice has been measured to have long radio attenuation lengths due to low temperature and relatively high purity (Barwick and others, Reference Barwick, Besson, Gorham and Saltzberg2005; Besson and others, Reference Besson2008; Barrella and other, Reference Barrella, Barwick and Saltzberg2011; Allison and others, Reference Allison2012; Avva and others, Reference Avva, Kovac, Miki, Saltzberg and Vieregg2015; Hanson and others, Reference Hanson2015). This, combined with the volume of glacial ice available in Greenland and Antarctica, makes polar ice sheets attractive sites for the construction of a radio neutrino detector.

RNO-G is one such experiment based on radio detection of UHEN in glacial ice (Aguilar and others, Reference Aguilar2021), among others in Antarctica (Anker and others, Reference Anker2019; Kravchenko and others, Reference Kravchenko2003; Gorham and others, Reference Gorham2009, Reference Gorham2019; Allison and others, Reference Allison2012, Reference Allison2020). RNO-G is being constructed near NSF's Summit Station at the highest point of the Greenland Ice Sheet. The planned detector will ultimately be composed of 35 autonomous stations separated by 1.25 km in a grid pattern. Each station is instrumented with radio antennas, with good response over the range 100–600 MHz, deployed both just below the surface and at depths down to 100 m in boreholes. Construction of the detector began during the summer of 2021 with the installation of the first three stations.

Previous measurements of radio attenuation lengths in Antarctica (Barwick and others, Reference Barwick, Besson, Gorham and Saltzberg2005; Besson and others, Reference Besson2008; Barrella and other, Reference Barrella, Barwick and Saltzberg2011; Allison and others, Reference Allison2012; Hanson and others, Reference Hanson2015) and Greenland (Paden and others, Reference Paden2005; Avva and others, Reference Avva, Kovac, Miki, Saltzberg and Vieregg2015; MacGregor and others, Reference MacGregor2015) have demonstrated that radio attenuation lengths vary at different ice locations, due primarily to differences in ice temperature and impurity levels. Since electric field attenuation length is a primary determinant of the expected number of observed UHENs at energies greater than 1 EeV, a precise, in situ measurement is required at Summit Station to assess RNO-G's science potential.

Our work builds upon one previous in situ measurement of the bulk ice electric field attenuation length performed at Summit Station by Avva and others (Reference Avva, Kovac, Miki, Saltzberg and Vieregg2015); that effort reported a depth-averaged attenuation length $\langle L_\alpha \rangle = 947^{ + 92}_{-85}$ m at 75 MHz. We herein quantify the attenuation length at higher frequencies to better match RNO-G's frequency range of 100–600 MHz. In addition to that prior analysis, there have been several previous measurements of the radio ice properties at Summit Station, including radar attenuation length measurements from air-borne radio sounding in the Greenland Ice Sheet Project (MacGregor and others, Reference MacGregor2015) and in situ radio sounding to investigate layering in the ice (Paden and others, Reference Paden2005). We include a comparison of our reported attenuation with previous measurements at Summit Station and Antarctica in Section 4.8.

2. Experimental approach

Our approach is similar to previous work in the astro-particle physics field (Barwick and others, Reference Barwick, Besson, Gorham and Saltzberg2005; Besson and others, Reference Besson2008; Barrella and other, Reference Barrella, Barwick and Saltzberg2011; Avva and others, Reference Avva, Kovac, Miki, Saltzberg and Vieregg2015; Hanson and others, Reference Hanson2015).

We transmit an impulsive, broadband radio signal downwards into the ice via a wideband, directional antenna, and measure the return signal as a voltage versus time trace on a second, identical antenna. The transmitted signal propagates through the ice sheet, reflects off of the bedrock, and returns to the receiving antenna on the surface. After correcting for geometric path loss, bedrock coefficient of reflection and electric field amplification from the focusing effect of the firn, the remaining power loss is attributed to absorption and scattering in the ice. Note that, experimentally, we do not distinguish between the two – our quoted attenuation length implicitly includes both effects. (Possible dispersive effects at the bedrock are quantified in Section 4.2.) To reduce systematic uncertainties, we remove the system response of the electronics (initial impulse, antenna response, cables, amplifiers, filters) by normalizing against a second measurement run in air. We recreate the through-ice setup (same initial impulse, antenna polarization, cables, amplifiers and filters but additional attenuators) for two antennas transmitting over a short distance in air and thereby largely cancel dependence on the system response from the in-ice data run.

We take the reflection off the bedrock to be specular. Motivation for this decision is described in Section 4.2 below. Given a specular reflection, the radar range equation reduces to the Friis equation (Friis, Reference Friis1946), and the direct through-air transmission formalism is applicable to our observed bedrock echoes. In that case, following the notation used by Avva and others (Reference Avva, Kovac, Miki, Saltzberg and Vieregg2015), the ratio of recorded voltages for each configuration (through-air vs. through-ice) as a function of frequency ν is,

(1)$${V_{\nu, {\rm ice}}\over V_{\nu, {\rm air}}} = \sqrt{F_{\,f} R} T_{{\rm ratio}} {d_{{\rm air}}\over d_{{\rm ice}}} \exp\left(-{d_{{\rm ice}}\over \langle L_\alpha \rangle}\right),\; $$

where d air and d ice are the distances the ice and air signals travel between antennas, R is the power reflection coefficient of the bedrock, F f is a focusing factor from the changing index of refraction in the firn (Stockham and others, Reference Stockham, Macy and Besson2016; Stockham, Reference Stockham2018), T ratio corrects for the change of transmission coefficient at the antenna feed between the antenna operating in air and in ice, and $\langle L_\alpha \rangle$ is the depth-averaged electric field attenuation length over the entire depth of the ice. Solving for the attenuation gives,

(2)$$\langle L_\alpha \rangle = d_{{\rm ice}} / \ln \left(\sqrt{F_{\,f}R} T_{{\rm ratio}} {V_{\nu, {\rm air}}\over V_{\nu, {\rm ice}}}{d_{{\rm air}}\over d_{{\rm ice}}} \right).$$

This equation differs from that in Avva and others (Reference Avva, Kovac, Miki, Saltzberg and Vieregg2015) by the inclusion of the focusing factor, which arises from the amplification of field strength from propagation in the firn (Sodha and others, Reference Sodha, Sirohi and Chakravarti1969). The addition of the focusing factor modifies the Avva and others (Reference Avva, Kovac, Miki, Saltzberg and Vieregg2015) bulk attenuation result from $947^{ + 92}_{-85}$ to $913^{ + 85}_{-79}$ m. We assume that the firn has an index of refraction varying linearly with the ice density profile ρ(z) at Summit Station (Kovacs and others, Reference Kovacs, Gow and Morey1995). The density profile has been experimentally measured (Alley and Koci, Reference Alley and Koci1988; Alley and others, Reference Alley1997; Hawley and others, Reference Hawley, Morris and McConnell2008; Arthern and others, Reference Arthern, Corr, Gillet-Chaulet, Hawley and Morris2013) and fit to a double exponential (Deaconu and others, Reference Deaconu2018), leading to a refractive index varying between n ~ 1.4 at the surface and n ~ 1.78 in deep (< 100 m) ice. The changing index of refraction focuses power from the transmitter on the downwards path, such that the Fresnel zone radius at the bedrock is reduced relative to the constant refractive index case. After reflection at the bedrock, the signal is partially de-focused on the return path, however the electric field areal flux density at the surface receiver is still amplified compared to the n = constant case. The focusing factor is equivalent to a correction to the expected 1/R 2 geometric spread factor, mathematically formalized by Červený and others (Reference Červený, Langer and Pšenčík1974). The focusing factor is present in MacGregor and others (Reference MacGregor2015), Matsuoka and others (Reference Matsuoka, Uratsuka, Fujita and Nishio2004), and Stockham and others (Reference Stockham, Macy and Besson2016), among others. We have used the finite-different time domain electrodynamics simulation software MEEP (Oskooi and others, Reference Oskooi2010) to confirm this effect.

3. Experimental setup

A system diagram of the experiment is presented in Figure 1. This measurement was performed in August 2021 at Summit Station, Greenland, using a separation distance of 244 m between the transmitter (coordinates 72.5801° N, 38.4569° W) and receiver (coordinates 72.5786° N, 38.4527° W) sites. The large separation distance assured that direct propagation from the transmitting antenna to the receiving antenna did not saturate the receiving amplifier. All antennas used were commercially available Create CLP-5130-2N (Creative Design Corp, 1994) log-periodic dipole antennas (LPDA) with ~ 8 dBi in-air forward gain over the band 105–1300  MHz.

Fig. 1. Diagram of experimental setup for bedrock reflection. On the transmitting side, we use a self-triggering FID Technologies + 5 kV high voltage pulse generator connected to a buried log-periodic dipole antenna (LPDA); an AVTECH fast pulse generator triggered by the FID pulser is connected to an in-air LPDA. On the receiving side, both the buried downward-pointing, and the in-air LPDAs are connected to a + 59 dB low noise amplifier; those outputs are then recorded on a Tektronix digital oscilloscope, triggered by the in-air signal.

Due to the large distance between stations, two parallel electronics signal chains are used, for triggering and the bedrock echo measurement, respectively. The bedrock echo electronics signal chain starts with the self-triggered high voltage FID TechnologyFootnote 1 model FPG6-1PNK pulse generator, which delivers a + 5 kV signal to a 50Ω coaxial feed. After the FID output, we apply a 100 MHz high pass filter using a MinicircuitsFootnote 2 NHP-100 filter. Following the filter, the signal is conveyed over 12 m of LMR-400 50Ω coaxial cable to an LPDA buried in the ice and pointed vertically downwards toward the bedrock; the bedrock-reflected return signal is then measured by a similarly-buried, downwards-pointing, receiver LPDA. The receiving antenna, located 244 m away along the surface of the ice, is aligned with the antenna tines parallel and collinear to the ones of the transmitting antenna, so that each antenna is in the gain null of the other to minimize contamination from horizontal ray paths. After measurement in the receiving antenna, the signal travels over a 10 m LMR-400 cable, bandpass filtered from 145–575 MHz using Minicircuits VHF-145+ and VLF-575+ filters, and then amplified by a custom RNO-G design low-noise amplifier with +59 dB of gain over the band 80–750 MHz. After the amplifier, the signal is bandpass filtered again using Minicircuits NHP-200 and VLF-575 + filters, and then recorded on a 2 GHz-bandwidth Tektronix MSO5204B oscilloscope.

The oscilloscope is triggered by the second electronic signal chain, ensuring a stable trigger over the distance between transmitter and receiver. The second chain begins with an AVTECH AVIR-1-C pulse generator triggered by the FID pulse generator TRG OUT, producing an impulsive, O(1  ns) pulse. The pulse generator is connected over a 12 m LMR-400 cable to an elevated, in-air LPDA pointed at a similarly elevated receiver LPDA located 244 m away and viewing the transmitter on boresight. The received in-air signal is then attenuated by 20 dB to prevent saturation, bandpass filtered using Minicircuits VHF-145 + and VLF-575 + filters, amplified by + 59 dB using the RNO-G low-noise amplifier, bandpass filtered again using Minicircuits VHF-145+ and VLF-575+ filters and finally captured by the oscilloscope. This in-air signal was used to trigger the oscilloscope, and therefore provides the reference t 0 for our measurements.

The oscilloscope was set to collect data over a 50 μs window; 10 000 individual triggers are averaged to suppress incoherent noise contributions and that average is written to scope memory. Twenty 10 000-event runs were collected, and then again averaged in post-processing, bringing the total number of triggers to 200 000.

To perform the air $\to$ air normalization run, we swapped the cables for the in-ice antennas with those from the in-air antennas. On the receiving side, two modifications were made: we added a 46 dB attenuator to prevent amplifier saturation and, for this configuration, we self-triggered on the arriving signal.

4. Experimental results

4.1. Bedrock power reflection coefficient

The power reflection coefficient at the ice-rock interface is not well-known and constitutes the largest uncertainty in our measurement of attenuation length. Taking an approach similar to Avva and others (Reference Avva, Kovac, Miki, Saltzberg and Vieregg2015), we take the power reflection coefficient to have a mean value of 0.215, a typical value for ice-bedrock interfaces as derived from radio sounding experiments (Barwick and others, Reference Barwick, Besson, Gorham and Saltzberg2005; Christianson and others, Reference Christianson2016). For uncertainty analysis, we assume the reflection coefficient can be drawn from a probability density function, uniformly distributed in the log of the reflection coefficient over the range from 0.01 to 1.0, which represent plausible extrema for the interface, from a frozen bedrock with high water content to an underlying layer of water (Christianson and others, Reference Christianson2016).

The observed return echo, in principle, could include both coherent and also incoherent contributions. Whereas the former sum linearly with the number of average triggers, the latter will scale as the square root of the number of events averaged N avg (Paden, Reference Paden2006). We have explicitly verified that our final results are insensitive to N avg, consistent with the assumption that the observed specular return echo, after subtracting the contribution from noise, is dominated by coherent scattering.

4.2. Bedrock echo

The reflection from the bedrock is visible above thermal noise in the time domain voltage trace of the receiving antenna, at a signal onset time of 35.55 μs after the oscilloscope trigger (Fig. 2). The bedrock echo is observed to include two components: a predominantly-specular, sharp, faster impulse (of duration ~ 500 ns), and a long (> 2 μs) extended signal which we associate with more diffuse, multi-path reflections off irregular features, both on the surface of, and within, the underlying bed reflector. For the purposes of the bulk radio attenuation measurement, and since an extended tail is not present for the in-air normalization run, we restrict consideration to the fast, specular component.

Fig. 2. Top: Recorded voltage as a function of time for the receiving in-ice antenna. Bottom: Recorded power, integrated in a sliding window of 100 ns to account for the group delay of the LPDA antennas. The specular component of the bedrock echo ‘signal’ is highlighted in magenta. Sub-surface internal layer reflections are visible at times earlier than 22 μs, after which noise dominates up to the point at which the bedrock echo is evident.

The uncertainty in the time window of the specular reflection is of O(10 ns), dominated by noise fluctuations at the edges of the window. This uncertainty is neglected since it is sub-dominant relative to the other systematic uncertainties in our final measurement. The final window start and end times are therefore defined to be 35.55 and 36.05 μs, respectively.

To determine the impact of neglecting the diffuse component of the bedrock echo on our final measured value of bulk attenuation, we investigate the dependence of our numerical result on the window length used in our analysis. We expect the measured attenuation to increase with increased window length due to the extended integration of power returning from the bedrock; at frequencies below 250 MHz, we obtain a $\sim 10\percnt$ larger attenuation length, but with increased uncertainty, as seen in Figure 3. At frequencies above 250 MHz, there is a negligible increase in attenuation length with increased window length. We note that the additional attenuation length is within the systematic uncertainties of our stated result due to the large bedrock reflection coefficient uncertainty.

Fig. 3. Measurement of the depth-averaged electric field attenuation at 200 MHz as a function of the window length used to select the bedrock echo. We note that including the entire diffuse component of the reflection into the final attenuation calculation increases the final result by no more than 10%.

The relationship between bulk ice attenuation and received power is different for the specular vs. diffuse components. We define attenuation from the Friis transmission equation (implicit in Eqn (1)) (Friis, Reference Friis1946). The Friis transmission equation is applicable for a specular reflection as it assumes direct line-of-sight propagation without interference within the first Fresnel zone, leading to a geometric path loss $\propto d_{{\rm ice}}^2$. The radar range equation is more applicable to the diffuse component since it includes power contributions from a rough surface via the definition of a radar cross-section and a geometric path loss ∝ (d ice/2)4 (Balanis, Reference Balanis2016). We find that use of the radar range equation instead of the Friis equation over the combined specular and diffuse components also increases the measured attenuation length at lower frequencies by a maximum of 10%, albeit introducing more model dependence from the unknown value of bedrock radar cross-section.

Based on the relatively small increase in obtained attenuation length from including the diffuse component, we quote our final result based on the Friis equation; this choice is consistent with previous similar measurements (Paden and others, Reference Paden2005; Barwick and others, Reference Barwick, Besson, Gorham and Saltzberg2005; Besson and others, Reference Besson2008; Barrella and other, Reference Barrella, Barwick and Saltzberg2011; Hanson and others, Reference Hanson2015; Avva and others, Reference Avva, Kovac, Miki, Saltzberg and Vieregg2015).

4.3. Bedrock depth

The bedrock depth can be derived from the absolute time of flight of the transmitted pulse and a model for the index of refraction of the ice as a function of depth. To reduce systematic biases (from location extrapolation and bedrock radio properties) and also as a cross check of our absolute timing calibration, we measure the bedrock depth from these data, rather than relying on previous measurements. The relationship between time of flight (Δt) and bedrock depth (half of the total distance propagated by the transmitted signal, d ice) can be found by solving for d ice in the integral:

(3)$$\Delta t = {2\over c} \int_0^{d_{{\rm ice}} / 2} n( z) dz,\; $$

where n(z) is the model for the index of refraction as a function of depth [m]. Index of refraction is related to dielectric constant (ε′) via $n( z) = \sqrt {\epsilon '( z) }$. The value of ε′(z) is derived from its relationship with measured ice density (ρ,[kg/m3]) (Kovacs and others, Reference Kovacs, Gow and Morey1995; Barwick and others, Reference Barwick2018):

(4)$$\epsilon'( z) = \left(1 + 0.854 \rho( z) \right)^2.$$

The parameterization of the dependence of ice density on depth follows Deaconu and others (Reference Deaconu2018), who performed a double exponential fit with a critical density at a depth of 14.9 m (Herron and Langway, Reference Herron and Langway1980):

(5)$$\rho( z) = \left\{\matrix {0.917 - 0.594 e^{-z / 30.8} \hfill & z \leq 14.9\,{\rm m} \cr 0.917 - 0.367 e^{-(z - 14.9) / 40.5} & z > 14.9\,{\rm m}}}$$

The uncertainty on the depth determination arises primarily from the uncertainty in the asymptotic index of refraction of deep glacial ice, which we take to be n = 1.78 ± 0.03 (Bogorodsky and others, Reference Bogorodsky, Bentley and Gudmandsen1985) for ice below the firn (deeper than 100 m at Summit Station). Using this refractive index profile, we calculate the bedrock to be at a depth of $3004^{ + 50}_{-52}$ m, corresponding to $d_{{\rm ice}} = 6008^{ + 100}_{-104}$ m for the through-ice bedrock echo total travel distance.

We note that, while our transmitting and receiving antennas were separated by 244 m, the through-ice signal approximately propagates vertically and Eqn (3) holds true. Over the measured 6 km propagation distance, horizontal propagation results in 4 m extra path length.

This bedrock depth is consistent with previous measurements of the bedrock depth from Greenland Ice Core Project (1994), at 3053.5  in 1993, and from Avva and others (Reference Avva, Kovac, Miki, Saltzberg and Vieregg2015), at $3014^{ + 48}_{-50}$ m in 2015.

4.4. Antenna coupling

The antenna transmission coefficient is defined as the quantity of power transmitted by the antenna from an incident radio frequency signal on a 50 Ω transmission line at the antenna feed point (S 21 in the scattering matrix). The transmission coefficient depends upon the dielectric properties of the antenna's embedded environment (Barwick and others, Reference Barwick2015; Glaser and others, Reference Glaser2019). To increase the power transmitted into the ice, we buried our antennas so that all active conductors were at least ~ 20 cm below the surface, thereby embedding them in an environment of n ≈ 1.4. The antennas used for the normalization are in air, for which n ~ 1.0. To correct for the change in match, we calculate a T ratio from the measured reflection coefficient (S 11 as shown in Figure 4) of the four antennas, two in air and two in ice, taken in the field. Assuming that all power not reflected at the feed is transmitted, the ratio becomes, in terms of the reflection coefficient S 11 in dB,

(6)$$T_{{\rm ratio}} = {1 - 10^{S_{11, {\rm ice}} / 10}\over 1 - 10^{S_{11, {\rm air}} / 10}}.$$

Fig. 4. Measured S 11 of each antenna used in the experiment. The difference in the low-frequency cutoff of the antenna when it is embedded in the ice compared to in air is due to the different indices of refraction of the two environments.

The antennas were found to transmit nearly all incident power in the frequency range of interest (150–550 MHz) in both the in-air and in-ice cases, resulting in a small T ratio correction. Averaged over frequency in the range of interest, T ratio = 1.00 ± 0.05, with the uncertainty assessed empirically from the variance of the measured match over the band. Our result is consistent with T ratio measured by other groups using the same or similar antennas (Barrella and other, Reference Barrella, Barwick and Saltzberg2011; Avva and others, Reference Avva, Kovac, Miki, Saltzberg and Vieregg2015).

In situ measurements, as well as simulations (Barwick and others, Reference Barwick2015; Glaser and others, Reference Glaser2019), have shown that the frequency-dependent antenna gain G 0, measured in air, also changes when the antenna is embedded in a dielectric medium. This change can be modeled as a down-shift in frequency, by the index of refraction of the medium (G′(ν) = G 0(ν/n)). For the LPDA used in this work, the gain over the frequency band of interest is uniform to < 0.5 dBi (Creative Design Corp, 1994), rendering the shift between in-air and in-ice measurements a subdominant systematic bias. The down-shift in frequency will cause a corresponding shift in low-frequency cutoff both in the gain and in the S 11 (as shown in Fig. 4), but the cutoff in both environments is below the high-pass filter of our analysis.

We note that there will be different contributions from the ice surface in both the in-air and in-ice antenna responses. We neglect these effects because they are likely to be small so long as the directional antennas (front-to-back [F/B] ratio of the LPDA ~–15 dB) are pointed away from the surface (Barwick and others, Reference Barwick2015).

4.5. Firn focusing

There is a geometric amplification of the bedrock echo electric field at the surface of the ice due to the changing index of refraction of the firn. To calculate the power focusing factor, for a negligibly thin firn layer, straightforward application of Snell's law prescribes that the electric field flux density measured at the surface, after bedrock reflection, follows (Stockham and others, Reference Stockham, Macy and Besson2016):

(7)$$F_{\rm f} = \left({n( z = 0) \over n( z = d_{{\rm ice}} / 2) }\right)^2.$$

We have verified that this equation agrees with a ray tracing simulation and 3D FDTD simulation. The uncertainty on the focusing factor arises from the uncertainty in the index of refraction model. Using n = 1.78 ± 0.03 (Bogorodsky and others, Reference Bogorodsky, Bentley and Gudmandsen1985) (as described previously in Section 4.3), and n = 1.4 ± 0.1 for the surface ice (Bogorodsky and others, Reference Bogorodsky, Bentley and Gudmandsen1985), we obtain a final focusing factor of F f = 1.61 ± 0.24.

4.6. In-air normalization amplitude

The amplitude of the signal from the in-air normalization run can be systematically biased from reflections off of the ice surface, increasing or decreasing the recorded power observed from the direct line-of-sight signal. Given the antenna heights above the ice (1.5 m) and distance between antennas (244 m), the first Fresnel zone is comprised by a nearly uniform, planar surface ice reflector, at all frequencies of interest. This leads to potential interference from reflections, depending on geometry: direct rays will interfere destructively/constructively with rays at the center/periphery of the Fresnel zone. To quantify any possible systematic bias, we compare the data against the absolute amplitude expectation of the signal from simulation. The absolute amplitude is derived from a measurement of the FID pulse shape, amplifier response, filter response, free space path loss and two independent simulations of the LPDA antenna response. Our simulations use either the Method of Moments software WIPL-DFootnote 3 or Finite Difference Time Domain software xFDTDFootnote 4, and have been found to agree with anechoic chamber measurements to 10% uncertainty (Barwick and others, Reference Barwick2017, Reference Barwick2015). The comparison of the simulated result with our data, seen in Figure 5, demonstrates that any possible systematic bias is not greater than 10% in voltage, consistent with previous results (Barwick and others, Reference Barwick2017, Reference Barwick2015).

Fig. 5. Comparison of data from the in-air normalization run compared against an absolute amplitude expectation as derived from two separate LPDA simulations. No systematic bias is evident within the ±10% voltage uncertainty in the antenna model (Barwick and others, Reference Barwick2015, Reference Barwick2017), over the frequency band of this analysis. The sharp dip between 180 and 220 MHz seen in data and simulations is most likely due to fine details in tine length and separation, which may be difficult to accurately simulate (Barwick and others, Reference Barwick2017).

4.7. Result and error analysis

The measured depth-averaged field attenuation length is presented in Figure 6. The previous in situ measurement reported by Avva and others (Reference Avva, Kovac, Miki, Saltzberg and Vieregg2015) is included for comparison, adjusted to remove the systematic bias from firn focusing not previously included in that analysis, resulting in a correction of their published result (from $947^{ + 92}_{-85}$ to $913^{ + 85}_{-79}$ m). We report the measurement only within the bandpass limits of our system, over which we have the highest sensitivity and lowest systematic biases from antenna modeling and filter response. Beyond the system bandpass limits, we have checked that our procedure yields an attenuation length numerically consistent with zero as expected for a noise-dominated regime.

Fig. 6. Measurement of the depth-averaged electric field attenuation as a function of frequency at Summit Station, within the system bandpass (shown as red dashed lines). Definitions of the error bars (and displayed upper limits) are provided in the text. The corrected result from Avva and others (Reference Avva, Kovac, Miki, Saltzberg and Vieregg2015) is shown for comparison. The frequency ν in the fit has units of MHz.

Contributions of different sources of uncertainty are calculated using a Monte Carlo method. We numerically calculate the estimated probability density distribution (PDF) of the bulk electric field attenuation within each frequency bin by repeatedly drawing random values of each component of the final measurement from their respective PDFs. For systematic uncertainties, we assume that each measured quantity used in the calculation of bulk field attenuation is uncorrelated and has a PDF either of a normal distribution (as is the case for F f, T ratio and d ice) or the distribution already described in the text (as is the case for R). The main component of statistical uncertainty is due to fluctuations in the power contributions from thermal noise in the recorded oscilloscope trace. In the 150–300 MHz band, the uncertainty from noise statistical fluctuation is sub-dominant to systematic uncertainties, contributing less than $10\percnt$ to the quoted uncertainty of each frequency bin. The final measurement is reported as a central value with one standard deviation (statistical plus systematic) error bars for those frequency bins that yield statistically significant results. For all other frequency bins, we report a 95% confidence level upper limit.

It is important to note that the majority of uncertainties are correlated in each frequency bin of the final measurement, with the primary contributions to the uncertainty arising from finite noise statistics and small systematic biases from the difference in LPDA response in ice vs. air. For the linear fit presented below, the visibly high goodness of fit is due to this correlation of uncertainties between frequency bins.

The reported bulk attenuation length includes losses from layer scattering. While expected to be a subdominant effect for vertical propagation due to the low reflection coefficient of the observed layers (Paden and others, Reference Paden2005), quantifying effects due to layer scattering for the more horizontal neutrino geometries must account for the larger Fresnel reflection coefficients at more glancing layer incidence angles (Aguilar and others, Reference Aguilar2022).

Birefringence of the ice can result in rotation of the signal that is dependent on polarization, leading to apparent loss of power at the co-polarized receiver. Measurements of the crystal orientation at Summit Station indicate uniaxial fabric at all depths (Thorsteinsson, Reference Thorsteinsson1996), unlike South Pole (Matsuoka and others, Reference Matsuoka2003; Barwick and others, Reference Barwick, Besson, Gorham and Saltzberg2005), indicating that birefringence will matter less for the Greenland site, though this remains to be quantified.

4.8. Discussion and summary

We derive the electric field attenuation length as a function of depth using a model of the field attenuation's dependence on temperature and chemical impurities, as described in detail in MacGregor and others (Reference MacGregor2007), MacGregor and others (Reference MacGregor2015) and Wolff and others (Reference Wolff, Miners, Moore and Paren1997). For a medium with non-zero intrinsic conductivity, the attenuation length can be expressed as:

(8)$$L_{\alpha}( f,\; z) = A( f) {\epsilon_0 \sqrt{\epsilon'_r( z) } c\over \sigma_{\infty}( z) },\; $$

where c is the speed of light, ε 0 is the permittivity of free space, εr(z) is the real component of the relative permittivity at a given depth z (defined in Eqn (4)), σ (z) is the infinite frequency limit of the electrical conductivity at a given depth, and the parameter A(f) is extracted experimentally, by requiring that the integrated, depth-dependent attenuation match our measured value of full-path attenuation, at each frequency f. Fujita and others (Reference Fujita, Matsuoka, Ishida, Matsuoka and Mae2000) suggest that the infinite frequency limit conductivity is valid at radio frequencies, since the molar conductivity does not change from low frequency (the frequency at which reference conductivity was measured, 0.3–3 MHz (Wolff and others, Reference Wolff, Miners, Moore and Paren1997)) to our frequency band (VHF/UHF, 150-550 MHz). The infinity frequency conductivity is related to chemical impurities and temperature via,

(9)$$\matrix{\eqalign{\sigma_\infty( z) & = \sigma_{{\rm pure}} \exp \left[{E_{{\rm pure}}\over k} \left({1\over T_{\rm r}} - {1\over T( z) }\right)\right]\cr & \quad + \mu_{{\rm H}^ + } [ {\rm H}^ + ] ( z) \exp \left[{E_{{\rm H}^ + }\over k} \left({1\over T_{\rm r}} - {1\over T( z) }\right)\right]\cr & \quad + \mu_{{\rm Cl}^-} [ {\rm Cl}^-] ( z) \exp \left[{E_{{\rm Cl}^-}\over k} \left({1\over T_{\rm r}} - {1\over T( z) }\right)\right]\cr & \quad + \mu_{{\rm NH}_4^ + } [ {\rm NH}_4^ + ] ( z) \exp \left[{E_{{\rm NH}_4^ + }\over k} \left({1\over T_{\rm r}} - {1\over T( z) }\right)\right],\; }}$$

where k is Boltzmann's constant, T(z) is in-ice temperature at a given depth as measured at the GRIP borehole Greenland Ice Core Project (1994)Footnote 5, T r is a reference temperature, σ pure is the conductivity of pure ice, $\mu _{\rm {H}^ + }$, $\mu _{\rm {Cl}^-}$ and $\mu _{\rm {NH}_4^ + }$ are molar conductivities, E pure, $E_{\rm {H}^ + }$, $E_{\rm {Cl}^-}$ and $E_{\rm {NH}_4^ + }$ are activation energies, and [H+](z)Footnote 6, [Cl](z)Footnote 7 and $[ \rm {NH}_4^ + ] ( z)$Footnote 7 are depth-dependent molar concentrations as measured at the GRIP borehole (Legrand and others, Reference Legrand, De Angelis and Maupetit1993; Greenland Ice Core Project, 1994; Wolff and others, Reference Wolff1995). The values of molar conductivities, conductivity of free ice and activation energies as measured by MacGregor and others (Reference MacGregor2007) and MacGregor and others (Reference MacGregor2015) and used in this analysis are given in Table 1. Note that, in this formulation, the temperature dependence of the attenuation length is explicitly absorbed into the conductivity dependence on temperature – the zero conductivity limit would correspond to no absorption and, correspondingly, no attenuation, for any temperature.

Table 1. Values of parameters used in the conductivity model of ice at summit station.

Compiled from MacGregor and others (Reference MacGregor2007) and MacGregor and others (Reference MacGregor2015) and from the GRIP borehole (Greenland Ice Core Project, 1994). Molar concentrations and ice temperature are tables of data measured at the GRIP borehole, and are available at the corresponding links in the footnotes.

Using the model of electric field attenuation length as a function of depth, we then unfold and solve for the parameter A(f), requiring that the depth-integrated attenuation matches our measured depth-averaged attenuation. The result reported at 300 MHz is plotted in Figure 7. We note a conspicuous enhancement in the $L_\alpha$(z) profile at depth $\gtrsim 1600$ m; this feature tracks a similarly precipitous drop in the tabulated GRIP H + and $NH_4^ +$ molar concentrations at z ~ 1600 m, resulting in a corresponding enhanced radio-frequency transparency at those depths.

Fig. 7. Electric field attenuation length as a function of depth at 300 MHz as derived from the two models of ice attenuation, one derived from temperature and chemical impurities (MacGregor and others, Reference MacGregor2007, Reference MacGregor2015) (black) and another from temperature alone (Bogorodsky and others, Reference Bogorodsky, Bentley and Gudmandsen1985)(red). Hatched and filled regions denote ± 1σ.

Avva and others (Reference Avva, Kovac, Miki, Saltzberg and Vieregg2015) derived the electric field attenuation length as a function of depth using a simplified model of the field attenuation dependence on temperature alone. They assume a linear relationship between the log of the attenuation versus the temperature of the ice:

(10)$$L_{\alpha}( f,\; T( z) ) = A( f) \cdot 10^{m T( z) }.$$

The parameter m is taken to be the average of the two sites measured by Bogorodsky and others (Reference Bogorodsky, Bentley and Gudmandsen1985) and set equal to − 0.017 (° C)−1. The parameter A(f) is derived from measured value of bulk field attenuation at each frequency f and set equal to 200–280 m. For comparison, our result using the simplified model as reported at 300 MHz is plotted in Figure 7.

The average electric field attenuation length of the top 1500 m of ice is of particular interest to RNO-G as the majority of neutrino interactions detectable by the experiment occur in this region (Aguilar and others, Reference Aguilar2022). It can be extracted from the bulk result using the field attenuation versus depth relation defined by MacGregor and others (Reference MacGregor2007) and MacGregor and others (Reference MacGregor2015). The average field attenuation length for the top 1500 m of ice result is shown in Figure 8. For reference, the bulk attenuation measurement at 300 MHz is $756^{ + 71}_{-87}$ m ($-11.49^{ + 1.49}_{-0.99}$ dB/km) while the average attenuation measurement of the top 1500 m of ice is $926^{ + 107}_{-124}$ m ($-9.38^{ + 1.45}_{-0.97}$ dB/km).

Fig. 8. Measurement of the average electric field attenuation for the top 1500 m of the ice sheet, as a function of frequency at Summit Station, derived from the measured bulk field attenuation in Figure 6 and the relationship between attenuation and temperature. Overlaid is the ± 1σ confidence interval of a linear fit of the data. Parameters of the fit are described in the text. Frequency ν in the fit is in units of MHz.

Our measurement can be used in simulations that calculate RNO-G's sensitivity to astrophysical neutrinos. For those simulations, we include, for convenience, a linear fit to the average electric field attenuation length for the top 1500 m of ice, shown in Figure 8. We find a significant correlation between the two parameters of the linear fit, which yields a slope of − 0.81 ± 0.14 m/MHz (− 0.12 ± 0.02 dB/km/MHz), intercept of 1154 ± 121 m (− 7.53 ± 0.72 dB/km), and a correlation coefficient ρ = −0.95.

We compare our obtained result of the average electric field attenuation length for the top 1500 m at 300 MHz of $926^{ + 107}_{-124}$ m to other similar measurements. Avva and others (Reference Avva, Kovac, Miki, Saltzberg and Vieregg2015) extrapolated their results at 75–300 MHz and estimated $\langle L_\alpha \rangle = 1022^{ + 230}_{-253}$ m, consistent (within uncertainty) with the result presented herein. The electric field attenuation length at the South Pole has been measured, with the focusing factor included, to be $\langle L_\alpha \rangle = 1660^{ + 255}_{-120}$ m at 300 MHz for the top 1500 m of ice (Allison and others, Reference Allison2012), a factor 1.8 times longer than our measurement. A linear fit to the log of attenuation length versus temperature (only) from Bogorodsky and others (Reference Bogorodsky, Bentley and Gudmandsen1985), 20°C colder ice yields an expected attenuation length ~ 2.2 times longer at South Pole than at Summit, however, the higher concentration of Cl at South Pole (~ 4μM (MacGregor and others, Reference MacGregor2007) compared to ~ 0.5μM at Summit Station in the top 1500 m of ice) moderates this expectation.

Measurements of the bulk radio field attenuation length at Summit Station have also been performed using air-borne radio sounding data. MacGregor and others (Reference MacGregor2015) inferred the bulk radio field attenuation by comparing the relative strengths of internal reflectors, obtaining attenuation lengths of 750–850 m at 150–200 MHz around Summit Station, consistent (within uncertainty) with the result presented herein. Stockham and others (Reference Stockham, Macy and Besson2016) measured the bulk radio field attenuation using the relative strength of the radio echo from the snow surface and the bedrock, obtaining at an attenuation length of 546 ± 23 m at 150–200 MHz around Summit Station. This measurement is notably lower than our result and the measurement done by MacGregor and others (Reference MacGregor2015) potentially due to radiometric calibration issues (MacGregor and others, Reference MacGregor2015; Stockham and others, Reference Stockham, Macy and Besson2016).

The RNO-G experiment, currently under construction at Summit Station in Greenland, is set to be one of the world's largest particle detectors. Our measurement of the bulk electric field attenuation length at Summit Station is consistent with previous measurements, with reduced systematic uncertainties. Our measurement will ultimately increase the precision of RNO-G's UHEN sensitivity estimates, which will either better motivate upper limits in the case of a null result or decrease uncertainties on the measured flux of ultra-high energy neutrinos in the universe in the case of observation.

Acknowledgements

We would like to thank the staff at Summit Station and Polar Field Services for logistical support in every way possible and to our colleagues at the British Antarctic Survey for unending enthusiasm while building and operating the BigRAID borehole drill for RNO-G.

We would like to acknowledge our home institutions and funding agencies for supporting the RNO-G work; in particular the Belgian Funds for Scientific Research (FRS-FNRS and FWO) and the FWO programme for International Research Infrastructure (IRI), the National Science Foundation through the NSF Awards 2118315 and 2112352 and the IceCube EPSCoR Initiative (Award ID 2019597), the German research foundation (DFG, Grant NE 2031/2-1), the Helmholtz Association (Initiative and Networking Fund, W2/W3 Program), the University of Chicago Research Computing Center and the European Research Council under the European Unions Horizon 2020 research and innovation programme (grant agreement No 805486).

References

Aartsen, MG and 52 others (2018) Differential limit on the extremely-high-energy cosmic neutrino flux in the presence of astrophysical background from nine years of IceCube data. Physical Review D 98(6), 062003. doi: 10.1103/PhysRevD.98.062003CrossRefGoogle Scholar
Aguilar, JA and 54 others (2021) Design and sensitivity of the Radio Neutrino Observatory in Greenland (RNO-G). JINST 16(03), P03025. doi: 10.1088/1748-0221/16/03/P03025CrossRefGoogle Scholar
Aguilar, JA and 57 others (2022) Reconstructing the neutrino energy for in-ice radio detectors: a study for the Radio Neutrino Observatory Greenland (RNO-G). The European Physical Journal C 82(2), 147. doi: 10.1140/epjc/s10052-022-10034-4CrossRefGoogle Scholar
Alley, R and 11 others (1997) Visual-stratigraphic dating of the GISP2 ice core: basis, reproducibility, and application. Journal of Geophysical Research: Oceans 102(C12), 2636726381. doi: 10.1029/96JC03837CrossRefGoogle Scholar
Alley, R and Koci, B (1988) Ice-core analysis at Site A, Greenland: preliminary results. Annals of Glaciology 10, 14. doi: 10.3189/S0260305500004067CrossRefGoogle Scholar
Allison, P and 46 others (2012) Design and initial performance of the Askaryan radio array prototype EeV neutrino detector at the South Pole. Astroparticle Physics 35, 457477. doi: 10.1016/j.astropartphys.2011.11.010CrossRefGoogle Scholar
Allison, P and 72 others (2020) Constraints on the diffuse flux of ultrahigh energy neutrinos from four years of Askaryan Radio Array data in two stations. Physical Review D 102(4), 043021. doi: 10.1103/PhysRevD.102.043021CrossRefGoogle Scholar
Anker, A and 23 others (2019) Targeting ultra-high energy neutrinos with the ARIANNA experiment. Advances in Space Research 64, 25952609. doi: 10.1016/j.asr.2019.06.016CrossRefGoogle Scholar
Arthern, RJ, Corr, HFJ, Gillet-Chaulet, F, Hawley, RL and Morris, EM (2013) Inversion for the density-depth profile of polar firn using a stepped-frequency radar. Journal of Geophysical Research (Earth Surface) 118(3), 12571263. doi: 10.1002/jgrf.20089CrossRefGoogle Scholar
Askaryan, GA (1961) Excess negative charge of an electron-photon shower and the coherent radio emission from it. Soviet Physics JETP 41, 616618.Google Scholar
Avva, J, Kovac, JM, Miki, C, Saltzberg, D and Vieregg, AG (2015) An in situ measurement of the radio-frequency attenuation in ice at Summit Station, Greenland. Journal of Glaciology 61, 10051011. doi: 10.3189/2015JoG15J057CrossRefGoogle Scholar
Balanis, CA (2016) Antenna Theory: Analysis and Design, 4th ed. Hoboken, New Jersey: Wiley Interscience.Google Scholar
Barrella, T, Barwick, S and Saltzberg, D (2011) Ross Ice Shelf (Antarctica) in situ radio-frequency attenuation. Journal of Glaciology 57(201), 6166. doi: 10.3189/002214311795306691CrossRefGoogle Scholar
Barwick, SW, Besson, DZ, Gorham, PW and Saltzberg, DP (2005) South Polar in situ radio-frequency ice attenuation. Journal of Glaciology 51(173), 231238. doi: 10.3189/172756505781829467CrossRefGoogle Scholar
Barwick, SW and 15 others (2015) Time domain response of the ARIANNA detector. Astroparticle Physics 62, 139151. doi: 10.1016/j.astropartphys.2014.09.002CrossRefGoogle Scholar
Barwick, SW and 18 others (2017) Radio detection of air showers with the ARIANNA experiment on the Ross Ice Shelf. Astroparticle Physics 90, 5068. doi: 10.1016/j.astropartphys.2017.02.003CrossRefGoogle Scholar
Barwick, SW and 18 others (2018) Observation of classically ‘forbidden’ electromagnetic wave propagation and implications for neutrino detection. JCAP 07, 055. doi: 10.1088/1475-7516/2018/07/055CrossRefGoogle Scholar
Besson, DZ and 42 others (2008) In situ radioglaciological measurements near Taylor Dome, Antarctica and implications for ultra-high energy (UHE) neutrino astronomy. Astroparticle Physics 29(2), 130157. doi: 10.1016/j.astropartphys.2007.12.004.CrossRefGoogle Scholar
Bogorodsky, VV, Bentley, CR and Gudmandsen, PE (1985) Radioglaciology. Dordrecht, Holland: Reidel Publishing Co.CrossRefGoogle Scholar
Červený, V, Langer, J and Pšenčík, I (1974) Computation of geometric spreading of seismic body waves in laterally inhomogeneous media with curved interfaces. Geophysical Journal International 38(1), 919. doi: 10.1111/j.1365-246X.1974.tb04105.xCrossRefGoogle Scholar
Christianson, K and 6 others (2016) Basal conditions at the grounding zone of Whillans Ice Stream, West Antarctica, from ice-penetrating radar. Journal of Geophysical Research: Earth Surface 121(11), 19541983. doi: 10.1002/2015JF003806CrossRefGoogle Scholar
Creative Design Corp (1994) Instruction Manual CLP5130-2 Create Radio Antenna.Google Scholar
Deaconu, C and 8 others (2018) Measurements and modeling of near-surface radio propagation in glacial ice and implications for neutrino experiments. Physical Review D 98(4), 043010. doi: 10.1103/PhysRevD.98.043010CrossRefGoogle Scholar
Friis, H (1946) A note on a simple transmission formula. Proceedings of the IRE 34(5), 254256. doi: 10.1109/JRPROC.1946.234568CrossRefGoogle Scholar
Fujita, S, Matsuoka, T, Ishida, T, Matsuoka, K and Mae, S (2000) A summary of the complex dielectric permittivity of ice in the megahertz range and its applications for radar sounding of polar ice sheets. In Physics of ice core records, Hokkaido University Press, pp. 185–212.Google Scholar
Glaser, C and 8 others (2019) NuRadioReco: a reconstruction framework for radio neutrino detectors. The European Physical Journal C 79(6), 29. doi: 10.1140/epjc/s10052-019-6971-5.CrossRefGoogle Scholar
Gorham, PW and 42 others (2009) The Antarctic impulsive transient antenna ultra-high energy neutrino detector design, performance, and sensitivity for 2006–2007 balloon flight. Astroparticle Physics 32, 1041. doi: 10.1016/j.astropartphys.2009.05.003CrossRefGoogle Scholar
Gorham, PW and 59 others (2019) Constraints on the ultrahigh-energy cosmic neutrino flux from the fourth flight of ANITA. Physical Review D 99(12), 122001. doi: 10.1103/PhysRevD.99.122001CrossRefGoogle Scholar
Hanson, JC and 13 others (2015) Radar absorption, basal reflection, thickness and polarization measurements from the Ross Ice Shelf, Antarctica. Journal of Glaciology 61(227), 438446. doi: 10.3189/2015JoG14J214CrossRefGoogle Scholar
Hawley, RL, Morris, EM and McConnell, JR (2008) Rapid techniques for determining annual accumulation applied at Summit, Greenland. Journal of Glaciology 54(188), 839845. doi: 10.3189/002214308787779951CrossRefGoogle Scholar
Herron, MM and Langway, CC (1980) Firn densification: an empirical model. Journal of Glaciology 25(93), 373385. doi: 10.3189/S0022143000015239CrossRefGoogle Scholar
Kovacs, A, Gow, AJ and Morey, RM (1995) The in-situ dielectric constant of polar firn revisited. Cold Regions Science and Technology 23(3), 245256. doi: 10.1016/0165-232X(94)00016-QCrossRefGoogle Scholar
Kravchenko, I and 14 others (2003) Performance and simulation of the RICE detector. Astroparticle Physics 19, 1536. doi: 10.1016/S0927-6505(02)00194-9CrossRefGoogle Scholar
Legrand, M, De Angelis, M and Maupetit, F (1993) Field investigation of major and minor ions along summit (central greenland) ice cores by ion chromatography. Journal of Chromatography A 640(1), 251258. doi: 10.1016/0021-9673(93)80188-ECrossRefGoogle Scholar
MacGregor, JA and 5 others (2007) Modeling englacial radar attenuation at Siple Dome, West Antarctica, using ice chemistry and temperature data. Journal of Geophysical Research: Earth Surface 112(F3), 16. doi: 10.1029/2006JF000717.CrossRefGoogle Scholar
MacGregor, JA and 11 others (2015) Radar attenuation and temperature within the Greenland Ice Sheet. Journal of Geophysical Research: Earth Surface 120(6), 9831008. doi: 10.1002/2014JF003418CrossRefGoogle Scholar
Matsuoka, K and 6 others (2003) Crystal orientation fabrics within the Antarctic ice sheet revealed by a multipolarization plane and dual-frequency radar survey. Journal of Geophysical Research 108, 2499. doi: 10.1029/2003JB002425CrossRefGoogle Scholar
Matsuoka, K, Uratsuka, S, Fujita, S and Nishio, F (2004) Ice-flow-induced scattering zone within the antarctic ice sheet revealed by high-frequency airborne radar. Journal of Glaciology 50(170), 382388. doi: 10.3189/172756504781829891CrossRefGoogle Scholar
Oskooi, AF and others (2010) Meep: a flexible free-software package for electromagnetic simulations by the FDTD method. Computer Physics Communications 181(3), 687702. doi: 10.1016/j.cpc.2009.11.008CrossRefGoogle Scholar
Paden, J (2006) Synthetic aperture radar for imaging the basal conditions of the polar ice sheets. Ph.D. thesis, Department of Electrical Engineering and Computer Science, University of Kansas, Lawrence, Kansas.Google Scholar
Paden, J and 5 others (2005) Wideband measurements of ice sheet attenuation and basal scattering. IEEE Geoscience and Remote Sensing Letters 2(2), 164168. doi: 10.1109/LGRS.2004.842474CrossRefGoogle Scholar
Sodha, MS, Sirohi, RS and Chakravarti, AK (1969) Focusing factor in a stratified medium. British Journal of Applied Physics 2(2), 14571461. doi: 10.1088/0022-3727/2/4/317.Google Scholar
Stockham, M (2018) Bulk and surface radio-frequency response of ice. Ph.D. thesis, University of Kansas, Lawrence, Kansas.Google Scholar
Stockham, M, Macy, J and Besson, DZ (2016) Radio frequency ice dielectric permittivity measurements using CReSIS data. Radio Science 51, 194212. doi: 10.1002/2015RS005849CrossRefGoogle Scholar
Thorsteinsson, T (1996) Textures and fabrics in the GRIP ice core, in relation to climate history and ice deformation. Berichte zur Polarforschung (Reports on Polar Research), 205, 5271. doi: 10.1029/97JC00161.Google Scholar
Wolff, EW and 5 others (1995) Long-term changes in the acid and salt concentrations of the Greenland ice core project ice core from electrical stratigraphy. Journal of Geophysical Research: Atmospheres 100(D8), 1624916263. doi: 10.1029/95JD01174CrossRefGoogle Scholar
Wolff, EW, Miners, WD, Moore, JC and Paren, JG (1997) Factors controlling the electrical conductivity of ice from the polar regionsa summary. The Journal of Physical Chemistry B 101(32), 60906094. doi: 10.1021/jp9631543CrossRefGoogle Scholar
Figure 0

Fig. 1. Diagram of experimental setup for bedrock reflection. On the transmitting side, we use a self-triggering FID Technologies + 5 kV high voltage pulse generator connected to a buried log-periodic dipole antenna (LPDA); an AVTECH fast pulse generator triggered by the FID pulser is connected to an in-air LPDA. On the receiving side, both the buried downward-pointing, and the in-air LPDAs are connected to a + 59 dB low noise amplifier; those outputs are then recorded on a Tektronix digital oscilloscope, triggered by the in-air signal.

Figure 1

Fig. 2. Top: Recorded voltage as a function of time for the receiving in-ice antenna. Bottom: Recorded power, integrated in a sliding window of 100 ns to account for the group delay of the LPDA antennas. The specular component of the bedrock echo ‘signal’ is highlighted in magenta. Sub-surface internal layer reflections are visible at times earlier than 22 μs, after which noise dominates up to the point at which the bedrock echo is evident.

Figure 2

Fig. 3. Measurement of the depth-averaged electric field attenuation at 200 MHz as a function of the window length used to select the bedrock echo. We note that including the entire diffuse component of the reflection into the final attenuation calculation increases the final result by no more than 10%.

Figure 3

Fig. 4. Measured S11 of each antenna used in the experiment. The difference in the low-frequency cutoff of the antenna when it is embedded in the ice compared to in air is due to the different indices of refraction of the two environments.

Figure 4

Fig. 5. Comparison of data from the in-air normalization run compared against an absolute amplitude expectation as derived from two separate LPDA simulations. No systematic bias is evident within the ±10% voltage uncertainty in the antenna model (Barwick and others, 2015, 2017), over the frequency band of this analysis. The sharp dip between 180 and 220 MHz seen in data and simulations is most likely due to fine details in tine length and separation, which may be difficult to accurately simulate (Barwick and others, 2017).

Figure 5

Fig. 6. Measurement of the depth-averaged electric field attenuation as a function of frequency at Summit Station, within the system bandpass (shown as red dashed lines). Definitions of the error bars (and displayed upper limits) are provided in the text. The corrected result from Avva and others (2015) is shown for comparison. The frequency ν in the fit has units of MHz.

Figure 6

Table 1. Values of parameters used in the conductivity model of ice at summit station.

Figure 7

Fig. 7. Electric field attenuation length as a function of depth at 300 MHz as derived from the two models of ice attenuation, one derived from temperature and chemical impurities (MacGregor and others, 2007, 2015) (black) and another from temperature alone (Bogorodsky and others, 1985)(red). Hatched and filled regions denote ± 1σ.

Figure 8

Fig. 8. Measurement of the average electric field attenuation for the top 1500 m of the ice sheet, as a function of frequency at Summit Station, derived from the measured bulk field attenuation in Figure 6 and the relationship between attenuation and temperature. Overlaid is the ± 1σ confidence interval of a linear fit of the data. Parameters of the fit are described in the text. Frequency ν in the fit is in units of MHz.