Comet nuclei: Morphology and implied processes of surface modification

Surface morphology and related issues for nuclei of three comets: Halley, Borrelly and Wild 2, are considered in the paper. Joint consideration of publications and results of our analysis of the comets’ images led to conclusions, partly new, partly repeating conclusions published by other researchers. It was found that typical for all three nuclei is the presence of rather flat areas: floors of craters and other depressions, mesas and terraces. This implies that flattening surfaces or planation is a process typical for the comet nuclei. Planation seems to work through the sublimation-driven slope collapse and retreat. This requires effective sublimation so this process should work only when a comet is close to the Sun and if on the nucleus there are starting slopes, steep and high enough to support the “long-distance” avalanching of the collapsing material. If the surface had no starting slopes, then instead of planation, the formation of pitted-and-hilly surfaces should occur. An example of this could be the mottled terrain of the Borelly nucleus. Both ways of the sublimational evolution on the nucleus surface should lead to accumulation of cometary regolith. The thickness of the degassed regolith is not known, but it is obvious that in surface depressions, including the flat-floor ones, it should be larger compared with nondepression areas. This may have implications for the in situ study of comets by the Deep Impact and Rosetta missions.Our morphological analysis puts constraints on the applicability of the popular “rubble-pile comet nucleus” hypothesis (Weissman, 1986. Are cometery nuclei primordial rubble piles? Nature 320, 242-244.). We believe that the rubble pile hypothesis can be applicable to the blocky Halley nucleus. The Borelly and Wild 2 nuclei also could be rubble piles. But in these cases the “rubbles” have to be either smaller than 30-50 m (a requirement to keep lineament geometry close to ideal), or larger than 1-2 km (a requirement to form the rather extended smooth, flat surfaces of mesa tops and crater floors). Another option is that the Borelly and Wild 2 nuclei are not rubble piles.In relation to surface morphology we suggest that three end-member types of the comet nuclei may exist: (1) impact cratered “pristine” bodies, (2) non-cratered fragments of catastrophic disruption, and (3) highly Sun-ablated bodies. In this threefold classification, the Wild 2 nucleus is partially ablated primarily cratered body. Borrelly is significantly ablated and could be either primarily cratered or not-cratered fragment. Halley is certainly partially ablated but with the available images it is difficult to say if remnants of impact craters do exist on it.Recently published observations and early results of analysis of the Tempel 1 nucleus images taken by Deep Impact mission are in agreement with our conclusions on the processes responsible for the Halley, Borrelly and Wild 2 nuclei morphologies. In particular, we have now more grounds to suggest that decrease in crater numbers and increase of the role of smooth flat surfaces in the sequence Wild 2?Tempel 1?Borelli reflects a progress in the sublimational degradation of the nucleus surface during comet passages close to the Sun.     

Single star scidar: atmospheric parameters profiling using the simulated annealing algorithm

In this paper, we present a new method to estimate, for each turbulent layer labelled i , the horizontal wind speed   v ( h _i )  , the standard deviation of the horizontal wind speed fluctuations  σ _v ( h_i )  and the integrated value of   C ² _n   over the thickness  Δ h_i   of the turbulent layer   C ² _n ( h_i )Δ h_i   , where   h_i   is the altitude of the turbulent layer. These parameters are extracted from single star scintillation spatiotemporal cross-correlation functions of atmospheric speckles obtained within the generalized mode. This method is based on the simulated annealing algorithm to find the optimal solution required to solve the problem. Astrophysics parameters for adaptive optics are also calculated using   C ² _n ( h_i )  and   v ( h_i )  values. The results of other techniques support this new method.     

OSIRIS: AO-assisted integral-field spectroscopy at the Keck Observatory

