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算法能够有效地剔除粗差。     

Provenance of Holocene beach sand in the Western Iberian margin: the use of the Kolmogorov–Smirnov test for the deciphering of sediment recycling in a modern coastal system

Detrital zircons from Holocene beach sand and igneous zircons from the Cretaceous syenite forming Cape Sines (Western Iberian margin) were dated using laser ablation – inductively coupled plasma – mass spectrometry. The U–Pb ages obtained were used for comparison with previous radiometric data from Carboniferous greywacke, Pliocene–Pleistocene sand and Cretaceous syenite forming the sea cliff at Cape Sines and the contiguous coast. New U–Pb dating of igneous morphologically simple and complex zircons from the syenite of the Sines pluton suggests that the history of zircon crystallization was more extensive (ca 87 to 74 Ma), in contrast to the findings of previous geochronology studies (ca 76 to 74 Ma). The U–Pb ages obtained in Holocene sand revealed a wide interval, ranging from the Cretaceous to the Archean, with predominance of Cretaceous (37%), Palaeozoic (35%) and Neoproterozoic (19%) detrital‐zircon ages. The paucity of round to sub‐rounded grains seems to indicate a short transportation history for most of the Cretaceous zircons (ca 95 to 73 Ma) which are more abundant in the beach sand that was sampled south of Cape Sines. Comparative analysis using the Kolmogorov–Smirnov statistical method, analysing sub‐populations separately, suggests that the zircon populations of the Carboniferous and Cretaceous rocks forming the sea cliff were reproduced faithfully in Quaternary sand, indicating sediment recycling. The similarity of the pre‐Cretaceous ages (>ca 280 Ma) of detrital zircons found in Holocene sand, as compared with Carboniferous greywacke and Pliocene–Pleistocene sand, provides support for the hypothesis that detritus was reworked into the beach from older sedimentary rocks exposed along the sea cliff. The largest percentage of Cretaceous zircons (<ca 95 Ma) found in Holocene sand, as compared with Pliocene–Pleistocene sand (secondary recycled source), suggests that the Sines pluton was the one of the primary sources that became progressively more exposed to erosion during Quaternary uplift. This work highlights the application of the Kolmogorov–Smirnov method in comparison of zircon age populations used to identify provenance and sediment recycling in modern and ancient detrital sedimentary sequences.     

天山南麓雷暴日短期预报研究

利用阿克苏及邻近地区12个气象站1980—2013年雷暴资料,以及同期高空资料,统计了各站年均雷暴日数,对发生区域雷暴天气的环流形势进行分类,归纳出各型的入型指标。通过逐步回归法,建立阿克苏及邻近地区区域雷暴概率回归预报模型,并对2013年进行试预报。结果表明:(1)阿克苏及邻近地区区域雷暴的影响系统主要分为4类:巴湖低槽型、急流型、西北气流型和温度槽型。(2)对2002—2012年5—9月(共1683 d)历史资料进行判别,满足入型条件的样本数为876 d,消空率为48%;对2013年5—9月(共153 d)历史资料进行判别,入型样本数为80 d,消空率为48%。(3)对2002—2012年5—9月所有入型样本进行回代检验,平均准确率为72.0%(平均TS评分为30.1%);对2013年5—9月所有入型样本进行试预报,平均准确率为63.2%(平均TS评分为28.2%)。     

风云卫星数据服务评价指标的验证分析

简述了目前我国风云卫星数据服务以及气象服务效益评估业务的现状,利用层次分析法提出风云卫星数据服务评价的6个1级指标和19个2级指标。设计了风云卫星数据服务评测量表,并采用分半信度、因子分析等方法对量表作信度和效度分析,结果表明所设计的评测量表具有良好的信度和效度,确定了包括服务平台易用性、服务流程规范性、服务响应及时性、服务内容完整性、服务工作主动性、服务宣传广泛性等6个方面的风云卫星数据服务评价指标。     

部分变量误差模型的整体抗差最小二乘估计

