Lower-crustal rifting in the Rukwa Graben, East Africa

An M_w 5.9 earthquake occurred in the Lake Rukwa rift, Tanzania, on 1994 August 18, and was well recorded by 20 broad-band seismic stations at distances of 160 to 800 km and 21 broad-band stations at teleseismic distances. The regional and teleseismic waveforms have been used to investigate the source characteristics of the main shock, and also to locate aftershocks that occurred within three weeks of the main shock. Teleseismic body-wave modelling yields the following source parameters for the main shock: source depth of 25 ± 2 km, a normal fault orientation, with a horizontal tension axis striking NE-SW and an almost vertical pressure axis (Nodal Plane I: strike 126°–142°, dip 63°–66°, and rake 280°–290°; Nodal Plane II: strike 273°–289°, dip 28°–31°, and rake 235°–245°), a scalar moment of 4.1 times 10¹⁷ N m, and a 2 s impulsive source time function. Four of the largest aftershocks also nucleated at depths of 25 km, as deduced from regional sPmp–Pmp times. The nodal planes are broadly consistent with the orientation of both the Lupa and Ufipa faults, which bound the Rukwa rift to the northeast and southwest, respectively. The rupture radius of the main shock, assuming a circular fault, is estimated to be 4 km with a corresponding stress drop of 6.5 MPa. Published estimates of crustal thickness beneath the Rukwa rift indicate that the foci of the main shock and aftershocks lie at least 10 km above the Moho. The presence of lower-crustal seismicity beneath the Rukwa rift suggests that the pre-rift thermal structure of the rifted crust has not been strongly modified by the rifting, at least to depths of 25 km.     

Numerical inversion of deformation caused by pressure sources: application to Mount Etna (Italy)

The interpretation of geodetic data in volcanic areas is usually based on analytical deformation models. Although numerical finite element (FE) modelling allows realistic features such as topography and crustal heterogeneities to be included, the technique is not computationally convenient for solving inverse problems using classical methods. In this paper, we develop a general tool to perform inversions of geodetic data by means of 3-D FE models. The forward model is a library of numerical displacement solutions, where each entry of the library is the surface displacement due to a single stress component applied to an element of the grid. The final solution is a weighted combination of the six stress components applied to a single element-source. The pre-computed forward models are implemented in a global search algorithm, followed by an appraisal of the sampled solutions. After providing extended testing, we apply the method to model the 1993–1997 inflation phase at Mt Etna, documented by GPS and EDM measurements. We consider four different forward libraries, computed in models characterized by homogeneous/heterogeneous medium and flat/topographic free surface. Our results suggest that the elastic heterogeneities of the medium can significantly alter the position of the inferred source, while the topography has minor effect.     

具梯度项的次线性椭圆型方程解的分析

本文研究了1类具梯度项的次线性椭圆型方程大解和完全有界解的存在性问题。运用上下解方法和极值原理分别得到了Rn上方程存在完全大解的充分必要条件和存在完全有界解的充分条件,并且证明了该方程在Rn中光滑有界区域Ω上不存在大解。     

非滑动式与滑动式拉格朗日插值法在BDS精密星历内插中的比较分析

在高精度卫星定位数据处理中,对精密星历内插处理是一项基础工作。利用拉格朗日多项式插值法,从三方面对北斗导航系统精密星历进行了内插精度分析,分别是非滑动式拉格朗日插值法在不同插值阶数下的插值精度,滑动式拉格朗日插值法在不同插值阶数下的插值精度,以及非滑动式和滑动式拉格朗日插值法在同一插值阶数下的插值精度。算例结果表明,如果将插值阶数作为横坐标,且向右为正,插值精度为纵坐标,向上为正,那么非滑动式拉格朗日插值法的插值精度与插值阶数呈现一个开口向下的抛物线关系;滑动式拉格朗日插值法的插值阶数高于7后,插值精度处于稳定且高的状态,插值误差最大值为1. 1 mm左右,均方差在0. 4 mm左右;在同一插值阶数下,滑动式拉格朗日插值法的精度高于非滑动式拉格朗日插值法的精度。     

