有限元法与伪谱法混合求解弹性波动方程

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.     

On cross-correlating weak lensing surveys

The present generation of weak lensing surveys will be superseded by surveys run from space with much better sky coverage and high level of signal-to-noise ratio, such as the Supernova/Acceleration Probe ( SNAP ). However, removal of any systematics or noise will remain a major cause of concern for any weak lensing survey. One of the best ways of spotting any undetected source of systematic noise is to compare surveys that probe the same part of the sky. In this paper we study various measures that are useful in cross-correlating weak lensing surveys with diverse survey strategies. Using two different statistics – the shear components and the aperture mass – we construct a class of estimators which encode such cross-correlations. These techniques will also be useful in studies where the entire source population from a specific survey can be divided into various redshift bins to study cross-correlations among them. We perform a detailed study of the angular size dependence and redshift dependence of these observables and of their sensitivity to the background cosmology. We find that one-point and two-point statistics provide complementary tools which allow one to constrain cosmological parameters and to obtain a simple estimate of the noise of the survey.     

Principal components in active galactic nuclei variability data and the estimation of the flux contributions from different components

It has been found that the near-infrared flux variations of Seyfert galaxies satisfy relations of the form   F_i ≈α _{i j} +β _{i j} F_j   , where F_i , F_j are the fluxes in filters i and j ; and  α _{i , j} , β _{i , j}   are constants. These relations have been used to estimate the constant contributions of the non-variable underlying galaxies. The paper attempts a formal treatment of the estimation procedure, allowing for the possible presence of a third component, namely non-variable hot dust. In an analysis of a sample of 38 Seyfert galaxies, inclusion of the hot dust component improves the model fit in approximately half the cases. All derived dust temperatures are below 300 K, in the range 540–860 K or above 1300 K. A noteworthy feature is the estimation of confidence intervals for the component contributions: this is achieved by bootstrapping. It is also pointed out that the model implies that such data could be fruitfully analysed in terms of principal components.     

Investigating the mechanisms leading to the deglaciation of past continental northern hemisphere ice sheets with the CLIMBER–GREMLINS coupled model

A coupling procedure between a climate model of intermediate complexity (CLIMBER-2.3) and a 3-dimensional thermo-mechanical ice-sheet model (GREMLINS) has been elaborated. The resulting coupled model describes the evolution of atmosphere, ocean, biosphere, cryosphere and their mutual interactions. It is used to perform several simulations of the Last Deglaciation period to identify the physical mechanisms at the origin of the deglaciation process. Our baseline experiment, forced by insolation and atmospheric CO₂, produces almost complete deglaciation of past northern hemisphere continental ice sheets, although ice remains over the Cordilleran region at the end of the simulation and also in Alaska and Eastern Siberia. Results clearly demonstrate that, in this study, the melting of the North American ice sheet is critically dependent on the deglaciation of Fennoscandia through processes involving switches of the thermohaline circulation from a glacial mode to a modern one and associated warming of the northern hemisphere. A set of sensitivity experiments has been carried out to test the relative importance of both forcing factors and internal processes in the deglaciation mechanism. It appears that the deglaciation is primarily driven by insolation. However, the atmospheric CO₂ modulates the timing of the melting of the Fennoscandian ice sheet, and results relative to Laurentide illustrate the existence of threshold CO₂ values, that can be translated in terms of critical temperature, below which the deglaciation is impeded. Finally, we show that the beginning of the deglaciation process of the Laurentide ice sheet may be influenced by the time at which the shift of the thermohaline circulation from one mode to the other occurs.     

CATASTROPHIC ECO-ENVIRONMENTAL CHANGE IN THE SONGNEN PLAIN, NORTHEASTERN CHINA SINCE 1900S

Although the Songnen Plain in the northeastern China was developed relatively late in the temperate zone of the world, its eco-environment has changed greatly. This paper analyzes the changes of land cover and the rates and trends ofdesertification during the past 100 years in the Songnen Plain. According to the macroscopic analysis, we find that the eco-environment in the plain has reached to the threshold of catastrophic change since the 1950s. The Thorn Needle Catastrophic Model was used to determine and validate this conclusion. Human activities, including large-scale construction projects, such as huge dams and dikes, and excessive grazing were the primary factors contributing to regional eco-environmental catastrophe. And irrational reclamation of the wilderness also affected the eco-environmental change. The results reveal the complex human-land interactions.     

