skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.
Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, June 11 until 2:00 AM ET on Friday, June 12 due to maintenance. We apologize for the inconvenience.


Title: Data release for A search for extremely-high-energy neutrinos and first constraints on the ultra-high-energy cosmic-ray proton fraction with IceCube
<h1 id="summary">Summary</h1> <p>Title: Data Release for A search for extremely-high-energy neutrinos and first constraints on the ultra-high-energy cosmic-ray proton fraction with IceCube</p> <p>The IceCube observatory analyzed 12.6 years of data in search of extremely-high-energy (EHE) neutrinos above 5 PeV. The resultant limit of the search (Fig 1), and the effective area of the event selection (Fig 7), are provided in this data release.</p> <h1 id="contents">Contents</h1> <ul> <li><p>README file: this file</p> </li> <li><p><code>differential_limit_and_sensitivity.csv</code>: a comma separated value file, giving the observed experimental differential limit, and sensitivity, of the search as a function of neutrino energy. This is the content of Fig 1 in the paper. The first column is the neutrino energy in GeV. The second column is the limit in units of GeV/cm2/s/sr. The third column is the sensitivity in units of GeV/cm2/s/sr.</p> </li> <li><p><code>effective_area.csv</code>: a comma separated value file, giving the effective area of the search as a function of energy. This is the content of Fig 7 in the paper. The first column is the neutrino energy in GeV. The second column is the total effective area of the search, summed across neutrino flavors, and averaged across neutrinos and antineutrinos, in meters-squared. The third column is the effective area of the search for the average of electron neutrino and electron antineutrinos in units of meters-squared. The fourth column is the same as the third, but for muon-flavor neutrinos. The fifth column is the same as the third and fourth, but for tau-flavor neutrinos.</p> </li> <li><p><code>demo.py</code>: a short python script to demonstrate how to read the files. Run like <code>python demo.py</code>. A standard base python installation is sufficient, as the only dependencies are numpy and matplotlib.</p> </li> </ul> <h1 id="contacts">Contacts</h1> <p>For any questions about this data release, please write to analysis@icecube.wisc.edu</p>  more » « less
Award ID(s):
2411849
PAR ID:
10614228
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Harvard Dataverse
Date Published:
Subject(s) / Keyword(s):
Astronomy and Astrophysics Physics neutrino cosmic ray
Format(s):
Medium: X Size: 300; 4500; 1441; 1843 Other: text/tab-separated-values; text/tab-separated-values; text/x-python-script; text/markdown
Size(s):
300 4500 1441 1843
Right(s):
Creative Commons Zero v1.0 Universal
Sponsoring Org:
National Science Foundation
More Like this
  1. {"Abstract":["Results of the Search for eV-scale sterile neutrinos with the DeepCore "golden" sample in machine readable format. Allows reproduction of Figure 5 of the publication.\n\n<ul>\n<li> The file "wilks_contour_90pct.csv" contains the 90% contour line in CSV format. The columns "U\\mu4_sq" and "U\\tau4_sq" correspond to the _squared_ elements of the 3+1 mixing matrix, |U_{&mu;4}|^2 and |U_{&tau;4}|^2, respectively. The contour was drawn assuming that Wilks' theorem holds with two degrees of freedom. The points are ordered by nearest neighbors. </li>\n<li> The file "expected_sensitivity_90pct.csv" contains the expected sensitivity at 90% CL of the analysis for Asimov data produced at the null hypothesis. The format is identical to that of the result. </li>\n<li> The files "marginalized_umu4_sq.csv" and "marginalized_utau4_sq.csv" contain the marginalized scan of the test statistic over the matrix elements. The column "Metric Difference" refers to the delta chi-square metric used in the analysis. Since this metric is a delta chi-square (rather than a likelihood), it should _not_ be multiplied by two when calculating significances. They allow reproduction of the top and right panels of Figure 5 of the publication. </li>\n<li> The file "real_data_wilks_contours_reprod_mar_2025_best_fit_points.tab" contains the best-fit values of all free nuisance parameters in a |U_{&mu;4}|^2 and |U_{&tau;4}|^2 scan. </li>\n<li> The file "real_data_wilks_contours_reprod_mar_2025_metric_diff.tab" contains the \\chi^2 difference between the scan point and the global best-fit in a |U_{&mu;4}|^2 and |U_{&tau;4}|^2 scan. </li>\n</ul>"]} 
    more » « less
  2. {"Abstract":["<b>Data and Code Repository</b> <br>\nThis repository contains the anonymized dataset and analysis code associated with the paper titled "Promoting Sustainable Travel Modes Through Health and Active Lifestyle Messaging." <br>\n<b>Contents</b>\n<ul>\n <li>anonymized_data.csv: Contains the raw anonymized dataset (1 MB) used in the analysis.</li>\n <li>analysis_v1.ipynb: Python scripts for analysis, and visualization. (Last Modified July 7, 2025, 5:00PM)</li>\n <li>README.md: Description of the repository contents and usage.</li>\n <li>datadictionary.md: A detailed explanation of each variable in the final dataset. 68 unique variables, 4,840 observations.</li>\n</ul>\n<b>Requirements</b> <br>\nPackages required to run this analysis are pandas==2.0.3, numpy==1.24.1, statsmodel.api==0.14.1. This code was tested on Python 3.8.13 and 3.9.2 and on macOS Sequoia 15.5 and Google Colab CPUs. <br>\n<b>Structure of the code</b> <br>\nThe first code block loads the dataset and required packages. The second code block has helper function that generates dataframe for statistical analysis in the later blocks. The third code block has helper variables and functions to load model specifications and formatting model coefficients for analysis in the later blocks Code blocks four and above generate statistical results used in the paper. <br>\n<b>Output</b> <br>\nThis code package generates the necessary derived data consisting of odds ratios and uncertainty for Figure 1 and Figure 2 in the main document:\n<ul>\n <li>Main Document Fig. 1 Treatment effects of air quality impacts targeting bus transit and active lifestyle messaging targeting walking or biking.</li>\n <li>Main Document Fig. 2 Comparative treatment effects of health and active lifestyle messaging for all respondents and those with underlying health conditions.</li>\n</ul>\nThis code package generates the following tables:\n<ul>\n <li>Main Document Table 1 Heterogeneous treatment effects across key subgroups.</li>\n <li>Table S4. Treatment effects of air pollution exposure messaging in Experiment 1 (personal and community benefits targeting bus transit)</li>\n <li>Table S5. Treatment effects of air quality improvement messaging in Experiment 2 (personal and community benefits targeting bus transit)</li>\n <li>Table S6. Treatment effects of air quality improvement messaging in Experiment 2 (personal gain targeting walking)</li>\n <li>Table S7. Treatment effects of air quality improvement messaging in Experiment 2 (personal and community benefits targeting walking)</li>\n <li>Table S8. Treatment effects of air quality improvement messaging in Experiment 2 (personal gain targeting biking)</li>\n <li>Table S9. Treatment effects of air quality improvement messaging in Experiment 2 (personal and community benefits targeting biking)</li>\n <li>Table S10. Treatment effects of all active lifestyle messaging in Experiment 3 (personal and community benefits targeting walking and biking)</li>\n <li>Table S11. Treatment effects of all active lifestyle messaging in Experiment 3 (personal gain targeting biking)</li>\n <li>Table S12. Treatment effects of active lifestyle messaging in Experiment 3 (personal and community benefits targeting biking)</li>\n <li>Table S13. Treatment effects of step count messaging in Experiment 3 (personal and community benefits targeting walking)</li>\n <li>Table S14. Treatment effects of calories burned messaging in Experiment 3 (personal and community benefits targeting walking)</li>\n <li>Table S15. Treatment effects of heart health messaging in Experiment 3 (personal and community benefits targeting walking and biking)</li>\n <li>Table S16. Heterogeneous treatment effects of air pollution exposure messaging in Experiment 1 (personal gain targeting bus transit)</li>\n <li>Table S17. Heterogeneous treatment effects of air pollution exposure messaging in Experiment 1 (personal and community benefits targeting bus transit)</li>\n <li>Table S18. Heterogeneous treatment effects of air quality improvement messaging in Experiment 2 (personal gain targeting bus transit)</li>\n <li>Table S19. Heterogeneous treatment effects of air quality improvement messaging in Experiment 2 (personal and community benefits targeting bus transit)</li>\n <li>Table S20. Heterogeneous treatment effects of active lifestyle messaging in Experiment 3 (personal gain targeting walking)</li>\n <li>Table S21. Heterogeneous treatment effects of active lifestyle messaging in Experiment 3 (personal and community benefits targeting walking)</li>\n <li>Table S22. Heterogeneous treatment effects of different active lifestyle messaging in Experiment 3 targeting walking</li>\n <li>Table S23. Heterogeneous treatment effects of calories burned messaging for commuters with varying daily travel times in Experiment 3 (targeting walking)</li>\n</ul>\n<b>Replication</b> <br>\nSupporting replication code is also available here: https://github.com/asensio-lab/health-active-lifestyle. This code package was last replicated on July 7, 2025 by @YifanLiu0304 <br>\n<b>Declaration of generative AI and AI-assisted technologies in the coding process</b> <br>\nDuring the preparation of this work the authors used ChatGPT in order to debug code errors such as KeyError, syntax errors in python scripts that were used to generate statistical tables. After using this tool/service, the researchers reviewed, edited, and replicated all code."],"TechnicalInfo":["Python, 3.9.2"],"Other":["Supporting replication code is also available here: https://github.com/asensio-lab/health-active-lifestyle"]} 
    more » « less
  3. <p>NSF COLDEX performed two airborne campaigns from South Pole Station over the Southern Flank of Dome A and 2022-23 and 2023-24, searching for a potential site of a continuous ice core that could sample the mid-Pleistocene transition. Ice thickness data extracted from the MARFA radar system has allow for a new understanding of this region.</p> <p>Here we generate crustal scale maps of ice thickness, bed elevation, specularity content, subglacial RMS deviation and fractional basal ice thickness with 1 km sampling, and 10 km resolution. We include both masked and unmasked grids.</p> <p> The projection is in the SCAR standard ESPG:3031 polar stereographic projection with true scale at 71˚S.</p> <p>These geotiffs were generated using performed using GMT6.5 (<a href="https://doi.org/10.1029/2019GC008515">Wessel et al., 2019</a>) using the pygmt interface, by binning the raw data to 2.5 km cells, and using the <a href="https://github.com/sakov/nn-c"> nnbathy </a> program to apply natural neighbor interpolation to 1 km sampling. A 10 km Gaussian filter - representing typical lines spacings - was applied and then a mask was applied for all locations where the nearest data point was further than 8 km. </p> Ice thickness, bed elevation and RMS deviation @ 400 m length scale (<a href="http://dx.doi.org/10.1029/2000JE001429">roughness</a>) data includes the following datasets: <ul> <li> UTIG/CRESIS <a href="https://doi.org/10.18738/T8/J38CO5">NSF COLDEX Airborne MARFA data</a></li> <li> British Antarctic Survey <a href="https://doi.org/10.5285/0f6f5a45-d8af-4511-a264-b0b35ee34af6">AGAP-North</a></li> <li> LDEO <a href="https://doi.org/10.1594/IEDA/317765"> AGAP-South </a></li> <li> British Antarctic Survey <a href="https://doi.org/10.5270/esa-8ffoo3e">Polargap</a></li> <li> UTIG Support Office for Airborne Research <a href="https://doi.org/10.15784/601588">Pensacola-Pole Transect (PPT) </a></li> <li> NASA/CReSIS <a href="https://doi.org/10.5067/GDQ0CUCVTE2Q"> 2016 and 2018 Operation Ice Bridge </a> </li> <li> ICECAP/PRIC <a href="https://doi.org/10.15784/601437"> SPICECAP Titan Dome Survey </a> </ul> <p>Specularity content (<a href="https://doi.org/10.1109/LGRS.2014.2337878">Schroeder et al. 2014</a>) is compiled from <a href="https://doi.org/10.18738/T8/KHUT1U"> Young et al. 2025a </a> and <a href="https://doi.org/10.18738/T8/6T5JS6"> Young et al. 2025b</a>.</p> <p>Basal ice fractional thickness is complied from manual interpretation by Vega Gonzàlez, Yan and Singh. </p> <p>Code to generated these grids can be found at <a href="https://github.com/smudog/COLDEX_dichotomy_paper_2025"> at github.com </a></p> 
    more » « less
  4. {"Abstract":["Data and plotting scripts for reproducing plots from "Physics potential of the IceCube Upgrade for atmospheric neutrino oscillations". Please refer to the related publication for a detailed explanation of the sample and analysis.\n<br><br>\nContents: \n<ol>\n<li> <code>README.md</code> contains useful information about the contents of this data release.</li>\n<li> <code>example.py</code> Example script demonstrating how to load the csv files, plot the chi squared map, and extract the 90% sensitivity contours shown in Figures 11 and 15.</li>\n<li> <code>modchi2map_nufitwoSK.csv</code> chi2 map used to produce Figure 11 (left), in which the injected truth point is the best fit point from <a href="http://www.nu-fit.org/sites/default/files/v52.tbl-parameters.pdf">NuFit 5.2 w/o SK (upper panel)</a>, sin<sup>2</sup>(&theta;<sub>23</sub>)=0.572 and &Delta;m<sup>2</sup><sub>32</sub>=2.43&times;10<sup>-3</sup> eV<sup>2</sup>.</li>\n<li> <code>modchi2map_nufitwSK.csv</code> chi2 map used to produce Figure 11 (right), in which the injected truth point is the best fit point from <a href="http://www.nu-fit.org/sites/default/files/v52.tbl-parameters.pdf">NuFit 5.2 w/ SK (lower panel)</a>, sin<sup>2</sup>(&theta;<sub>23</sub>)=0.451 and &Delta;m<sup>2</sup><sub>32</sub>=2.43&times;10<sup>-3</sup> eV<sup>2</sup>.</li>\n<li> <code>modchi2map_icecube.csv</code> chi2 map used to produce Figure 15, in which the injected truth point is the best fit point from the latest published IceCube DeepCore oscillation result <a href="https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.134.091801">Phys. Rev. Lett. 134, 091801 (2025)</a>, sin<sup>2</sup>(&theta;<sub>23</sub>)=0.54 and &Delta;m<sup>2</sup><sub>32</sub>=2.40&times;10<sup>-3</sup> eV<sup>2</sup>.</li>\n</ol>\n\nPlease note: The CSV files are available for download through dataverse in either csv or tab format."]} 
    more » « less
  5. {"Abstract":["<b>NSF COLDEX CXA1 laser ranging, magnetic field strength and orientation, and positioning</b><br>\n<br>\nThis data collection represents georeferenced, time registered instrument measurements (Level 1B data) converted to SI units, and is of most interest to users who wish to reprocess the data. Users interested in geophysical observables (eg magnetic anomalies and surface elevation) should used the derived Level 2 datasets. The data format are space delimited ASCII files, following the formats used for UTIG predecessor ICECAP/OIB project at NASA's NSIDC DAAC. Fields are described in the # delimited detailed header for each granule.<br>\n<br>\n<b>DATASET ORGANIZATION</b><br>\nTransects were collected in preplanned systems with the following parameters:\n<ul>\n<li>CLX radials (CLX/MKB##/R###), attempting to emulate flow lines from Dome A and radiating (in the EPSG:3031 polar stereographic projection) from easting 965 km northing 385 km, with a separation of 0.25 degrees.</li>\n<li>CLX corridor (CLX/MKB##/X###) rotated from the EPSG:3031 polar stereographic projection at -150 degrees and separated by 10 km in the Y direction and 3.75 km in the X direction.</li>\n<li>CLX2 corridor (CLX2/MKB##/X###) rotated from the EPSG:3031 polar stereographic projection at -150 degrees and separated by 2.5 km in its Y direction and 2.5 km in its X direction.</li>\n<li>SAD corridor (SAD/MKB##/X###|Y####) designed to characterize the Saddle region near South Pole approximately perpendicular to the flow lines, rooted from the EPSG:3031 polar stereographic projection at -73.8 degrees and separated by 2.5 km in its Y direction and 2.5 km in the its X direction.</li>\n<li>Untargeted transit lines used the name of the expedition (CXA1) as the project, and used the flight and the increment within the flight to name the Transect (eg (CXA1/MKB2n/F10T02a).</li>\n</ul>\n\n<br>\n<b>POSITIONING</b><br>\nGPS provides the positioning (and timing) for all other data streams:<br>\n<ul>\n<li>IPUTG1B (ICECAP Positioning/UTIG GPS Level 1B): real time CA-code 1 Hz GPS position data</li>\n</ul>\n<br>\n<b>MAGNETICS</b><br>\nMagnetic data provides constraints on the depth to crystalline rock, the temperature structure of the crust, and the underlying tectonics<br>\n<ul>\n<li>IMGEO1B (ICECAP Magnetics: Geometrics Level 1B): georeferenced total magnetic field data; 10 Hz ASCII</li>\n<li>IMFGM1B (ICECAP Magnetics/Flux Gate Magnetometer Level 1B): georeferenced aircraft referenced magnetic field vectors; 10 Hz ASCII</li>\n</ul>\n<br>\n<b>ALTIMETRY</b><br>\nLaser altimetry provides for surface shape and altimetry validation. See <a NSF href="https://doi.org/10.18738/T8/99IEOG"> COLDEX 2022-23 Riegl Laser Altimeter Level 2 Geolocated Surface Elevation Triplets</a> for derived surface elevation data.<br>\n<ul>\n<li>ILUTP1B (ICECAP Lidar/UTIG Profiler Level 1B): georeferenced laser range data, no orientation corrections; 3.5 Hz ASCII</li>\n</ul>\n<br>\nThis work was supported by the U.S. National Science Foundation Center for Oldest Ice Exploration (NSF COLDEX), an NSF Science and Technology Center (NSF 2019719). We thank the NSF Office of Polar Programs, the NSF Office of Integrative Activities and the NSF Antarctic Infrastructure and Logistics Program, Kenn Borek Airlines and the Antarctic Support Contractor for logistical support. Additional support came from the Open Polar Radar project (NSF grant 2127606) award to the University of Texas at Austin and a gift to UTIG from the G. Unger Vetlesen foundation."]} 
    more » « less