Solar Phase Angle Research Articles

Combined spin orientation and phase function of asteroids

Context. Large sky surveys provide numerous non-targeted observations of small bodies of the Solar System. The upcoming LSST of the Vera C. Rubin Observatory will be the largest source of small body photometry in the next decade. With non-coordinated epochs of observation, colors – and therefore taxonomy and composition – can only be computed by comparing absolute magnitudes obtained in each filter by solving the phase function (evolution of brightness of the small body against the solar phase angle). Current models in use in the community (HG, HG12*, and HG1G2), however, fail to reproduce the long-term photometry of many targets due to the change in the aspect angle between apparitions. Aims. We aim to derive a generic yet simple phase function model accounting for the variable geometry of the small bodies over multiple apparitions. Methods. As a spinoff of the HG1 G2 model, we propose the sHG1G2 phase function model in which we introduce a term describing the brightness changes due to spin orientation and polar oblateness. We applied this new model to 13 245 908 observations of 122 675 Solar System objects (SSOs). These observations were acquired in the g and r filters with the Zwicky Transient Facility between November 1, 2019 and December 1, 2023. We retrieved them and implemented the new sHG1G2 model in FINK, a broker of alerts designed for the LSST. Results. The sHG1G2 model leads to smaller residuals than other phase function models, providing a better description of the photometry of asteroids. We determined the absolute magnitude, H, and phase function coefficients (G1, G2) in each filter, the spin orientation (α0, δ0), and the polar-to-equatorial oblateness, R, for 95 593 SSOs, which constitutes about a tenfold increase in the number of characterized objects compared to the current census. Conclusions. The application of the sHG1 G2 model to ZTF alert data using the FINK broker shows that the model is appropriate for extracting physical properties of asteroids from multi-band and sparse photometry, such as the forthcoming LSST survey.

Multicolour photometry of LEO mega-constellations Starlink and OneWeb

ABSTRACT The development of low earth orbit (LEO) mega-constellation fundamentally threatens ground-based optical astronomical observations. To study the photometric properties of the LEO mega-constellations, we used the Xinglong 50 cm telescope to conduct a large-sample, high-precision, and multicolour target-tracking photometry of two typical LEO mega-constellations: Starlink and OneWeb. Over a three-month observation period starting on 2022 January 1, we collected 1447 light curves of 404 satellites in four typical versions: Starlink v1.0, DarkSat, VisorSat, Starlink v1.5, and OneWeb. According to data statistics, Starlink v1.0 has the smallest median magnitude at clear and Sloan Digital Sky Survey gri band, and OneWeb is the dimmest bus. The brightness of Starlink v1.5 is slightly brighter than VisorSat. We construct a detailed photometric model with solar phase angle variations by calculating the illumination-visibility geometry based on the orbital parameters. Our data analysis shows that the solar phase angle is the significant characteristic which influencing Starlink satellites’ brightness, but it is not sensitive to OneWeb satellites. VisorSat and Starlink v1.5 versions, which are equipped with deployable visors, have significantly reduced scattered light compared to the previous Starlink v1.0 version. The multiband LOWESS and colour index are analysed in characterizing the energy and colour features of LEO mega-constellation satellites. This work found that the proportion of scattered sunlight mitigation achieved with VisorSat and Starlink v1.5 was 55.1 and 40.4 per cent, respectively. The colour index of different buses shows an evident clustering feature. Our observation and analysis could provide valuable quantitative data and photometric models, which can contribute to assessing the impact of LEO mega-constellations on astronomical observations.

3 μm Phase Curves of Main-belt Asteroids from NEOWISE Photometry

