基于相位一致性的遥感影像道路特征检测方法

在空间域中进行道路特征检测容易受图像亮度、对比度和噪声的影响,为了克服这一缺点,该文采用在频率域中进行道路特征检测。基于频率域的相位一致性道路特征检测方法具有亮度和对比度不变性,特征在频率域中表现出最大的相位一致性,通过计算影像6个方向的相位一致性值得到相位一致性值之和的灰度图,对灰度图进行二值化,然后使用数学形态学方法进行细化,去除噪声,检测出道路特征。实验表明,该方法对对比度、亮度变化不敏感,能够精确地检测出与人眼视觉相一致的道路特征。     

南极格罗夫山普通球粒陨石的矿物岩石学特征

本文研究了第19次中国南极科学考察队收集的97块普通球粒陨石样品的矿物岩石学特征,在此基础上划分了陨石的岩石学类型、化学群、冲击变质和风化类型。共发现5块非平衡型普通球粒陨石:GRV 020011(L3),GRV 020035(H3),GRV 020054(H3),GRV 020137(L3)和GRV020169(L3)。化学群分类的结果是H群陨石24块,L群70块,LL群3块。多数样品的冲击变质程度较低,但也有少量样品具有S5型的冲击变质特征。样品的风化程度普遍较低,仅GRV020055(H5),GRV022123(L5)和GRV022462(LL6)3块陨石达到了W3型。     

Trajectory analysis of the lower Brent Group (Jurassic), Northern North Sea: contrasting depositional patterns during the advance of a major deltaic system

Trajectory analysis is an alternative approach to systems tract analysis in unravelling the sequence stratigraphic development of sedimentary successions. Whereas the latter anticipates a succession of the depositional history in terms of a given order of systems tracts, trajectory analysis combines trajectory classes in any order, thus providing a more flexible interpretation of the depositional evolution with fewer a priori assumptions about the type or the nature of the mechanisms driving sequence development. The overall regressive part of the Brent Delta (Middle Jurassic, Northern North Sea) has been analysed using this approach. The distribution, thicknesses and stacking patterns of facies associations have been analysed to unravel the trajectorial behaviour of the system. In proximal areas (Oseberg domain), thin shoreface/foreshore packages associated with a prograding strandplain are overlain by upper delta-plain (floodplain) and distributary channel deposits. Flat or descending regressive trajectories can explain the stratigraphic development in this area. A short distance to the north (Huldra domain), the presence of thicker shoreface/foreshore packages and lower delta-plain sediments suggests a low-angle ascending regressive trajectory. In more distal areas (Gullfaks and Visund domains), a higher rate of aggradation leads to the development of even thicker shoreface/foreshore packages and the development of lagoons and bays in the lower delta-plain realm. Alternating high- and low-angle ascending regressive trajectories can explain the distal development.     

基于多重分割关联子的高分辨率遥感场景分类

高分辨率遥感影像提供了丰富的外观信息和空间结构信息,广泛应用于土地利用分类当中,源于文章领域的视觉词袋（Bag-of-Visual-Words,BoVW）模型现已成功应用于图像分类领域。传统的BoVW模型忽略了特征之间的空间布局信息和像素一致性信息,提出多重分割关联子特征,融合图像的外观信息、空间布局信息和像素一致性信息,实验表明该方法能够获取优于许多经典的遥感图像特征的性能。     

菲尔德斯半岛黑背海鸥鸣声特征初步分析

本文对南极菲尔德斯半岛黑背海鸥在某些生态行为下发出的鸣声波形结构和参数进行了分析。得出 :幼黑背海鸥发出“吱”声 ,频带在 1 6 1 0至 483 1 Hz,谱峰在 1 .7、2 .5和 2 .9k Hz;黑背海鸥发出“哦”声 ,频带在 1 1 71至 43 4 3 Hz,谱峰约在 1 .5 k Hz;黑背海鸥发出“哦啊”声 ,频带在1 0 72至 3 93 2 Hz,谱峰在 1 .2、1 .7和 3 k Hz,“咯”声 ,频带在 1 3 6 0至 41 5 3 Hz,谱峰在 1 .3、2和 2 .9k Hz     

