A machine learning approach for predicting computational intensity and domain decomposition in parallel geoprocessing

ABSTRACT

High performance computing is required for fast geoprocessing of geospatial big data. Using spatial domains to represent computational intensity (CIT) and domain decomposition for parallelism are prominent strategies when designing parallel geoprocessing applications. Traditional domain decomposition is limited in evaluating the computational intensity, which often results in load imbalance and poor parallel performance. From the data science perspective, machine learning from Artificial Intelligence (AI) shows promise for better CIT evaluation. This paper proposes a machine learning approach for predicting computational intensity, followed by an optimized domain decomposition, which divides the spatial domain into balanced subdivisions based on the predicted CIT to achieve better parallel performance. The approach provides a reference framework on how various machine learning methods including feature selection and model training can be used in predicting computational intensity and optimizing parallel geoprocessing against different cases. Some comparative experiments between the approach and traditional methods were performed using the two cases, DEM generation from point clouds and spatial intersection on vector data. The results not only demonstrate the advantage of the approach, but also provide hints on how traditional GIS computation can be improved by the AI machine learning.     

哀牢山-红河左旋走滑剪切带中新世抬升的时间序列

在哀牢山-红河剪切带的元阳到者龙的6个不同地点,对18个变质岩样品进行系统分析钾长石的~(40)Ar-~(39)Ar年龄。发现用MDD模式得到的18条冷却曲线(温度范围在400-150℃之间)均存在一次(或两次)快速冷却过程(冷却速度>100°C/Ma),这种快速冷却(抬升)过程不仅反映了该剪切带在大约25-17Ma时,韧性走滑作用结束和正断作用的开始,而且它的起始时间还有自东南向西北逐渐变新的趋势。利用6个不同地点快速冷却(抬升)的起始时间与相对距离作图时,得到了较好的线性函数关系(斜率为34mm/a)     

洪山泉流量动态系统分析

为查明洪山泉域的开采程度,在前人研究基础上,建立了现阶段泉流量-降雨量的线性模拟模型。通过对模型的系统辨识,对系统的"慢时变"现象作出了合理解释,并发现系统在1994年以后发生了"快时变",说明开采量从1995年开始迅猛增加。同时,由多年来泉流量序列的自相关分析,发现其宏观周期到2000年基本保持不变,说明虽然第二输入(即人为开采量)对该系统的作用已不可忽略,但尚未达到足以改变其宏观频率的程度。     

三峡库区四方碑滑坡稳定性与变形趋势预测

三峡水库建成后,库水位周期性涨落和暴雨产生的渗流作用导致大量古滑坡的复活或新滑坡的发生。以库区近水平层状结构的四方碑滑坡为例,依据库水位实际调动,将水位从175 m至145 m不同降速与50年一遇暴雨进行工况组合,计算4种工况下滑坡的稳定性及破坏概率。然后采用Geo-studio软件的Sigma模块对滑坡进行变形模拟,运用R/S分析方法判断滑坡的变形持续性,并结合野外调查情况,综合评价分析四方碑滑坡的稳定性。结果表明:滑坡在各工况下整体均处于基本稳定状态,具有低危险性;变形模拟结果显示滑坡前缘位移最大,与野外调查情况一致;各监测点Hurst指数均介于0.5~1,表明时间序列具有正持续性,在研究的时间限度内滑坡的局部破坏增强,应在汛期加强对滑坡前缘的巡查和预警。     

特提斯构造域东段叠合盆地演化和油气成藏规律特征初论

通过分析特提斯构造域东段区域地质和含油气盆地勘探开发基础数据,从板块构造演化入手,系统编制特提斯构造域东段沉积构造演化剖面图和生储盖组合剖面图,研究盆地演化阶段、叠合特征、油气成藏条件及油气藏类型,揭示中亚和中国西部前陆盆地演化和油气富集规律异同。研究表明：古亚洲洋、古特提斯洋和新特提斯洋控制了特提斯构造域东段的区域构造分带、盆地演化、盆地类型及油气成藏模式。根据古洋壳缝合线可分为北、中、南3个构造带,古生代以来多期微板块的拼贴,导致特提斯构造域东段含油气盆地演化分为3个演化阶段,早古生代伸展、晚古生代挤压、早中生代伸展和新生代挤压构造作用控制了研究区盆地的叠合演化,发育下古生界、上古生界和中生界3套区域分布的优质烃源岩和下古生界、上古生界、中生界和新生界4套储盖组合,形成多种类型的油气藏。     

青藏高原东边缘山区极值风速推算方法的研究

利用1971—2001年青藏高原东边缘地区12个气象台站风资料,分析了影响极值风速的主要因子,通过对地形的参数化处理,建立了极值风速随海拔高度和地形参数变化的拟合模型,可推算该地山区的极值风速,并利用山区临时观测资料对拟合模型进行了检验。结果表明,利用引入地形参数所建立的极值风速拟合模型来推算山区极值风速具有一定的可行性。     

