利用探地雷达频谱反演层状介质几何与电性参数

通过对地下层状介质探地雷达(GPR)回波广义反射系数的奇偶分解,建立了联系GPR反射系数序列频谱与介质几何参数、电性参数的代价函数,这些介质参数包括地下反射面的深度、层厚度以及各层的介电常数和电导率,从而提出了一种由GPR频谱同时估算地下介质多参数的全局优化反演方法.为了对多参数全局优化算法给出一个合理的参数初值,研究了不同参数对反射系数序列频谱属性的影响规律,提出了利用不同频谱属性分别估算不同参数的分步反演方法.以分步反演方法得到的结果作为多参数全局优化反演的初值,可以极大地提高反演计算的效率和反演结果的可靠性.用理论模型合成数据和GPR公路检测数据对本文方法进行了测试,结果表明本文方法效果良好,具有较高的分辨率,能较好的给出厚度小于调谐厚度的薄层的深度、厚度和介电常数等参数.     

混合ADI-FDTD亚网格技术在探地雷达频散媒质中的高效正演

提出混合ADI-FDTD亚网格技术开展频散介质GPR正演,即在物性参数变化剧烈局部区域采用细网格剖分ADI-FDTD计算,其他的区域采用粗网格剖分常规FDTD计算,ADI-FDTD突破了CFL条件的限制,可选取与粗网格一致的大时间步长,有效地提高了计算效率.本文首先基于Debye方程,推导了粗网格FDTD及细网格ADI-FDTD频散介质差分格式,着重对粗细两种网格结合的场值交换方式进行了深入探讨,给出了该算法的计算流程.然后以一个薄层模型为例,分别应用粗网格、细网格、混合ADI-FDTD亚网格算法对该模型进行正演,计算资源的占用及模拟精度说明了混合ADI-FDTD亚网格算法的优势.最后,建立频散介质与非频散介质的组合模型,应用3种方法对该模型进行正演,对比3种方法优劣,分析雷达剖面中非频散介质及频散介质中波形特征,有效地指导雷达资料的精确解释.     

消除探地雷达数据的子波衰减和频散的反滤波方法

消除探地雷达数据的子波衰减和频散可以很好地提高探地雷达的勘探深度和勘探分辨率.常用的消除探地雷达数据的子波衰减和频散方法为反Q滤波方法.该方法需要利用地下介质的Q参数,但是正确求取地下介质的Q参数很困难.针对这一问题,本文提出了一种消除探地雷达数据的子波衰减和频散的反滤波方法.该方法以地下介质反射系数是随机数为前提,利用地下介质等效滤波器具有最小相位这个特性,通过求取等效滤波器的振幅谱来求取等效滤波器的反滤波器.最后,利用该反滤波器对探地雷达数据进行反滤波,实现消除探地雷达数据的子波衰减和频散.     

基于PML边界条件的二阶电磁波动方程GPR时域有限元模拟

完全匹配层(PML)作为一种稳定高效的吸收边界条件,广泛应用于基于一阶电磁波动方程的探地雷达(GPR)数值模拟中.为解决基于二阶电磁波动方程的GPR数值模拟的吸收边界问题,本文借鉴二阶弹性波动方程的PML边界条件构建思想,提出了一种适合二阶电磁波动方程GPR时域有限元模拟的PML边界条件.从二阶电磁波动方程出发,基于复拉伸坐标变换,推导了PML算法的频域表达式;通过合理构造辅助微分方程,得到了PML算法的时域表达式,并以变分形式(弱形式)加载到GPR时域有限元方程中,实现了PML边界条件在二阶电磁波动方程GPR时域有限元模拟中的应用.在此基础上,对比了无边界条件、Sarma边界条件和PML边界条件下均匀模型的波场快照、单道波形、时域反射误差和能量衰减曲线,结果表明:PML边界条件的吸收效果要远优于Sarma边界条件,具有近似零反射系数.一个复杂介质模型的正演模拟验证了PML边界条件在非均匀地电结构中电磁波传播模拟的良好吸收效果.     

A strategy for an accurate estimation of the basal permittivity in the Martian North Polar Layered Deposits

