Triangulated spherical splines for geopotential reconstruction

We present an alternate mathematical technique than contemporary spherical harmonics to approximate the geopotential based on triangulated spherical spline functions, which are smooth piecewise spherical harmonic polynomials over spherical triangulations. The new method is capable of multi-spatial resolution modeling and could thus enhance spatial resolutions for regional gravity field inversion using data from space gravimetry missions such as CHAMP, GRACE or GOCE. First, we propose to use the minimal energy spherical spline interpolation to find a good approximation of the geopotential at the orbital altitude of the satellite. Then we explain how to solve Laplace’s equation on the Earth’s exterior to compute a spherical spline to approximate the geopotential at the Earth’s surface. We propose a domain decomposition technique, which can compute an approximation of the minimal energy spherical spline interpolation on the orbital altitude and a multiple star technique to compute the spherical spline approximation by the collocation method. We prove that the spherical spline constructed by means of the domain decomposition technique converges to the minimal energy spline interpolation. We also prove that the modeled spline geopotential is continuous from the satellite altitude down to the Earth’s surface. We have implemented the two computational algorithms and applied them in a numerical experiment using simulated CHAMP geopotential observations computed at satellite altitude (450 km) assuming EGM96 (n _max = 90) is the truth model. We then validate our approach by comparing the computed geopotential values using the resulting spherical spline model down to the Earth’s surface, with the truth EGM96 values over several study regions. Our numerical evidence demonstrates that the algorithms produce a viable alternative of regional gravity field solution potentially exploiting the full accuracy of data from space gravimetry missions. The major advantage of our method is that it allows us to compute the geopotential over the regions of interest as well as enhancing the spatial resolution commensurable with the characteristics of satellite coverage, which could not be done using a global spherical harmonic representation. The results in this paper are based on the research supported by the National Science Foundation under the grant no. 0327577.     

遥感影像混合像元分解中的端元选择方法综述

端元选择是进行遥感影像混合像元分解的首要步骤,也是最关键的步骤,其直接影响混合像元分解的精度。该文对遥感影像混合像元分解中端元产生的特定背景、当前存在的端元选择途径和端元选择方法进行综述,并通过对当前端元选择方法的分析,提出了选择或构造端元选择方法应遵循的几个原则。     

Influence of pore size and geometry on peat unsaturated hydraulic conductivity computed from 3D computed tomography image analysis

In organic soils, hydraulic conductivity is related to the degree of decomposition and soil compression, which reduce the effective pore diameter and consequently restrict water flow. This study investigates how the size distribution and geometry of air‐filled pores control the unsaturated hydraulic conductivity of peat soils using high‐resolution (45 µm) three‐dimensional (3D) X‐ray computed tomography (CT) and digital image processing of four peat sub‐samples from varying depths under a constant soil water pressure head. Pore structure and configuration in peat were found to be irregular, with volume and cross‐sectional area showing fractal behaviour that suggests pores having smaller values of the fractal dimension in deeper, more decomposed peat, have higher tortuosity and lower connectivity, which influences hydraulic conductivity. The image analysis showed that the large reduction of unsaturated hydraulic conductivity with depth is essentially controlled by air‐filled pore hydraulic radius, tortuosity, air‐filled pore density and the fractal dimension due to degree of decomposition and compression of the organic matter. The comparisons between unsaturated hydraulic conductivity computed from the air‐filled pore size and geometric distribution showed satisfactory agreement with direct measurements using the permeameter method. This understanding is important in characterizing peat properties and its heterogeneity for monitoring the progress of complex flow processes at the field scale in peatlands. Copyright © 2010 John Wiley & Sons, Ltd.     

热解石墨坩埚的应用

本文介绍本院研制的热解石墨坩埚的特性和应用情况。该坩埚具有纯度高、致密、耐高温、耐强酸强碱腐蚀、易洗涤等特点,已广泛应用于矿物、建材、玻璃等方面试样的前处理。与铂、银、镍、刚玉、石英等坩埚比较,具有引入杂质少成本低等优点。     

ELIMINATING COMPLEX EIGENVECTORS AND EIGENVALUES IN MULTIWAY ANALYSES USING THE DIRECT TRILINEAR DECOMPOSITION METHOD

