天文学研究发展趋势分析--恒星与恒星系统

对在1981-2000年世界上所发表的和中国学者所发表的有关恒星与恒星系统的论文作统计发现：此期间世界上这一领域的发展较平稳，而我国的发展快速．这反映了改革开放后，我国基础学科研究大有进展．从各分支所占的比重和发展来看，我国在恒星与恒星系统的研究与世界同期有几乎相同的分布，因此总体上我国在这一领域的发展基本正常．当然有些分支发展较快，如有关超新星及其遗迹、星际介质和恒星形成区、化学丰度的研究等，这和一些较强的研究团组形成有关；在双星研究方面，我国则与世界发展一致，双星研究始终是恒星研究领域的重点；而在世界范围内较突出的关于银河系的研究，在我国却相对较弱．恒星和恒星系统这一研究领域20年的论文数统计显示，我国学者所发表的论文只占世界总论文数的1．3％，虽然在最后5年有大幅上升，但也只占2．0％，这与我国IAU会员数所占比例相比是偏少的．就世界整体而言，恒星领域的研究进展与整个天文学领域相比是较慢的，显然这与一批能做深空探测和高能波段观测的设备投入有关．因此，除了对恒星及恒星系统领域作统计分析外，对整个天文学领域各大分支作分析可能对制定今后我国天文学发展计划更有利。     

Agricultural impacts on ecosystem functioning in temperate areas of North and South America

Land use has a large impact on ecosystem functioning, though evidences of these impacts at the regional scale are scarce. The objective of this paper was to analyze the impacts of agricultural land use on ecosystem functioning (radiation interception and carbon uptake) in temperate areas of North and South America. From land cover maps generated using high-resolution satellite images we selected sites dominated by row crops (RC), small grain crops (SG), pastures (PA), and rangelands (RA) in the Central Plains of USA and the Pampas of Argentina. These two regions share climatic characteristics and the agricultural conditions (crop types) are also very similar. Both areas were originally dominated by temperate grasslands. In these sites we extracted the temporal series of the normalized difference vegetation index (NDVI) from the NOAA satellites for the period 1989–1998 and calculated the mean seasonal NDVI curve for each site. Additionally, we calculated the mean annual NDVI, the maximum NDVI, the date of the year when the max NDVI was recorded and the interannual variability of these three attributes. We compared the mean values of each NDVI-derived attribute between land cover types and between continents. The NDVI seasonal patterns for each land cover type were roughly similar between the Central Plains and the Pampas during the growing season. The largest differences were observed during the winter and spring, when the NDVI of all land cover types in the Central Plains remained at lower values than in the Pampas. This was probably caused by the high annual thermal amplitude in the Central Plains that results in a much more restricted growing season. As a result of these differences in the shape of the NDVI curve, the mean annual NDVI in the Central Plains was lower than in the Pampas for all land cover types but the maximum NDVI did not differ importantly. In both regions, row crops delayed the date of the NDVI peak, small grain crops advanced it and pastures did not change it importantly, compared with rangelands. The interannual variability of the NDVI attributes was higher for small grains than for row crops in both regions. However, small grains crops were consistently more variable between years in the Central Plains than in the Pampas. The opposite occurred with pastures and rangelands, which were more variable in the Pampas than in the Central Plains. This paper confirms and generalizes previous findings that showed important imprints of land use on ecosystem functioning in temperate ecosystems. Our results support the idea that the changes in land cover that have occurred in the Central Plains and the Pampas leaded to similar changes in the way that ecosystems absorb solar radiation and in the patterns of carbon uptake.     

Unusual stripes in emission and absorption in solar radio bursts: Ropes of fibers in the meter wave band

