Excitation of Chandler Wobble by Pacific,Indian and Atlantic Oceans from 1980 to 2005

In addition to the atmosphere, the oceans play important roles in the excitation of the Chandler wobble. The contributions made by the Pacific Ocean, the Indian Ocean and the Atlantic Ocean from 1980 to 2005 to the excitation of the Chandler wobble are first and systematically researched by taking advantage of the data of the current velocity field and ocean floor pressure provided by the marine circulation model of the Estimating the Circulation and Climate of the Ocean (ECCO). Studies show that the contributions of the three oceans to the excitation of the Chandler wobble are different from one another: the excitation energy of the Pacific Ocean makes up about 22.2% of the observational excitation energy, the largest one among the three oceans, that of the Indian Ocean accounts for about 12.7% and that of the Atlantic Ocean amounts to about 7.1%, the smallest among the three great oceans. The remarkable increase in the excitation energy of the Chandler wobble by the Pacific Ocean may be possibly due to the effect of the strong ENSO event which occurred from 1982 to 1983.     

光学系统参数确定带宽约束的天文图像盲反卷积

大气湍流极大限制了地基大口径望远镜观测天文目标图像的空间分辨率.根据最大似然估计原理,提出了用实际光学带宽约束的可有效减小天文观测图像中大气湍流影响的盲反卷积方法,通过共轭梯度优化算法使卷积误差函数趋向最小.建立了望远镜光学系统参数和图像频域带宽的关系,采用变量正性约束、点扩散函数带宽有限约束,提高算法的收敛性.为避免图像处理中有效傅立叶变换频率超出截止频率,要求采集望远镜焦面图像时单个成像单元(如CCD像素单元)应小于四分之一衍射斑直径.算法中未用目标支持域约束,所提出的方法适用于全视场天文图像恢复.用计算机模拟和对实际天文目标双鱼座图像数据的恢复结果验证了所提出方法的有效性.     

一般回路形式下的暗条电流平衡高度

Ｋｉｐｐｅｎｈａｈｎ－Ｓｃｈｌｕｔｅｒ模型回路和Ｋｕｐｐｅｒｕｓ－Ｒａａｄｕ模型回路是暗条电流回路的两种极端形式，本文则从一般回路形式来研究暗条电流与暗条高度的关系。结果发现：当暗条电流与活动区磁场作用的Ｌｏｒｅｎｔｚ力方向与暗条重力方向相同时，流入活动区光球层的电流越大，其在暗条电流中所占的比例越大，暗条越稳定。当二力方向相反时，流入活动区光球层的电流越大，其在暗条电流中所占的比例越大，暗条越不稳定。     

Spectral and thermodynamic constraints on the existence of gypsum at the Juventae Chasma on Mars

The spectral imaging of the Mars obtained with the Mars Express/OMEGA experiment demonstrates that a majority of the sulfates-rich regions are associated with the interior light-toned layered deposits within the canyon system in the equatorial zone of the planet. While all sulfates-rich deposits inside the canyons are characterized by the presence of the kieserite and hydrated magnesium sulfates, the spectral features of gypsum were detected only in the Juventae Chasma and the Iani Chaos. The detection of gypsum in the upper part of the layered deposits, stacking the erosional remnant on the floor of the Juventae Chasma (above the spectral signature of the kieserite and polyhydrated sulfates detected on the flanks of the remnant) represents a more intriguing case. To clarify the question of the presence of gypsum in the Juventae Chasma, we present reanalyzed OMEGA spectra within that area and performed the chemical equilibrium modelling of sulfates precipitation sequence at the freezing and the evaporation of a hypothetical aqueous solution which could have existed within the Chasma in the past. Our results did not confirm the presence of distinct spectral signatures of gypsum. The results of equilibrium modelling also exclude significant precipitation of gypsum during the latest stage of the aqueous sedimentation, responsible for the formation of the upper part of the erosional remnant.     

Circular binary star and an interstellar interloper: The analytical solution

