在地震预测研究中使用GPS和震源机制资料的问题探讨——与王秀文等商榷

针对"山西地区应力场变化与地震的关系"一文,本文指出:(1)P波初动存在误读的可能,出现矛盾符号是正常的;(2)地震仪器极性也可能反向,不校核可能会导致错误的震源机制解.还讨论如何正确地表示震源机制随时间的变化.分析认为GPS观测结果若不认真考虑观测误差、模型误差以及应变信息的层次,所得结果可能会与震源机制解反演的区域...     

通用测震数据获取软件包的设计与实现

在掌握数据源及相关数据格式的基础上,实现了读取SSS实时流数据、数据库连续波形、事件波形数据以及SEED文件的通用获取测震数据软件包.该软件包可使应用程序通过统一数据接口获得测震数据.经测试该软件包功能完整、高效易用,提高了地震科技人员的工作效率.     

Integration of space geodesy: A US National Geodetic Observatory

In the interest of improving the performance and efficiency of space geodesy a diverse group in the US, in collaboration with IGGOS, has begun to establish a unified National Geodetic Observatory (NGO). To launch this effort an international team will conduct a multi-year program of research into the technical issues of integrating SLR, VLBI, and GPS geodesy to produce a unified set of global geodetic products. The goal is to improve measurement accuracy by up to an order of magnitude while lowering the cost to current sponsors. A secondary goal is to expand and diversify international sponsorship of space geodesy. Principal benefits will be to open new vistas of research in geodynamics and surface change while freeing scarce NASA funds for scientific studies. NGO will proceed in partnership with, and under the auspices of, the International Association of Geodesy (IAG) as an element of the Integrated Global Geodetic Observation System project. The collaboration will be conducted within, and will make full use of, the IAG's existing international services: the IGS, IVS, ILRS, and IERS. Seed funding for organizational activities and technical analysis will come from NASA's Solid Earth and Natural Hazards Program. Additional funds to develop an integrated geodetic data system known as Inter-service Data Integration for Geodetic Operations (INDIGO), will come from a separate NASA program in Earth science information technology. INDIGO will offer ready access to the full variety of NASA's space geodetic data and will extend the GPS Seamless Archive (GSAC) philosophy to all space geodetic data types.     

A NATIONAL SYSTEM FOR EXCHANGE OF WATER DATA

ABSTRACT

This paper describes a system proposed by a Federal Interagency Work Group for providing access to large quantities of hydrological data. The National Water Data Exchange (NAWDEX) will be a confederation of data-acquisition organizations that are willing to make their data available to other data users. NAWDEX will take full advantage of on-going data handling activities and will provide the means to present data in a standard format. An ultimate goal of the system will be standardization and quality control of all of the data. Management of NAWDEX will be lodged in a Systems Central, which will assure that data are made available according to the time frame and the format required by the user. Access to the data will be through an index maintained by Systems Central. Full implementation of NAWDEX will require 4 to 5 years. A technical advisory group, consisting of members representing the Federal and non-Federal sectors, and topical work groups, established in response to specific needs, will provide input during the development of NAWDEX.     

21世纪重力与磁法勘探的展望

对21世纪重力与磁法勘探的仪器，数据处理技术，解释理论与方法，应用领域等方面的发展方向进行了分析与展望，发展航空标量，矢量，梯度重力测量和航空全梯度磁力测量，三分量磁力测量，提高综合信息采集能力；开展卫星重磁测量，综合卫星，航空，地面重磁测量资料研究地球结构与构造；发展高精度数据处理技术；重磁异常弱信号的提取，不同深度重磁异常的划分，低纬底变倾角化磁极以及位场面延拓；发展复杂条件下三维重磁场多参数综合反演可视化技术以及快速自动反演技术；探索磁性多参数的应用新领域，充分发挥磁法在环境污染调查中的作用并开拓应用新领域。     

各向异性扩散在电磁测深曲线平滑中的应用

电磁类地球物理方法由于极易受到各类噪声的干扰,使得估计的视电阻率曲线或相位曲线发生畸变,严重地影响了反演解译精度,如何对这类观测曲线进行合理有效地平滑,是目前数据处理中的重点,也是难点.本文把各向异性扩散（anisotropic diffusion AD）引入到曲线平滑中,提出了各向异性扩散的曲线平滑方法,以平滑点梯度值的降函数作为扩散速度,在梯度较大的位置予以较小的扩散速度以保护特征,在梯度较小的位置予以较大的扩散速度平滑噪声引起的扰动.同时为了减小“飞点”(outliers)对平滑的干扰,采用局部单调扩散进行预处理.实验结果显示,本文提出的平滑方法有效地平滑了噪声与“飞点”的干扰,恢复了曲线的基本形态,有效地保持了曲线基本特征.     