The shapes of asteroid phase curves are influenced by the physical properties of asteroid surfaces. The variation of an asteroid’s brightness as a function of the solar phase angle can tell us about surface properties such as grain size distribution, roughness, porosity, and composition. Phase curves are traditionally derived from photometric observations at visible wavelengths, but phase curves using infrared data can also provide useful information about an asteroid surface. Using photometric observations centered near ∼3.4 μm from the W1 band of the Near-Earth Object Wide-field Infrared Survey Explorer mission, we construct thermally and rotationally corrected infrared phase curves for a sample of main-belt asteroids, which includes asteroids observed by the AKARI satellite, as well as subsets of the Themis and Flora dynamical families. We calculate the linear slope of the phase curves as a measure of their shape and compare W1 phase slopes to band depths of absorption features associated with hydrated materials, spectral slopes, visible albedos, W1 albedos, and diameters. We observe a steepening of the W1 phase slope of C-type asteroids with increasing 2.7 μm band depth but little correlation between the phase slope and 3 μm band depth or 3 μm spectral slope. The C-types in our sample exhibit steeper average W1 phase slopes than M- or S-types, similar to visible-light phase slopes. We also observe steeper W1 phase slopes for smaller-diameter objects within the Themis family and explore comparisons to Jupiter-family comets in phase slope versus albedo space.

Optical Imaging Polarimetry of Comet 21P/Giacobini–Zinner during Its 2018 Apparition

We conducted polarimetric observations of comet 21P/Giacobini–Zinner in the optical wavelength region using the polarimetric imager PICO equipped with the 50 cm telescope for Public Outreach at the Mitaka Campus of the National Astronomical Observatory of Japan on UT 2018 September 16, at a solar phase angle of 77° for the comet. We used Gunn’s i′-band filter to cover the wavelength region where sunlight reflected by dust grains dominates cometary spectra. The intensity map of comet 21P/Giacobini–Zinner showed an elongated coma structure along the east–west direction. The linear polarization degree map was almost uniform and showed no evidence of jet or arc structures in the inner coma. The linear polarization degree was 22% in the coma, which is consistent with other optical polarization observations of 21P/Giacobini–Zinner for this and previous apparitions reported in the literature. The measured polarization angles (179°) were nearly perpendicular to the position angle of the scattering plane of 271°. There is no evidence of significant changes in the dust properties (composition, porosity, and size distribution) in the comet from its 1985 to 2018 apparitions. The spectral gradient of the linear polarization degree for the comet was nearly flat in the wavelength range of ∼510–830 nm at a solar phase angle of ∼77° during the 2018 apparition.

Multicolor and multi-spot observations of Starlink’s Visorsat

Abstract This study provides the results of simultaneous multicolor observations for the first Visorsat (STARLINK-1436) and the ordinary Starlink satellite (STARLINK-1113) in the U, B, V, g′, r, i, RC, IC, $z$, J, H, and Ks bands to quantitatively investigate the extent to which Visorsat reduces its reflected light. Our results are as follows: (1) in most cases, Visorsat is fainter than STARLINK-1113, and the sunshade on Visorsat therefore contributes to the reduction of the reflected sunlight; (2) the magnitude at 550 km altitude (normalized magnitude) of both satellites often reaches the naked-eye limiting magnitude (<6.0); (3) from a blackbody radiation model of the reflected flux, the peak of the reflected components of both satellites is around the $z$ band; and (4) the albedo of the near-infrared range is larger than that of the optical range. Under the assumption that Visorsat and STARLINK-1113 have the same reflectivity, we estimate the covering factor, Cf, of the sunshade on Visorsat, using the blackbody radiation model: the covering factor ranges from 0.18 ≤ Cf ≤ 0.92. From the multivariable analysis of the solar phase angle (Sun–target–observer), the normalized magnitude, and the covering factor, the phase angle versus covering factor distribution presents a moderate anti-correlation between them, suggesting that the magnitudes of Visorsat depend not only on the phase angle but also on the orientation of the sunshade along our line of sight. However, the impact on astronomical observations from Visorsat-designed satellites remains serious. Thus, new countermeasures are necessary for the Starlink satellites to further reduce reflected sunlight.

A new constraint on HCl abundance at the cloud top of Venus

