The Alice UV spectrograph aboard NASA's New Horizons mission is sensitive to MeV electrons that penetrate the instrument's thin aluminum housing and interact with its microchannel plate detector. We have searched for penetrating electrons at heliocentric distance of 2-45 AU, finding no evidence of discrete events outside of the Jovian magnetosphere. However, we do find a gradual long-term increase in the Alice instrument's global dark count rate at a rate of ~1.5% per year, which may be caused by a heightened gamma-ray background from aging of the spacecraft's radioisotope thermoelectric generator fuel. If this hypothesis is correct, then the Alice instrument's global dark count rate should flatten and then decrease over the next 5-10 years.

The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes.

The Cold Classical Kuiper Belt, a class of small bodies in undisturbed orbits beyond Neptune, are primitive objects preserving information about Solar System formation. The New Horizons spacecraft flew past one of these objects, the 36 km long contact binary (486958) Arrokoth (2014 MU69), in January 2019. Images from the flyby show that Arrokoth has no detectable rings, and no satellites (larger than 180 meters diameter) within a radius of 8000 km, and has a lightly-cratered smooth surface with complex geological features, unlike those on previously visited Solar System bodies. The density of impact craters indicates the surface dates from the formation of the Solar System. The two lobes of the contact binary have closely aligned poles and equators, constraining their accretion mechanism.

The New Horizons spacecraft's encounter with the cold classical Kuiper belt object (486958) Arrokoth (formerly 2014 MU69) revealed a contact-binary planetesimal. We investigate how it formed, finding it is the product of a gentle, low-speed merger in the early Solar System. Its two lenticular lobes suggest low-velocity accumulation of numerous smaller planetesimals within a gravitationally collapsing, solid particle cloud. The geometric alignment of the lobes indicates the lobes were a co-orbiting binary that experienced angular momentum loss and subsequent merger, possibly due to dynamical friction and collisions within the cloud or later gas drag. Arrokoth's contact-binary shape was preserved by the benign dynamical and collisional environment of the cold classical Kuiper belt, and so informs the accretion processes that operated in the early Solar System.

On January 1st 2019, the New Horizons spacecraft flew by the classical Kuiper belt object (486958) Arrokoth (provisionally designated 2014 MU69), possibly the most primitive object ever explored by a spacecraft. The I/F of Arrokoth is analyzed and fit with a photometric function that is a linear combination of the Lommel-Seeliger (lunar) and Lambert photometric functions. Arrokoth has a geometric albedo of p_V = 0.21_(-0.04)^(+0.05) at a wavelength of 550 nm and ~0.24 at 610 nm. Arrokoth's geometric albedo is greater than the median but consistent with a distribution of cold classical Kuiper belt objects whose geometric albedos were determined by fitting a thermal model to radiometric observations. Thus, Arrokoth's geometric albedo adds to the orbital and spectral evidence that it is a cold classical Kuiper belt object. Maps of the normal reflectance and hemispherical albedo of Arrokoth are presented. The normal reflectance of Arrokoth's surface varies with location, ranging from ~0.10-0.40 at 610 nm with an approximately Gaussian distribution. Both Arrokoth's extrema dark and extrema bright surfaces are correlated to topographic depressions. Arrokoth has a bilobate shape and the two lobes have similar normal reflectance distributions: both are approximately Gaussian, peak at ~0.25 at 610 nm, and range from ~0.10-0.40, which is consistent with co-formation and co-evolution of the two lobes. The hemispherical albedo of Arrokoth varies substantially with both incidence angle and location, the average hemispherical albedo at 610 nm is 0.063 +/- 0.015. The Bond albedo of Arrokoth at 610 nm is 0.062 +/- 0.015.

The outer Solar System object (486958) Arrokoth (provisional designation 2014 MU$_{69}$) has been largely undisturbed since its formation. We study its surface composition using data collected by the New Horizons spacecraft. Methanol ice is present along with organic material, which may have formed through radiation of simple molecules. Water ice was not detected. This composition indicates hydrogenation of carbon monoxide-rich ice and/ or energetic processing of methane condensed on water ice grains in the cold, outer edge of the early Solar System. There are only small regional variations in color and spectra across the surface, suggesting Arrokoth formed from a homogeneous or well-mixed reservoir of solids. Microwave thermal emission from the winter night side is consistent with a mean brightness temperature of 29$\pm$5 K.