We consider an interstellar interloper moving at a relatively large distance from a circular binary star. We use the analytical method of separating rapid and slow subsystems, the rapid subsystem being the binary and the slow subsystem being the interstellar interloper. We show that due to the higher than geometrical symmetry of the problem, in addition to the conservation of the energy and the projection of the angular momentum on the axis of the rotation of the binary, the square of the angular momentum is also conserved. In the course of the time evolution, the vector of the angular momentum rotates about that axis at the constant angle to the axis. After obtaining this general counterintuitive result, we focus at the case where the interstellar interloper is coplanar with the binary. We provide an explicit equation of the motion of the interloper. Then we calculate analytically the angle of deflection of the interloper from the straight line. We analyze the difference in the angle of deflection between this three-body problem and the corresponding two-body problem: we show that this difference remains almost constant (a negative constant) at the range of the eccentricities of the interloper trajectory relatively close to unity and linearly increases (by the absolute value, remaining negative) with the eccentricity as the latter becomes much greater than unity.     

原子钟时域频率稳定度分析中等效自由度的计算与分析

该文基于标准化自协方差函数的组合算法计算了阿仑方差、哈达玛方差的重叠和非重叠形式以及改进阿仑方差、时间方差的等效自由度,并用经验公式计算了总方法估计的等效自由度,在此基础上对各种方差估计不同噪声情况下的等效自由度进行了比较分析,得出如下结论,总方法估计的等效自由度大于相应非总方法估计的等效自由度,重叠和非重叠阿仑方差估计的等效自由度大于相应哈达玛方差估计的等效自由度;对非重叠估计而言,噪声的频率越低,其等效自由度越大;而相应的重叠估计和总方法估计与之相反,噪声的频率越高,其等效自由度越大,这也说明重叠估计和相应的总方法估计对高频噪声等效自由度的提高优于低频噪声.     

Double Plasma Resonance for Decimetric dot-Emissions

The mechanism of the double plasma resonance has been applied to the dot-like emission structures observed in the decimetric band. In particular, the existence of dot-emissions on definite trajectories in the dynamic spectra offers a rare opportunity to delineate the details of the double resonance mechanism as well as the diagnostics of the radio source and the ambient corona. The modeling of the chains of dots points to the excitation of the electron Bernstein modes predominantly at the second cyclotron harmonic and their subsequent coalescence to produce the observed dot-emissions. The coalescence process doubles the cyclotron harmonic number and thereby reduces the cyclotron absorption significantly. The discrete nature of the emissions results from the finite extent of the resonant region due to the presence of inhomogeneities. The determination of the relative scales of spatial variation of the magnetic field and the electron density in the region of the decimetric emission is particularly desirable as it has important consequences for the energy release, acceleration processes and the configuration of the participating plasma structures.     

Distributions of tangential discontinuity in the corotating solar wind structure

Distributions of the tangential discontinuity (TD) in the solar wind sector structure are investigated on the basis of the magnetic field data and the ion plasma parameters from the Explorer 33 satellite from 23 January to 23 March 1968. The TD is separated from the observed field fluctuations by calculating the direction of the plasma flow and also the direction of the minimum field fluctuation with respect to the ambient magnetic field direction.It is found that the TD is formed by the thin layered field-aligned currents (the current sheets), and that the TD is predominantly built up in the leading edge of the solar wind where the compression of the plasma and the magnetic field takes place.It is suggested that the current sheets might be locally generated in the leading edge in the turbulent conditions arising from collisions between the fast- and the slow-stream of the solar sector structure.     

The Oscillations of Coronal Loops Including the Shell

We investigate the MHD waves in a double magnetic flux tube embedded in a uniform external magnetic field. The tube consists of a dense hot cylindrical cord surrounded by a co-axial shell. The plasma and the magnetic field are taken to be uniform inside the cord and also inside the shell. Two slow and two fast magnetosonic modes can exist in the thin double tube. The first slow mode is trapped by the cord, the other is trapped by the shell. The oscillations of the second mode have opposite phases inside the cord and shell. The speeds of the slow modes propagating along the tube are close to the tube speeds inside the cord and the shell. The behavior of the fast modes depends on the magnitude of Alfvén speed inside the shell. If it is less than the Alfvén speed inside the cord and in the environment, then the fast mode is trapped by the shell and the other may be trapped under the certain conditions. In the opposite case when the Alfvén speed in the shell is greater than those inside the cord and in the environment, then the fast mode is radiated by the tube and the other may also be radiated under certain conditions. The oscillation of the cord and the shell with opposite phases is the distinctive feature of the process. The proposed model allows to explain the basic phenomena connected to the coronal oscillations: i) the damping of oscillations stipulated in the double tube model by the radiative loss, ii) the presence of two different modes of perturbations propagating along the loop with close speeds, iii) the opposite phases of oscillations of modulated radio emission, coming from the near coronal sources having sharply different densities.     

