Modeling spatially non-stationary land use/cover change in the lower Connecticut River Basin by combining geographically weighted logistic regression and the CA-Markov model   总被引：2，自引：0，他引：2  

Hui Wang  Scott R. Stephenson  Shijin Qu 《International journal of geographical information science》2019,33(7):1313-1334

Simulating land use/cover change (LUCC) and determining its transition rules have been a focus of research for several decades. Previous studies used ordinary logistic regression (OLR) to determine transition rules in cellular automata (CA) modeling of LUCC, which often neglected the spatially non-stationary relationships between driving factors and land use/cover categories. We use an integrated geographically weighted logistic regression (GWLR) CA-Markov method to simulate LUCC from 2001–2011 over 29 towns in the Connecticut River Basin. Results are compared with those obtained from the OLR-CA-Markov method, and the sensitivity of LUCC simulated by the GWLR-CA-Markov method to the spatial non-stationarity-based suitability map is investigated. Analysis of residuals indicates better goodness of fit in model calibration for geographically weighted regression (GWR) than OLR. Coefficients of driving factors indicate that GWLR outperforms OLR in depicting the local suitability of land use/cover categories. Kappa statistics of the simulated maps indicate high agreement with observed land use/cover for both OLR-CA-Markov and GWLR-CA-Markov methods. Similarity in simulation accuracy between the methods suggests that the sensitivity of simulated LUCC to suitability inputs is low with respect to spatial non-stationarity. Therefore, this study provides critical insight on the role of spatial non-stationarity throughout the process of LUCC simulation.  相似文献   

全球地磁感应测深数据三维反演   总被引：2，自引：0，他引：2       下载免费PDF全文

李世文  翁爱华  张艳辉  李建平  杨悦  唐裕  邹宗霖  李春成 《地球物理学报》2019,62(5):1908-1920

全球地磁感应测深能获得地幔转换带及下地幔上部的导电结构.但目前稀疏的地磁台站分布及部分台站的观测数据稳定性较差,影响了三维反演对地下电性结构的分辨力和反演可靠性.为此,区别于传统的L2-范数反演方法,本文提出并实现了基于L1-范数的地磁测深响应三维反演技术.在反演中,利用L1-范数度量数据预测误差,降低"飞点"数据的影响,将相关系数较小的C-响应估计也纳入反演数据中.三维正演模拟采用球坐标系下的交错网格有限差分法,反演采用有限内存拟牛顿法.文中利用指数概率密度分布函数构造非高斯噪声的合成数据,并采用棋盘模型对反演方法的可靠性进行了验证.之后,我们将本文提出的三维反演方法用于全球129个地磁观测台站的C-响应数据反演,结果表明在地幔转换带深部,中国东北地区为高导电异常,南欧和北非则均为高阻;夏威夷在900km以下为高导;菲律宾海及以东地区在转换带表现为明显的高阻,这些结果与前人研究结果一致.由于采用了更多的台站数据,我们的反演结果还发现一些新的异常:南美洲南端,转换带表现为明显的高导;澳大利亚东南部,地幔转换带深部,也存在一个明显的高导异常,这些异常分布和地震层析成像的低速区一致.因此,L1-范数三维反演能够充分利用全球C-响应数据信息,提高地磁测深对地球深部电性结构的分辨能力,更好的研究全球地幔电性结构.  相似文献   

VIIRS/DNB夜间灯光月度产品插补方法对比——以北京为例     

陈慕琳  蔡红艳 《地理科学进展》2019,38(1):126-138

