Interpolation by regularized spline with tension: I. Theory and implementation

Bivariate and trivariate functions for interpolation from scattered data are derived. They are constructed by explicit minimization of a general smoothness functional, and they include a tension parameter that controls the character of the interpolation function (e.g., for bivariate case the surface can be tuned from a membrane to a thin steel plate), Tension can be applied also in a chosen direction, for modeling of phenomena with a simple type of anisotropy. The functions have regular derivatives of all orders everywhere. This makes them suitable for analysis of surface geometry and for direct application in models where derivatives are necessary. For processing of large datasets (thousands of data points), which are now common in geosciences, a segmentation algorithm with a flexible size of overlapping neighborhood is presented. Simple examples demonstrating flexibility and accuracy of the functions are presented.On leave from the Department of Physical Geography and Cartography, Comenius University, Mlynská dolina, Bratislava, Czechoslovakia.On leave from the Institute of Physics, Dúbravská cesta 9, Bratislava, Czechoslovakia.     

变参考慢度Born近似傅氏偏移

针对常规Born近似傅氏偏移方法对于剧烈横向变速介质不能精确成像的状况而提出了变参考慢度Born近似傅氏偏移，理论上解决了任意速度变化地质模型的偏移成像问题。此外，为进一步提高复杂地层的成像精度和波场延拓算子的稳定性，对散射波场的计算公式作了改进。将改进的方法运用于盐丘模型的正演和偏移试验，并与常规Born近似偏移方法相比较，可明显看出变参考慢度Born近似傅氏偏移方法在效果上要优于后者，其处理速度横向变化的能力大大增强。     

济宁许厂煤矿伴生分散元素铼的富集成矿探析

铼是一种分散金属,目前尚没有在煤层中富集成矿的案例。该文依据山东省济宁煤田许厂煤矿16上、17煤层光谱分析成果,结合铼元素的性质、物质来源和石炭-二叠纪太原组的聚煤环境等因素进行综合分析,论证了其富集可能性。简要阐述了铼富集成矿的意义,建议煤田地质勘查工作中加强对煤中伴生铼元素的分析和研究。     

散乱数据插值的HASM方法

本文运用高精度曲面建模(HASM)方法,研究了空间散乱数据插值算法,并以陕西咸阳彬县大佛寺煤矿区的实测数据为案例,使用交叉统计检验方法比较分析HASM方法与地理信息系统(GIS)领域中常用传统插值方法的空间插值结果。结果表明,HASM方法具有较高的精度。本文还分析了传统常用插值方法的一些统计结果,得出一些结论,这些结论对于使用常用传统GIS空间插值方法的研究者有很好的参考作用。     

迭代法井间地震衍射层析成像

从各向同性非均匀介质中弹性波波动方程出发，导出了在弱不均匀条件下，频率－波数域衍射场方程，给出了格林函数解析解。在电磁波地球物理层析成像数据合成的基础上，用类“体电流法（ＶＣＭ）”，将解析解线性化，得到未作Ｂｏｒｎ近似或Ｒｙｔｏｖ近似的衍射场数值解，便于计算合成数据。采用迭代法进行反演，每步迭代中，采用等比Ｎｅｕｍａｎｎ级数解；反问题用ＳＶＤ法求解，探讨了阻尼因子的取法。数值模型试验表明，这种算法是有效和成功的，成像的分辨率较高，对介质的均匀性要求有所降低，为实际应用创造了进一步条件。     

Effect of relative humidity on the angular distribution of scattered light intensity for rural and urban aerosols

The change in atmospheric relative humidity affects the physical and optical properties of aerosol particles. It would be interesting to study the effect of an increase in relative humidity on the angular scattering of light by aerosols (by incorporating the changes due to it as the complex refractive index and the parameters of the size distribution function). In the present paper we have computed the angular scattering function for rural and urban aerosols, for light of wavelength 0.55 μn. The results obtained for these two models, representative of different environments, are interesting and show some discriminating features.     

The effect of single vegetation elements on wind speed and sediment transport in the Sahelian zone of Burkina Faso

