On Deriving Plasma Velocity Information from CDS/NIS Observations: Application to the Dynamics of Blinkers

Using standard instrument software and two independently developed data reduction and analysis procedures, we re-examine the accuracy of plasma velocity information derived from data obtained by the Solar and Heliospheric Observatory (SOHO)–Coronal Diagnostic Spectrometer (CDS). We discuss only the Ov 629 Å line data obtained by the Normal Incidence Spectrometer (NIS) and analyse a quiet Sun (QS) and active region (AR) dataset. Using the QS data, we demonstrate that the well-known North-South tilt in wavelength along the NIS slit varies significantly with time, which is not accounted for in the standard CDS correction procedures. In addition, when residual N – S trends exist in the data after processing, they may not be detected, nor removed, using the standard analysis software. This underscores the need for careful analysis of velocity results for individual datasets when using standard correction procedures. Furthermore, even when the results obtained by the two independent methods are well correlated (coefficients greater than 0.9), discrepancies in the values of the derived Doppler velocities can remain (95% within ±5 km s^?1). Therefore, we apply the results to examine the velocities obtained for EUV blinkers by previous authors. It is found that a strong correlation exists in the patterns of variation of the blinker velocities (> 0.98), even though there may be differences in their magnitudes. That is, in a clear majority of cases, the methods agree that a blinker is red-shifted or blue-shifted, although the uncertainty in the absolute velocity may be large.     

On helioseismic tests of basic physics

An important goal of helioseismology is to provide information about the basic physics and parameters that determine the structure of the solar interior. Here we discuss the procedures applied in such analyses, using as an example attempts to obtain significant constraints on the value of Newton's gravitational constant G from helioseismology. The analysis is based on complete direct and inverse helioseismic analysis of a set of accurate observed acoustic frequencies. We confirm, as found by previous investigations based on different approaches, that the actual level of precision of the helioseismic inferences does not allow us to constrain G with a precision better than that which can be reached with direct experimental measurements. The conclusion emphasizes the importance in helioseismic inferences of considering not only the accuracy with which solar oscillations are measured, but also the effect of uncertainties in other aspects of the model computation and helioseismic analysis.     

Ionospheric refraction in radio source observations at long radio wavelengths

Ionospheric refraction effects encountered in radio source observations in the 30 to 75 MHz range with the Clark Lake TPT telescope are discussed. It is found that simple calibration procedures are sufficient to provide positions of unknown sources with an accuracy of approximately one arcmin. Observations made near sunrise, or during disturbed ionospheric conditions must be discarded. If no corrections are applied, RMS errors of a few arcmin are to be expected. On leave at Radiosterrenwacht Dwingeloo, The Netherlands     

An accuracy estimation of the World CCD asteroid observations in the years 1999-2005

The MPC database of the asteroid observations (each position from near 20 millions) was used in analysis of observational accuracy for more than 300 active world observatories both professional and amateur. The values of the “Mean error of a single observation” σ (for α,δ) were derived based on the Pulkovo method of accuracy estimation. These values may be used for observatory weight assignment in the orbital improvement procedures. The accuracy of the best amateur observations is proved to be comparable with professional one (σ=±0^″.20). The detailed results in electronic format are accessible from the first author.     

Estimation of unmodeled forces on a lunar satellite

In previous investigations, a procedure for sequentially estimating the state of a lunar orbiting space vehicle acted upon by unmodeled terms in the lunar potential has been developed. Results obtained by processing tracking data from the Apollo 10 and 11 missions indicate that the algorithm provides more precise estimates of the vehicle state than conventional orbit determination procedures and, hence, provides an accurate input for navigation purposes. The question of the agreement of the estimates with the actual unmodeled accelerations has not been established. This investigation considers the question of the accuracy with which the algorithm can estimate the acceleration due to unmodeled lunar surface mascons. It is shown that an accurate estimate of the time history of the unmodeled acceleration can be obtained. The investigation also considers the effects of the magnitude and location of the mascons, as well as the effect of the observation accuracy.     

射电天文及卫星测地中大气折射的改正及其精度

