Improved spatial delineation of streambed properties and water fluxes using distributed temperature sensing

A new method was developed for analysing and delineating streambed water fluxes, flow conditions and hydraulic properties using coiled fibre‐optic distributed temperature sensing or closely spaced discrete temperature sensors. This method allows for a thorough treatment of the spatial information embedded in temperature data by creating a matrix visualization of all possible sensor pairs. Application of the method to a 5‐day field dataset reveals the complexity of shallow streambed thermal regimes. To understand how velocity estimates are affected by violations of assumptions of one‐dimensional, saturated, homogeneous flow and to aid in the interpretation of field observations, the method was also applied to temperature data generated by numerical models of common field conditions: horizontal layering, presence of lateral flow and variable streambed saturation. The results show that each condition creates a distinct signature visible in the triangular matrices. The matrices are used to perform a comparison of the behaviour of one‐dimensional analytical heat‐tracing models. The results show that the amplitude ratio‐based method of velocity calculation leads to the most reliable estimates. The minimum sensor spacing required to obtain reliable velocity estimates with discrete sensors is also investigated using field data. The developed method will aid future heat‐tracing studies by providing a technique for visualizing and comparing results from fibre‐optic distributed temperature sensing installations and testing the robustness of analytical heat‐tracing models. Copyright © 2016 John Wiley & Sons, Ltd.     

基于NDVI时间序列数据的江西省水稻种植制度变化研究

利用Savitzky-Golay滤波对覆盖江西省范围的SPOT VGT NDVI时间序列数据进行平滑处理的基础上,结合坡度数据,通过非监督分类的方法提取了江西省2000、2005和2010年水稻种植范围,并根据NDVI的年内动态变化,从水稻种植范围、水稻生长季起始时间、水稻复种指数和NDVI最大振幅等分析了江西省水稻种植和生长情况,探讨2000~2010年江西省水稻生产的变化。     

L1范数与IQR统计量组合的GNSS坐标序列粗差探测算法

GNSS坐标时间序列中不可避免地含有粗差,未剔除的粗差将会导致参数估计有偏。因此,粗差探测与剔除是GNSS坐标序列分析中一项重要的数据预处理工作。针对GNSS坐标时间序列特点,提出了一种将L1范数(L1-norm)估计与四分位距统计量IQR(interquartile range)组合的移动开窗粗差探测算法,称之为L1_Mod IQR。该方法的主要思想是,首先利用L1范数估计得到较"真实"的残差,然后再对残差采用IQR统计量进行粗差探测。将L1_Mod IQR法与"3σ"法、基于最小二乘的τ检验法等粗差探测算法进行了模拟计算与对比,验证了该算法的有效性。进一步采用L1_Mod IQR算法对中国区域10个IGS站的高程时间序列进行了分析,结果表明中国区域IGS站高程序列的粗差剔除率最小为0.1%,最大为2.6%。并且以WUHN站为例与SOPAC提供的结果进行了对比,结果表明SOPAC提供的"Clean"数据仍含有大量的粗差,而L1_Mod IQR算法能够有效地剔除粗差。     

Significance of shallow core transects for reef models and sea‐level curves,Heron Reef,Great Barrier Reef

A sequence of shallow reef cores from Heron Reef, Great Barrier Reef, provides new insights into Holocene reef growth models. Isochron analysis of a leeward core transect suggests that the north‐western end of Heron Reef reached current sea‐level by ca 6·5 kyr bp and then prograded leeward at a rate of ca 19·6 m/kyr between 5·1 kyr and 4·1 kyr bp (pre‐1950) to the present reef margin. A single short core on the opposing margin of the reef is consistent with greater and more recent progradation there. Further to the east, one windward core reached modern sea‐level by ca 6·3 kyr bp , suggesting near ‘keep‐up’ behaviour at that location, but the opposing leeward margin behind the lagoon reached sea‐level much more recently. Hence, Heron Reef exhibited significantly different reef growth behaviour on different parts of the same margin. Mean reef accretion rates calculated from within 20 m of one another in the leeward core transect varied between ca 2·9 m and 4·7 m/kyr depending on relative position in the prograding wedge. These cores serve as a warning regarding the use of isolated cores to inform reef growth rates because apparent aggradation at any given location on a reef varies depending on its location relative to a prograding margin. Only transects of closely spaced cores can document reef behaviour adequately so as to inform reef growth models and sea‐level curves. The cores also emphasize potential problems in U‐series dates for corals within a shallow (ca 1·5 m) zone beneath the reef flat. Apparent age inversions restricted to that active diagenetic zone may reflect remobilization and concentration of Th in irregularly distributed microbialites or biofilms that were missed during sample vetting. Importantly, the Th‐containing contaminant causes ages to appear too old, rather than too young, as would be expected from younger cement.     

