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.     

冀北中低级变质表壳岩的年代学特征及其对华北克拉通构造演化的约束

华北克拉通在新太古代末期发生克拉通化,形成了现今规模的古陆,大量的太古宙岩石均经历了~2500Ma左右的区域高级变质作用(高角闪岩相-麻粒岩相)。而华北克拉通北部冀北地区出露一套中低级变质(绿片岩相-角闪岩相)的火山-沉积岩系,主要包括胡麻营地区红旗营子表壳岩和大阴山地区单塔子表壳岩中变质程度较低的部分。胡麻营地区红旗营子表壳岩系主要岩石组合为变基性火山岩、绿帘角闪岩、斜长角闪岩、含石榴石斜长角闪岩、角闪斜长片麻岩、黑云斜长片麻岩、黑云角闪斜长片麻岩、黑云二长片麻岩、石英片岩、磁铁石英岩等,SIMS锆石U-Pb定年结果表明斜长角闪岩形成于2486±18Ma(MSWD=1.4),而黑云斜长片麻岩形成于2507±37Ma(MSWD=2.0)。大阴山地区单塔子中低级变质表壳岩系主要由浅变质火山岩、云母石英片岩、斜长角闪岩、磁铁石英岩和大理岩等组成,SHRIMP锆石U-Pb定年结果显示,浅变质火山岩中的变玄武岩形成于2490±19Ma(MSWD=2.0),而变英安岩形成于2502±8Ma(MSWD=0.83)。因此,冀北中低级变质的表壳岩系主要形成于新太古代末期,形成年龄为2507~2486Ma;结合冀东青龙地区新太古代末期(2511~2503Ma)的浅变质火山-沉积岩系(青龙表壳岩),我们认为新太古代末期,中低级变质表壳岩系广泛分布于华北克拉通的核部和边缘地区,此套岩系覆盖在太古宙高级变质杂岩之上,代表华北克拉通化之后的稳定盖层,是克拉通化的主要标志之一。     

大兴安岭中段塔尔气地区早白垩世花岗岩成因及形成构造环境

通过大兴安岭中段塔尔气地区三道桥和桑多尔岩体全岩地球化学、年代学和锆石Lu-Hf同位素资料,探讨了其形成时代、岩石成因、源区性质等方面问题。三道桥和桑多尔岩体分别由正长花岗岩和二长花岗岩组成。LA-ICP-MS锆石U-Pb定年结果表明三道桥和桑多尔岩体均形成于早白垩世,其年龄分别为143 Ma和141 Ma。岩石地球化学分析结果显示,三道桥和桑多尔岩体均为高钾钙碱性的I型花岗岩,具有相似的岩浆源区。上述两岩体的锆石ε_(Hf)(t)值为正值(分别为4.9~9.3和2.1~8.1),具较年轻的二阶段Hf模式年龄(t_(DM2)=598~882 Ma,676~1 062 Ma),暗示其源于新元古代—显生宙期间新增生地壳物质的部分熔融。结合区域研究资料,认为三道桥和桑多尔岩体的形成与蒙古-鄂霍茨克洋闭合后的岩石圈伸展密切相关。     

钟姑地区典型矿床的重磁找矿模式

钟姑地区位于宁芜中生代火山岩盆地南段,已发现有白象山铁矿等多个大型铁矿,是长江中下游成矿带以玢岩型铁矿为主的矿集区。本文以钟姑地区实测岩石物性参数为桥梁,以实测重磁数据为基础,分析了该区白象山铁矿、钟九铁矿和云楼铁矿的重磁场分布特征,并根据典型矿区已知地质条件,运用成熟的2.5D重磁联合反演技术进行计算,从定量角度认识铁矿体深部发育形态、位置与重磁异常之间对应关系,建立钟姑地区白象山铁矿等典型矿床的重磁找矿模式,为该区进一步找矿突破提供地球物理场信息。     

北祁连地区中寒武世岩相古地理研究

现今地处造山带的北祁连地区在中寒武世时期拉张的构造背景下海域进一步拓宽。通过对北祁连地区中寒武统地层野外调查及室内分析研究,在"构造控盆、盆控相"的思路控制下,以沉积岩性、结构-构造、古生物组合等为基础,在优势相原则的指导下,共划分出浅海陆棚相、半深海相、深海盆地相3种沉积相类型。结合前人的研究资料,编制北祁连地区中寒武世的沉积构造格架及沉积盆地类型图,分析其空间配置关系。在此基础上,采用造山带岩相古地理编图方法,恢复与重建研究区的古地理面貌,编制相应的古地理图。     

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.     

Tectono-sedimentary evolution of the Miocene Oued Dayr Formation (Ghomaride Complex,Internal Domain of the Maghrebian Rif Belt,Morocco)

This work deals with sedimentological, petrographic, and structural analyses of a middle Miocene late-orogenic sedimentary cycle, denoted Oued Dayr Formation, recognized in the Rifian sector of the Maghrebian Chain (Morocco). The analyzed Formation (75 m thick) starts with 15–20 m of light colored polymict conglomerates, with minor sandstone beds, lying on the Paleozoic basement and Mesozoic cover of the Ghomaride Nappe. Facies analysis indicates a fining-upward deposition in a marine environment characterized by increasing deepening, reflecting a subsidence rate that exceeds sedimentary supply. Petrographic analysis points out that sandstones are represented by litharenites originated by erosion of recycled orogen. The conglomerates pebbles and cobbles consist of Alpine low- to high-grade metamorphic rocks as metarenites, phyllites, mylonitic quartzites, micaschists, augen gneisses deriving from the exhumed deep metamorphic basement, the overlying metasedimentary of the Sebtide Nappes and of sedimentary rocks as sandstones, jaspes, limestones, and shales deriving from the Ghomaride Nappes and their sedimentary cover. Data reveal mixed provenance indicating that the Oued Dayr Formation was fed by the Internal Nappes stack of the Maghrebian Chain. Structural analysis shows that the Oued Dayr Formation accumulated in a Thrust-Top basin, during an early extension (D₀ phase), recorded by synsedimentary normal faults within middle Langhian deposits on the rear of the Internal Nappes stack. Subsequent ductile and brittle compressional (D₁, D₂, D₃) and extensional (D₄) deformation phases occurred during and/or after the stacking, exhumation, and early unroofing of Sebtide Complex coeval with the opening of the western Mediterranean back-arc basins since middle Miocene time.