We provided a new constraint on HCl abundance at the cloud top of Venus by infrared spectroscopy using a cross-dispersed high-resolution echelle spectrograph, iSHELL, mounted on the NASA Infrared Telescope Facility (IRTF). This study investigated the inconsistency in HCl abundance between previous ground-based observations and Venus Express solar occultation measurements. Venusian dayside observations at a solar phase angle of ∼90° were conducted during August 6–7, 2018, and August 18–20, 2020 (UT), when the Venusian afternoon and morning sides were visible, respectively. The high spectral resolving power of ∼80,000 and large Doppler shift (∼13 km/s) enabled the measurement of Venusian lines with less contamination by terrestrial lines. We analyzed the H35Cl P(5) and H37Cl P(6) lines at 2775.8 and 2750.1 cm−1, respectively, in the 1–0 band, together with 16O12C18O P- and R-branch lines of the 20001–00001 band, which fell in the same spectral orders as the HCl lines. The 16O12C18O lines were used to derive the cloud top altitude because the upper clouds significantly impacted the retrieval of HCl abundance. The cloud tops had an equatorially symmetric structure. The average altitude was 70.8 ± 0.6 km in the region equatorward of ±30° and decreased toward higher latitudes. The HCl volume mixing ratio was derived as 0.379 ± 0.013 ppm at a probing altitude of 70.6 ± 1.1 km and showed no significant latitudinal dependence within the range of ±70°. A difference of ∼0.02 ppm between 2018 and 2020 would result mainly from temporal variation. The H35Cl/H37Cl abundance ratio was 3.01 ± 0.16, with no prominent latitudinal dependence. The obtained HCl volume mixing ratio agreed with the results of previous ground-based measurements, which were approximately one order of magnitude larger than those derived from Venus Express solar occultation measurements. The systematic uncertainties in our retrieval analysis cannot explain this significant inconsistency. The impact of diffuse light produced from aerosol scattering on the retrieval method for solar occultation measurements should be investigated further to solve this issue.

The Diverse Shapes of Dwarf Planet and Large KBO Phase Curves Observed from New Horizons

Observing dwarf planets and other large Kuiper Belt objects (KBOs) from vantage points between 8 and 47 au from the Sun, NASA’s New Horizons spacecraft has found diversity in the shapes of their solar phase curves. Here we extend solar phase angle coverage of dwarf planets (136199) Eris, (136472) Makemake, and (136108) Haumea; large KBOs (28978) Ixion, (50000) Quaoar, (307261) 2002 MS4, and (556416) 2014 OE394; and Neptune’s satellite Triton to phase angles as high as α = 94° using New Horizons data and fit the resulting solar phase curves to the Hapke photometric model. When accounting for sparse α sampling, these fits yield large uncertainties in the Hapke parameters; however, opposition effect parameters are generally well constrained and suggest a significant range of regolith maturation ages among these bodies. The expanded range in α enables evaluation of Bond albedos, phase integrals, rotation curves at high α, and comparisons of the surface scattering properties of these objects with those of others in the solar system. The dwarf planets with surface compositions dominated by hypervolatiles, Eris and Makemake, and Triton (a likely former KBO) have shallower solar phase curve slopes (i.e., lower phase coefficients, higher phase integrals, and Bond albedos) than objects with volatile-poor surfaces. The total amplitude of Haumea’s rotation curve at α = 48° is Δm = 0.6 ± 0.2 mag, nearly twice that of its rotation curve measured from Earth at low phase angles. Bond albedos range from 0.037 ± 0.007 for Ixion to for Eris.

Photometry and model of near-Earth asteroid 2021 DW1 from one apparition