Based on data from the spectrographs of IZMIRAN and Tremsdorf station (Astronomical Institute, Potsdam), we analyze the ropes of narrow-band fibers in the spectra of solar radio bursts in the meter wave band by invoking events of satellite data (SOHO/LASCO, EIT, MDI) for the analysis. We consider in detail basic properties of the ropes in four events in comparison with previously known data. The fibers in ropes are more commonly observed with an overlap in time and frequency, but occasionally (more often at the end of the ropes) they can follow with a separation in time. The fiber duration and recurrence period seldom remain stable and, in general, increase from 0.3–0.5 s at the beginning to several seconds at the end of the rope. The relative values of the instantaneous and total fiber frequency bandwidths change only slightly in different events; δ f / f ≈ 0.003–0.005 and Δf / f ≈ 0.02–0.03. Most of the ropes exhibit a low-frequency absorption. The fibers in ropes are similar to ordinary intermediate drift bursts (fiber bursts), but they drift in a narrow frequency band and have a more frequent recurrence in some events. The ropes of fibers are usually observed in the time interval when the shock front catches up with the leading edge of a coronal mass ejection. Under the condition of a unified approach to interpreting the ropes of fibers in all events, their basic properties can be explained in terms of the model of fiber bursts. The connection of fibers with the developed zebra pattern is shown within the framework of a unified approach to the formation theory of stripes in emission and absorption in the model on whistlers.     

Uncertainties in climate change prediction: El Niño-Southern Oscillation and monsoons

The sensitivity of climate phenomena in the low latitudes to enhanced greenhouse conditions is a scientific issue of high relevance to billions of people in the poorest countries of the globe. So far, most studies dealt with individual model results. In the present analysis, we refer to 79 coupled ocean–atmosphere simulations from 12 different climate models under 6 different IPCC scenarios. The basic question is as to what extent various state-of-the-art climate models agree in predicting changes in the main features of El Niño-Southern Oscillation (ENSO) and the monsoon climates in South Asia and West Africa. The individual model runs are compared with observational data in order to judge whether the spatio-temporal characteristics of ENSO are well reproduced. The model experiments can be grouped into multi-model ensembles. Thus, climate change signals in the classical index time series, in the principal components and in the time series of interannual variability can be evaluated against the background of internal variability and model uncertainty.There are large differences between the individual model predictions until the end of the 21st century, especially in terms of monsoon rainfall and the Southern Oscillation index (SOI). The majority of the models tends to project La Niña-like anomalies in the SOI and an intensification of the summer monsoon precipitation in India and West Africa. However, the response barely exceeds the level of natural variability and the systematic intermodel variations are larger than the impact of different IPCC scenarios. Nonetheless, there is one prominent climate change signal, which stands out from model variations and internal noise: All forced model experiments agree in predicting a substantial warming in the eastern tropical Pacific. This oceanic heating does not necessarily lead to a modification of ENSO towards more frequent El Niño and/or La Niña events. It simply represents a change in the background state of ENSO. Indeed, we did not find convincing multi-model evidence for a modification of the wavelet spectra in terms of ENSO or the monsoons. Some models suggest an intensification of the annual cycle but this signal is fairly model-dependent. Thus, large model uncertainty still exists with respect to the future behaviour of climate in the low latitudes. This has to be taken into account when addressing climate change signals in individual model experiments and ensembles.     

Impacts of deforestation and afforestation in the Mediterranean region as simulated by the MPI atmospheric GCM

For two reasons it is important to study the sensitivity of the global climate to changes in the vegetation cover over land. First, in the real world, changes in the vegetation cover may have regional and global implications. Second, in numerical simulations, the sensitivity of the simulated climate may depend on the specific parameterization schemes employed in the model and on the model's large-scale systematic errors. The Max-Planck-Institute's global general circulation model ECHAM4 has been used to study the sensitivity of the local and global climate during a full annual cycle to deforestation and afforestation in the Mediterranean region. The deforestation represents an extreme desertification scenario for this region. The changes in the afforestation experiment are based on the pattern of the vegetation cover 2000 years before present when the climate in the Mediterranean was more humid. The comparison of the deforestation integration to the control shows a slight cooling at the surface and reduced precipitation during the summer as a result of less evapotranspiration of plants and less evaporation from the assumption of eroded soils. There is no significant signal during the winter season due to the stronger influence of the mid-latitude baroclinic disturbances. In general, the results of the afforestation experiment are opposite to those of the deforestation case. A significant response was found in the vicinity of grid points where the land surface characteristics were modified. The response in the Sahara in the afforestation experiment is in agreement with the results from other general circulation model studies.     

