青藏高原拉萨地块西部念青唐古拉岩群的地球化学特征及构造意义

拉萨地块西部呈断块状沿狮泉河-申扎-嘉黎蛇绿混杂岩带附近分布的念青唐古拉岩群被认为是前寒武纪变质基底。本文对念青唐古拉岩群进行了系统的岩石学、地球化学、同位素年代学及构造地质学研究。研究结果表明片岩-片麻岩-变粒岩含十字石、石榴子石等特征变质矿物,遵循粒度分异规律,其原岩可能为来自冈瓦纳古陆核北缘中新元古代弧盆体系的活动大陆边缘浊积岩。斜长角闪岩具低硅、高铁镁、富钙的基性岩特征,其原岩为岛弧型基性火山岩。念青唐古拉岩群中的花岗伟晶岩锆石LA-ICPMS U-Pb年龄为1150±13Ma,具过铝质S型花岗岩地球化学特征,可能为中元古代(1150±13 Ma)以前就开始沉积的念青唐古拉岩群基底岩石通过部分熔融形成。与花岗伟晶岩渐变过渡接触的二云斜长片麻岩第一组变质重结晶锆石U-Pb年龄为701±15 Ma,结合十字石特征变质矿物,暗示了该地区中温高压变质作用的峰期变质,变质程度达角闪岩相;第二组热液流体锆石UPb年龄为301±8.4 Ma,可能与冈瓦纳大陆北缘古特提斯洋演化过程中的岩浆热液作用有关。     

Initiation of large-volume silicic centers in the Yellowstone hotspot track: insights from H<Subscript>2</Subscript>O- and F-rich quartz-hosted rhyolitic melt inclusions in the Arbon Valley Tuff of the Snake River Plain

During the onset of caldera cluster volcanism at a new location in the Snake River Plain (SRP), there is an increase in basalt fluxing into the crust and diverse silicic volcanic products are generated. The SRP contains abundant and compositionally diverse hot, dry, and often low-δ¹⁸O silicic volcanic rocks produced through time during the formation of individual caldera clusters, but more H₂O-rich eruptive products are rare. We report analyses of quartz-hosted melt inclusions from pumice clasts from the upper and lower Arbon Valley Tuff (AVT) to gain insight into the initiation of caldera cluster volcanism. The AVT, a voluminous, caldera-forming rhyolite, represents the commencement of volcanism (10.44 Ma) at the Picabo volcanic field of the Yellowstone hotspot track. This is a normal δ¹⁸O rhyolite consisting of early and late erupted members (lower and upper AVT, respectively) with extremely radiogenic Sr isotopes and unradiogenic Nd isotopes, requiring that ~50 % of the mass of these elements is derived from melts of Archean upper crust. Our data reveal distinctive features of the early erupted lower AVT melt including: variable F concentrations up to 1.4 wt%, homogenous and low Cl concentrations (~0.08 wt%), H₂O contents ranging from 2.3 to 6.4 wt%, CO₂ contents ranging from 79 to 410 ppm, and enrichment of incompatible elements compared to the late erupted AVT, subsequent Picabo rhyolites, SRP rhyolites, and melt inclusions from other metaluminous rhyolites (e.g., Bishop Tuff, Mesa Falls Tuff). We couple melt inclusion data with Ti measurements and cathodoluminescence (CL) imaging of the host quartz phenocrysts to elucidate the petrogenetic evolution of the AVT rhyolitic magma. We observe complex and multistage CL zoning patterns, the most critical being multiple truncations indicative of several dissolution–reprecipitation episodes with bright CL cores (higher Ti) and occasional bright CL rims (higher Ti). We interpret the high H₂O, F, F/Cl, and incompatible trace element concentrations in the context of a model involving melting of Archean crust and mixing of the crustal melt with basaltic differentiates, followed by multiple stages of fractional crystallization, remelting, and melt extraction. This multistage process, which we refer to as distillation, is further supported by the complex CL zoning patterns in quartz. We interpret new Δ¹⁸O(Qz-Mt) isotope measurements, demonstrating a 0.4 ‰ or ~180 °C temperature difference, and strong Sr isotopic and chemical differences between the upper and lower AVT to represent two separate eruptions. Similarities between the AVT and the first caldera-forming eruptions of other caldera clusters in the SRP (Yellowstone, Heise and Bruneau Jarbidge) suggest that the more evolved, lower-temperature, more H₂O-rich rhyolites of the SRP are important in the initiation of a caldera cluster during the onset of plume impingement.     