IPv6环境下地震观测网络系统构建

现有的地震观测网络系统都是基于IPv4协议开发的, 由于IPv4地址的匮乏, 无法对各个台站直接进行数据访问, 只能进行台网级数据访问; 另外, 现有的地震观测网络系统只是简单的通过用户名和密码进行验证, 而且传输时没有对用户名和密码进行加密, 安全性低。 通过研制基于IPv6协议的地震观测网络系统, 开发出了新型的IPv6环境下地震观测系统组网协议, 实现了台站级数据共享的目标; 使用各种类型实时数据加密方法, 保证了数据的完整性和机密性; 建立了新型的证书管理和认证机制; 在实验的基础上, 开发出了现有基于IPv4的地震观测系统向IPv6的过渡技术, 为下一步IPv6技术在地震行业的推广应用打下了良好的基础。     

基于“剪切+张裂”一般位错模型频率域求解微震震源机制

研究微地震的震源机制,获得压裂区域的破裂方向、尺度和应力状态等信息,在非常规油气开采过程中具有重要意义.对于微震,通常采用剪切位错或者矩张量模型对震源进行描述.本文从其实际发震机制出发,使用了"剪切+张裂"的一般位错点源模型,并基于此模型发展了一种利用全波形信息,通过波形振幅谱相关和初至约束,在频率域求解微震震源机制的方法.该方法适用于地面和井中观测,能够在得到常规震源参数（断层走向、倾角和滑动角）的同时给出裂缝断层剪切和张裂错动的距离信息,更直观体现破裂程度.理论数值测试证明方法有效、可行,在未滤波的情况下,实际数据的波形拟和结果仍较为一致,同时还发现错动距离与应力降等常规破裂参数并不严格相关,说明剪切、张裂错距可作为独立的新参数来定量评估水压致裂效果,指导工程开发进行.     

WATEQB – BASICA REVISION OF WATEQF FOR IBM PERSONAL COMPUTERS

Abstract. WATEQB is a BASICA revision for IBM PC of WATEQF, a geochemical model for chemical equilibrium and saturation of natural waters. This model computes the speciation of inorganic ions and complex species in solution for a given water analysis.
WATEQB is identical to the FORTRAN IV version from the standpoint of chemical process; however, the structure of the program differs from the original version. The structure of the program is not complex and allows the user to try some new modifications, if necessary.
Due to the limited memory of personal computers and in contrast to the large volume of the original version of the program, the author has created four data files, one of which is of random access type. Reading the thermochemical data is possible through these four data files.
This paper demonstrates how to set up the input data to execute WATEQB and describes the program rather than the model. The users of WATEQB will increase as personal computers become more common.     

Variations in the altitude of the F2 peak associated with trough-formation processes

A novel approach is described which can help to determine, from ground-based data, which of the possible production mechanisms for the mid-latitude F-region ionospheric trough is dominant during a particular event. This approach involves numerically modelling the possible causal mechanisms of the mid-latitude trough to see how each will affect the altitude of the F2-layer electron-concentration peak (hmF2), and then comparing these predictions with the observed variation of hmF2 during trough formation. The modelling work predicts that, if the neutral-wind velocity does not vary, hmF2 will remain almost constant if the trough is formed via stagnation, but will rise if it is formed as a result of high ion velocities or neutral upwelling. Observations made at Halley (76°S, 27°W, L=4.2), Antarctica, show that most frequently the only changes in hmF2 during trough formation are those expected due to variations in the neutral wind, which suggests that stagnation is the most common production mechanism. During the most geomagnetically active night studied, on which Ap varied between 18 and 32, there was a rise in hmF2 that cannot be explained by changes in the neutral wind. On this night the plasma also decayed faster, and the poleward edge of the trough was seen earlier than on other nights. These differences, together with the fact that the ion velocities remained relatively low, suggest the trough was caused by a change in neutral composition, possibly advected into the observing area.     

裸眼井中软地层多极子源激发的弹性波的传播

