Comparison of centroid computation algorithms in a Shack–Hartmann sensor

Analytical theory is combined with extensive numerical simulations to compare different flavours of centroiding algorithms: thresholding, weighted centroid, correlation, quad cell (QC). For each method, optimal parameters are defined in function of photon flux, readout noise and turbulence level. We find that at very low flux the noise of QC and weighted centroid leads the best result, but the latter method can provide linear and optimal response if the weight follows spot displacements. Both methods can work with average flux as low as 10 photons per subaperture under a readout noise of three electrons. At high-flux levels, the dominant errors come from non-linearity of response, from spot truncations and distortions and from detector pixel sampling. It is shown that at high flux, centre of gravity approaches and correlation methods are equivalent (and provide better results than QC estimator) as soon as their parameters are optimized. Finally, examples of applications are given to illustrate the results obtained in the paper.     

Searching for water ice on 47171 1999 TC36, 1998 SG35, and 2000 QC243: ESO large program on TNOs and centaurs

Transneptunian objects and Centaurs are supposed to be among the most pristine bodies of the Solar System. To investigate their physical properties and their surface composition, an ESO large program at the Very Large Telescope was carried out. In this paper we present photometric and spectroscopic near-infrared data of two Centaurs (1998 SG35 and 2000 QC243) and one transneptunian object (47171 1999 TC36). For 47171 1999 TC36 and 1998 SG35 visible photometry is also presented. Models of the surface composition of these objects are presented and discussed. By including a small percentage of water ice in our geographical mixtures, we obtain a better agreement with the observations in the H and K bands.     

The Greenwich Photo-heliographic Results (1874 – 1976): Initial Corrections to the Printed Publications

A new sunspot and faculae digital dataset for the interval 1874?–?1955 has been prepared under the auspices of the NOAA National Geophysical Data Center (NGDC). This digital dataset contains measurements of the positions and areas of both sunspots and faculae published initially by the Royal Observatory, Greenwich, and subsequently by the Royal Greenwich Observatory (RGO), under the title Greenwich Photo-heliographic Results (GPR), 1874?–?1976. Quality control (QC) procedures based on logical consistency have been used to identify the more obvious errors in the RGO publications. Typical examples of identifiable errors are North versus South errors in specifying heliographic latitude, errors in specifying heliographic (Carrington) longitude, errors in the dates and times, errors in sunspot group numbers, arithmetic errors in the summation process, and the occasional omission of solar ephemerides. Although the number of errors in the RGO publications is remarkably small, an initial table of necessary corrections is provided for the interval 1874?–?1917. Moreover, as noted in the preceding companion papers, the existence of two independently prepared digital datasets, which both contain information on sunspot positions and areas, makes it possible to outline a preliminary strategy for the development of an even more accurate digital dataset. Further work is in progress to generate an extremely reliable sunspot digital dataset, based on the long programme of solar observations supported first by the Royal Observatory, Greenwich, and then by the Royal Greenwich Observatory.     

Photometric study of the major satellites of Uranus

In this paper, we analyze the results of ground-based and space-born photometric observations of the major satellites of Uranus—Miranda, Ariel, Umbriel, Titania, and Oberon. All sets of photometric observations of the satellites available in the literature were examined for uniformity and systematic differences and summarized to a unified set by wavelength ranging from 0.25 to 2.4 μm. This set covers the interval of phase angles from 0.034° to 35°. The compound phase curves of brightness of the satellites in the spectral bands at 0.25, 0.41, 0.48, 0.56, 0.75, 0.91, 1.4, and 1.8 μm, which include a pronounced opposition surge and linear part, were constructed. For each satellite, the geometric albedo was found in different spectral bands taking into account the brightness opposition effect, and its spectral dependence was studied. It has been shown that the reflectance of the satellites linearly depends on the wavelength at different phase angles, but has different spectral gradients. The parameters of the phase functions of brightness, including the amplitude and the angular width of the brightness opposition surge, the phase coefficient, and the phase angle at which the nonlinear increase in brightness starts, were determined and their dependences on wavelength and geometric albedo were analyzed. Our investigations show that, in their optical properties, the satellites Miranda and Ariel, Titania and Oberon, and Umbriel present three types of surfaces. The observed parameters of the brightness opposition effect for the Uranian satellites, some ice satellites of Jupiter and Saturn, and the E-and S-type asteroids are analyzed and compared within the framework of the coherent backscattering and mutual shadowing mechanisms.     

Report of the IAU/IAG Working Group on cartographic coordinates and rotational elements: 2006

