Star forming regions (SFRs) have been postulated as possible sources of cosmic rays (CRs) in our galaxy. One example of a gamma-ray source associated with an SFR is the Fermi-LAT cocoon, an extended region of gamma-ray emission in the Cygnus X region and attributed to a possible superbubble with freshly accelerated CRs. Because the emission region is surrounded by ionization fronts, it has been named the "Cygnus cocoon". CRs in the cocoon could have originated in the OB2 association and been accelerated at the interaction sites of stellar winds of massive O type stars. So far, there is no clear association at TeV energies. Spectral and morphological studies of TeV gamma-ray emission detected by the High Altitude Water Cherenkov (HAWC) observatory at the 2HWC J2031+415 region reveal that the spectral energy distribution of the cocoon extends from GeV to at least tens of TeV. Using HAWC data, we are able to study the acceleration of particles to highest energies in the Cygnus OB2 SFR.

We describe plans for the development of the Southern Wide-field Gamma-ray Observatory (SWGO), a next-generation instrument with sensitivity to the very-high-energy (VHE) band to be constructed in the Southern Hemisphere. SWGO will provide wide-field coverage of a large portion of the southern sky, effectively complementing current and future instruments in the global multi-messenger effort to understand extreme astrophysical phenomena throughout the universe. A detailed description of science topics addressed by SWGO is available in the science case white paper [1]. The development of SWGO will draw on extensive experience within the community in designing, constructing, and successfully operating wide-field instruments using observations of extensive air showers. The detector will consist of a compact inner array of particle detection units surrounded by a sparser outer array. A key advantage of the design of SWGO is that it can be constructed using current, already proven technology. We estimate a construction cost of 54M USD and a cost of 7.5M USD for 5 years of operation, with an anticipated US contribution of 20M USD ensuring that the US will be a driving force for the SWGO effort. The recently formed SWGO collaboration will conduct site selection and detector optimization studies prior to construction, with full operations foreseen to begin in 2026. Throughout this document, references to science white papers submitted to the Astro2020 Decadal Survey with particular relevance to the key science goals of SWGO, which include unveiling Galactic particle accelerators [2-10], exploring the dynamic universe [11-21], and probing physics beyond the Standard Model [22-25], are highlighted in red boldface.

TeV Halos, extended regions of TeV gamma-ray emission around middle-aged pulsars, have recently been established as a new source class in gamma-ray astronomy. These halos have been attributed to relativistic electrons and positrons that have left the acceleration region close to the pulsar and are diffusing in the surrounding medium. Measuring the morphology of TeV Halos enables for the first time a direct measurement of the electron diffusion on scales of tens of parsecs. There are hints that the presence of relativistic particles affects the diffusion rate in the pulsars' surroundings. Understanding electron diffusion is necessary to constrain the origins of the apparent `excess' of cosmic-ray positrons at tens of GeV. TeV Halos can also be used to find mis-aligned pulsars, as well as study certain properties of the Galaxy's pulsar population. Future VHE gamma-ray instruments will detect more of those TeV Halos and determine how much pulsars contribute to the observed cosmic-ray electron and positron fluxes, and how they affect diffusion in their environments.

We outline the science motivation for SGSO, the Southern Gamma-Ray Survey Observatory. SGSO will be a next-generation wide field-of-view gamma-ray survey instrument, sensitive to gamma-rays in the energy range from 100 GeV to hundreds of TeV. Its science topics include unveiling galactic and extragalactic particle accelerators, monitoring the transient sky at very high energies, probing particle physics beyond the Standard Model, and the characterization of the cosmic ray flux. SGSO will consist of an air shower detector array, located in South America. Due to its location and large field of view, SGSO will be complementary to other current and planned gamma-ray observatories such as HAWC, LHAASO, and CTA.

Only five binary systems have been found to emit at TeV energies. Each of these systems is composed of a massive O or B type star and a compact object (black hole or a pulsar). The type of compact object and the origin of the gamma-ray emission is unknown for most of these systems. Extending spectral observations to higher energies can help disentangle the nature of the compact object as well as the particle acceleration mechanisms present. Interestingly, the TeV emission from these systems does not always coincide with their emission in GeV or X-ray, which is how many such systems have been originally discovered. Increased coverage of these systems may allow HAWC to see precisely when in the orbit the TeV emission begins and ends. The HAWC Observatory detects TeV gamma-rays with high sensitivity, covering over two-thirds of the overhead sky every day. Applying a stacking method to known TeV binary systems can help HAWC enhance the signal from TeV binaries above the steady background from other sources in the galaxy. We will present results from this stacking analysis using 760 days of HAWC data.