We present the results from four stellar occultations by (486958) Arrokoth, the flyby target of the New Horizons extended mission. Three of the four efforts led to positive detections of the body, and all constrained the presence of rings and other debris, finding none. Twenty-five mobile stations were deployed for 2017 June 3 and augmented by fixed telescopes. There were no positive detections from this effort. The event on 2017 July 10 was observed by SOFIA with one very short chord. Twenty-four deployed stations on 2017 July 17 resulted in five chords that clearly showed a complicated shape consistent with a contact binary with rough dimensions of 20 by 30 km for the overall outline. A visible albedo of 10% was derived from these data. Twenty-two systems were deployed for the fourth event on 2018 Aug 4 and resulted in two chords. The combination of the occultation data and the flyby results provides a significant refinement of the rotation period, now estimated to be 15.9380 $\pm$ 0.0005 hours. The occultation data also provided high-precision astrometric constraints on the position of the object that were crucial for supporting the navigation for the New Horizons flyby. This work demonstrates an effective method for obtaining detailed size and shape information and probing for rings and dust on distant Kuiper Belt objects as well as being an important source of positional data that can aid in spacecraft navigation that is particularly useful for small and distant bodies.

The flyby of Pluto and Charon by the New Horizons spacecraft provided high-resolution images of cratered surfaces embedded in the Kuiper belt, an extensive region of bodies orbiting beyond Neptune. Impact craters on Pluto and Charon were formed by collisions with other Kuiper belt objects (KBOs) with diameters from ~40 kilometers to ~300 meters, smaller than most KBOs observed directly by telescopes. We find a relative paucity of small craters less than approximately 13 kilometers in diameter, which cannot be explained solely by geological resurfacing. This implies a deficit of small KBOs (less than 1 to 2 kilometers in diameter). Some surfaces on Pluto and Charon are likely greater than 4 billion years old, thus their crater records provide information on the size-frequency distribution of KBOs in the early Solar System.

The deep nitrogen-covered Sputnik Planitia (SP; informal name) basin on Pluto is located very close to the longitude of Pluto's tidal axis[1] and may be an impact feature [2], by analogy with other large basins in the solar system[3,4]. Reorientation[5-7] due to tidal and rotational torques can explain SP's location, but requires it to be a positive gravity anomaly[7], despite its negative topography. Here we argue that if SP formed via impact and if Pluto possesses a subsurface ocean, a positive gravity anomaly would naturally result because of shell thinning and ocean uplift, followed by later modest N2 deposition. Without a subsurface ocean a positive gravity anomaly requires an implausibly thick N2 layer (greater than 40 km). A rigid, conductive ice shell is required to prolong such an ocean's lifetime to the present day[8] and maintain ocean uplift. Because N2 deposition is latitude-dependent[9], nitrogen loading and reorientation may have exhibited complex feedbacks[7].

The vast, deep, volatile-ice-filled basin informally named Sputnik Planum is central to Pluto's geological activity[1,2]. Composed of molecular nitrogen, methane, and carbon monoxide ices[3], but dominated by N2-ice, this ice layer is organized into cells or polygons, typically ~10-40 km across, that resemble the surface manifestation of solid state convection[1,2]. Here we report, based on available rheological measurements[4], that solid layers of N2 ice approximately greater than 1 km thick should convect for estimated present-day heat flow conditions on Pluto. More importantly, we show numerically that convective overturn in a several-km-thick layer of solid nitrogen can explain the great lateral width of the cells. The temperature dependence of N2-ice viscosity implies that the SP ice layer convects in the so-called sluggish lid regime[5], a unique convective mode heretofore not definitively observed in the Solar System. Average surface horizontal velocities of a few cm/yr imply surface transport or renewal times of ~500,000 years, well under the 10 Myr upper limit crater retention age for Sputnik Planum[2]. Similar convective surface renewal may also occur on other dwarf planets in the Kuiper belt, which may help explain the high albedos of some of them.