部分变量误差模型(partial EIV model)的加权整体最小二乘(weighted total least-squares,WTLS)估计不具备抵御粗差的能力。鉴于粗差可能同时出现在观测值和系数矩阵中,本文在提出部分变量误差模型WTLS估计的两步迭代解法的基础上,运用抗差M估计的等价权方法,发展了一种整体抗差最小二乘(TRLS)估计方法,并采用一致最大功效统计量确定降权因子。针对WTLS估计两步迭代解法的特点,设计了两个不同的降权方案:第1个方案是在估计系数矩阵元素时,不对观测值降权,仅对系数矩阵降权;第2个方案是在估计系数矩阵元素时,既对系数矩阵降权,同时也对观测值降权。通过对模拟2D仿射变换和线性拟合实例进行计算和分析,结果表明第1方案优于第2方案,并且优于基于残差和验后单位权方差的抗差估计和现有的变量误差模型抗差估计。     

空间层次聚类显著性判别的重排检验方法

同时顾及空间邻近与专题属性相似的空间层次聚类是挖掘空间分布模式的一种有效手段。空间层次聚类方法虽然可以获得多层次的聚集结构,但聚类结果显著性的统计判别依然是一个尚未解决的难题。为此,本文提出了一种空间层次聚类结果显著性的统计判别方法,用于确定空间层次聚类的停止准则,减少聚类过程对参数设置的依赖。通过试验分析与比较发现,该方法能够有效判别空间层次聚类结果的显著性和确定层次聚类合并过程的停止条件,同时具有很好的抗噪性,避免随机结构的干扰。     

非量测相机影像三维模型构建及精度检验

针对三维模型构建投资大、周期长、技术难度较高的问题,该文提出了基于Agisoft PhotoScan的非量测相机影像三维模型构建及精度检验方法。在对非量测数码相机数据进行预处理的基础上,利用Agisoft PhotoScan软件进行数据定向、点云提取、建立不规则三角网及数字高程模型和数字正射影像的制作;根据生成的数字高程模型和正射影像建立研究区域的三维场景模型。实验结果证明,该方法可快速建立高精度三维模型,缩短建模时间并降低投资成本。     

Contamination,bioaccessibility and human health risk of heavy metals in exposed-lawn soils from 28 urban parks in southern China's largest city,Guangzhou

The total concentrations and oral bioaccessibility of heavy metals in surface-exposed lawn soils from 28 urban parks in Guangzhou were investigated, and the health risks posed to humans were evaluated. The descending order of total heavy metal concentrations was Fe > Mn > Pb > Zn > Cu > Cr > Ni > Cd, but Cd showed the highest percentage bioaccessibility (75.96%). Principal component analysis showed that Grouped Cd, Pb, Cr, Ni, Cu and Zn, and grouped Cr and Mn could be controlled two different types of human sources. Whereas, Ni and Fe were controlled by both anthropogenic and natural sources. The carcinogenic risk probabilities for Pb and Cr to children and adults were under the acceptable level (<1 × 10⁻⁴). Hazard Quotient value for each metal and Hazard Index values for all metals studied indicated no significant risk of non-carcinogenic effects to children and adults in Guangzhou urban park soils.     

阿根廷米纳毕戈塔地区水系沉积物采样粒度试验研究

中国地质调查局与阿根廷地质矿产调查局合作在阿根廷西北部米纳毕戈塔地区进行1︰25万水系沉积物测量,该区属于干旱-半干旱高寒山区,为突出找矿效果,消除或减少风积物干扰,需要确定适合于该区的水系沉积物采样粒度。为此,在该区选择一个有已知矿床的1︰5万图幅进行采样粒度试验,分别选择10~60目、-60目、60~80目、-80目4个粒度级进行粒度试验;根据我国区域地球化学勘查规范,每个粒度分析39种元素。结果表明,大部分元素在4种粒度水系沉积物中的分布形态基本一致,都能够较好地反映出区内已知矿床,说明采用的采样方法可以有效地避免风成砂干扰;金、银、铜、锡等成矿元素在4种粒度中的分布略有差异,10~60目的金、铜和锡异常对已知矿体的反应更清晰准确。因此,本区水系沉积物地球化学测量的最佳采样粒度是10~60目。根据本次试验结果,中-阿地调局在本区联合开展1︰25万水系沉积物地球化学测量时采用了10~60目水系沉积物作为采样介质,取得了很好的效果。这是本区第一次进行水系沉积粒度试验,对本区将来的地球化学调查和研究具有指导意义。