The direct trilinear decomposition method(DTDM)is an algorithm for performing quantitative curveresolution of three-dimensional data that follow the so-called trilinear model,e.g.chromatography-spectroscopy or emission-excitation fluorescence.Under certain conditions complexeigenvalues and eigenvectors emerge when the generalized eigenproblem is solved in DTDM.Previouspublications never treated those cases.In this paper we show how similarity transformations can be usedto eliminate the imaginary part of the complex eigenvalues and eigenvectors,thereby increasing theusefulness of DTDM in practical applications.The similarity transformation technique was first used byour laboratory to solve the similar problem in the generalized rank annihilation method(GRAM).Because unique elution profiles and spectra can be derived by using data matrices from three or moresamples simultaneously,DTDM with similarity transformations is more efficient than GRAM in the casewhere there are many samples to be investigated.     

第二代航磁黄沙坪异常的分解及其矿产新空间

本文采用同时具有识別、分解、解释异常三重作用的三元解释法■,在黄沙坪异常负值区有效地提取了由未知的深源磁性体引起的剩余异常。该异常与已知的中、浅源磁性体分处于极为类同的成矿地质背景中。后者与隐伏岩体、铅锌多金属盲矿体之间有着密切的成生联系。这为前者间接找矿提供了有力的依据。定量解释结果中,两者间的磁矩等效关系,论证了已知矿床西侧深部空间里,仍存在找到又一个大型矿床的前景。     

Characteristics and analysis of the geomagnetic variations in regions around the Qiongzhou Strait

Ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ａｎｄ　ａｎａｌｙｓｉｓ　ｏｆ　ｔｈｅ　ｇｅｏｍａｇｎｅｔｉｃ　ｖａｒｉａｔｉｏｎｓ　ｉｎ　ｒｅｇｉｏｎｓ　ａｒｏｕｎｄ　ｔｈｅ　Ｑｉｏｎｇｚｈｏｕ　Ｓｔｒａｉｔ（范国华）（姚同起）（顾左文）（朱克佳）（陈伯舫）（冯戬...     

Field measurements of NO and NO2 emissions from fertilized and unfertilized soils

Field measurements of NO and NO₂ emissions from soils have been performed in Finthen near Mainz (F.R.G.) and in Utrera near Seville (Spain). The applied method employed a flow box coupled with a chemiluminescent NO _x detector allowing the determination of minimum flux rates of 2 g N m^-2 h^-1 for NO and 3 g m^-2 h^-1 for NO₂.The NO and NO₂ flux rates were found to be strongly dependent on soil surface temperatures and showed strong daily variations with maximum values during the early afternoon and minimum values during the early morning. Between the daily variation patterns of NO and NO₂, there was a time lag of about 2 h which seem to be due to the different physico-chemical properties of NO and NO₂. The apparent activation energy of NO emission calculated from the Arrhenius equation ranged between 44 and 103 kJ per mole. The NO and NO₂ emission rates were positively correlated with soil moisture in the upper soil layer.The measurements carried out in August in Finthen clearly indicate the establishment of NO and NO₂ equilibrium mixing ratios which appeared to be on the order of 20 ppbv for NO and 10 ppbv for NO₂. The soil acted as a net sink for ambient air NO and NO₂ mixing ratios higher than the equilibrium values and a net source for NO and NO₂ mixing ratios lower than the equilibrium values. This behaviour as well as the observation of equilibrium mixing ratios clearly indicate that NO and NO₂ are formed and destroyed concurrently in the soil.Average flux rates measured on bare unfertilized soils were about 10 g N m^-2 h^-1 for NO₂ and 8 g N m^-2 h^-1 for NO. The NO and NO₂ flux rates were significantly reduced on plant covered soil plots. In some cases, the flux rates of both gases became negative indicating that the vegetation may act as a sink for atmospheric NO and NO₂.Application of mineral fertilizers increased the NO and NO₂ emission rates. Highest emission rates were observed for urea followed by NH₄Cl, NH₄NO₃ and NaNO₃. The fertilizer loss rates ranged from 0.1% for NaNO₃ to 5.4% for urea. Vegetation cover substantially reduced the fertilizer loss rate.The total NO _x emission from soil is estimated to be 11 Tg N yr^-1. This figure is an upper limit and includes the emission of 7 Tg N yr^-1 from natural unfertilized soils, 2 Tg N yr^-1 from fertilized soils as well as 2 Tg N yr^-1 from animal excreta. Despite its speculative character, this estimation indicates that NO _x emission by soil is important for tropospheric chemistry especially in remote areas where the NO _x production by other sources is comparatively small.