Aims. Very small asteroids (VSAs, objects with diameters smaller than about 150 m) can be spun up by the YORP effect to rotation periods as short as tens of seconds. This effect has been observed for many of them. It is also hypothesised, that in the same process their spin axes are asymptotically drawn to the position perpendicular to the orbital plane. So far this effect has been observed only for one VSA and needs further verification. For that, spin axes of several other VSAs should be determined by observing their brightness variations at many different positions on the sky. Methods. On 4 March 2021 at 9 UTC a 30-m in diameter near-Earth asteroid 2021 DW1 passed the Earth at a distance of 570 000 km, reaching the maximum brightness of V = 14.6 mag. We observed it photometrically from 2 March, when it was visible at V = 16.5 mag, until 7 March (V = 18.2 mag). During that time 2021 DW1 swept a 170° long arc in the northern sky, spanning solar phase angles in the range from 36° to 86°. This made it an excellent target for physical characterisation, including spin axis and shape derivation. Results. Convex inversion of the asteroid lightcurves gives a sidereal period of rotation Psid = 0.013760 ± 0.000001 h, and two solutions for the spin axis ecliptic coordinates: (A) λ1 = 57° ± 10°, β1 = 29° ± 10° and (B) λ2 = 67° ± 10°, β2 = −40° ± 10°. The magnitude-phase curve can be fitted with a standard H, G function with H = 24.8 ± 0.5 mag and an assumed G = 0.24. The asteroid colour indices are g − i = 0.79 ± 0.01 mag, and i − z = 0.01 ± 0.02 mag which indicates an S taxonomic class, with an average geometric albedo pV = 0.23 ± 0.02. The asteroid effective diameter, derived from H and pV, is Deff = 30 ± 10 m. Conclusions. It was found that the inclination of the spin axis of 2021 DW1 is not perpendicular to the orbital plane (obliquity ϵ = 54° ± 10° or ϵ = 123° ± 10°). More spin axes of VSAs should be determined to check, if 2021 DW1 is an exception or a typical case.

Polarimetric properties of the near-Sun asteroid (155140) 2005 UD in comparison with other asteroids and meteoritic samples

ABSTRACTThe investigation of asteroids near the Sun is important for understanding the final evolutionary stage of primitive Solar system objects. A near-Sun asteroid (NSA), (155140) 2005 UD, has orbital elements similar to those of (3200) Phaethon (the target asteroid for the JAXA’s DESTINY+ mission). We conducted photometric and polarimetric observations of 2005 UD and found that this asteroid exhibits a polarization phase curve similar to that of Phaethon over a wide range of observed solar phase angles (α = 20–105°) but different from those of (101955) Bennu and (162173) Ryugu (asteroids composed of hydrated carbonaceous materials). At a low phase angle (α ≲ 30°), the polarimetric properties of these NSAs (2005 UD and Phaethon) are consistent with anhydrous carbonaceous chondrites, while the properties of Bennu are consistent with hydrous carbonaceous chondrites. We derived the geometric albedo, pV ∼ 0.1 (in the range of 0.088–0.109); mean V-band absolute magnitude, HV = 17.54 ± 0.02; synodic rotational period, $T_\mathrm{rot} = 5.2388 \pm 0.0022 \, \mathrm{h}$ (the two-peaked solution is assumed); and effective mean diameter, $D_\mathrm{eff} = 1.32 \pm 0.06 \, \mathrm{km}$. At large phase angles (α ≳ 80°), the polarization phase curve are likely explained by the dominance of large grains and the paucity of small micron-sized grains. We conclude that the polarimetric similarity of these NSAs can be attributed to the intense solar heating of carbonaceous materials around their perihelia, where large anhydrous particles with small porosity could be produced by sintering.

On Charon’s Far-ultraviolet Surface Reflectance

Abstract We present the first measurements of Charon’s far-ultraviolet (FUV) surface reflectance, obtained by the Alice spectrograph on New Horizons. We find no measurable flux shortward of 1650 Å, and Charon’s geometric albedo is <0.019 (3σ) at 1600 Å. From 1650 to 1725 Å, Charon’s geometric albedo increases to 0.166 ± 0.068 and remains nearly constant until 1850 Å. As this spectral shape is characteristic of H2O ice absorption, Charon is the first Kuiper Belt object with a H2O ice surface to be detected in the FUV. Charon’s geometric albedo is ∼3.7 times lower than Enceladus’s at these wavelengths but has a very similar spectral shape. We attribute this to similarities in their surface compositions and the difference in absolute reflectivity to a high concentration or more-absorbing contaminants on Charon’s surface. Finally, we find that Charon has different solar phase behavior in the FUV than Enceladus, Mimas, Tethys, and Dione, with a stronger opposition surge than Enceladus and a shallower decline at intermediate solar phase angles than any of these Saturnian satellites.