The Cygnus arm of our galaxy is a source-rich and complex region hosting multiple gamma-ray source types such as pulsar wind nebulae (PWN), supernova remnants, binary systems, and star clusters. The High Altitude Water Cherenkov (HAWC) observatory has been collecting data continuously since 2015 and has reported five sources within the Cygnus region. Several other instruments have also observed gamma-ray sources in this region. For instance, Fermi-LAT found gamma-ray emission at GeV energies due to a Cocoon of freshly accelerated cosmic rays, which is co-located with a known PWN TeV 2032+4130 seen by several TeV gamma-ray observatories. TeV J2032+4130 is likely powered by the pulsar PSR J2032+4127 based on the multi-wavelength observation and asymmetric morphology reported by VERITAS. The study of HAWC data will provide more information regarding the morphology, emission origin, and the correlation with the GeV emission. This presentation will discuss the analysis of data collected with the HAWC instrument and the Fermi-LAT and the results obtained to provide a deeper understanding of the Cygnus Cocoon across five decades of energy range.

Gamma-ray emission from large structures is useful for tracing the propagation and distribution of cosmic rays throughout our Galaxy. For example, the search for gamma-ray emission from Giant Molecular Clouds may allow us to probe the flux of cosmic rays in distant galactic regions and to compare it with the flux measured at Earth. Also, the composition of the cosmic rays can be measured by separating the gamma-ray emission from hadronic or leptonic processes. In the case of emission from the Fermi Bubbles specifically, constraining the mechanism of gamma-ray production can point to their origin. HAWC possesses a large field of view and good sensitivity to spatially extended sources, which currently makes it the best suited ground-based observatory to detect extended regions. We will present preliminary results on the search of gamma-ray emission from Molecular Clouds, as well as upper limits on the differential flux from the Fermi Bubbles.

Studying gamma-ray emission by Galactic objects is key to understanding the origins and acceleration mechanisms of Galactic cosmic ray electrons and hadrons. The HAWC observatory provides an unprecedented view of the gamma-ray sky at TeV energies and is particularly suited for the study of Galactic objects. However, the interpretation of the measured data poses several challenges. The high density of sources and source candidates can cause source confusion and make it harder to disentangle the origin of the emission. The relatively low angular resolution of HAWC, compared to instruments in optical or radio wavelengths, can further cause the emission of neighboring sources to bleed into each other or even make them look like one extended source. On the other hand, with its wide field of view, HAWC is uniquely suited for the study of extended sources. However, this requires the simultaneous modeling of both their morphology and emission spectrum. Joint likelihood fits to data taken over a larger range of energies can help overcome these challenges and achieve the full potential of the HAWC detector. In this presentation, we will discuss how systematic uncertainties related to joint likelihood fits can affect the measurements.

A number of Galactic sources emit GeV-TeV gamma rays that are produced through leptonic and/or hadronic mechanisms. Spectral analysis in this energy range is crucial in order to understand the emission mechanisms. The HAWC Gamma-Ray Observatory, with a large field of view and location at $19^\circ$ N latitude, is surveying the Galactic Plane from high Galactic longitudes down to near the Galactic Center. Data taken with partially-constructed HAWC array in 2013-2014 exhibit TeV gamma-ray emission along the Galactic Plane. A high-level analysis likelihood framework for HAWC, also presented at this meeting, has been developed concurrently with the Multi-Mission Maximum Likelihood (3ML) architecture to deconvolve the Galactic sources and to perform multi-instrument analysis. It has been tested on early HAWC data and the same method will be applied on HAWC data with the full array. I will present preliminary results on Galactic sources from TeV observations with HAWC and from joint analysis on Fermi and HAWC data in GeV-TeV energy range.