Channel bed steps and pool shapes along Soda Creek, Three Sisters Wilderness, Oregon

Field study of bedrock step–pool systems along the upper reaches of Soda Creek in the Three Sisters Wilderness of Oregon shows strong correlation between several form variables (shape) and channel slope. Although step height and step length showed no regular spacing and variable correlation with channel slope, length to height ratios demonstrated strong negative correlations: steep slopes (20% to 80%) featured greater step height and shorter pool lengths than did flatter channel slopes. Correlations between step height to length ratios and channel slope varied between three lithologies. Explained variations ranged from 0.984 for the oldest channel steps developed in basalt, to 0.982 for steps of intermediate age developed in andesite, to 0.964 for the youngest steps developed in dacite. Sample size was 57, 40, and 33, respectively. The frequency of pool shape classes did not vary by lithology, but specific shape classes developed under differing slope conditions by rock type. All pool classes have adjusted (developed) their form to maximize resistance to flow H/L/S, and they have done so in remarkably uniform fashion.     

Polarization analysis in the wavelet domain based on the adaptive covariance method

We present an improved method for computing polarization attributes of particle motion from multicomponent seismic recordings in the time–frequency domain by using the continuous wavelet transform. This method is based on the analysis of the covariance matrix. We use an approximate analytical formula to compute the elements of the covariance matrix for a time window which is derived from an averaged instantaneous frequency of the multicomponent record. The length of the time-window is automatically and adaptively set to match the dominant period of the analysing wavelet at each time–frequency point. Then the eigenparameters are estimated for each time–frequency point without interpolation. With these key features, our method provides a suitable approach for polarization analysis of dispersive signals or overlapping seismic arrivals in multicomponent seismic data. For polarization analysis in the time domain, we show that the proposed method is consistent with existing polarization analysis methods. We apply the method to real data sets from exploration and earthquake seismology to illustrate some filtering applications and wave type characterizations.     

The correlation of microtremors: empirical limits and relations between results in frequency and time domains

The use of the correlation of microtremor records is on its way to develop into a common tool to estimate local shear wave velocity structure. For this reason, the establishment of the conditions for the correct use of this method and its limitations when applied to real data is becoming increasingly important. In addition to the use of frequency domain spatial correlation technique [the Spatial Autocorrelation (SPAC) method], the use of time domain correlation to obtain the Green's function of the medium is rapidly gaining presence. We explore the use of microtremor correlation techniques in the time domain to determine local velocity structure and compare with previous results obtained with the same data using SPAC. Our data come from three experiments carried out in Parkway and Wainuiomata valleys in New Zealand, using broad-band portable stations. Interstation distances range from 5 m to 2.1 km, and our results are useful in the frequency band from 0.1 to 7 Hz. Frequency domain correlation requires an isotropic microtremor field, a condition that need not be satisfied in the time domain. Two station correlations provide useful results due to the temporal stationarity and isotropy, in average, of the microtremor wavefield. This manifests itself in the symmetry of the temporal correlation functions with respect to zero time. Our results show that the local velocity structure and the interstation distance are the key factors conditioning the frequency range where surface wave dispersion can be correctly measured either in frequency or time domains. We confirm that, when the interstation distance becomes much larger than the dominant wavelengths, only the correlation in time domain is useful. All our results indicate that the signal obtained in the correlation of vertical component microtremors is due to the fundamental mode of Rayleigh waves, which appears as the most stable propagation mode, without any indication of body waves.     

STARTER: a statistical GIS-based model for the prediction of snow avalanche susceptibility using terrain features—application to Alta Val Badia, Italian Dolomites

