近10年西藏高原雪线时空变化及其与气象因素关系分析

本文利用2000年3月-2011年2月西藏地区的MODIS雪盖产品数据、DEM数据以及地面气象观测数据,结合GIS空间分析方法,分析了西藏地区不同自然区划地带下雪线的时空变化特征及其与气象因素的关系。研究表明：西藏及各区域年平均雪线波动变化比较平稳,全区年平均雪线为4848.6m,呈微弱上升趋势,线性倾向率为6.54m/10a;各季节平均雪线中,秋季雪线的变化对年平均贡献最大,二者相关系数达0.796。冬季雪线呈下降趋势（相关系数为-0.625）,其余三季则均表现为上升趋势,但均不显著;除东喜马拉雅南翼山地雪线逐月变化波动明显外（标准差为60.3m）,其余均表现为平缓波动形势;西藏地区的雪线空间分布基本上表现为由东南向西北方向逐步升高的态势,其中东南部和西北部雪线分布密集且复杂。中部雪线则相对较稀疏,其高、低值区分别与山脉和河谷分布相对应;整体上,西藏雪线与气温正相关,与降水量负相关,但是各区域四季雪线与气温、降水量之间又存在差异。雪线是积雪各要素特征变化最为敏感的指示器,研究西藏高原雪线的时空分布特征及其与气象因素之间的关系,对了解西藏高原乃至整个青藏高原的气候变化具有重要的意义。     

A new method to calculate the resonant response of a sunspot model atmosphere to magneto-atmospheric waves

In order to understand the observed oscillations in sunspots we present a new method for calculating the resonant response of a realistic semi-empirical model of the sunspot umbral atmosphere and subphotosphere to magneto-atmospheric waves in a vertical magnetic field. The depth dependence of both the adiabatic coefficient and the turbulent pressure is taken into account. This requires an extension of the wave equations by Ferraro & Plumpton (1958). We compare the coefficients of wave transmission, re flection, and conversion between fast mode and slow mode waves for different assumptions, compare the results with those from earlier modelling efforts, and point out possible sources of mistakes. The depth dependence of the adiabatic coefficient strongly influences the resulting spectrum of resonance frequencies. The condition of a conservation of wave flux is violated if the depth dependence of the turbulent pressure is not properly considered.     

Scattering of internal waves in a two-layer fluid flowing through a channel with small undulations

Using two-dimensional linear water wave theory, we consider the problem of normal water wave (internal wave) propagation over small undulations in a channel flow consisting of a two-layer fluid in which the upper layer is bounded by a fixed wall, an approximation to the free surface, and the lower one is bounded by a bottom surface that has small undulations. The effects of surface tension at the surface of separation is neglected. Assuming irrotational motion, a perturbation analysis is employed to calculate the first-order corrections to the velocity potentials in the two-layer fluid by using Green’s integral theorem in a suitable manner and the reflection and transmission coefficients in terms of integrals involving the shape function c(x) representing the bottom undulation. Two special forms of the shape function are considered for which explicit expressions for reflection and transmission coefficients are evaluated. For the specific case of a patch of sinusoidal ripples having the same wave number throughout, the reflection coefficient up to the first order is an oscillatory function in the quotient of twice the interface wave number and the ripple wave number. When this quotient approaches one, the theory predicts a resonant interaction between the bed and the interface, and the reflection coefficient becomes a multiple of the number of ripples. High reflection of the incident wave energy occurs if this number is large. Again, when a patch of sinusoidal ripples having two different wave numbers for two consecutive stretches is considered, the interaction between the bed and the interface near resonance attains in the neighborhood of two (singular) points along the x-axis (when the ripple wave number of the bottom undulation become approximately twice as large as the interface wave number). The theoretical observations are presented in graphical form.     

Style of fluid flow and deformation in and around an ancient out-of-sequence thrust: An example from the Nobeoka Tectonic Line in the Shimanto accretionary complex, Southwest Japan

To understand the characteristics of deformation of an out-of-sequence thrust (OST) and the style of fluid flow along it, we investigated the Nobeoka Tectonic Line, which has been interpreted as a deep OST (7–9 km), in the Shimanto accretionary complex, Southwest Japan. The shear zone in the footwall differs significantly in the along-strike direction not only in thickness, which varied from 100 to 300 m, but also in lithology and mineral vein development. These variations might reflect primarily differences in lithology; that is, the sandstone-dominant shear zone with a large amount of mineral veins precipitated in microcracks is relatively thick, whereas the shale-dominant shear zone with a small amount of veins and with textures indicating highly pressurized pore fluid, is thinner. By comparison with characteristics of a shallow OST (3–5 km), we conclude that the shallow OST has experienced repeated brittle failure with rapid slip and focused fluid flow whereas the deep OST has experienced both brittle and ductile deformation, followed by fluid flow of various styles, depending on the lithology.     

Geodynamic framework of large volcanic fields highlighted by SRTM DEMs: Method evaluation and perspectives exampled on three areas from the Cameroon Volcanic Line

This study is a part of a wider investigation to evaluate how much can be learnt by using low-cost methods such as systematic moderate-resolution remote sensing to map and quantify geological structures at the regional scale on very large volcanic provinces only partly studied in the field. Volcanic-centre and cinder-cone distribution, faults and structural lineaments are mapped combining Shuttle Radar Topography Mission (SRTM), Digital Elevation Models (DEMs) and Landsat satellite images. As an example of the method, we present the interpretation of structural data and morphological features of three contrasted areas from the Cameroon Volcanic Line (Tombel graben, Upper Benue valley, and Ngaoundéré area) for which local field studies are available for comparison. At a local scale, this remote-sensing method of mapping displays good to excellent correlations with previously published data and, by itself, it allows one to constrain the structural setting of each area. Numerical treatment of vent and cinder-cone localisation can be related to tension fractures (T direction), whereas numerical treatment of the lineaments constrains the associated fault system to a single transtensional (strike-slip + extension) Riedel type fracture network. The first results on the Cameroon Volcanic Line are promising and could be used at a larger scale on numerous volcanic provinces for which field data are not yet available.     