Every three years the IAU/IAG Working Group on Cartographic Coordinates and Rotational Elements revises tables giving the directions of the poles of rotation and the prime meridians of the planets, satellites, minor planets, and comets. This report introduces improved values for the pole and rotation rate of Pluto, Charon, and Phoebe, the pole of Jupiter, the sizes and shapes of Saturn satellites and Charon, and the poles, rotation rates, and sizes of some minor planets and comets. A high precision realization for the pole and rotation rate of the Moon is provided. The expression for the Sun’s rotation has been changed to be consistent with the planets and to account for light travel time     

Peripheral Structures of Planetary Systems

We analyze the conditions for the formation and time evolution of peripheral comet structures of solar-type planetary systems. In the Solar system, these include the Kuiper belt, the Oort cloud, the comet spear, and the Galactic comet ring that marks the Galactic orbit of the Sun. We consider the role of the viscosity of a protoplanetary gas–dust disk, major planets, field stars, globular clusters, giant molecular clouds, and the Galactic gravitational field in the formation of these peripheral structures marked by comets and asteroids. We give a list of the closest past and future passages of neighboring stars through the solar Oort cloud that perturb the motion of its comets and, thus, contribute to the enhancement of its cometary activity, on the one hand, and to the replenishment of the solar comet spear with new members, on the other hand.     

The analysis of the characteristics of variation and change patterns of sunspot cycles

The paper focus on the variation character of sunspot number and solar cycles based on the new version sunspot number (SSN) data. According to seven main variables describing solar cycles, including peak value, the length of cycle, the length of ascending phase, the ratio of the ascending time to the descending time, slope, half width, and area under the curve of solar cycle, clustering, principal component and factor analysis, are applied to analyze variation characteristic and patterns of the 24 solar cycles. We cluster these 24 cycles to find groups in these solar cycles, and search for the main factor determining strength, length and occurrence time of the peak, and the furthest cycle from the average. The cycles within a cluster will be similar or related to one another and different from or unrelated to the cycles in other clusters. These results could help us search for similar cycles conveniently, obtain the understanding of the characteristics of solar cycle variation and analysis of sunspot number change and evolution characteristics, and analyze the origin and the variation mechanism of solar cycle.     

The international planetary patrol program: An assessment of the first three years

Day-to-day and hour-to-hour changes in the large-scale atmospheric and surface features of the planets can now be studied more effectively than previously possible. Since 1969 a network of observatories has obtained almost uninterrupted photographic coverage during all apparitions of Jupiter and Mars, plus some of Venus. Patrol films and catalogues are available to the scientific community. Recent or current analyses include the distribution and motion of clouds on Mars, the development and decay of Martian dust storms, the seasonal, diurnal and random fluctuations in contrast between adjacent light and dark regions on Mars, the detection of vertical shear in the Jovian atmosphere, the longitudinal oscillation of the Red Spot, the dependence of rotation period on xenographic latitude and on time, the eruption and spread of SEB disturbances, and the retrograde circulation of the Venus cloud deck.     

Assessing the performance of molecular gas clump identification algorithms

The detection of clumps(cores) in molecular clouds is an important issue in sub-millimetre astronomy. However, the completeness of the identification and the accuracy of the returned parameters of the automated clump identification algorithms are still not clear. In this work, we test the performance and bias of the GaussClumps, ClumpFind, FellWalker, Reinhold, and Dendrograms algorithms in identifying simulated clumps. By designing the simulated clumps with various sizes, peak brightness, and crowdedness, we investigate the characteristics of the algorithms and their performance. In the aspect of detection completeness, FellWalker, Dendrograms, and GaussClumps are the first, second, and third best algorithms, respectively. The numbers of correct identifications of the six algorithms gradually increase as the size and signal-to-noise ratio(SNRs) of the simulated clumps increase and they decrease as the crowdedness increases. In the aspect of the accuracy of retrieved parameters, FellWalker and Dendrograms exhibit better performance than the other algorithms. The average deviations in clump parameters for all algorithms gradually increase as the size and SNR of clumps increase. Most of the algorithms except FellWalker exhibit significant deviation in extracting the total flux of clumps. Taken together, FellWalker, GaussClumps,and Dendrograms exhibit the best performance in detection completeness and extracting parameters. The deviation in virial parameter for the six algorithms is relatively low. When applying the six algorithms to the clump identification for the Rosette molecular cloud, ClumpFind1994, ClumpFind2006, GaussClumps,FellWalker, and Reinhold exhibit performance that is consistent with the results from the simulated test.     

从KARST到FAST---中国天眼的概念起源与决策过程

