基于方差递推法的Kalman滤波在钟差预报中的应用

卫星导航系统中星载原子钟的钟差预报对于导航、定位及授时具有重要的作用。为了提高卫星钟差预报的精度,设计了一种两步确定卫星钟噪声协方差矩阵的Kalman滤波钟差预报模型。该方法首先基于Hadamard总方差确定卫星钟噪声协方差矩阵的初值,然后,使用方差递推法得到滤波过程中卫星钟的噪声协方差矩阵。使用GPS系统的星载铷钟数据进行短期预报,并与常用的二次多项式模型、灰色模型进行对比,结果表明:本文中提出的方法可以实现高精度的卫星钟差预报且预报效果优于两种常用模型,同时,该方法能够在一定程度上弥补预报误差随预报时间增加而不断变大的不足。     

航空物探勘查金属矿产时提取弱信息的解释方法

摘 要　 结构逻辑法是一种提取航空物探局部弱信息的新的解释方法。在介绍该方法原理的 基础上‚探讨了利用这种方法提取航空物探局部弱信息的效果及其在矿产资源快速勘查评价 中的作用。     

主应力轴旋转对剪切带形成的影响分析

软土的剪切带形成对于研究土的工程性质、土体稳定性评价有重要的意义。 采用可以考虑主应力轴旋转影响的土的弹塑性本构模型———关口-太田模型, 分析主应力轴旋转对剪切带形成的影响, 得出最可能发生局部剪切变形的方向同主应力轴的旋转角有关, 但主轴旋转对激发变形局部化失稳的影响不大的结论。     

Comparison of ordinary and disjunctive kriging in an Indian iron ore deposit

This paper presents the results of disjunctive kriging applied to a supergene iron ore deposit of Bailadila Range of India. Disjunctive kriging is applied firstly to compare estimates of the blocks by ordinary kriging and secondly to estimate benchwise local recoverable reserves of the orebody. Good agreement exists between block estimates by ordinary kriging and disjunctive kriging except for peripheral blocks with less borehole information. Estimation of benchwise reserves shows that the behavior of the distribution of grades is different in various benches. The study shows that disjunctive kriging can be applied successfully for estimation of local recoverable reserves in the case of a good grade hematite iron ore deposit.     

Strain distribution within a km-scale,mid-crustal shear zone: The Kuckaus Mylonite Zone,Namibia

The subvertical Kuckaus Mylonite Zone (KMZ) is a km-wide, crustal-scale, Proterozoic, dextral strike-slip shear zone in the Aus granulite terrain, SW Namibia. The KMZ was active under retrograde, amphibolite to greenschist facies conditions, and deformed felsic (and minor mafic) gneisses which had previously experienced granulite facies metamorphism during the Namaqua Orogeny. Lenses of pre- to syn-tectonic leucogranite bodies are also deformed in the shear zone. Pre-KMZ deformation (D₁) is preserved as moderately dipping gneissic foliations and tightly folded migmatitic layering. Shear strain within the KMZ is heterogeneous, and the shear zone comprises anastomosing high strain ultramylonite zones wrapping around less deformed to nearly undeformed lozenges. Strain is localized along the edge of leucogranites and between gneissic lozenges preserving D₁ migmatitic foliations. Strain localization appears controlled by pre-existing foliations, grain size, and compositional anisotropy between leucogranite and granulite. The local presence of retrograde minerals indicate that fluid infiltration occurred in places, but most ultramylonite in the KMZ is free of retrograde minerals. In particular, rock composition and D₁ fabric heterogeneity are highlighted as major contributors to the strain distribution in time and space, with deformation localization along planes of rheological contrast and along pre-existing foliations. Therefore, the spatial distribution of strain in crustal-scale ductile shear zones may be highly dependent on lithology and the orientation of pre-existing fabric elements. In addition, foliation development and grain size reduction in high strain zones further localizes strain during progressive shear, maintaining the anastomosing shear zone network established by the pre-existing heterogeneity.     

Air-sea carbon-dioxide flux estimated by eddy covariance method from a buoy observation

Precise measurements of the CO2 gas transfer across the air-sea interface provide a better understanding of the global carbon cycle.The air-sea CO2 fluxes are obtained by the eddy covariance method and the bulk method from a buoy observation in the northern Huanghai sea.The effects of buoy motion on flux calculated by the eddy covariance method are demonstrated.The research shows that a motion correction can improve the correlation coefficient between the CO2 fluxes estimated from two different levels.Without the CO2-H2 O cross-correlation correction which is termed as PKT correction,the air-sea CO2 fluxes estimated by eddy covariance method using the motion corrected data are nearly an order of magnitude larger than those estimated by the bulk method.After the CO2-H2 O cross-correlation correction,some eddy covariance CO2 fluxes indeed become closer to the bulk CO2 flux,whereas some are overcorrected which are in response to small water vapor flux.     

