Testing variograms for positive-definiteness

A method based on Bochner's theorem is described for demonstrating the positive-definiteness of variogram models and for generating classes of valid variogram functions.Work performed while on leave at the Centre de Géostatistique et de Morphologie Mathématique.     

二维梯形沉积盆地放大效应的数值模拟研究

盆地效应作为场地效应的一种特殊类型,其对地震动的影响仍在持续研究当中。建立二维梯形沉积盆地模型,模型边界采用透射人工边界,利用显式有限元法作为求解方法。根据Buckingham提出的π定理定义了多个影响盆地效应的无量纲参数,包括盆地深宽比、介质波速比、盆地深度与入射波长之比、盆地开口宽度与入射波长之比。以Ricker子波作为输入波,研究SH波垂直入射盆地模型时,以上无量纲参数对盆地放大效应的影响。研究发现盆地深宽比越大盆地诱发的聚焦效应越明显,盆地深宽比越小盆地边缘效应越明显。在盆地的深宽比和入射波确定时,随着盆地开口宽度与入射波长之比不断增大盆地的最值放大系数逐渐从盆地中心移向盆地边缘,表现出越来越明显的边缘放大效应。其中波速比差异对放大系数的影响最大,波速比差异越大介质的波阻抗比越大,放大效应越显著。通过将影响盆地放大效应的参数无量纲化,使梯形盆地对地震动放大效应的研究的结论更具有普遍意义。     

The evolution of zircon during low‐P partial melting of metapelitic rocks: theoretical predictions and a case study from Mt Stafford,central Australia

Palaeoproterozoic metasedimentary migmatite reflects the highest temperature parts of a regional aureole at Mt Stafford, central Australia, comprising rocks that experienced 500–800 °C at ≈3 kbar. Whole‐rock major element concentrations are correlated with Zr content, psammitic compositions having nearly twice the Zr content of pelitic compositions. Zirconium is concentrated in mesosome compared with leucosome. Zircon is largely detrital, mostly lacking any overgrowth contemporary with migmatite formation. Comparatively small proportions of micro‐zircon (<10 μm) in sub‐solidus rocks are mostly hosted by quartz and plagioclase. Much higher proportions (three to five times) of micro‐zircon in migmatite are hosted by prograde K‐feldspar, cordierite and biotite. T–X and P–T NCKFMASHTZr pseudosections constructed using thermocalc model the distribution of Zr between solid and silicate liquid phases. Half of the detrital zircon (~100 ppm Zr) is predicted to be dissolved into silicate liquid at ≈800 °C and all dissolved by 850 °C, if all zircon is involved in the equilibration volume. Melt segregation at relatively low temperature is predicted to enrich the residuum in Zr, consistent with the observed distribution of Zr between mesosome and leucosome. The limited development of metamorphic zircon rims or overgrowths at Mt Stafford is explained by three concurrent processes: (i) Zr liberated during prograde metamorphism formed micro‐zircon, rather than following the prediction that Zr will partition into silicate liquid; (ii) some detrital zircon was probably armoured by other rock‐forming minerals, reducing Zr content in the effective bulk rock composition; and (iii) small proportions of melt loss during migmatization removed Zr that otherwise would have been available to form metamorphic rims.     

Solution to the Stokes Boundary-Value Problem on an Ellipsoid of Revolution

We have constructed Green's function to Stokes's boundary-value problem with the gravity data distributed over an ellipsoid of revolution. We show that the problem has a unique solution provided that the first eccentricity e₀ of the ellipsoid of revolution is less than 0·65041. The ellipsoidal Stokes function describing the effect of ellipticity of the boundary is expressed in the E-approximation as a finite sum of elementary functions which describe analytically the behaviour of the ellipsoidal Stokes function at the singular point = 0. We prove that the degree of singularity of the ellipsoidal Stokes function in the vicinity of its singular point is the same as that of the spherical Stokes function.     

Characterizing mixed-use buildings based on multi-source big data

To-date few research has successfully integrated big data from multiple sources to characterize urban mixed-use buildings. In this paper, we introduce a probabilistic model to integrate multi-source and geospatial big data (social network data, taxi trajectories, Points of Interest and remote sensing images) to characterize urban mixed-use buildings. The usefulness of our model is demonstrated with a case study of the Tianhe District in megacity Guangzhou, China. The model predicted building functions at 85% accuracy based on ground truth data from field surveys. We further explored the spatial patterns of the identified building functions. Most mixed-use buildings are located along major streets. Our proposed model can identify mixed-use buildings in a city; information is useful for planning evaluation and urban policymaking.     

