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

新一代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(插补效果偏低、计算速度快)进行插补。     

Mid-late Holocene anthropogenic and natural variations in geochemistry and sedimentation in the Yazoo Basin: Clark Lake,Sharkey County,Mississippi

A 4000-yr history at Clark Lake reveals natural variations in sediment geochemistry and nutrient levels and anthropogenic influence on twentieth century sedimentation rates. Sediment texture and ridge and swale topography indicate that the channel system creating Clark Lake is now occupied by the Yazoo River. Between 1.24 and 0.60 m (2522–865 cal yr BP), total carbon percentages and the C/N ratio average 42% and 17, respectively. The base of the 1650-yr interval and onset of high organic activity in the lake coincides with the abandonment of the Yazoo Meander Belt. The top of the interval is marked by a drop in C/N ratio to 11 and a geochemical transition zone. No change in the rate of sediment accumulation or clay mineralogy was observed despite increased cumulative percentages of Al, Fe, K, Mg, Mn, Ca, and Na. The decrease in organic activity at 865 cal yr BP is attributed to unfavorable growth conditions related to the entrenchment of the Yazoo River and changes in the hydroperiod of the area. Significant floodplain alteration was completed in 1978, but the only major effect was increased sedimentation in the late twentieth century due to the completion of the Whittington Auxiliary Channel.     

1980—2015年风蚀影响下中国北方土壤有机质与养分流失时空特征

中国北方干旱/半干旱区是全球主要的沙尘源区之一,风蚀造成的沙尘排放可导致大量的土壤有机质（SOM）与养分流失,并通过传输与沉降过程对其进行空间再分配,对空气质量、气候变化、植被生长及生物地球化学过程等具有重要影响。本文利用WRF/Chem（Weather Research Forecasting with Chemistry）v3.7.1大气化学传输模型,对1980—2015年间中国北方沙尘排放及其引起的SOM、全氮（TN）与全磷（TP）的时空变化过程进行了精细化模拟,探究了中国北方风蚀引起的SOM、TN与TP养分流失的时空变化特征。结果表明：① 1980—2015年来平均每年约有66.59 Tg的沙尘颗粒排放至大气;② 沙尘排放具有较大的时空差异,沙尘排放源区主要集中在新疆东部、内蒙古西部的巴丹吉林沙漠和腾格里沙漠等地区;③ 每年因沙尘排放引起的SOM、TN和TP流失量分别约为0.07 Tg、0.004 Tg和0.005 Tg;④ 1980—2015年间中国北方沙尘排放及其引起的SOM、TN与TP的流失量具有较强的年际波动,未呈现显著的趋势性变化特征。本研究对于认知风蚀在碳循环与养分循环过程中的作用,以及对深入了解中国北方的土地退化机理具有重要意义。     

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.     

再探温跃层有孔虫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遥感进行地表冻融特性监测的理论依据。