Pluto in Glory: Discovery of Its Huge Opposition Surge

AbstractNear‐infrared observations of the Pluto Charon system were captured with the Palomar High Angular Resolution Observer (PHARO) adaptive optics (AO) system on the 200‐inch Hale telescope during the historically small solar phase angles in 2018–2019. Both objects exhibit large opposition surges of ∼30%–35% in the last half degree of solar phase angle, which is among the largest observed for icy moons and other Kuiper Belt Objects. In addition, Pluto's surge is exceptionally steep. Pluto's unusual phase curve may be due to an unusual surface texture caused by seasonal volatile transport and active geologic processes. These observations enable accurate determination of Pluto's geometric albedo in the JHK filter system, which we find to be 0.86 ± 0.04, 0.59 ± 0.05, and 0.39 ± 0.04 for Pluto, respectively, and 0.68 ± 0.06 for Charon in the J filter.

Pluto’s Haze Abundance and Size Distribution from Limb Scatter Observations by MVIC

Abstract The New Horizons spacecraft observed Pluto and Charon at solar-phase angles between 16° and 169°. In this work, we use the Multispectral Visible Imaging Camera (MVIC) observations to construct multiwavelength phase curves of Pluto’s atmosphere, using the limb scatter technique. Observational artifacts and biases were removed using Charon as a representative airless body. The size and distribution of the haze particles were constrained using a Titan fractal aggregate phase function. We find that monodispersed and log-normal populations cannot simultaneously describe the observed steep forward scattering, indicative of wavelength-scale particles, and the non-negligible backscattering indicative of particles much smaller than the wavelength. Instead, we find it necessary to use bimodal or power-law distributions, especially below ∼200 km, to properly describe the MVIC observations. Above 200 km, where the atmosphere is isotropically scattering, a monodisperse, log-normal, or a bimodal/power law approximating a monodispersed population is able to fit the phase curves well. As compared to the results of previously published articles, we find that Pluto’s atmosphere must contain haze particle number densities an order of magnitude greater for small (∼10 nm) and large (∼1 μm) radii, and relatively fewer intermediate sizes (∼100 nm). These conclusions support a lower aggregate aerosol growth rate than that found by Gao et al., indicating a higher charge-to-radius ratio, upwards of 60e − μm−1. In order to generate large particles with a lower growth rate, the atmosphere must also have a lower sedimentation velocity (<∼0.01 m s−1 at 200 km), which is possible with a fractal dimension of less than 2.

Optical-to-NIR magnitude measurements of the Starlink LEO Darksat satellite and effectiveness of the darkening treatment

Aims. We aim to measure the Sloan r′, Sloan i′, J, and Ks magnitudes of Starlink’s STARLINK-1130 (Darksat) and STARLINK-1113 low Earth orbit (LEO) communication satellites and determine the effectiveness of the Darksat darkening treatment from the optical to the near-infrared (NIR). Methods. Four observations of Starlink’s LEO communication satellites, Darksat and STARLINK-1113, were conducted on two nights with two telescopes. The Chakana 0.6 m telescope at the Ckoirama observatory (Chile) observed both satellites on 5 Mar. 2020 (UTC) and 7 Mar. 2020 (UTC) using a Sloan r′ and Sloan i′ filter, respectively. The ESO VISTA 4.1 m telescope with the VIRCAM instrument observed both satellites on 5 Mar. 2020 (UTC) and 7 Mar. 2020 (UTC) in the NIR J-band and Ks-band, respectively. Results. The calibration, image processing, and analysis of the Darksat images give r ≈ 5.6 mag, i ≈ 5.0 mag, J ≈ 4.2 mag, and Ks ≈ 4.0 mag when scaled to a range of 550 km (airmass = 1) and corrected for the solar incidence and observer phase angles. In comparison, the STARLINK-1113 images give r ≈ 4.9 mag, i ≈ 4.4 mag, J ≈ 3.8 mag, and Ks ≈ 3.6 mag when corrected for range, solar incidence, and observer phase angles. The data and results presented in this work show that the special darkening coating used by Starlink for Darksat has darkened the Sloan r′ magnitude by 50%, Sloan i′ magnitude by 42%, NIR J magnitude by 32%, and NIR Ks magnitude by 28%. Conclusions. The results show that both satellites increase in reflective brightness with increasing wavelength and that the effectiveness of the darkening treatment is reduced at longer wavelengths. This shows that the mitigation strategies being developed by Starlink and other LEO satellite operators need to take into account other wavelengths, not just the optical. This work highlights the continued importance of obtaining multi-wavelength observations of many different LEO satellites in order to characterise their reflective properties and to aid the community in developing impact simulations and developing mitigation tools.

