Bias estimation and correction for triangle-based surface area calculations

The calculation of surface area is meaningful for a variety of space-filling phenomena, e.g., the packing of plants or animals within an area of land. With Digital Elevation Model (DEM) data we can calculate the surface area by using a continuous surface model, such as by the Triangulated Irregular Network (TIN). However, just as the triangle-based surface area discussed in this paper, the surface area is generally biased because it is a nonlinear mapping about the DEM data which contain measurement errors. To reduce the bias in the surface area, we propose a second-order bias correction by applying nonlinear error propagation to the triangle-based surface area. This process reveals that the random errors in the DEM data result in a bias in the triangle-based surface area while the systematic errors in the DEM data can be reduced by using the height differences. The bias is theoretically given by a probability integral which can be approximated by numerical approaches including the numerical integral and the Monte Carlo method; but these approaches need a theoretical distribution assumption about the DEM measurement errors, and have a very high computational cost. In most cases, we only have variance information on the measurement errors; thus, a bias estimation based on nonlinear error propagation is proposed. Based on the second-order bias estimation proposed, the variance of the surface area can be improved immediately by removing the bias from the original variance estimation. The main results are verified by the Monte Carlo method and by the numerical integral. They show that an unbiased surface area can be obtained by removing the proposed bias estimation from the triangle-based surface area originally calculated from the DEM data.     

A fine-granularity scheduling algorithm for parallel XDraw viewshed analysis

Viewshed analysis is widely used in many terrain applications such as siting problem, path planning problem, and etc. But viewshed computation is very time-consuming, in particular for applications with large-scale terrain data. Parallel computing as a mainstream technique with the tremendous potential has been introduced to enhance the computation performance of viewshed analysis. This paper presents a revised parallel viewshed computation approach based on the existing serial XDraw algorithm in a distributed parallel computing environment. A layered data-dependent model for processing data dependency in the XDraw algorithm is built to explore scheduling strategy so that a fine-granularity scheduling strategy on the process-level and thread-level parallel computing model can be accepted to improve the efficiency of the viewshed computation. And a parallel computing algorithm, XDraw-L, is designed and implemented taken into account this scheduling strategy. The experimental results demonstrate a distinct improvement of computation performance of the XDraw-L algorithm in this paper compared with the coarse-partition algorithm, like XDraw-E which is presented by Song et al. (Earth Sci Inf 10(5):511–523, 2016), and XDraw-B that is the basic algorithm of serial XDraw. Our fine-granularity scheduling algorithm can greatly improve the scheduling performance of the grid cells between the layers within a triangle region.     

一种道路网三维建模方法

针对道路网三维建模中成本高、周期长等问题,该文介绍了一种基于规则的道路网三维建模方法。在规则文件驱动下,该方法充分考虑道路网的结构化特征和地理环境,利用抽象化的道路网实体数据和属性信息,实现批量化道路网建模,并在云平台上发布,进而实现道路网三维景观的浏览、查询和分析。采用上述方法实现了一段包含立体交叉在内的道路网的三维模型的快速搭建。结果表明,该方法的规则复用性强、建模效率高、模型调整速度快,为在三维GIS环境下构建道路网模型提供了高效的技术手段。     

Identifying urban–rural differences in social vulnerability to natural hazards: a case study of China

Identifying and analyzing the urban–rural differences of social vulnerability to natural hazards is imperative to ensure that urbanization develops in a way that lessens the impacts of disasters and generate building resilient livelihoods in China. Using data from the 2000 and 2010 population censuses, this study conducted an assessment of the social vulnerability index (SVI) by applying the projection pursuit cluster model. The temporal and spatial changes of social vulnerability in urban and rural areas were then examined during China’s rapid urbanization period. An index of urban–rural differences in social vulnerability (SVI_D) was derived, and the global and local Moran’s I of the SVI_D were calculated to assess the spatial variation and association between the urban and rural SVI. In order to fully determine the impacts of urbanization in relation to social vulnerability, a spatial autoregressive model and Bivariate Moran’s I between urbanization and SVI were both calculated. The urban and rural SVI both displayed a steadily decreasing trend from 2000 to 2010, although the urban SVI was always larger than the rural SVI in the same year. In 17.5% of the prefectures, the rural SVI was larger than the urban SVI in 2000, but was smaller than the urban SVI in 2010. About 12.6% of the urban areas in the prefectures became less vulnerable than rural areas over the study period, while in more than 51.73% of the prefectures the urban–rural SVI gap decreased over the same period. The SVI_D values in all prefectures had a significantly positive spatial autocorrelation and spatial clusters were apparent. Over time, social vulnerability to natural hazards at the prefecture-level displayed a gathering–scattering pattern across China. Though a regional variation of social vulnerability developed during China’s rapid urbanization, the overall trend was for a steady reduction in social vulnerability in both urban and rural areas.

     

昌宁—孟连结合带北端癞痢头山岩体地球化学特征及其对成岩成矿作用的启示

