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.     

变化环境下渭河流域水文干旱演变特征剖析

环境变化影响区域水资源的可持续开发利用,导致水文过程出现非平稳特征,需发展非平稳水文干旱评估方法。选取渭河流域为研究区,依据流域内2个水文站、62个雨量站和24个气象站1961-2013年数据,基于可变下渗容量模型定量分离气候变化和人类活动对径流衰减的贡献;采用标准化径流指数（Standardized Runoff Index, SRI）剖析水文干旱时空演变特征;提出多种SRI参数化方案,对比评定各方案表征非平稳干旱的合理性以及环境变化对干旱演变的影响作用。结果表明:自1991年以来渭河流域年径流量呈显著衰减趋势,人类活动是径流演变的主要因素,对咸阳和华县站径流量变化的贡献率分别为-66.7%和-71.0%;时变参数方案计算的干旱指数能合理重建历史水文干旱序列;人类活动是渭河流域1991年以来短历时水文干旱发生的主导因素,气候变化主要影响长历时旱涝的演变趋势。     

Study on the Global Carbon Assimilation System Based on Multisource Remote Sensing Data

The Paris agreement signed in April, 2016 aims to balance global anthropogenic carbon emissions and terrestrial carbon sinks by the middle of the 21^st century. To fulfill this goal, it is necessary to calculate carbon fluxes of different regions reliably. The global carbon assimilation system is an effective technique for achieving this goal. The Ministry of Science and Technology of China supports the project entitled as study on the global carbon assimilation system based on multisource remote sensing data through the national key research and development programs for global change and adaptation during the thirteen-five period. This project will develop synergic inversion techniques for retrieving key parameters of biological and atmospheric cycles and for assimilating multisource remote sensing and ground based data. Then, the high resolution global carbon assimilation system coupled with an ecological model will be constructed. This system is able to assimilate jointly multisource observation data and to optimize key model parameters, photosynthesis and respiration carbon fluxes of global terrestrial ecosystems, and anthropogenic carbon emission fluxes of key regions. This system will be used to study quantitatively the spatial and temporal patterns of carbon fluxes of global terrestrial ecosystems and anthropogenic carbon emission fluxes of key regions and to identify the mechanisms driving the global terrestrial carbon sinks and sources. The outputs of this study will be helpful for the fulfillment of the key research and development programs for global change and adaptation and provide valuable data and technical support for the decision-making in China.     

伊通盆地鹿乡断陷C区块含油饱和度与含油性分析

综合研究伊通盆地鹿乡断陷双二段储层的岩性、含油性和含油饱和度,确定了适用于本区的含油饱和度计算方法,并与实际含油性进行了对比,分析了沉积和构造对含油性的影响。研究区双二段含油储层的岩性以粉砂岩、细砂岩和砂砾岩为主,含油级别为荧光级以上。根据阿尔奇公式计算出的含油饱和度大部分与实际的含油性相符合,并将含油饱和度70%、50%和38%值作为本区气层、油层、油水层和水层的区分标准。本区储层的含油性受到沉积和构造两方面的控制,沉积作用主要控制储层的岩石结构和物性,构造控制油气的聚集,即在构造高部位,实际含油饱和度比计算结果要高,而在构造低部位,实际含油饱和度比计算结果要低。     

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.     

Sensitivity analysis for the GALDIT method based on the assessment of vulnerability to pollution in the northern Sfax coastal aquifer,Tunisia

Coastal aquifer of northern Sfax (Tunisia) suffers from the high risk to seawater intrusion and the water quality degradation due to the overexploitation. Hence, assessing the study area vulnerability to pollution is highly crucial so as to protect the groundwater resources. The assessment has been performed by applying the GALDIT method using Geographic Information System (GIS) software and multi-criteria evaluation techniques, and the sensitivity analysis approach to evaluate the effect of each GALDIT parameter on the vulnerability assessment. The GALDIT vulnerability map classifies the study area into three vulnerability classes: low vulnerability (30–50), moderate vulnerability (50–70), and high vulnerability (70–90), which represent 5, 30, and 65 % of the study area, respectively. The map illustrates that the coastal zones of the aquifer are the most threatened areas. The sensitivity analysis results show that the aquifer hydraulic conductivity (A) and the thickness of the aquifer (T) represent the determining factors in the modified vulnerability model. The real weight was used to elaborate the modified GALDIT model which was correlated with resistivity values for validation. This study could serve as a scientific basis for sustainable land planning and groundwater management in the study area.