Photometry of asteroid (101955) Bennu with OVIRS on OSIRIS-REx

NASA's OSIRIS-REx spacecraft arrived at its sampling target, asteroid (101955) Bennu, in December 2018 and started a series of global observation campaigns. Here we investigate the global photometric properties of Bennu as observed by the OSIRIS-REx Visible and InfraRed Spectrometer (OVIRS) over the time period December 9, 2018, to September 26, 2019. In this study we used observations obtained over wavelengths ranging from 0.4 to 3.7 μm, with a solar phase angle range of 5.3° to 132.6°. Our aim is to characterize the global average disk-resolved photometric properties of Bennu with multiple models. The best-fit model is a McEwen model with an exponential phase function and an exponential polynomial partition function. We use this model to correct the OVIRS spectra of Bennu to a standard reference viewing and illumination geometry at visible to infrared wavelengths for the purposes of global spectral mapping. We derive a bolometric Bond albedo map in which Bennu's surface values range from 0.021 to 0.027. We find a phase reddening effect, and our model is effective at removing this phase reddening. Our average model albedo shows a blueish spectrum with a > 10% absorption feature centered at 2.74 μm. Of all comparisons with previously visited asteroids and comets, only 28P/Neujmin, 2P/Encke, and (162173) Ryugu are darker than Bennu. We find that Bennu is a few percent brighter than Ryugu in the wavelengths respectively observed by the OSIRIS-REx and Hayabusa2 missions (from 0.48 to 0.86 μm). We also compare our spectroscopic photometry of Bennu with the OSIRIS-REx imaging photometry and with ground-based predictions.

A study of photometric observations of (3200) phaethon from Lulin Observatory

The near-Earth asteroid (3200) Phaethon is regarded as the parent body of the Geminid meteor stream. It is also an active asteroid with a very blue spectrum which is similar to that of the members of the Pallas family. For obtaining this asteroid's surface property in more detail, we conducted a series of multicolor photometry of this asteroid over a wide range of solar phase angle during its close approach to the Earth in the end of 2017. Our observations confirmed that Phaethon is a C-complex asteroid as previous studies have claimed. We also obtained the phase curve for this asteroid in solar phase angle of 22.67 ∼ 114.48°. From this result, we estimate that the absolute magnitude (H) of this asteroid is 14.450 ± 0.017 mag, and its phase slope is (β) ∼0.0422 ± 0.0005. Using the absolute magnitude that we estimated, together with an estimate of Phaethon's albedo (0.122 ± 0.03) that previous literature derived, we can derived the Phaethon's diameter, 4.90−0.54+0.79 km. The diameter estimate is roughly consistent with what was observed through the radar ranging measurement conducted as the same period as our observation.

Examining the effects of on-orbit aging of SL-12 rocket bodies using visible band spectra with the MMT telescope