利用全波场和分波场算法对裸眼井软地层中由单极子源和多极子源（偶极子源和四极子源）激发的弹性波进行了数值研究,给出了声波全波列理论波形、纯纵波和纯横渡理论波形及纵横波的频率传输特性曲线.结果表明,纯纵波和纯横波是无频散几何衰减波;除了在v_s>（1-ρ₁/ρ₂）^1/2v₁的井眼中由单极子源激发的横波无传播特征频率外,纵横波都具有各自的传播特征频率,纵波的传播特征频率有一系列离散值,横波的传播特征频率只有一个且等于非对称模式波的截止频率;多极子源激发的纵波传播特征频率的最小值大于横波的传播特征频率,随着极子数增多,二者的特征频率值皆增高,差别也越大,因此,当声源发射主频较低,接近横波传播特征频率而远离纵波传播特征频率时,低频多极子源将压制纵波并激发出较强的横波,分析结果还表明,通常所说的低频多极子源激发出的横波波群实际上是纯横波与截止频率附近的非对称模式波叠加的结果,在其它条件不变时,二老对横波波群相对贡献的大小取决于声源激发主频的高低及其带宽的大小.     

Recent progress in understanding Pc1 pearl formation

We discuss recent progress in understanding mechanisms of formation of Pc1 pearl emissions. This problem remains unsettled in spite of many years of experimental and theoretical studies. Modern satellite observations by e.g. Polar and Cluster still do not reveal the whole picture experimentally since they do not stay long enough in the generation region to give a full account of all the spatio-temporal structure belonging to a single Pc1 event. Ground-based observations have also been extremely helpful, especially in combination with spacecraft data, but still existing experimental information allows one to advocate different scenario of Pc1 nature. On the other hand, a complete self-consistent theory taking into account all factors significant for Pc1 generation remains to be developed. Several mechanisms are discussed with respect to formation of Pc1 pearl spectrum, among them are nonlinear modification of the ionospheric reflection by precipitating energetic protons, modulation of ion-cyclotron instability by long-period (Pc3/4) pulsations, reflection of waves from layers of heavy-ion gyroresonances, and nonlinearities of wave generation process. We show that each of these mechanisms have their attractive features and explains certain part experimental data but any of them, if taken alone, meets some difficulties when compared to observations. We conclude that development of a refined nonlinear theory and further correlated analysis of modern satellite and ground-based data is needed to solve this very intriguing problem.     

Elastic full waveform inversion with probabilistic petrophysical clustering

Full waveform inversion aims to use all information provided by seismic data to deliver high-resolution models of subsurface parameters. However, multiparameter full waveform inversion suffers from an inherent trade-off between parameters and from ill-posedness due to the highly non-linear nature of full waveform inversion. Also, the models recovered using elastic full waveform inversion are subject to local minima if the initial models are far from the optimal solution. In addition, an objective function purely based on the misfit between recorded and modelled data may honour the seismic data, but disregard the geological context. Hence, the inverted models may be geologically inconsistent, and not represent feasible lithological units. We propose that all the aforementioned difficulties can be alleviated by explicitly incorporating petrophysical information into the inversion through a penalty function based on multiple probability density functions, where each probability density function represents a different lithology with distinct properties. We treat lithological units as clusters and use unsupervised K-means clustering to separate the petrophysical information into different units of distinct lithologies that are not easily distinguishable. Through several synthetic examples, we demonstrate that the proposed framework leads full waveform inversion to elastic models that are superior to models obtained either without incorporating petrophysical information, or with a probabilistic penalty function based on a single probability density function.     

DRAC: Dose Rate and Age Calculator for trapped charge dating

Accurate calculation of the environmental radiation dose rate (Ḋ) is an essential part of trapped charge dating methods, such as luminescence and electron spin resonance dating. Although the calculation of Ḋ is not mathematically complex, the incorporation of multiple variables and the propagation of uncertainties can be challenging. The Dose Rate and Age Calculator (DRAC) is an open access, web-based program which enables rapid Ḋ calculation for trapped charge dating applications. Users can select from recently published attenuation and conversion factors to make mathematically robust, reproducible Ḋ calculations. Comparison of DRAC calculated Ḋ values against the published Ḋ determinations of 422 samples from 32 studies results in a reproducibility ratio of 1.01 ± 0.05. It is anticipated that DRAC will facilitate easier inter-laboratory comparisons and will provide greater transparency for Ḋ calculations. DRAC will be updated to reflect the latest advances in Ḋ calculation and is freely accessible at www.aber.ac.uk/alrl/drac. The code for DRAC is available from github at https://github.com/DRAC-calculator/DRAC-calculator.     

The use of low frequencies in a full‐waveform inversion and impedance inversion land seismic case study