在大科学时代,大科学工程的建设日益复杂,设计者需综合考虑设备性能、技术储备、经费、风险、环境等因素才能进行合理决策.在国内外天文学发展的背景中,梳理了500 m口径球面射电望远镜(Five-hundredmeter Aperture Spherical radio Telescope, FAST)概念形成的过程以及FAST团队在其中所做的决策,包括大望远镜中国方案的设想、先导单元的提出、主动反射面技术的采用等,并对决策过程进行了分析探讨. FAST诞生于我国与国际天文学发展的互动与融合进程,实现了从跟进到占据先机的转变,可为在相关基础薄弱的领域建设大科学工程提供参考.     

Multispectral observations complementary to the study of high-energy solar phenomena

Non-thermal phenomena on the Sun are characterized by the transient acceleration of electrons and ions to energies ranging from several keV to tens of GeV, and the impulsive heating of plasma to temperatures exceeding 5 × 10⁷ K. These energetic processes result in the emission of a broad spectrum of electromagnetic radiation and of high-energy neutrons, as well as the escape of high energy electrons and ions from the acceleration region. The determination of the energy spectrum, polarization, and spatial distribution of these emissions, which contain detailed information on the acceleration and heating process, and the conditions at the sites at which this energy is generated and dissipated, is the principal objective of high-energy solar studies.The study of the evolution of magnetic structures in the solar convection zone and atmosphere which underlie the metastable conditions which precede these energetic processes, of the conditions that trigger the release of energy, and of the impact of the energy released on the solar atmosphere, is most effectively carried out by observations of thermal and quasi-thermal phenomena which precede, coincide with, and follow the impulsive acceleration and heating event itself. Multispectral observations of the phenomena associated with non-thermal events on the Sun are reviewed, and the requirements for visible, ultraviolet, extreme ultraviolet, and soft X-ray observations which are necessary for future advances are briefly described.     

CELIAS - Charge, Element and Isotope Analysis System for SOHO

The CELIAS experiment on SOHO is designed to measure the mass, ionic charge and energy of the low and high speed solar wind, of suprathermal ions, and of low energy flare particles. Through analysis of the elemental and isotopic abundances, the ionic charge state, and the velocity distributions of ions originating in the solar atmosphere, the investigation focuses on the plasma processes on various temporal and spatial scales in the solar chromosphere, transition zone, and corona. CELIAS includes 3 mass- and charge-discriminating sensors based on the time-of-flight technique: CTOF for the elemental, charge and velocity distribution of the solar wind, MTOF for the elemental and isotopic composition of the solar wind, and STOF for the mass, charge and energy distribution of suprathermal ions. The instrument will provide detailed in situ diagnostics of the solar wind and of accelerated particles, which will complement the optical and spectroscopic investigations of the solar atmosphere on SOHO. CELIAS also contains a Solar Extreme Ultraviolet Monitor, SEM, which continously measures the EUV flux in a wide band of 17 – 70 nm, and a narrow band around the 30.4 nm He II line.Principal-InvestigatorPrincipal-Investigator for data phase     

Virial oscillations of celestial bodies: V. The structure of the potential and kinetic energies of a celestial body as a record of its creation history

The physical meaning of the terms of the potential and kinetic energy expressions, expanded by means of the density variation function for a nonuniform self-gravitating sphere, is discussed. The terms of the expansions represent the energy and the moment of inertia of the uniform sphere, the energy and the moment of inertia of the nonuniformities interacting with the uniform sphere, and the energy of the nonuniformities interacting with each other. It follows from the physical meaning of the above components of the energy structure, and also from the observational fact of the expansion of the Universe that the phase transition, notably, fusion of particles and nuclei and condensation of liquid and solid phases of the expanded matter accompanied by release of energy, must be the physical cause of initial thermal and gravitational instability of the matter. The released kinetic energy being constrained by the general motion of the expansion, develops regional and local turbulent (cyclonic) motion of the matter, which should be the second physical effect responsible for the creation of celestial bodies and their rotation.     

Thermal radiation from impact plumes