The characterization of deep space debris has posed a significant challenge in the study of space objects. To be most effective, characterization must be performed quickly and under non-ideal operational conditions, generally using non-resolved techniques. Multi-color photometry and the resultant color indices offer the potential to rapidly discriminate between debris and intact space objects such as rocket bodies and satellites. However, these studies are not well informed by high resolution spectra of these same objects, due to the lack of prior measurements with large astronomical telescopes. We present a spectroscopic study of orbital debris aimed at characterizing changes over time in object properties, using ground-based spectroscopy of several similar rocket bodies in geosynchronous orbit, launched at different times. High and moderate resolution spectroscopy is not routinely collected by SSA resources. Nonetheless, several researchers have collected satellite spectra for research purposes. Several researchers have also noted the progressive reddening of spacecraft surfaces.In this study, we have collected high resolution spectra of five Russian SL-12 rocket bodies in geosynchronous orbit. The spectra were collected with the Blue Channel Spectrograph on the 6.5 m MMT telescope at Mt. Hopkins. The measurements were taken using the 300-line grating, which is blazed for the red, and can cover a 520 nm range at dispersion 0.196 nm/pixel. The large collecting aperture of the MMT allowed the rapid collection of multiple high signal-to-noise spectra with only 2 minutes per exposure. This short exposure allowed us to have confidence the solar phase angle was near constant during each collection, but that the spectra were averaged over the rotation of the rocket body. These spectra allow analysis of both the variation in albedo over a large wavelength range, and searches for discrete absorption features. The SL-12 (also called the "Proton K") was a mainstay Russian four-stage to GEO launch vehicle that was used from 1974 to 2012. The SL-12 fourth stage rocket bodies (henceforth referred to as "SL-12 RB") offer a convenient ensemble of objects for which photometric techniques can be developed and tested. For this study, spectra of five SL-12s with a range of years-on-orbit (YOO) ranging from 23-35 years were collected, allowing a comparative study of the evolution of the spectra over a 12-year difference in age. Additionally, all these objects have been previously observed in the near-IR with the UKIRT WFCAM. The spectra are analyzed for evidence of the effects of on-orbit reddening and other changes over time.

First observations and magnitude measurement of Starlink’s Darksat

Aims. We measured the Sloan g′ magnitudes of the Starlink’s STARLINK-1130 (Darksat) and 1113 low Earth orbit (LEO) communication satellites to determine the effectiveness of the Darksat darkening treatment at 475.4 nm. Methods. Two observations of the Starlink’s Darksat LEO communication satellite were conducted on 2020/02/08 and 2020/03/06 using Sloan r′ and g′ filters, respectively. A second satellite, STARLINK-1113, was observed on 2020/03/06 using a Sloan g′ filter. The initial observation on 2020/02/08 was a test observation conducted when Darksat was still in the process of manoeuvring to its nominal orbit and orientation. Based on the successful test observation, the first main observation took place on 2020/03/06, along with an observation of the second Starlink satellite. Results. The calibration, image processing, and analysis of the Darksat Sloan g′ image gives an estimated Sloan g′ magnitude of 7.46 ± 0.04 at a range of 976.50 km. For STARLINK-1113, an estimated Sloan g′ magnitude of 6.59 ± 0.05 at a range of 941.62 km was found. When scaled to a range of 550 km and corrected for the solar and observer phase angles, a reduction by a factor of two is seen in the reflected solar flux between Darksat and STARLINK-1113. Conclusions. The data and results presented in this work demonstrate that the special darkening coating used by Starlink for Darksat has darkened the Sloan g’ magnitude by 0.77 ± 0.05 mag when the range is equal to a nominal orbital height (550 km). This result will serve members of the astronomical community who are actively modelling the satellite mega-constellations to ascertain their actual impact on both amateur and professional astronomical observations. Both concurrent and subsequent observations are planned to cover the full optical and NIR spectrum using an ensemble of instruments, telescopes, and observatories.

Enhancing Debris Tracking Missions in Geostationary Orbit with Advanced Solar Photon Thrusters