The kinetic nature of turbulence at short scales in the solar wind

The analysis of the transition from the large-scale fluid regime to the short-scale kinetic range of wavelengths in the development of the turbulent cascade of energy is nowadays subject of fervent discussion in the space plasmas scientific community. We make use of Hybrid Vlasov-Maxwell simulations where the full kinetic dynamics of ions is taken into account, while electrons are treated as a fluid. We investigate the development of turbulence in the solar wind, in 1D-3V phase space configuration and in the frequency range across the ion cyclotron frequency. These simulations allow for the analysis of the role of kinetic effects in the short-scale region of the energy spectra in the direction parallel to the background magnetic field. Our numerical results show the presence of a significant electrostatic activity at small wavelengths, triggered by the resonant interaction of ions with longitudinal waves. Our model does not allow to take into account the evolution of the turbulent spectra in the plane perpendicular to the ambient field, due to limited dimensionality in phase space. On the other hand, this model permits to isolate and study the possibility of transferring the electromagnetic large-scale energy on the small-scale kinetic electrostatic component of the spectrum. Peculiar features observed in the spacecraft data in the solar wind are qualitatively reproduced within the hybrid-Vlasov model, such as the generation of perpendicular temperature anisotropy and accelerated longitudinal beams of ions in the distribution of particle velocities as well as the appearance of a marked peak of electrostatic activity in the short-scale termination of the turbulent spectra.     

Numerical modelling of permafrost in bedrock in northern Fennoscandia during the Holocene

The occurrence of permafrost in bedrock in northern Fennoscandia and its dependence on past and presently ongoing climatic variations was investigated with one- (1D) and two-dimensional (2D) numerical models by solving the transient heat conduction equation with latent heat effects included. The study area is characterized by discontinuous permafrost occurrences such as palsa mires and local mountain permafrost. The ground temperature changes during the Holocene were constructed using climatic proxy data. This variation was used as a forcing function at the ground surface in the calculations. Several versions of the present ground temperature were applied, resulting in different subsurface freezing–thawing conditions in the past depending on the assumed porosity and geothermal conditions.Our results suggest that in high altitude areas with a cold climate (present mean annual ground temperature between 0°C and −3°C), there may have been considerable variations in permafrost thickness (ranging from 0 to 150 m), as well as periods of no permafrost at all. The higher is the porosity of bedrock filled with ice, the stronger is the retarding effect of permafrost against climatic variations.Two-dimensional models including topographic effects with altitude-dependent ground temperatures and slope orientation and inclination dependent solar radiation were applied to a case of mountain permafrost in Ylläs, western Finnish Lapland, where bedrock permafrost is known to occur in boreholes to a depth of about 60 m. Modelling suggests complicated changes in permafrost thickness with time as well as contrasting situations on southern and northern slopes of the mountain.Extrapolating the climatic warming of the last 200 years to the end of the next century when the anticipated increase in the annual average air temperature is expected to be about 2 K indicates that the permafrost occurrences in bedrock in northern Fennoscandia would be thawing rapidly in low-porosity formations. However, already a porosity of 5% filled with ice would retard the thawing considerably.     

Dynamo model with a small number of modes and magnetic activity of T Tauri stars

We suggest a model based on the representation of the stellar magnetic field as a superposition of a finite number of poloidal and toroidal free decay modes to describe the dynamo action in fully convective stars. For the adopted law of stellar differential rotation, we determined the dynamo number in exceeding which the generation of a cyclically varying magnetic field is possible in stars without a radiative core and derived an expression for the period of the cycle. The dynamo cycles in fully convective stars and in stars with thin convective envelopes are shown to differ qualitatively: first, the distributions of spots in latitude during the cycle are different for these two types of stars and, second, the model predicts a great weakening of the spot formation in fully convective stars at certain phases of the cycle. To compare the theory with observations, we have analyzed the historical light curve for the weak-line T Tauri star V410 Tau and found that its long-term activity is not a well-defined cycle with a definite period—its activity is more likely quasi-cyclic with a characteristic time of ～4 yr and with a chaotic component superimposed. we have also concluded that a redistribution of spots in longitude is responsible for the secular brightness variations in the star. This does not allow the results of photometric observations to be directly compared with predictions of ourmodel, in which, for simplicity, we assumed a symmetry in longitude and investigated the temporal evolution of the spot distribution in latitude. Therefore, we discuss the questions of what and how observations can be compared with predictions of the dynamo theory.     