本文综述了电离层、对流层中电波折射引起的射电天文观测及卫星测地中的各种误差及各种改正方法和它们的精度。 对流层影响的主要改正方法是实测大气温度、压力等参数,用数学模型计算。电离层影响的改正目前有三种方法:一是实时测量电离层主要参数——电子总含量的变化,然后用数学模型方法改正。二是采用双频同时观测的手段来消除电离层折射的影响。三是采用自校准方法。文中还比较了两种不同的自校准方法——常规自校准方法和多余量自校准方法。     

Accurate computation of highly eccentric satellite orbits

For computing highly eccentric (e0.9) Earth satellite orbits with special perturbation methods, a comparison is made between different schemes, namely the direct integration of the equations of motion in Cartesian coordinates, changes of the independent variable, use of a time element, stabilization and use of regular elements. A one-step and a multi-step integration are also compared.It is shown that stabilization and regularization procedures are very helpful for non or smoothly perturbed orbits. In practical cases for space research where all perturbations are considered, these procedures are no longer so efficient. The recommended method in these cases is a multi-step integration of the Cartesian coordinates with a change of the independent variable defining an analytical step size regulation. However, the use of a time element and a stabilization procedure for the equations of motion improves the accuracy, except when a small step size is chosen.     

Isotropic and anisotropic pointing models

This paper describes an alternative approach for generating pointing models for telescopes equipped with serial kinematics, esp. equatorial or alt-az mounts. Our model construction does not exploit any assumption for the underlying physical constraints of the mount, however, one can assign various effects to the respective components of the equations. In order to recover the pointing model parameters, classical linear least squares fitting procedures can be applied. This parameterization also lacks any kind of parametric singularity. We demonstrate the efficiency of this type of model on real measurements with meter-class telescopes where the results provide a root mean square accuracy of 1.5?2 arcseconds.     

The speckle polarimeter of the 2.5-m telescope: Design and calibration

The speckle polarimeter is a facility instrument of the 2.5-mSAIMSU telescope that combines the features of a speckle interferometer and a polarimeter. The speckle polarimeter is designed for observations in several visible bands in the following modes: speckle interferometry, polarimetry, speckle polarimetry, and polaroastrometry. In this paper we describe the instrument design and the procedures for determining the angular scale of the camera and the position angle of the camera and the polarimeter. Our measurements of the parameters for the binary star HD 9165 are used as an example to demonstrate the technique of speckle interferometry. For bright objects the accuracy of astrometry is limited by the error of the correction for the distortion caused by the atmospheric dispersion compensator. At zenith distances less than 45? the additional relative measurement error of the separation is 0.7%, while the additional error of the position angle is 0.3°. In the absence of a dispersion compensator the accuracy of astrometry is limited by the uncertainty in the scale and position angle of the camera, which are 0.15% and 0.06°, respectively. We have performed polarimetric measurements of unpolarized stars and polarization standards. The instrumental polarization at the Cassegrain focus in the V band does not exceed 0.01%. The instrumental polarization for the Nasmyth focus varies between 2 and 4% within the visible range; we have constructed its model and give a method for its elimination from the measurements. For stars with an intrinsic polarization of less than 0.2% during observations at the Cassegrain focus the error is determined mainly by the photon and readout noises and can reach 5 × 10^?5.     

Two Novel Approaches for Photometric Redshift Estimation based on SDSS and 2MASS

We investigate two training-set methods: support vector machines (SVMs) and Kernel Regression (KR) for photometric redshift estimation with the data from the databases of Sloan Digital Sky Survey Data Release 5 and Two Micron All Sky Survey. We probe the performances of SVMs and KR for different input patterns. Our experiments show that with more parameters considered, the accuracy does not always increase, and only when appropriate parameters are chosen, the accuracy can improve. For different approaches, the best input pattern is different. With different parameters as input, the optimal bandwidth is dissimilar for KR. The rms errors of photometric redshifts based on SVM and KR methods are less than 0.03 and 0.02, respectively. Strengths and weaknesses of the two approaches are summarized. Compared to other methods of estimating photometric redshifts, they show their superiorities, especially KR, in terms of accuracy.     

A New Analysis of Fireball Data from the Meteorite Observation and Recovery Project (MORP)

Sixty fireball cameras operated in Western Canada from 1971 to 1985. Over one thousand (1016) fireballs were recorded at more than one station, but only 367 were reduced, of which 285 have been published, including that of the Innisfree meteorite. Digitization of all the data is underway, and procedures are being developed which will allow the automatic reduction of events not previously examined. The results of the analysis of 80 fireballs reduced but not previously published are presented. When the new analysis is complete, the MORP archive will be a valuable source of information on meteoroid orbits.     

