球面四元三角网的三拓扑数计算

讨论了任意球面三角格网p的三拓扑数计算.任意球面三角格网p的三拓扑数是指和该球面三角格网互为三邻近的目标球面三角格网的个数,它是讨论球面栅格区域局部拓扑不变量的重要参数,也是描述和推断球面栅格拓扑关系首先要解决的问题.     

GGOS对流层延迟产品精度分析及在PPP中的应用

对流层延迟是卫星导航定位的主要误差源,GNSS广域增强需要高精度的对流层延迟产品进行误差修正。对流层延迟可通过GNSS进行实时估计,也可通过融合多源数据的数值气象预报模型获取。IGS发布的全球对流层天顶延迟产品由GNSS解算,其精度可达4mm,时间分辨率为5min,但其分布不均匀,在广袤的海洋区域无数据覆盖。GGOS Atmosphere基于ECMWF 40年再分析资料,可提供1979年以来时间分辨率为6h、空间分辨率为2.5°×2°的全球天顶对流层总延迟格网数据。本文通过2015年全球IGS测站的ZTD资料对GGOS的ZTD产品进行了评估,研究了GGOS Atmosphere对流层延迟产品与IGS发布ZTD资料之间的系统差,通过线性拟合估计出每个测站GGOS-ZTD与IGSZTD系统差系数(包括比例误差a和固定误差b),然后对比例误差a、固定误差b进行球谐展开,建立了两种ZTD数据源之间的系统差模型。选取IGS测站和陆态网测站,对附加系统偏差改正后的GGOSZTD产品对PPP的收敛速度的影响进行研究。本文研究结果表明:GGOS-ZTD与IGS-ZTD间存在系统偏差,其bias平均为-0.54cm;两者之间较差的RMS平均为1.31cm,说明GGOS-ZTD产品足以满足广大GNSS导航定位用户对对流层延迟改正的需要。将改正了系统差后的GGOS-ZTD产品用于ALBH、DEAR、ISPA测站、PALM测站、ADIS测站、YNMH测站、WUHN测站进行PPP试验,发现可明显提高定位收敛速度,尤其是在U方向上,收敛速度分别提高10.58%、31.68%、15.96%、43.89%、51.46%、14.69%、18.40%。     

A statistical analysis of the distribution of cordierite and biotite in hornfels from the Bugaboo contact aureole: implications for the kinetics of porphyroblast crystallization

The three‐dimensional disposition of cordierite and biotite crystals in a hornfels from the contact aureole of the Bugaboo Batholith is quantified using high‐resolution X‐ray micro‐computed tomography and global as well as scale‐dependent pattern statistics. The results demonstrate a random distribution of cordierite and biotite crystal sizes for all scales across the entire rock volume studied indicative of interface‐controlled prograde metamorphic reaction kinetics. The reaction considered responsible for the mineral assemblage and the formation of cordierite and biotite in the hornfels is Ms + Chl + Qtz = Crd + And + Bt + . Rock‐specific phase equilibria point to metamorphic conditions of ~520 –550 °C and 3 kbar for this reaction. The common approach to approximate the shape of crystals as spherical underestimates the influence of the Strauss hard‐core process on rock texture and may be misinterpreted to reflect ordering of crystal sizes by inhibition of nucleation and growth commonly associated with diffusion‐controlled reaction kinetics. According to our findings, Strauss hard‐core ordering develops at length scales equal to and less than the average major axis of the crystal population. This is significantly larger than what is obtained if a spherical crystal geometry would be assumed, and increases with deviation from sphericity. For the cordierite and biotite populations investigated in this research, Strauss hard‐core ordering developed at length scales of up to ～2.2 and 1.25 mm, respectively, which is almost 1 mm longer than the scales that would be obtained if a spherical geometry would have been assumed. Our results highlight the importance of a critical assessment of the geometrical model assumptions commonly applied in the three‐dimensional analysis of crystal size distributions, and underline the need for a quantitative understanding of interface processes in order to appreciate their role in the kinetics of contact metamorphic reactions and rock texture formation.     

A global spectral element model for poisson equations and advective flow over a sphere

A global spherical Fourier–Legendre spectral element method is proposed to solve Poisson equations and advective flow over a sphere. In the meridional direction, Legendre polynomials are used and the region is divided into several elements. In order to avoid coordinate singularities at the north and south poles in the meridional direction, Legendre–Gauss–Radau points are chosen at the elements involving the two poles. Fourier polynomials are applied in the zonal direction for its periodicity,with only one element. Then, the partial differential equations are solved on the longitude–latitude meshes without coordinate transformation between spherical and Cartesian coordinates. For verification of the proposed method, a few Poisson equations and advective flows are tested. Firstly, the method is found to be valid for test cases with smooth solution. The results of the Poisson equations demonstrate that the present method exhibits high accuracy and exponential convergence. Highprecision solutions are also obtained with near negligible numerical diffusion during the time evolution for advective flow with smooth shape. Secondly, the results of advective flow with non-smooth shape and deformational flow are also shown to be reasonable and effective. As a result, the present method is proved to be capable of solving flow through different types of elements, and thereby a desirable method with reliability and high accuracy for solving partial differential equations over a sphere.     

GNSS多尺度应变场估计及区域形变检测