A complex analytic solution for the attitude motion of a near-symmetric rigid body under body-fixed torques

Although analytic solutions for the attitude motion of a rigid body are available for several special cases, a comprehensive theory does not exist in the literature for the more complicated problems found in spacecraft dynamics. In the present paper, analytic solutions in complex form are derived for the attitude motion of a near-symmetric rigid body under the influence of constant body-fixed torques. The solution is very compact, which enables efficient and rapid machine computation. Numerical simulations reveal that the solution is very accurate when applied to typical spinning spacecraft problems.     

土压力滑楔理论的统一解

基于滑楔理论计算土压力的基本原理,运用考虑中间主应力效应的统一强度理论,对传统的土压力公式不能反映土体材料强度的中间主应力效应的公式进行了改进,推导了基于滑楔理论的土压力统一解公式,Coulomb和Rankine土压力公式为其特例.通过变换统一强度参数,可获得一系列的解.在工程运用中可根据土质和工程的实际情况,适当选择统一强度参数和应力状态来确定土压力的大小.     

Restricted Three-Body Problem: An Approximation of its General Solution – Part One – the Manifold of Symmetric Periodic Solutions

This paper gives the results of a programme attempting to exploit ‘la seule bréche’ (Poincaré, 1892, p. 82) of non-integrable systems, namely to develop an approximate general solution for the three out of its four component-solutions of the planar restricted three-body problem. This is accomplished by computing a large number of families of ‘solutions précieuses’ (periodic solutions) covering densely the space of initial conditions of this problem. More specifically, we calculated numerically and only for μ = 0.4, all families of symmetric periodic solutions (1st component of the general solution) existing in the domain D:(x ₀ ∊ [−2,2],C ∊ [−2,5]) of the (x ₀, C) space and consisting of symmetric solutions re-entering after 1 up to 50 revolutions (see graph in Fig. 4). Then we tested the parts of the domain D that is void of such families and established that they belong to the category of escape motions (2nd component of the general solution). The approximation of the 3rd component (asymmetric solutions) we shall present in a future publication. The 4th component of the general solution of the problem, namely the one consisting of the bounded non-periodic solutions, is considered as approximated by those of the 1st or the 2nd component on account of the `Last Geometric Theorem of Poincaré' (Birkhoff, 1913). The results obtained provoked interest to repeat the same work inside the larger closed domain D:(x ₀ ∊ [−6,2], C ∊ [−5,5]) and the results are presented in Fig. 15. A test run of the programme developed led to reproduction of the results presented by Hénon (1965) with better accuracy and many additional families not included in the sited paper. Pointer directions construed from the main body of results led to the definition of useful concepts of the basic family of order n, n = 1, 2,… and the completeness criterion of the solution inside a compact sub-domain of the (x ₀, C) space. The same results inspired the ‘partition theorem’, which conjectures the possibility of partitioning an initial conditions domain D into a finite set of sub-domains D i that fulfill the completeness criterion and allow complete approximation of the general solution of this problem by computing a relatively small number of family curves. The numerical results of this project include a large number of families that were computed in detail covering their natural termination, the morphology, and stability of their member solutions. Zooming into sub-domains of D permitted clear presentation of the families of symmetric solutions contained in them. Such zooming was made for various values of the parameter N, which defines the re-entrance revolutions number, which was selected to be from 50 to 500. The areas generating escape solutions have being investigated. In Appendix A we present families of symmetric solutions terminating at asymptotic solutions, and in Appendix B the morphology of large period symmetric solutions though examples of orbits that re-enter after from 8 to 500 revolutions. The paper concludes that approximations of the general solution of the planar restricted problem is possible and presents such approximations, only for some sub-domains that fulfill the completeness criterion, on the basis of sufficiently large number of families.