This article describes a methodology to localise areas with high potential towards natural snowpack instability under particular meteorological conditions, at the scale of an Alpine valley. Localisation is based on statistically relating known release areas of past avalanche events to maps of: (1) slope inclination, (2) slope orientation (aspect), (3) elevation, (4) distance from crest lines, (5) terrain roughness and (6) concavities/convexities. The maps have been built using two different GIS softwares while the statistical analyses have been performed with a specific software handling also Fuzzy Set theory algorithms. The results of the statistical analyses have been verified on test release—areas which have not been used as input data for the statistical analyses. Verification allowed to quantify how reliably the susceptibility values were calculated, to compare the values obtained using different combinations of terrain features and to finally decide on the most efficient combination. The susceptibility maps were calculated and verified for three different meteorological scenarios (given by three classes of snow depth). Verification has shown that the accuracy of the susceptibility maps was between 67% and 82%. The three susceptibility maps show a remarkable difference in the spatial pattern of the highest susceptibility pixels suggesting that for different meteorological scenarios different classes of terrain features need to be considered.The possibility to make combinations of terrain features and to assess and verify their statistical relationship with release areas of past avalanche events is the major original step made by STARTER. Linking those release areas to meteorological scenarios is an attempt to include in the analysis the combined influence of terrain features and meteorological conditions towards snowpack instability.     

Frequency response modelling of seismic waves using finite difference time domain with phase sensitive detection (TD–PSD)

This paper describes an efficient approach for computing the frequency response of seismic waves propagating in 2- and 3-D earth models within which the magnitude and phase are required at many locations. The approach consists of running an explicit finite difference time domain (TD) code with a time harmonic source out to steady-state. The magnitudes and phases at locations in the model are computed using phase sensitive detection (PSD). PSD does not require storage of time-series (unlike a fast Fourier transform), reducing its memory requirements. Additionally, the response from multiple sources can be obtained from a single finite difference run by encoding each source with a different frequency. For 2-D models with many sources, this time domain phase sensitive detection (TD–PSD) approach has a higher arithmetic complexity than direct solution of the finite difference frequency domain (FD) equations using nested dissection re-ordering (FD–ND). The storage requirements for 2-D finite difference TD–PSD are lower than FD–ND. For 3-D finite difference models, TD–PSD has significantly lower arithmetic complexity and storage requirements than FD–ND, and therefore, may prove useful for computing the frequency response of large 3-D earth models.     

Effects of topography and crustal heterogeneities on the source estimation of LP event at Kilauea volcano

The main goal of this study is to improve the modelling of the source mechanism associated with the generation of long period (LP) signals in volcanic areas. Our intent is to evaluate the effects that detailed structural features of the volcanic models play in the generation of LP signal and the consequent retrieval of LP source characteristics. In particular, effects associated with the presence of topography and crustal heterogeneities are here studied in detail. We focus our study on a LP event observed at Kilauea volcano, Hawaii, in 2001 May. A detailed analysis of this event and its source modelling is accompanied by a set of synthetic tests, which aim to evaluate the effects of topography and the presence of low velocity shallow layers in the source region. The forward problem of Green's function generation is solved numerically following a pseudo-spectral approach, assuming different 3-D models. The inversion is done in the frequency domain and the resulting source mechanism is represented by the sum of two time-dependent terms: a full moment tensor and a single force. Synthetic tests show how characteristic velocity structures, associated with shallow sources, may be partially responsible for the generation of the observed long-lasting ringing waveforms. When applying the inversion technique to Kilauea LP data set, inversions carried out for different crustal models led to very similar source geometries, indicating a subhorizontal cracks. On the other hand, the source time function and its duration are significantly different for different models. These results support the indication of a strong influence of crustal layering on the generation of the LP signal, while the assumption of homogeneous velocity model may bring to misleading results.     

南极格罗夫山普通球粒陨石的矿物岩石学特征

本文研究了第19次中国南极科学考察队收集的97块普通球粒陨石样品的矿物岩石学特征,在此基础上划分了陨石的岩石学类型、化学群、冲击变质和风化类型。共发现5块非平衡型普通球粒陨石:GRV 020011(L3),GRV 020035(H3),GRV 020054(H3),GRV 020137(L3)和GRV020169(L3)。化学群分类的结果是H群陨石24块,L群70块,LL群3块。多数样品的冲击变质程度较低,但也有少量样品具有S5型的冲击变质特征。样品的风化程度普遍较低,仅GRV020055(H5),GRV022123(L5)和GRV022462(LL6)3块陨石达到了W3型。     

2009/2010—2011/2012年度夏秋季南奥克尼群岛水域南极磷虾捕捞群体年龄结构时空变化

