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.     

Assessment and Prediction of Liquefaction Potential Using Different Artificial Neural Network Models: A Case Study

Soil liquefaction as a transformation of granular material from solid to liquid state is a type of ground failure commonly associated with moderate to large earthquakes and refers to the loss of strength in saturated, cohesionless soils due to the build-up of pore water pressures and reduction of the effective stress during dynamic loading. In this paper, assessment and prediction of liquefaction potential of soils subjected to earthquake using two different artificial neural network models based on mechanical and geotechnical related parameters (model A) and earthquake related parameters (model B) have been proposed. In model A the depth, unit weight, SPT-N value, shear wave velocity, soil type and fine contents and in model B the depth, stress reduction factor, cyclic stress ratio, cyclic resistance ratio, pore pressure, total and effective vertical stress were considered as network inputs. Among the numerous tested models, the 6-4-4-2-1 structure correspond to model A and 7-5-4-6-1 for model B due to minimum network root mean square errors were selected as optimized network architecture models in this study. The performance of the network models were controlled approved and evaluated using several statistical criteria, regression analysis as well as detailed comparison with known accepted procedures. The results represented that the model A satisfied almost all the employed criteria and showed better performance than model B. The sensitivity analysis in this study showed that depth, shear wave velocity and SPT-N value for model A and cyclic resistance ratio, cyclic stress ratio and effective vertical stress for model B are the three most effective parameters on liquefaction potential analysis. Moreover, the calculated absolute error for model A represented better performance than model B. The reasonable agreement of network output in comparison with the results from previously accepted methods indicate satisfactory network performance for prediction of liquefaction potential analysis.     

华北克拉通大青山-集宁地区古元古代变质沉积岩的锆石氧同位素组成:SHRIMP微区原位分析

在已获得锆石U-Pb年龄基础上,我们首次对孔兹岩带典型出露区大青山和集宁土贵乌拉地区古元古代变质沉积岩进行了锆石SHRIMP氧同位素研究。锆石具有复杂的内部结构和年龄分布。大青山地区古元古代早期变质沉积岩4个样品碎屑锆石的δ18O为5.52‰~7.11‰,部分重结晶锆石的δ18O为7.22‰~7.90‰,变质新生锆石的δ18O为6.37‰~8.31‰。大青山地区古元古代晚期变质沉积岩2个样品的锆石O同位素组成特征与古元古代早期的类似,另外2个样品的锆石O同位素组成与之不同,碎屑锆石、部分重结晶锆石和变质新生锆石的δ18O分别为6.26‰~10.80‰、9.00‰~11.20‰和9.66‰~11.90‰。集宁土贵乌拉地区古元古代晚期超高温变质沉积岩4个样品不存在碎屑锆石,变质锆石的δ18O变化范围为11.41‰~13.57‰。主要认识如下:1)大青山地区碎屑沉积物主要来自新太古代晚期-古元古代早期成熟度不高的TTG花岗质岩石物源区,与之相比,集宁土贵乌拉地区古元古代晚期超高温岩石的变质原岩成熟度更高;2)不同类型变质沉积岩变质新生锆石的δ18O和变质新生锆石与碎屑锆石的Δ18O存在明显区别,主要反映了岩石体系和变质流体的O同位素组成不同;3)大青山地区高角闪岩相-麻粒岩相变质沉积岩,重结晶锆石的O同位素完全重置,但U-Th-Pb体系未完全重置,集宁土贵乌拉地区超高温变质沉积岩的重结晶锆石O和U-Th-Pb同位素体系都完全重置。变质作用强度不同是主要原因。     

阿尔金南缘中段清水泉斜长花岗岩同位素年龄及成因研究

阿尔金南缘清水泉地区与基性-超基性岩伴生的花岗岩为斜长花岗岩。岩石地球化学显示该花岗岩高硅、富铝和钠,低镁和钾;轻稀土富集,具有Eu的正异常(δEu为1.01~2.01)。岩石富Rb、Ba,特别高Sr(779×10-6~864×10-6),低Y(1.17×10-6~1.51×10-6)及Yb(0.15×10-6~0.20×10-6),强烈亏损Nb、Ta等。斜长花岗岩锆石振荡环带清晰,Th/U和Nb/Ta比值分别为0.38~0.52,2.92~5.04;具有明显的Ce正异常和Eu负异常,为典型的岩浆锆石,利用LA-ICP-MS微区原位定年获得该花岗岩206Pb/238U-207Pb/235U谐和年龄为465Ma,206Pb/238U加权平均年龄为451±4Ma。锆石饱和温度计和锆石Ti温度计演算结果显示锆石的结晶温度分别为783~811℃和693~821℃。推测花岗岩源区压力范围为1.8~2.0GPa,形成深度在60km以上。综合分析清水泉花岗岩主、微量元素地球化学特征,并结合区域地质,认为该花岗岩属"I"型花岗岩,由地幔基性岩浆上侵分异形成,产于伸展环境。     

南水北调中线突发水污染事件的快速预测

