Use of variograms for field magnetometry analysis in Upper Silesia Industrial Region

Measurements of topsoil magnetic susceptibility are often used for quick assessment of soil contamination of anthropogenic origin, with heavy metals or other pollutants. However, because of complicated correlations between low-field magnetic susceptibility (shortened to magnetic susceptibility) of topsoil and soil pollution, the outcome of a field magnetometry survey can not be related directly to soil pollution. For each case study, the results should be interpreted on their own taking into account not only the type of pollution but also pedogenic, biogenic and environmental factors. In practice, it is very difficult to measure and consider all these factors. Here we illustrate the merit of geostatistical methods, which are focused on the spatial variability of a phenomenon, in the interpretation of soil magnetometry results. This article presents the analysis of spatial variability of top soil layers magnetic susceptibility-within the Upper Silesia Industrial Region (USIR)-using semivariance analysis. It also explains how to adjust the sampling density of field magnetometry measurements to spatial variability of the soil pollution as well as to the spatial scale of the investigated area. For this purpose, the values of magnetic susceptibility have been measured by using various sampling densities at areas of different size located within USIR. This enabled to determine the main scales of magnetic susceptibility spatial variability of soils within USIR using semivariance. A few distinct scales of variability were found from the site scale to a more regional scale. Variability ranges of 30 km, 12 km, and 5 km refer to the large regional scale, whereas smaller ranges of few hundreds down to a few tens of meters, can be attributed to the local (site) scale. In addition, the precision of the measuring campaigns, performed within USIR with different sampling densities, was compared through the analysis of the spatial variability of the soil magnetic susceptibility signal by using ordinary kriging. jarek97@yahoo.com, piotr.fabijanczyk@is.pw.edu.pl     

Two-dimensional power series solution for non-axisymmetrical transport in a radially convergent tracer test with scale-dependent dispersion

It has been known for many years that dispersivity increases with solute travel distance in a subsurface environment. The increase of dispersivity with solute travel distance results from the significant variation of hydraulic properties of heterogeneous media and was identified in the literature as scale-dependent dispersion. This study presents an analytical solution for describing two-dimensional non-axisymmetrical solute transport in a radially convergent flow tracer test with scale-dependent dispersion. The power series technique coupling with the Laplace and finite Fourier cosine transform has been applied to yield the analytical solution to the two-dimensional, scale-dependent advection–dispersion equation in cylindrical coordinates with variable-dependent coefficients. Comparison between the breakthrough curves of the power series solution and the numerical solutions shows excellent agreement at different observation points and for various ranges of scale-related transport parameters of interest. The developed power series solution facilitates fast prediction of the breakthrough curves at any observation point.     

Modeling doline populations with logistic growth functions

Initiation and evolution of a large karst depression population consisting of parent and daughter dolines were spatially and temporally modeled using logistic growth functions. Logistic growth models are well suited for analyzing doline population initiation and evolution because they reflect the density‐dependent growth mechanisms present in the evolution of karst depressions. Seven assumptions based upon previous studies were refined into mathematical statements and tested using more than 2000 dolines from a subpopulation of both parent and daughter karst depressions on the Western Highland Rim and Pennyroyal Plain of Tennessee and Kentucky. Logistic growth models quantify the initiation and evolution of doline populations and interface well with recent models describing the evolution of three‐dimensional conduit systems. Logistic growth models should apply to modeling other doline populations in karst terranes with hydraulically efficient, three‐dimensional conduit systems. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical experiments on the generation of long ocean waves in coastal waters of the Buenos Aires province,Argentina

A numerical model (two horizontal dimensions, vertically integrated) is used to investigate the generation of long ocean waves, ranging from 20 min to almost 2 h, at Buenos Aires continental shelf. The domain includes the Río de la Plata estuary and the continental shelf together and extends from 33.5° to 40.5°S latitude, and from 51° to 63°W longitude. Sea-level oscillations are modeled by forcing with passage of atmospheric cold fronts and atmospheric gravity waves. Both forcing mechanisms, which have been present during high activity lapses of long ocean waves, are mathematically implemented. After several numerical simulations, it is concluded that the pressure and wind fields associated to cold fronts do not generate long ocean waves in the area, though they do produce disturbances with periods longer than the tidal ones. On the other hand, it is so concluded that atmospheric gravity waves are an effective mechanism to force long ocean waves. Results obtained show that generation of long ocean waves is highly sensitive depending on the propagation direction and the phase speed of the atmospheric gravity waves. The long ocean wave event detected during the large-amplitude gravity-wave event of 13 October 1985 is successfully simulated. Finally, all our results suggest that atmospheric gravity waves are a highly effective mechanism forcing for the generation of long ocean waves in Buenos Aires coastal waters.     

