Modeling the land surface reflectance for optical remote sensing data in rugged terrain

A model for topographic correction and land surface reflectance estimation for optical remote sensing data in rugged terrian is presented.Considering a directional-directional reflectance that is used for direct solar irradiance correction and a hemispheric-directional reflectance that is used for atmospheric diffuse irradiance and terrain background reflected irradiance correction respectively,the directional reflectance-based model for topographic effects removing and land surface reflectance calculation is developed by deducing the directional reflectance with topographic effects and using a radiative transfer model.A canopy reflectance simulated by GOMS model and Landsat/TM raw data covering Jiangxi rugged area were taken to validate the performance of the model presented in the paper.The validation results show that the model presented here has a remarkable ability to correct topography and estimate land surface reflectance and also provides a technique method for sequently quantitative remote sensing application in terrain area.     

Inversion of magnetotelluric data in the presence of strong bathymetry/topography

We describe an algorithm for inverting magnetotelluric (MT) data in the presence of strong bathymetry or topography. Instead of correcting distortions due to bathymetry or topography we incorporate them directly into the inversion. To achieve a high accuracy in computing MT responses we use finite‐difference approximations that permit sloping discrete boundaries inside model elements. The same approach is applicable to any seafloor electromagnetic calculation and can also be used to incorporate steep topography on land. We test our approach on various topographic features and compare our results to that of a finite element approach. Finally, we present inversion examples that illustrate the effectiveness of our inversion algorithm in recovering true subsurface structures in the presence of strong bathymetry and topography.     

地形效应影响下地震动参数与斜坡稳定性的相关性研究

地形对地震动的影响比较复杂, 考虑地形放大效应的地震滑坡稳定性分析需要选择合适的地震动参数. 本文使用自贡地形影响台阵记录到的2008年汶川M_S8.0地震主震加速度记录, 分析了地震动峰值加速度、 阿里亚斯烈度以及90%能量持时随地形高度的变化, 探讨了地形效应作用下峰值加速度和阿里亚斯烈度与地震动作用下斜坡稳定性的相关性. 结果表明: ① 地形场地对峰值加速度和阿里亚斯烈度均有显著的放大效应. 地形放大效应较为复杂, 其整体上随台站高度的增加而增大, 水平向的放大效应大于竖直向. 水平向峰值加速度的放大系数为1.1—1.8, 阿里亚斯烈度的放大系数为1.2—3.3; 竖直向相应放大系数分别为1.1—1.3和1.2—1.7. ② 地形对地震动持时也有一定的放大效应, 但不同高度、 不同分量的放大效应没有显著差异, 其放大系数均约为1.3. ③ 阿里亚斯烈度和峰值加速度均能很好地表征地形对地震动的影响, 与地震动对斜坡稳定性的影响具有很强的相关性. 与峰值加速度相比, 阿里亚斯烈度综合了地震动的多方面特征, 可以更好地表征地形对地震动的影响, 与地震动作用下斜坡稳定性的相关性更强.      

考虑电流型畸变的大地电磁三维反演

大地电磁测深(MT)的观测数据易受到由近地表小尺度非均匀体或地形起伏引起的电流型畸变干扰,消除或压制这种干扰对获取可靠的深部电性结构至关重要.当区域结构为二维时,电流型畸变可采用张量分解等方法予以消除或压制.当区域结构为三维时,畸变问题更加复杂和严重,传统张量分解方法往往效果不佳或无效,严重地制约了MT三维反演技术的实用性.对此,本文提出一种考虑电流型畸变的MT三维反演算法,将完整的电流型畸变参数引入到目标函数,并采用非线性共轭梯度法与电阻率参数同时反演,从而达到压制畸变的目的.该算法有两个关键点:一是通过分析实测数据所遭受畸变的分布特征,在目标函数中对其进行有效约束;二是在迭代过程中,通过自适应地调整双正则化因子保障算法的稳定和效率.理论模型测试结果显示,常规三维反演算法不能合理解释数据中的畸变成分,而只能通过引入虚假异常体强制地拟合受畸变数据,从而造成电阻率模型严重失真.与之相比,本文算法能够在反演中自动求解各测点所受到的畸变,获得更接近真实的电阻率模型.     

Chrome-rich clinopyroxene in orthopyroxenite from Maowu, Dabie Mountains, central China: A second record and its implications for petrogenesis

