Fast computation of Hankel Transform using orthonormal exponential approximation of complex kernel function

The computation of electromagnetic (EM) fields, for 1-D layered earth model, requires evaluation of Hankel Transform (HT) of the EM kernel function. The digital filtering is the most widely used technique to evaluate HT integrals. However, it has some obvious shortcomings. We present an alternative scheme, based on an orthonormal exponential approximation of the kernel function, for evaluating HT integrals. This approximation of the kernel function was chosen because the analytical solution of HT of an exponential function is readily available in literature. This expansion reduces the integral to a simple algebraic sum. The implementation of such a scheme requires that the weights and the exponents of the exponential function be estimated. The exponents were estimated through a guided search algorithm while the weights were obtained using Marquardt matrix inversion method. The algorithm was tested on analytical HT pairs available in literature. The results are compared with those obtained using the digital filtering technique with Anderson filters. The field curves for four types (A-, K-, H-and Q-type) of 3-layer earth models are generated using the present scheme and compared with the corresponding curves obtained using the Anderson sc heme. It is concluded that the present scheme is more accurate than the Anderson scheme     

Seismic Hazards Assessment for Radioactive Waste Disposal Sites in Regions of Low Seismic Activity

A comparative tectonic quiescence and lack of earthquakes make the stable centres of continents attractive for siting long-term radioactive waste storage facilities. The low rates of deformation in such regions, however, make it difficult to characterize their long-term seismotectonic behaviour, leading to uncertain estimates for the very low probability hazard estimates required by society. In an attempt to overcome the deficiency of both contemporary seismicity and paleoseismic data in central Canada, we have used earthquake histories from regions with similar seismotectonic characteristics from around the world. Substituting space for time, we estimate a long-term rate per 10⁶ km² of 0.004 magnitude ≥6 earthquakes per annum, of which 33–100% might rupture to the surface.     

我国西部1∶5万地形图空白区地物解译研究——以西藏那曲其香错地区为例

我国西部部分地区1∶5万地形图数据尚未完成,给区域和国家经济发展带来极大不便。空间技术、计算机技术和信息技术的高速发展为遥感技术应用提供了技术支持,高分辨率影像获取能力的提高为大比例尺地物特征的提取提供了数据可能。本文采用SPOT5遥感影像,以1∶5万地形图图式标准为依据,进行青藏高原其香错地区地物判读,结果表明SPOT5能基本满足1∶5万地形图制图要求。为进一步提高判读结果,可采用多时相、多分辨率影像进行1∶5万地形图地物判读。     

地学自由曲面径向基函数网络重建的详细算法

径向基函数(RadialBasisFunction,简称RBF)神经网络是一种理想的地学离散数据网格化工具,能够适应各种不同分布形式和边界条件的数据,收敛速度较快,可以逼近任何复杂曲面。这里详细介绍了RBF神经网络的算法。适当径向基函数的形式和偏差系数是使用RBF神经网络作地学曲面重建的关键。大量的实际数据验证结果表明,当选用Gauss型径向基函数时,一般可获得比较理想的网格化效果,同时具备"曲面平滑"和"拟合度高"的特点。     

Analysis and classification of hyperspectral data for mapping land degradation: An application in southern Spain

Desertification is a severe stage of land degradation, manifested by “desert-like” conditions in dryland areas. Climatic conditions together with geomorphologic processes help to mould desert-like soil surface features in arid zones. The identification of these soil features serves as a useful input for understanding the desertification process and land degradation as a whole. In the present study, imaging spectrometer data were used to detect and map desert-like surface features. Absorption feature parameters in the spectral region between 0.4 and 2.5 μm wavelengths were analysed and correlated with soil properties, such as soil colour, soil salinity, gypsum content, etc. Soil groupings were made based on their similarities and their spectral reflectance curves were studied. Distinct differences in the reflectance curves throughout the spectrum were exhibited between groups. Although the samples belonging to the same group shared common properties, the curves still showed differences within the same group.Characteristic reflectance curves of soil surface features were derived from spectral measurements both in the field and in the laboratory, and mean reflectance values derived from image pixels representing known features. Linear unmixing and spectral angle matching techniques were applied to assess their suitability in mapping surface features for land degradation studies. The study showed that linear unmixing provided more realistic results for mapping “desert-like” surface features than the spectral angle matching technique.     

土壤热异常对地表能量平衡影响初探

将来自土壤深部的热通量引入off line的陆面过程模式 (NCAR—LSM ) ,通过长达 2a的数值试验对比分析了它对各层次土壤温度和地表能量平衡的影响。　　在土壤底部引入 5W /m2 的热通量使底层土壤显著升温 ,但升温随着接近表层而迅速衰减。积分 3个月后 ,由地下进入地表的热流量增幅可达 1W/m2 以上 ,并持续增大到 5W /m2 ,地表最大升温约 0 .5K ,同时地表感热、蒸发潜热及长波辐射通量均有 1W /m2 左右的正异常 ;若将土壤热传导系数放大一个量级以加速热量交换 ,则地表升温提高到 1K以上 ,长波辐射增加 3W /m2 以上 ,超过了气溶胶全球平均的辐射效应。结果表明 :一定量值的土壤热异常对地表能量平衡和短期气候变化 (10 -1～ 10 1a)有着不可忽略的影响。同时 ,深入的资料分析、完善的陆面过程模式以及它与大气模式的耦合试验也是亟待进行的相关工作。     