Soil loss caused by wind erosion is a widespread phenomenon in the Sahelian zone of West Africa. According to Sahelian farmers, scattered vegetation standing in amongst the crop has the potential for a wind erosion control strategy. This study was conducted to study the effect of single vegetation elements on the pattern of average wind speed and sediment transport. This was done by two experiments that were carried out during the rainy seasons of 2002 and 2003 in north Burkina Faso, West Africa. Wind speeds were measured using three sonic anemometers, at a sampling frequency of 16 Hz. Sediment transport was determined by calculating the mass fluxes from 17 MWAC catchers. In this study, a shrub was defined as a vegetation element with branches until ground and a tree as a vegetation element with a distinctive trunk below a canopy. Behind shrubs wind speed near the soil surface was reduced up to approximately seven times the height of the shrub. The observed reduction in wind speed in the area where wind speed was reduced was 15 per cent on average. At the sides of the shrub, wind speed was increased, by on average 6 per cent. As the area of increase in wind speed is one‐third of the area of decrease in wind speed, the net effect of a shrub is a reduction in wind speed. A similar pattern was visible for the pattern of sediment transport around a shrub. Downwind of a shrub, sediment transport was diminished up to seven times the height of the shrub. Probably most of this material was trapped by the shrub. Trees showed a local increase of wind around the trunk, which is expected to relate to an increase in sediment transport around the trunk. Mass flux measurements of sediment transport were not made, but visual observations in the field substantiate this. Behind the canopy of a tree, a tree acts similarly to a shrub regarding its effects on average wind speed, but as a tree is generally a larger obstacle than a shrub the extent of this effect is larger than for shrubs. Thus, whereas shrubs are more effective than trees regarding their direct effect on soil loss by trapping sand particles near the soil surface, trees are more effective in affecting soil loss indirectly by reducing the wind speed downwind more effectively than shrubs. Therefore, to reduce soil loss in an area, the presence of both trees and shrubs is crucial. Copyright © 2007 John Wiley & Sons, Ltd.     

Sources of Inaccuracy in Boron Isotope Measurement Using LA-MC-ICP-MS

Laser ablation multi-collector-inductively coupled plasma-mass spectrometry (LA-MC-ICP-MS) has become a valuable tool for the in situ measurement of the boron isotope composition of geological samples at high (tens to hundreds of μm) spatial resolution. That said, this application suffers from significant analytical challenges. We focus in this study on the underlying processes of two of the main causes for inaccuracies using this technique. We provide empirical evidence that not only Ca ions (Sadekov et al. 2019, Standish et al. 2019, Evans et al. 2021) but also Ar ions, that are reflected within the flight tube of the mass spectrometer, are the source for previously reported issues with spectral baselines. We also address the impact of plasma conditions on the instrumental mass fractionation as a source for matrix- and mass-load-related analytical biases. Comparing experimental data with the results of a dedicated release and diffusion model (RDM) we estimate that a close to complete (~ 97%) release of boron from the sample aerosol is needed to allow for consistently accurate LA boron isotope measurement results without the need for corrections.     

On The Origin of The High-Perihelion Scattered Disk: The Role of The Kozai Mechanism And Mean Motion Resonances

We study the transfer process from the scattered disk (SD) to the high-perihelion scattered disk (HPSD) (defined as the population with perihelion distances q > 40 AU and semimajor axes a>50 AU) by means of two different models. One model (Model 1) assumes that SD objects (SDOs) were formed closer to the Sun and driven outwards by resonant coupling with the accreting Neptune during the stage of outward migration (Gomes 2003b, Earth, Moon, Planets 92, 29–42.). The other model (Model 2) considers the observed population of SDOs plus clones that try to compensate for observational discovery bias (Fernández et al. 2004, Icarus , in press). We find that the Kozai mechanism (coupling between the argument of perihelion, eccentricity, and inclination), associated with a mean motion resonance (MMR), is the main responsible for raising both the perihelion distance and the inclination of SDOs. The highest perihelion distance for a body of our samples was found to be q = 69.2 AU. This shows that bodies can be temporarily detached from the planetary region by dynamical interactions with the planets. This phenomenon is temporary since the same coupling of Kozai with a MMR will at some point bring the bodies back to states of lower-q values. However, the dynamical time scale in high-q states may be very long, up to several Gyr. For Model 1, about 10% of the bodies driven away by Neptune get trapped into the HPSD when the resonant coupling Kozai-MMR is disrupted by Neptune’s migration. Therefore, Model 1 also supplies a fossil HPSD, whose bodies remain in non-resonant orbits and thus stable for the age of the solar system, in addition to the HPSD formed by temporary captures of SDOs after the giant planets reached their current orbits. We find that about 12 – 15% of the surviving bodies of our samples are incorporated into the HPSD after about 4 – 5 Gyr, and that a large fraction of the captures occur for up to the 1:8 MMR (a ⋍ 120 AU), although we record captures up to the 1:24 MMR (a ≃ 260 AU). Because of the Kozai mechanism, HPSD objects have on average inclinations about 25°–50°, which are higher than those of the classical Edgeworth–Kuiper (EK) belt or the SD. Our results suggest that Sedna belongs to a dynamically distinct population from the HPSD, possibly being a member of the inner core of the Oort cloud. As regards to 2000 CR₁₀₅ , it is marginally within the region occupied by HPSD objects in the parametric planes (q,a) and (a,i), so it is not ruled out that it might be a member of the HPSD, though it might as well belong to the inner core.