利用GNSS数据,基于球面小波构建了GNSS多尺度应变场的估计模型,给出了球面小波位置和尺度的选取、模型正则化因子和参数的估计方法。为检验所构建模型的正确性,根据负位错理论正演了震间闭锁断层区域地表速度场作为模拟数据源,估计并分析了闭锁状态下走滑断层区域应变场分布特征,结果与实际吻合。模拟不同空间影响范围的闭锁逆冲断层形变,开展了多尺度应变场检测不同空间尺度形变信息的试验。结果表明,对于50km影响范围的小断层形变信号,在大尺度(第4—7尺度)应变场中并没有体现。而在小尺度(第8尺度)应变场中表现得非常明显。对于150km影响范围的大断层形变,在小尺度(第8尺度)应变场中仅出现了一小部分信息,而在大尺度(第4—7尺度)应变场中表现得更加完整和明显。由此得出,不同空间影响范围的地壳形变信息会在相应尺度的应变场中得以体现;小尺度应变场具有检测局部地壳形变信息的优势。     

起伏地表地震波场定向照明方法

基于惠更斯-菲涅尔原理和组合震源理论,提出了起伏地表情况下地震波场定向照明分析方法。该方法根据地表起伏形态和给定的照明方向进行坐标系变换,在新坐标系中计算组合震源的激发延时并模拟起伏地表定向的地震波场,从而得出地震波场源-检定向照明度。简单的二维起伏速度模型和起伏的Marmousi模型数值试验验证了该方法的有效性,表明通过坐标变换进行地震波场定向照明分析在地表起伏的勘探区域具有较大的适用性和灵活性,可以用来分析地下复杂构造对照明能量分布的贡献。     

Element chemostratigraphy of the Devonian/Carboniferous boundary – A compositional approach

The Devonian/Carboniferous (D/C) boundary is a critical interval in the Phanerozoic history, which is associated with vigorous climatic perturbations, continental glaciation, global sea-level fall and rapidly increased extinction rates in marine realms. In many sections world-wide, these global changes left a marked lithological signature, in particular the Hangenberg black shale (products of deep-shelf anoxia) and the overlying Hangenberg sandstone (sudden siliciclastic influx into predominantly carbonate depositional environments). Both layers bear a distinct geochemical signature. Even though either or both of these two lithologies are absent at many sections, their correlative counterparts can be indicated by subtle geochemical markers. We studied elemental geochemistry of fourteen D/C boundary sections in six key areas across Europe with the aim to select globally correlatable elemental proxy for the D/C boundary. Analysis of raw/log-transformed geochemical data (EDXRF, c.p.s. units), presenting the standard approach here, indicates that concentrations of terrigenous elements (Al, K, Rb, Ti and Zr) are mainly controlled by diluted Ca (carried by marine calcium carbonate) in limestone facies and, accordingly, their variations can be related to carbonate production in the sea rather than to terrigenous input from continent. Nevertheless, due to the relative nature of geochemical observations, reliance solely on statistical processing of raw data might lead to incomplete picture of multivariate data structure and/or biased interpretations. For this reason, the aim of this contribution is to discuss the logratio alternatives of the standard statistical methods, which may better reflect the relative nature of the data. For this purpose, principal component analysis was employed to reveal main geochemical patterns and while the geochemical signature of the D/C boundary was further analysed using Q-mode clustering that leads to predicative orthonormal logratio coordinates – balances. The comprehensive picture of the multivariate data structure provided by these statistical tools makes them a primary choice for exploratory compositional data analysis. At the same time, it turns out that the standard and compositional approaches have synergic effects. This fact can be extensively used in further geochemical studies.     

Diffusion-based cartogram on spheres

A planar cartogram is a two-dimensional map, on which the area of each closed region is in direct proportion to a chosen extensive property. To date, various algorithms have been proposed to construct planar cartograms. This work extends the two-dimensional, diffusion-based, topologically invariant cartogram algorithm proposed by Gastner and Newman onto spheres. Unlike its planar counterpart, the spherical formulation does not require boundary conditions and is invariant to the rotation of input data on the sphere. An implementation of this spherical cartogram transformation is designed to generate readable topology-preserving cartograms on spheres. Lastly, the method is illustrated with applications to global data such as worldwide human population, gross domestic product (purchasing power parity), carbon dioxide emissions and regional data such as the Electoral College of the United States presidential election of 2016.     

变水头入渗试验推求垂向渗透系数的计算方法

在Philip与Nestingen土壤水运动计算方法基础上,推导出描述单环变水头入渗试验中水深变化过程的计算公式,提出了确定土壤垂向渗透系数的MP方法.利用土壤非饱和数值模拟及室内土槽入渗试验得到的单环水位变化过程,验证了MP方法推求土壤垂向渗透系数的可靠性.结果表明:MP方法推求的土壤垂向渗透系数比目前采用的Philip方法、Nestingen方法及Green-Ampt方法精度更高,具有广泛的应用前景.     

全球对地观测球体幕显示中图像处理技术

覆盖全球的气象信息在球体幕上进行显示具有很强的立体感和纵深感,在多个行业具有应用价值。基于中国气象局的气象资料,在国内首次研制了全球对地观测球体幕显示系统软件并能够稳定运行。本文主要介绍了系统的核心技术图像处理的实现。图像处理包括两个步骤:图像增强和图像变形。图像增强处理增强了包含大量气象信息的图像显示效果,图像变形处理通过几何校正、图像拼接等技术实现将备投图像从平面显示正确投影到球体幕上。经过系统图像处理后,覆盖全球的观测资料能在球体幕上正确、无缝、无变形、无失真显示,其利用模拟真实地球的球体屏幕,让观众从太空的角度直观形象地体验全球气象信息展示。