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

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.     

An efficient and robust approach to predict ship self-propulsion coefficients

Achieving a reliable and accurate numerical prediction of the self-propulsion performance of a ship is still an open problem that poses some relevant issues. Several CFD methods, ranging from boundary element methods (BEM) to higher-fidelity viscous Reynolds averaged Navier–Stokes (RANS) based solvers, can be used to accurately analyze the separate problems, i.e. the open water propeller and the hull calm water resistance. However, when the fully-coupled self-propulsion problem is considered, i.e. the hull advancing at uniform speed propelled by its own propulsion system, several complexities rise up. Typical flow simplifications adopted to speed-up the simulations of the single analysis (hull and propeller separately) lose their validity requiring a more complex solver to tackle the fully-coupled problem. The complexity rises up further when considering a maneuver condition. This aspect increases the computational burden and, consequently, the required time which becomes prohibitive in a preliminary ship design stage.The majority of the simplified methods proposed in literature to include propeller effects, without directly solve the propeller flow, in a high-fidelity viscous solver are not able to provide all the commonly required self-propulsion coefficients. In this work, a new method to enrich the results from a body force based approach is proposed and investigated, with the aim to reduce as much as possible the computational burden without losing any useful result. This procedure is tested for validation on the KCS hull form in self-propulsion and maneuver conditions.     

A stochastic method for the prediction of icebreaker bow extreme stresses

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.     

The Cabo Frio upwelling overrides geographical patterns in the population dynamics of the shrimp Artemesia longinaris Spence Bate, 1888 (Decapoda: Penaeidae)

Upwelling areas are among the most productive ecosystems on the planet, influencing the biology of marine organisms. This study investigated the population dynamics of the shrimp Artemesia longinaris in two regions in southeastern Brazil, one inside (Macaé—Rio de Janeiro State) and one outside (Ubatuba—State of São Paulo) the Cabo Frio upwelling area. The aim was to verify the influence of the upwelling phenomenon on the abundance, growth, longevity, size of sexual maturity, and reproductive period of the species. In total, 188,902 individuals were captured at Macaé and 3,461 at Ubatuba. Individuals captured at Macaé showed larger maximum size, higher longevity, and slower growth rate, besides reaching sexual maturity at larger sizes than at Ubatuba. Continuous reproduction was observed in both regions, with juvenile recruitment peaks in spring and summer. Local conditions observed at Macaé were influenced by the Cabo Frio upwelling zone, characterized by productive and cooler waters that are around 20°C during most of the year. The upwelling phenomenon is probably the main factor influencing the population parameters studied here, changing the geographic patterns previously observed for the variation of these parameters in A. longinaris.     

再探温跃层有孔虫B/Ca与温度和[CO32-]的依赖性及其在24 ka以来热带西太平洋次表层碳酸盐系统重建中的应用

The B/Ca ratio of planktonic foraminifer shells has been used as a proxy for reconstructing past ocean carbonate chemistry. However, recent studies have revealed significant uncertainties associated with this proxy, such as whether seawater temperature or [ CO_3~(2-)] is the dominant control on the partition coefficient(K_D) of planktonic foraminiferal B/Ca. To address these uncertainties and thus improve our understanding of the planktonic foraminiferal B/Ca proxy, we analysed B/Ca ratios in the tests of Neogloboquadrina dutertrei(300–355 μm) and Pulleniatina obliquiloculata(355–400 μm) in surface sediment samples from the tropical western Pacific and South China Sea. The relationship between these B/Ca ratios and bottom water calcite saturation states(Δ[ CO_3~(2-)]) is weak, thus suggesting only a small dissolution effect on the B/Ca of the two species. The correlation coefficients(R~2) between the B/Ca ratios of N. dutertrei and P. obliquiloculata and environmental parameters(e.g., temperature, salinity, phosphate, DIC and ALK) in the tropical western Pacific and South China Sea are not high enough to justify using B/Ca ratios as a palaeoenvironmental proxy in the study areas. The significant correlation between K_D values of N. dutertrei and P. obliquiloculata and carbonate system parameters(e.g.,[ CO_3~(2-)], DIC, ALK, pH and [ HCO_3~-]) in the study area reflect chemical links between the K_D denominator and these variables. Based on our surface sediment calibration, an empirical relationship between the K_D of N.dutertrei and temperature is proposed in the tropical western Pacific. We also generated a record of B/Ca ratios in N. dutertrei(300–355 μm) from Core MD06-3052 in the tropical western Pacific over the past 24 ka to evaluate the application of the revised B/Ca proxy method. Based on the reconstructed empirical relationship for B/Ca and subsurface seawater ALK, we estimated subsurface seawater carbonate system parameters in the tropical western Pacific since 24 ka. In general, the estimated subsurface seawater pH and [ CO_3~(2-)] show an increase with time, and the record of subsurface seawater pCO_2 shows a decrease with time, in the tropical western Pacific over the past 24 ka. The consistent trends in subsurface seawater pCO_2 and opal flux during deglaciation may imply that the reported increase in subsurface water pCO_2 in the study area was promoted by enhanced upwelling in the Southern Ocean.     

日照豪迈码头港池布局对泥沙输移影响研究