A unique feature of Pluto's large satellite Charon is its dark red northern polar cap. Similar colours on Pluto's surface have been attributed to organic macromolecules produced by energetic radiation processing of hydrocarbons. The polar location of this material on Charon implicates the temperature extremes that result from Charon's high obliquity and long seasons. The escape of Pluto's atmosphere provides a potential feed stock for production of complex chemistry. Gas from Pluto that is transiently cold-trapped and processed at Charon's winter pole was proposed as an explanation on the basis of an image of Charon's northern hemisphere, but not modelled quantitatively. Here we report images of the southern hemisphere illuminated by Pluto-shine and also images taken during the approach phase showing the northern polar cap over a range of longitudes. We model the surface thermal environment on Charon, the supply and temporary cold-trapping of material escaping from Pluto, and, while cold-trapped, its photolytic processing into more complex and less volatile molecules. The model results are consistent with the proposed mechanism producing the observed colour pattern on Charon.

Pluto's atmospheric haze settles out rapidly compared with geological timescales. It needs to be accounted for as a surface material, distinct from Pluto's icy bedrock and from the volatile ices that migrate via sublimation and condensation on seasonal timescales. This paper explores how a steady supply of atmospheric haze might affect three distinct provinces on Pluto. We pose the question of why they each look so different from one another if the same haze material is settling out onto all of them. Cthulhu is a more ancient region with comparatively little present-day geological activity, where the haze appears to simply accumulate over time. Sputnik Planitia is a very active region where glacial convection, as well as sublimation and condensation rapidly refresh the surface, hiding recently deposited haze from view. Lowell Regio is a region of intermediate age featuring very distinct coloration from the rest of Pluto. Using a simple model haze particle as a colorant, we are not able to match the colors in both Lowell Regio and Cthulhu. To account for their distinct colors, we propose that after arrival at Pluto's surface, haze particles may be less inert than might be supposed from the low surface temperatures. They must either interact with local materials and environments to produce distinct products in different regions, or else the supply of haze must be non-uniform in time and/or location, such that different products are delivered to different places.

In preparation for the Jan 1/2019 encounter between the New Horizons spacecraft and the Kuiper Belt object 2014 MU69, we provide estimates of the expected impact crater surface density on the Kuiper Belt object. Using the observed crater fields on Charon and Pluto down to the resolution limit of the 2015 New Horizons flyby of those bodies and estimates of the orbital distribution of the crater-forming projectiles, we calculate the number of craters per unit area formed as a function of the time a surface on 2014 MU69 has been exposed to bombardment. We find that if the shallow crater-size distribution from 1-15 km exhibited on Pluto and Charon is indeed due to the sizes of Kuiper Belt projectiles, 2014 MU69 should exhibit a surface that is only lightly cratered below 1 km scale, despite being bombarded for $\sim4$~billion years. Its surface should therefore be more clearly indicative of its accretionary environment. In addition, this object may the first observed for which the majority of the bombardment is from exogenic projectiles moving less than or near the speed of sound in the target materials, implying morphologies more akin to secondary craters elsewhere in the Solar System. Lastly, if the shallow Kuiper Belt size distribution implied from the Pluto and Charon imaging is confirmed at MU69, then we conclude that this size distribution is a preserved relic of its state $\simeq$4.5~Gyr ago and provides a direct constraint on the planetesimal formation process itself.

With the number of confirmed rocky exoplanets increasing steadily, their characterisation and the search for exoplanetary biospheres is becoming an increasingly urgent issue in astrobiology. We aim to investigate the possibility of characterising an exoplanet (in terms of habitability, geology, presence of life etc.) by studying material ejected from the surface during an impact event. For given parameters characterising the impact event, we estimate the escaping mass and assess its subsequent collisional evolution in a circumstellar orbit, assuming a Sun-like host star. We calculate the fractional luminosity of the dust as a function of time after the impact event and study its detectability with current and future instrumentation. We consider the possibility to constrain the dust composition, giving information on the geology or the presence of a biosphere. As examples, we investigate whether calcite, silica or ejected microorganisms could be detected. For a 20 km diameter impactor, we find that the dust mass escaping the exoplanet is roughly comparable to the zodiacal dust. The collisional evolution is best modelled by considering two independent dust populations, a spalled population consisting of non-melted ejecta evolving on timescales of millions of years, and dust recondensed from melt or vapour evolving on much shorter timescales. While the presence of dust can potentially be inferred with current telescopes, studying its composition requires advanced instrumentation not yet available. The direct detection of biological matter turns out to be extremely challenging. Despite considerable difficulties (small dust masses, noise such as exozodiacal dust etc.), studying dusty material ejected from an exoplanetary surface might become an interesting complement to atmospheric studies in the future.

