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.     

ECOLOGICAL SERIES OF SOIL ANIMALS IN DARLIDAI MOUNTAIN

The ecological series of soil animals under the broad-leaved and pine mixed forest in Darlidai Mountainwas studied. Seven sample plots were selected according to different altitude gradients, which belong to different vegeta-tion types. By investigating and analyzing soil animals in every sample plot it is found that there are 45 groups and 1956individuals, which axe involved in 3 phylums, 7 classes, 16 orders, respectively. The altitude is a key factor which af-fects ecological series of soil animals. Both the groups and individuals of soil animals increase with altitude increasingunder certain conditions, which contrastes with ordinary cases, resulting from special micro-climate in studied area. Thegroups and individuls of soil animals are the most under the broad-leaved and pine forest on the top of the mountain, andthe least under Picea-Abies forest in the foot of the mountain.     

对青藏东北缘现今块体划分、运动及变形的初步研究

利用2维非连续变形分析方法（DDA），以位移代替围压作为边界约束力，研究青藏东北缘现今块体划分及其运动变形。根据该地区地质构造及地震活动，以GPS点测量位移作为模拟结果约束点，得出了较合理的块体划分模型和随时间演化的主应变分布图，并把应变高值区与近几年来发生的5级以上地震作对比，得出了研究区内地震危险性可能较大的区域。另外，对模拟的甘青块体与阿拉善块体的边缘带断裂左旋运动做了大概计算。     

Revisiting Sustainable Development of Dry Valleys in Hengduan Mountains Region

Dry valleys are a striking geographic landscape in Hengduan Mountains Region and are characterized by low rainfall, desert type of vegetation and fragile environment. Past efforts and resources have been concentrated mainly on rehabilitation of degraded ecosystem and fragile environment, particularly reforestation, while socio-economic development has been largely overlooked. Despite successes in pocket areas, the overall trend of unsustainability and environmental deterioration are continuing. It is important to understand that uplift of the Tibetan Plateau is the root cause of development of dry valleys, and development and formation of dry valleys is a natural process. Human intervention has played a secondary role in development of dry valleys and degradation of dry valleys though human intervention in many cases has speeded up environmental degradation of the dry valleys. It is important to understand that dry valleys are climatic enclaves and an integrated approach that combines rehabilitation of degraded ecosystems and socio-economic development should be adopted if the overall goal of sustainable development of dry valleys is to be achieved. Promotion of niche-based cash crops, rural energy including hydropower, solar energy, biogas and fuelwood plantation is recommended as the priority activities.     

CPLD技术及其在现代时统中的应用

对CPLD(复杂可编程逻辑器件)技术的基本特征和发展趋势作了简要介绍,揭示了该技术在现代数字系统中的重要地位及作用.利用CPLD对时统设备IRIG-B码产生器进行集成,其实验结果表明,集成了的B码产生器不但简单、可靠,而且便于调试,克服了以往硬件电路复杂的缺点.     

Effect of spatial variation characteristics on contouring of design storm depth

The ordinary kriging method, a geostatistical interpolation technique, was applied for developing contour maps of design storm depth in northern Taiwan using intensity–duration–frequency (IDF) data. Results of variogram modelling on design storm depths indicate that the design storms can be categorized into two distinct storm types: (i) storms of short duration and high spatial variation and (ii) storms of long duration and less spatial variation. For storms of the first category, the influence range of rainfall depth decreases when the recurrence interval increases, owing to the increasing degree of their spatial independence. However, for storms of the second category, the influence range of rainfall depth does not change significantly and has an average of approximately 72 km. For very extreme events, such as events of short duration and long recurrence interval, we do not recommend usage of the established design storm contours, because most of the interstation distances exceed the influence ranges. Our study concludes that the influence range of the design storm depth is dependent on the design duration and recurrence interval and is a key factor in developing design storm contours. Copyright © 2003 John Wiley & Sons, Ltd.     

富氏谱分析法模拟地壳垂直运动的地学依据及其策应力模型

应用地壳波浪与镶嵌构造学说对富氏谱分析法提取地壳垂直形变信息的科学性做了地质学意义上的阐释 ,并提出了根据多期形变资料提取特定波段上构造策应力的数学模型     

Optical Flash of GRB 990123: Constraints on the Physical Parameters of the Reverse Shock

