Preliminary Study on the Translunar Halo Orbits of the Real Earth–Moon System

The computation of translunar Halo orbits of the real Earth–Moon system (REMS) has been an open problem for a long time, but now, it is possible to compute Halo orbits of the REMS in a systematic way. In this paper, we describe the method used for the numerical computation of Halo orbits for a time span longer than 41 years. Halo orbits of the REMS are computed from quasi-periodic Halo orbits of the quasi-bicircular problem (QBCP). The QBCP is a model for the dynamics of a spacecraft in the Earth–Moon–Sun system. It is a Hamiltonian system with three degrees of freedom and depending periodically on time. In this model, Earth, Moon and Sun are moving in a self-consistent motion close to bicircular. The computed Halo orbits of the REMS are compared with the family of Halo orbits of the QBCP. The results show that the QBCP is a good model to understand the main features of the Halo family of the REMS.     

P-III分布参数的概率权重矩法S函数计算

概率权重矩法是一种估计统计分布参数的方法.本文根据不完全Γ函数在无限区间积分,推导了P-Ⅲ分布参数的S函数的计算公式.通过现有计算公式比较,其计算结果具有较高的计算精度,避免了的大量的数值积分计算.文中公式只要借助于计算编程进行求解,给定超几何函数项一定的截断误差,其运算具有较高的运行速度.文中计算公式是一种P-Ⅲ分布参数S函数的计算途径.     

Fitting the Linear Model of Coregionalization by Generalized Least Squares

In geostatistical studies, the fitting of the linear model of coregionalization (LMC) to direct and cross experimental semivariograms is usually performed with a weighted least-squares (WLS) procedure based on the number of pairs of observations at each lag. So far, no study has investigated the efficiency of other least-squares procedures, such as ordinary least squares (OLS), generalized least squares (GLS), and WLS with other weighing functions, in the context of the LMC. In this article, we compare the statistical properties of the sill estimators obtained with eight least-squares procedures for fitting the LMC: OLS, four WLS, and three GLS. The WLS procedures are based on approximations of the variance of semivariogram estimates at each distance lag. The GLS procedures use a variance–covariance matrix of semivariogram estimates that is (i) estimated using the fourth-order moments with sill estimates (GLS₁), (ii) calculated using the fourth-order moments with the theoretical sills (GLS₂), and (iii) based on an approximation using the correlation between semivariogram estimates in the case of spatial independence of the observations (GLS₃). The current algorithm for fitting the LMC by WLS while ensuring the positive semidefiniteness of sill matrix estimates is modified to include any least-squares procedure. A Monte Carlo study is performed for 16 scenarios corresponding to different combinations of the number of variables, number of spatial structures, values of ranges, and scale dependence of the correlations among variables. Simulation results show that the mean square error is accounted for mostly by the variance of the sill estimators instead of their squared bias. Overall, the estimated GLS₁ and theoretical GLS₂ are the most efficient, followed by the WLS procedure that is based on the number of pairs of observations and the average distance at each lag. On that basis, GLS₁ can be recommended for future studies using the LMC.     

Frequency analysis of nonstationary annual maximum flood series using the time‐varying two‐component mixture distributions

The most popular practice for analysing nonstationarity of flood series is to use a fixed single‐type probability distribution incorporated with the time‐varying moments. However, the type of probability distribution could be both complex because of distinct flood populations and time‐varying under changing environments. To allow the investigation of this complex nature, the time‐varying two‐component mixture distributions (TTMD) method is proposed in this study by considering the time variations of not only the moments of its component distributions but also the weighting coefficients. Having identified the existence of mixed flood populations based on circular statistics, the proposed TTMD was applied to model the annual maximum flood series of two stations in the Weihe River basin, with the model parameters calibrated by the meta‐heuristic maximum likelihood method. The performance of TTMD was evaluated by different diagnostic plots and indexes and compared with stationary single‐type distributions, stationary mixture distributions and time‐varying single‐type distributions. The results highlighted the advantages of TTMD with physically‐based covariates for both stations. Besides, the optimal TTMD models were considered to be capable of settling the issue of nonstationarity and capturing the mixed flood populations satisfactorily. Copyright © 2016 John Wiley & Sons, Ltd.     

跟踪和预报雷达强回波运动初探

利用数字化天气雷达资料，针对强对流天气下的回波块，提取面积Ｓ、形状描绘子Ｒ和矩不变量Ｍ３个特征量，作为判别因子，对强回波块进行跟踪，然后作出回波移动预报。对７次降水过程的实测雷达资料试验结果表明：该方法的跟踪和预报准确率是令人满意的。     