A second occurrence of chrome-rich clinopyroxene has been discovered as inclusions in orthopyroxene in orthopyroxenite from Maowu, the Dabie Mountains, Central China. The average formula for chrome-rich clinopyroxene can be expressed as (Na_0.39Ca_0.54)_0.93(Mg_0.57Fe²⁺_0.06Fe³⁺_0.01Cr_0.24Al_0.15)_1.03Si_2.02O₆, with a maximum amount of kosmochlor component of 28.52 mol%. The unit cell parameters obtained from a single-crystal are a  = 9.614 Å, b  = 8.800 Å, c =  5.240 Å, β = 106.59°, space group C2 / c . The indices of refraction are α = 1.697, β = 1.704, γ = 1.726. Chrome-rich clinopyroxene, which coexists with chromite, chromian rutile and chromian pyrope, crystallized at a temperature of 1025 °C and very high pressure, and therefore represents a mantle relic. Together with the appearance of low-pressure inclusion mineral assemblage and the estimation of physical–chemical conditions for matrix minerals, the Maowu eclogite–ultramafic complex is considered to be formed during ultrahigh-pressure metamorphism from the mantle-derived protolith.     

Gravimetric survey terrain correction using linear analytical approximation

Various methods for computing the terrain correction in a high‐precision gravity survey are currently available. The present paper suggests a new method that uses linear analytical terrain approximations. In this method, digital terrain models for the near‐station topographic masses are obtained by vectorizing scan images of large‐scaled topographic maps, and the terrain correction computation is carried out using a Fourier series approximation of discrete height values. Distant topography data are represented with the help of digital GTOPO30 and Shuttle Radar Topography Mission cartographic information. We formulate linear analytical approximations of terrain corrections for the whole region using harmonic functions as the basis of our computational algorithm. Stochastic modelling allows effective assessment of the accuracy of terrain correction computation. The Perm Krai case study has shown that our method makes full use of all the terrain data available from topographic maps and digital terrain models and delivers a digital terrain correction computed to a priori precision. Our computer methodology can be successfully applied for the terrain correction computation in different survey areas.     

复杂地形条件下航空伽玛能谱地形改正方法探讨

应用矩形辐射体航空伽玛辐射场理论,研究应用于复杂地形条件下的航空伽玛能谱资料的地形改正方法,并针对该方法进行验证. 本文利用航空物探测量过程中获得的DTM数据（达到了地形改正所需的地形起伏数据精度）,应用矩形辐射体航空伽玛辐射场理论,根据地面辐射体与航空伽玛场分布之间的正演关系,对航空伽玛能谱解释方法的原理进行了正演分析和反演推导. 探索出一种按影响角进行地形改正的方法. 该方法特点是适合于任意飞行方式(缓地形和水平飞行均可)和任意地形条件航空伽玛能谱的逐点地形改正. 结果表明该地形改正方法能够基本消除航空伽玛能谱测量中的地形起伏产生的影响,经地形修正后的航空伽玛能谱异常能较正确地反映地面辐射体的真实情况.     

A geoid solution for airborne gravity data

Airborne gravity data is usually attached with satellite positioning of data points, which allow for the direct determination of the gravity disturbance at flight level. Assuming a suitable gridding of such data, Hotine’s modified integral formula can be combined with an Earth Gravity Model for the computation of the disturbing potential (T) at flight level. Based on T and the gravity disturbance data, we directly downward continue T to the geoid, and we present the final solution for the geoid height, including topographic corrections. It can be proved that the Taylor expansion of T converges if the flight level is at least twice the height of the topography, and the terrain potential will not contribute to the topographic correction. Hence, the simple topographic bias of the Bouguer shell yields the only topographic correction. Some numerical results demonstrate the technique used for downward continuation and topographic correction.     

USING UNMANNED AERIAL VEHICLE PHOTOGRAMMETRY TECHNOLOGY TO OBTAIN QUANTITATIVE PARAMETERS OF ACTIVE TECTONICS

With the development of photogrammetry technology and the popularity of unmanned aerial vehicles (UAVs)technology in recent years, using UAV photogrammetry technology to rapidly acquire high precision and high resolution topographic and geomorphic data on the fault zone has gradually become an important technical means. This paper first summarizes the basic principle and workflow of a new digital photogrammetry technology, SfM (Structure from Motion), which is simple, efficient and low cost. Using this technology, we conducted aerial image acquisition and data processing for a typical fault landform on the northern of Caka Basin in Qinghai. The digital elevation model (DEM)with 6.1cm/pix resolution is generated and the density of point cloud is as high as 273 points/m². The coverage area is 0.463km². Further, the terrain and slope data parallel to the fault direction are extracted by topographic analysis method, and combined with the contour map and the slope diagram generated by the DEM, a fine interpretation and quantitative study of complex multilevel geomorphic surfaces is carried out. Finally, based on the results of sophisticated interpretation of geomorphology, we got the vertical displacements of the T₁ terrace to the T₃ terrace as (1.01±0.06)m, (1.37±0.13)m and (3.10±0.11)m, and the minimum vertical displacements of the T₄ terrace and the T₅ terrace as (3.77±0.14)m and (5.46±0.26)m, respectively, through the topographic profile data extracted by DEM. Such vertical displacement parameters are difficult to obtain directly by traditional remote sensing images, which shows the great application prospect of UAV photogrammetry technology in the quantitative study of active tectonics.     

