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.     

Characteristics of Spatial and Temporal Variations of Monthly Mean Surface Air Temperature over Qinghai-Tibet Plateau

1 Introduction The Qinghai-Tibet Plateau, known as the highest plateau with the most complex topography in the world, covers an area of more than 200km2, with a mean elevation of more than 4000m a.s.l. (Ye and Gao, 1979). Surrounded by the Earth’s highest mountains, such as the Himalayas, Pamir, Kunlun Mountains, the plateau plays a significant role in climate change in China even in the world, thus attracted great attention of researchers. Up to now, many achievements have been gained by…     

用遗传算法解算病态方程

对应用遗传算法解决病态方程问题进行了探讨。利用拟合法而不是通过法方程求解参数，从而避免了法方程系数求逆，使病态方程的解答有了较好的结果。通过模拟计算并和其他方法进行比较，证明该方法是可行的和有效的。     

地壳演化与成矿作用——以川滇地洼系“四层楼”铜矿床序列为例

提出川滇地洼系“四层楼”铜矿床序列的形成与陆壳演化的成生联系,是与本区陆壳由前地槽—地槽—地台—地洼演化各阶段与之相匹配的成矿作用的产物.与此同时,并总结了本区“四层楼”铜矿床序列的成矿作用具有明显的继承性、新生性、旋回性及层控性四大特点和多因复成矿床的成矿模式.     

Fluid entrapment and reequilibration during subduction and exhumation: A case study from the high-grade Nestos shear zone, Central Rhodope, Greece

A fluid inclusion study on metamorphic minerals of successive growth stages was performed on highly deformed paragneisses from the Nestos Shear Zone at Xanthi (Central Rhodope), in which microdiamonds provide unequivocal evidence for ultrahigh-pressure (UHP) metamorphism. The correlation of fluid inclusion density isochores and fluid inclusion reequilibration textures with geothermobarometric data and the relative chronology of micro- and macro-scale deformation stages allow a better understanding of both the fluid and metamorphic evolution along the P–T–d path. Textural evidence for subduction towards the NE is recorded by the orientation of intragranular NE-oriented fluid inclusion planes and the presence of single, annular fluid inclusion decrepitation textures. These textures occur within quartz “foam” structures enclosed in an earlier generation of garnets with prolate geometries and rarely within recrystallized matrix quartz, and reequilibrated both in composition and density during later stages of exhumation. No fluid inclusions pertaining to the postulated ultrahigh-pressure stage for microdiamond-bearing garnet–kyanite–gneisses have yet been found. The prolate shape of garnets developed during the earliest stages of exhumation that is recorded structurally by (L  S) tectonites, which subsequently accommodated progressive ductile SW shearing and folding up to shallow crustal levels. The majority of matrix kyanite and a later generation of garnet were formed during SW-directed shear under plane-strain conditions. Fluid inclusions entrapped in quartz during this stage of deformation underwent density loss and transformed to almost pure CO₂ inclusions by preferential loss of H₂O. Those inclusions armoured within garnet retained their primary 3-phase H₂O–CO₂ compositions. Reequilibration of fluid inclusions in quartz aggregates is most likely the result of recrystallization along with stress-induced, preferential H₂O leakage along dislocations and planar lattice defects which results in the predominance of CO₂ inclusions with supercritical densities. Carbonic fluid inclusions from adjacent kyanite–corundum-bearing pegmatoids and, the presence of shear-plane-parallel fluid inclusion planes within late quartz boudin structures consisting of pure CO₂-fluid inclusions with negative crystal shapes, bear witness of the latest stage of deformation by NE-directed extensional shear.This study shows that the textures of early fluid inclusions that formed already during the prograde metamorphic path can be preserved and used to derive information about the kinematics of subduction that is difficult to obtain from other sources. The textures of early inclusions, together with later generations of unaltered primary and secondary inclusions in metamorphic index minerals that can be linked to specific deformation stages and even P–T conditions, are a welcome supplement for the reconstruction of a rather detailed P–T–d path.     

Extensional Late Paleozoic deformation on the western margin of Pangea, Patlanoaya area, Acatlán Complex, southern Mexico