Velocity model building and impedance inversion generally suffer from a lack of intermediate wavenumber content in seismic data. Intermediate wavenumbers may be retrieved directly from seismic data sets if enough low frequencies are recorded. Over the past years, improvements in acquisition have allowed us to obtain seismic data with a broader frequency spectrum. To illustrate the benefits of broadband acquisition, notably the recording of low frequencies, we discuss the inversion of land seismic data acquired in Inner Mongolia, China. This data set contains frequencies from 1.5–80 Hz. We show that the velocity estimate based on an acoustic full‐waveform inversion approach is superior to one obtained from reflection traveltime inversion because after full‐waveform inversion the background velocity conforms to geology. We also illustrate the added value of low frequencies in an impedance estimate.     

The greenhouse effect

The greenhouse effect on the Earth is identified by the difference between the effective radiating temperature of the planet and its surface temperature. The difference between the energy emitted by the surface and that emitted upward to space by the upper atmosphere quantifies it; it can therefore be defined as the long wave energy trapped in the atmosphere. Climate forcing and the response of the climate system within which climate feedback mechanisms are contained, will be defined in this review. Quantitative examples will illustrate what could happen if the greenhouse effect is perturbed by the human activities, in particular if atmospheric CO₂ concentrations would double in the future. Recent measurements by satellites of the greenhouse effect will be given. The net cooling effect of clouds on the Earth and whether or not there will be less cooling by clouds as the planet warms, are discussed following a series of papers recently published by Ramanathan and his collaborators.     

Adjusting chlorophyll-a estimates through temporal weighting based on the seasonal development of phytobiomass

North-temperate lakes have been shown to progress through a general cycle of chlorophyll-a phenology. Because estimates of seasonal means are often based on only a few samples collected against this variable background, these estimates can be biased or uncertain. Our goal was to reduce the seasonal uncertainty and thereby produce more accurate estimates of chlorophyll concentration by defining a correction for phenological development. Time-series data from 149 lake-years were used to develop equations from which chlorophyll values could be corrected to the seasonal mean in relation to their particular date of measurement. However, we found the seasonal correction to be ineffectual in reducing uncertainty about nutrient-response regressions. After reviewing a number of hypotheses, we conclude that the correction derived from the average response for many lakes will be inadequate to adjust for the seasonal pattern occurring within any particular lake. This occurs because the temporal weighting correction, generated through repeated averaging, underestimates the seasonal variability which exists among individual lakes. An effective correction, if it is developed will have to be based on patterns within single lakes or possibly different lakes within a single region.     

Simulating the temporal and spatial dynamics of the North Sea using the new model GETM (general estuarine transport model)

Here we present results of a 1-year realistic North Sea simulation from the new model GETM (general estuarine transport model) and assess the capabilities of this model by comparing them to model results from the well-known HAMSOM (Hamburg shelf sea and ocean model) model, in situ data from the North Sea project and satellite-derived sea-surface temperature data. The annual cycle and the spatial variability of stratification and mixing in the North Sea is simulated. It is shown that the new model is successful in reproducing the general temporal and spatial dynamics of the North Sea. The major advantages of GETM for achieving improved results in this simulation are the implementation of general vertical coordinates, of a state-of-the-art turbulence model and of higher-order advection schemes. By exploiting the full capabilities of these features a more realistic simulation could be achieved. We found that the greatest differences in the model results are produced by applying advection schemes of different complexity. Here we are able to demonstrate that better advection schemes lead to stronger horizontal gradients and stronger vertical stratification during summer. When comparing these results to measurements from the North Sea project and to satellite data, we find that these stronger gradients are more realistic. Therefore, we consider it as essential to use such high-order advection schemes if the spatial variability of estuarine or shelf seas like the North Sea is to be resolved adequately. The advanced turbulence closure scheme also contributed to more realistic simulation of the vertical stratification. Finally, general vertical coordinates better resolve the shallow regions, but are also useful for the deeper regions, as they allow a better estimation of sea-surface temperature compared to traditional coordinates.Responsible Editor: Phil Dyke     

使用MATLAB获取地震前兆数据库中观测数据

地震前兆数据库中的观测数据具有极高的科研及应用价值.MATLAB作为具有强大功能的科学计算及分析平台是地震前兆数据分析的得力工具.虽然如此,但MATLAB提供的数据库工具箱并不能完全操控基于ORACLE的前兆数据库.因此,直接通过MATLAB获取前兆观测数据进行后续的前兆统计及分析工作受到了影响.本文首先介绍了主流的数据库访问机制以及在MATLAB中如何实现以这些机制对前兆数据库进行访问并指出了其中的缺陷;然后,给出了通过组件对象模型(COM)的访问方式获取前兆观测数据的方法;最后,给出了使用该方法进行函数封装后,在实际应用中的实例的全部源代码.通过本文提供的方法,台站工作人员或科研人员可以快速、方便的通过MATLAB获取前兆数据库中的相关信息进行科研及日常管理工作.