Integrating lysimeter drainage and eddy covariance flux measurements in a groundwater recharge model

Abstract

Field-scale water balance is difficult to characterize because controls exerted by soils and vegetation are mostly inferred from local-scale measurements with relatively small support volumes. Eddy covariance flux and lysimeters have been used to infer and evaluate field-scale water balances because they have larger footprint areas than local soil moisture measurements. This study quantifies heterogeneity of soil deep drainage (D) in four 12.5-m² repacked lysimeters, compares evapotranspiration from eddy covariance (ET_EC) and mass balance residuals of lysimeters (ETwb_Lys), and models D to estimate groundwater recharge. Variation in measured D was attributed to redirection of snowmelt infiltration and differences in lysimeter hydraulic properties caused by surface soil treatment. During the growing seasons of 2010, 2011 and 2012, ETwb_Lys (278, 289 and 269 mm, respectively) was in good agreement with ET_EC (298, 301 and 335 mm). Annual recharge estimated from modelled D was 486, 624 and 613 mm for three calendar years 2010, 2011 and 2012, respectively. In summary, lysimeter D and ET_EC can be integrated to estimate and model groundwater recharge.

Editor D. Koutsoyiannis     

Anatomy of a subduction complex: architecture of the Franciscan Complex,California, at multiple length and time scales

The Franciscan Complex of California records over 150 million years of continuous E-dipping subduction that terminated with conversion to a dextral transform plate boundary. The Franciscan comprises mélange and coherent units forming a stack of thrust nappes, with significant along-strike variability, and downward-decreasing metamorphic grade and accretion ages. The Franciscan records progressive subduction, accretion, metamorphism, and exhumation, spanning the extended period of subduction, rather than events superimposed on pre-existing stratigraphy. High-pressure (HP) metamorphic rocks lack a thermal overprint, indicating continuity of subduction from subduction initiation at ca. 165 Ma to termination at ca. 25 Ma. Accretionary periods may have alternated with episodes of subduction erosion that removed some previously accreted material, but the complex collectively reflects a net addition of material to the upper plate. Mélanges (serpentinite and siliciclastic matrix) with exotic blocks have sedimentary origins as submarine mass transport deposits, whereas mélanges formed by tectonism comprise disrupted ocean plate stratigraphy and lack exotic blocks. The former are interbedded with and grade into coherent siliciclastic units. Palaeomegathrust horizons, separating nappes accreted at different times, appear restricted to narrow zones of <100 m thickness. Exhumation of Franciscan units, both coherent and mélange, was accommodated by significant extension of the hanging wall and cross-sectional extrusion. The amount of total exhumation, as well as exhumation since subduction termination, needs to be considered when comparing Franciscan architecture to modern and ancient subduction complexes. Equal dextral separation of folded Franciscan nappes and late Cenozoic (post-subduction) units across strands of the (post-subduction) San Andreas fault system shows that the folding of nappes took place prior to subduction termination. Dextral separation of similar clastic sedimentary suites in the Franciscan and the coeval Great Valley Group forearc basin is approximately that of the San Andreas fault system, precluding major syn-subduction strike-slip displacement within the Franciscan.     

Multifractal modeling and spatial statistics

In general, the multifractal model provides more information about measurements on spatial objects than a fractal model. It also results in mathematical equations for the covariance function and semivariogram in spatial statistics which are determined primarily by the second-order mass exponent. However, these equations can be approximated by power-law relations which are comparable directly to equations based on fractal modeling. The multifractal approach is used to describe the underlying spatial structure of De Wijs 's example of zinc values from a sphalerite-bearing quartz vein near Pulacayo, Bolivia. It is shown that these data are multifractal instead of fractal, and that the second-order mass exponent (=0.979±0.011 for the example) can be used in spatial statistical analysis.     

Influence of the intermediate principal stress and principal stress direction on the mechanical behavior of cohesionless soils using the discrete element method

In this paper, the Discrete Element Method (DEM) is employed to numerically explore the response of hollow cylinder specimens of granular soils under complex stress paths. Two series of numerical tests are conducted to clarify the effects of the principal stress direction α and the intermediate principal stress through the b-value on the mechanical response of granular materials. The effects of α and b-value on the non-coaxiality of the principal stress and the principal plastic strain increment directions are investigated. It is observed that b-value and α significantly affect the non-coaxial behavior of granular materials. Finally, the results are discussed and compared with those obtained from physical laboratory tests.