P-T-t-deformation paths recorded by kinzigites during diapirism in the western Variscan belt (Golfe du Morbihan, southern Brittany, France)

Abstract The high-grade metamorphic rocks of southern Brittany underwent a complex tectonic evolution under various P-T conditions (high-P, high-T), related to stacking of nappes during Palaeozoic continentcontinent collision. The east to west thrusting observed in the whole belt is strongly perturbed by vertical movements attributed to the ascent of anatectic granites in the high-T area. The field reconstruction of subvertical, closed elliptical structures in gneisses and migmatites, associated with the subhorizontal, doubly radial pattern of stretching lineation in the mica schists, suggests the existence of an elliptical diapiric body buried at depth beneath the present erosion level. Deformation is associated with a complex P-T evolution partly recorded in aluminous gneisses (kinzigites, e.g. morbihanites). A chronology of successive episodes of mineral growth at different compositions is established by detailed studies of the mineral-microstructure relationships in X-Z sections, using the deformation-partitioning concept (low- and high-strain zones). Several thermometric and barometric calibrations are applied to mineral pairs either in contact or not in contact but in equivalent microstructiiral positions with respect to the deformation history. This methodology provides a continuous microstructural control of P-T variations through time and leads to three P-T-t-d paths constructed from numerous successive P-T estimations. Path 1 is a clockwise retrograde path preserved in low-strain zones, which records general exhumation movements after crustal thickening. Paths 2 and 3 are clockwise prograde/retrograde paths from high-strain zones; they are interpreted and discussed in the light of models of crustal anatexis and upward movement of magma (diapirism). Deformation and P-T effects induced by diapirism can be distinguished from the general deformation-metamorphic history of a belt, and would seem to be produced during a late stage of its history. The present microstructural-petrological approach to defining successive mineral equilibria in relation to progressive deformation steps provides a far more accurate evaluation of the metamorphic evolution than is possible by ‘standard’thermobarometry.     

Parameter identification of inelastic structures under dynamic loads

The generalized model of differential hysteresis contains 13 control parameters with which it can curve‐fit practically any hysteretic trace. Three identification algorithms are developed to estimate the control parameters of hysteresis for different classes of inelastic structures. These algorithms are based upon the simplex, extended Kalman filter, and generalized reduced gradient methods. Novel techniques such as global search and internal constraints are incorporated to facilitate convergence and stability. Effectiveness of the devised algorithms is demonstrated through simulations of two inelastic systems with both pinching and degradation characteristics in their hysteretic traces. Owing to very modest computing requirements, these identification algorithms may become acceptable as a design tool for mapping the hysteretic traces of inelastic structures. Copyright © 2002 John Wiley & Sons, Ltd.     

锚固站对BDS自主导航精度的影响分析

介绍了自主导航的轨道确定及时间同步观测方程。以北斗仿真全星座为对象,通过采用仿真星间及卫星与地面锚固站间观测值,进行了60 d自主导航解算,分别探讨了锚固站数量及锚固站观测连续性对北斗卫星导航系统（BDS）3类卫星自主导航精度的影响。结果表明：锚固站数量及观测连续性对RERR及CERR无影响;加入1个锚固站即可显著改进URE结果精度,继续增加锚固站数量虽然可进一步提高URE精度但其改进效果有限;锚固站观测中断时间越长,其对应自主导航精度越低。因此,在BDS自主导航运行模式下应保持较高的锚固站观测频次以保证自主导航精度;另外,锚固站数量及观测连续性对北斗系统3类卫星自主导航精度的影响无显著差异。     

集合卡尔曼滤波同化中雷达位置的敏感性研究

