The mean of 830 visual magnitudes adjusted to a distance of 550 km (the operational altitude) is 4.63 +/-0.02. The data on DarkSat, the low-albedo satellite, indicate that it is fainter than the others by 1.6 magnitudes or 78%. However, there is considerable uncertainty in this value due to the small number of observations. Some satellites were observed to flare by 10,000 times their normal brightness. These statistics can serve as a baseline for assessing the reduced brightness of the VisorSat design for future Starlink satellites.

The Starlink satellites are shaped like flat panels. The flat sides face zenith and nadir during normal operations. Their brightness is determined by the product of the solar illumination on the downward facing side of the panel multiplied by the area of that side projected toward the observer on Earth. This geometry leads to a unique brightness function that is not shared by other satellites. For example, the observed brightness is very sensitive to the solar elevation angle. There are circumstances where sunlight only illuminates the upward facing side of the satellite rendering it invisible to Earth-based observers. A brightness model depending on the solar aspect and the observer aspect of the flat panel, in addition to the satellite distance, is described. Absolute brightness is the only free parameter of the model, and it is taken to be that at a distance of 1,000 km when the solar and observer factors are unity. This model has been successfully fitted to a set of observed magnitudes. The absolute visual magnitude of a Starlink satellite as determined from this fitting is 4.1 +/- 0.1. The model could be used to determine the absolute magnitude of the Starlink satellite known as Dark Sat which has a special low-albedo coating.

Improved equations for computing planetary magnitudes are reported. These formulas model V-band observations acquired from the time of the earliest filter photometry in the 1950s up to the present era. The new equations incorporate several terms that have not previously been used for generating physical ephemerides. These include the rotation and revolution angles of Mars, the sub-solar and sub-Earth latitudes of Uranus, and the secular time dependence of Neptune. Formulas for use in The Astronomical Almanac cover the planetary phase angles visible from Earth. Supplementary equations cover those phase angles beyond the geocentric limits. Geocentric magnitudes were computed over a span of at least 50 years and the results were statistically analyzed. The mean, variation and extreme magnitudes for each planet are reported. Other bands besides V on the Johnson-Cousins and Sloan photometric systems are briefly discussed. The planetary magnitude data products available from the U.S. Naval Observatory are also listed. An appendix describes source code and test data sets that are available on-line for computing planetary magnitudes according to the equations and circumstances given in this paper. The files are posted as supplementary material for this paper. They are also available at SourceForge under project https://sourceforge.net/projects/planetary-magnitudes/ under the 'Files' tab in the folder 'Ap_Mag_Current_Version'.

A new version of the Catalog containing Sloan magnitudes for the brightest stars is presented. The accuracy of the data indicates that the Catalog is a reliable source of comparison star magnitudes for astronomical photometry. Version 2 complements the APASS database of fainter stars.

Mid-IR albedo values of Neptune are derived from Spitzer Space Telescope measurements reported by Stauffer et al. (2016). The method of this derivation is described and the results indicate that the geometric albedo was about 1% or less at the time of the observations in 2016. Short-term mid-IR variability of Neptune, evidenced by the Spitzer observations themselves, indicates an albedo at 3.6 microns ranging from 0.2% to 0.6% with a mean of 0.4%. The corresponding albedos at 4.5 microns are 0.7%, 1.3% and 0.9%. Furthermore, the 60-year history of visible-light brightness variations, which show that Neptune was significantly fainter a few decades ago, suggests that the mid-IR albedo during that earlier period of time may have been much less than 1%. The albedo values reported here can have implications for models of Neptune's atmosphere. However, the mid-IR brightness of Neptune cannot contribute very strongly to its total albedo because the Sun emits only about 2% of its flux long-ward of 3 micron. By contrast, the Sun emits 42% of its flux at visible and near-UV wavelengths where the planet's albedo is in the tens of percents.

Published reflectance data covering several different wavelength intervals has been combined and analyzed in order to determine the spherical bolometric albedo of Mercury. The resulting value of 0.088 +/- 0.003 spans wavelengths from 0 to 4 {\mu}m which includes over 99% of the solar flux. This bolometric result is greater than the value determined between 0.43 and 1.01 {\mu}m by Domingue et al. (2011, Planet. Space Sci., 59, 1853-1872). The difference is due to higher reflectivity at wavelengths beyond 1.01 {\mu}m. The average effective blackbody temperature of Mercury corresponding to the newly determined albedo is 436.3 K. This temperature takes into account the eccentricity of the planet's orbit (M\'endez and Rivera-Valet\'in. 2017. ApJL, 837, L1).

