Evaluating empirical regression equations for Vs and estimating Vs30 in northeastern Taiwan

The shear wave velocity of shallow sediments is very important in seismic wave amplification and thus Vs30 is a well-known parameter for site classification. Based on the relationship between Vs and soil indexes, the empirical regression equations were evaluated using more than 600 data sets from the Ilan area and the Taipei Basin. Multivariable analysis, which can increase the accuracy of regression equations, was used in this study. Three extrapolations were compared, and then the most accurate bottom-constant extrapolation was adopted to estimate Vs30 at 16 boreholes, none of which reached a depth of 30 m. Ultimately Vs30 was derived for 110 free-field strong motion stations, and the stations in northeastern Taiwan were reclassified according to the Vs30-based NEHRP provisions. Regression equations of Vs and the extrapolation of Vs30 were also applied to boreholes with only N-values and corresponding depths of less than 30 m for assessing the Vs30 and NEHRP class.     

居里等温面深度的探测方法

地壳居里等温面的研究对确定深部磁体结构和了解深部热状态及结构有重要意义.本文对探测居里等温面深度的方法进行了科学地总结,分类和评估,并主要介绍了各种方法的基本原理.较直接的探测居里等温面的方法是根据大地热流值,地温\|深度关系或岩石矿物实验来确定,但是这种方法在实际中并不常用.最常用的方法是应用航磁和卫星磁数据.根据所研究的目标的不同,磁法又可以被分为两大类：单体磁异常法和组合磁异常法.这两大类下又可以有各种不同的应用技术,其中功率谱分析方法最为流行.虽然没有任何一种方法是完美的,针对组合磁异常的统计功率谱分析方法被认为是较成熟的一项技术,并得到了很多的应用及讨论.组合磁异常的功率谱分析方法的理论基础是Spector\|Grant 统计模型,这个统计模型很适用于对区域磁异常组合的分析.为了提高分辨率,无论单体磁异常法或组合磁异常法都需要考虑很大的一片研究区域,这样的结果会降低所获得的居里等温面的横向分辨率,由此所获得的居里等温面深度通常只是某一区域平均居里点深度的一级近似.一项比较特别的技术是Mayhew的最初为分析卫星磁数据而建立起来的等层模型方法.这一方法的基本原理是,根据一些校准点,将一等厚度层上所观测到的不同的磁场强度转换为等磁场强度的深度差异.这一技术与功率谱分析方法的一个重要差别在于它不需要对磁体顶部和底部所引起的功率谱进行区分,而这种区分是功率谱分析方法中所必需的.由于对居里等温面深度的探测是一种反演技术并包含了反演中的非唯一性和数学上的不稳定性,在各种方法的应用中应采取谨慎的态度.近年来,自相似和分形理论的应用为这一领域的研究提供了很多崭新的研究机会和科学思想.自相似和分形模型是对Spector\|Grant统计模型的完善和发展,因为自相似模型考虑了更真实的磁场强度差异的分形分布.另外,这一新理论也揭示了以前很多种方法中的局限性.     

An attempt to define Curie point depths in Greece from aeromagnetic and heat flow data

The objective of this study is to understand the nature and extent of the regional geothermal system at depth beneath the area of Greece by constructing the Curie isotherms.Spectral analysis of aeromagnetic data in conjunction with heat flow information revealed an almost inverse linear relation between heat flow and Curie depths and was used to construct the Curie isotherms from the existing heat flow data.The results showed that Curie depths in the area range from about 20 km in western Greece, up to 1 km beneath the Hellenic volcanic arc. These results are consistent with the existing geothermal and geotectonic regime in the area.     

A temperature prediction-correction method for estimating surface soil heat flux from soil temperature and moisture data

Surface soil heat flux is a component of surface energy budget and its estimation is needed in land-atmosphere interaction studies. This paper develops a new simple method to estimate soil heat flux from soil temperature and moisture observations. It gives soil temperature profile with the thermal diffusion equation and, then, adjusts the temperature profile with differences between observed and computed soil temperatures. The soil flux is obtained through integrating the soil temperature profile. Compared with previous methods, the new method does not require accurate thermal conductivity. Case studies based on observations, synthetic data, and sensitivity analyses show that the new method is preferable and the results obtained with it are not sensitive to the availability of temperature data in the topsoil. In addition, we pointed out that the soil heat flux measured with a heat-plate can be quite erroneous in magnitude though its phase is accurate.     

