Spatial-temporal characterization of the San Andreas Fault by fault-zone trapped waves at seismic experiment site,Parkfield, California

In this article, we review our previous research for spatial and temporal characterizations of the San Andreas Fault (SAF) at Parkfield, using the fault-zone trapped wave (FZTW) since the middle 1980s. Parkfield, California has been taken as a scientific seismic experimental site in the USA since the 1970s, and the SAF is the target fault to investigate earthquake physics and forecasting. More than ten types of field experiments (including seismic, geophysical, geochemical, geodetic and so on) have been carried out at this experimental site since then. In the fall of 2003, a pair of scientific wells were drilled at the San Andreas Fault Observatory at Depth (SAFOD) site; the main-hole (MH) passed a ～200-m-wide low-velocity zone (LVZ) with highly fractured rocks of the SAF at a depth of ～3.2 km below the wellhead on the ground level (Hickman et al., 2005; Zoback, 2007; Lockner et al., 2011). Borehole seismographs were installed in the SAFOD MH in 2004, which were located within the LVZ of the fault at ～3-km depth to probe the internal structure and physical properties of the SAF. On September 282004, a M6 earthquake occurred ～15 km southeast of the town of Parkfield. The data recorded in the field experiments before and after the 2004 M6 earthquake provided a unique opportunity to monitor the co-mainshock damage and post-seismic heal of the SAF associated with this strong earthquake. This retrospective review of the results from a sequence of our previous experiments at the Parkfield SAF, California, will be valuable for other researchers who are carrying out seismic experiments at the active faults to develop the community seismic wave velocity models, the fault models and the earthquake forecasting models in global seismogenic regions.     

基于MATLAB的基本矩阵随机采样鲁棒估计

基本矩阵的估算可以分为线性算法、非线性优化迭代算法、鲁棒估计算法三大类,而以鲁棒估计效率、效果最为突出,但是实现比较复杂。本文通过数码相机获取的两张未标定照片,采样互相关系数法和RANSAC算法,在MATLAB中实现了基本矩阵的自动拟合,实验表明拟合精度较高、速度快,为下一步三维场景的自动重建和量测奠定了基础。     

配置法在局部大地水准面精化中的实例研究

从最小二乘配置方法的基本原理出发,以我国某地区范围内1km分辨率的大地水准面高模型数据为例,根据实用公式计算了试验区大地水准面高的协方差值后,采用多项式函数模型和高斯函数模型分别拟合了该地区大地水准面高的局部协方差函数,并对试验区内18个检核点做了推估计算。根据推估值(Nfit)与实测值(NGPSL)的比较分析表明,虽然多项式协方差函数模型略优于高斯协方差函数模型,但它们都能以厘米级的精度拟合局部大地水准面,这表明了配置法用于精化厘米级大地水准面的有效性。     

面向离散点的空间权重矩阵生成算法与实证研究

采用阈值法和k-近邻法度量空间上离散点间的空间邻接关系,针对不同的距离计算方式(欧式距离和曼哈顿距离)设计了面向离散点的空间权重矩阵生成算法,使用C#语言在计算机上实现。用该算法对收集的8 367个常州市地价样点构建了不同土地用途地价样点的空间权重矩阵,并计算出分用途的常州市城市地价空间自相关指数。     

国产磁通门磁力仪稳定性分析

国产磁通门磁力仪(GM-3磁通门磁力仪)在全国各地磁台站推广应用,两年来,该仪器在武汉地震台工作比较稳定。本文对该仪器的基线值进行了多层线性回归分析。结果表明仪器性能比较稳定。     

应用分块矩阵求逆解算约束秩亏间接平差模型

在通常的四分块矩阵求广义逆矩阵和凯利逆矩阵公式基础上,分析左上角子矩阵为秩亏、右下角子矩阵为零的特殊四分块矩阵的凯利逆矩阵存在条件,应用广义逆矩阵法和矩阵变换法推导该类特殊四分块矩阵的凯利逆矩阵显性表达公式,并用于解算约束秩亏间接平差模型的参数估计.实验数据表明,当满足存在性条件时,应用分块矩阵求逆公式解算约束秩亏间接...     

图像重建块迭代算法中松弛参数选取的研究