Seismic refraction traveltime inversion for static corrections in a glaciated shield rock environment: a case study

A traveltime inversion technique is applied to model the upper ∼40 m of the subsurface of a glaciated shield rock area in order to calculate static corrections for a multi-azimuth multi-depth walk-away vertical seismic profile and a surface seismic reflection profile. First break information from a seismic refraction survey is used in conjunction with a ray-tracing program and an iterative damped least-squares inversion algorithm to create a two-dimensional model of the subsurface. The layout of the seismic survey required crooked seismic lines and substantial gaps in the source and receiver coverage to be accounted for. Additionally, there is substantial topographical variation and a complex geology consisting of glaciofluvial sediment and glacial till overlying a crystalline bedrock. The resolution and reliability of the models is measured through a parameter perturbation technique, normalized χ² values, root means square traveltime residuals and comparison to known geology.     

Experiments on the dynamic correction of magnetotelluric and magnetovariational data

Pathways of further development of the dynamic correction technique suggested by E. Fainberg are discussed. This technique allows the correction of magnetotelluric (MT) and magnetovariational (MV) data distorted by subsurface irregularities within a broad range of frequencies beyond the static shift interval. This paper addresses new algorithms for the dynamic correction of tensor MT and MV response functions. Moreover, experiments on the dynamic correction of MT and MV data calculated for models containing heterogeneities in sedimentary cover and infinite or half-infinite conductive crustal prism are described. Model experiments show that the dynamic correction can ensure sufficient accuracy of the input data for the reliable solution of the inverse problem of magnetotellurics.     

基于三维正演的音频大地电磁阻抗相位不变量校正技术

复杂地形、地质条件的大地电磁数据解释容易出现假象,采用三维正演技术模拟地形和地表不均匀体的背景响应,对实测数据阻抗相位不变量进行校正,实现更准确的定性分析;对三维异常体模型的合成数据进行一维、二维多参数反演试算,以确定地形剧变区选择反演技术的最佳方案.合成数据的试反演结果显示一维反演水平切片假异常较多,二维反演能压制测向假异常,但不能压制走向的假异常,水平切片多出现测向条带.使用本文提出的阻抗相位不变量校正法扣除地形、地表背景响应,结合一维、二维反演,能使实际资料解释成果更加可靠.     

Origins of hydrocarbons in the Sagara oil field, central Japan

Abstract   We collected free-gas and in situ fluid samples up to a depth of 200.6 m from the Sagara oil field, central Japan (34°44'N, 138°15'E), during the Sagara Drilling Program (SDP) and measured the concentrations and stable carbon isotopic compositions of CH₄ and C₂H₆ in the samples. A combination of the CH₄/C₂H₆ ratios with the carbon isotope ratios of methane indicates that the hydrocarbon gases are predominantly of thermogenic origin at all depths. The isotope signature of hydrocarbon gases of δ¹³      < δ¹³     suggests that these gases in the Sagara oil field are not generated by polymerization, but by the decomposition of organic materials.     

Geographical distribution of helium isotope ratios in northeastern Japan

In order to study the precise geographical distribution of helium isotope ratios in northeastern Japan and compare it with geophysical data, we collected 43 gas and water samples from hot and mineral springs in the region where the ratio had never been reported, and measured the ³He/⁴He and ⁴He/²⁰Ne ratios of these samples. It was found that the ³He/⁴He ratios show clear contrasts between the forearc and the back-arc regions in the Tohoku district in northeastern Japan. In the forearc region, the ratios are smaller than 1 R_A (1 R_A = 1.4 × 10⁻⁶; R_A means the ³He/⁴He ratio of the atmosphere). On the other hand, those along the volcanic front and in the back-arc region are apparently higher. Moreover, we found a variation in the ³He/⁴He ratios along the volcanic front. In Miyagi Prefecture (38–39°N), the ratios range from 2 to 5 R_A. On the other hand, the ratios are less than 1 R_A in and around the southern border between Iwate and Akita Prefectures (39–39.5°N). Comparing the distribution of helium isotope ratios to results of recent geophysical studies, we found that the features in geographical distribution of helium isotope ratios are similar to those of seismic low-velocity zone distributions and high Qp⁻¹ distributions in the uppermost mantle. These observations strongly suggest that the helium isotope ratios reflect the distribution of melts in the uppermost mantle and are a useful tool for investigating the origin, behavior, and distribution of deep fluids and melts.     