The use of mixing models and chemical geothermometers for estimating underground temperatures in geothermal systems

Application of various chemical geothermometers and mixing models indicate underground temperatures of 260°C, 280°C and 265°C in the Geysir, Hveravellir and Landmannalaugar geothermal fields in Iceland, respectively. Mixing of the hot water with cold water occurs in the upflow zones of all these geothermal systems. Linear relations between chloride, boron and δ¹⁸O constitute the main evidence for mixing, which is further substantiated by chloride, silica and sulphate relations in the Geysir and Hveravellir fields.A new carbonate-silica mixing model is proposed which is useful in distinguishing boiled and non-boiled geothermal waters. This model can also be used to estimate underground temperatures using data from warm springs. This model, as well as the chloride-enthalpy model and the Na-Li, and CO₂-gas geothermometers, invariably yield similar results as the quartz geothermometer sometimes also does. By contrast, the Na-K and the Na-K-Ca geothermometers yield low values in the case of boiling hot springs, largely due to loss of potassium from solution in the upflow. The results of these geothermometers are unreliable for mixed waters due to leaching subsequent to mixing.     

Curie point depth from spectral analysis of aeromagnetic data from Cerro Prieto geothermal area,Baja California,México

Using aeromagnetic data acquired in the area from the Cerro Prieto geothermal field, we estimated the depth to the Curie point isotherm, interpreted as the base of the magnetic sources, following statistical spectral-based techniques. According to our results the Curie point isotherm is located at a depths ranging from 14 to 17 km. Our result is somewhat deeper than that obtained previously based only in 2-D and 3-D forward modeling of previous low-quality data. However, our results are supported by independent information comprising geothermal gradients, seismicity distribution in the crust, and gravity determined crustal thickness. Our results imply a high thermal gradient (ranging between 33 and 38 °C/km) and high heat flow (of about 100 mW/m²) for the study area. The thermal regime for the area is inferred to be similar to that from the Salton trough.     

A spreadsheet method of estimating best-fit hydraulic gradients using head data from multiple wells

Hydraulic gradients from planar water tables, or piezometric surfaces, and horizontal flow regimes can be quickly and conveniently calculated from data sets involving numerous wells. The matrix-solving functions of a modem spreadsheet program (Excel) were used to determine the equation of a water-table plane, Ax + By + Cz - D = 0, and the equation coefficients were then used to determine the magnitude of the hydraulic gradient, according to gradient = square root of A2 + B2/C2, and its direction, according to alpha = arctan B/A, where alpha is the angle measured from the x-axis. A pre-prepared Excel file constructed to handle data from up to 20 wells at once is available for free downloading at www.geo.ku.edu/hydro/KUHydro.html.     

Whole genome data for omics-based research on the self-fertilizing fish Kryptolebias marmoratus

Genome resources have advantages for understanding diverse areas such as biological patterns and functioning of organisms. Omics platforms are useful approaches for the study of organs and organisms. These approaches can be powerful screening tools for whole genome, proteome, and metabolome profiling, and can be used to understand molecular changes in response to internal and external stimuli. This methodology has been applied successfully in freshwater model fish such as the zebrafish Danio rerio and the Japanese medaka Oryzias latipes in research areas such as basic physiology, developmental biology, genetics, and environmental biology. However, information is still scarce about model fish that inhabit brackish water or seawater. To develop the self-fertilizing killifish Kryptolebias marmoratus as a potential model species with unique characteristics and research merits, we obtained genomic information about K. marmoratus. We address ways to use these data for genome-based molecular mechanistic studies. We review the current state of genome information on K. marmoratus to initiate omics approaches. We evaluate the potential applications of integrated omics platforms for future studies in environmental science, developmental biology, and biomedical research. We conclude that information about the K. marmoratus genome will provide a better understanding of the molecular functions of genes, proteins, and metabolites that are involved in the biological functions of this species. Omics platforms, particularly combined technologies that make effective use of bioinformatics, will provide powerful tools for hypothesis-driven investigations and discovery-driven discussions on diverse aspects of this species and on fish and vertebrates in general.     