任务轨道设计

本文首先说明太空任务与轨道设计的关系 ,接着介绍轨道的基本性质。从地球重力势的观点看各种常用的绕地轨道 ,包括地球和太阳同步轨道及Molniya轨道。从扰动的观点看常用的星际轨道 ,包括LISA、ASTROD、SOHO轨道。最后对星际轨道设计 ,说明二点边界值问题的数值解法、飞掠星体的应用、最佳化的考虑 ,并用以设计 2 0 1 5年发射的ASTROD初步任务轨道。     

Parameter information from non-linear cosmological fields

We develop a general formalism for analysing parameter information from non-Gaussian cosmic fields. The method can be adapted to include the non-linear effects in galaxy redshift surveys, weak lensing surveys and cosmic velocity field surveys as part of parameter estimation. It can also be used as a test of non-Gaussianity of the cosmic microwave background. Generalizing maximum-likelihood analysis to second order, we calculate the non-linear Fisher information matrix and likelihood surfaces in parameter space. To this order we find that the information content is always increased by including non-linearity. Our methods are applied to a realistic model of a galaxy redshift survey, including non-linear evolution, galaxy bias, shot-noise and redshift-space distortions to second order. We find that including non-linearities allows all of the degeneracies between parameters to be lifted. Marginalized parameter uncertainties of a few per cent will then be obtainable using forthcoming galaxy redshift surveys.     

On the r-mode spectrum of relativistic stars: the inclusion of the radiation reaction

We consider both mode calculations and time-evolutions of axial r modes for relativistic uniformly rotating non-barotropic neutron stars, using the slow-rotation formalism, in which rotational corrections are considered up to linear order in the angular velocity Ω. We study various stellar models, such as uniform density models, polytropic models with different polytropic indices n , and some models based on realistic equations of state. For weakly relativistic uniform density models and polytropes with small values of n , we can recover the growth times predicted from Newtonian theory when standard multipole formulae for the gravitational radiation are used. However, for more compact models, we find that relativistic linear perturbation theory predicts a weakening of the instability compared to the Newtonian results. When turning to polytropic equations of state, we find that for certain ranges of the polytropic index n , the r mode disappears, and instead of a growth, the time-evolutions show a rapid decay of the amplitude. This is clearly at variance with the Newtonian predictions. It is, however, fully consistent with our previous results obtained in the low-frequency approximation.     

An analytical expression for apsidal superhump precession and comparisons with numerical simulations and dwarf nova observations

We compare analytical expressions of precession rates from apsidal (positive) superhumps in close binary systems with numerical disc simulation results and observed values. In the analytical expressions, we include both the dynamical effects on the precession of the disc and effects caused by pressure forces that have been theorized to provide a retrograde effect (i.e. slowing) on the prograde disc precession. We establish new limits on density wave pitch angle to a normalized disc sound speed 60≥Ω_orb  d  tan  i / c >2.214 . Using average values for the density wave pitch angle i and speed of sound c , we find good correlation between numerical simulations and the analytical expression for the apsidal superhump period excess, which includes both the prograde and retrograde effects, for mass ratios of 0.025≤ q ≤0.33 . We also show good correlations with the four known eclipsing systems, OY Car, Z Cha, HT Cas, and WZ Sge. Our analytical expression for apsidal superhump period excess as a function of orbital period is consistent with the trend found in observed systems.     

THE CORPORATE INTERVIEW AS A RESEARCH METHOD IN ECONOMIC GEOGRAPHY*

The open-ended corporate interview as a qualitative research method is proposed as a valuable component of an evidentiary strategy in economic geography. It is argued to be more sensitive than other survey methods to historical, institutional, and strategic complexity. The corporate interview method is particularly appropriate in periods of economic and social change that challenge traditional analytical categories and theoretical principles. Some problems inherent to the method, and strategies for minimizing their impact on the research project, are described.