在图像重建中,Landweber迭代算法是图像重建算法中的重要方法.本文将针对Landweber分块迭代算法中松弛参数的选取进行研究.在重建过程中采取对投影矩阵按投影角度分块的方法,选取特定的松弛参数.通过数值实验得出结论:对于按角度分块的块迭代算法,松弛参数选取为λ乘以块矩阵与其共轭转置矩阵乘积的最大特征值分之一,当采集完全投影数据,且λ接近(1/6)～(1/7)时效果最好.另外,本文按角度分块的做法和松弛参数的选取方法对于有限角度图像重建问题也是可行的.     

In-season crop classification using elements of the Kennaugh matrix derived from polarimetric RADARSAT-2 SAR data

Accurate spatio-temporal classification of crops is of prime importance for in-season crop monitoring. Synthetic Aperture Radar (SAR) data provides diverse physical information about crop morphology. In the present work, we propose a day-wise and a time-series approach for crop classification using full-polarimetric SAR data. In this context, the 4 × 4 real Kennaugh matrix representation of a full-polarimetric SAR data is utilized, which can provide valuable information about various morphological and dielectric attributes of a scatterer. The elements of the Kennaugh matrix are used as the parameters for the classification of crop types using the random forest and the extreme gradient boosting classifiers.The time-series approach uses data patterns throughout the whole growth period, while the day-wise approach analyzes the PolSAR data from each acquisition into a single data stack for training and validation. The main advantage of this approach is the possibility of generating an intermediate crop map, whenever a SAR acquisition is available for any particular day. Besides, the day-wise approach has the least climatic influence as compared to the time series approach. However, as time-series data retains the crop growth signature in the entire growth cycle, the classification accuracy is usually higher than the day-wise data.Within the Joint Experiment for Crop Assessment and Monitoring (JECAM) initiative, in situ measurements collected over the Canadian and Indian test sites and C-band full-polarimetric RADARSAT-2 data are used for the training and validation of the classifiers. Besides, the sensitivity of the Kennaugh matrix elements to crop morphology is apparent in this study. The overall classification accuracies of 87.75% and 80.41% are achieved for the time-series data over the Indian and Canadian test sites, respectively. However, for the day-wise data, a ∼6% decrease in the overall accuracy is observed for both the classifiers.     

Local Oceanic Vertical Deflection Determination with Gravity Data Along a Profile

Deflections of the vertical (DOVs) over oceans cannot be directly measured, which restricts their applications. A local covariance function of anomalous potential is put forward in this paper in conjunction with the least-squares collocation (LSC) method to compute the oceanic DOVs utilization of oceanic gravity data along a profile. The covariance functions of gravity field quantities have been derived directly as functions of x, y and z without the need to introduce coordinate transformations corresponding to along- or cross-profile components. In the proposed methodology, gravity data along a profile were used to calculate the residual gravity anomaly using the remove-compute-restore technique. The residual gravity anomaly was used to calculate the parameters of the proposed covariance function of the local anomalous gravity field, which was used in the LSC to compute the residual DOVs along the profile. The residual DOVs added model DOVs to recover the DOVs along the profile. The results of a simulation experiment prove that the proposed methodology is feasible and effective.     

基于经验模态分解和小波分解估算海气通量涡相关计算中的截断时间尺度

基于时长38天的海表风场实测数据,应用经验模态分解(EmpiricalModeDecomposition,EMD)和小波分解(Wavelet Decomposition, WD)这两种数据处理方法首先对涡相关法中的截断时间尺度(CutoffTimescale,CTS)进行估算,结果显示:基于EMD与WD方法估算出的CTS一般都在40秒左右(EMD的结果略小),远远小于传统涡相关法中CTS的取值(固定为10分钟),且EMD和WD的使用使得每一段数据都能够根据自身的湍流特点而获得合适的CTS; EMD方法和WD方法有效的去除了计算结果中的非湍部分,且对通量传输方向的刻画也更加合理,极大提高了通量的计算精度,所得通量与传统方法计算的通量偏差平均值高达45%;研究还对EMD和WD的优缺点进行了对比分析,结果表明EMD相比于WD有更高的自主性,而WD对信号的分离程度则更高。