Validation of a multidecadal RCM hindcast over East Asia

The precipitation and low-level air temperature in East Asia from a regional climate model (RCM) hindcast for the 22-year period 1979–2000 is evaluated against observational data in preparation for the model use in regional climate change research. Emphasis of the evaluation is placed on the RCM capability in capturing the temporal and spatial variability of precipitation and low-level temperature, especially in conjunction with important climatological events such as, ENSO and East Asian monsoon, at three spatial scales of continental, subcontinental, and river basins.Spatial anomaly correlation time series of geopotential height and temperature show that the simulated upper-air fields remain consistent with the driving large-scale fields, NCEP Reanalysis 2 (R2), throughout the period. The simulated seasonal shifts in 850 hPa winds also agree well with R2 over eastern China and the western Pacific Ocean although the magnitudes of the shifts are overestimated, especially over the eastern slope of the Tibetan Plateau and in northern Manchuria. The simulated precipitation climatology agrees reasonably with that from two analysis datasets based on station- and remote-sensing data. Outstanding characteristics of precipitation including the location of the main rainband, climatological means, and the spatiotemporal variability in association with East Asian Monsoon, ENSO, and extreme events, are well represented in the hindcast. The most notable bias in the simulated precipitation is an overestimation of winter rainfall in southwestern coast of China, near the border with Vietnam. The simulation overestimates the interannual variability of seasonal precipitation especially in southern China, however, the corresponding coefficients of variation agree reasonably with observations except in very dry regions. This suggests that climate sensitivity of scaled precipitation can be useful for projecting climate change signals. The simulated low-level temperature climatology agrees reasonably with observational data as well. The most noticeable biases in the simulated low-level temperature are the warm (cold) biases in southern Siberia (northeastern China) during winter (summer) and the systematic underestimation of low-level temperature in the Tibetan Plateau for all seasons. The daily maximum temperature is underestimated for all seasons by 2−3 K with the largest biases in spring and fall except in the northwestern Mongolia region where it has been overestimated during winter. The daily minimum temperature biases ranges from 0.3 K in spring to 2 K in winter, and are much smaller than those in daily maximum temperature. The evaluation of the multidecadal hindcast shows that model errors mostly confined in the region near the lateral boundaries of the model domain with only minor biases in eastern China. This allows us to be cautiously optimistic about the RCM usefulness for studies of precipitation and low-level temperature changes in East Asia induced by increased emissions of greenhouse gases.     

Glacier retreat and climatic variability in the eastern Terskey–Alatoo, inner Tien Shan between the middle of the 19th century and beginning of the 21st century

Changes in the extent of glaciers and rates of glacier termini retreat in the eastern Terskey–Alatoo Range, the Tien Shan Mountains, Central Asia have been evaluated using the remote sensing techniques. Changes in the extent of 335 glaciers between the end of the Little Ice Age (LIA; mid-19th century), 1990 and 2003 have been estimated through the delineation of glacier outlines and the LIA moraine positions on the Landsat TM and ASTER imagery for 1990 and 2003 respectively. By 2003, the glacier surface area had decreased by 19% of the LIA value, which constitutes a 76 km² reduction in glacier surface area. Mapping of 109 glaciers using the 1965 1:25,000 maps revealed that glacier surface area decreased by 12.6% of the 1965 value between 1965 and 2003. Detailed mapping of 10 glaciers using historical maps and aerial photographs from the 1943–1977 period, has enabled glacier extent variations over the 20th century to be identified with a higher temporal resolution. Glacial retreat was slow in the early 20th century but increased considerably between 1943 and 1956 and then again after 1977. The post-1990 period has been marked by the most rapid glacier retreat since the end of the LIA. The observed changes in the extent of glaciers are in line with the observed climatic warming. The regional weather stations have revealed a strong climatic warming during the ablation season since the 1950s at a rate of 0.02–0.03 °C a^− 1. At the higher elevations in the study area represented by the Tien Shan meteorological station, the summer warming was accompanied by negative anomalies in annual precipitation in the 1990s enhancing glacier retreat. However, trends in precipitation in the post-1997 period cannot be evaluated due to the change in observational practices at this station. Neither station in the study area exhibits significant long-term trends in precipitation.     