OSIRIS (OH-Suppressing Infra-Red Integral-field Spectrograph) is a new facility instrument for the Keck Observatory. After seeing first light in February 2005, OSIRIS is currently undergoing commissioning. OSIRIS provides the capability of performing three-dimensional spectroscopy in the near-infrared z, J, H, and K bands at the resolution limit of the Keck II telescope, which is equipped with adaptive optics and a laser guide star. The science case for OSIRIS is summarized, and the instrument and associated data reduction software are described.     

Wavelet Analysis of Space Solar Telescope Images

The scientific satellite SST (Space Solar Telescope) is an important research project strongly supported by the Chinese Academy of Sciences. Every day, SST acquires 50 GB of data (after processing) but only 10GB can be transmitted to the ground because of limited time of satellite passage and limited channel volume. Therefore, the data must be compressed before transmission. Wavelets analysis is a new technique developed over the last 10 years, with great potential of application. We start with a brief introduction to the essential principles of wavelet analysis, and then describe the main idea of embedded zerotree wavelet coding, used for compressing the SST images. The results show that this coding is adequate for the job.     

一次西南低涡过程的地面场中尺度分析

本文运用带通滤波分析方法,对四川省1981年7月由西南低涡等天气系统引起的大暴雨过程作了地面气象要素场的中尺度诊断分析。结果表明,带通滤波后的场具有明显的中尺度特征,其中一些物理量场的适当配置与6小时后的暴雨中心有较好的对应关系。说明了此带通滤波方法在实际业务预报中具有一定的使用价值。     

Modelling Gas And Dust Release From Comet Hale–Bopp

Numerical simulations of the evolving activity of comet Hale-Bopp are presented, assuming a porous, spherical nucleus, 20 km in radius, made of dust and gas-laden amorphous ice. The main effects included are: crystallization of amorphous ice and release of occluded gas, condensation, sublimation and flow of gases through the pores, changing pore sizes, and flow of dust grains. The model parameters, such as initial pore size and porosity, emissivity, dust grain size, are varied in order to match the observed activity. In all cases, a sharp rise in the activity of the nucleus occurs at a large heliocentric distance pre-perihelion, marked by a few orders of magnitude increase in the CO and the CO₂ fluxes and in the rate of dust emission. This is due to the onset of crystallization, advancing down to a few meters below the surface, accompanied by release of the trapped gases. A period of sustained, but variable, activity ensues. The emission of water molecules is found to surpass that of CO at a heliocentric distance of 3 AU. Thereafter the activity is largely determined by the behaviour of the dust. If a dust mantle is allowed to build up, the water production rate does not increase dramatically towards perihelion; if most of the dust is ejected, the surface activity increases rapidly, producing a very bright comet. This revised version was published online in July 2006 with corrections to the Cover Date.     

台湾海峡西岸重力场变化梯度场图像分析

讨论了应用曲面拟合技术和利用离散的重力观测点的重力变化对特定研究区的重力场变化进行数值模拟以获取该区重力变化空间梯度场图像的具体方法;并应用该方法获得了2004年3月福建安溪3.6级地震前的重力变化空间梯度场图像,发现空间梯度场图像能清晰地反映出地震前重力场群体性变化的差异性;利用此方法对2005年6月~2006年3月间台湾海峡西岸重力场变化进行数值模拟,获得其空间梯度场图像。     

Correcting systematic effects in a large set of photometric light curves

We suggest a new algorithm to remove systematic effects in a large set of light curves obtained by a photometric survey. The algorithm can remove systematic effects, such as those associated with atmospheric extinction, detector efficiency, or point spread function changes over the detector. The algorithm works without any prior knowledge of the effects, as long as they linearly appear in many stars of the sample. The approach, which was originally developed to remove atmospheric extinction effects, is based on a lower rank approximation of matrices, an approach which has already been suggested and used in chemometrics, for example. The proposed algorithm is especially useful in cases where the uncertainties of the measurements are unequal. For equal uncertainties, the algorithm reduces to the Principal Component Analysis (PCA) algorithm. We present a simulation to demonstrate the effectiveness of the proposed algorithm and we point out its potential, in the search for transit candidates in particular.