The High Altitude Water Cherenkov (HAWC) observatory is a wide field-of-view detector sensitive to gamma rays of 100 GeV to a few hundred TeV. Located in central Mexico at 19 degrees North latitude and 4100 m above sea level, HAWC will observe gamma rays and cosmic rays with an array of water Cherenkov detectors. The full HAWC array is scheduled to be operational in Spring 2015. In this paper, we study the HAWC sensitivity to the gamma-ray signatures of high-mass (multi- TeV) dark matter annihilation. The HAWC observatory will be sensitive to diverse searches for dark matter annihilation, including annihilation from extended dark matter sources, the diffuse gamma-ray emission from dark matter annihilation, and gamma-ray emission from non-luminous dark matter subhalos. Here we consider the HAWC sensitivity to a subset of these sources, including dwarf galaxies, the M31 galaxy, the Virgo cluster, and the Galactic center. We simulate the HAWC response to gamma rays from these sources in several well-motivated dark matter annihilation channels. If no gamma-ray excess is observed, we show the limits HAWC can place on the dark matter cross-section from these sources. In particular, in the case of dark matter annihilation into gauge bosons, HAWC will be able to detect a narrow range of dark matter masses to cross-sections below thermal. HAWC should also be sensitive to non-thermal cross-sections for masses up to nearly 1000 TeV. The constraints placed by HAWC on the dark matter cross-section from known sources should be competitive with current limits in the mass range where HAWC has similar sensitivity. HAWC can additionally explore higher dark matter masses than are currently constrained.

The Cygnus region is a very bright and complex portion of the TeV sky, host to unidentified sources and a diffuse excess with respect to conventional cosmic-ray propagation models. Two of the brightest TeV sources, MGRO J2019+37 and MGRO J2031+41, are analyzed using Milagro data with a new technique, and their emission is tested under two different spectral assumptions: a power law and a power law with an exponential cutoff. The new analysis technique is based on an energy estimator that uses the fraction of photomultiplier tubes in the observatory that detect the extensive air shower. The photon spectrum is measured in the range 1 to 200 TeV using the last 3 years of Milagro data (2005-2008), with the detector in its final configuration. MGRO J2019+37 is detected with a significance of 12.3 standard deviations ($\sigma$), and is better fit by a power law with an exponential cutoff than by a simple power law, with a probability $>98$% (F-test). The best-fitting parameters for the power law with exponential cutoff model are a normalization at 10 TeV of $7^{+5}_{-2}\times10^{-10}$ $\mathrm{s^{-1}\: m^{-2}\: TeV^{-1}}$, a spectral index of $2.0^{+0.5}_{-1.0}$ and a cutoff energy of $29^{+50}_{-16}$ TeV. MGRO J2031+41 is detected with a significance of 7.3$\sigma$, with no evidence of a cutoff. The best-fitting parameters for a power law are a normalization of $2.4^{+0.6}_{-0.5}\times10^{-10}$ $\mathrm{s^{-1}\: m^{-2}\: TeV^{-1}}$ and a spectral index of $3.08^{+0.19}_{-0.17}$. The overall flux is subject to an $\sim$30% systematic uncertainty. The systematic uncertainty on the power law indices is $\sim$0.1. A comparison with previous results from TeV J2032+4130, MGRO J2031+41 and MGRO J2019+37 is also presented.

The Crab Nebula was detected with the Milagro experiment at a statistical significance of 17 standard deviations over the lifetime of the experiment. The experiment was sensitive to approximately 100 GeV - 100 TeV gamma ray air showers by observing the particle footprint reaching the ground. The fraction of detectors recording signals from photons at the ground is a suitable proxy for the energy of the primary particle and has been used to measure the photon energy spectrum of the Crab Nebula between ~1 and ~100 TeV. The TeV emission is believed to be caused by inverse-Compton up-scattering scattering of ambient photons by an energetic electron population. The location of a TeV steepening or cutoff in the energy spectrum reveals important details about the underlying electron population. We describe the experiment and the technique for distinguishing gamma-ray events from the much more-abundant hadronic events. We describe the calculation of the significance of the excess from the Crab and how the energy spectrum is fit. The fit is consistent with values measured by IACTs between 1 and 20 TeV. Fixing the spectral index to values that have been measured below 1 TeV by IACT experiments (2.4 to 2.6), the fit to the Milagro data suggests that Crab exhibits a spectral steepening or cutoff between about 20 to 40 TeV.