Derivation and relative performance of strings of line elements for modeling (un)confined and semi-confined flow

In the analytic element method, strings of line-sinks may be used to model streams and strings of line-doublets may be used to model impermeable walls or boundaries of inhomogeneities. The resulting solutions are analytic, but the boundary conditions are met approximately. Equations for line elements may be derived in two ways: through integration of point elements (the integral solution) and through application of separation of variables in elliptical coordinates (the elliptical solution). Using both approaches, two sets of line elements are presented for four flow problems: line-sinks and line-doublets in (un)confined flow, and line-sinks and line-doublets in semi-confined flow. Elliptical line elements have the advantage that they do not need a far-field expansion for accurate evaluation far away from the element. The derivation of elliptical line elements is new and applicable to both (un)confined flow and semi-confined flow; only the resulting expressions for elliptical line elements for semi-confined flow have not been found in the current groundwater literature. Existing solutions for elliptical line elements for (un)confined flow were intended for the modeling of isolated features. Four examples are presented, one for each flow problem, to demonstrate that strings of elliptical line elements may be used to obtain accurate solutions; elliptical line-doublets for semi-confined flow can only be strung together in combination with two integral line-doublets. For a zigzag canal in (un)confined flow, a string of elliptical line-sinks performed better than a string of integral line-sinks of the same order when the smallest angle between two adjacent segments is less than 130°. Elliptical line-doublets performed better than integral line-doublets for a square inhomogeneity in a uniform, confined flow field; the difference was smaller for an octagonal inhomogeneity. For semi-confined flow, the difference between the integral and elliptical line-sinks was small when modeling a zigzag canal.     

基于LIC的海流动态模拟与可视化

线积分卷积(LIC)方法是一种针对矢量场数据的可视化技术,本文主要研究了LIC方法应用于海洋流场的可视化质量改进问题。主要包括以下几点改进:根据海流运动方向改变稀疏背景噪声的指向以使流线变得更加细长明亮;针对稀疏噪声中海流变化剧烈的区域流线会出现分叉的问题,采用忽略稀疏背景噪声流线追踪时非主流线部分的纹理点的积分计算方法来解决此问题,同时可以省去不必要纹理点的积分计算的开支,加快显示速度;采用赋予海流随机的初始相位的方法来解决稀疏噪声循环动画突兀的问题;通过流线的运动速度快慢或者流线的颜色来表示海流速度大小;通过对海流数据进行插值并改善边界流线追踪条件的方法可形成完整细致的全球海流图像。实验结果表明,该方法可以生成纹理对比强烈、细节清晰、动画连贯的高质量海流运动图像,具有良好的可视化效果。     

Scattering of water wave by a submerged horizontal plate and a submerged permeable breakwater

Based on a two-dimensional linear water wave theory, this study develops the boundary element method (BEM) to examine normally incident wave scattering by a fixed, submerged, horizontal, impermeable plate and a submerged permeable breakwater in water of finite depth. Numerical results for the transmission coefficients are also presented. In addition, the numerical technique's accuracy is demonstrated by comparing the numerical results with previously published numerical and experimental ones. According to that comparison, the transmission coefficient relies not only on the submergence of the horizontal impermeable plate and the height of the permeable breakwater, but also on the distance between horizontal plate and permeable breakwater. Results presented herein confirm that the transmission coefficient is minimum for the distance approximately equal to four times the water depth.     

含裂缝地层的过套管电阻率测井响应特征研究

由于石油压裂开采等开发措施的实施,会使在套管周围存在着许多被油气等高阻流体填充的微裂缝.尽管这些裂缝或孔隙可能很小,但对套管井电阻率测量会有非常大的影响,甚至会改变原地层的电阻率测井特征,因此裂缝测井响应的计算及考察对过套管电阻率测井十分重要.为解决裂缝测井响应的计算问题,本文提出了计算等效电阻的电流通量管模型,利用该电流通量管模型给出了含垂直裂缝地层横向电阻的计算方法,基于传输线方程法及地层电阻的过套管测量方法实现了含裂缝地层的过套管电阻率测井响应的数值计算,通过计算实例考察了地层裂缝对过套管电阻率测井响应的影响.算例表明:裂缝中的高阻流体对地层视电阻率测量结果会有较大的影响;环形裂隙比垂直裂缝有更大的电阻率测井响应.本文的研究为解决微裂缝过套管电阻率测井响应的计算这一关键技术问题提供了一种可行的计算与考察方法.     

中国气溶胶分布的地理学和气候学特征

中国人口地理分界线——胡焕庸线根据人口、地理,气候和经济等特点把中国(不包括港、澳、台地区)分为东、西两部分。用2000~2010年MODIS大气气溶胶光学厚度(AOD)资料,分析气溶胶分布的地理学和气候学特征后发现,胡焕庸线还可被视为中国气溶胶地理学的分界线,在其两侧气溶胶的性质和浓度都有明显差别。在人口稠密和海拔较低的东部,由人类活动产生的气溶胶为主,年平均AOD约为0.45;在西部,自然过程释放的气溶胶主导的AOD约为0.25。近10 a来东部AOD的年际间变化呈现增加趋势,西部AOD出现微弱减少的趋势。东部人为气溶胶年际间变化受亚洲季风影响。西部自然气溶胶年际间变化主要受沙漠地区沙尘气溶胶排放源的影响,沙尘天气过程主要控制其气溶胶的释放。