Wright et al. 2018 (W2018) responded to certain criticisms of the NEOWISE results raised by Myhrvold (2018a, 2018b) in a preprint submitted for publication to Icarus. This paper responds to W2018 and has been updated to address comments from two rounds of peer review at Icarus. Points raised by Wright et al. during peer review and in recent publications are also addressed. The response here presents issues of broad interest to the planetary science community. W2018 acknowledges some of the most serious issues identified by Myhrvold, including the importance of Kirchhoff's law and the misrepresentation of diameters copied from previous radar, occultation, and spacecraft observations as NEOWISE thermal modeling results. W2018 confirms copying ROS diameters for 117 asteroids but does not report actual modeled diameters or the identity of the asteroids involved. W2018 confirms that many NEOWISE "fit" curves do not pass near the data points they claimed to fit. W2018 reveals more serious problems with NEOWISE results than had previously been understood. Wright et al. belatedly disclose a software bug, found and corrected in 2011, that apparently corrupted the vast majority of NEOWISE results, which have yet to be corrected. W2018 asserts that the bug had rare and small impacts on results but elsewhere provides analysis that directly contradicts this claim. We learn from W2018 that H values used by NEOWISE were not taken from the literature, as claimed, but rather derived by fitting data, a method inconsistent with prior NEOWISE papers and practice in thermal modeling. The systematic underestimation of WISE observation errors shown in Myhrvold 2018b is confirmed with a new analysis of repeated observation of calibration stars. Each argument in W2018 is shown here to be mistaken. Ultimately, W2018 reveals why the results of the empirical examination that Myhrvold published were important.

We publicly release ATM, a Python package designed to model asteroid flux measurements to estimate an asteroid's size, surface temperature distribution, and emissivity. The full multi-dimensional posterior pdf is found using Markov Chain Monte Carlo. Data files with $\sim$ 2.5 million WISE flux measurements for $\sim$ 150,000 asteroids and additional MPC data are also included with the package, as well as Python Jupyter Notebooks with examples of analysis. The entirety of the analysis presented here, including all the figures, tables, and catalogs, can be easily reproduced with these publicly released Notebooks. We show that ATM can match the best-fit size estimates for well-observed asteroids published in 2016 by the NEOWISE team (Mainzer et al. 2016) with a sub-percent bias and a scatter of only 6%. We estimate that the accuracy of WISE-based asteroid size estimates is approximately in the range of 15-20% for most objects. We also study optical data collected by the Sloan Digital Sky Survey (SDSS) and show that correlations of optical colors and WISE-based best-fit model parameters indicate robustness of the latter. Our analysis also gives support to the claim by Harris & Drube (2014) that candidate metallic asteroids can be selected using the best-fit temperature parameter and infrared albedo. We investigate a correlation between SDSS colors and optical albedo derived using WISE-based size estimates and show that this correlation can be used to estimate asteroid sizes with optical data alone, with a precision of about 21% relative to WISE-based size estimates. After accounting for systematic errors, the difference in accuracy between infrared and optical color-based size estimates becomes less than a factor of two. (abridged)

This study addresses thermal modeling of asteroids with a new derivation of the Near Earth Asteroid Thermal (NEATM) model which correctly accounts for the presence of reflected sunlight in short wave IR bands. Kirchhoff's law of thermal radiation applies to this case and has important implications. New insight is provided into the eta parameter in the NEATM model and it is extended to thermal models besides NEATM. The role of surface material properties on eta is examined using laboratory spectra of meteorites and other asteroid compositional proxies; the common assumption that emissivity e = 0.9 in asteroid thermal models may not be justified and can lead to misestimating physical parameters. In addition, indeterminacy in thermal modeling can limit its ability to uniquely determine temperature and other physical properties. A new curve fitting approach allows thermal modeling to be done independent of visible band observational parameters such as the absolute magnitude H. These new thermal modeling techniques are applied to observational data for selected asteroids from the WISE/NEOWISE mission. The previous NEOWISE analysis assumes Kirchhoff's law does not apply. It also deviates strongly from established statistical practice and systematically underestimates the sampling error inherent in observing potentially irregular asteroids from a finite sample of observations. As a result, the new analysis finds asteroid diameter and other physical properties that have large differences from published NEOWISE results, with greatly increased error estimates. NEOWISE results have a claimed +/-10% accuracy for diameter estimates, but this is unsupported. ABSTRACT CONTINUED IN PDF...

Multiple terrestrial and space-based telescopes have been proposed for detecting and tracking near-Earth objects (NEOs). Detailed simulations of the search performance of these systems have used complex computer codes that are not widely available, which hinders accurate cross-comparison of the proposals and obscures whether they have consistent assumptions. Moreover, some proposed instruments would survey infrared (IR) bands, whereas others would operate in the visible band, and differences among asteroid thermal and visible light models used in the simulations further complicate like-to-like comparisons. I use simple physical principles to estimate basic performance metrics for the ground-based Large Synoptic Survey Telescope and three space-based instruments - Sentinel, NEOCam, and a Cubesat constellation. The performance is measured against two different NEO distributions, the Bottke et al. distribution of general NEOs, and the Veres et al. distribution of Earth impacting NEO. The results of the comparison show simplified relative performance metrics, including the expected number of NEOs visible in the search volumes and the initial detection rates expected for each system. Although these simplified comparisons do not capture all of the details, they give considerable insight into the physical factors limiting performance. Multiple asteroid thermal models are considered, including FRM, NEATM, and a new generalized form of FRM (GFRM). I describe issues with how IR albedo and emissivity have been estimated in previous studies, which may render them inaccurate. A thermal model for tumbling asteroids is also developed and suggests that tumbling asteroids may be surprisingly difficult for IR telescopes to observe.