On the robustness of the toroidal velocity theorem

In dynamo theory the toroidal velocity theorem in its classical version (Elsasser, Phys. Rev. 1946, vol. 69, pp. 106–116, Bullard and Gellman, Phil. Trans. R. Soc. Lond. 1954, vol. A247, pp. 213–278) rules out dynamo action in a spherical conducting volume provided that the fluid is incompressible, the conductivity is uniform, and the velocity field is purely toroidal. We prove in this note that this result is robust in the sense that slight compressibility of the fluid, small non-radial variations and even large radial variations in conductivity, and the presence of a small non-toroidal velocity component do not invalidate the theorem. Moreover, by proper choice of the conductivity distribution modelling the conducting volume, small deviations from spherical symmetry of the conductor can also be taken into account.     

透射边界高频耦合失稳机理及稳定实现——SH波动

就大型近场波动的高效数值模拟而言,稳定实现高阶人工边界是一个尚未圆满解决的问题.本文针对使用多次透射公式的SH波动集中质量有限元模拟,依据GKS定理的群速度解释,进一步阐明了人工边界与内域离散格式耦合所导致高频失稳的机理,即两者支持群速度指向内域的外行高频平面谐波,波动能量自发地从人工边界进入內域,从而导致失稳,而这类谐波是由集中质量有限元离散引入的.本文提出了消除此种耦合失稳的一种方法:通过修改有限元刚度阵来改变内域离散格式,并保证修改格式的精度不低于原有格式的精度.理论分析和数值实验表明此法能稳定实现透射边界.本文研究结果具有推广应用前景.     

基于贝叶斯定理的地震危险性概率预测研究

由于地震孕育过程的复杂性和观测技术的局限性,不同地震观测资料表现出异常变化与后续较大地震的对应关系存在不确定性,因此对预测意见进行概率表达是一种科学恰当的做法。本文基于泊松分布的危险区背景地震概率预测和单项预测方法(包括测震、流体、形变、电磁等学科)的历史预测效能,采用贝叶斯定理计算得到单项预测方法的短期或年度地震危险概率预测结果,进而采用综合概率方法,给出基于多种单项预测方法的短期或年度地震危险概率预测结果。短期概率预测初步结果表明,2018年2～9月,中国大陆72%的5级以上地震都位于相对高概率预测区域。     

Homogeneous planar and two-dimensional mean-field antidynamo theorems with zero mean flow

In an electrically conducting fluid, two types of turbulence with a preferred direction are distinguished: planar turbulence, in which every velocity in the turbulent ensemble of flows has no component in the given direction; and two-dimensional turbulence, in which every velocity in the turbulent ensemble is invariant under translation in the preferred direction. Under the additional assumptions of two-scale and homogeneous turbulence with zero mean flow, the associated magnetohydrodynamic alpha- and beta-effects are derived in the second-order correlation approximation (SOCA) when the electrically conducting fluid occupies all space. Limitations of the SOCA are well known, but alpha- and beta-effects of a turbulent flow are useful in interpreting the dynamo effects of the turbulence. Two antidynamo theorems, which establish necessary conditions for dynamo action, are shown to follow from the special structures of these alpha- and beta-effects. The theorems, which are analogues of the laminar planar velocity and two-dimensional antidynamo theorems, apply to all turbulent ensembles with the prescribed alpha- and beta-effects, not just the planar and two-dimensional ensembles. The mean magnetic field is general in the planar theorem but only two-dimensional in the two-dimensional theorem. The two theorems relax the previous restriction to turbulence which is both two-dimensional and planar. The laminar theorems imply decay of the total magnetic field for any velocity of the associated turbulent ensemble. However, the mean-field theorems are not fully consistent with the laminar theorems because further conditions beyond those arising from the turbulence must be imposed on the beta-effect to establish decay of the mean magnetic field. In particular, negative turbulent magnetic diffusivities must be restricted. It is interesting that there is no inconsistency in the alpha-effects. The failure of the SOCA with the two-scale approximation to simply preserve the laminar antidynamo theorems at the beta-effect level is a further demonstration of the restricted validity of the theory and shows that negative diffusivity effects derived by approximation methods must be treated cautiously.