面向GIS增量数据的订阅与发布机制研究

在分析已有数据分发策略的基础上,提出了一种基于"订阅与发布"机制的GIS增量数据快速分发方法。该方法首先利用用户的注册信息记录用户感兴趣的数据范围或类型,当服务器接收到增量数据时,根据不同的用户对增量数据进行筛选,并提出了一种基于联想模式的数据关联查询策略,再将筛选出来的数据发送给相应的客户端用户。经试验论证,该方法能够针对不同的用户需求,将GIS增量数据快速准确地发送给客户端用户,有助于数据的有效利用和管理。     

奇异值分解法在病态问题中的应用

用截断奇异值分解法和修正奇异值分解法对大地测量病态问题进行了处理,并与最小二乘的结果进行了比较,最后将两者方法进行结合同样对病态方程进行了处理,得到了结合奇异值分解的解,并与截断奇异值、修正奇异值分解的解以及真值进行比较,发现结合奇异值分解的解即修正奇异值截断法比截断奇异值和修正奇异值的解更加靠近于真值,修正奇异值截断法相比于截断奇异值和修正奇异值法对于病态方程抗干扰能力更强,更具有实际意义。     

北斗三频伪距相关随机模型单站建模方法

为了提供可靠的伪距随机模型,基于3个线性无关的北斗三频伪距/载波无几何无电离层(GIF)组合,研究了一种利用单站数据估计北斗三频伪距相关随机模型算法。该算法首先利用低阶多项式对GIF组合进行拟合,以尽可能消除非伪距噪声以外的其他常数和误差,然后利用多元线性回归分析实现对3个线性无关的GIF组合随机噪声同时建模,最后再由线性组合关系变换得到原始北斗三频伪距相关随机模型。经北斗三频实测数据验证结果表明,该算法可实现北斗非差三频伪距相关随机模型的单站解算,有利于为导航定位以及完好性监测提供精准随机模型。     

一种快速二维 Delaunay 三角网点定位算法

在构建二维Delaunay三角网的逐点插入法中,定位待插点所在三角形的快慢是影响整个算法构网速度的关键因素。针对目前已有算法存在的搜索路径长、搜索路径求解计算量大等问题,结合三角形重心的几何性质,对点定位算法进行改进,避免求三角形重心和相交边的过程。实验结果表明,文中算法较目前其他点定位算法能够有效地缩短搜索路径,减少点定位的计算时间,提高Delaunay三角网构网过程中点定位的效率。     

面向应用的航天遥感科学论证研究

为适应中国民用航天遥感从科学试验向业务服务模式转变,更好地探索、了解与解决应用需求与航天遥感系统对接等方面遇到的技术问题,促进航天遥感统筹协调可持续发展,中国适时于2004年成立了国家航天局航天遥感论证中心。10余年来,论证中心以航天遥感系统为研究对象,系统开展了面向应用的航天遥感科学论证概念、理论方法、技术工程与应用研究。本文是论证中心团队长期从事航天遥感科学论证研究与实践的系统总结,介绍了遥感论证初步认知、遥感论证关注问题、遥感论证理论体系与模型方法集、遥感论证能力建设及遥感论证实践等方面内容,给出了遥感论证定义并详细分析了研究范围和内容,提出了由知识维、进程维和逻辑维所组成的遥感论证作用域3维空间结构,指出社会发展加快和信息化水平提高,带动整个航天遥感数据信息链向更大规模、更短响应时间周期、更综合数据集成、更高数据质量、更加智能化方向发展,航天遥感系统将进入新的"智慧遥感"发展阶段。得益于十余年来中国民用航天快速发展,我们经历了风云三号新型载荷校飞、多角度多光谱偏振遥感器论证、环境星应用工程论证等实践,取得了多方面理论方法的突破,并应用到2030民用航天发展规划、高分辨率对地观测系统、国家自然灾害空间信息基础设施、国家民用空间基础设施中长期发展、2030中国综合地球观测系统规划等论证当中。经过不断实践,快速迭代,形成了遥感论证理论体系及相应的十大模型方法,包括遥感信息流模型、遥感信息特征模型、遥感信息应用模型、遥感信息量分析模型、遥感数据工程模型、航天遥感系统结构模型、航天遥感系统状态描述模型、航天遥感系统质量模型、航天遥感系统发展动力模型及能力体系模型。这些模型方法全面反映了航天遥感系统特征、结构、状态、发展动力、条件等,可广泛用于对航天遥感系统进行顶层设计、规划、考察、分析、评价、预测,并开展实践探索。     

