The diffuse-to-global and diffuse-to-direct-beam spectral irradiance ratios as turbidity indexes in an urban environment

Continuous measurements of solar spectral radiation using the Multi-filter Rotating Shadow Band Radiometer (MFRSR) are performed at the Actinometric Station of the National Observatory of Athens (ASNOA). The present study utilizes three clear-sky days of continuous observations, from local sunrise to local sunset, in order to investigate the daily variation of the radiation components (diffuse, global and direct-beam) as well as their ratios (diffuse-to-global, DGR, and diffuse-to-direct-beam, DDR) under different atmospheric conditions. Both ratios have received great scientific interest, especially for investigating solar irradiance modifications under various atmospheric conditions, aerosol load and optical properties. Apart from this, the present study shows that the DDR can constitute a measure of atmospheric turbidity when it is determined at longer wavelengths, while the DGR cannot. The effect of the solar zenith angle (SZA) on both ratios is significant at the shorter wavelengths with varying sensitivities depending on the aerosol field and sun elevation. The present study confirms the results obtained by previous solar irradiance measurements in Athens and also those computed via radiative transfer codes and sheds light on the scientific knowledge of the use of spectral DDR as an atmospheric turbidity index.     

The reference spectral noise ratio method to evaluate the seismic response of a site

The seismic response of several sites in the Durban area of South Africa has been investigated using a refined version of the spectral ratio (H/V) generated by cultural seismic noise. Two samples of noise separated by an interval of several minutes for the vertical and the two orthogonal horizontal components were taken at each site. Though the two samples differed noticeably in the time as well as in the frequency domain, when the components are considered separately, the ratios of the horizontal to the vertical spectra were very similar for most of the pairs. This similitude was obvious not only in the frequencies of the peaks, but also in the range of amplifications involved. These observations have been used to introduce a refinement to the Nakamura Method. This refinement consists in using as a reference the H/V spectral ratio of a well-known, hard rock site. All the other soft sites of interest are then compared with this standard reference site. This comparison provides valuable information on the frequencies and amplification factors introduced by the critical soft sediments. This refined method will be referred here as the ‘Reference spectral noise ratio’.This method, is compared to a similar one devised to investigate the response of horizontal layers to longitudinal seismic waves, using the ratio of the spectra generated by long- period waves.     

A data variance technique for automated despiking of magnetotelluric data with a remote reference

The magnetotelluric method employs co‐located surface measurements of electric and magnetic fields to infer the local electrical structure of the earth. The frequency dependent ‘apparent resistivity’ curves can be inaccurate at long periods if input data are contaminated – even when robust remote reference techniques are employed. Data despiking prior to processing can result in significantly more reliable estimates of long period apparent resistivities. This paper outlines a two‐step method of automatic identification and replacement for spike‐like contamination of magnetotelluric data; based on the simultaneity of natural electric and magnetic field variations at distant sites. This simultaneity is exploited both to identify windows in time when the array data are compromised as well as to generate synthetic data that replace observed transient noise spikes. In the first step windows in data time series that contain spikes are identified according to an intersite comparison of channel ‘activity’– such as the variance of differenced data within each window. In the second step, plausible data for replacement of flagged windows are calculated by Wiener filtering coincident data in clean channels. The Wiener filters – which express the time‐domain relationship between various array channels – are computed using an uncontaminated segment of array training data. Examples are shown where the algorithm is applied to artificially contaminated data and to real field data. In both cases all spikes are successfully identified. In the case of implanted artificial noise, the synthetic replacement time series are very similar to the original recording. In all cases, apparent resistivity and phase curves obtained by processing the despiked data are much improved over curves obtained from raw data.     

分频混色技术在高精度地震解释中的应用

分频解释技术已经越来越多的应用到地震储层预测和油藏描述中,成为地震解释的一种常规技术手段,分频解释技术中分频RGB颜色融合显示方法因其独有的视觉效果也开始受到关注.本文介绍了分频混色技术在高精度地震解释中的应用以及在应用过程中应注意的几点问题,包括优选高精度频谱分解方法、RGB颜色融合技术应用中频率选择的原则等.在此基础上,进一步探讨了分频数据RGB混色方法能有效解决地震沉积学中地层切片精度不满足薄储层预测的问题.RGB颜色融合显示方法能充分发掘高精度地震数据宽频的优势,得到精度更高的储层预测结果.     