Triangulated spherical splines for geopotential reconstruction

We present an alternate mathematical technique than contemporary spherical harmonics to approximate the geopotential based on triangulated spherical spline functions, which are smooth piecewise spherical harmonic polynomials over spherical triangulations. The new method is capable of multi-spatial resolution modeling and could thus enhance spatial resolutions for regional gravity field inversion using data from space gravimetry missions such as CHAMP, GRACE or GOCE. First, we propose to use the minimal energy spherical spline interpolation to find a good approximation of the geopotential at the orbital altitude of the satellite. Then we explain how to solve Laplace’s equation on the Earth’s exterior to compute a spherical spline to approximate the geopotential at the Earth’s surface. We propose a domain decomposition technique, which can compute an approximation of the minimal energy spherical spline interpolation on the orbital altitude and a multiple star technique to compute the spherical spline approximation by the collocation method. We prove that the spherical spline constructed by means of the domain decomposition technique converges to the minimal energy spline interpolation. We also prove that the modeled spline geopotential is continuous from the satellite altitude down to the Earth’s surface. We have implemented the two computational algorithms and applied them in a numerical experiment using simulated CHAMP geopotential observations computed at satellite altitude (450 km) assuming EGM96 (n _max = 90) is the truth model. We then validate our approach by comparing the computed geopotential values using the resulting spherical spline model down to the Earth’s surface, with the truth EGM96 values over several study regions. Our numerical evidence demonstrates that the algorithms produce a viable alternative of regional gravity field solution potentially exploiting the full accuracy of data from space gravimetry missions. The major advantage of our method is that it allows us to compute the geopotential over the regions of interest as well as enhancing the spatial resolution commensurable with the characteristics of satellite coverage, which could not be done using a global spherical harmonic representation. The results in this paper are based on the research supported by the National Science Foundation under the grant no. 0327577.     

一个有限差分大气环流模式的并行效率

适合于大规模并行处理机（ＭＰＰ）的中国科学院大气物理研究所新一代九层大气环流模式（ＩＡＰ ９Ｌ ＡＧＣＭ）被发展成功。本文描述了并行模式的设计原则，给出了该模式的并行版本在我国几种主要的并行机上的并行执行效率数据。采用消息传递方式（ＭＰＩ）的区域分解方法用于实现程序的并行，结果表明，对于该模式，只有一维区域分解方法才能取得较好的并行可扩展性，即随着处理器的增加而能得到相应较好的并行加速比和并行效率。     

精密时钟单钟幂律谱系数的取值分析

精密时钟噪声可用五个幂律谱系数来表征其强弱。由于不存在实际可用的理想时钟作为参考,用以估计幂律谱系数的单钟噪声无法直接测出,而只能利用钟组内各钟之间的钟差。然而这是一个不适定问题．本文探讨了单钟幂律谱系数的取值范围与优化估计方法。     

Clay mineral,geochemical and Sr–Nd isotopic fingerprinting of sediments in the Murray–Darling fluvial system,southeast Australia

The Willyama Supergroup of the Broken Hill region in southern Australia consists of supracrustal sedimentary and magmatic rocks, formed between 1810 and 1600 Ma. A statistical analysis of nearly 2000 SHRIMP U–Pb zircon spot ages, compiled from published and unpublished sources, provides evidence for three distinct tectonostratigraphic successions and four magmatic events during this interval. Succession 1 includes Redan Geophysical Zone gneisses and the lower part of the Thackaringa Group (Cues Formation). These rocks were deposited after 1810 Ma and host granite sills of the first magmatic event (1710–1700 Ma). Succession 2 includes the upper Thackaringa Group (Himalaya Formation), the Broken Hill Group and the Sundown Group and was deposited between 1710 and 1660 Ma. These rocks all contain detrital zircons from the first magmatic event (1710–1700 Ma) and in some cases from the second magmatic event (1690–1680 Ma). The second magmatic event (1690–1680 Ma) was bimodal, resulted from crustal extension, and was coeval with deposition of the Broken Hill Group and deepening of the basin. With this event a mafic sill swarm focused in the Broken Hill Domain. Mafic sills lack any trace of inheritance, unlike the granitoids that commonly contain inherited zircons typical of the supracrustal sediments. Succession 3, the Paragon Group and equivalents were deposited after 1660 Ma, but before a regional metamorphic event at 1600 Ma. Metamorphism was closely followed by inversion of the succession into a fold‐and‐thrust belt, accompanied by a fourth late to post‐orogenic magmatic event (ca 1580 Ma) characterised by granite intrusion and regional acid volcanism (the local equivalents of the Gawler Range Volcanics in South Australia).