Cosmic rays (CRs) and magnetic fields may be dynamically important in driving large-scale galactic outflows from rapidly star-forming galaxies. We construct two-dimensional axisymmetric models of the local starburst and super-wind galaxy M82 using the CR propagation code GALPROP. Using prescribed gas density and magnetic field distributions, wind profiles, CR injection rates, and stellar radiation fields, we simultaneously fit both the integrated gamma-ray emission and the spatially-resolved multi-frequency radio emission extended along M82's minor axis. We explore the resulting constraints on the gas density, magnetic field strength, CR energy density, and the assumed CR advection profile. In accord with earlier one-zone studies, we generically find low central CR pressures, strong secondary electron/positron production, and an important role for relativistic bremsstrahlung losses in shaping the synchrotron spectrum. We find that the relatively low central CR density produces CR pressure gradients that are weak compared to gravity, strongly limiting the role of CRs in driving M82's fast and mass-loaded galactic outflow. Our models require strong magnetic fields and advection speeds of order ~1000 km/s on kpc scales along the minor axis in order to reproduce the extended radio emission. Degeneracies between the controlling physical parameters of the model and caveats to these findings are discussed.

Milagro observations have found bright, diffuse TeV emission concentrated along the galactic plane of the Milky Way. The intensity and spectrum of this emission is difficult to explain with current models where gamma-ray production is dominated by hadronic mechanisms, and has been named the "TeV excess". We show that TeV emission from pulsars naturally explains this excess. In particular, recent observations have detected "TeV halos" surrounding pulsars that are either nearby or particularly luminous. Here, we show that the full population of Milky Way pulsars will produce diffuse TeV emission concentrated along the Milky Way plane. The total gamma-ray flux from TeV halos is expected to exceed the hadronic gamma-ray flux at energies above ~500 GeV. Moreover, the spectrum and intensity of TeV halo emission naturally matches the TeV excess. If this scenario is common to all galaxies, it will decrease the contribution of star-forming galaxies to the IceCube neutrino flux. Finally, we show that upcoming HAWC observations will resolve a significant fraction of the TeV excess into individual TeV halos, conclusively confirming, or ruling out, this model.