位于昌宁—孟连结合带北端的癞痢头山岩体,主要由中粗晶花岗闪长岩(MGG)和细晶花岗闪长岩(FG)组成,形成于中三叠世(232 Ma).MGG和FG岩相学特征相似,无明显空间分带关系,呈过渡渐变关系.它们具高硅(w(SiO2) =66.2％～68.2％,平均67.0％),富钾(K2O/Na2O=1.04～2.59,平均1.64),贫碱(w(Na2O+K2O)=5.45％～7.31％,平均5.99％)、低钠(w(Na2O)=1.54％～3.58％,平均2.38％),高镁值(Mg#=47.3～58.7,平均51.5)和过铝质(A/CNK=1.41～2.17,平均1.65)特征,属过铝质高钾钙碱性系列.二者轻、重稀土元素含量比值介于6.91～12.1,表现为轻稀土元素相对富集、重稀土元素相对亏损.岩石富集w(Rb) (54.7～176×10-6),w(U)(1.10～4.77×10-6),w(Th) (9.33～35.30×10-6)等大离子亲石元素,亏损w(Ta) (1.22～1.51×10-6),w(Nb)(14.9～18.5×10-6),w(Zr)(78.4～171.0×10-6)等高场强元素,以及δEu(0.35～0.57)和"TNT"负异常特征,表明初始岩浆以陆壳成分为主并有幔源组分的加入.癞痢头山岩体样品具有较高的A/CNK值(1.41～2.17)和刚玉分子含量(1.87％～8.67％),以及Nb/Ta比值(11.03～13.76)和Zr/Hf比值(23.13～47.85)分别与平均地壳和壳源岩石相应比值一致,显S型花岗岩特征.综合分析,认为癞痢头山岩体属于S型花岗岩,起源于古老地壳物质(贫黏土的古老杂砂岩)部分熔融岩浆和底侵作用带入的幔源岩浆的混合岩浆,是保山—思茅地块陆陆后碰撞阶段力学性质由挤压向伸展转化动力学背景下的产物,具备成矿作用发生的物质基础,有一定的找矿潜力.     

中国大陆强震孕育深处的变形模式

地震预测面临着地球内部的"不可入性"和大地震"非频发性"的问题,是极具挑战性且尚待解决的世界性科学难题之一(陈运泰,2009) .地震活动是断裂活动的直接体现,精定位的震源深度客观反映了断裂深部发生地震滑动变形的起始位置.断裂深部发生的慢滑动和重复地震等对揭示断裂深部行为具有重要的指示意义(Harris,2017).在...     

海底孔内水合物观测系统与钻井实施

通过耦合布设海底表面和孔内水合物观测系统获取物理、化学、微生物等数据是了解掌握海洋水合物环境动态变化、碳循环规律和资源开发技术的有效手段,据此初步提出了一套海底孔内水合物观测系统,并指出了系统的关键技术难点。而在钻井布设该孔内观测系统时,要注意防止井内安全事故的发生,充分评估安放位置处含水合物地层的地质力学稳定性并采取适当的钻井和完井方式稳定或强化含水合物地层以支撑孔内观察仪器,实现孔内长期监测。最后讨论了在我国南海实施孔内水合物观测系统的必要性,并提出了初步想法。     

Summer fast ice evolution off Zhongshan Station,Antarctica

Based on the field data acquired in the program of fast ice observation off Zhongshan Station, Prydz Bay, East Antarctica during the austral summer 2005/ 2006, physical properties evolution of fast ice during the ice ablation season is analyzed in detail. Results show that the annual maximum ice thickness in 2005 occurred in later November, and then ice started to reek, and the ablation duration was 62 days; sea water under the ice became warmer synchronously; corresponding to the warming sea ice temperature, a ＂relative cold mid-layer＂ appeared in sea ice; the fast ice marginal line recoiled back to the shore observably, and the recoil distance was 20.9 km from 18 December 2005 through 14 January 2006. In addition, based on the data of sea ice thickness survey along the investigation course of MV Xuelong on December 18 of 2005, the ice thickness distribution paten in the marginal ice zone have been described ： sea ice thickness increased, but the diversity of floe ice thick-ness decreased from open water to fast ice zone distinctly.     

鱼类摄食生态研究的理论及方法

综合分析了鱼类摄食生态研究的理论，方法及其应用中的诸多问题：胃含物分析法，营养生境及食物重叠，食物链与食物网，营养结构与营养级，摄食器官的形态特征与摄食行为的关系，食物选择以及生理生太因子对鱼类摄食的影响；提出了食物排序的概念、优势度指数Ｄ及其判定公式，并以此比较客观地建立了摄食生态相似的鱼类群落主次食物的划分标准。     

Interdecadal Variability of the East Asian Summer Monsoon and Associated Atmospheric Circulations

Based on the National Centers for Environmental Prediction and National Center for Atmospheric Research （NCEP/NCAR） reanalysis data from 1950-1999, interdecadal variability of the East Asian Summer Monsoon （EASM） and its associated atmospheric circulations are investigated. The EASM exhibits a distinct interdecadal variation, with stronger （weaker） summer monsoon maintained from 1950-1964 （1976-1997）. In the former case, there is an enhanced Walker cell in the eastern Pacific and an anti-Walker cell in the western Pacific. The associated ascending motion resides in the central Pacific, which flows eastward and westward in the upper troposphere, descending in the eastern and western ends of the Pacific basin. At the same time, an anomalous East Asian Hadley Cell （EAHC） is found to connect the low-latitude and mid-latitude systems in East Asia, which strengthens the EASM. The descending branch of the EAHC lies in the west part of the anti-Walker cell, flowing northward in the lower troposphere and then ascending at the south of Lake Baikal （40°-50°N, 95°- 115°E） before returning to low latitudes in the upper troposphere, thus strengthening the EASM. The relationship between the EASM and SST in the eastern tropical Pacific is also discussed. A possible mechanism is proposed to link interdecadal variation of the EASM with the eastern tropical Pacific SST. A warmer sea surface temperature anomaly （SSTA） therein induces anomalous ascending motion in the eastern Pacific, resulting in a weaker Walker cell, and at the same time inducing an anomalous Walker cell in the western Pacific and an enhanced EAHC, leading to a weaker EASM. Furthermore, the interdecadal variation of summer precipitation over North China is found to be the south of Lake Baikal through enhancing and reducing strongly regulated by the velocity potential over the regional vertical motions.