The optical flash accompanying GRB 990123 is believed to be powered by the reverse shock of a thin shell. With the best-fit physical parameters for GRB 990123 and the assumption that the parameters in the optical flash are the same as in the afterglow, we show that: 1) the shell is thick rather than thin, and we have provided the light curve for the thick shell case which coincides with the observation; 2) the theoretical peak flux of the optical flash accounts for only 3×10~-4 of the observed. In order to remove this discrepancy, the physical parameters, the electron energy and magnetic ratios, εe and εB, should be 0.61 and 0.39, which are very different from their values for the late afterglow.     

Silicon limitation on primary production and its destiny in Jiaozhou Bay, China——Ⅳ：Study on cross-bay transect from estuary to ocean

The authors analyzed the data collected in the Ecological Station Jiaozhou Bay from May 1991 to November 1994, including 12 seasonal investigations, to determine the characteristics, dynamic cycles and variation trends of the silicate in the bay. The results indicated that the rivers around Jiaozhou Bay provided abundant supply of silicate to the bay. The silicate concentration there depended on river flow variation. The horizontal variation of silicate concentration on the transect showed that the silicate concentration decreased with distance from shorelines. The vertical variation of it showed that silicate sank and deposited on the sea bottom by phytoplankton uptake and death, and zooplankton excretion. In this way, silicon would endlessly be transferred from terrestrial sources to the sea bottom. The silicon took up by phytoplankton and by other biogeochemical processes led to insufficient silicon supply for phytoplankton growth. In this paper, a 2D dynamic model of river flow versus silicate concentration was established by which silicate concentrations of 0.028–0.062 μmol/L in seawater was yielded by inputting certain seasonal unit river flows (m³/s), or in other words, the silicate supply rate; and when the unit river flow was set to zero, meaning no river input, the silicate concentrations were between 0.05–0.69 μmol/L in the bay. In terms of the silicate supply rate, Jiaozhou Bay was divided into three parts. The division shows a given river flow could generate several different silicon levels in corresponding regions, so as to the silicon-limitation levels to the phytoplankton in these regions. Another dynamic model of river flow versus primary production was set up by which the phytoplankton primary production of 5.21–15.55 (mgC/m²·d)/(m³/s) were obtained in our case at unit river flow values via silicate concentration or primary production conversion rate. Similarly, the values of primary production of 121.98–195.33 (mgC/m²·d) were achieved at zero unit river flow condition. A primary production conversion rate reflects the sensitivity to silicon depletion so as to different phytoplankton primary production and silicon requirements by different phytoplankton assemblages in different marine areas. In addition, the authors differentiated two equations (Eqs. 1 and 2) in the models to obtain the river flow variation that determines the silicate concentration variation, and in turn, the variation of primary production. These results proved further that nutrient silicon is a limiting factor for phytoplankton growth. This study was funded by NSFC (No. 40036010), and the Director's Fund of the Beihai Sea Monitoring Center, the State Oceanic Administration.     

Is Pre-Eruptive Null Point Reconnection Required for Triggering Eruptions?

We study the magnetic field evolution and topology of the active region NOAA 10486 before the 3B/X1.2 flare of October 26, 2003, using observational data from the French–Italian THEMIS telescope, the Michelson Doppler Imager (MDI) onboard Solar and Heliospheric Observatory (SOHO), the Solar Magnetic Field Telescope (SMFT) at Huairou Solar Observation Station (HSOS), and the Transition Region and Coronal Explorer (TRACE). Three dimensional (3D) extrapolation of photospheric magnetic field, assuming a potential field configuration, reveals the existence of two magnetic null points in the corona above the active region. We look at their role in the triggering of the main flare, by using the bright patches observed in TRACE 1600 Å images as tracers at the solar surface of energy release associated with magnetic reconnection at the null points. All the bright patches observed before the flare correspond to the low-altitude null point. They have no direct relationship with the X1.2 flare because the related separatrix is located far from the eruptive site. No bright patch corresponds to the high-altitude null point before the flare. We conclude that eruptions can be triggered without pre-eruptive coronal null point reconnection, and the presence of null points is not a sufficient condition for the occurrence of flares. We propose that this eruptive flare results from the loss of equilibrium due to persistent flux emergence, continuous photospheric motion and strong shear along the magnetic neutral line. The opening of the coronal field lines above the active region should be a byproduct of the large 3B/X1.2 flare rather than its trigger.