This paper deals with the investigation of the Mars subsurface by means of data collected by the Mars Advanced Radar for Subsurface and Ionosphere Sounding working at few megahertz frequencies. A data processing strategy, which combines a simple inversion model and an accurate procedure for data selection is presented. This strategy permits to mitigate the theoretical and practical difficulties of the inverse problem arising because of the inaccurate knowledge of the parameters regarding both the scenario under investigation and the radiated electromagnetic field impinging on the Mars surface. The results presented in this paper show that it is possible to reliably retrieve the electromagnetic properties of deeper structures if such strategy is accurately applied. An example is given here, where the analysis of the data collected on Gemina Lingula, a region of the North Polar layer deposits, allowed us to retrieve permittivity values for the basal unit in agreement with those usually associated to the Earth basaltic rocks.     

Inverse determination of soil density and stress state using dispersion wave measurements and cone penetration tests in a non-layered soil

The determination of the void ratio and stress state distribution in undisturbed soil is still an ambitious aim which cannot be reached by sampling from bore holes. Therefore, an alternative method is proposed to determine soil density and stress state using dispersion wave measurements and cone penetration tests. Dispersion wave measurements and resonant-column tests are carried out to measure the shear wave velocity and shear modulus distribution with depth. Using finite element calculations a relationship between cone penetration resistance, stress state and void ratio is derived. From these results the void ratio distribution and the stress state can be calculated inversely. The inverse method is applied to real test data. The results of the inverse parameter determination are shown and assessed regarding the possibilities and the limitations of the presented method.     

Ground penetrating radar inversion in 1-D: an approach for the estimation of electrical conductivity, dielectric permittivity and magnetic permeability

This paper presents a method for inverting ground penetrating radargrams in terms of one-dimensional profiles. We resort to a special type of linearization of the damped E-field wave equation to solve the inverse problem. The numerical algorithm for the inversion is iterative and requires the solution of several forward problems, which we evaluate using the matrix propagation approach. Analytical expressions for the derivatives with respect to physical properties are obtained using the self-adjoint Green's function method. We consider three physical properties of materials; namely dielectrical permittivity, magnetic permeability and electrical conductivity. The inverse problem is solved minimizing the quadratic norm of the residuals using quadratic programming optimization. In the iterative process to speed up convergence we use the Levenberg–Mardquardt method. The special type of linearization is based on an integral equation that involves derivatives of the electric field with respect to magnetic permeability, electrical conductivity and dielectric permittivity; this equation is the result of analyzing the implication of the scaling properties of the electromagnetic field. The ground is modeled using thin horizontal layers to approximate general variations of the physical properties. We show that standard synthetic radargrams due to dielectric permittivity contrasts can be matched using electrical conductivity or magnetic permeability variations. The results indicate that it is impossible to differentiate one property from the other using GPR data.     

A dripper-TDR method for in situ determination of hydraulic conductivity and chemical transport properties of surface soils

Field determined hydraulic and chemical transport properties can be useful for the protection of groundwater resources from land-applied chemicals. Most field methods to determine flow and transport parameters are either time or energy consuming and/or they provide a single measurement for a given time period. In this study, we present a dripper-TDR field method that allows measurement of hydraulic conductivity and chemical transport parameters at multiple field locations within a short time period. Specifically, the dripper-TDR determines saturated hydraulic conductivity (K_s), macroscopic capillary length (λ_c), immobile water fraction (θ_im/θ), mass exchange coefficient (α) and dispersion coefficient (D_m). Multiple dripper lines were positioned over five crop rows in a field. Background and step solutions were applied through drippers to determine surface hydraulic conductivity parameters at 44 locations and surface transport properties at 38 locations. The hydraulic conductivity parameters (K_s, λ_c) were determined by application of three discharge rates from the drippers and measurements of the resultant steady-state flux densities at the soil surface beneath each dripper. Time domain reflectometry (TDR) was used to measure the bulk electrical conductivity of the soil during steady infiltration of a salt solution. Breakthrough curves (BTCs) for all sites were determined from the TDR measurements. The K_s and λ_c values were found to be lognormally distributed with average values of 31.4 cm h⁻¹ and 6.0 cm, respectively. BTC analysis produced chemical properties, θ_im/θ, α, and D_m with average values of 0.23, 0.0036 h⁻¹, and 1220 cm² h⁻¹, respectively. The estimated values of the flow and transport parameters were found to be within the ranges of values reported by previous studies conducted at nearby field locations. The dripper TDR method is a rapid and useful technique for in situ measurements of hydraulic conductivity and solute transport properties. The measurements reported in this study give clear evidence to the occurrence of non-equilibrium water and chemical movement in surface soil. The method allows for quantification of non-equilibrium model parameters and preferential flow. Quantifying the parameters is a necessary step toward determining the influences of surface properties on infiltration, runoff, and vadose zone transport.     