On the Formation of Planetesimals: Radial Contraction of the Dust Layer Interacting with the Protoplanetary Disk Gas

Radial contraction of the dust layer in the midplane of a gas–dust protoplanetary disk that consists of large dust aggregates is modeled. Sizes of aggregates vary from centimeters to meters assuming the monodispersion of the layer. The highly nonlinear continuity equation for the solid phase of the dust layer is solved numerically. The purpose of the study is to identify the conditions under which the solid matter is accumulated in the layer, which contributes to the formation of planetesimals as a result of gravitational instability of the dust phase of the layer. We consider the collective interaction of the layer with the surrounding gas of the protoplanetary disk: shear stresses act on the gas in the dust layer that has a higher orbital velocity than the gas outside the layer, this leads to a loss of angular momentum and a radial drift of the layer. The stress magnitude is determined by the turbulent viscosity, which is represented as the sum of the α-viscosity associated with global turbulence in the disk and the viscosity associated with turbulence that is localized in a thin equatorial region comprising the dust layer and is caused by the Kelvin–Helmholtz instability. The evaporation of water ice and the continuity of the mass flux of the nonvolatile component on the ice line is also taken into account. It is shown that the accumulation of solid matter on either side of the ice line and in other regions of the disk is determined primarily by the ratio of the radii of dust aggregates on either side of the ice line. If after the ice evaporation the sizes (or density) of dust aggregates decrease by an order of magnitude or more, the density of the solid phase of the layer’s matter in the annular zone adjacent to the ice line from the inside increases sharply. If, however, the sizes of the aggregates on the inner side of the ice line are only a few times smaller than behind the ice line, then in the same zone there is a deficit of mass at the place of the modern asteroid belt. We have obtained constraints on the parameters at which the layer compaction is possible: the global turbulence viscosity parameter (α < 10^?5), the initial radial distribution of the surface density of the dust layer, and the distribution of the gas surface density in the disk. Restrictions on the surface density depend on the size of dust aggregates. It is shown that the timescale of radial contraction of a dust layer consisting of meter-sized bodies is two orders of magnitude and that of decimeter ones, an order of magnitude greater than the timescale of the radial drift of individual particles if there is no dust layer.     

The evolution of loop structures in flux rings within the solar convection zone

Choudhuri and Gilman (1987) considered certain implications of the hypothesis that the magnetic flux within the Sun is generated at the bottom of the convection zone and then rises through it. Taking flux rings symmetric around the rotation axis and using reasonable values of different parameters, they found that the Coriolis force deflects these flux rings into trajectories parallel to the rotation axis so that they emerge at rather high latitudes. This paper looks into the question of whether the action of the Coriolis force is subdued when the initial configuration of the flux ring has non-axisymmetries in the form of loop structures. The results depend dramatically on whether the flux ring with the loops lies completely within the convection zone or whether the lower parts of it are embedded in the stable layers underneath the convection zone. In the first case, the Coriolis force supresses the non-axisymmetric perturbations so that the flux ring tends to remain symmetric and the trajectories are very similar to those of Choudhuri and Gilman (1987). In the second case, however, the lower parts of the flux ring may remain anchored underneath the bottom of the convection zone, but the upper parts of the loops still tend to move parallel to the rotation axis and emerge at high latitudes. Thus the problem of the magnetic flux not being able to come out at the sunspot latitudes still persists after the non-axisymmetries in the flux rings are taken into account.National Center for Atmospheric Research is sponsored by the National Science Foundation.     

On the accuracy of the orbit of asteroid (99942) APOPHIS at the time of its encounter with the Earth in 2029

