Modeling urban growth in a metropolitan area based on bidirectional flows,an improved gravitational field model,and partitioned cellular automata

Simulating urban landscape dynamics in metropolitan areas has attracted much attention lately, but the difficulty remains. Although large-scale urban simulation studies consider spatial interaction as an important factor, spatial interaction cannot be accurately measured based on a single element flow, and its effects may not strictly follow a distance decay function. Furthermore, different cities may require different transition rules. In this study, we combined bidirectional flows of population and information and an improved gravitational field model to model the urban spatial interaction, and we then integrated a partitioned cellular automata (CA) model to simulate the urban growth for different cities in the Yangtze River middle reaches megalopolis. It was found that the simulation results generated by the CA model considering spatial interaction are significantly improved. Furthermore, partitioned conversion thresholds can effectively improve the model performance. The proposed model showed a much better performance in the simulation of subordinate cities surrounding the core cities, than for the core cities and fringe cities. We suggest that large-scale urban simulation should pay more attention to the development of partitioned transition rules. The effects of intercity urban flows should also be considered in the simulation of small- and medium-sized cities near the regional cores.     

台风影响下长江口羽流的响应与淡水输运机制

河口区域冲淡水锋面及物质输运会对台风产生快速而复杂的响应,并产生显著的生态效应及沉积过程,且不同类型的台风会对水动力及物质输运产生不同影响。本文以长江河口为研究区域,利用盐度锋面附近的浮标站点观测数据,发现在登陆型超强台风“利奇马”的影响下,旋转流特征消失,站点全水深出现持续近两天的北向流;而在转向型强台风“巴威”的影响下,出现持续约三天的全水深南向流。应用覆盖长江口及邻近海域的有限体积海洋模型（FVCOM）对两次台风过程进行模拟,实验对比台风过境前后流场、盐度场并且计算淡水通量,发现台风过境时通过破坏夏季典型的向海扩张的表层平流型羽流结构,从而促进盐度垂向混合,导致大量淡水堆积在靠岸一侧,加强了沿岸的淡水输运,进而形成底部捕获型羽流。对于登陆型台风“利奇马”的影响,淡水输运主要沿苏北海岸向北,而在转向型台风“巴威”偏北风的影响下,淡水则显著沿岸向浙闽海域流动。这两种台风过后,表层淡水开始向海扩展,垂向盐度分层再一次呈现,在2～3天内完全恢复为表层平流型羽流结构。     

A minimum-volume oriented bounding box strategy for improving the performance of urban cellular automata based on vectorization and parallel computing technology

ABSTRACT

As an effective tool for simulating spatiotemporal urban processes in the real world, urban cellular automata (CA) models involve multiple data layers and complicated calibration algorithms, which make their computational capability become a bottleneck. Numerous approaches and techniques have been applied to the development of high-performance urban CA models, among which the integration of vectorization and parallel computing has broad application prospects due to its powerful computational ability and scalability. Unfortunately, this hybrid algorithm becomes inefficient when the axis-aligned bounding box (AABB) of study areas contains many unavailable cells. This paper presents a minimum-volume oriented bounding box (OBB) strategy to solve the above problem. Specifically, geometric transformation (i.e. translation and rotation) is applied to find the OBB of the study area before implementing the hybrid algorithm, and a set of functions are established to describe the spatial coordinate relationship between the AABB and OBB layers. Experiments conducted in this study demonstrate that the OBB strategy can further reduce the computational time of urban CA models after vectorization and parallelism. For example, when the cell size is 15 m and the neighborhood size is 3 × 3, an approximately 10-fold speedup in computational time can result from vectorization in the MATLAB environment, followed by an 18-fold speedup after implementing parallel computing in a quad-core processor and, finally, a speedup of 25-fold by further using an OBB strategy. We thus argue that OBB strategy can make the integration of vectorization and parallel computing more efficient and may provide scalable solutions for significantly improving the applicability of urban CA models.     

Jurassic Subduction of the Paleo-Pacific Ocean in NE China: Zircon U–Pb and Sr-Nd-Hf isotopic Evidence from the Huashan and Taili Monzogranites

The Qilian orogen along the NE edge of the Tibet‐Qinghai Plateau records the evolution of Proto‐Tethyan Ocean that closed through subduction along the southern margin of the North China block during the Early Paleozoic. The South Qilian belt is the southern unit of this orogen and dominated by Cambrian‐Ordovician volcano‐sedimentary rocks and Neoproteozoic Hualong complex that contains similar rock assemblages of the Central Qilian block. Our recent geological mapping and petrologic results demonstrate that volcano‐sedimentary rocks show typical rock assembles of a Cambrian‐early Ordovician arc‐trench system in Lajishan Mts. along the northern margin of the Hualong Complex. Island arc rocks including basalt, andesite, dacite, rhyolite, and breccia is in fault contact with ophiolite complex consisting of mantle peridotite, serpentinite, gabbro, dolerite, plagiogranite, and basalt. Accretionary complexes are tectonically separated from the ophiolite‐arc rocks, with various rock assemblages spatially. They consist of pillow basalt, basalt breccia, tuff, chert, and limestone blocks with a seamount origin within the scaly shale in Dingmaoshan and Donggoumeikuang areas, and basalt, chert, and sandstone blocks within muddy shale matrix and mélange at Lajishankou area. Abundant radiolarians occur in red chert, and trilobite, brachiopod, and coral fossils occur within Dingmaoshan limestone blocks. Although partial basalt or chert blocks are highly disrupted, duplex, thrust fault, rootless intrafolial fold, tight fold, and penetrative foliation are well‐developed at Donggoumeikuang area. Spatially, accretionary complexes lie structurally beneath ophiolite complex and above the turbidites of the Central Qilian block. Ophiolite and accretionary complexes are also overlapped by late Ordovician molasse deposits sourced from Cambrian arc‐trench system and the Central Qilian block. These observations demonstrate that a Cambrian‐early Ordovician trench‐arc system within the South Qilian belt formed during the early Paleozoic southward subduction of the South Qilian Ocean collided with the Central Qilian block prior to the late Ordovician.     