Modeling of Stress Changes at Mayon Volcano,Philippines

In this paper, we analyze the stress change associated with the inverted volcanic source models at Mayon volcano, Philippines, where there are gravity changes without significant deformation after the 1993 eruption. We detail the applicable data and the associated inversion techniques and models prior to calculating the appropriate stress changes. It is determined that the stress change associated with the central magmatic source produces compressional stress changes at a secondary source to the northwest, prompting a change in the local water storage in the underlying fractured rock.     

对重力势二阶垂直偏导数的一种新的取样方法

从目前的研究来看，二阶垂直导数（X和Y两个方向）的求解方法较简单，它采用的五个因子的数字滤波的最小二乘法来优化数据。但是在较小区域面积时（对工程和环境应用中的重力数据），重力异常很容易放大和被区域重力场掩盖，我们认为产生的缘由是表示重力场的x和y方向上的二次多项式引起的。因此，为了逼近重力水平方向的偏导可以设计一个好的滤波器来拟合在xy平面的重力场值，再得出想要的重力垂直方向偏导，文中最后对提出的方法在前人经典模型上和野外实践进行了验证。     

几种Q矢量的比较

利用 1 991 - 0 7- 0 5T2 0— 0 6T2 0一次江淮梅雨锋暴雨过程实况资料 ,结合地面降水分布 ,从定性 (矢量场 ,即 Q)、定量 (散度场 ,即 2 · Q)的角度细致、具体地比较分析了准地转 Q矢量、半地转 Q矢量、非地转 Q矢量及湿 Q矢量的诊断特性。结果表明 :半地转 Q矢量诊断能力优越于准地转 Q矢量。非地转 Q矢量、湿 Q矢量诊断能力明显优越于准地转 Q矢量、半地转 Q矢量。准地转 Q矢量诊断能力最差 ;从定性的角度分析 ,非地转 Q矢量与湿 Q矢量诊断能力相差不大。但进一步定量分析发现 ,在整个梅雨锋暴雨过程中 ,湿 Q矢量的诊断能力大于其他 Q矢量 ;70 0 h Pa高度上各 Q矢量的矢量、散度辐合区较其在 850 h Pa和 50 0 h Pa上对降水反映更好 ,尤其是70 0 h Pa湿 Q矢量散度辐合区与降水区有非常好的对应关系。     

湛江港邻近海域台风浪的模拟研究

对第三代近岸海浪数值模式SWAN及包含的物理过程进行了简要介绍，利用该模式对影响南海湛江港海域的二次台风浪过程进行了模拟研究：由藤田台风风场模型同化相应时刻的台风要素、NCAR/NCEP网格点资料、单站观测资料后，提供模式所需风场；利用自嵌套的方式，提供模式波谱边界条件；两次模拟结果与实际海浪观测资料相符较好，可以为该海域台风浪的模拟预报提供较为重要的参考。     

多部多普勒雷达同步探测三维风场反演系统

“我国重大气候和天气灾害的形成机理与预测理论研究”项目首次启用双多普勒雷达组网同步观测，如何从同步观测的暴雨资料中反演三维场，具有重要的理论和现实意义，根据项目需要，建立了国内首套多部多普勒天气雷达同步探测三维风场反演软件系统，并反演出2001年7月13日由双多普勒雷达同步观测到的一个中β尺度暴雨系统的三维风场，结果显示，该三维风场结构合理，系统能满足外场试验的要求。     

弹道扰动引力的有限元逼近算法

为了克服多项式逼近弹道扰动引力的缺点，根据有限元插值的原理，采用了对弹道周围空间区域进行有限元剖分的方法，利用剖分单元各顶点的扰动引力分量内插出弹道点对应的扰动引力分量值。结果表明，文中提出的逼近算法能够快速精确可靠地逼近弹道扰动引力，是一种具有应用价值的方法。