We estimate the effect of trajectory measurement errors on the orbital parameters of asteroid Apophis determined from improvements. For this purpose, based on all of the optical and radar observations available to date, we have computed a nominal orbit of the asteroid. The scatter ellipsoid of the initial conditions of motion has been obtained by two methods. In the first, universally accepted method, the scatter ellipsoid is calculated by assuming a linear dependence of the errors in the parameters being determined on observational errors. In the second method, the scatter region of the orbital parameters around the nominal-orbit parameters is determined by the Monte Carlo method. We show that the region determined by the latter method at the initial epoch differs only slightly from the scatter ellipsoid for the linear approximation. We estimate the sizes of the projections of the corresponding regions onto the target plane at the time of the closest encounter of the asteroid with the Earth in 2029. The projections are approximated by ellipses. Our computations have shown that the ellipse has the following sizes: 389.6 km for the semimajor axis and 16.4 km for the semiminor axis in the linear case and 330.0 and 11.1 km, respectively, in the nonlinear case.     

Shock jump relations for a dusty gas atmosphere

This paper presents simplified forms of jump relations for one dimensional shock waves propagating in a dusty gas. The dusty gas is assumed to be a mixture of a perfect gas and spherically small solid particles, in which solid particles are continuously distributed. The simplified jump relations for the pressure, the temperature, the density, the velocity of the mixture and the speed of sound have been derived in terms of the upstream Mach number. The expressions for the adiabatic compressibility of the mixture and the change-in-entropy across the shock front have also been derived in terms of the upstream Mach number. Further, the handy forms of shock jump relations have been obtained in terms of the initial volume fraction of small solid particles and the ratio of specific heats of the mixture, simultaneously for the two cases viz., (i) when the shock is weak and, (ii) when it is strong. The simplified shock jump relations reduce to the Rankine-Hugoniot conditions for shock waves in an ideal gas when the mass fraction (concentration) of solid particles in the mixture becomes zero. Finally, the effects due to the mass fraction of solid particles in the mixture, and the ratio of the density of solid particles to the initial density of the gas are studied on the pressure, the temperature, the density, the velocity of the mixture, the speed of sound, the adiabatic compressibility of the mixture and the change-in-entropy across the shock front. The results provided a clear picture of whether and how the presence of dust particles affects the flow field behind the shock front. The aim of this paper is to contribute to the understanding of how the shock waves behave in the gas-solid particle two-phase flows.     

Sediment budget and tectonic evolution of the Meuse catchment in the Ardennes and the Roer Valley Rift System

The Meuse river system is located in the northeastern part of the Paris Basin, the Ardennes, and the Roer Valley Rift System (RVRS). The Meuse river system developed during the uplift of the Ardennes since the Eocene and it was affected by renewed rifting of the RVRS starting in the Late Oligocene. In response to the uplift of the Ardennes, the river system incised and a terrace sequence developed during the Plio–Pleistocene. The sediments generated by erosion in the catchment were transported into the RVRS and further to the north, into the Zuiderzee Basin and the North Sea Basin. Using a digital terrain model, the amount of eroded rock volume versus time for the Meuse catchment has been computed using the Paleogene and older planation surfaces and the fluvial terraces. Comparison of the amount of eroded material with the volume of sediment preserved in the RVRS for the early Middle Pleistocene shows that about 17.5% of the sediment volume transported into the RVRS remained there, the rest being transported further into the Zuiderzee Basin and the North Sea Basin. The Quaternary tectonic uplift of the Ardennes inferred from the incision history of the Meuse river system is characterized by a long-term uplift, on which a Middle Pleistocene acceleration is superimposed. The accelerated uplift is contemporaneous with an uplift event in the RVRS and in the neighbouring Eifel area, and with the onset of the youngest phase of volcanism in the Eifel area. The areal distribution of this uplift is characterized by a dome shape centered around the Eifel area.     

Spectroscopic observations of the exoplanet WASP-32b transit