Large-Scale Transport of Magnetic Flux on the sun

The structure of the large-scale background magnetic field evolves in time and space. The large-scale horizontal transport velocity field of the magnetic flux patterns was inferred over the whole solar photosphere in the course of two solar activity cycles from year 1976 to 1999. The method of velocity determination and the testing procedures of the velocity accuracy are presented. The non-axially symmetric component of the horizontal velocity was found and both zonal and meridional velocity regions were described. The horizontal large-scale transport velocity regions vary in shape and the intensity during different phases of the 11-year solar activity cycle. The total horizontal transport velocity is characterized by the presence of variable amounts of the vector field vortices with symmetric orientation relative to the solar equator. The zonal velocity regions, distributed inside of the zonal belt limited by latitudes ± 35°, are persistent for about 4 Carrington rotations. Recurrent structures of similar velocity distributions are not coherent over the whole solar photosphere.     

Mapping monthly quantiles of ionospheric electron content and foF2

Monthly quantiles of ionization parameters are important information for system planning and to assess one aspect of ionospheric variability. For ionosonde parameters the data base is probably large enough to allow mapping of deciles. For electron content the situation is not so good: at least for Europe the data base only allows to construct maps for the upper and and the lower monthly quartiles.We present mapping procedures and maps for monthly quartiles of electron content and of foF2 based on experiences with mapping procedures for monthly medians. It turned out that for the lower quartiles it is necessary to to use more coefficients for the quartile maps than for the median maps.We discuss possibilities to gain decile information for electron content from foF2. For a few selected months and locations deciles can be gained from the Faraday effect on the signals of geostationary satellites. These selected cases are used to check on procedures to use foF2 experience to construct electron content deciles.     

Magnetofrictional Extrapolations of Low and Lou’s Force-Free Equilibria

We present a careful investigation of the magnetofrictional relaxation and extrapolation technique applied to the reconstruction of two test fields. These fields are taken from the family of nonlinear force-free magnetic equilibria constructed by Low and Lou (Astrophys. J. 352, 343, 1990), which have emerged as standard tests for extrapolation techniques in recent years. For the practically relevant case that only the field values in the bottom plane of the considered volume (vector magnetogram) are used as input information (i.e., not including the knowledge about the test field at the side and top boundaries), the test field is reconstructed to a higher accuracy than obtained previously. Detailed diagnostics of the reconstruction accuracy show that the implementation of fourth-order spatial discretization was essential to reach this accuracy for the given test fields and to achieve near machine precision in satisfying the solenoidal condition. Different variants of boundary conditions are tested, which all yield comparable accuracy. In its present implementation, the technique yields a scaling of computing time with total number of grid points only slightly below N ^5/3, which is too steep for applications to large (≥1024²) magnetograms, except on supercomputers. Directions for improvement are outlined.     

Sequential filter design for precision orbit determination and physical constant refinement

Earth-based spacecraft tracking data have historically been processed with classical least squares filtering techniques both for navigation purposes and for physical constant determination. The small, stochastic non-gravitational forces acting on the spacecraft are described to motivate the use of sequential estimation as an alternative to the least squares fitting procedures. The stochastic forces are investigated both in terms of their effect on the tracking data and their influence on estimation accuracy. A flexiible sequential filter design which leaves the existing trajectory, variational equations, data observable and partial computations undisturbed is described. A detailed filter design is presented that meets the precision demands and flexibility requirements of deep space navigation and scientific problems, one which provides a high degree of numerical integrity and numerical analysis capability, facilitates the efficient computation of multiple solutions, and makes few demands on the supporting computational structure.This paper presents the results of one phase of research carried out in part at the Jet Propulsion Laboratory, California Institute of Technology, under Contract No. NAS 7-100, sponsored by NASA.     

Can Dopplergrams be corrected for straylight?