Saturn's geometric magnitude and albedo were determined by Mallama et al. (2017, Icarus 282, 19-33). That analysis depended largely on photometric data obtained when the ring system interfered with observations of the globe. In this study the magnitude and albedo are re-evaluated based on spectrophotometry obtained during a ring-plane crossing (Karkoschka, 1998, Icarus, 133, 134-146). A better correction for the illumination phase angle of Saturn's globe is also derived and it is applied to the spectroscopy. The resulting values for most of the 12 Johnson-Cousins and Sloan wavelength bands agree reasonably well with those reported by Mallama et al. but some significant differences are noted. The revised albedo value for the Rc band reduces the discrepancy with that inferred from the work of Dyudina et al. (2005, ApJ 618, 973-986).

Complete sets of reference magnitudes in all 7 Johnson-Cousins bands (U, B, V, R, I, Rc and Ic) and the 5 principal Sloan bands (u', g', r', i', and z') are presented for the 8 planets. These data are accompanied by illumination phase functions and other formulas which characterize the instantaneous brightness of the planets. The main source of Johnson-Cousins magnitudes is a series of individualized photometric studies reported in recent years. Gaps in that dataset were filled with magnitudes synthesized in this study from published spectrophotometry. The planetary Sloan magnitudes, which are established here for the first time, are an average of newly recorded Sloan filter photometry, synthetic magnitudes and values transformed from the Johnson-Cousins system. Geometric albedos derived from these two sets of magnitudes are consistent within each photometric system and between the systems for all planets and in all bands. This consistency validates the albedos themselves as well as the magnitudes from which they were derived. In addition, a quantity termed the delta stellar magnitude is introduced to indicate the difference between the magnitude of a planet and that of its parent star. A table of these delta values for exo-planets possessing a range of physical characteristics is presented. The delta magnitudes are for phase angle 90 degrees where a planet is near the greatest apparent separation from its star. This quantity may be useful in exo-planet detection and observation strategies when an estimate of the signal-to-noise ratio is needed. Likewise, the phase curves presented in this paper can be used for characterizing exo-planets. Finally, magnitudes for the proposed Planet Nine are estimated, and we note that P9 may be especially faint at red and near-IR wavelengths.

Neptune has brightened by more than 10% during the past several decades. We report on the analysis of published Johnson-Cousins B and V magnitudes dating back to 1954 along with new U, B, V, R, Rc, I and Ic photometry that we recorded during the past 24 years. Characteristic magnitudes, colors and albedos in all seven band-passes are derived from the ensemble of data. Additionally, 25 spectra spanning 26 hours of observation on 5 nights are analyzed. The spectrophotometry demonstrates that planetary flux and albedo is inversely related to the equivalent widths of methane bands. We attribute the changes in band strength, flux and albedo to the high altitude clouds which rotate across the planet's visible disk. Bright clouds increase albedo and flux while reducing methane absorption. Synthetic V magnitudes derived from the spectroscopy also agree closely with the photometric quantities, which cross-validates the two techniques. The spectroscopic and photometric results are discussed within the framework of the secular brightness variations of Neptune.

Uranus is fainter when the Sun and Earth are near its equatorial plane than when they are near the projection of its poles. The average of the absolute values of the sub-Earth and sub-Sun latitudes (referred to as the sub-latitude here) is used to quantify this dependency. The rates of change of magnitude with sub-latitude for four of the Johnson-Cousins band-passes are B-band, -0.48 +/- 0.11 milli-magnitudes per degree; V-band, -0.84 +/- 0.04 ; R-band, -5.33 +/- 0.30; and I-band -2.79 +/- 0.41. Evaluated over the range of observed sub-latitudes, the blue flux changes by a modest 3% while the red flux varies by a much more substantial 30%. These disk-integrated variations are consistent with the published brightness characteristics of the North and South Polar Regions, with the latitudinal distribution of methane and with a planetary hemispheric asymmetry. Reference magnitudes and colors are also reported along with geometric albedos for the seven Johnson-Cousins band-passes.

The brightness of cloud bands on Jupiter and Saturn as a function of latitude is reported. Bright Jovian bands near the equator are located in regions of anti-cyclonic circulation of the atmosphere. By contrast, bright equatorial bands on Saturn are associated with cyclonic motion. Modern definitions of the cloud band terms 'zone' and 'belt' are distinguished from their old meanings.