利用正演模拟实现面波衰减

瑞雷波是分布在自由表面附近的P-SV面波,具有强振幅、低频、频散的特征,常用于工程物探,而在油气勘探及深部探测中常作为一种"干扰"波,会影响地震资料的信噪比,需要将面波净化处理。为实现面波衰减,本文将工程物探方法和石油物探技术相结合,通过面波频散曲线反演获得高精度的近地表模型,在三维近地表横波速度模型基础上进行面波正演,利用面波模型自适应匹配相减的方法将其减去,在不伤及有效波的情况下将面波的影响降至最低,一定程度上避免了常规面波衰减方法伤及有效波、损失低频弱信息、导致资料有效频带变窄、对地震构造解释及反演造成负面影响等缺点。最后应用实际数据验证了本研究方法的有效性。     

岩石圈和上地幔的热—化学结构概率反演——以LitMod1D为例

本文介绍了一种一维联合多种地球物理数据的概率反演方法 Lit Mod1D的基本原理及应用。该联合反演方法结合了地球物理学与岩石学的正演模拟,将岩石学、矿物物理学、地球化学、以及地球物理学等观测数据限定在一个自适应的热—动力学框架下,通过解热传导、热力学、重力势、流变学和地壳均衡等方程,获得岩石圈及上地幔深部(地表到410 km不连续面)的热、组分及地震波速度等结构。Lit Mod1D基于贝叶斯框架进行反演,使用马尔可夫链蒙特卡洛(Markov Chain Monte Carlo,MCMC)随机搜索方法,最终得到的岩石圈及上地幔深部的热—化学结构不仅能够很好的拟合各种地球物理数据,同时,由于不同的地球物理观测数据对于不同深度的组分或/和热结构异常有着不同的敏感度,因此相比于单独反演降低了反演的不确定性。我们选取了华北克拉通的阳原和狼山测点作为介绍该方法的应用实例,其中克拉通东部的阳原测点之下显示了一个热的、薄的岩石圈地幔,而西部的狼山测点则显示了一个冷的、相对较厚的岩石圈地幔,对应了华北克拉通岩石圈地幔减薄的空间变化。反演的Mg#值与相应地区的捕虏体有很好的一致性,两处均不存在亏损岩石圈地幔,说明其岩石圈地幔已被新生饱满地幔置换了。     

基于MODIS数据地表反照率时空变化特征及影响因子研究

地表反照率表征地球表面对太阳辐射的反射能力,是影响地表辐射能量收支平衡的关键参数。本文以淮河流域为例,利用MODIS（MODerate resolution Imaging Spectroradiometer）数据,采用网格趋势分析、异常变化分析、相关分析和灰色关联度分析等方法,分析了淮河流域2005～2015年地表反照率的时空变化规律,以及土地利用类型、地形因子、地表参数和气候等影响因子。结果表明:淮河流域年平均地表反照率整体呈“北高南低、东高西低”的空间分布规律,变化在0.043～0.223,平均值为 0.145。低值区主要集中于水体密集和山区丘陵地带,且标准差相对较小;高值区主要集中于流域中部及东部平原地带,且标准差较大。61.5%的区域地表反照率呈增加趋势,且存在季节性差异,夏季平均地表反照率最大,春季次之,秋季最小,冬季则由于降雪覆盖和农田利用的影响波动幅度较大。淮河流域地表反照率与归一化植被指数（Normalized Difference Vegetation Index,NDVI）、地表温度、气温和降水在大部分区域呈正相关,面积占比分别达到90.23%、82.32%、85.41%和93.70%。灰色关联度分析表明,不同土地利用类型（水体除外）下年均地表反照率受各因子影响排序为:NDVI>气温>地表温度>降水,空间变化受各因子影响排序为:NDVI>降水>地表温度>气温>高程。     

基于GIS的岩溶区土地利用的地形梯度分布研究——以酉阳县板溪乡为例

基于 RS 和 GIS 技术,采用地形位指数、分布指数和土地利用程度综合指数,探讨了2002、2011两年重庆市酉阳县板溪乡不同地形梯度土地利用结构的变化特征和定量分析了土地利用程度。结果显示：建设用地、耕地和水域对地形条件的选择性较强;建设用地主要聚集在海拔〈850 m、坡度〈15°176;的区域,耕地主要聚集在海拔〈1050 m 和坡度〈25°176;的区域,5和7级地形位为制约水域分布的地形条件;2011年,灌木林地在海拔650-1250 m、坡度〉15°176;的区域集中增加。总体上,该研究区的土地利用程度较高,但仍然存在耕地分布不合理的现象。