New mapping in the northern part of the Paleozoic Acatlán Complex (Patlanoaya area) records several ductile shear zones and brittle faults with normal kinematics (previously thought to be thrusts). These movement zones separate a variety of units that pass structurally upwards from: (i) blueschist-eclogitic metamorphic rocks (Piaxtla Suite) and mylonitic megacrystic granites (Columpio del Diablo granite ≡ Ordovician granites elsewhere in the complex); (ii) a gently E-dipping, listric, normal shear zone with top to the east kinematic indicators that formed under upper greenschist to lower amphibolite conditions; (iii) the Middle–Late Ordovician Las Minas quartzite (upper greenschist facies psammites with minor interbedded pelites intruded by mafic dikes and a leucogranite dike from the Columpio del Diablo granite) unconformably overlain by the Otate meta-arenite (lower greenschist facies psammites and pelites): roughly temporal equivalents are the Middle–Late Ordovician Mal Paso and Ojo de Agua units (interbedded metasandstone and slate, and metapelite and mafic minor intrusions, respectively) — some of these units are intruded by the massive, 461 ± 2 Ma, Palo Liso megacrystic granite: decussate, contact metamorphic muscovite yielded a ⁴⁰Ar/³⁹Ar plateau age of 440 ± 4 Ma; (iv) a steeply-moderately, E-dipping normal fault; (v) latest Devonian–Middle Permian sedimentary rocks (Patlanoaya Group: here elevated from formation status). The upward decrease in metamorphic grade is paralleled by a decrease in the number of penetrative fabrics, which varies from (i) three in the Piaxtla Suite, through (ii) two in the Las Minas unit (E-trending sheath folds deformed by NE-trending, subhorizontal folds with top to the southeast asymmetry, both associated with a solution cleavage), (iii) one in the Otate, Mal Paso, and Ojo de Agua units (steeply SE-dipping, NE–SW plunging, open-close folds), to (iv) none in the Patlanoaya Group. ⁴⁰Ar/³⁹Ar analyses of muscovite from the earliest cleavage in the Las Minas unit yielded a plateau age of 347 ± 3 Ma and show low temperature ages of  260 Ma. Post-dating all of these structures and the Patlanoaya Group are NE-plunging, subvertical folds and kink bands. An E–W, vertical normal fault juxtaposes the low-grade rocks against the Anacahuite amphibolite that is cut by megacrystic granite sheets, both of which were deformed by two penetrative fabrics. Amphibole from this unit has yielded a ⁴⁰Ar/³⁹Ar plateau age of 299 ± 6 Ma, which records cooling through  490 °C and is probably related to a Permo-Carboniferous reheating event during exhumation. The extensional deformation is inferred to have started in the latest Devonian ( 360 Ma) during deposition of the basal Patlanoaya Group, lasting through the rapid exhumation of the Piaxtla Suite at  350–340 Ma synchronous with cleavage development in the Las Minas unit, deposition of the Patlanoaya Group with active fault-related exhumation suggested by Mississippian and Early Permian conglomerates ( 340 and 300 Ma, respectively), and continuing at least into the Middle Permian (≡ 260 Ma muscovite ages). The continuity of Mid-Continent Mississippian fauna from the USA to southern Mexico suggests that this extensional deformation occurred on the western margin of Pangea after closure of the Rheic Ocean.     

Magnetotelluric investigations for imaging electrical structure of Garhwal Himalayan corridor,Uttarakhand, India

Magnetotelluric investigations have been carried out in the Garhwal Himalayan corridor to delineate the electrical structure of the crust along a profile extending from Indo-Gangetic Plain to Higher Himalayan region in Uttarakhand, India. The profile passing through major Himalayan thrusts: Himalayan Frontal Thrust (HFF), Main Boundary Thrust (MBT) and Main Central Thrust (MCT), is nearly perpendicular to the regional geological strike. Data processing and impedance analysis indicate that out of 44 stations MT data recorded, only 27 stations data show in general, the validity of 2D assumption. The average geoelectric strike, N70°W, was estimated for the profile using tensor decomposition. 2D smooth geoelectrical model has been presented, which provides the electrical image of the shallow and deeper crustal structure. The major features of the model are (i) a low resistivity (<50Ωm), shallow feature interpreted as sediments of Siwalik and Indo-Gangetic Plain, (ii) highly resistive (> 1000Ωm) zone below the sediments at a depth of 6 km, interpreted as the top surface of the Indian plate, (iii) a low resistivity (< 10Ωm) below the depth of 6 km near MCT zone coincides with the intense micro-seismic activity in the region. The zone is interpreted as the partial melting or fluid phase at mid crustal depth. Sensitivity test indicates that the major features of the geoelectrical model are relevant and desired by the MT data.     

Recharge source and hydrogeochemical evolution of shallow groundwater in a complex alluvial fan system,southwest of North China Plain