Methods of calculating petrophysical properties from lattice preferred orientation data

We consider the theoretical problems of calculating the physical properties of an aggregate from the constituent crystal properties and the lattice preferred orientation. The notion of a macroscopically homogeneous sample with an internally varing distribution of stress and strain fields is introduced to explain why further efforts have to be made to improve on the physically based Voigt and Reuss bounds. It is shown that the Voigt and Reuss bounds become increasingly separated with inceasing anisotropy, emphasising the need for better methods. The problem of highly anisotropic minerals is illustrated with polycrystals of plagioclase feldspar and biotite. Biotite is used to illustrate the mean velocity, the geometric mean and the self-consistent methods. The self-consistent method, which is generally accepted to give the best estimate, is almost identical to geometric mean recently introduced by Matthies and Humbert (1993) and similar to the arithmetic mean of the Voigt and Reuss bounds (Hill, 1952). The geometric mean has the powerful physical condition that the aggregate mean is equal to the mean of the inverse property (e.g. mean elastic stiffness and compliance). Despite its lack of theoretical justification the Hill average remains a useful estimate.     

基于磁力资料的南海及邻区磁性基底和居里面深度特征研究

磁性基底和居里面是研究地壳和岩石圈的地质构造和热演化过程的两个重要磁性界面.为了研究南海及邻区磁性基底和居里面所反映的深部构造及其热活动的地质效应,本文在对磁异常进行化极处理的基础上,采用最小曲率位场分离方法,获得了磁性基底和居里面引起的化极磁异常,利用双界面模型快速反演方法,反演了南海及邻区的磁性基底和居里面深度,研究了磁性基底、居里面深度及其分布特征,讨论了磁性基底、居里面与新生界深度之间相关性特征及其地质意义.研究表明,磁性基底深度5~20 km,洋盆南北两侧磁性基底走向分别以NE、NEE向为主,中南半岛周缘磁性基底呈NW、NNW走向.居里面深度15~32 km,宏观表现为"洋壳浅、周缘深"及周缘"北浅南深"的特征,洋盆地区居里面深度呈现"西南浅、东部深",洋壳与陆壳接触带在居里面深度上表现为梯级带特征.新生界深度与磁性基底深度相关性（Correlation between the depth of magnetic basement and Cenozoic,CDMBC）多以不规则形状分布,在盆地的沉积中心呈现正相关;新生界深度与居里面深度相关性（Correlation between the depth of Curie surface and Cenozoic,CDCSC）多呈NE、NEE向带状正相关分布,走向与盆地走向一致;莺歌海盆地、琼东南盆地、万安盆地南部和曾母盆地CDMBC呈正相关、CDCSC呈负相关,莺歌海相关性特征推测为：居里面随岩石圈变形隆起而抬升,磁性基底张裂下沉,发生大规模沉降引起;琼东南盆地相关性特征推测为：居里面随岩石圈变形下坳而下降,沉积中心与磁性基底下沉方向一致;万安盆地和曾母盆地相关性特征推测为：深部流体沿南海西缘断裂直接进入地壳,引起该处居里面深度变浅.     

弹性波被动源数据一次波估计方法研究

与声波被动源数据相比,弹性波被动源数据含有更加丰富的地下地质信息.弹性波被动源数据合成的虚拟炮记录中,同样存在表面相关多次波,影响数据的后续处理以及对地下地质构造的准确判断.对弹性波被动源数据合成的虚拟炮记录进行多次波的预测与匹配相减,将比声波数据更加困难.本文提出弹性波被动源数据一次波估计方法,针对各向同性介质条件下的弹性波被动源数据进行研究,依据声波被动源数据一次波估计方法的研究思路,推导弹性波被动源数据一次波估计方法的理论过程,充分考虑弹性波被动源数据所包含的波场成分,包括纵波、横波以及在传播过程中产生的转换波.在传播过程中,由于能量的损失,略去多次传播的高阶项,进而给出相应的目标函数.利用数据的稀疏性,结合L₁范数约束的最优化问题求解方法,避免对弹性波虚拟炮记录进行多次波的预测与匹配相减,解决多次波影响问题.本文通过合成的弹性波被动源数据,进行一次波估计的测试,并与主动源和合成的虚拟炮记录进行对比,验证本文提出方法的有效性与准确性.

     

