Dynamic land use change simulation using cellular automata with spatially nonstationary transition rules

The dynamic relationships between land use change and its driving forces vary spatially and can be identified by geographically weighted regression (GWR). We present a novel cellular automata (GWR-CA) model that incorporates GWR-derived spatially varying relationships to simulate land use change. Our GWR-CA model is characterized by spatially nonstationary transition rules that fully address local interactions in land use change. More importantly, each driving factor in our GWR model contains effects that both promote and resist land use change. We applied GWR-CA to simulate rapid land use change in Suzhou City on the Yangtze River Delta from 2000 to 2015. The GWR coefficients were visualized to highlight their spatial patterns and local variation, which are closely associated with their effects on land use change. The transition rules indicate low land conversion potential in the city’s center and outer suburbs, but higher land conversion potential in the inner near suburbs along the belt expressway. Residual statistics show that GWR fits the input data better than logistic regression (LR). Compared with an LR-based CA model, GWR-CA improves overall accuracy by 4.1% and captures 5.5% more urban growth, suggesting that GWR-CA may be superior in modeling land use change. Our results demonstrate that the GWR-CA model is effective in capturing spatially varying land transition rules to produce more realistic results, and is suitable for simulating land use change and urban expansion in rapidly urbanizing regions.     

东北地区冬半年积雪与气温对冻土的影响

利用东北地区121个气象站逐日冻土深度、积雪深度、平均气温、地表平均气温及降水量数据,分析了1964—2017年冬半年冻土的变化特征及气象要素对冻土的影响。结果表明：东北地区积雪深度、平均气温、地表平均气温与冻土深度相关系数较高,降水量相关性不大。20世纪60年代平均气温、地表平均气温及负积温最低,最大冻土深度为历年代最深;随着气候变暖,最大冻土深度以6.15 cm?（10a）^-1的速率显著减小。冬半年平均最大冻土深度为123 cm,呈显著纬向分布,自辽东半岛向大兴安岭北部递增;随纬度和海拔高度的增加,平均气温和地表平均气温降低,负积温增加,且由北向南地气温差增大。最大冻土深度全区有90%以上的站点减少,减少速率以0.1~10 cm?（10a）^-1为主。冻土持续时间随纬度升高而增加,月最大冻土深度和积雪深度最大值分别出现在3月和1月,最大冻土深度的增加要滞后于积雪深度的增加。由于积雪对地温的保温作用,积雪深度较浅时,冻土深度增加较明显,随着积雪深度的增加,冻土深度变化较小,积雪对冻土起到了保温的作用。对于高纬度地区站点,30 cm左右为积雪的保温界限值;对于沿海站点,积雪保温的界限值在5 cm左右;在相同地形下,冻土深度较浅区域积雪的保温值因海拔高度、气候特点而异。最大冻土深度对地表平均气温升温的响应更为显著,地表平均气温和平均气温每升高1 ℃,最大冻土深度将减小8.4 cm和10.6 cm,负积温每减少100 ℃?d,最大冻土深度减少4.9 cm。     

The Impact of Climatic and Environmental Factors on n-Alkanes Indices in Southwestern Tibetan Plateau

n-Alkanes are widely used in paleoenvironmental reconstructions.However,our understanding of changes in the distribution of n-alkanes with climatic and environmental factors remains unclear in arid/semi-arid regions.We sampled 26 surface sediments from three climatic zones across the southwestern Tibetan Plateau to evaluate the sensitivity of chain length distributions of n-alkanes to climatic and environmental parameters.Our observations demonstrate that average chain length (ACL), proportion of aquatic macrophyte (Paq),carbon preference index (CPI) and ratio of the contents of nC27 and nC31 (nC27/nC31) are all sensitive to hydroclimatic conditions.In contrast to commonly-adopted assumptions,the correlations between these indices and hydrological parameters are not always good,which indicates that the interpretation of n-alkane indices is special on the southwestern Tibetan Plateau.These might be related to the vegetation characteristics and seasonality of biological activity,and need to be considered in paleoclimatic reconstruction.The impact of seasonal precipitation on n-alkanes indices was also evaluated.     

普朗斑岩型铜矿床成矿元素多重分形特征及其矿化强度指示

成矿元素品位的变化受控于长期地质作用过程,具有较强的非线性特征,定量刻画其变化特征有助于深入理解成矿过程,并为找矿预测提供理论依据。运用多重分形去趋势移动均值（MFDMA） 法,分析云南普朗斑岩型铜矿床4号勘探线钻孔Cu元素品位的多重分形特征及其对矿化强度的指示意义。结果显示：所有钻孔Cu元素品位具有多重分形特征,且不同矿化等级的元素分布存在局部奇异性差异,多重分形强度随矿化强度等级减弱而增加;这些特征指示在钾硅化带—带内侧强矿化钻孔的高品位相对聚集,而带外侧角岩化—青磐岩化带的高品位相对分散;通过数据随机重排技术推测Cu品位多重分形结构奇异性差异是由其分布概率密度和长相关性共同作用引起,后者的作用贡献大于前者。     

辽宁省沈抚新区区域地壳稳定性评价