新一代VIIRS/DNB(Visible Infrared Imaging Radiometer Suite Day/Night Band)夜间灯光数据因其具备更精细的时空分辨率、数据不存在饱和现象、不同年份数据可比性强等优点而迅速代替DMSP/OLS(Defense Meteorological Satellite Program's Operational Linescan System)夜间灯光数据成为新的研究热点。但由于杂散光的污染,VIIRS/DNB夏季数据数值缺失严重,致使数据在空间和时间上不连续,因此,插补缺失数据成为后续应用的前提。鉴于此,论文从插补结果异常值、与参照值对比、计算性能三个方面,系统比较了三次样条插值(样条)、三次Hermite插值(Hermite)、灰色预测模型(GM)、三次指数平滑法(指数)4种插补方法的适用性,以期为插补VIIRS/DNB夜间灯光数据提供方法选择的依据。研究结果表明：①异常值比较方面,Hermite法未出现异常值,另外3种算法仅出现少量异常值(0.02%~1.34%);②与参照值的对比方面,Hermite法与参考值接近程度最高,GM接近程度最低,样条法和指数法介于两者中间;③算法性能比较方面,4种方法都具备计算简单、容易编程的特点,但指数法的算法时长是另外3种方法的10倍以上。因此,综合多方表现,当插补月份前后两侧均有足够长的原始数据时,插补效果好、计算速度快、不会出现过冲现象的Hermite法最适宜,样条法次之;当插补月份仅单侧有足够长的数据时,适宜采用指数法(插补效果好、计算速度较慢)或GM(插补效果偏低、计算速度快)进行插补。  相似文献   

频率域海洋可控源电磁垂直各向异性三维反演   总被引：1，自引：0，他引：1       下载免费PDF全文

彭荣华  胡祥云  李建慧  刘亚军 《地球物理学报》2019,62(6):2165-2175

地层宏观电性各向异性会对可控源电磁响应产生重要影响.由于海底地层电性结构常表现为电导率各向异性，若仅对海洋可控源电磁（MCSEM）数据进行常规各向同性反演，有可能无法获得准确的反演解释结果，从而削弱MCSEM技术的可靠性.本文实现了电导率垂直各向异性（VTI）条件下频率域海洋可控源电磁数据三维反演算法.其中，三维正演采用基于二次场控制方程的交错网格有限体积法，并利用直接矩阵分解技术来求解离散所得的大型线性方程组，有利于快速计算多场源的响应.反演采用具有近似二次收敛性的高斯牛顿算法对目标函数进行最优化.最后，对具有VTI电性各向异性特征的盐丘构造模型的MCSEM合成数据分别进行了电导率各向同性和垂直各向异性三维反演，结果表明：各向同性三维反演算法无法对受VTI介质影响的MCSEM数据进行正确的反演解释，而垂直各向异性三维反演能够获得更为可靠的地下电阻率结构和异常体分布，展现出对海底电性各向异性结构更为优良的反演解释能力.  相似文献   

球坐标系下三维大地电磁正演研究   总被引：1，自引：0，他引：1       下载免费PDF全文

罗威  王绪本  王堃鹏  张刚  李德伟 《地球物理学报》2019,62(10):3885-3897

大地电磁正演理论研究热点一直以来主要集中在如何提高计算效率和精度，但在剖面足够长、探测深度足够大的情况下，传统的笛卡尔坐标系数值模拟方式难以准确拟合地球曲率形态.本文研究了基于球坐标系的三维大地电磁正演，推导了交错网格有限差分三维正演公式，与一维解析解和三维标准模型测试对比，验证了正演算法的正确性.通过理论模型计算，对比分析球坐标和笛卡尔坐标系正演结果表明：球坐标系模拟更合理，避免了传统笛卡尔坐标拉伸投影所引入的误差，可代替目前的笛卡尔坐标模拟方法.基于球坐标和笛卡尔坐标系的三维大地电磁正演响应值随着频率变低差异越明显.球坐标和笛卡尔坐标计算结果差异度与频率、模型结构和电阻率有关.本文模型计算结果在数万秒周期处已出现接近10%的差异，对于较大尺度的长周期大地电磁，地球曲率的影响不能忽略.  相似文献   

基于稀疏约束反演的三维混采数据分离          下载免费PDF全文

魏亚杰  张盼  许卓 《地球物理学报》2019,62(10):4000-4009

混合震源采集技术相对于传统的地震数据采集，在极大提高采集效率的同时引入了混叠噪声，很大程度上影响了成像结果的精度.二维混采数据中，我们通常利用混叠噪声在非共炮域呈非相干分布这一特点来压制混叠噪声，从而实现混合震源数据分离.相对于二维混采数据，三维混采数据具有数据量巨大，构建混合震源算子困难，混合度的增加引入了高强度混叠噪声的特点.针对上述问题，本文采用稀疏约束反演方法在Radon域实现混采数据分离，混叠噪声强度比较大的情况下，稀疏约束反演方法能够得到更高精度的分离结果；利用震源激发的GPS时间通过长记录的方式在共接收点道集对上一次迭代分离结果做混合、伪分离，实现了单个共接收点道集自身混合、伪分离，避免了对整个数据做运算，同时不需要构建混合震源算子.通过模拟数据和实际数据计算来验证上述方法的适用性.  相似文献   

非迁移潮作用下的电离层赤道异常区单峰现象及其经度分布          下载免费PDF全文

田耀宇  郝永强  张东和  肖佐 《地球物理学报》2019,62(11):4067-4081

当喷泉效应较弱而双峰结构发展不充分的时候，可能在赤道异常区仅能够观测到一个电子密度的峰值，称之为单峰现象.本文利用CHAllenging Minisatellite Payload卫星在2001-2010年的电子密度数据给出了单峰的发生规律，单峰在地方时早上08：00-10：00和下午16：00-19：00发生率高，发生位置在经度上呈现多波数分布，尤其在10：00-18：00明显：在分季时多呈现四波，而在冬至季时以三波为主.单峰发生多的经度，正好对应着双峰的结构特征较弱之处.究其原因，是非迁移潮的DE2和DE3分量调制了背景风场和大气发电机电场，在电场和喷泉效应减弱的经度，双峰结构难以形成时，就会表现为单峰结构.本文扩展了对单峰现象的地方时、季节和经度分布等规律的了解，明确了非迁移潮在其中施加的影响，由此，单峰同双峰现象一样可以用于研究非迁移潮对热层-电离层的作用.  相似文献   

GPS/PWV时间序列特征提取方法的研究     

胡广保  叶世榕  张彦祥  夏朋飞  夏凤雨 《大地测量与地球动力学》2019,39(1):36-40

针对GAMIT/GLOBK软件解算得到的4 a GPS/PWV时间序列的特征提取问题，提出基于滤波辅助的部分集成经验模态分解（PEEMD）与Hilbert谱分析相结合的特征提取方法。首先，在PEEMD方法的基础上,结合滤波辅助的PEEMD方法与Hilbert谱分析，建立GPS/PWV时间序列特征提取模型；然后，将所提出的方法应用于TNML测站4 a的GPS/PWV长时间序列和7 d的GPS/PWV短时间序列分析中，并将滤波辅助的PEEMD结果与传统的小波分解结果进行对比。结果表明，该特征提取方法能准确有效地提取出 GPS/PWV时间序列中的周年周期和日周期特征分量，滤波辅助的PEEMD分解结果与小波分解结果一致，且提取的特征分量与原始信号更加吻合。  相似文献   

A stochastic method for the prediction of icebreaker bow extreme stresses     

《Applied Ocean Research》2019

It is well established that the ship-ice interaction process is quite complex and associated ice loads on the icebreaker hull is a stochastic process. Obviously, novel accurate statistical methods and models should be developed and applied to estimate extreme bow stresses.This paper studies icebreaker bow stresses based on measured distribution of ice thickness in the Arctic Ocean on the way to and from the North Pole. Since the vessel route was carefully selected searching for easier ice conditions, the Arctic Ocean crossing was not a straight linear but a meandering path. Thus, the specific ship route data was biased with respect to general ice statistics in the region, but true with respect to the route specific ice data encountered by a ship navigating in that region. Therefore the route specific ice thickness data is directly needed for ship design and navigation analysis. It is assumed that captains are competent and knowledgeable, and therefore will select a route that provides the most favourable ice conditions.This paper contributes to study of the newest Chinese self-designed polar icebreaker, serving the purpose of enhancing icebreaker operational reliability. Finite Element Method software package ANSYS/LS-DYNA has been employed to simulate bow stress pattern for a particular icebreaker operating in the Arctic Ocean. Extreme bow stresses were estimated using Naess-Gaidai method. The latter is a first application of Naess-Gaidai method to a distribution with lower bound. Thus this paper aims at introducing an efficient method of estimating route-specific icebreaker extreme bow stresses.  相似文献   

Evaluation of aquatic ecological systems through dynamics of ichthyofaunal diversity in a Himalayan torrential river,Murti     

《Limnologica》2020

The biodiversity hotspot region of the Eastern Himalayas consists of a vast freshwater network enriched with species diversity. Many small-scale torrential rivers and water reaches contribute to the species pool of all the major rivers by converging downstream. These reaches are most likely to be degraded at a faster rate as compared to the large-scale rivers following an increased rate of urbanization, habitat alterations, and changing climatic conditions. Therefore, this study aims to explore River Murti, which is a representative small scale river system characterized by a large altitudinal gradient and a diverse watershed area. Ichthyofaunal diversity (i.e., diversity, evenness & richness) and 21 environmental variables are measured through a tri-seasonal sampling effort conducted along 14 selected locations. A total of 41 fish species (including species belonging to 4 Near Threatened, 8 Vulnerable, and 1 Endangered) are found inhabiting this river. Ichthyofaunal assemblage is found to be primarily modulated by habitat diversity and landscape variables. Three Aquatic Ecological Systems (AES) have been identified along this river in a top-down approach based on recorded environmental variables. We have calculated an observed/expected ratio for each diversity indices along 14 locations based on predicted temporal variability using boosted regression (BRT) models. The evaluation of diversity status has been kept at 0.5 to account for a 50% loss or deviation from observed (O/E50). This evaluation has been successfully used to delineate AES1 with majorly “Impaired” status and thus ensures its importance in terms of species conservation. Our study indicates the contribution of 11 major environmental drivers modulating the species assemblage patterns in these AES. Amongst them, altitude, substrate coarseness, river morphology, and shelter availability are strongly associated with species diversity as per the BRT models. These underlying factors are also correlated with “basin pressure,” suggesting that anthropogenic disturbances, as well as the changing climate, might play an important role in the gradual change in environmental conditions, which in turn could cause a shift in species assemblage structure.  相似文献