Methods for estimating rotational components of seismic ground motion and their numerical comparisons*

Rotational components play an important role in natural earthquake research, engineering seismic investigation, building monitoring, seismic exploration and other fields. Traditional researches mainly focus on three translational components, but less on rotational ones. As the precision of rotational sensing techniques has increased, many scholars have paid more attention to the seismic rotational motions. Because the rotational observations are not very popular before and now, approximately converting the translational components into rotational components is utilized in rotation analysis. Based on numerical six-component seismic data with the finite difference method, we compare three different conversion methods, the travelling-wave, frequency-domain and the difference method, to analyze their characteristics and feasibilities when they are applied to estimate rotational components with translational observations.     

Spatial rainfall model using a pattern classifier for estimating missing daily rainfall data

Missing data in daily rainfall records are very common in water engineering practice. However, they must be replaced by proper estimates to be reliably used in hydrologic models. Presented herein is an effort to develop a new spatial daily rainfall model that is specifically intended to fill in gaps in a daily rainfall dataset. The proposed model is different from a convectional daily rainfall generation scheme in that it takes advantage of concurrent measurements at the nearby sites to increase the accuracy of estimation. The model is based on a two-step approach to handle the occurrence and the amount of daily rainfalls separately. This study tested four neural network classifiers for a rainfall occurrence processor, and two regression techniques for a rainfall amount processor. The test results revealed that a probabilistic neural network approach is preferred for determining the occurrence of daily rainfalls, and a stepwise regression with a log-transformation is recommended for estimating daily rainfall amounts.     

叠前地震资料提高分辨率处理方法及应用

对于隐蔽性油气藏和复杂构造油气藏,地震反射资料难以满足精度要求,需要开展地震资料提高分辨率处理.为此,使用反Q滤波、广义S变换、谱整形技术及谱白化技术等方法处理了涠西南地区叠前地震资料,并且分析了四种提高分辨率方法的应用效果.结果表明,四种方法都能提高主频、拓宽频带,其中反Q滤波的效果最好,广义S变换效果次之.反Q滤波结果中波组的低频和中高频能量分配合理,并且在中高频段反Q滤波结果的信噪比明显高于其它三种方法.此外,使用不同Q值对反Q滤波结果影响比较大.随着Q值增加,主频和带宽都稍有降低的趋势,但是信噪比提高了.当Q值大于200.0时信噪比改善的空间已经很小,这时再提高Q值,换来的只是主频降低和带宽缩小.     

Curie Point Depths of Albania Inferred from Ground Total Field Magnetic Data

Ground total magnetic field data of Albania were used to produce estimates of the Curie point isotherm. The strategy followed was to estimate the depth to the bottom of the deepest magnetic sources. Firstly, the average depth to the top of the deepest crustal block, z_t, was computed by linear fitting to the second lowest-frequency segment of the azimuthally averaged power spectrum of the total magnetic field data. Then, the depth to the centroid of the deepest crustal block, z₀, was computed by linear fitting to the lowest-frequency segment of the azimuthally averaged power spectrum of a distribution of magnetic dipoles. Finally, the depth to the bottom, the inferred Curie point depth, z_b, was calculated from z_b=2z₀−z_t. Curie depth estimates for Albania vary from about 17 to 25.5 km (below sea level). These results are consistent with the depths inferred by extrapolating geothermal gradient and heat-flow values, suggesting that the Curie point depth analysis is useful to estimate the regional thermal structure. It also suggests that the approach was valid and that ground total magnetic field data can be used for this purpose.     

Methods of analyzing radio polarization data

Abstract