Straylight effects in solar intensity images have been framed in a well-established mathematical formulation (Chae et al., 1998a, b). Many procedures can be applied to recover the point spread function (PSF) and correct the atmospheric and instrumental distortions. Concerning velocity data, some procedures to make corrections for straylight have been addressed. Many experiments infer the velocity as the mean displacement of a spectral line and build a Dopplergram as a nonlinear combination of the intensity images sampled at different wavelengths. Unfortunately, often only the final products are available and it is not possible to apply an appropriate correction to the individual intensity images. A standard procedure to correct for the straylight has been applied to a time-series of velocity data and the results have been controlled in the single Dopplergram (that seems to be good) and in a power spectrum at the atmospheric frequencies (where any correction should reduce the power). These frequencies are characteristic of the observing site and typically in the range of the g modes: this means that no long run from a network would decrease their contribution to the noise. The results show that the correction procedures, based on the intensity and velocity gradients, are not reliable if controlled in the power spectrum. This conclusion has led to a change of the acquisition procedure at the Kanzelhöhe Solar Observatory, where now the individual intensity images that build the Dopplergram are stored.     

中高轨卫星广播星历精度分析

GPS广播星历参数具有物理意义明确、参数少、精度高等优点,可以考虑将它应用于其他卫星导航系统。但是GPS系统的卫星构成比较单一,而其他卫星导航系统可能包含中地球轨道 (MEO)、倾斜地球同步轨道(IGSO)和地球静止轨道(GEO)等多种不同类型的中高轨卫星。分析了采用GPS广播星历参数时,MEO、IGSO和GEO卫星的广播星历拟合精度,特别讨论了轨道倾角接近于0的GEO卫星的广播星历拟合精度,并给出了相应的改进措施。计算表明,对于 MEO卫星,2 h的广播星历拟合精度(三维位置)可达厘米级;对于IGSO卫星和轨道倾角较大的GEO卫星,4 h的广播星历拟合精度约为0．1 m,径向位置误差在厘米量级;而对于轨道倾角接近于0的GEO卫星,若不采取特殊措施,由于轨道倾角和升交点经度统计相关,其广播星历拟合精度很差,为此提出了一种坐标转换方法。采用此方法后的广播星历拟合精度可达0．1 m,径向位置误差为厘米量级。     

稳定化方法与后稳定化方法的比较

对BaumgaLrte的稳定化和Chin的后稳定化进行了详尽讨论与数值比较．用经典数值方法并结合这两种稳定化方式都能提高数值精度和改善数值稳定性．在最佳稳定参数下稳定化精度一般不等价于后稳定化．两者精度优劣并无常定．考虑到Baumgarte的稳定化使得数值积分的右函数更复杂和增加计算耗费,尤其是存在稳定参数最佳选取的麻烦,故推荐后稳定化投入实算．但值得注意的是用后稳定化与没有经过稳定化处理的经典积分器来比不宜扩大积分步长．     

Photometric redshifts of galaxies from SDSS and 2MASS

In order to find the physical parameters which determine the accuracy of pho- tometric redshifts, we compare the spectroscopic and photometric redshifts (photo-z's) for a large sample of ～ 80 000 SDSS-2MASS galaxies. Photo-z's in this paper are es- timated by using the artificial neural network photometric redshift method (ANNz). For a subset of～40000 randomly selected galaxies, we find that the photometric redshift recovers the spectroscopic redshifi distribution very well with rms of 0.016. Our main results are as follows: (1) Using magnitudes directly as input parameters produces more accurate photo-z's than using colors; (2) The inclusion of 2MASS (3, H, Ks) bands does not improve photo-z's significantly, which indicates that near infrared data might not be important for the low-redshift sample; (3) Adding the concentration index (essentially the steepness of the galaxy brightness profile) as an extra input can improve the photo-z's estimation up to～10 percent; (4) Dividing the sample into early- and late-type galaxies by using the concentration index, normal and abnormal galaxies by using the emission line flux ratios, and red and blue galaxies by using color index (g - r), we can improve the accuracy of photo-z's significantly; (5) Our analysis shows that the outliers (where there is a big difference between the spectroscopic and photometric redshifts) are mainly correlated with galaxy types, e.g., most outliers are late-type (blue) galaxies.     

卫星双向时间比对及其误差分析

介绍了卫星双向时间比对(TwsTT)方法的发展过程及其基本原理，讨论了TwsTT的误差源，分析了各误差源对比对精度的影响。得出影响TwSTT精度的主要误差源为设备时延误差和路径传播时延误差，而由于卫星和地面站运动引起的误差相对较小，从而给出了实际应用时不同比对精度下需要考虑的误差源及对误差源的精度要求。