Comparison of spectral ratio techniques for estimation of site effects using microtremor data and earthquake motions recorded at the surface and in boreholes

The authors apply soil response estimation techniques employing accelerograms for fifteen earthquakes recorded at the Yokohama Strong Motion Network and its vertical array of nine sites, plus microtremor data recorded at all 150 sites. Assessment of the reliability of surface to reference site spectral ratios, and horizontal to vertical spectral ratios of S‐waves, coda, and microtremors, relied on cross‐validation with both surface to borehole spectral ratios, and theoretical soil response functions for vertically incident SH plane waves calculated from the vertical array recordings and logging data, respectively. The results indicate the following: (1) the reference site method provides more accurate soil response estimations; (2) the H/V of microtremors provides relatively better soil response estimations than its S‐wave and coda counterparts; and (3) the H/V of microtremors provides coincident soil response estimations at two‐thirds of the sites when considering no more than ±30% mismatch in the determination of the predominant periods. The matching of parameters obtained from the analysis of microtremors, combined with their relatively low economic cost and the possibility of recording without strict spatial or time restrictions turns microtremors into a particularly appealing approach. Copyright © 2003 John Wiley & Sons, Ltd.     

Bayesian spectral density approach for modal updating using ambient data

The problem of identification of the modal parameters of a structural model using measured ambient response time histories is addressed. A Bayesian spectral density approach (BSDA) for modal updating is presented which uses the statistical properties of a spectral density estimator to obtain not only the optimal values of the updated modal parameters but also their associated uncertainties by calculating the posterior joint probability distribution of these parameters. Calculation of the uncertainties of the identified modal parameters is very important if one plans to proceed with the updating of a theoretical finite element model based on modal estimates. It is found that the updated PDF of the modal parameters can be well approximated by a Gaussian distribution centred at the optimal parameters at which the posterior PDF is maximized. Examples using simulated data are presented to illustrate the proposed method. Copyright © 2001 John Wiley & Sons, Ltd.     

最小曲率位场分离方法研究

位场分离是位场数据处理和解释中的重点和难点之一.本文给出了单步长非原位和原位两种最小曲率位场分离差分迭代格式,并利用Fourier频谱分析理论研究了这两种迭代格式的收敛性.通过研究表明,单步长非原位迭代格式不收敛,只有单步长原位迭代格式收敛,但单步长原位迭代格式受迭代方向选择的影响,随着迭代次数的增大其影响逐渐消失.根据单步长非原位迭代格式的频谱特点,提出了叠加步长非原位和原位迭代格式,同样利用Fourier频谱分析理论研究了叠加步长非原位和原位迭代格式的收敛性.通过研究认为,一维叠加步长非原位迭代格式收敛,但二维叠加步长非原位迭代格式不收敛;不论是一维或二维,其原位迭代格式均收敛.进一步的理论研究表明,非原位迭代格式的频率响应是一个实偶函数,而原位迭代格式的频率响应是一个复函数;单步长迭代格式的频率响应具有一定的周期性,而叠加步长迭代格式的频率响应无周期性特征;叠加步长迭代格式比单步长迭代格式的收敛性好.     

Comparison of spectral and spatial methods for a Moho recovery from gravity and vertical gravity-gradient data