Plumes produced by the impacts of asteroids and comets consist of rock vapor and heated air. They emit visible light, ultraviolet, and infrared radiation, which can greatly affect the environment. We have carried out numerical simulations of the impacts of stony and cometary bodies with a diameter of 0.3, 1, and 3 km, which enter the atmosphere at various angles, using a hydrodynamic model supplemented by radiation transfer. We assumed that the cosmic object has no strength, and deforms, fragments, and vaporizes in the atmosphere. After the impact on the ground, the formation of craters and plumes was simulated, taking the internal friction of destroyed rocks and the trail formed in the atmosphere into account. The equation of radiative transfer, added to the equations of gas dynamics, was used in the approximation of radiative heat conduction or, if the Rosseland optical depth of a radiating volume of gas and vapor was less than unity, in the volume‐emission approximation. We used temperature and density distributions obtained in these simulations to calculate radiation fluxes on the Earth's surface by integrating the equation of radiative transfer along rays passing through a luminous region. We used tables of the equation of state of dunite and quartz (for stony impactors and a target) and air, as well as tables of absorption coefficients of air, vapor of ordinary chondrite, and vapor of cometary material. We have calculated the radiation impulse on the ground and the impact radiation efficiency (a ratio of thermal radiation energy incident on the ground to the kinetic energy of a body), which ranges from ~0.5% to ~9%, depending on the impactor size and the angle of entry into the atmosphere. Direct thermal radiation from fireballs and impact plumes, poses a great danger to people, animals, plants, and economic objects. After the impacts of asteroids at a speed of 20 km s^?1 at an angle of 45°, a fire can occur at a distance of 250 km if the asteroid has a diameter of 0.3 km, and at a distance of 2000 km if the diameter is 3 km.     

The composition,structure, temperature and dynamics of the upper thermosphere in the polar regions during October to December 1981

During the period October to December 1981, the Dynamics Explorer-2 (DE-2) spacecraft successively observed the South polar and the North polar regions, and recorded the temperature, composition and dynamical structure of the upper thermosphere. In October 1981, perigee was about 310 km altitude, in the vicinity of the South Pole, with the satellite orbit in the 09.00–21.00 L.T. plane. During late November and December, the perigee had precessed to the region of the North Pole, with the spacecraft sampling the upper thermosphere in the 06.00 18.00 L.T. plane. DE-2 observed the meridional wind with a Fabry-Perot interferometer (FPI), the zonal wind with the wind and temperature spectrometer (WATS), the neutral temperature with the FPI, and the neutral atmosphere composition and density with the neutral atmosphere composition spectrometer (NACS). A comparison between the South (summer) Pole and the North (winter) Pole data shows considerable seasonal differences in all neutral atmosphere parameters. The region of the summer pole, under similar geomagnetic and solar activity conditions, and at a level of about 300 km, is about 300 K warmer than that of the winter pole, and the density of atomic oxygen is strongly depleted (and nitrogen enhanced) around the summer pole (compared with the winter pole). Only part of the differences in temperature and composition structure can be related to the seasonal variation of solar insolation, however, and both polar regions display structural variations (with latitude and Universal Time) which are unmistakeable characteristics of strong magnetospheric forcing. The magnitude of the neutral atmosphere perturbations in winds, temperature, density and composition within both summer and winter polar regions all increase with increasing levels of geomagnetic activity.The UCL 3-dimensional time dependent global model has been used to simulate the diurnal, seasonal and geomagnetic response of the neutral thermosphere, attempting to follow the major features of the solar and geomagnetic inputs to the thermosphere which were present during the late 1981 period.In the UCL model, geomagnetic forcing is characterized by semi-empirical models of the polar electric field which show a dependence on the Y component of the Interplanetary Magnetic Field, due to Heppner and Maynard (1983), It is possible to obtain an overall agreement, in both summer and winter hemispheres, with the thermospheric wind structure at high latitudes, and to explain the geomagnetic control of the combined thermal and compositional structure both qualitatively and quantitatively. To obtain such agreement, however, it is essential to enhance the polar ionosphere as a consequence of magnetospheric particle precipitation, reflecting both widespread auroral (kilovolt) electrons, and “soft” cusp and polar cap sources. Geomagnetic forcing of the high latitude thermosphere cannot be explained purely by a polar convective electric field, and the thermal as well as ionising properties of these polar and auroral electron sources are crucial components of the total geomagnetic input.     

First Results from SOHO

SOHO, the Solar and Heliospheric Observatory, is a project of international cooperation between ESA and NASA to study the Sun, from its deep core to the outer corona, and the solar wind. Three helioseismology instruments are providing unique data for the study of the structure and dynamics of the solar interior, from the very deep core to the outermost layers of the convection zone. A set of five complementary remote sensing instruments, consisting of EUV, UV and visible light imagers, spectrographs and coronagraphs, give us our first comprehensive view of the outer solar atmosphere and corona, leading to a better understanding of the enigmatic coronal heating and solar wind acceleration processes. Finally, three experiments complement the remote sensing observations by making in- situ measurements of the composition and energy of the solar wind and charged energetic particles, and another instrument maps the neutral hydrogen in the heliosphere and its dynamic change by the solar wind. This paper reports some of the first results from the SOHO mission. This revised version was published online in July 2006 with corrections to the Cover Date.     