Methods of estimating the strength and direction of galactic magnetic fields from radio polarization measurements are reviewed. Particular attention is paid to the analysis of the Faraday rotation in order to derive the large scale magnetic field structure. Ways in which an axisymmetric spiral field structure can be observationally distinguished from a bisymmetric spiral structure are described, as are the ways in which field symmetric with respect to the galactic plane can be distinguished from those that are antisymmetric.     

Evaluation of crystallization temperatures of mineral paragenesis in consolidated crust and uppermost mantle from explosion seismology data

In terms of the general endogeneous evolution of the lithosphere, the continental crystalline crust and the uppermost mantle, formed by regional metamorphic and magmatic processes, show mineral paragenesis stratification, expressed by a regular mineral sequence. The continuous macrolayering of mineral paragenesis through lithospheric depth profile is caused by phase transformations and variations in composition of complex natural systems, and affects the vertical distribution of seismic velocities,V _p,V _s, and other physical parameters.To evaluate palaeotemperatures (crystallization temperatures of mineral paragenesis), consistentV _pandV _s (Z) velocity models for the consolidated crust of two regionally separated areas of different geological structures — Precambrian shield (Voronezh Massif) and a young fold-mountain structure in the central part of the Transasian orogeneous belt (Himalaya) — are used as starting data.The velocity models are recalculated into (Z) and (Z) profiles (Z) being the seismic parameter. (Z) the Debye temperature). These, according to Debye theory, allow the determination of variations in entropy, thermodynamic and temperature gradients at the time of crustal generation.For the two regions chosen, palaeotemperature distributions are eventually calculated for the depth intervals given by velocity profiles. Crystallization temperatures calculated from seismic data show good agrrement with the values obtained from mineralogical thermobarometry.     

Correction factors for estimating suspended sediment export from loess catchments

This study focuses on problems related to estimating the exact exported suspended sediment quantity and quality from flow‐proportional samples taken at a fixed height at the outlet of a small, agricultural catchment in central Belgium. First, a correction factor is introduced to estimate the real suspended sediment concentration. This factor depends on the stream depth and the grain size distribution of the sediment and is necessary because of the presence of a concentration gradient in the stream. Secondly, water depth is corrected for the sediment volume, and thirdly, the grain size distribution of the exported sediment is corrected. The application of these correction factors leads to a reduction of the estimated total sediment export by 18 per cent (from 4810 Mg over a period of three years to 3930 Mg). Furthermore, the importance of small grain size fractions increases, leading to higher enrichment ratios for the clay fraction and lower ratios for the sand fraction. Suspended sediment export is thus more selective than was previously thought. Due to the overestimation of the sediment export, total phosphorus losses from the catchment are also overestimated by 17 per cent. This paper therefore emphasizes that point sampling of suspended sediment is biased and that a method should be incorporated to correct the sediment export in terms of quantity and quality. Copyright © 2001 John Wiley & Sons, Ltd.     

Comparison of methods for estimating earth resistivity from airborne electromagnetic measurements

Earth resistivity estimates from frequency domain airborne electromagnetic data can vary over more than two orders of magnitude depending on the half-space estimation method used. Lookup tables are fast methods for estimating half-space resistivities, and can be based on in-phase and quadrature measurements for a specified frequency, or on quadrature and sensor height. Inverse methods are slower, but allow sensor height to be incorporated more directly. Extreme topographic relief can affect estimates from each of the methods, particularly if the portion of the line over the topographic feature is not at a constant height above ground level. Quadrature–sensor height lookup table estimates are generally too low over narrow valleys. The other methods are also affected, but behave less predictably. Over very good conductors, quadrature–sensor height tables can yield resistivity estimates that are too high. In-phase–quadrature tables and inverse methods yield resistivity estimates that are too high when the earth has high magnetic susceptibility, whereas quadrature–sensor height methods are unaffected. However, it is possible to incorporate magnetic susceptibility into the in-phase–quadrature lookup table. In-phase–quadrature lookup tables can give different results according to whether the tables are ordered according to the in-phase component or the quadrature component. The rules for handling negative in-phase measurements are particularly critical. Although resistivity maps produced from the different methods tend to be similar, details can vary considerably, calling into question the ability to make detailed interpretations based on half-space models.