Recursive algorithm for fast GNSS orbit fitting

Gaussian elimination is an efficient and numerically stable algorithm for estimating parameters and their precision. However, before estimating the parameters, it is often prudent to perform statistical tests to achieve the best fitting model. We use Gaussian elimination to select the best fitting model among candidate models. A succinct relationship between the weighted sum of squared residuals and the previous one is revealed by a volume formula. For quick parameter estimation and determination of weighted sum of squared residuals, a recursive elimination algorithm is proposed in the context of Gaussian elimination. In order to improve the model selection efficiency, the parameter estimation and the determination of the weighted sum of squared residuals are carried out in parallel using the proposed recursive elimination algorithm in which the improvement at each recursive stage is judged by the Bayesian information criterion. Ultimately, the computational complexity and numerical stability of the recursive elimination proposed are briefly discussed, and a GNSS orbit interpolation example is used to verify the results. It shows that the proposed recursive elimination algorithm inherits the numerical stability of the Gaussian elimination, and this algorithm can be used to examine the gain from the newly introduced parameter, dynamically assess the fitting model, and fix the optimal model efficiently. The optimal fitting model with the lowest information is very close to the real situation verified by checkpoints.     

A new method for specular and diffuse pseudorange multipath error extraction using wavelet analysis

Multipath remains one of the major challenges in Global Navigation Satellite System (GNSS) positioning because it is considered the dominant source of ranging errors, which can be classified into specular and diffuse types. We present a new method using wavelets to extract the pseudorange multipath in the time domain and breaking it down into the two components. The main idea is an analysis-reconstruction approach based on application of both continuous wavelet transform (CWT) and discrete wavelet transform (DWT). The proposed procedure involves the use of L1 code-minus-carrier (CMC) observable where higher-frequency terms are isolated as residuals. CMC residuals are analyzed by applying the CWT, and we propose the scalogram as a technique for discerning time–frequency variations of the multipath signal. Unlike Fourier transform, the potential of the CWT scalogram for examining the non-stationary and multifrequency nature of the multipath is confirmed as it simultaneously allows fine detection and time localization of the most representative frequencies of the signal. This interpretation of the CWT scalogram is relevant when choosing the levels of reconstruction with DWT, allowing accurate time domain extraction of both the specular and diffuse multipath. The performance and robustness of the method and its boundary applicability are assessed. The experiment was carried out using a receiver of Campania GNSS Network. The results are given in which specular multipath error is achieved using DWT level 7 approximation component and diffuse multipath error is achieved using DWT level 6 denoised detail component.     

Optimization and Parametrical Investigation to Assess the Reconstitution of Different Types of Indian Sand Using Portable Travelling Pluviator

By using the air pluviation technique, it is aimed to achieve the desired relative density with uniform void ratio throughout the specimen in order to maintain the homogeneity and to avoid the spatial variability. Further, in order to achieve the maximum deposition intensity, a systematic optimization study has been carried out rigorously in a test tank to determine the diameter of the orifice to be employed for the sieve plates of different porosity and the number of sieve plates to be installed in the diffuser sieve sets. The study has been conducted with four different patterns of sieves with different porosity to achieve a wide range of relative densities for four different uniformly graded Indian sands. The dynamic penetrometer which is considered to be one of the cost effective instruments has been efficiently used to determine the soil resistance at various locations of the test tank for every given height of fall in order to check the uniformity of placement density throughout the sand bed. The study reveals that the sand beds of different relative densities could be achieved using different patterns of diffuser sieves at optimum sand flow rate without compromising the uniformity. The effect of height of fall as well as porosity of diffuser sieves on the relative density of different sands has been studied in detail. The deposition intensity and the relative density obtained from the present study are compared with the values available in the literature.