年最大值抽样估算的暴雨频率设计值的低估问题探讨

由年最大值抽样（AMS）和年超大值抽样（AES）的基本理论及重现期（RP）的定义可知,AMS并不符合以“事件”为基础的重现期的定义.以美国西南半干旱区1438个雨量站和太湖流域96个雨量站的降雨资料为例,通过经验频率与超过概率的比较,发现AMS估算的暴雨频率设计值偏小,尤其是对常遇频率降雨设计值的影响更加显著.美国的降雨量资料站点多、系列长,实际资料验证与理论分析一致.通过对太湖流域AMS资料的分布形态进行分析的结果表明：太湖流域的站点不多,资料长度不够,且大部分站点在雨量大值区数据稀少,使得频率直方图不连续,是造成我国太湖流域的资料验证效果不理想的可能原因.     

慢特征分析法在气象上的应用进展

气候系统具有非平稳特征,根本原因在于其外强迫随时间发生改变,因此外部驱动力的分析对于理解气候系统的动力学特征至关重要,而如何有效提取系统外部驱动信息是一个亟待解决的前沿科学问题。最近几年,在生物神经学领域中应用的一种提取非平稳信号中外强迫信息的方法——慢特征分析法(Slow Feature Analysis,SFA),在气象领域中也得到了初步成功的尝试,结果显示出此方法对气候系统的外强迫信息分析及有关动力学机制的探究有较好的应用前景。本文主要介绍SFA方法的理论思想及实施步骤,并通过一个理想的非平稳时间序列检验其提取外强迫信息的能力,结果证明在衰减的Logistic模型中,可利用SFA算法提取出模型中的外强迫,且与真实外强迫的相关系数可达0.99;此外,还介绍将该方法应用于Arosa臭氧时间序列,分析其提取的外强迫信息的动力学特征;并介绍了在气候时间序列建模中引入外强迫因子的预测效果。     

中国东部夏季暴雨的年代际跃变及其大尺度环流背景

本文利用1960~2011年中国东部地面测站的逐日降水资料和JRA-55再分析资料探讨了夏季暴雨分布的年代际跃变及其相关联的大尺度环流异常特征。基于暴雨频数和占比(夏季暴雨占比是指5~8月暴雨降水量对总降水量的贡献百分比)的分析结果表明:中国东部夏季暴雨分布在20世纪70年代末和90年代初经历两次反相的经向"三极子"跃变。中国东部夏季暴雨的年代际演变过程可分为三个时段:1960~1979年为华南和华北暴雨偏多、江淮流域暴雨偏少的经向"三极子"分布;1980~1991年为南方和华北暴雨偏少、江淮流域暴雨偏多的经向三极子"分布;1992~2011年为南方暴雨显著偏多、华北暴雨持续偏少,逐渐形成经向"偶极子"分布,并导致近十多年我国夏季"南涝北旱"的整体格局。1970年代末(1990年代初)跃变相关联的大尺度环流异常配置:东亚夏季风的减弱(增强),西太平洋副高的增强西伸但南撤(北抬),南亚高压的减弱南缩(增强东扩),以及蒙古高原中低层的气旋式(反气旋式)环流异常。与此同时,低层局地环流也发生调整:华北和黄淮地区以及华南和江南地区均为反气旋式(气旋式)环流异常,而江淮流域和四川盆地受控于风场切变式辐合(辐散)异常;涡度场发生相应变化,南北方大部分地区的负(正)涡度异常不(有)利于低涡的发展,而江淮流域和四川盆地的正(负)涡度异常有(不)利于低涡的发展,进而引发江南和华南暴雨减少(增加)、江淮流域和四川盆地暴雨增加(减少)、黄淮和华北暴雨减少(增加)的经向"三极子"跃变。     