天文地球动力学中的相对论效应

天文地球动力学利用空间与地面观测手段 ,监测和研究地球整体与各圈层的物质运动以及它们间的相互作用 ,这都离不开广义相对论涉及的时间与空间。随着空间对地观测精度的提高 ,为了充分利用高精度观测提供的信息 ,在天文地球动力学的研究中必须考虑相对论效应。所涉及的相对论效应包括 :( 1 )相对论参考系的建立 ,( 2 )在恰当的参考系中对观测者和被观测对象的相对论运动方程 (平动和自转 )的描述 ,( 3 )观测者和被观测对象间的电磁信号传播 ,( 4 )依赖于坐标选择的结果与具有物理意义的可观测量间的转换 ,( 5)某些基本概念与定义在广义相对论框架下的重新确认。本文对天文地球动力学中的这些相对论效应作了简要的评述。     

时间频率量的特征及其对时频系统建设的影响

分析了物理量测量中时间频率量的特点，主要有：时间的流逝性；其基准是自然基准；时间和频率既密切相关又有区别；时间频率具有最高的测量精细度(分辨率)与准确度；其计量标准可通过电磁波发播；其测量精确度与测量时间有关。另外，从基准、守时、授时、时间频率设备的研制、生产和队伍建设等方面阐发了这些特点对时间频率系统建设的影响。     

Mars geodesy,rotation and gravity

This review provides explanations of how geodesy, rotation and gravity can be addressed using radioscience data of an orbiter around a planet or of the lander on its surface.The planet Mars is the center of the discussion.The information one can get from orbitography and radioscience in general concerns the global static gravitational field, the time variation of the gravitational field induced by mass exchange between the atmosphere and the ice caps, the time variation of the gravitational field induced by the tides, the secular changes in the spacecraft's orbit induced by the little moons of Mars named Phobos and Deimos, the gravity induced by particular targets, the Martian ephemerides, and Mars' rotation and orientation.The paper addresses as well the determination of the geophysical parameters of Mars and, in particular, the state of Mars' core and its size, which is important for understanding the planet's evolution.Indeed, the state and dimension of the core determined from the moment of inertia and nutation depend in turn on the percentage of light elements in the core as well as on the core temperature, which is related to heat transport in the mantle.For example, the radius of the core has implications for possible mantle convection scenarios and, in particular, for the presence of a perovskite phase transition at the bottom of the mantle.This is also important for our understanding of the large volcanic province Tharsis on the surface of Mars.     

RARE EARTH AND OTHER ELEMENTS IN COMPONENTS OF THE ABEE ENSTATITE CHONDRITE

Abee clast samples, a matrix sample, a dark inclusion, magnetic and nonmagnetic samples, and bulk samples were analyzed by neutron activation analysis (NAA). The REE were determined by radiochemical NAA. Na, K, Sc, Cr, Mn, Fe, Co, Ni, Zn, As, Se, Sm, Ir, Au were determined by instrumental NAA. High abundances of As, Ir, and Au in the magnetic separate and the correlation of their abundances with the metal content of the clasts indicate that As, Ir, and Au chiefly occur in the metal. Correlations for Zn and Sc indicate that they chiefly occur in niningerite, but a significant amount of Sc may also occur in oldhamite. The dark inclusions do not follow the As and Zn correlations, suggesting that the dark inclusions and clasts are not equilibrated with each other. Correlation of the REE and oldhamite abundances for both the clasts and dark inclusions indicates that the REE chiefly occur in oldhamite. In view of the INAA and mineralogical evidence for non-equilibration among the clasts and a dark inclusion (Sears et al., 1981), the similar REE patterns for clasts (3,3) and dark inclusion (5,1), and the similar mineral composition of oldhamite in clast (3,3) and dark inclusion (5,1), suggest that the oldhamite in the clasts and dark inclusions is of a common origin, which Sears et al. (1981) showed could be formed by condensation. A Tb anomaly of a factor of 2 was found in sample (2,9 and 9,2), and a La anomaly of a factor of 2 was found in clast (3,3). The only other REE anomaly in Abee, a factor of 3.5 for Yb, was found by Nakamura and Masuda(1973). In view of the evidence for equilibration among the clasts, this anomaly must have been introduced shortly before the brecciation process and indicates that no significant reheating has occurred. This concurs with the findings of Sugiura and Strangeway (1981) and Bogard et al. (1982).