The New Horizons encounter with the cold classical Kuiper Belt object (KBO) 2014 MU69 (informally named 'Ultima Thule,' hereafter Ultima) on 1 January 2019 will be the first time a spacecraft has ever closely observed one of the free-orbiting small denizens of the Kuiper Belt. Related to but not thought to have formed in the same region of the Solar System as the comets that been explored so far, it will also be the largest, most distant, and most primitive body yet visited by spacecraft. In this letter we begin with a brief overview of cold classical KBOs, of which Ultima is a prime example. We give a short preview of our encounter plans. We note what is currently known about Ultima from earth-based observations. We then review our expectations and capabilities to evaluate Ultima's composition, surface geology, structure, near space environment, small moons, rings, and the search for activity.

Pluto's landscape is shaped by the cycles of the volatile ices covering its surface. In particular, the Sputnik Planitia (SP) ice sheet displays a large diversity of terrains, with bright and dark plains, pits, topographic depressions and evidences of recent and past glacial flows. Outside SP, New Horizons also revealed numerous N2 ice deposits, in Tombaugh Regio and at mid-northern latitudes. These observations suggest a complex history involving volatile and glacial processes on different timescales. We present numerical simulations of volatile transport on Pluto performed with a model able to simulate the N2 cycle over millions of years (Myrs), taking into account the changes of obliquity and orbital parameters as experienced by Pluto. Results show that over one obliquity cycle, the latitudes of SP between 25{\deg}S-30{\deg}N are dominated by N2 condensation, while the latitudes between 30-50{\deg}N are dominated by N2 sublimation. We find that a net amount of 1 km of ice has sublimed at the northern edge of SP during the last 2 Myrs. By comparing these results with the observed geology of SP, we can relate the formation of the pits and the brightness of the ice to the ice flux occurring at the annual timescale, while the glacial flows at its eastern edge and the erosion of the water ice mountains all around the ice sheet are related to the astronomical timescale. We also perform simulations with a glacial flow scheme which shows that SP is currently at its minimum extent. We also explore the stability of N2 ice outside SP. Results show that it is not stable at the poles but rather in the equatorial regions, in particular in depressions, where thick deposits may persist over tens of Myrs, before being trapped in SP. Finally, another key result is that the minimum and maximum surface pressures obtained over the simulated Myrs remain in the range of mm-Pa and Pa, respectively.

Several trace constituents can reach saturation vapor pressure in Pluto's present day atmosphere. As such, we describe a search for discrete could features in Pluto's atmosphere using New Horizons data obtained on 14-15 July 2015, during the Pluto flyby closest approach. We report that Pluto's present day atmosphere is at least largely (>99% by surface area) free of discrete clouds. We also report a handful of features that may plausibly be clouds, all of which were detected obliquely and at high phase observing geometry. No cloud candidates were identified away from the terminators or in low phase (backscattering geometry) images.

The exploration of the Pluto-Charon system by the New Horizons spacecraft represents the first opportunity to understand the distribution of albedo and other photometric properties of the surfaces of objects in the Solar System's "Third Zone" of distant ice-rich bodies. Images of the entire illuminated surface of Pluto and Charon obtained by the Long Range Reconnaissance Imager (LORRI) camera provide a global map of Pluto that reveals surface albedo variegations larger than any other Solar System world except for Saturn's moon Iapetus. Normal reflectances on Pluto range from 0.08-1.0, and the low-albedo areas of Pluto are darker than any region of Charon. Charon exhibits a much blander surface with normal reflectances ranging from 0.20-0.73. Pluto's albedo features are well-correlated with geologic features, although some exogenous low-albedo dust may be responsible for features seen to the west of the area informally named Tombaugh Regio. The albedo patterns of both Pluto and Charon are latitudinally organized, with the exception of Tombaugh Regio, with darker regions concentrated at the Pluto's equator and Charon's northern pole The phase curve of Pluto is similar to that of Triton, the large moon of Neptune believed to be a captured Kuiper Belt Object (KBO), while Charon's is similar to that of the Moon. Preliminary Bond albedos are 0.25+/-0.03 for Charon and 0.72+/-0.07 for Pluto. Maps of an approximation to the Bond albedo for both Pluto and Charon are presented for the first time. Our work shows a connection between very high albedo (near unity) and planetary activity, a result that suggests the KBO Eris may be currently active.