We report studies of ultra-high energy cosmic ray composition via analysis of depth of airshower maximum (Xmax), for airshower events collected by the High Resolution Fly's Eye (HiRes) observatory. The HiRes data are consistent with a constant elongation rate d<Xmax>/d(log(E)) of 47.9 +- 6.0 (stat.) +- 3.2 (syst.)g/cm^2/decade for energies between 1.6 EeV and 63 EeV, and are consistent with a predominantly protonic composition of cosmic rays when interpreted via the QGSJET01 and QGSJET-II high-energy hadronic interaction models. These measurements constrain models in which the galactic-to-extragalactic transition is the cause of the energy spectrum "ankle' at 4 EeV.

The High-Altitude Water Cherenkov (HAWC) Experiment is a second-generation highsensitivity gamma-ray and cosmic-ray detector that builds on the experience and technology of the Milagro observatory. Like Milagro, HAWC utilizes the water Cherenkov technique to measure extensive air showers. Instead of a pond filled with water (as in Milagro) an array of closely packed water tanks is used. The event direction will be reconstructed using the times when the PMTs in each tank are triggered. Therefore, the timing calibration will be crucial for reaching an angular resolution as low as 0.25 degrees.We propose to use a laser calibration system, patterned after the calibration system in Milagro. Like Milagro, the HAWC optical calibration system will use ~1 ns laser light pulses. Unlike Milagro, the PMTs are optically isolated and require their own optical fiber calibration. For HAWC the laser light pulses will be directed through a series of optical fan-outs and fibers to illuminate the PMTs in approximately one half of the tanks on any given pulse. Time slewing corrections will be made using neutraldensity filters to control the light intensity over 4 orders of magnitude. This system is envisioned to run continuously at a low rate and will be controlled remotely. In this paper, we present the design of the calibration system and first measurements of its performance.

We present the result of a search of the Milagro sky map for spatial correlations with sources from a subset of the recent Fermi Bright Source List (BSL). The BSL consists of the 205 most significant sources detected above 100 MeV by the Fermi Large Area Telescope. We select sources based on their categorization in the BSL, taking all confirmed or possible Galactic sources in the field of view of Milagro. Of the 34 Fermi sources selected, 14 are observed by Milagro at a significance of 3 standard deviations or more. We conduct this search with a new analysis which employs newly-optimized gamma-hadron separation and utilizes the full 8-year Milagro dataset. Milagro is sensitive to gamma rays with energy from 1 to 100 TeV with a peak sensitivity from 10-50 TeV depending on the source spectrum and declination. These results extend the observation of these sources far above the Fermi energy band. With the new analysis and additional data, multi-TeV emission is definitively observed associated with the Fermi pulsar, J2229.0+6114, in the Boomerang Pulsar Wind Nebula (PWN). Furthermore, an extended region of multi-TeV emission is associated with the Fermi pulsar, J0634.0+1745, the Geminga pulsar.

Air-fluorescence detectors such as the High Resolution Fly's Eye (HiRes) detector are very sensitive to upward-going, Earth-skimming ultrahigh energy electron-neutrino-induced showers. This is due to the relatively large interaction cross sections of these high-energy neutrinos and to the Landau-Pomeranchuk-Migdal (LPM) effect. The LPM effect causes a significant decrease in the cross sections for bremsstrahlung and pair production, allowing charged-current electron-neutrino-induced showers occurring deep in the Earth's crust to be detectable as they exit the Earth into the atmosphere. A search for upward-going neutrino-induced showers in the HiRes-II monocular dataset has yielded a null result. From an LPM calculation of the energy spectrum of charged particles as a function of primary energy and depth for electron-induced showers in rock, we calculate the shape of the resulting profile of these showers in air. We describe a full detector Monte Carlo simulation to determine the detector response to upward-going electron-neutrino-induced cascades and present an upper limit on the flux of electron-neutrinos.

Measurements are reported of the yield and spectrum of fluorescence, excited by a 28.5 GeV electron beam, in air at a range of pressures of interest to ultra-high energy cosmic ray detectors. The wavelength range was 300 - 420 nm. System calibration has been performed using Rayleigh scattering of a nitrogen laser beam. In atmospheric pressure dry air at 304 K the yield is 20.8 +/- 1.6 photons per MeV.

Measurements are reported on the fluorescence of air as a function of depth in electromagnetic showers initiated by bunches of 28.5 GeV electrons. The light yield is compared with the expected and observed depth profiles of ionization in the showers. It validates the use of atmospheric fluorescence profiles in measuring ultra high energy cosmic rays.