In global studies investigating the Earth’s lithospheric structure, the spectral expressions for the gravimetric forward and inverse modeling of the global gravitational and crustal structure models are preferably used, because of their numerical efficiency. In regional studies, the applied numerical schemes typically utilize the expressions in spatial form. Since the gravity-gradient observations have a more localized support than the gravity measurements, the gravity-gradient data (such as products from the Gravity field and steady-state Ocean Circulation Explorer - GOCE - gravity-gradiometry satellite mission) could preferably be used in regional studies, because of reducing significantly the spatial data-coverage required for a regional inversion or interpretation. In this study, we investigate this aspect in context of a regional Moho recovery. In particular, we compare the numerical performance of solving the Vening Meinesz-Moritz’s (VMM) inverse problem of isostasy in spectral and spatial domains from the gravity and (vertical) gravity-gradient data. We demonstrate that the VMM spectral solutions from the gravity and gravity-gradient data are (almost) the same, while the VMM spatial solutions differ from the corresponding spectral solutions, especially when using the gravity-gradient data. The validation of the VMM solutions, however, reveals that the VMM spatial solution from the gravity-gradient data has a slightly better agreement with seismic models. A more detailed numerical analysis shows that the VMM spatial solution formulated for the gravity gradient is very sensitive to horizontal spatial variations of the vertical gravity gradient, especially in vicinity of the computation point. Consequently, this solution provides better results in regions with a relatively well-known crustal structure, while suppressing errors caused by crustal model uncertainties from distant zones. Based on these findings we argue that the gravity-gradient data are more suitable than the gravity data for a regional Moho recovery.     

测震台网数据传输补数技术

针对测震台网波形数据传输中断导致丢数现象,开发测震台网自动补数系统SDARS,自动检测数据丢失现象,并从地震台数据采集器下载缺失数据;设计手动辅助补数工具软件SeisDataUtil,提供不同格式的地震数据文件的手动补数功能。通过此套补数系统,有效提高测震台网波形数据完整率,最小化减少波形数据丢失。     

多分量地震资料处理解释技术研究

随着油田勘探开发的不断深入,多波多分量地震资料的应用也越来越广泛.本文结合实际资料,对多波多分量地震资料的处理和解释方法进行了探讨和研究,主要针对多波资料的静校正、叠加成像、转换波资料的共转换点的求取方法等技术进行了分析研究,开发了综合利用多波资料进行地震属性和各向异性研究的方法,形成了一套较完善的多分量地震资料的处理、解释技术系列.经过在胜利油田罗家地区的初步应用,表明了预测结果与钻井资料相吻合.     

Kernel estimation as a basic tool for geomorphological data analysis

Kernel estimation, based on the convolution of a probability density function with a set of magnitudes or event dates, provides tuneable smooth pictures of probability density functions and event intensity functions. Such pictures are in several respects superior to those provided by histograms, box plots, cumulative distributions or raw plots. They permit examination of broad features and fine structure, are readily produced with modest computational effort and are essentially free of artefacts arising from binning. Examples are given using data on cirque lengths, limestone pavements, glacier areas and dated flood deposits. The technique deserves widespread use in geomorphology and allied sciences. Copyright © 2007 John Wiley & Sons, Ltd.     

低信噪比转换波地震资料静校正(英文)

Converted waves have slow velocity and low signal-to-noise ratio. It is also difficult to pick first-breaks and bin the common-conversion-points （CCP）. Some statics methods, which work well for P-wave data, can＇t be effectively used for solving convertedwave statics problems. This has become the main obstacle to breakthroughs in convertedwave data processing. To improve converted-wave static corrections, first, a statics method based on the common-receiver-point （CRP） stack is used for the initial receiver static correction to enhance the coherency of the CRP stack. Second, a stack-power-maximization static correction which improves the continuity of the CCP stack is used for detailed receiver statics. Finally, a non-surface-consistent residual moveout correction of the CCP gathers is used to enhance the stack power of reflection signals from different depths. Converted-wave statics are solved by the joint use of the three correction methods.     

Laboratory technique for measurement of spectral induced polarization response of soil sampies1