On July 14th 2015, NASA's New Horizons mission gave us an unprecedented detailed view of the Pluto system. The complex compositional diversity of Pluto's encounter hemisphere was revealed by the Ralph/LEISA infrared spectrometer on board of New Horizons. We present compositional maps of Pluto defining the spatial distribution of the abundance and textural properties of the volatiles methane and nitrogen ices and non-volatiles water ice and tholin. These results are obtained by applying a pixel-by-pixel Hapke radiative transfer model to the LEISA scans. Our analysis focuses mainly on the large scale latitudinal variations of methane and nitrogen ices and aims at setting observational constraints to volatile transport models. Specifically, we find three latitudinal bands: the first, enriched in methane, extends from the pole to 55deg N, the second dominated by nitrogen, continues south to 35deg N, and the third, composed again mainly of methane, reaches 20deg N. We demonstrate that the distribution of volatiles across these surface units can be explained by differences in insolation over the past few decades. The latitudinal pattern is broken by Sputnik Planitia, a large reservoir of volatiles, with nitrogen playing the most important role. The physical properties of methane and nitrogen in this region are suggestive of the presence of a cold trap or possible volatile stratification. Furthermore our modeling results point to a possible sublimation transport of nitrogen from the northwest edge of Sputnik Planitia toward the south.

NASA's New Horizons spacecraft has revealed the complex geology of Pluto and Charon. Pluto's encounter hemisphere shows ongoing surface geological activity centered on a vast basin containing a thick layer of volatile ices that appears to be involved in convection and advection, with a crater retention age no greater than $\approx$10 Ma. Surrounding terrains show active glacial flow, apparent transport and rotation of large buoyant water-ice crustal blocks, and pitting, likely by sublimation erosion and/or collapse. More enigmatic features include tall mounds with central depressions that are conceivably cryovolcanic, and ridges with complex bladed textures. Pluto also has ancient cratered terrains up to ~4 Ga old that are extensionally fractured and extensively mantled and perhaps eroded by glacial or other processes. Charon does not appear to be currently active, but experienced major extensional tectonism and resurfacing (probably cryovolcanic) nearly 4 billion years ago. Impact crater populations on Pluto and Charon are not consistent with the steepest proposed impactor size-frequency distributions proposed for the Kuiper belt.

Observations made during the New Horizons flyby provide a detailed snapshot of the current state of Pluto's atmosphere. While the lower atmosphere (at altitudes <200 km) is consistent with ground-based stellar occultations, the upper atmosphere is much colder and more compact than indicated by pre-encounter models. Molecular nitrogen (N$_2$) dominates the atmosphere (at altitudes <1800 km or so), while methane (CH$_4$), acetylene (C$_2$H$_2$), ethylene (C$_2$H$_4$), and ethane (C$_2$H$_6$) are abundant minor species, and likely feed the production of an extensive haze which encompasses Pluto. The cold upper atmosphere shuts off the anticipated enhanced-Jeans, hydrodynamic-like escape of Pluto's atmosphere to space. It is unclear whether the current state of Pluto's atmosphere is representative of its average state--over seasonal or geologic time scales.

The New Horizons mission has provided resolved measurements of Pluto's moons Styx, Nix, Kerberos, and Hydra. All four are small, with equivalent spherical diameters of $\approx$40 km for Nix and Hydra and ~10 km for Styx and Kerberos. They are also highly elongated, with maximum to minimum axis ratios of $\approx$2. All four moons have high albedos ( $\approx$50-90 %) suggestive of a water-ice surface composition. Crater densities on Nix and Hydra imply surface ages $\gtrsim$ 4 Ga. The small moons rotate much faster than synchronous, with rotational poles clustered nearly orthogonal to the common pole directions of Pluto and Charon. These results reinforce the hypothesis that the small moons formed in the aftermath of a collision that produced the Pluto-Charon binary.