A method for averaging the perturbing function in Hill’s problem

A modified method for averaging the perturbing function in Hill’s problem is suggested. The averaging is performed in the revolution period of the satellite over the mean anomaly of its motion with a full allowance for a variation in the position of the perturbing body. At its fixed position, the semimajor axis of the satellite orbit during the revolution of the satellite is constant in view of the evolution equations, while the remaining orbital elements undergo secular and long-period perturbations. Therefore, when the motion of the perturbing body is taken into account, the semimajor axis of the satellite orbit undergoes the strongest perturbations. The suggested approach generalizes the averaging method in which only the linear (in time) term is included in the perturbing function. This method requires no expansion in powers of time. The described method is illustrated by calculating the perturbations of the semimajor axes for two distant satellites of Saturn, S/2000 S 1 and S/2000 S5. An approximate analytic solution is compared with the results of numerical integration of the averaged system of equations of motion for these satellites.     

Determining the rotation law of the solar convection zone

The condition of minimum total dissipation is used to derive stationary rotation and azimuthal magnetic field distributions in the bulk of the solar convection zone with an upper boundary at which the relative radius is r/R=0.95. General equilibrium con figurations with symmetric and antisymmetric (about the equator) angular-velocity and field components are determined. The calculated rotation law matches the observed one in general parameters, but the decrease in angular velocity at high latitudes in theory is larger than that in observations. Besides, there are additional sharp variations in the rotation and field distributions in the theoretical curves near the generation zone of solar torsional waves. The possible cause of the latter discrepancy is discussed. The change in equilibrium distributions due to the presence of an inverse molecular-weight gradient at the base of the convection zone is also studied. This gradient is known to be produced by accelerated gravitational helium settling in the convection zone.     

The brightest OH maser in the sky: A flare of emission in W75 N

A flare of maser radio emission in the 1665-MHz OH line with a flux density of about 1000 Jy was discovered in the star-forming region W75 N in 2003. At the time of its observations, it was the strongest OH maser in the entire history of research since the discovery of cosmic OH masers in 1965. The linear polarization of the flare emission reached 100%. A weaker flare with a flux density of 145 Jy was observed in this source in 2000–2001; this was probably a precursor of the intense flare. The intensity of two other spectral features decreased when the flare emerged. This change in the intensity of the emission from maser condensations (a brightening of some of them and a weakening of others) can be explained by the passage of a magnetohydrodynamic shock through regions of enhanced gas concentration.     

类星体发射线和吸收线的统计分析

本文根据文[1，2]三个类星体表的3941个类星体作统计分析，发现在发射红移Zem1．4—3．2范围内吸收线较多，尤其在Zem1．8—2．4范围内明显增加；同时在Zem1．4—2．8范围内出现吸收线红移Zab＞Zem的现象．结果表明，在Zem1．8—2．2间隔内的类星体活动性加剧，似可以说反映了类星体与环境相互作用的演化效应．     

Can gas dynamics in centres of galaxies reveal orbiting massive black holes?

If supermassive black holes in centres of galaxies form by merging of black hole remnants of massive Population III stars, then there should be a few black holes of mass one or two orders of magnitude smaller than that of the central ones, orbiting around the centre of a typical galaxy. These black holes constitute a weak perturbation in the gravitational potential, which can generate wave phenomena in gas within a disc close to the centre of the galaxy. Here, we show that a single orbiting black hole generates a three-arm spiral pattern in the central gaseous disc. The density excess in the spiral arms in the disc reaches values of 3–12 per cent when the orbiting black hole is about 10 times less massive than the central black hole. Therefore, the observed density pattern in gas can be used as a signature in detecting the most massive orbiting black holes.     

Detecting protoplanets with ALMA