针对对流尺度集合卡尔曼滤波（EnKF）雷达资料同化中雷达位置对同化的影响进行研究。为了考察强对流出现在雷达不同方位时集合卡尔曼滤波同化雷达资料的能力,以一个理想风暴为例,设计了8个均匀分布在模拟区域周围的模拟雷达进行试验。单雷达同化试验中,初期同化对雷达位置较敏感,而十几个循环后对雷达方位的敏感性降低。造成初期同化效果较差的雷达观测位于模拟区域正南和正北方向,这两部雷达与模拟区域中心的连线垂直于风暴移动方向（即环境气流的方向）。双雷达试验的结果表明,正东、正南、正西和正北方向的雷达组合观测会使同化初期误差较大,这说明并不是所有与风暴连线成90°的雷达组合都能在短时同化中得到合理的分析结果,还需要都处于模拟区域对角线上（即与环境气流成45°夹角）,同化效果才较好。短时同化后的确定性预报结果表明,较大分析误差也会导致较大预报误差。这些分析误差主要是由于同化初期不准确的集合平均场驱动出的不合理的背景误差协方差造成的。当背景场随着同化循环得到改进后,驱动出的合理的背景误差协方差使得不同位置雷达同化造成的差异逐步减小。基于上述结果,引入迭代集合均方根滤波（iEnSRF）算法,结果显示使用该算法后,雷达位置对同化效果的影响减小,同化不同位置的雷达资料均能有效降低分析和预报误差。      

利用地震噪声准实时监测短周期面波波速变化

依据由噪声信号提取面波格林函数的原理,利用福建省地震台网25个宽频带台站2007年7月2日—8月29日的脉动观测资料,得到了瑞雷波群速度在福建地区的分布,并将该结果作为下一步相对变化动态成像的背景速度分布．分布结果表明,福建地区周期为3—5s的瑞雷波群速度大致在2.9—3.1km/s之间,平均速度为3.0km/s,瑞雷波群速度分布呈北高南低的现象,这与福建地区北部多山、南部多平原盆地的地理环境有很好的吻合．而且,该群速度分布图在漳州盆地地区表现出一个非常明显的低速,这主要是因为受到盆地沉积层的影响．通过滑动窗（窗长为20d,步长1d）技术得到了观测区内周期约为3-5s的瑞雷波波速分布变化的连续图像;再利用扣除背景影响的技术,得到了2007年8月14日—2008年7月1日福建地区瑞雷波波速的相对变化时空动态图像．通过分析相对变化时空动态图像与该时间范围内发生的地震的对应关系,表明福建地区瑞雷波波速在多次网内中等强度的地震或震群（ML＞3.0）中均表现出震前波速升高,震后下降恢复的变化趋势．初步分析认为,这可能与震前整个地区受到的应力增大震后应力得到释放所导致的介质变化有关．      

快速有效的转换波共转换点叠加技术

转换波共转换点(Common Conversion Point简称CCP)叠加的关键在于CCP抽道集和非双曲线正常时差校正NMO(Normal MoveOut). 目前方法的精度限制了其在中-浅层或大炮检距情况下的应用. 我们对CCP叠加技术进行了系统研究,导出了新的CCP位置计算公式和非双曲线时距关系式,并给出了具体的CCP抽道集方法. 理论模型试验和实验资料处理表明本文方法精度高、简便易行,特别是对于中-浅层和大炮检距情况也能得到良好的叠加效果.     

基于NDVI时间序列数据的江西省水稻种植制度变化研究

利用Savitzky-Golay滤波对覆盖江西省范围的SPOT VGT NDVI时间序列数据进行平滑处理的基础上,结合坡度数据,通过非监督分类的方法提取了江西省2000、2005和2010年水稻种植范围,并根据NDVI的年内动态变化,从水稻种植范围、水稻生长季起始时间、水稻复种指数和NDVI最大振幅等分析了江西省水稻种植和生长情况,探讨2000~2010年江西省水稻生产的变化。     

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

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

Underwater Navigation Based on Real-Time Simultaneous Sound Speed Profile Correction

Precise sound speed profile (SSP) information is critical when using sonar for underwater terrain navigation. Nevertheless, acquiring SSP information in real-time is difficult, especially in underwater navigation environment. To account for this problem, this paper presents an underwater navigation method that applies real-time SSP correction to the terrain measurements. The method uses a probe to measure the surface sound speed simultaneously. Meanwhile, redundancy of topographic measurement data is exploited to derive the equivalent SSP information. The acquired equivalent SSP is updated continuously with particle filter algorithm. The terrain measurements can be corrected in real-time with the equivalent SSP to improve the performance of underwater terrain navigation. By removing the impact of inaccurate SSP from the terrain measurements, the proposed method can achieve precise and robust underwater navigation result without using an underway-profiling instrument. Simulated results confirm the good performance of the proposed method.