Many cities around the world are developed at alluvial fans. With economic and industrial development and increase in population, quality and quantity of groundwater are often damaged by over-exploitation in these areas. In order to realistically assess these groundwater resources and their sustainability, it is vital to understand the recharge sources and hydrogeochemical evolution of groundwater in alluvial fans. In March 2006, groundwater and surface water were sampled for major element analysis and stable isotope (oxygen-18 and deuterium) compositions in Xinxiang, which is located at a complex alluvial fan system composed of a mountainous area, Taihang Mt. alluvial fan and Yellow River alluvial fan. In the Taihang mountainous area, the groundwater was recharged by precipitation and was characterized by Ca–HCO₃ type water with depleted δ¹⁸O and δD (mean value of −8.8‰ δ¹⁸O). Along the flow path from the mountainous area to Taihang Mt. alluvial fan, the groundwater became geochemically complex (Ca–Na–Mg–HCO₃–Cl–SO₄ type), and heavier δ¹⁸O and δD were observed (around −8‰ δ¹⁸O). Before the surface water with mean δ¹⁸O of −8.7‰ recharged to groundwater, it underwent isotopic enrichment in Taihang Mt. alluvial fan. Chemical mixture and ion exchange are expected to be responsible for the chemical evolution of groundwater in Yellow River alluvial fan. Transferred water from the Yellow River is the main source of the groundwater in the Yellow River alluvial fan in the south of the study area, and stable isotopic compositions of the groundwater (mean value of −8.8‰ δ¹⁸O) were similar to those of transferred water (−8.9‰), increasing from the southern boundary of the study area to the distal end of the fan. The groundwater underwent chemical evolution from Ca–HCO₃, Na–HCO₃, to Na–SO₄. A conceptual model, integrating stiff diagrams, is used to describe the spatial variation of recharge sources, chemical evolution, and groundwater flow paths in the complex alluvial fan aquifer system.     

华北克拉通北缘隆化地区S型花岗岩的独居石年龄图谱

位于华北克拉通北缘中段的隆化S型花岗岩由石榴石黑云母花岗岩、石榴石花岗岩以及片麻理化的黑云母花岗岩组成。其主体岩性石榴石黑云母花岗岩SiO_2和Al_2O_3含量分别为64.09%～69.6%以及14.6%～16.13%,K_2O/Na_2O>1.0,A/CNK>1,0,Mg~#在20.76～34.89之间变化,具有明显的Nb、Ta、P、Ti和Sr亏损以及Rb、K和Th富集。石榴石黑云母花岗岩(样品JB6031-1)采用独居石电子探针U-Th-Pb化学法进行测年,获得了2553±120Ma、2180±42Ma和1854±24Ma三个年龄峰值。一颗独居石内部成分分带上6个分析点定年结果构成2553±120Ma的峰值年龄,这一年龄与我们最新获得的2506±7Ma和2541±8Ma(继承锆石年龄)LA-ICP-MS锆石U-Pb同位素年龄相似,我们将这一独居石年龄解释为继承独居石的年龄,表明在赤城-隆化断裂以北存在太古宙陆块,并且在后期构造-热事件中发生部分熔融形成S型花岗岩。该独居石颗粒幔部成分分带上10个分析点的测年结果揭示的峰值年龄为2181±42Ma,该年龄也是出现频率最高的年龄值,我们建议2181±42Ma为S型花岗岩的结晶年龄,反映了S型花岗岩的侵位时代。独居石颗粒外部成分分带上8个分析点的测年结果构成1854±24Ma的峰值年龄,该年龄与华北克拉通中部带的变质年龄接近,我们将其解释为S型花岗岩的变质年龄,表明华北克拉通北缘的构造演化与中部带的构造演化密切相关。     

热带地区积云对流的长期变化特征

使用１９７９年１月至１９８４年１２月向外长波辐射（ＯＬＲ）资料，对热带地区积云对流的长期变化特征进行了研究。结果表明：热带地区积云对流活动存在显著的季节变化，冬季积云对流区主要是东西向，位于南印度洋和西太平洋的近赤道地区；夏季则北移至北印度洋和菲律宾附近的西太平洋地区。低纬地区积云对流活动存在明显的季节性位移，北印度洋地区的积云对流活动主要集中在５—１０月，７—８月位置最北；北半球热带西太平洋地区的积云对流活动则主要集中在６—１１月，８—９月位置最北。标准差分析表明，冬季北半球热带西太平洋、赤道中太平洋及热带印度洋东部地区积云对流的年际变化最明显。经验正交函数（ＥＯＦ）分析的主要空间型反映了赤道中太平洋、热带西太平洋、阿拉伯海和副热带西太平洋地区的积云对流活动存在一定的关系。结合遥相关计算还表明秋云对流存在４种遥相关型，即２种东西向偶极型涛动型、西太平洋型和北印度洋型。