Computation of Large Anisotropic Seismic Heterogeneities (CLASH)

A general tomographic technique is designed in order (i) to operate in anisotropic media; (ii) to account for the uneven seismic sampling and (iii) to handle massive data sets in a reasonable computing time. One modus operandi to compute a 3-D body wave velocity model relies on surface wave phase velocity measurements. An intermediate step, shared by other approaches, consists in translating, for each period of a given mode branch, the phase velocities integrated along ray paths into local velocity perturbations. To this end, we develop a method, which accounts for the azimuthal anisotropy in its comprehensive form. The weakly non-linear forward problem allows to use a conjugate gradient optimization. The Earth's surface is regularly discretized and the partial derivatives are assigned to the individual grid points. Possible lack of lateral resolution, due to the inescapable uneven ray path coverage, is taken into account through the a priori covariances on parameters with laterally variable correlation lengths. This method allows to efficiently separate the 2ψ and the 4ψ anisotropic effects from the isotropic perturbations. Fundamental mode and overtone phase velocity maps, derived with real Rayleigh wave data sets, are presented and compared with previous maps. The isotropic models concur well with the results of Trampert & Woodhouse. Large 4ψ heterogeneities are located in the tectonically active regions and over the continental lithospheres such as North America, Antarctica or Australia. At various periods, a significant 4ψ signature is correlated with the Hawaii hotspot track. Finally, concurring with the conclusions of Trampert & Woodhouse, our phase velocity maps show that Rayleigh wave data sets do need both 2ψ and 4ψ anisotropic terms.     

Spectral Simulation and Conditioning to Local Continuity Trends in Geologic Modeling

Spectral simulation has gained application in building geologic models due to the advantage of better honoring the spatial continuity of petrophysical properties, such as reservoir porosity and shale volume. Distinct from sequential simulation methods, spectral simulation is a global algorithm in the sense that a global density spectrum is calculated once and the inverse Fourier transform is performed on the Fourier coefficient also only once to generate a simulation realization. The generated realizations honor the spatial continuity structure globally over the whole field instead of only within a search neighborhood, as with sequential simulation algorithms. However, the disadvantage of global spectral simulation is that it traditionally cannot account for the local information such as the local continuity trends, which are often observed in reservoirs and hence are important to be accounted for in geologic models. This disadvantage has limited wider application of spectral simulation in building geologic models. In this paper, we present ways of conditioning geologic models to the relevant local information. To account for the local continuity trends, we first scale different frequency components of the original model with local-amplitude spectrum ratios that are specific to the local trend. The sum of these scaled frequency components renders a new model that displays the desired local continuity trend. The implementation details of this new method are discussed and examples are provided to illustrate the algorithm.     

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万地形图地物判读。     

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.     

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.