Empirical Assessment of Ground Water–to-Indoor Air Attenuation Factors for the CDOT-MTL Denver Site

The area surrounding the Colorado Department of Transportation Materials Testing Laboratory in Denver was the subject of intense investigation, involving the collection of thousands of ground water, soil-gas, and indoor air samples in order to investigate indoor air impacts associated with a subsurface release of chlorinated solvents. The preremediation portion of that data set is analyzed and reduced in this work to ground water–to-indoor air attenuation factors (α_gw= the ratio of the measured indoor air concentration to the soil-gas concentration predicted to be in equilibrium with the local ground water concentration). The empirical α_gw values for this site range from about 10⁻⁶ to 10⁻⁴ with an overall average of 3 × 10⁻⁵ (μg/L indoor air)/(μg/L soil gas). The analysis of this data set highlights the need for a thorough data review and data screening when using large data sets to derive empirical relationships between subsurface concentrations and indoor air. More specifically, it is necessary to identify those parts of the data that contain a strong vapor intrusion pathway signal, which generally will require concentrations well above reported detection levels combined with spatial or temporal correlation of subsurface and indoor concentrations.     

基于非结构化网格的任意复杂2D RMT有限元模拟

射频大地电磁法(RMT)是以潜艇天线发射的射线源等作为场源的一种地球物理勘探方法,目前被广泛应用于近地表工程和环境地球物理勘探.RMT数据解释常采用基于静态假设的大地电磁法(MT)程序,往往会反演出不真实的浅层目标体.为解决这一长期困扰RMT资料解释的问题,本文实现了考虑位移电流的RMT有限元数值模拟.为了处理任意复杂模型,如起伏地形,非结构的三角形网格被用于离散RMT模型.首先通过算例对比,验证了程序的正确性和可靠性.通过Dike模型讨论了空气层厚度对RMT数值解的影响,结果表明当空气层厚度大于1/4波长即可满足精度要求.以山脊模型为例计算了位移电流对RMT响应的影响,表明位移电流的影响会随着山脊高程的增加而增大.最后通过舒家店实际模型进一步验证了位移电流在RMT中的重要性.     

Azimuth-dependent AVO in reservoirs containing non-orthogonal fracture sets

Azimuthal anisotropy in rocks can result from the presence of one or more sets of partially aligned fractures with orientations determined by the stress history of the rock. The symmetry of a rock with horizontal bedding that contains two or more non-orthogonal sets of vertical fractures may be approximated as monoclinic with a horizontal plane of mirror symmetry. For offsets that are small compared with the depth of the reflector, the azimuthal variation in P-wave AVO gradient for such a medium varies with azimuth as     where φ is the azimuth measured with respect to the fast polarization direction for a vertically polarized shear wave. φ ₂ depends on both the normal compliance B _N and the shear compliance B _T of the fractures and may differ from zero if B _N B _T varies significantly between fracture sets. If B _N B _T is the same for all fractures,     and the principal axes of the azimuthal variation in P-wave AVO for fixed offset are determined by the polarization directions of a vertically propagating shear wave. At larger offsets, terms in     and     are required to describe the azimuthal variation in AVO accurately. φ ₄ and φ ₆ also depend on B _N B _T. For gas-filled open fractures     but a lower value of B _N B _T may result from the presence of a fluid with non-zero bulk modulus.     

非平坦地形条件下电阻率三维反演

本文实现了非平坦地形条件下电阻率三维反演,讨论了几种消除反演中地形影响的方法. 结果表明,只有将地形直接带入反演算法中,进行带地形电阻率三维反演才能有效消除地形影响及其对反演结果的偏差,得到与地下电性结构相符的反演结果.     

黄土高原复杂地形上高质量湍流通量数据获取方法

利用兰州大学半干旱气候与环境观测站(Semi-Arid Climate and Environment Observatory of Lanzhou University,SACOL)湍流观测资料,分析了二次坐标旋转(double rotation,DR)、平面拟合(planar fit,PF)和分风区平面拟合(fetch planar fit,FPF)在复杂地形上的适用性,总结出一套适用于SACOL的总体湍流特征参数化方案.经过超声虚温订正、坐标旋转、空气密度脉动订正以及平稳性检验、总体湍流特征检验、总体质量分级处理,摩擦速度(u*)、感热通量、潜热通量、CO_2通量高质量数据所占比例分别为45%~62%、66%~68%、62%~65%、52%~54%.采用DR得到的高质量数据比例与采用PF相比,u*提高了17%,后三种通量略降低2%~3%.PF和FPF两种结果的差别主要体现在u*上,只考虑主导风向数据DR得到的u*质量仍最好.综合兼顾数据质量和计算工作量,在复杂地形上处理湍流观测资料的最优坐标旋转方法是DR.