在广泛收集资料和野外调查工作的基础上,经综合分析控制和影响区域地壳稳定性的主要因素和内外动力地质的耦合作用,基于ArcGIS平台采用多因素加权叠加分析方法开展了沈抚新区区域地壳稳定性评价研究。评价因素主要选取活动断裂、地震活动性和深部地球物理以及工程岩组、地形地貌与地表地质灾害等共12个影响因子,并依据评价因素的关联性和重要性进行分类赋值,进而建立了区域地壳稳定性评价模型。对沈抚新区的区域地壳稳定性进行了定量化评价。结果表明,沈抚新区以稳定区和次稳定区为主,面积为728.9 km2,占总面积的81.07%,有利于规划和建设。建议在地下工程穿越不稳定区和次不稳定区时,要采取必要的工程措施,确保基础和边坡的稳定性。规划的地铁9号线主要穿越次不稳定区和次稳定区。施工带在埋深15~25 m之间,该段地层主要为卵砾石堆积层、泥包砾堆积层与变质片麻岩强风化带及其交接部位,施工带处于地下水浸没带或变动带,其支护比较困难,建议尽早确定合理的施工方案,保障施工顺利和安全,此外还需要采取一定降水防腐措施。      

中国地区集合预报产品自适应递减平均偏差订正法的改进研究

基于中国地区T213集合预报产品2 m温度预报数据,采用卡尔曼滤波类型的自适应递减平均法进行偏差订正处理,原方案在剧烈降温天气订正效果表现不理想。通过对递减平均参数w的重新构建得到改进的订正方案w(i,p)(i为站点信息,p为天气过程信息),在此基础上进一步优化对历史信息的有效提取,得到改进的方案w(i,p)相似法和w(i,p)统计法,并进行效果检验。结果表明:改进为包含空间和天气过程信息的函数w(i,p)后方案的订正效果得到不同程度的提高,其中24 h剧烈降温预报各成员预报均方根误差平均减小了0. 15℃;而进一步改进的w(i,p)统计法在当前几种剧烈降温预报中订正效果最优,其集合平均偏差与w(i,p)方案相比减小2. 54℃。     

GRAPES区域集合预报两种侧边界扰动方法对比试验

中国气象局数值预报中心自2014年建立了区域集合预报业务系统,其使用的侧边界扰动由全球集合预报系统动力降尺度得到。为深入了解侧边界扰动对区域集合预报的影响,基于15 km水平分辨率的区域集合预报模式,使用动力降尺度方法和尺度化滞后平均法(scaled lagged average forecasting,SLAF)设计构造了两种侧边界扰动方案,并开展了2015年7月共6天的集合预报试验,利用集合均方根误差、集合离散度、连续分级概率评分、离群值、Brier Score及相对作用特征曲线面积等概率预报检验方法进行了多方面检验,分析了两种侧边界扰动方案对区域集合预报质量的影响。结果表明:动力降尺度侧边界扰动方案(DOWN)的扰动总能量在各垂直层次均大于SLAF方案,使得边界上前者的离散度大于后者,集合扰动增长更为合理;对于等压面要素和地面要素,DOWN方案的离散度、Outlier、CRPS等评分优于SLAF方案,反映了DOWN方案构造的侧边界扰动更加合理;在降水概率预报技巧方面,SLAF方案在评分上具有一定优势,但评分的提高没有通过显著性水平检验,因此认为两种方案对降水预报的改进基本相当。     

四川地形扰动对降水分布影响

引入一维加权平均的谱分析方法定量研究四川地形强迫对该区域降水分布的影响。结果表明:纬向地形和冬季降水谱峰锁相于同一波长（475.8 km）,呈共振关系,地形与其他季节降水呈漂移关系,这与经向和纬向上环流变动有关,即冬季纬向环流占主导,纬向地形触发的大气波动对冬季降水策动作用大;夏季降水是各种不同尺度系统相互作用的结果,地形是重要因素之一。经向和纬向地形特征尺度分别为296.8 km和475.8 km,反映了地形强迫的中尺度特征,且纬向地形谱峰比经向大1个数量级,纬向强迫更明显。夏季降水谱峰比冬季大2个数量级,降水系统纬向特征尺度比冬季小约150 km,说明夏季在纬向地形强迫下,降水系统尺度减小的同时其强度大大增加,这在一定程度上可以解释中尺度对流性降水在夏季偏多。四川夏季最大降水位于雅安地区,其地形扰动比四川整体扰动更明显,故产生的降水也更大。夏季降水和经向地形锁相于同一波长（37.1 km）,经向地形对雅安夏季强降水起关键作用。     

Energy average formula of photon gas rederived by using the generalised Hermann—Feynman theorem

By virtue of the generalised Hermann--Feynmam theorem we re-derive the energy average formula of photon gas. This is another useful application of the theorem.     

Application of the weighted residual method to diffraction by 2-D canyons of arbitrary shape: I. Incident SH waves

A solution for the two-dimensional scattering and diffraction of plane SH waves by canyons of arbitrary shape in an elastic half space is presented. The wave field for arbitrary geometry in this paper is computed numerically by the method of weighted residues (moment method). The wave displacement field computed by the present residual method for the case of a semi-circular canyon was shown to agree analytically and numerically with that computed by the exact closed form series solution. The same observations about ground amplifications, their dependence on frequencies and orientations of the incident waves, can be stated here for canyons of arbitrary shape as previously made for circular canyons.