Appraisal of regional and index flood quantile estimators

Our results illustrate the performance of at-site and regional GEV/PWM flood quantile estimators in regions with different coefficients of variation, degrees of regional heterogeneity, record lengths, and number of sites. Analytic approximations of bias and variance are employed. For realistic GEV distributions and short records, the index-flood quantile estimator performs better than a 2-parameter GEV/PWM quantile estimator with a regional shape parameter, or a 3-parameter at-site GEV/PWM quantile estimator, in both humid and especially in arid regions, as long as the degree of regional heterogeneity is moderate. As regional heterogeneity or record lengths increases, 2-parameter estimators quickly dominate. Flood frequency models that assign probabilities larger than 2% to negative flows are unrealistic; experiments employing such distributions provide questionable results. This appraisal generally demonstrates the value of regionalizing estimators of the shape of a flood distribution, and sometimes the coefficient of variation.     

Low Energy Transfer to the Moon

In 1991, the Japanese Hiten mission used a low energy transfer with a ballistic capture at the Moon which required less Vthan a standard Hohmann transfer. In this paper, we apply the dynamical systems techniques developed in our earlier work to reproduce systematically a Hiten-like mission. We approximate the Sun–Earth–Moon-spacecraft 4-body system as two 3-body systems. Using the invariant manifold structures of the Lagrange points of the 3-body systems, we are able to construct low energy transfer trajectories from the Earth which execute ballistic capture at the Moon. The techniques used in the design and construction of this trajectory may be applied in many situations.     

Approaches for modeling longitudinal dispersion in pore-networks

The purpose of this study is to quantify the dispersivity in the longitudinal direction by upscaling pore scale mixing over a network domain and to verify the dispersivity with that obtained through the more rigorous upscaling technique, the Brownian particle tracking model (BPTM). We model a porous medium with a network of pore-units that are comprised of pore bodies and bonds of finite volume. Such a pore-unit is assumed to be a mixing cell with the steady state flow condition for a single fluid. Dispersivity can be obtained by solving the mixing cell model (MCM) for the concentration in each pore-unit and by averaging the concentrations for a large number of pore units (as a function of time and space). A minimal size of network that ascertains an asymptotic value of dispersivity was determined and verified with large size pore networks. This numerically computed dispersivity is compared with the results from the BPTM for the same porous medium and flow conditions. We show that the dispersivity obtained from the MCM is equally reliable for the heterogeneous pore-networks and can be estimated as a function of pore size heterogeneity. For homogeneous networks with the MCM, the iteration time step plays an important role. On the other hand, for networks with the BPTM, the assumption of intra-bond velocity profile affects the results.     

Walsh列率域中多维分形模型与GIS环境下地球物理信号处理

地球物理信号通常在多个尺度段表现尺度不变性，这些不变性起因于不同的地质、地球物理或成矿过程的自相似性. 利用这种在多个尺度段的尺度不变性可以设计多维分形滤波器，滤波所得信号表征了具尺度不变性的地质地球物理或成矿过程，可以用于成矿预测或环境评价. 本文研究了Walsh变换列率空间地球物理信号的列率功率谱密度与列率之间的分形与多维分形关系, 试验证实了大洋钻探、石油以及煤系地层地球物理测井资料在Walsh域的多维分形性质，提出了用于分解地球物理场,提取有用信号并用于矿产资源勘探或环境评价的多维分形W-A模型. 利用波列率域中的多维分形关系构造了W A图解(W-A Plot). 借助W-A图解可以确定最小平方误差（LS）意义下Walsh功率谱变化的不同自相似性的频率分界点,从而用于设计W-A分形滤波器. 这种滤波器可将地球物理场分解成具有不同自相似性的局部场,并且保留原场的各向异性结构. 与通常使用的基于Fourier变换的滤波技术相比,W-A模型具有许多优点：W-A适用于检测地球物理信号中的突变、线性、环状、局部与纹理结构等弱信号. 同时，由于Walsh变换中只有简单的变号(加法与减法)，其计算速度远快于建立在复数乘法之上的Fourier变换，所以W-A计算速度远快于Fourier域的滤波方法，可以用于地球物理信号的现场实时处理. 用加拿大Nova-cotia省西南地区的布格重力异常进行了W-A方法的试算，处理结果反映了地质、矿产分布规律，能够很好地进行矿产预测.