Laboratory measurements of soil samples are necessary to assess the effect of mineralogy, grain size distribution, moisture content, and electrolyte composition on the resistivity spectrum of soil material. Laboratory results are also required for the interpretation of field data. Induced polarization phenomena in glacial soils are poorly understood and so far no convenient laboratory techniques are available for its measurement. Coarse grain size and the need to measure unsaturated samples and to monitor the homogeneity of the sample require a sample holder–electrode construction that differs from those presented in Clay mineral Studies. This study presents a spectral induced polarization laboratory system that is suitable for measuring fine- and coarse-grained and both saturated and unsaturated soil samples. The noise caused by the electrode–electrolyte interface is studied in detail. It is shown that easy-to-use platinum or acid-free steel potential electrodes are convenient over a broad frequency band ranging from 0.016Hz up to more than 1000 Hz. The laboratory experiments and comparisons between laboratory and field results also indicate that sampling and sample packing procedures have only a minor influence on the phase spectrum of glacial soils.     

Evapotranspiration modelling from climatic data using a neural computing technique

Evapotranspiration is one of the basic components of the hydrologic cycle and essential for estimating irrigation water requirements. This paper investigates the modelling of evapotranspiration using the feed‐forward artificial neural network (ANN) technique with the Levenberg–Marquardt (LM) training algorithm. The LM algorithm has never been used in evapotranspiration estimation before. The LM is used for the optimization of network weights, since this algorithm is more powerful and faster than the conventional gradient descent. Various combinations of daily climatic data, i.e. wind speed, air temperature, relative humidity and solar radiation, from three stations in Los Angeles, USA, are used as inputs to the ANN so as to evaluate the degree of effect of each of these variables on evapotranspiration. A comparison is made between the estimates provided by the ANN and those of the following empirical models: Penman, Hargreaves, Turc. Mean square error, mean absolute error and determination coefficient statistics are used as comparing criteria for the evaluation of the models' performances. Based on the comparisons, it was found that the neural computing technique could be employed successfully in modelling evapotranspiration process from the available climatic data. The results also indicate that the Hargreaves method provides better performance than the Penman and Turc methods in estimation of the evapotranspiration. The accuracy of the ANN technique in evapotranspiration estimation using nearby station data was also investigated. Copyright © 2007 John Wiley & Sons, Ltd.     

Horizontal to vertical spectral ratio and geological conditions: an analysis of strong motion data from Greece and Taiwan (SMART-1)

Nowadays, several research efforts have been directed towards the evaluation of the horizontal to vertical, (h/v), spectral ratio technique, following the pioneering work of Nakamura originally proposed for microtremor measurements. In the present study, an extension of this technique based on Greek and Taiwanese (SMART-1) strong motion data is attempted. Whenever possible comparisons with results obtained by the standard spectral ratio technique for the same sites are made. Regarding the Greek data it is concluded that (h/v) spectral ratio clearly distinguishes the range of resonant peaks between ‘rock’ and ‘alluvium’ sites. On the other hand, nonlinear effects that take place beneath SMART-1 array most probably obscure the expected consistency between (h/v) and standard spectral ratio technique.     

Gravity profiling as a technique for determining the thickness of glacier ice

Summary Detailed gravity measurements recently carried out on the Gorner glacier, Switzerland, are used to determine the variation of thickness across the glacier ice. The Gorner glacier was chosen as a test site because seismic control was available. The glacier ice at a profile near the Monte Rosa massif is associated with a relative gravity low of about –23 mgal. Model oalculations yield a corresponding ice thickness of about 400 m at the central part of the profile. A comparison of the derived residual gravity anomaly with the calculated effect of the 3-D ice model based on seismic information is made. It is shown that the regional field determined for the Gorner glacier is appropriate and gives the correct residual anomaly associated with the glacier ice. Therefore, the proposed gravity technique for determining variations of the thickness of glacier ice appears to be a valuable and rather inexpensive method for surveying glaciers.Institut für Geophysik, ETH-Zürich, Contribution No 145.     

叠前提高西部山地资料信噪比

中国西部山地地震资料非常复杂,地表静校正复杂和低信噪比是困扰处理人员的两大难题。本文在分析对比各种压制倾斜干扰的技术方法的基础上,提出采用静校正和叠前提高信噪比交互迭代的方法,有效地解决了这两个难题。对于各种干扰进行有针对性地压制与消除;在提高信噪比的基础上提高了资料的分辨率。