Theoretical investigations show that planet-disk interactions cause structures in circumstellar disks, which are usually much larger in size than the planet itself and thus more easily detectable. The specific result of planet-disk interactions depends on the evolutionary stage of the disk. Exemplary signatures of planets embedded in disks are gaps and spiral density waves in the case of young, gas-rich protoplanetary disks and characteristic asymmetric density patterns in debris disks. Numerical simulations convincingly demonstrate that high-resolution imaging performed with observational facilities which are already available or will become available in the near future will allow to trace these “fingerprints” of planets in protoplanetary and debris disks. These observations will provide a deep insight into specific phases of the formation and early evolution of planets in circumstellar disks. In this context, the Atacama Large Millimeter Array (ALMA) will play a crucial role by allowing to trace features in disks which are indicative for various stages of the formation and early evolution of planets in circumstellar disks.     

星系团的分层效应

对星系团各类分层效应的有关问题做了概要的评述,包括成员星系在位置空间和(或)速度空间中的形态分层、光度(质量)分层和元素丰度分层的表现形式和探测途径,分层效应可能的形成机制及其对星系和星系团的结构和演化的影响。     

Electrostatic Potentials in Supernova Remnant Shocks

Recent advances in the understanding of the properties of supernova remnant shocks have been precipitated by theChandra and XMM X-ray Observatories, and the HESS Atmospheric Čerenkov Telescope in the TeV band. A critical problem for this field is the understanding of the relative degree of dissipative heating/energization of electrons and ions in the shock layer. This impacts the interpretation of X-ray observations, and moreover influences the efficiency of injection into the acceleration process, which in turn feeds back into the thermal shock layer energetics and dynamics. This paper outlines the first stages of our exploration of the role of charge separation potentials in non-relativistic electron-ion shocks where the inertial gyro-scales are widely disparate, using results from a Monte Carlo simulation. Charge density spatial profiles were obtained in the linear regime, sampling the inertial scales for both ions and electrons, for different magnetic field obliquities. These were readily integrated to acquire electric field profiles in the absence of self-consistent, spatial readjustments between the electrons and the ions. It was found that while diffusion plays little role in modulating the linear field structure in highly oblique and perpendicular shocks, in quasi-parallel shocks, where charge separations induced by gyrations are small, and shock-layer electric fields are predominantly generated on diffusive scales.     

On the Dynamics of Dust Grains in a Hierarchical Environment

Most main sequence stars are binaries or higher multiplicity Systems and it appears that at birth most stars have circumstellar disks. It is commonly accepted that planetary systems arise from the material of these disks; consequently, binary and multiple systems may have a main role in planet formation. In this paper, we study the stage of planetary formation during which the particulate material is still dispersed as centimetre-to-metre sized primordial aggregates. We investigate the response of the particles, in a protoplanetary disk with radius R_D = 100 AU around a solar-like star, to the gravitational field of bound perturbing companions in a moderately wide (300–1600 AU) orbit. For this purpose, we have carried out a series of simulations of coplanar hierarchical configurations using a direct integration code that models gravitational and viscous forces. The massive protoplanetary disk is around one of the components of the binary. The evolution in time of the dust sub-disk depends mainly on the nature (prograde or retrograde) of the relative revolution of the stellar companion, and on the temperature and mass of the circumstellar disk. Our results show that for binary companions near the limit of tidal truncation of the disk, the perturbation leads to an enhanced accretion rate onto the primary, decreasing the lifetime of the particles in the protoplanetary disk with respect to the case of a single star. As a consequence of an enhanced accretion rate the mass of the disk decreases faster, which leads to a longer resultant lifetime for particles in the disk. On the other hand, binary companions may induce tidal arms in the dust phase of protoplanetary disks. Spiral perturbations with m = 1 may increase in a factor 10 or more the dust surface density in the neighbourhood of the arm, facilitating the growth of the particles. Moreover, in a massive disk (0.01M⊙) the survival time of particles is significantly shorter than in a less massive nebula (0.001M⊙) and the temperature of the disk severely influences the spiral-in time of particles. The rapid evolution of the dust component found in post T Tauri stars can be explained as a result of their binary nature. Binarity may also influence the evolution of circumpulsar disks. This revised version was published online in July 2006 with corrections to the Cover Date.