土壤热异常影响地表能量平衡的个例分析和数值模拟

The statistical relationship between soil thermal anomaly and short-term climate change is presented based on a typical case study. Furthermore, possible physical mechanisms behind the relationship are revealed through using an off-line land surface model with a reasonable soil thermal forcing at the bottom of the soil layer.In the first experiment, the given heat flux is 5 W m-2 at the bottom of the soil layer (in depth of 6.3 m)for 3 months, while only a positive ground temperature anomaly of 0.06℃ can be found compared to the control run. The anomaly, however, could reach 0.65℃ if the soil thermal conductivity was one order of magnitude larger. It could be even as large as 0.81℃ assuming the heat flux at bottom is 10 W m-2. Meanwhile, an increase of about 10 W m-2 was detected both for heat flux in soil and sensible heat on land surface, which is not neglectable to the short-term climate change. The results show that considerable response in land surface energy budget could be expected when the soil thermal forcing reaches a certain spatial-tem poral scale. Therefore, land surface models should not ignore the upward heat flux from the bottom of the soil layer. Moreover, integration for a longer period of time and coupled land-atmosphere model are also necessary for the better understanding of this issue.     

近百年中国东部夏季降水的时空变率

利用中国东部25°N以北28个站1880-1999年夏季季降水序列,用旋转复经验正交函数（RCEOF）方法,研究了中国东部地区百年干湿的时空演变规律。结果表明,夏季降水空间变率大值区依次为:长江中下游地区、淮河流域、江南、华北、西南及东北。除西南外的5个关键区大体上反映了从6月到8月夏季雨带自南向北椎进所滞留的地区。旋转空间位相分布揭示了长江中下游地区、江南、东北的旱涝异常主要表现为驻波振动特征;而淮河流域、华北、西南地区显示出降水异常信号具有部分的行波特征。尤其第4空间模显示出旱涝异常信号从东北南部可沿着黄淮下游传到长江下游地区。对于近百年中国东部地区夏季于湿变化,长江中下游地区、淮河流域、华北及东北四个地区都存在20-25年时间尺度的周期振荡;长江中下游地区及华北地区都存在准60年时间尺度的振荡周期;东北地区主要表现出36年时间尺度的振荡周期;淮河流域存在明显的70-80年时间尺度的振荡周期;华北地区存在的11年时间尺度的振荡周期恰好与太阳黑子活动的11年周期相一致。在年代际时间尺度（包括次年代际时间尺度）上,长江中下游、淮河流域及华北地区的夏季降水的变化与太阳活动有显著的正相关。     

Atmospheric Surface-Layer Processes During the Total Solar Eclipse Of 11 August 1999

The central path of the total solar eclipse (TSE) of 11 August 1999 mostly passedthrough a region of active monsoon in India, with the eclipse ending around localsunset. Measurements in the surface layer (SL) were carried out close to the eclipseaxis at Akola (20°42' N, 77°2' E) in central India. The globalradiation flux reaching the surface vanishes around totality at 1803:24 (LT), followedby a small recovery before again dipping to zero at sunset. The temperatures in the SL, and subsoil at 50-mm depth, show a local minimum with a lag of about 10 min after the second contact, whereas the lag appears to vanish when the temperature series is detrended. The SL exhibits near-neutral, though generally stable, conditions from about 1500 hr itself due to heavy cloud cover followed by the eclipse-induced cooling of the surface. The wind component across the eclipse axis vanishes at totality, the wind vector aligning with the azimuth of the traversing moon shadow. The deceleration of the mean flow can be due to the combined effect of the colder surface and downward heat flux, where the locally altered horizontal temperature gradients may cause, as in this instance, the cessation of the cross flow.The oscillations in temperature and wind that show significant peaks, around 90–100 min as well as the semi-period of the eclipse near 60 min, persist for several hours past the eclipse event. A fresh analysis of the published data on the TSE of 7 March 1970 shows spectral peaks in the temperature nearly coincident with those already reported from the surface pressure records. The oscillations in the SL fields during the two TSE events are very similar implying that the source mechanism is also likely to be the same in both the cases.     

The Production of Methanol by Flowering Plants and the Global Cycle of Methanol

Methanol has been recognised as an important constituent of the background atmosphere, but little is known about its overall cycle in the biosphere/atmosphere system. A model is proposed for the production and emission to the atmosphere of methanol by flowering plants based on plant structure and metabolic properties, particularly the demethylation of pectin in the primary cell walls. This model provides a framework to extend seven sets of measurements of methanol emission rates to the global terrestrial biosphere. A global rate of release of methanol from plants to the atmosphere of 100 Tg y^–1 is calculated. A separate model of the global cycle of methanol is constructed involving emissions from plant growth and decay, atmospheric and oceanic chemical production, biomass burning and industrial production. Removal processes occur through hydroxyl radical attack in the atmosphere, in clouds and oceans, and wet and dry deposition. The model successfully reproduces the methanol concentrations in the continental boundary-layer and the free atmosphere, including the inter-hemispheric gradient in the free atmosphere. The model demonstrates a new concept in global biogeochemistry, the coupling of plant cell growth with the global atmospheric concentration of methanol. The model indicates that the ocean provides a storage reservoir capable of holding at least 66 times more methanol than the atmosphere. The ocean surface layer reservoir essentially buffers the atmospheric concentration of methanol, providing a physically based smoothing mechanism with a time constant of the order of one year.