The advanced solar photon thruster (ASPT) concept was first proposed, developed, and analyzed by Dachwald and Wurm. Since its debut as an alternative to the flat solar sail (FSS) for interplanetary missions, it has been deemed inferior for its technical complexity, poor scalability, and time-optimal transfer trajectory performance. A novel application of the ASPT to the space object tracking and detection mission has revealed several key performance benefits illuminating small debris in geostationary orbit (GEO). A numerical procedure to optimize the ASPT average solar reflected intensity was developed, and the optimized performance results were compared with those of an FSS with an equal total mass. Results demonstrate that although the FSS offers slightly better acceleration performance, in agreement with previously published literature, the ASPT can provide larger reflected solar intensities over a wider range of applicable phase angles, resulting in a 200% detection improvement when tracking a 10-cm-diam object. In addition, the vehicle slewing requirements for the ASPT are reduced to just one revolution per day as the surface directing the sun light onto small objects is decoupled from the large solar collector pointing at the sun. More specifically, a ground observer could continually detect small debris when illuminated by the ASPT, whereas the FSS design allowed for its detection when solar phase angles were 60 deg or larger. The ASPT offers substantial performance benefits to the GEO small debris tracking mission when considering flight safety, anomaly characterization, local area clearing around high-value assets, and space object catalog maintenance activities.

Dayside cloud top structure of Venus retrieved from Akatsuki IR2 observations

We describe the dayside cloud top structure of Venus as retrieved from 93 images acquired at a wide variety of solar phase angles (0–120°) using the 2.02-μm channel of the 2-μm camera (IR2) onboard the Venus orbiter, Akatsuki, from April 4 to May 25, 2016. Since the 2.02-μm channel is located in a CO2 absorption band, the sunlight reflected from Venus allowed us to determine the cloud top altitude corresponding to unit aerosol optical depth at 2.02 μm. First, the observed solar phase angle dependence and the center-to-limb variation of the reflected sunlight in the region equatorward of 30° were used to construct a spatially averaged cloud top structure characterized by cloud top altitude zc, Mode 2 modal radius rg,2, and cloud scale height H, which were 70.4 km, 1.06 μm, and 5.3 km, respectively. Second, cloud top altitudes at individual locations were retrieved on a pixel-by-pixel basis with an assumption that rg,2 and H were uniform for the entire planet. The latitudinal structure of the cloud top altitude was symmetric with respect to the equator. The average cloud top altitude was 70.5 km in the equatorial region and showed a gradual decrease of ~2 km by the 45° latitude. It rapidly dropped at latitudes of 50–60° and reached 61 km in latitudes of 70–75°. The average cloud top altitude in the region equatorward of 30° showed negligible local time dependence, with changes up to 1 km at most. Local variations in cloud top altitude, including stationary gravity wave features, occurred within several hundreds of meters. Although long zonal or tilted streaky features poleward of ~45° were clearly identifiable, features in the low and middle latitudes were usually subtle. These did not necessarily appear as local variations at the cloud top level, where mottled and patchy UV patterns were observed, suggestive of convection and turbulence at the cloud top level.

Multi wavelength optical broad band photometric properties of a representative sample of geostationary satellites

We present the results of multi-wavelength optical observations, using B, V, R and I band filters, of a sample of 61 active geostationary communications satellites, on 2015 July 9–11. These observations were done outside the glinting season, with a longitudinal solar angle smaller than 6° and a latitudinal solar phase angle of ~16.5°. This approach enabled us to determine the distribution of quiescent magnitudes and colors of these targets during a period of the night when they are expected to be the brightest. We present the distribution of magnitudes and colors of the sample. We find that the colors of these satellites are in general redder than the Solar spectrum, with a wide distribution of values (0.62<B-V<1.49,1.00<B-R<2.28,1.31<B-I<3.07,0.34<V-R<0.90,0.63<V-I<1.69,0.27<R-I<0.79). Satellites from different manufacturers are found to cluster around specific colors in color-color plots, indicating that this information can be used to identify different satellites and resolve cross tagging issues. We also investigate whether the colors of satellites change with time and find that in most cases there is no significant change to these properties as the satellites age.