依据南水北调中线干渠资料,开展了正常输水情况下串联明渠内可溶污染物浓度分布规律数值模拟研究。采用数值模拟、数学归纳和统计分析方法,提出表征污染物输移扩散特征的峰值输移距离、污染带长度和峰值浓度的快速预测公式;通过示范工程验证了快速预测公式的可行性。结果表明:①串联明渠内,峰值输移距离随渠道流速减小而减小,并且污染带长度增加值随明渠内流速减小而减小,但是峰值浓度随明渠流速减小而增加;②快速预测公式计算结果与现场试验实测结果的误差均不到15%,证明了快速预测公式的合理性和可行性。这些研究结果为南水北调中线工程突发可溶性水污染事件应急预警方案的制定提供了科学依据。     

湖南大义山岩体白沙子岭矿区钨锡矿深部找矿探索

大义山花岗岩体受"大义山式"北西向张扭性断裂构造控制,具多次侵位的特点,与岩体有关的矿产主要为锡多金属矿。本次补勘的白沙子岭锡矿区有云英岩脉型、蚀变岩体型、矽卡岩型等3类锡多金属矿类型。白沙子岭矿段平行排列的矿脉群在平面上整体呈一椭圆形展布,脉群从中间向两侧单脉的长度、厚度、品位逐渐变小,单个脉体从中间向两端厚度、品位逐渐变低。矿脉两端顶底板断层面上见到产状完全相反的擦痕,推断深部存在着一个椭圆形隐伏母岩体,成矿母岩的上部会形成蚀变岩体型矿体。为此本次补勘工作施工了一坑内深钻,上部见到了蚀变岩体型富锡矿体,下部见钨矿,预测深部找矿有较大的找矿前景。     

冰层取心回转钻进模拟实验台的设计与测试

极地钻探实践表明,冰盖底部冰岩交界附近地质情况异常复杂,不但可能存在暖冰、基底融水,甚至还存在厚度不等的冰岩夹层,取心钻探异常困难,而优选钻头类型、确定合理的钻进参数是保证其安全、快速钻进的重要因素。本文设计了一套能够模拟冰层回转钻进的实验台,其技术参数为：钻压、转速分别可在0～10 kN、0～300 r/min范围内调节,最大扭矩约100 N·m。该实验台能够测量钻头切削具温度、钻孔深度及钻进速度等参数,为深入研究钻头类型及结构参数、钻进参数对扭矩、钻速和切削温度的影响规律提供了手段。采用PDC复合片钻头进行了冰钻实验,结果表明,实验台能够准确调节钻压和转速,可满足实验要求。     

Reconstruction of Holocene paleohydrological and sedimentological characteristics of Umm Ar-Rimam basin (northern Kuwait)

The paleohydrological and sedimentological characteristics of a playa lake in northern Kuwait (Arabian Gulf) are reconstructed using sedimentological, geochemical, and isotopic techniques. The sequence consists of up to 8 cycles of S-poor, alluvial sediments capped by a thin organic soil interbedded with gravity-fall calcrete sediments. The succession is locally derived from mainly Quaternary sediments and is regressive with upsection filling of the subsiding basin by cycles of sheetwash flow in response to climatic change. There is no natural, open-water lake water as indicated by low total organic carbon (TOC) data, but the presence of incised calcrete yardangs suggests that more extensive open-water conditions were operative in the past. Stable isotope (δ¹⁸O‰ and δ¹³C‰) values of the authigenic carbonates indicate the following three distinct processes: evaporation, meteoric fluid infiltration, and rapid per-descensum flow (rapid downward movement of water and playa sediment through pipes) through a porous, clastic sequence. Because evaporites are scarce, other factors besides evaporation action control chemical and isotopic compositions of the per-decensum lake fluids. Consequently, the isotopic composition cannot be interpreted exclusively as an indicator of salinity or evaporation ratio. The degassing of CO₂ during groundwater discharge may explain the enriched carbon isotope values for the authigenic carbonates precipitated in the sediments. Hydrologically closed lake water bodies tend to show low negative carbonate oxygen and carbon isotopic signatures. Isotopically negative δ¹³C values imply a strong input of soil-zone carbon to the groundwater of the top 60 cm of the sediment. Lakes that are hydrologically closed and evaporate or equilibrate with atmospheric CO₂ will tend to have low negative δ¹⁸O and δ¹³C values in the carbonates as reported by Talbot (Chem Geol: Isotope Geosci Sect 80(4):261–279, 1990). Biologically active lakes will tend toward lower δ¹³C of dissolved carbon due to the photosynthetic effects of ¹²C withdrawal as reported by Dunagan and Driese (J Sed Res 69:772–783, 1999). Increased biological activity during sedimentation may account for low carbon isotope values where plants were abundant, but in shrinkage-dominated systems (those of clay-rich soil subjecting to wet-dry conditions), carbon isotopes will be largely inherited from the calcretic limestones in the land extending landward of the coast and not influenced by coastal processes (known as hinterland), such as Umm Ar-Rimam depression. This basin does not fit the classic shallow playa-type basins of the Arabian Peninsula but rather the recharge playas of the southwestern USA.