利用南极磷虾渔业科学观察员收集的南极磷虾生物学数据,在将研究区域划分成60'×30'精细尺度单元的基础上,分析了南奥克尼群岛南极磷虾种群捕捞群体年龄结构的时空变化。结果表明,2009/2010年度,南极磷虾体长范围为33.2—69.0 mm,平均体长为51.8 mm,优势体长组为52—60 mm(52.3%),另在44—48 mm有一个次峰区;2010/2011年度,体长范围为22.9—59.3 mm,平均体长为44.5 mm,优势体长组为42—50 mm(56.8%);2011/2012年度,体长范围为18.5—59.4 mm,平均体长为41.7 mm,优势体长组为36—46 mm(67.3%)。三个年度间南极磷虾体长频次分布呈显著性差异(p0.001),且不具相似性(PSI=41.25—73.42)。2009/2010—2011/2012年度各旬别磷虾体长组可分成3个组,其中2011年2月上旬和2012年3月中旬可划分为一组,2010年1月下旬至2月中旬以及2011年2月中旬可归为一组,这几个体长组均呈双峰分布,其他的体长组划归为一组。1+和2+龄虾大多出现在离南奥克尼群岛最近的水域,而离南奥克尼群岛较远的单元中则多为4+龄以上的磷虾个体。     

A probabilistic approach to the inversion of data from a seismic array and its application to volcanic signals

Array techniques are particularly well‐suited for detecting and quantifying the complex seismic wavefields associated with volcanic activity such as volcanic tremor and long‐period events. The methods based on the analysis of the signal in the frequency domain, or spectral methods, have the main advantages of both resolving closely spaced sources and reducing the necessary computer time, but may severely fail in the analysis of monochromatic, non‐stationary signals. Conversely, the time‐domain methods, based on the maximization of a multichannel coherence estimate, can be applied even for short‐duration pulses. However, for both the time and the frequency domain approaches, an exhaustive definition of the errors associated with the slowness vector estimate is not yet available. Such a definition become crucial once the slowness vector estimates are used to infer source location and extent. In this work we develop a method based on a probabilistic formalism, which allows for a complete definition of the uncertainties associated with the estimate of frequency–slowness power spectra from measurement of the zero‐lag cross‐correlation. The method is based on the estimate of the theoretical frequency–slowness power spectrum, which is expressed as the convolution of the true signal slowness with the array response pattern. Using a Bayesian formalism, the a posteriori probability density function for signal slowness is expressed as the difference, in the least‐squares sense, between the model spectrum and that derived from application of the zero‐lag cross‐correlation technique. The method is tested using synthetic waveforms resembling the quasi‐monochromatic signals often associated with the volcanic activity. Examples of application to data from Stromboli volcano, Italy, allow for the estimate of source location and extent of the explosive activity.     

A World Model of Chemical Denudation in Karst Terrains

In this paper, I develop global karst chemical denudation models using chemical equilibrium equations. Theoretical karst chemical denudation rates are calculated as soil carbon dioxide concentration varies from 10^?3.5 atm to 10^?1.0 atm, temperature varies from 0°C to 30°C, and annual runoff varies from 500 mm to 3,000 mm. Both open and closed karst solution systems are examined. The Drake (1980) and the Brook‐Folkoff‐Box (1983) soil pCO² equations are used to develop chemical denudation models for different carbonate rocks (limestone and dolomite), climate (tropical, temperate, and arctic/alpine terrains), and karst solution type (open and closed systems). The major conclusion is that the karst solution type, least known in the past karst studies, is an important factor in controlling chemical denudation rates.     

A World Model of Chemical Denudation in Karst Terrains

In this paper, I develop global karst chemical denudation models using chemical equilibrium equations. Theoretical karst chemical denudation rates are calculated as soil carbon dioxide concentration varies from 10^?3.5 atm to 10^?1.0 atm, temperature varies from 0°C to 30°C, and annual runoff varies from 500 mm to 3,000 mm. Both open and closed karst solution systems are examined. The Drake (1980) and the Brook-Folkoff-Box (1983) soil pCO² equations are used to develop chemical denudation models for different carbonate rocks (limestone and dolomite), climate (tropical, temperate, and arctic/alpine terrains), and karst solution type (open and closed systems). The major conclusion is that the karst solution type, least known in the past karst studies, is an important factor in controlling chemical denudation rates.