We present first results of spectroscopic observations of transiting exoplanets in the Special Astrophysical Observatory of the Russian Academy of Sciences with the Main Stellar Spectrograph of the 6-m BTA telescope. For the exoplanetWASP-32b, we detected a significant variation of intensity and equivalent width in the Hα spectral line of the parent star at the time of a transit. The equivalent width of the line during transit is by 8–10% larger than outside the planet passage. Residual intensity in the core of the line reveals the following tendency: the line is by 10–15% deeper inside transit than outside it. Observations with the long-slit spectrograph of the Crimean Astrophysical Observatory at the 2.6-m ZTSh telescope also showed a transit event in the Hα line, although, with a smaller amplitude and shape inverted in relation to the data from the 6-m telescope. While in the observations with the BTA the Hα line becomes deeper during the transit, in the ZTSh observations, the residual intensity of the Hα line decreases during the transit. Reducing and analysis of the archive data of WASP-32b observations with the HARPS spectrograph also confirm the Hα line modulation at the time of the transit. The observed data give evidence of the envelope in WASP-32b filling the Roche lobe and a comet-like tail of changing geometry and orientation relative to the observer. These changes determine different depths and shapes of the Hα spectral line at the time of transits.     

Modifications of the magnetospheric plasma due to ponderomotive forces

By means of two simple examples it is demonstrated that the ponderomotive force induced by the geomagnetic pulsations significantly modifies the plasma distribution in the Earth's magnetosphere. The first example considers the quasi-static plasma equilibrium at closed magnetic shells. Ponderomotive forces, averaged over the geomagnetic pulsation period, increase the plasma density at the equator of the oscillating magnetic shell. If the oscillation amplitude exceeds a certain critical level, the quasi-static solution predicts a maximum of the plasma density at the equator of the magnetic shell. The second example considers the quasi-stationary plasma flow along the open field lines stretching into the geomagnetic tail (the polar wind). The inclusion of the ponderomotive force may shift the critical point for the flow transition through the supersonic barrier closer to the Earth.     

On some problems concerning the calculation of the form of the magnetopause and the electric field on the magnetopause in the framework of the continuum medium model

In order to understand the reason of the existence of the electric field in the magnetosphere, and for the theoretical evaluation of its value, it is necessary to find the solution of the problem of determination of the magnetosphere boundary form in the frameworks of the continuum medium model which takes into account part of the magnetospheric plasma movement in supporting the magnetospheric boundary equilibrium. A number of problems for finding the distribution of the pressure, the density, the magnetic field and the electric field on the particular tangential discontinuity is considered in the case when the form of discontinuity is set (the direct problem) and a number of problems for finding the form of the discontinuity and the distribution of the above-mentioned physical quantities on the discontinuity is considered when the law of the change of the external pressure along the boundary is set (for example, with the help of the approximate Newton equation). The problem which is considered here, which deals with the calculation of the boundary form and with the calculation of the distribution of the corresponding physical quantities on the discontinuity of the 1st kind for the compressible fluid with the magnetic field with field lines which are perpendicular to the plane of the flow in question, concerns the last sort of problems. The comparison of the results of the calculation with the data in the equatorial cross-section of the magnetosphere demonstrates that the calculated form of the boundary, the value of the velocity of the return flow and the value of the electric field on the magnetopause, agree satisfactorily with the observational data.     

Grain formation in the ejecta of a supernova: the nucleation rate of grains heated by radiation

In the ejecta of a supernova, the temperature of the created grains differs from that of the gas due to the radiation from the star. We investigated the grain formation in the supernova using a developed new nucleation rate where the temperature difference between the gas and the grains is taken into account. If the temperature of the grains is higher than that of the gas, the nucleation process does not occur when the gaseous temperature attains the condensation temperature. As a result we found that the temperature difference between the gas and the grains in SN1987A is about 50–200K which leads that the nucleation is delayed for about 20–100 days.     

Simple Model of a Stochastically Excited Solar Dynamo

The aim of this paper is to investigate the dynamical nature of the complexity observed in the time evolution of the sunspot number. We report a detailed analysis of the sunspot number time series, and use the daily records to build the phase space of the underlying dynamical system. The observed features of the phase space prompted us to describe the global behavior of the solar cycle in terms of a noise-driven relaxation oscillator. We find the equations whose solutions best fit the observed series, which adequately describe the shape of the peaks and the oscillations of the system. The system of equations obtained from this fitting procedure is shown to be equivalent to a truncation of the dynamo equations. A linear transformation maps the phase space of these equations into the phase space reconstructed from the observations. The irregularities of the solar cycle were modeled through the introduction of a stochastic parameter in the equations to simulate the randomness arising in the process of eruption of magnetic flow to the solar surface. The mean values and deviations obtained for the periods, rise times and peak values, are in good agreement with the values obtained from the sunspot time series.