Study of high temperature targets identification and temperature retrieval experimental model in SWIR remote sensing based Landsat8

For surface features in short-wave infrared (SWIR, 1.3–3.0 μm) in remote sensing imagery, pixel values depict the total energy including reflection and emission. For surface features at normal temperature in SWIR band, emission energy can be ignored. While for surface features at high temperature in SWIR band, emission energy is equal to or even higher than the reflection energy. So remote sensing imagery of SWIR band can be used to separate emission and reflection energy as well as to realize temperature retrieval of high temperature targets. In this study, the seventh band (SWIR band) of Landsat8 OLI remote sensing imagery is used to perform the theoretical model research for temperature retrieval of high temperature targets. In the meantime, it is also used with the corresponding observation experiment of synchronization satellite to check the theoretical model. The result shows that the radiant flux density for mixed pixels with high temperature targets is higher than adjacent pixels without high temperature targets. Thus, the high temperature pixels can be identified in SWIR band. The retrieval results of temperature and fractional area for high temperature targets are consistent with reality. In the study, the result illustrates that it is effective to identify high temperature targets in remote sensing imagery of SWIR band and the model is appropriate for temperature retrieval use.     

Regional assessment of pan-Pacific urban environments over 25 years using annual gap free Landsat data

Urbanization and the associated change in land cover has been intensifying across the globe in recent decades. Regional studies on the rate and amount of urban expansion are critical for understanding how patterns of change differ within and among cities with varying structure and development characteristics. Yet spatially consistent and timely information on urban development is difficult to access particularly across international jurisdictions. Remote sensing based technologies offer a unique perspective on urban land cover with the data offering significant potential to urban studies due to its consistent and ubiquitous nature. In this research we applied a pixel-based image composite technique to generate annual gap-free surface reflectance Landsat composites from 1984 to 2012 for 25 urban environments across 12 countries in the Pacific Rim. Using time series composites, spectral indices were calculated and compared using a hexagonal grid ring model to assess changes in vegetative and urban patterns. Trajectories were then clustered to further investigate the spatio-temporal dynamics and relationships among the 25 cities. Performance of the clustering analyses varied depended on the temporal and spatial metrics however overall clustering results indicated relatively strong spatio-temporal similarities among a number of key cities. Three pairs of cities—Melbourne and Sydney; Tianjin and Manila; and Singapore City and Kuala Lumpur were found to be highly similar in their urban and vegetation dynamics temporally and spatially. In contrast Vancouver and Las Vegas had no similar analogous. This work demonstrates the value of utilising annual Landsat time series composites for assessing urban vegetation and urban dynamics at regional scales and potential use in sustainable urban planning, resources allocation, and policy making.     

基于 DTW 距离的时序 NDVI 数据植被信息提取--以秦巴山区为例

秦巴山区是我国重要的生态屏障,对该区的植被信息提取开展研究,可为区内生态服务功能及自然资源开发利用提供基础数据。通过加窗处理改进DTW距离相似性算法,结合临近度模糊分类方法对2005—2014年的MODIS NDVI时序数据进行植被信息提取。首先利用S-G滤波对MODIS NDVI时序数据进行重建;再利用2013年的采样数据构建各类植被的标准NDVI时序曲线,逐像元计算与标准NDVI时序曲线的加窗DTW距离,利用临近度模糊分类实现植被信息提取;最后验证提取精度。结果表明,算法具有较高的运行效率,可避免错误匹配,以较高的精度(总体精度83.8%,kappa系数0.77)实现长时间序列的植被信息提取。