An evaluation of the universal soil loss equation and field techniques for assessing soil erosion on a tropical alfisol in western Nigeria

Four techniques for soil erosion assessment were compared over two consecutive seasons for bare-fallow plots and a maize-cowpea sequence in 1985 at IITA, Ibadan, Nigeria. The techniques used were: tracer (aluminium paint), nails (16 and 25), the rill method, and the Universal Soil Loss Equation (USLE). Soil loss estimated by these techniques was compared with that determined using the runoff plot technique. There was significantly more soil loss (P < 0·01) in bare-fallow than in plots under maize (Zea mays) or cowpea (Vigna unguiculata). In the first season, soil loss from plots sown to maize was 40·2 Mg ha^?1 compared with 153·3 Mg ha^?1 from bare-fallow plots. In the second season, bare-fallow plots lost 87·5 Mg ha^?1 against 39·4 Mg ha^?1 lost from plots growing cowpea. The techniques used for assessing erosion had no influence on the magnitude of soil erosion and did not interfere with the processes of erosion. There was no significant difference (P < 0·05) between soil erosion determined by the nails and the runoff plot technique. Soil loss determined on six plots (three under maize, three bare-fallow) by the rill technique, at the end of the season, was significantly lower (P < 0·05) than that determined by the runoff plot technique. The soil loss estimated by the rill method was 143·2, 108·8 and 121·9 Mg ha^?1 for 11, 11, and 8 per cent slopes respectively, in comparison with 201·5, 162·0, and 166·4 Mg ha^?1 measured by the runoff plot method. Soil loss measured on three bare-fallow plots on 10 different dates by the rill technique was also significantly lower (P < 0·01) than that measured by the runoff plot. In the first season the USLE significantly underestimated soil loss. On 11, 11, and 8 per cent slopes, respectively, soil loss determined by the USLE was 77, 92, and 63 per cent of that measured by the runoff plot. However, in the second season there was no significant difference between soil loss determined by the USLE and that determined by the conventional runoff plot technique.     

The effects of fire on the surface properties and seed germination in two shallow monoliths from a rangeland soil subjected to simulated raindrop impact and water erosion

This paper describes the methods used and some preliminary results of simulated erosion studies on soils with cryptogamic crusts from a semiarid rangeland environment. Two 0·3 m² shallow monoliths were collected from the upper 20 cm of a Typic Haplargid from the semiarid Australian rangelands and subjected to a range of rainfall intensities and durations representing potentially erosive summer and winter rainfall events. One of the monoliths was cleared of vegetation by a simulated low intensity bushfire. Macro- and micromorphological properties of the surface, as well as runoff and erosion losses, were measured during the experiment. Runoff and erosion losses were, as expected, greater for all conditions on the burned than on the unburned monolith. Intensive rainfall damaged the cryptogamic crust unprotected by vegetation by widening and deepening desiccation cracks around the cryptogams, and breaking away and dispersing larger soil fragments from the crack margins. The burned and eroded surfaces provided a much poorer environment for seed entrapment, germination, and growth than did the unburned surface.     

Variations in the depth distribution of phosphorus in soil profiles and implications for model-based catchment-scale predictions of phosphorus delivery to surface waters

The PSYCHIC process-based model for predicting sediment and phosphorus (P) transfer within catchments uses spatial data on soil-P derived from the National Soil Inventory (NSI) data set. These soil-P values are based on bulked 0–15 cm depth and do not account for variations in soil-P with depth. We describe the depth distribution of soil-P (total and Olsen) in grassland and arable soils for the dominant soil types in the two PSYCHIC study catchments: the Avon and the Wye, UK. There were clear variations in soil-P (particularly Olsen-P) concentrations with depth in untilled grassland soils while concentrations of total-P were broadly constant within the plough layer of arable soils. Concentrations of Olsen-P in arable soils, however, exhibited maximum values near the soil surface reflecting surface applications of fertilisers and manures between consecutive ploughing events. When the soil-P concentrations for the surface soil (0–5 cm average) were compared to both the profile-averaged (0–15 cm) and the NSI (0–15 cm) values, those for the surface soil were considerably greater than those for the average 0–15 cm depth. Modelled estimates of P loss using the depth-weighted average soil-P concentrations for the 0–5 cm depth layer were up to 14% greater than those based on the NSI data set due to the preferential accumulation of P at the soil surface. These findings have important implications for the use of soil-P data (and other data) in models to predict P losses from land to water and the interpretation of these predictions for river basin management.