We report the recovery of returning Halley-type comet 12P/Pons-Brooks using the 4.3 m Lowell Discovery Telescope, at a heliocentric distance of 11.89 au. Comparative analysis with a dust model suggests that the comet may have been active since $\sim30$ au from the Sun. We derive a nucleus radius of $17\pm6$ km from the nucleus photometry, though this number is likely an overestimation due to the contamination from dust and gas. Continuing monitoring is encouraged in anticipation of the comet's forthcoming perihelion in 2024 April.

We present a study of interstellar comet 2I/2019 Q4 (Borisov) using both preperihelion and postperihelion observations spanning late September 2019 through late January 2020. The intrinsic brightness of the comet was observed to continuously decline throughout the timespan, likely due to the decreasing effective scattering cross-section as a result of volatile sublimation with a slope of $-0.43 \pm 0.02$ km$^{2}$ d$^{-1}$. We witnessed no significant change in the slightly reddish colour of the comet, with mean values of $\left \langle g - r \right \rangle = 0.68 \pm 0.04$, $\left \langle r - i \right \rangle = 0.23 \pm 0.03$, and the normalised reflectivity gradient across the $g$ and $i$ bands $\overline{S'} \left(g,i\right) = \left(10.6 \pm 1.4\right)$ % per $10^3$ \AA, all unremarkable in the context of solar system comets. Using the available astrometric observations, we confidently detect the nongravitational acceleration of the comet following a shallow heliocentric distance dependency of $r_{\rm H}^{-1 \pm 1}$. Accordingly, we estimate that the nucleus is most likely $\lesssim$0.4 km in radius, and that a fraction of $\gtrsim$0.2% of the total nucleus in mass has been eroded due to the sublimation activity since the earliest observation of the comet in 2018 December by the time of perihelion. Our morphology simulation suggests that the dust ejection speed increased from $\sim$4 m s$^{-1}$ in September 2019 to $\sim$7 m s$^{-1}$ around perihelion for the optically dominant dust grains of $\beta \sim 0.01$, and that the observable dust grains are no smaller than micron size.

We present a study of C/2019 Y4 (ATLAS) using Sloan $gri$ observations from mid-January to early April 2020. During this timespan, the comet brightened with a growth in the effective cross-section of $\left(2.0 \pm 0.1 \right) \times 10^{2}$ m$^{2}$ s$^{-1}$ from the beginning to $\sim$70 d preperihelion in late March 2020, followed by a brightness fade and the comet gradually losing the central condensation. Meanwhile, the comet became progressively bluer, and was even bluer than the Sun ($g - r \approx 0.2$) when the brightness peaked, likely due to activation of subterranean fresh volatiles exposed to sunlight. With the tailward-bias corrected astrometry we found an enormous radial nongravitational parameter, $A_{1} = \left(+2.25 \pm 0.13\right) \times 10^{-7}$ au d$^{-2}$ in the heliocentric motion of the comet. Taking all of these finds into consideration, we conclude that the comet has disintegrated since mid-March 2020. By no means was the split new to the comet, as we quantified that the comet had undergone another split event around last perihelion $\sim$5 kyr ago, during which its sibling C/1844 Y1 (Great Comet) was produced, with the in-plane component of the separation velocity $\gtrsim$1 m s$^{-1}$. We constrained that the nucleus of C/2019 Y4 before disintegration was $\gtrsim$60 m in radius, and has been protractedly ejecting dust grains of $\sim$10-40 $\mu$m (assuming dust bulk density 0.5 g cm$^{-3}$) with ejection speed $\sim$30 m s$^{-1}$ in early March 2020 and increased to $\sim$80 m s$^{-1}$ towards the end of the month for grains of $\sim$10 $\mu$m.

We present a preliminary analysis of comet P/2019 LM$_4$ (Palomar) as observed by the Zwicky Transient Facility (ZTF) survey in 2019 and 2020. We find that the discovery of the comet in 2019 and the recovery in 2020 is largely attributed to two separate outbursts that are $\gtrsim2$ and $\gtrsim3.9$ mag in strength. The outbursts occurred around the end of April to early May of 2019 as well as between 2020 May 8.31 and 9.52 UTC, respectively.

Granvik et al. (2016) report an absence of asteroids on orbits with perihelia near the Sun that they attribute to the 'supercatastrophic disruption' of these bodies. Here we investigate whether there is evidence for this process among other bodies with similarly low perihelia: near-Earth asteroids, SOHO comets, as well as meter-sized and millimeter-sized meteoroids. We determine no known near-Earth asteroids have past (last 10^4 years) histories residing significantly inside the Granvik et al. (2016) limit, indirectly supporting the disruption hypothesis. The exception is asteroid (467372) 2004 LG which spent 2500 years within this limit, and thus presents a challenge to that theory. Phaethon has a perihelion distance hovering just above the limit and may be undergoing slow disruption, which may be the source of its dust complex. We find that the rate at which ungrouped SOHO comets are observed is consistent with expected rates for the injection of small (25 m) class asteroids into the near-Sun region and suggest that this fraction of the SOHO-observed comet population may in fact be asteroidal in origin. We also find that there is an absence of meter-sized bodies with near-Sun perihelia but an excess of millimeter-sized meteoroids. This implies that if near-Sun asteroids disrupt, they do not simply fragment into meter-sized chunks but disintegrate ultimately into millimeter-sized particles. We propose that the disruption of near-Sun asteroids as well as the anomalous brightening and destruction processes that affect SOHO comets occur through meteoroid erosion, that is, the removal of material through impacts by high-speed near-Sun meteoroids.

The discovery of the second interstellar object 2I/Borisov on 2019 August 30 raises the question of whether it was ejected recently from a nearby stellar system. Here we compute the asymptotic incoming trajectory of 2I/Borisov, based on both recent and pre-discovery data extending back to December 2018, using a range of force models that account for cometary outgassing. From Gaia DR2 astrometry and radial velocities, we trace back in time the Galactic orbits of 7.4 million stars to look for close encounters with 2I/Borisov. The closest encounter we find took place 910kyr ago with the M0V star Ross 573, at a separation of 0.068pc (90% confidence interval of 0.053-0.09pc) with a relative velocity of 23km/s. This encounter is nine times closer than the closest past encounter identified for the first interstellar object 1I/'Oumuamua. Ejection of 2I/Borisov via a three-body encounter in a binary or planetary system is possible, although such a large ejection velocity is unlikely to be obtained and Ross 573 shows no signs of binarity. We also identify and discuss some other recent close encounters, recognizing that if 2I/Borisov is more than about 10Myr old, our search would be unlikely to find its parent system.

Comet 2I/Borisov, the first unambiguous interstellar comet ever found, was discovered in August 2019 at $\sim3$ au from the Sun on its inbound leg. No pre-discovery detection beyond 3 au has yet been reported, mostly due to the comet's proximity to the Sun as seen from the Earth. Here we present a search for pre-discovery detections of comet Borisov using images taken by the Catalina Sky Survey (CSS), Pan-STARRS and Zwicky Transient Facility (ZTF), with a further comprehensive follow-up campaign being presented in \citet{Bolin2019}. We identified comet Borisov in ZTF images taken in May 2019 and use these data to update its orbit. This allowed us to identify the comet in images acquired as far back as December 2018, when it was 7.8 au from the Sun. The comet was not detected in November 2018 when it was 8.6 au from the Sun, possibly implying an onset of activity around this time. This suggests that the activity of the comet is either driven by a more volatile species other than H$_2$O, such as CO or CO$_2$, or by exothermic crystallization of amorphous ice. We derive the radius of the nucleus to be $<7$ km using the non-detection in November 2018, and estimate an area of $\sim0.5$---$10 \mathrm{km^2}$ has been active between December 2018 and September 2019, though this number is model-dependent and is highly uncertain. The behavior of comet Borisov during its inbound leg is observationally consistent with dynamically new comets observed in our solar system, suggesting some similarities between the two.

Near-Earth Objects (NEOs) that orbit the Sun on or within Earth's orbit are tricky to detect for Earth-based observers due to their proximity to the Sun in the sky. These small bodies hold clues to the dynamical history of the inner solar system as well as the physical evolution of planetesimals in extreme environments. Populations in this region include the Atira and Vatira asteroids, as well as Venus and Earth co-orbital asteroids. Here we present a twilight search for these small bodies, conducted using the 1.2-m Oschin Schmidt and the Zwicky Transient Facility (ZTF) camera at Palomar Observatory. The ZTF twilight survey operates at solar elongations down to $35^\circ$ with limiting magnitude of $r=19.5$. During a total of 40 evening sessions and 62 morning sessions conducted between 2018 November 15 and 2019 June 23, we detected 6 Atiras, including 2 new discoveries 2019 AQ$_3$ and 2019 LF$_6$, but no Vatiras or Earth/Venus co-orbital asteroids. NEO population models show that these new discoveries are likely only the tip of the iceberg, with the bulk of the population yet to be found. The population models also suggest that we have only detected 5--$7\%$ of the $H<20$ Atira population over the 7-month survey. Co-orbital asteroids are smaller in diameters and require deeper surveys. A systematic and efficient survey of the near-Sun region will require deeper searches and/or facilities that can operate at small solar elongations.

Comets are primitive objects that formed in the protoplanetary disk, and have been largely preserved over the history of the Solar System. However, they are not pristine, and surfaces of cometary nuclei do evolve. In order to understand the extent of their primitive nature, we must define the mechanisms that affect their surfaces and comae. We examine the lightcurve of comet 240P/NEAT over three consecutive orbits, and investigate three events of significant brightening ($\Delta m \sim -2$ mag). Unlike typical cometary outbursts, each of the three events are long-lived, with enhanced activity for at least 3 to 6 months. The third event, observed by the Zwicky Transient Facility, occurred in at least two stages. The anomalous behavior appears to have started after the comet was perturbed by Jupiter in 2007, reducing its perihelion distance from 2.53 to 2.12 au. We suggest that the brightening events are temporary transitions to a higher baseline activity level, brought on by the increased insolation, which has warmed previously insulated sub-surface layers. The new activity is isolated to one or two locations on the nucleus, indicating that the surface or immediate sub-surface is heterogeneous. Further study of this phenomenon may provide insight into cometary outbursts, the structure of the near-surface nucleus, and cometary nucleus mantling.

We present DeepStreaks, a convolutional-neural-network, deep-learning system designed to efficiently identify streaking fast-moving near-Earth objects that are detected in the data of the Zwicky Transient Facility (ZTF), a wide-field, time-domain survey using a dedicated 47 sq. deg camera attached to the Samuel Oschin 48-inch Telescope at the Palomar Observatory in California, United States. The system demonstrates a 96-98% true positive rate, depending on the night, while keeping the false positive rate below 1%. The sensitivity of DeepStreaks is quantified by the performance on the test data sets as well as using known near-Earth objects observed by ZTF. The system is deployed and adapted for usage within the ZTF Solar-System framework and has significantly reduced human involvement in the streak identification process, from several hours to typically under 10 minutes per day.

The Zwicky Transit Factory (ZTF) is a powerful time domain survey facility with a large field of view. We apply the synthetic tracking technique to integrate a ZTF's long-dwell dataset, which consists of 133 nominal 30-second exposure frames spanning about 1.5 hours, to search for slowly moving asteroids down to approximately 23rd magnitude. We found more than one thousand objects from searching 40 CCD-quadrant subfields, each of which covers a field size of $\sim$0.73 deg$^2$. While most of the objects are main belt asteroids, there are asteroids belonging to families of Trojan, Hilda, Hungaria, Phocaea, and near-Earth-asteroids. Such an approach is effective and productive. Here we report the data process and results.

Jupiter-family comet 289P/Blanpain was first discovered in 1819 and was then lost for $\sim200$ years, only to be rediscovered in 2003 as a small, weakly active comet. The comet is associated with the Phoenicids, an otherwise minor meteor shower that produced significant outbursts in 1956 and 2014. The shower points to the existence of significant mass-loss events of P/Blanpain in recent history. P/Blanpain was recovered during an apparent large outburst in 2013 July at an appreciable heliocentric distance of 3.9 au, with brightness increase of 9 mag, making it one of the largest comet outbursts ever observed. Here we present an analysis of archival data taken by several telescopes. We find that the 2013 outburst has produced $\sim10^{8}$ kg of dust, which accounts for a modest fraction ($\sim1\%$) of the mass of P/Blanpain's nucleus as measured in 2004. Based on analysis of long-term lightcurve and modeling of coma morphology, we conclude that the 2013 outburst was most likely driven by the crystallization of amorphous water ice triggered by a spin-up disruption of the nucleus. Dust dynamical model shows that a small fraction of the dust ejecta will reach the Earth in 2036 and 2041, but are only expected to produce minor enhancements to the Phoenicid meteor shower. The 2013 outburst of P/Blanpain, though remarkable for a comet at small sizes, does not necessary imply a catastrophic disruption of the nucleus. The upcoming close encounter of P/Blanpain in 2020 January will provide an opportunity to examine the current state of the comet.

We report on a multi-observatory campaign to examine asteroid 3200 Phaethon during its December 2017 close approach to Earth, in order to improve our measurements of its fundamental parameters, and to search for surface variations, cometary activity and fragmentation. The mean colors of Phaethon are B-V = 0.702 +/- 0.004, V-R = 0.309 +/- 0.003, R-I = 0.266 +/- 0.004, neutral to slightly blue, consistent with previous classifications of Phaethon as a F-type or B-type asteroid. Variations in Phaethon's B-V colors (but not V-R or R-I) with observer sub-latitude are seen and may be associated with craters observed by the Arecibo radar. High cadence photometry over phases from 20 to 100 degrees allows a fit to the values of the HG photometric parameters; H = 14.57 +/- 0.02, 13.63 +/- 0.02, 13.28 +/- 0.02, 13.07 +/- 0.02; G = 0.00 +/- 0.01, -0.09 +/- 0.01, -0.10 +/- 0.01, -0.08 +/- 0.01 in the BVRI filters respectively; the negative G values are consistent with other observations of F type asteroids. Light curve variations were seen that are also consistent with concavities reported by Arecibo, indicative of large craters on Phaethon's surface whose ejecta may be the source of the Geminid meteoroid stream. A search for gas/dust production set an upper limit of 0.06 +/- 0.02 kg/s when Phaethon was 1.449 AU from the Sun, and 0.2 +/- 0.1 kg/s at 1.067 AU. A search for meter-class fragments accompanying Phaethon did not find any whose on-sky motion was not also consistent with background main belt asteroids.

We present a novel algorithm for scheduling the observations of time-domain imaging surveys. Our Integer Linear Programming approach optimizes an observing plan for an entire night by assigning targets to temporal blocks, enabling strict control of the number of exposures obtained per field and minimizing filter changes. A subsequent optimization step minimizes slew times between each observation. Our optimization metric self-consistently weights contributions from time-varying airmass, seeing, and sky brightness to maximize the transient discovery rate. We describe the implementation of this algorithm on the surveys of the Zwicky Transient Facility and present its on-sky performance.

We describe ZStreak, a semi-real-time pipeline specialized in detecting small, fast-moving near-Earth asteroids (NEAs) that is currently operating on the data from the newly-commissioned Zwicky Transient Facility (ZTF) survey. Based on a prototype originally developed by Waszczak et al. (2017) for the Palomar Transient Factory (PTF), the predecessor of ZTF, ZStreak features an improved machine-learning model that can cope with the $10\times$ data rate increment between PTF and ZTF. Since its first discovery on 2018 February 5 (2018 CL), ZTF/ZStreak has discovered $45$ confirmed new NEAs over a total of 232 observable nights until 2018 December 31. Most of the discoveries are small NEAs, with diameters less than $\sim100$ m. By analyzing the discovery circumstances, we find that objects having the first to last detection time interval under 2 hr are at risk of being lost. We will further improve real-time follow-up capabilities, and work on suppressing false positives using deep learning.

Main-belt asteroid (6478) Gault unexpectedly sprouted two tails in late 2018 and early 2019, identifying it as a new active asteroid. Here we present observations obtained by the 1.2-m Zwicky Transient Facility survey telescope that provide detailed time-series coverage of the onset and evolution of Gault's activity. Gault exhibited two brightening events, with the first one starting on 2018 Oct. 18$\pm5$ days and a second one starting on 2018 Dec. 24$\pm1$ days. The amounts of mass released are $2\times10^7$ kg and $1\times10^6$ kg, respectively. Based on photometric measurements, each event persisted for about a month. Gault's color has not changed appreciably over time, with a pre-outburst color of $g_\mathrm{PS1}-r_\mathrm{PS1}=0.50\pm0.04$ and $g_\mathrm{PS1}-r_\mathrm{PS1}=0.46\pm0.04$ during the two outbursts. Simulations of dust dynamics shows that the ejecta consists of dust grains of up to 10 $\mu$m in size that are ejected at low velocities below $1~\mathrm{m~s^{-1}}$ regardless of particle sizes. This is consistent with non-sublimation-driven ejection events. The size distribution of the dust exhibits a broken power-law, with particles at 10--20 $\mu$m following a power-law of $-2.5$ to $-3.0$, while larger particles follow a steeper slope of $-4.0$. The derived properties can be explained by either rotational excitation of the nucleus or a merger of a near-contact binary, with the latter scenario to be statistically more likely.

The under-abundance of asteroids on orbits with small perihelion distances suggests that thermally-driven disruption may be an important process in the removal of rocky bodies in the Solar System. Here we report our study of how the debris streams arise from possible thermally-driven disruptions in the near-Sun region. We calculate that a small body with a diameter $\gtrsim0.5$ km can produce a sufficient amount of material to allow the detection of the debris at the Earth as meteor showers, and that bodies at such sizes thermally disrupt every $\sim2$ kyrs. We also find that objects from the inner parts of the asteroid belt are more likely to become Sun-approacher than those from the outer parts. We simulate the formation and evolution of the debris streams produced from a set of synthetic disrupting asteroids drawn from Granvik et al. (2016)'s near-Earth object population model, and find that they evolve 10--70 times faster than streams produced at ordinary solar distances. We compare the simulation results to a catalog of known meteor showers on Sun-approaching orbits. We show that there is a clear overabundance of Sun-approaching meteor showers, which is best explained by a combining effect of comet contamination and an extended disintegration phase that lasts up to a few kyrs. We suggest that a few asteroid-like Sun-approaching objects that brighten significantly at their perihelion passages could, in fact, be disrupting asteroids. An extended period of thermal disruption may also explain the widespread detection of transiting debris in exoplanetary systems. Data and codes that generate the figures and main results of this work are publicly available on 10.5281/zenodo.2547298 and https://github.com/Yeqzids/near-sun-disruptions.

The Zwicky Transient Facility is a large optical survey in multiple filters producing hundreds of thousands of transient alerts per night. We describe here various machine learning (ML) implementations and plans to make the maximal use of the large data set by taking advantage of the temporal nature of the data, and further combining it with other data sets. We start with the initial steps of separating bogus candidates from real ones, separating stars and galaxies, and go on to the classification of real objects into various classes. Besides the usual methods (e.g., based on features extracted from light curves) we also describe early plans for alternate methods including the use of domain adaptation, and deep learning. In a similar fashion we describe efforts to detect fast moving asteroids. We also describe the use of the Zooniverse platform for helping with classifications through the creation of training samples, and active learning. Finally we mention the synergistic aspects of ZTF and LSST from the ML perspective.

The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48-inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg$^2$ field of view and 8 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF's public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope.

The Zwicky Transient Facility (ZTF) is a new robotic time-domain survey currently in progress using the Palomar 48-inch Schmidt Telescope. ZTF uses a 47 square degree field with a 600 megapixel camera to scan the entire northern visible sky at rates of ~3760 square degrees/hour to median depths of g ~ 20.8 and r ~ 20.6 mag (AB, 5sigma in 30 sec). We describe the Science Data System that is housed at IPAC, Caltech. This comprises the data-processing pipelines, alert production system, data archive, and user interfaces for accessing and analyzing the products. The realtime pipeline employs a novel image-differencing algorithm, optimized for the detection of point source transient events. These events are vetted for reliability using a machine-learned classifier and combined with contextual information to generate data-rich alert packets. The packets become available for distribution typically within 13 minutes (95th percentile) of observation. Detected events are also linked to generate candidate moving-object tracks using a novel algorithm. Objects that move fast enough to streak in the individual exposures are also extracted and vetted. The reconstructed astrometric accuracy per science image with respect to Gaia is typically 45 to 85 milliarcsec. This is the RMS per axis on the sky for sources extracted with photometric S/N >= 10. The derived photometric precision (repeatability) at bright unsaturated fluxes varies between 8 and 25 millimag. Photometric calibration accuracy with respect to Pan-STARRS1 is generally better than 2%. The products support a broad range of scientific applications: fast and young supernovae, rare flux transients, variable stars, eclipsing binaries, variability from active galactic nuclei, counterparts to gravitational wave sources, a more complete census of Type Ia supernovae, and Solar System objects.

The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single 1-m class survey telescope. The public surveys will cover the observable northern sky every three nights in g and r filters and the visible Galactic plane every night in g and r. Alerts generated by these surveys are sent in real time to brokers. A consortium of universities which provided funding ("partnership") are undertaking several boutique surveys. The combination of these surveys producing one million alerts per night allows for exploration of transient and variable astrophysical phenomena brighter than r $\sim$ 20.5 on timescales of minutes to years. We describe the primary science objectives driving ZTF including the physics of supernovae and relativistic explosions, multi-messenger astrophysics, supernova cosmology, active galactic nuclei and tidal disruption events, stellar variability, and Solar System objects.

We propose a LSST Solar System near-Sun Survey, to be implemented during twilight hours, that extends the seasonal reach of LSST to its maximum as fresh sky is uncovered at about 50 square degrees per night (1500 sq. deg. per lunation) in the morning eastern sky, and surveyable sky is lost at the same rate to the western evening sky due to the Earth's synodic motion. By establishing near-horizon fence post picket lines to the far west and far east we address Solar System science use cases (including Near Earth Objects, Interior Earth Objects, Potentially Hazardous Asteroids, Earth Trojans, near-Sun asteroids, sun-grazing comets, and dormant comets) as well as provide the first look and last look that LSST will have at the transient and variable objects within each survey field. This proposed near-Sun Survey will also maximize the overlap with the field of regard of the proposed NEOCam spacecraft that will be stationed at the Earth's L1 Lagrange point and survey near quadrature with the Sun. This will allow LSST to incidently follow-up NEOCam targets and vice-versa (as well as targets from missions such as Euclid), and will roughly correspond to the Earth's L4 and L5 regions.

Near-Earth asteroid (3200) Phaethon is notable for its association to a strong annual meteor shower, the Geminids, indicative of one or more episodes of mass ejection in the past. The mechanism of Phaethon's past activity is not yet understood. Here we present a Hubble Space Telescope (HST) search of meter-sized fragments in the vicinity of Phaethon, carried out during Phaethon's historic approach to the Earth in mid-December of 2017. Numerical simulation conducted to guide HST's pointing also show that the dynamical evolution of Phaethon-originated particles is quick, as ejected materials take no longer than $\sim250$ yr to spread to the entire orbit of Phaethon. Our search was completed down to 4-meter-class limit (assuming Phaethon-like albedo) and was expected to detect 0.035% particles ejected by Phaethon in the last several decades. The negative result of our search capped the total mass loss of Phaethon over the past few dozen orbits to be $10^{12}$ kg at $3\sigma$ level, taking the best estimates of size power-law from meteor observations and spacecraft data. Our result also implies a millimeter-sized dust flux of $<10^{-12} \mathrm{m^{-2} s^{-1}}$ within 0.1 au of Phaethon, suggesting that any Phaethon-bound mission is unlikely to encounter dense dust clouds.

Small bodies in the solar system are conventionally classified into asteroids and comets. However, it is recently found that a small number of objects can exhibit properties of both asteroids and comets. Some are more consistent with asteroids despite episodic ejections and are labeled as "active asteroids", while some might be aging comets with depleting volatiles. Ejecta produced by active asteroids and/or dormant comets are potentially detectable as meteor showers at the Earth if they are in Earth-crossing orbits, allowing us to retrieve information about the historic activities of these objects. Meteor showers from small bodies with low and/or intermittent activities are usually weak, making shower confirmation and parent association challenging. We show that statistical tests are useful for identifying likely parent-shower pairs. Comprehensive analyses of physical and dynamical properties of meteor showers can lead to deepen understanding on the history of their parents. Meteor outbursts can trace to recent episodic ejections from the parents, and "orphan" showers may point to historic disintegration events. The flourish of NEO and meteor surveys during the past decade has produced a number of high-confidence parent-shower associations, most have not been studied in detail. More work is needed to understand the formation and evolution of these parent-shower pairs.

At 12:07 UT (8:07 pm China Standard Time) on 2017 Oct 4, a bright bolide was widely observed in the Shangri-La region in the Province of Yunnan, China (Figure 1). The event was well observed by the general public as it took place on the night of the Mid Autumn Festival which associates with moon gazing. Sonic booms and ground shaking were reported in an area about a thousand square kilometers wide to the northwest of the Shangri-La City. Data from the U.S. government sensor suggested that the impact energy of the event is approximately 0.54 kt TNT equivalent, with the terminus of the bolide positioned at $28.1^\circ$ N, $99.4^\circ$ E. This is the largest observed bolide event overland since the bolide event took place in Mauritania on 2016 Jun 27 (1.2 kt). In this Research Note we present a preliminary analysis of this event.

Jupiter-family Comet D/1770 L1 (Lexell) was the first discovered Near-Earth Object (NEO), and passed the Earth on 1770 Jul 1 at a recorded distance of 0.015 au. The comet was subsequently lost due to unfavorable observing circumstances during its next apparition followed by a close encounter with Jupiter in 1779. Since then, the fate of D/Lexell has attracted interest from the scientific community, and now we revisit this long-standing question. We investigate the dynamical evolution of D/Lexell based on a set of orbits recalculated using the observations made by Charles Messier, the comet's discoverer, and find that there is a $98\%$ chance that D/Lexell remains in the Solar System by the year of 2000. This finding remains valid even if a moderate non-gravitational effect is imposed. Messier's observations also suggest that the comet is one of the largest known near-Earth comets, with a nucleus of $\gtrsim 10$ km in diameter. This implies that the comet should have been detected by contemporary NEO surveys regardless of its activity level if it has remained in the inner Solar System. We identify asteroid 2010 JL$_{33}$ as a possible descendant of D/Lexell, with a $0.8\%$ probability of chance alignment, but a direct orbital linkage of the two bodies has not been successfully accomplished. We also use the recalculated orbit to investigate the meteors potentially originating from D/Lexell. While no associated meteors have been unambiguously detected, we show that meteor observations can be used to better constrain the orbit of D/Lexell despite the comet being long lost.

The Large Synoptic Survey Telescope (LSST) is uniquely equipped to search for Solar System bodies due to its unprecedented combination of depth and wide field coverage. Over a ten-year period starting in 2022, LSST will generate the largest catalog of Solar System objects to date. The main goal of the LSST Solar System Science Collaboration (SSSC) is to facilitate the efforts of the planetary community to study the planets and small body populations residing within our Solar System using LSST data. To prepare for future survey cadence decisions and ensure that interesting and novel Solar System science is achievable with LSST, the SSSC has identified and prioritized key Solar System research areas for investigation with LSST in this roadmap. The ranked science priorities highlighted in this living document will inform LSST survey cadence decisions and aid in identifying software tools and pipelines needed to be developed by the planetary community as added value products and resources before the planned start of LSST science operations.

1I/2017 U1 (`Oumuamua), a recently discovered asteroid in a hyperbolic orbit, is likely the first macroscopic object of extrasolar origin identified in the solar system. Here, we present imaging and spectroscopic observations of \textquoteleft Oumuamua using the Palomar Hale Telescope as well as a search of meteor activity potentially linked to this object using the Canadian Meteor Orbit Radar. We find that \textquoteleft Oumuamua exhibits a moderate spectral gradient of $10\%\pm6\%~(100~\mathrm{nm})^{-1}$, a value significantly lower than that of outer solar system bodies, indicative of a formation and/or previous residence in a warmer environment. Imaging observation and spectral line analysis show no evidence that \textquoteleft Oumuamua is presently active. Negative meteor observation is as expected, since ejection driven by sublimation of commonly-known cometary species such as CO requires an extreme ejection speed of $\sim40$ m s$^{-1}$ at $\sim100$ au in order to reach the Earth. No obvious candidate stars are proposed as the point of origin for \textquoteleft Oumuamua. Given a mean free path of $\sim10^9$ ly in the solar neighborhood, \textquoteleft Oumuamua has likely spent a very long time in the interstellar space before encountering the solar system.

Citizen science refers to scientific research conducted or participated by non-professional scientists (such as hobbyists or members from the general public). Citizen astronomy is a classic example of citizen science. Citizen astronomers has benefited from technological advancements in the recent decades as they fill the scientific gaps left by professional astronomers, in the areas such as time domain observations, visual classification and data mining. Chinese citizen astronomers have made a visible contribution in the discoveries of new objects; however, comparing to their counterparts in the western world, they appear to be less interested in researches that do not involve making new discovery, such as visual classification, long-term monitoring of objects, and data mining. From a questionnaire survey that aimed to investigate the motivation of Chinese citizen astronomers, we find that this population is predominantly male ($92\%$) who mostly reside in economically developed provinces. A large fraction ($69\%$) of the respondents are students and young professionals younger than the age of 25, which differs significantly from the occupation and age distribution of typical Chinese Internet users as well as the user distribution of large international citizen science projects such as the Galaxy Zoo. This suggests that youth generation in China is more willing to participate citizen astronomy research than average generation. Additionally, we find that interests in astronomy, desire to learn new knowledges, have a fun experience and meet new friends in the community are all important driving factors for Chinese citizen astronomers to participate research. This also differs from their counterparts in western countries. (full abstract in PDF)

(144977) 2005 EC127 is an V-/A-type inner-main-belt asteroid with a diameter of 0.6 +- 0.1 km. Asteroids of this size are believed to have rubble-pile structure, and, therefore, cannot have a rotation period shorter than 2.2 hours. However, our measurements show that asteroid 2005 EC127 completes one rotation in 1.65 +- 0.01 hours with a peak-to-peak light-curve variation of ~0.5 mag. Therefore, this asteroid is identified as a large super-fast rotator. Either a rubble-pile asteroid with a bulk density of ~6 g cm^-3 or an asteroid with an internal cohesion of 47 +- 30 Pa can explain 2005 EC127. However, the scenario of high bulk density is very unlikely for asteroids. To date, only six large super-fast rotators, including 2005 EC127, have been reported, and this number is very small when compared with the much more numerous fast rotators. We also note that none of the six reported large SFRs are classified as C-type asteroids.

Dormant or near-dormant short-period comets can unexpectedly regain ability to eject dust. In many known cases, the resurrection is short-lived and lasts less than one orbit. However, it is possible that some resurrected comets can remain active in latter perihelion passages. We search the archival images of various facilities to look for these "reactivated" comets. We identify two candidates, 297P/Beshore and 332P/Ikeya-Murakami, both of which were found to be inactive or weakly active in the previous orbit before their discovery. We derive a reactivation rate of $\sim0.007~\mathrm{comet}^{-1}~\mathrm{orbit}^{-1}$, which implies that typical short-period comets only become temporary dormant for less than a few times. Smaller comets are prone to rotational instability and may undergo temporary dormancy more frequently. Next generation high-cadence surveys may find more reactivation events of these comets.

Encke-type comet 332P/Ikeya-Murakami is experiencing cascading fragmentation events during its 2016 apparition. It is likely the first splitting Encke-type comet ever being observed. A nongravitational solution to the astrometry reveals a statistical detection of the radial and transverse nongravitational parameters, $A_{1} = \left(1.54 \pm 0.39\right) \times 10^{-8}$ AU day$^{-2}$, and $A_{2} = \left(7.19 \pm 1.92\right) \times 10^{-9}$ AU day$^{-2}$, respectively, which implies a nucleus erosion rate of $\left(0.91 \pm 0.17\right)$% per orbital revolution. The mass-loss rate likely has to be supported by a much larger fraction of an active surface area than known cases of short-period comets; it may be relevant to the ongoing fragmentation. We failed to detect any serendipitous pre-discovery observations of the comet in archival data from major sky surveys, whereby we infer that 332P used to be largely inactive, and is perhaps among the few short-period comets which have been reactivated from weakly active or dormant states. We therefore constrain an upper limit to the nucleus size as $2.0 \pm 0.2$ km in radius. A search for small bodies in similar orbits to that of 332P reveals comet P/2010 B2 (WISE) as the best candidate. From an empirical generalised Jupiter-family (Encke-type included) comet population model, we estimate the likelihood of chance alignment of the 332P--P/2010 B2 pair to be 1 in 33, a small number indicative of a genetic linkage between the two comets on a statistical basis. The pair possibly originated from a common progenitor which underwent a disintegration event well before the twentieth century.

We describe 2016 January to April observations of the fragments of 332P/Ikeya-Murakami, a comet earlier observed in a 2010 October outburst (Ishiguro et al 2014). We present photometry of the fragments, and perform simulations to infer the time of breakup. We argue that the eastern-most rapidly brightening fragment ($F4$) best corresponds to the original nucleus, rather than the initial bright fragment $F1$. We compute radial and tangential non-gravitational parameters, $A_1 = (1.5 \pm 0.4) \times 10^{-8}$ AU day$^{-2}$ and $(7.2 \pm 1.9) \times 10^{-9}$ AU day$^{-2}$; both are consistent with zero at the $4\sigma$ level. Monte Carlo simulations indicate that the fragments were emitted on the outbound journey well after the 2010 outburst, with bright fragment $F1$ splitting in mid--2013 and the fainter fragments within months of the 2016 January recovery. Western fragment $F7$ is the oldest, dating from 2011. We suggest that the delayed onset of the splitting is consistent with a self-propagating crystallization of water ice.

Dormant comets in the near-Earth object (NEO) population are thought to be involved in the terrestrial accretion of water and organic materials. Identification of dormant comets is difficult as they are observationally indistinguishable from their asteroidal counterparts, however they may have produced dust during their final active stages which potentially are detectable today as weak meteor showers at the Earth. Here we present the result of a reconnaissance survey looking for dormant comets using 13~567~542 meteor orbits measured by the Canadian Meteor Orbit Radar (CMOR). We simulate the dynamical evolution of the hypothetical meteoroid streams originated from 407 near-Earth asteroids in cometary orbits (NEACOs) that resemble orbital characteristics of Jupiter-family comets (JFCs). Out of the 44 hypothetical showers that are predicted to be detectable by CMOR, we identify 5 positive detections that are statistically unlikely to be chance associations, including 3 previously known associations. This translates to a lower limit to the dormant comet fraction of $2.0\pm1.7\%$ in the NEO population and a dormancy rate of $\sim 10^{-5}~\mathrm{yr^{-1}}$ per comet. The low dormancy rate confirms disruption and dynamical removal as the dominant end state for near-Earth JFCs. We also predict the existence of a significant number of meteoroid streams whose parents have already been disrupted or dynamically removed.

Previous studies have revealed Jupiter-family comet 252P/LINEAR as a comet that was recently transported into the near-Earth object (NEO) region in $\sim1800$~AD yet only being weakly active. In this Letter, we examine the "formed (almost) dead" hypothesis for 252P/LINEAR using both dynamical and observational approaches. By statistically examining the dynamical evolution of 252P/LINEAR over a period of $10^7$~years, we find the median elapsed residency in the NEO region to be $4\times10^2$~years which highlights the likelihood of 252P/LINEAR as an (almost) first-time NEO. With available cometary and meteor observations, we find the dust production rate of 252P/LINEAR to be at the order of $10^6$~kg per orbit since its entry to the NEO region. These two lines of evidence support the hypothesis that the comet was likely to have formed in a volatile-poor environment. Cometary and meteor observations during the comet's unprecedented close approach to the Earth around 2016 Mar. 21 would be useful for the understanding of the surface and evolutionary properties of this unique comet.

The Ngari (Ali) prefecture of Tibet, one of the highest areas in the world, has recently emerged as a promising site for future astronomical observation. Here we use 31 years of reanalysis data from the Climate Forecast System Reanalysis (CFSR) to examine the astroclimatology of Ngari, using the recently-erected Ali Observatory at Shiquanhe (5~047~m above mean sea level) as the representative site. We find the percentage of photometric night, median atmospheric seeing and median precipitable water vapor (PWV) of the Shiquanhe site to be $57\%$, $0.8"$ and 2.5~mm, comparable some of the world's best astronomical observatories. Additional calculation supports the Shiquanhe region as one of the better sites for astronomical observations over the Tibetan Plateau. Based on the studies taken at comparable environment at Atacama, extraordinary observing condition may be possible at the few vehicle-accessible 6~000~m heights in the Shiquanhe region. Such possibility should be thoroughly investigated in future.

Jupiter-family comet 15P/Finlay has been reportedly quiet in activity for over a century but has harbored two outbursts during its 2014/2015 perihelion passage. Here we present an analysis of these two outbursts using a set of cometary observations. The outbursts took place between 2014 Dec. 15.4--16.0 UT and 2015 Jan. 15.5--16.0 UT as constrained by ground-based and spacecraft observations. We find a characteristic ejection speed of $V_0=300$ to $650 \mathrm{m \cdot s^{-1}}$ for the ejecta of the first outburst and $V_0=550$ to $750 \mathrm{m \cdot s^{-1}}$ for that of the second outburst using a Monte Carlo dust model. The mass of the ejecta is calculated to be $M_\mathrm{d}=2$ to $3\times10^5 \mathrm{kg}$ for the first outburst and $M_\mathrm{d}=4$ to $5\times10^5 \mathrm{kg}$ for the second outburst, corresponds to less than $10^{-7}$ of the nucleus mass. The specific energy of the two outbursts is found to be $0.3$ to $2\times10^5 \mathrm{J \cdot kg^{-1}}$. We also revisit the long-standing puzzle of the non-detection of the hypothetical Finlayid meteor shower by performing a cued search using the 13-year data from the Canadian Meteor Orbit Radar, which does not reveal any positives. The Earth will pass the 2014/2015 outburst ejecta around 2021 Oct. 6 at 22 h UT to Oct. 7 at 1 h UT, with a chance for some significant meteor activity in the radio range, which may provide further clues to the Finlayid puzzle. A southerly radiant in the constellation of Ara will favor the observers in the southern tip of Africa.

We present studies of C/2015 D1 (SOHO), the first sunskirting comet ever seen from ground stations over the past half century. The Solar and Heliospheric Observatory (SOHO) witnessed its peculiar light curve with a huge dip followed by a flareup around perihelion: the dip was likely caused by sublimation of olivines, directly evidenced by a coincident temporary disappearance of the tail. The flareup likely reflects a disintegration event, which we suggest was triggered by intense thermal stress established within the nucleus interior. Photometric data reveal an increasingly dusty coma, indicative of volatile depletion. A catastrophic mass loss rate of $\sim$10$^{5}$ kg s$^{-1}$ around perihelion was seen. Ground-based Xingming Observatory spotted the post-perihelion debris cloud. Our morphological simulations of post-perihelion images find newly released dust grains of size $a \gtrsim 10$ $\mu$m in radius, however, a temporal increase in $a_{\min}$ was also witnessed, possibly due to swift dispersions of smaller grains swept away by radiation forces without replenishment. Together with the fading profile of the light curve, a power law dust size distribution with index $\gamma = 3.2 \pm 0.1$ is derived. We detected no active remaining cometary nuclei over $\sim$0.1 km in radius in post-perihelion images acquired at Lowell Observatory. Applying radial non-gravitational parameter, $\mathcal{A}_{1} = \left(1.209 \pm 0.118 \right) \times 10^{-6}$ AU day$^{-2}$, from an isothermal water-ice sublimation model to the SOHO astrometry significantly reduces residuals and sinusoidal trends in the orbit determination. The nucleus mass $\sim$10$^{8}$--10$^{9}$ kg, and the radius $\sim$50--150 m (bulk density $\rho_{\mathrm{d}} = 0.4$ g cm$^{-3}$ assumed) before the disintegration are deduced from the photometric data; consistent results were determined from the non-gravitational effects.

It is speculated that some weakly active comets may be transitional objects between active and dormant comets. These objects are at a unique stage of the evolution of cometary nuclei, as they are still identifiable as active comets, in contrast to inactive comets that are observationally indistinguishable from low albedo asteroids. In this paper, we present a synthesis of comet and meteor observations of Jupiter-family comet 209P/LINEAR, one of the most weakly active comets recorded to-date. Images taken by the Xingming 0.35-m telescope and the Gemini Flamingo-2 camera are modeled by a Monte Carlo dust model, which yields a low dust ejection speed ($1/10$ of that of moderately active comets), dominance of large dust grains, and a low dust production of $0.4~\mathrm{kg \cdot s^{-1}}$ at 19~d after the 2014 perihelion passage. We also find a reddish nucleus of 209P/LINEAR that is similar to D-type asteroids and most Trojan asteroids. Meteor observations with the Canadian Meteor Orbit Radar (CMOR), coupled with meteoroid stream modeling, suggest a low dust production of the parent over the past few hundred orbits, although there are hints of a some temporary increase in activity in the 18th century. Dynamical simulations indicate 209P/LINEAR may have resided in a stable near-Earth orbit for $\sim 10^4$~yr, which is significantly longer than typical JFCs. All these lines of evidence imply that 209P/LINEAR as an aging comet quietly exhausting its remaining near surface volatiles. We also compare 209P/LINEAR to other low activity comets, where evidence for a diversity of the origin of low activity is seen.

The Kreutz family of sungrazing comets contains over 2~000 known members, many of which are believed to be under $\sim 100$~m sizes (\textit{mini~comets}) and have only been studied at small heliocentric distances ($r_\mathrm{H}$) with space-based SOHO/STEREO spacecraft. To understand the brightening process of mini Kreutz comets, we conducted a survey using CFHT/MegaCam at moderate $r_\mathrm{H}$ guided by SOHO/STEREO observations. We identify two comets that should be in our search area but are not detected, indicating that the comets have either followed a steeper brightening rate within the previously-reported rapid brighten stage (the \textit{brightening~burst}), or the brightening burst starts earlier than expected. We present a composite analysis of pre-perihelion light-curves of five Kreutz comets that covered to $\sim 1$~AU. We observe a significant diversity in the light-curves that can be used to grossly classify them into two types: C/Ikeya-Seki and C/SWAN follow the canonical $r_\mathrm{H}^{-4}$ while the others follow $r_\mathrm{H}^{-7}$. In particular, C/SWAN seems to have undergone an outburst ($\Delta m>5$~mag) or a rapid brightening ($n\gtrsim11$) between $r_\mathrm{H}=1.06$~AU to $0.52$~AU, and shows hints of structural/compositional differences compared to other bright Kreutz comets. We also find evidence that the Kreutz comets as a population lose their mass less efficiently than the dynamically new comet, C/ISON, and are relatively devoid of species that drive C/ISON's activity at large $r_\mathrm{H}$. Concurrent observations of C/STEREO in different wavelengths also suggest that a blue-ward species such as CN may be the main driver for brightening burst instead of previously-thought sodium.

The dynamically new comet, C/2013 A1 (Siding Spring), is to make a close approach to Mars on 2014 October 19 at 18:30 UT at a distance of 40+/-1 Martian radius. Such extremely rare event offers a precious opportunity for the spacecrafts on Mars to closely study a dynamically new comet itself as well as the planet-comet interaction. Meanwhile, the high speed meteoroids released from C/Siding Spring also pose a threat to physically damage the spacecrafts. Here we present our observations and modeling results of C/Siding Spring to characterize the comet and assess the risk posed to the spacecrafts on Mars. We find that the optical tail of C/Siding Spring is dominated by larger particles at the time of the observation. Synchrone simulation suggests that the comet was already active in late 2012 when it was more than 7 AU from the Sun. By parameterizing the dust activity with a semi-analytic model, we find that the ejection speed of C/Siding Spring is comparable to comets such as the target of the Rosetta mission, 67P/Churyumov-Gerasimenko. Under nominal situation, the simulated dust cone will miss the planet by about 20 Martian radius. At the extreme ends of uncertainties, the simulated dust cone will engulf Mars, but the meteoric influx at Mars is still comparable to the nominal sporadic influx, seemly indicating that intense and enduring meteoroid bombardment due to C/Siding Spring is unlikely. Further simulation also suggests that gravitational disruption of the dust tail may be significant enough to be observable at Earth.

An unexpected intense outburst of the Draconid meteor shower was detected by the Canadian Meteor Orbit Radar (CMOR) on October 8, 2012. The peak flux occurred at ~16:40 UT on October 8 with a maximum of 2.4 +/- 0.3 hr-1 km-2 (appropriate to meteoroid mass larger than 10-7 kg), equivalent to a ZHRmax = 9000 +/- 1000 using 5-minute intervals, using a mass distribution index of s = 1.88 +/- 0.01 as determined from the amplitude distribution of underdense Draconid echoes. This makes the out- burst among the strongest Draconid returns since 1946 and the highest flux shower since the 1966 Leonid meteor storm, assuming a constant power-law distribution holds from radar to visual meteoroid sizes. The weighted mean geocentric radiant in the time interval of 15-19h UT, Oct 8, 2012 was {\alpha}g = 262.4 +/- 0.1 deg, {\delta}g = 55.7 +/- 0.1 deg (epoch J2000.0). Visual observers also reported increased activity around the peak time, but with a much lower rate (ZHR 200), suggesting that the magnitude-cumulative num- ber relationship is not a simple power-law. Ablation modeling of the observed meteors as a population does not yield a unique solution for the grain size and distribution of Draconid meteoroids, but is consistent with a typical Draconid meteoroid of mtotal between 10-6 to 10-4 kg being composed of 10 - 100 grains. Dynamical simulations indicate that the outburst was caused by dust particles released during the 1966 per- ihelion passage of the parent comet, 21P/Giacobini-Zinner, although there are dis- crepancies between the modelled and observed timing of the encounter, presumably caused by approaches of the comet to Jupiter during 1966-1972. Based on the results of our dynamical simulation, we predict possible increased activity of the Draconid meteor shower in 2018, 2019, 2021 and 2025.

Previous studies have suggested that Comet 209P/LINEAR may produce strong me- teor activity on Earth on 2014 May 24; however, exact timing and activity level is difficult to estimate due to the limited physical observations of the comet. Here we reanalyze the optical observations of 209P/LINEAR obtained during its 2009 appari- tion. We find that the comet is relatively depleted in dust production, with Af {\rho} at 1 cm level within eight months around its perihelion. This feature suggested that this comet may be currently transitioning from typical comet to a dormant comet. Syn- dyne simulation shows that the optical cometary tail is dominated by larger particles with {\beta} ~ 0.003. Numerical simulations of the cometary dust trails confirm the arrival of particles on 2014 May 24 from some of the 1798-1979 trails. The nominal radiant is at RA 122 +/- 1 deg, Dec 79 +/- 1 deg (J2000) in the constellation of Camelopardalis. Given that the comet is found to be depleted in dust production, we concluded that a meteor storm (ZHR>=1000) may be unlikely. However, our simulation also shows that the size distribution of the arrived particles is skewed strongly to larger particles. Coupling with the result of syndyne simulation, we think that the event, if detectable, may be dominated by bright meteors. We encourage observers to monitor the expected meteor event as it will provide us with rare direct information on the dynamical history of 209P/LINEAR which is otherwise irretrievably lost.

A strong outburst of the October Draconid meteor shower was predicted for October 8, 2011. Here we present the observations obtained by the Canadian Meteor Orbit Radar (CMOR) during the 2011 outburst. CMOR recorded 61 multi-station Draconid echoes and 179 single-station overdense Draconid echoes (covering the magnitude range of +3 <= MV <= +7) between 16-20h UT on October 8, 2011. The mean radiant for the outburst was determined to be a_g = 261.9o +/- 0.3o, d_g = +55.3o +/- 0.3o (J2000) from observations of the underdense multi-station echoes. This radiant location agrees with model predictions to ~1o . The determined geocentric velocity was found to be ~10-15% lower than the model value (17.0 - 19.1 km s^-1 versus 20.4 km s^-1), a discrepancy we attribute to undercorrection for atmospheric deceleration of low density Draconid meteoroids as well as to poor radar radiant geometry during the outburst peak. The mass index at the time of the outburst was determined to be ~ 1.75 using the amplitude distribution of underdense echoes...

To explore the issue of performing a non-interactive numerical weather forecast with an operational global model in assist of astronomical observing, we use the Xu-Randall cloud scheme and the Trinquet-Vernin AXP seeing model with the global numerical output from the Global Forecast System to generate 3-72h forecasts for cloud coverage and atmospheric seeing, and compare them with sequence observations from 9 sites from different regions of the world with different climatic background in the period of January 2008 to December 2009. The evaluation shows that the proportion of prefect forecast of cloud cover forecast varies from ~50% to ~85%. The probability of cloud detection is estimated to be around ~30% to ~90%, while the false alarm rate is generally moderate and is much lower than the probability of detection in most cases. The seeing forecast has a moderate mean difference (absolute mean difference <0.3" in most cases) and root-mean-square-error or RMSE (0.2"-0.4" in most cases) comparing with the observation. The probability of forecast with <30% error varies between 40% to 60% for entire atmosphere forecast and 40% to 50% for free atmosphere forecast for almost all sites, which being placed in the better cluster among major seeing models. However, the forecast errors are quite large for a few particular sites. Further analysis suggests that the error might primarily be caused by the poor capability of GFS/AXP model to simulate the effect of turbulence near ground and on sub-kilometer scale. In all, although the quality of the GFS model forecast may not be comparable with the human-participated forecast at this moment, our study has illustrated its suitability for basic observing reference, and has proposed its potential to gain better performance with additional efforts on model refinement.

We present the results of BVRI photometry and classification of 53 unusual asteroids, including 35 near-Earth asteroids (NEAs), 6 high eccentricity/inclination asteroids, and 12 recently-identified asteroid-pair candidates. Most of these asteroids were not reportedly classified prior to this work. For the few asteroids that have been previously studied, the results are generally in rough agreement. In addition, we merge the results from several photometric/spectroscopic surveys to create a largest-ever sample with 449 spectrally classified NEAs for statistical analysis. We identify a "transition point" of the relative number of C/X-like and S-like NEAs at H~18<=>D~1km with confidence level at ~95% or higher. We find that the C/X-like:S-like ratio for 18<=H<22 is about two times higher than that of H<18 (0.33+/-0.04 versus 0.17+/-0.02), virtually supporting the hypothesis that smaller NEAs generally have less weathered surface (therefore, less reddish appearance) caused by younger collision ages.