The New Horizons spacecraft mapped colors and infrared spectra across the encounter hemispheres of Pluto and Charon. The volatile ices CH$_4$, CO, and N$_2$, that dominate Pluto's surface, have complicated spatial distributions resulting from sublimation, condensation, and glacial flow acting over seasonal and geological timescales. Pluto's H$_2$O ice "bedrock" is also mapped, with isolated outcrops occurring in a variety of settings. Pluto's surface exhibits complex regional color diversity associated with its distinct provinces. Charon's color pattern is simpler, dominated by neutral low latitudes and a reddish northern polar region. Charon near infrared spectra reveal highly localized areas with strong NH$_3$ absorption tied to small craters with relatively fresh-appearing impact ejecta.

We discuss two semi-independent calibration techniques used to determine the in-flight radiometric calibration for the New Horizons' Multi-spectral Visible Imaging Camera (MVIC). The first calibration technique compares the observed stellar flux to modeled values. The difference between the two provides a calibration factor that allows the observed flux to be adjusted to the expected levels for all observations, for each detector. The second calibration technique is a channel-wise relative radiometric calibration for MVIC's blue, near-infrared and methane color channels using observations of Charon and scaling from the red channel stellar calibration. Both calibration techniques produce very similar results (better than 7% agreement), providing strong validation for the techniques used. Since the stellar calibration can be performed without a color target in the field of view and covers all of MVIC's detectors, this calibration was used to provide the radiometric keywords delivered by the New Horizons project to the Planetary Data System (PDS). These keywords allow each observation to be converted from counts to physical units; a description of how these keywords were generated is included. Finally, mitigation techniques adopted for the gain drift observed in the near-infrared detector and one of the panchromatic framing cameras is also discussed.

The Pluto system was recently explored by NASA's New Horizons spacecraft, making closest approach on 14 July 2015. Pluto's surface displays diverse landforms, terrain ages, albedos, colors, and composition gradients. Evidence is found for a water-ice crust, geologically young surface units, surface ice convection, wind streaks, volatile transport, and glacial flow. Pluto's atmosphere is highly extended, with trace hydrocarbons, a global haze layer, and a surface pressure near 10 microbars. Pluto's diverse surface geology and long-term activity raise fundamental questions about how small planets remain active many billions of years after formation. Pluto's large moon Charon displays tectonics and evidence for a heterogeneous crustal composition, its north pole displays puzzling dark terrain. Small satellites Hydra and Nix have higher albedos than expected.

N2 is abundant in Pluto's atmosphere and on its surface, but the supply is depleted by prodigious atmospheric escape. We demonstrate that cometary impacts could not have delivered enough N mass to resupply Pluto's N2 atmospheric escape over time; thus Pluto's N2 is likely endogenous, and therefore was either acquired early in its history or created by chemistry inside/on Pluto. We find that cratering could excavate a considerable amount of N2 to resupply the atmosphere against escape if the near-surface N2 reservoir is deep. However, we find that this process likely falls short of that necessary to resupply the atmosphere against escape by at least an order of magnitude. We conclude that either the escape of N2 from Pluto's atmosphere was on average much lower than the predictions for the current epoch, or that internal activity could be necessary to bring N2 to the surface and resupply escape losses. Observations made by the New Horizons spacecraft in mid-2015 will yield further constraints on the provenance and evolution of Pluto's surface and atmospheric N2, and could reveal evidence for past or present internal activity.

The identification of impact craters on planetary surfaces provides important information about their geological history. Most studies have relied on individual analysts who map and identify craters and interpret crater statistics. However, little work has been done to determine how the counts vary as a function of technique, terrain, or between researchers. Furthermore, several novel internet-based projects ask volunteers with little to no training to identify craters, and it was unclear how their results compare against the typical professional researcher. To better understand the variation among experts and to compare with volunteers, eight professional researchers have identified impact features in two separate regions of the moon. Small craters (diameters ranging from 10 m to 500 m) were measured on a lunar mare region and larger craters (100s m to a few km in diameter) were measured on both lunar highlands and maria. Volunteer data were collected for the small craters on the mare. Our comparison shows that the level of agreement among experts depends on crater diameter, number of craters per diameter bin, and terrain type, with differences of up to $\sim\pm45%$. We also found artifacts near the minimum crater diameter that was studied. These results indicate that caution must be used in most cases when interpreting small variations in crater size-frequency distributions and for craters $\le10$ pixels across. Because of the natural variability found, projects that emphasize many people identifying craters on the same area and using a consensus result are likely to yield the most consistent and robust information.