Cosmic structure growth and dark energy

Dark energy has a dramatic effect on the dynamics of the Universe, causing the recently discovered acceleration of the expansion. The dynamics are also central to the behaviour of the growth of large-scale structure, offering the possibility that observations of structure formation provide a sensitive probe of the cosmology and dark energy characteristics. In particular, dark energy with a time-varying equation of state can have an influence on structure formation stretching back well into the matter-dominated epoch. We analyse this impact, first calculating the linear perturbation results, including those for weak gravitational lensing. These dynamical models possess definite observable differences from constant equation of state models. Then we present a large-scale numerical simulation of structure formation, including the largest volume to date involving a time-varying equation of state. We find the halo mass function is well described by the Jenkins et al. mass function formula. We also show how to interpret modifications of the Friedmann equation in terms of a time-variable equation of state. The results presented here provide steps toward realistic computation of the effect of dark energy in cosmological probes involving large-scale structure, such as cluster counts, the Sunyaev–Zel'dovich effect or weak gravitational lensing.     

Numerical Stabilization of Orbital Motion

Mainly, the author focuses on Baumgarte's method and its applications in satellite, asteroid, stellar and planetary problems. In the paper arguments are given for the use of energy relations for stabilization in the elliptical two-body problem. Stabilizing properties of Baumgarte's equations and others are discussed. A simple approach is proposed for stabilizing the equations of almost circular motion. By using Baumgarte's technique, the author derives stabilized equations of perturbed restricted three-body problem. It is shown experimentally that stabilization in the problems mentioned above can raise the accuracy of numerical integration by several orders.     

Coronal Holes and Icosahedral Symmetry

The study of the expansion of the solar wind out of a system of coronal holes is continued. To this end, we consider the numerical integration of partial differential equations for problems with icosahedral symmetry, in general. First, employing Weyl theory, orbifold coordinates are introduced. Second, the icosahedral coordinates are discussed in detail. Third, following an analysis of the properties of these coordinates and the derivation of a few expressions useful for grid construction, various alternatives for the distribution of lattice points required for numerical integration are considered. A comparison of these numerical grids motivates the choice of a specific grid optimized for the numerical integration carried out in the accompanying paper by Kalish et al.(2002). This revised version was published online in July 2006 with corrections to the Cover Date.     

介质横向不均匀对震后地表位移场的影响

利用三维有限元模拟两层地球模型(包括弹性层和粘弹性层)介质的横向不均匀和断层错动的复杂性对震后地表位移场的影响，其数值结果包括弹性解和粘弹性解。结果表明，震后位移场的分布特征不仅与粘弹性层的流变特性有关，还与弹性层介质的剪切模量以及断层的同震错动有关。     

Seismic rays in an ellipsoidal earth model: theoretical analysis and estimation of deviations due to the ellipticity

Summary. A first-order form of the Euler's equations for rays in an ellipsoidal model of the Earth is obtained. The conditions affecting the velocity law for a monotonic increase, with respect to the arc length, in the angular distance to the epicentre, and in the angle of incidence, are the same in the ellipsoidal and spherical models. It is therefore possible to trace rays and to compute travel times directly in an ellipsoidal earth as in the spherical model. Thus comparison with the rays of the same coordinates in a spherical earth provides an estimate of the various deviations of these rays due to the Earth's flattening, and the corresponding travel-time differences, for mantle P -waves and for shallow earthquakes. All these deviations are functions both of the latitude and of the epicentral distance. The difference in the distance to the Earth's centre at points with the same geocentric latitude on rays in the ellipsoidal and in the spherical model may reach several kilometres. Directly related to the deformation of the isovelocity surfaces, this difference is the only cause of significant perturbation in travel times. Other differences, such as that corresponding to the ray torsion, are of the first order in ellipticity, and may exceed 1 km. They induce only small differences in travel time, less than 0.01s. Thus, we show that the ellipticity correction obtained by Jeffreys (1935) and Bullen (1937) by a perturbational method can be recovered by a direct evaluation of the travel times in an ellipsoidal model of the Earth. Moreover, as stated by Dziewonski & Gilbert (1976), we verify the non-dependence of this correction on the choice of the velocity law.     

首都圈地区现今断层活动的数学模拟及地震危险性讨论

以首都圈地区现今活动断层上近２０年的位移测量资料为依据，用二维线弹性有限元对该区断层的活动特征进行了拟合，并结合有关资料讨论了该区的地震危险性。结果表明：１９７７－１９８６年间该区主压应力优势方位与华北较一致，约为Ｎ４５°－８０°Ｅ；１９８６－１９９０年其主压应力方向向北偏转，约为Ｎ５°－６０°Ｅ；１９９０年至现在其主压应力又向东偏转，角度大于第一阶段，约为Ｎ８０°－９５°Ｅ。该区的张家口－延庆一带近期有发生中强震的可能；丰镇－阳高－大同地区和凉城－古营盘地区应力较高，也应引起注意。     

密度和压缩系数的散射层析成像法

本文在速度成像的基础上研究了同时对密度和压缩系数成像的散射波层析成像法.对不同散射角度的计算可以得到一系列反演图像,拟合这些图像,从而可以有效地达到对密度和压缩系数（或速度）成像的目的.与单纯的速度成像相比,增加了反演的难度.首先是对资料的方位性要求增加;其次是对资料的利用率下降.即便如此,从对较少量的炮点和检波点资料的数值计算来看,仍取得了满意的成像结果.我们对组成字母“A”的散射体结构进行了成像计算,结果能够同时再现密度和压缩系数,成像清晰,表明了方法的可行性,并能应用于复杂结构的成像问题.     

青藏高原及邻近区域的S波三维速度结构

本文收集了WWSSN台网和我国台网中13个地震台站的长周期地震记录,用140条10-90s瑞利波频散曲线和作者提出的Tarantola-Backus面波频散层析成象方法,作了青藏高原及邻区的速度反演,得出该地区岩石层速度结构的三维图象.结果表明,1.在10-110km深度范围内,速度结构出现与大地构造特征相关的分区性,显示出四个构造单元:青藏块体、柴达木-巴颜喀拉-三江块体、塔里木块体和印度块体.2.高原内部,深度为10-70km内速度较低,莫霍界面呈不对称盆形分布,藏北那曲附近地壳厚度超过70km,高原边缘壳厚为45-50km,90-110km为高速异常,表明高原内部存在上地幔盖层.3.高原北部的班公湖断裂和东部的三江断裂系是该区重要的分界线,是岩石层结构存在明显差异的重要接触部位,可能是冈瓦纳古陆与欧亚古陆的缝合带.4.柴达木-巴颜喀拉-三江块体内部速度分布不均匀,地壳厚度由北向南从45km加深到60km;在深度90-110km存在一低速层.5.塔里木地块内速度随深度均匀增加,从地壳到上地幔110km内没有发现低速层.地壳厚度约50km.