港池的布局会改变其周边海域水动力条件及其泥沙冲淤状况,为了进一步了解其影响程度,本文以山东日照豪迈重工临港厂区运输码头的改造工程为研究对象,利用平面二维数值模型MIKE21/3 Integrated Models,建立了潮流和波浪耦合作用下的泥沙输移数值模型,对该工程附近海域进行了波浪、潮流和泥沙输移的数值模拟。同时采用实测数据对数值模型进行了验证,数值模拟结果与实测资料拟合较好,表明MIKE21能有效地模拟运输码头及其周围海域潮流的变化过程。以日照豪迈重工临港厂区运输码头为依托,基于数值模拟结果,分析不同改造方案下的水动力条件和泥沙冲淤状况。结果表明:港池布局对水动力条件影响甚微,对地形地貌冲淤影响较大;模拟结果确定港池南向开口为最佳改造方案,为工程的规划和设计提供科学依据。     

GNSS-R/IR监测地表冻融状态对延迟多普勒波形和多路径数据影响分析

将GNSS-R/IR技术的应用领域拓展到地表冻融状态的监测中,本文利用冻融土混合介质介电常数模型计算土壤介电常数,采用双站全极化相干反射率模型和随机粗糙面散射模型,分别计算了经冻融土反射的GPS相干反射量的镜像反射率,以及GPS非相干反射分量的漫散射特性。模拟分析了冻融转换时,GPS多路径信息(GNSS-IR)以及包含漫散射信号的延迟多普勒图(GNSS-R)的变化特征。理论研究表明冻融转换过程中,地表介电常数的变化导致GPS多路径信息和延迟多普勒图的明显变化。本文从散射机理上揭示了利用GNSS-R和GNSS-IR遥感进行地表冻融特性监测的理论依据。     

面向双/多变量的连续面域拓扑图可视化方法

面域拓扑图是一种利用区域面积大小定量表达区域属性信息的可视化方法。由于其区域面积本身已经表示某一变量,因此这更有利于双/多变量的制图表达。针对目前基于面域拓扑图的双/多变量表达方法中存在的难以表达相邻区域之间基本状况和不利于不同地理现象的空间分布规律及差异表达的问题,本文提出一种面向双/多变量的连续面域拓扑图可视化方法。首先通过格网密度补偿和积分步长逐步试探的方法对基于扩散模型的连续面域拓扑图生成算法进行部分优化,完成基本变量的表达,然后分别通过空间内插和符号扩展完成第2和第3变量在连续面域拓扑图中的表达。最后以慕尼黑市人口密度和银行/ATM分布(双变量)数据以及奥格斯堡市人口密度数据、幼儿园分布以及规模数据(多变量)为试验数据进行可视化,并通过实证分析验证了该方法的有效性和优越性。     

机载WIDAS数据的Landsat卫星反照率初步验证

随着精细化监测的需求,中高空间分辨率的地表反照率产品逐渐成为气候模型的主要输入。目前,中高空间分辨率反照率产品的验证主要基于地表站点的通量塔观测数据,区域机载飞行数据的验证依然相对较少。因此,本文基于区域机载数据验证Landsat反照率产品。针对内蒙古自治区根河森林试验区所获取的机载红外广角双模式成像仪(WIDAS)多角度反射率数据,应用BRDF原型反演算法估算其反照率,分析了应用机载数据验证中高空间分辨率反照率产品的潜力。2016年内蒙古根河森林试验区机载WIDAS飞行多角度观测的可用多角度范围为25°,以前的研究表明BRDF原型反照率反演算法表现出对小观测角度的反照率反演结果的鲁棒性。因此,机载WIDAS反照率在一定程度可用于星载反照率的验证。首先,基于核驱动模型和各向异性平整指数(AFX)提取了试验区4种MODIS二向性反射分布函数(BRDF)原型;然后,将其作为先验知识应用到根河森林WIDAS机载数据的反照率反演中;最后,用WIDAS反照率和单个地面通量塔观测的反照率对Landsat卫星数据的反照率进行初步验证。验证结果表明Landsat反照率与WIDAS反照率结果较为一致,但略有低估,总体均方根误差(RMSE)约为0.02,偏差为0.0057。在多角度观测范围较小时,BRDF原型的反照率反演算法可为星载地表反照率的验证提供了一种有效的验证手段。     

DMSP/OLS夜间灯光数据的四川省GDP空间化分析

针对DMSP/OLS夜间灯光数据模拟省级格网GDP问题,该文选取2000、2006以及2010年DMSP/OLS稳定夜间灯光数据、辐射定标夜间灯光数据和四川省第二、三产业社会经济GDP数据,运用相关性分析、对比分析以及回归分析的方法,模拟出四川省GDP空间化密度图并分析与高程之间的关系。结果表明:在省级范围内,GDP与辐射定标灯光数据的相关性明显高于与稳定灯光数据的相关性,且前者的灯光指数CNLI与GDP相关性最高;运用3年的辐射定标灯光数据模拟得到四川省GDP密度图与GDP统计数据的拟合精度均达到0.99以上;2000-2010年四川省GDP呈现由成都市向四周不断扩张的趋势并随高程变化显著,四川省GDP高值区主要集中在高程1 000m以下的成都平原及周边区域,西部高原和山地地区GDP值较低。