Integral inversion of SST data of type GRACE

Satellite missions CHAMP and GRACE dedicated to global mapping of the Earth’s gravity field yield accurate satellite-to-satellite tracking (SST) data used for recovery of global geopotential models usually in a form of a finite set of Stokes’s coefficients. The US-German Gravity Recovery And Climate Experiment (GRACE) yields SST data in both the high-low and low-low mode. Observed satellite positions and changes in the intersatellite range can be inverted through the Newtonian equation of motion into values of the unknown geopotential. The geopotential is usually approximated in observation equations by a truncated harmonic series with unknown coefficients. An alternative approach based on integral inversion of the SST data of type GRACE into discrete values of the geopotential at a geocentric sphere is discussed in this article. In this approach, observation equations have a form of Green’s surface integrals with scalar-valued integral kernels. Despite their higher complexity, the kernel functions exhibit features typical for other integral kernels used in geodesy for inversion of gravity field data. The two approaches are discussed and compared based on their relative advantages and intended applications. The combination of heterogeneous gravity data through integral equations is also outlined in the article. panovak@kma.zcu.cz     

中国新一代高精度、高分辨率大地水准面的研究和实施

采用移去恢复技术，利用我国高分辨率DTM和重力资料推算我国大陆重力大地水准面;然后再和我国GPS水准所构成的高程异常控制网拟合,推算了具有dm级精度、15′×15′分辨率的我国大陆大地水准面。利用全国地壳运动监测网络的80余个高精度GPS水准点进行外部检核，检核结果证实和原设计精度完全一致，即该大陆大地水准面的绝对精度,在东经102°以东高于±0.3m，在东经102°以西、北纬36°以北为±0.4m，36°以南为±(0.4～0.6)m。利用卫星测高数据计算垂线偏差，反解我国海域大地水准面。为了检核，由测高垂线偏差反演为重力异常，与海上万余点船测重力值进行了外部检核；同时用上述反演的重力异常推算大地水准面，与直接解得的相应结果进行比较作为内部检核。由重力和GPS水准数据推算的上述大陆大地水准面，和主要由卫星测高数据确定的海洋大地水准面，二者之间一般都存在以系统误差为主的拼接差。顾及这一现象并结合我国在陆海大地水准面拼接区重力资料稀疏的实际，研究提出了扩展拼接技术，即在沿海选取部分陆海毗邻的局部地区，在这局部地区内，陆地用实测平均重力格网数据，海洋用测高平均重力格网数据，统一推算这陆海局部重力大地水准面。然后利用这一局部大地水准面的陆地部分和已经用GPS水准校正的陆地大地水准面进行拟合。最后用拟合参数校正中国全部海域的测高重力大地水准面，从而保持陆地部分大地水准面不变，最大限度地削弱拼接点和测高海洋大地水准面的系统误差。     

高程现代化问题

海拔高在工程应用中有不同可代替的作用。测量海拔高的经典方法是几何水准，它的主要缺点是劳动强度大、效率低、花费大、实时性差。GPS海拔高测量在很大程度上克服了这些缺点，而且可以满足许多应用的高程需求。GPS海拔高测量代替水准测量，以GPS海拔高测量为主的“高程现代化”代表大地测量的一个发展方向。     

球面小波分析技术及在全球重力模型中的应用

在介绍球面小波理论的基础上,推导和比较了几种球面小波,分析了最新地球重力模型-EGM96,以此为依据,把球面小波多分辨分析用于计算全球自由空气异常及重力大地水准面,并对处理结果做出解释。     

Characterisation of tobelite (NH4)Al2[AlSi3O10](OH)2 and ND4-tobelite (ND4)Al2[AlSi3O10](OD)2 using IR spectroscopy and Rietveld refinement of XRD spectra

Tobelite (NH₄) Al₂ [AlSi₃O₁₀] (OH)₂, the ammonium analogue of muscovite, and its deuterated form ND₄-tobelite (ND₄) Al₂ [AlSi₃O₁₀] (OD)₂ have been synthesised at 600?°C and 200 and 500 Mpa using a well homogenised, stoichiometric SiO₂-Al₂O₃ oxide mix with Al₂O₃ in excess of 5 mol% and a 25% NH₃ solution whose relative abundance was such that the amount of NH₄ ⁺ stoichiometrically available was in excess of 50%. Characterisation of both tobelite and ND₄-tobelite using IR-spectroscopy, Rietveld refinement of X-ray powder diffraction data, and electron microprobe analysis indicate that, similar to K⁺ in muscovite, the NH₄ ⁺ or ND₄ ⁺ molecule occupies the interlayer site. IR absorption bands caused by NH₄ ⁺ and ND₄ ⁺ can be explained, to a very good approximation, on the basis of T_d symmetry. Nevertheless, substantial line broadening and the occurrence of shoulders indicate a deviation from ideal T_d symmetry. However, even at 77?K, no discrete splitting of the degenerate states could be confirmed. The OH stretching frequencies observed for synthetic tobelite are quite similar to those for muscovite, indicating that the replacement of K⁺ by NH₄ ⁺ has no effect. The low FWHH of the OH bands indicate that the hydroxyl groups are well ordered within the structure. Rietveld refinement of tobelite and ND₄-tobelite indicates that all samples synthesised consist of the 3 different mica polytypes which are typical of muscovite – namely 1M (C2/m), 2M ₁ (C2/c) and 2M ₂ (C2/c). Tobelite and ND₄-tobelite synthesised at 500 Mpa principally contain the 1M polytype, whereas the principle polytype for ND₄-tobelite synthesised at 200 Mpa, is 2M ₂. Rietveld refinement of X-ray diffraction spectra for tobelite synthesised at 200 Mpa was problematic due to the very broad FWHH of the X-ray peaks indicating poor crystallinity. In comparision to synthetic muscovite, the cell dimensions observed for tobelite and its deuterated analogue are quite similar except for the lattice constant c. Due to the larger radius of NH₄ ⁺ or ND₄ ⁺ compared to K⁺ cation, the c-direction is expanded form 10.275 Å in muscovite to approximately 10.540 Å in tobelite and ND₄-tobelite.     

卫星测高研究应用新进展

卫星测高是利用空间技术进行海洋观测的崭新技术手段。GEOSAT卫星的升空又使得这一研究向前大大地推进。本文简要介绍卫星测高在确定大地水准面精度、研究岩石圈特性、探测海底深部构造、研究海洋动力环境以及其经济效益和社会效益的最新进展。     

中国东南陆缘大地水准面的确定

利用我国东南地区１１８２２个重力数据、３０″×３０″数值地形模型和ＷＤＭ９４重力位模型,以Ｓｔｏｋｅｓ理论为基础,采用低通滤波方法处理远区影响,计算了中国东南陆缘２．５′×２．５′格网重力大地水准面,其结果与福建、山东的ＧＰＳ水准相比较,精度分别为±０．１２４ｍ和±０．０６２ｍ。     

Crystal chemistry of tremolite–tschermakite solid solutions

Tremolite–tschermakite solid solutions have been synthesized between 700 and 850 °C and 200 and 2000 MPa. The starting materials were oxide–hydroxide mixtures and an additional 0.1–1.8 molal CaBr₂ solution. The run products were characterized using SEM, HRTEM, EMP, XRD and FTIR. The synthesized Al tremolites formed needles and lath-shaped crystals of up to 300 × 20 μm. HRTEM investigations showed that the majority of the amphiboles were well ordered. The EMP analysis revealed that the Al tremolites were solid solutions in the ternary tremolite–tschermakite–cummingtonite. The highest observed Al content was close to the composition of magnesiohornblende (X_ts=0.54). Different cummingtonite concentrations (X_cum=0.00–0.18) were observed, which generally increased with Al content. Rietveld refinements of the lattice constants showed a linear decrease of the cell parameters a and b with increasing Al content, whereas c and β increased. Small deviations from the linear behaviour were caused by variable amounts of the cummingtonite component. For pure tschermakite lattice parameters of a=9.7438(11) Å, b=17.936(14) Å, c=5.2995(3) Å, β=105.68(9)° and V=891.7 ± 1.4 ų were extrapolated by least-squares regression. Using the a and β lattice parameters for tremolite, tschermakite and cummingtonite, it was possible to derive amphibole compositions using powder XRD. IR spectra of the Al tremolites showed a total of 12 individual bands. The FWHMs of all bands increased with increasing Al content. According to their FWHMs, these bands were grouped into three band systems at 3664–3676 cm^?1 (I), 3633–3664 cm^?1 (II) and 3526–3633 cm^?1 (III). Assuming ^[6]Al substitution at M2 and/or M3 and ^[4]Al at T1, three principal different configurational groups could be assigned as local environments for the proton. I: only Si⁴⁺ at T1 and one or two Al³⁺ at M2 and/or M3_far, II: one Al³⁺ at T1 and one to three Al³⁺ at M2 and/or at M3_far, III: either Al³⁺ on M3_near and/or two Al³⁺ on T1 and additional one to four Al³⁺ at M2. It is assumed that these three configurational groups correspond to the three groups of observed bands. This was quantitativly supported by Monte-Carlo simulations. A model with random distribution at M2 and M3 including Al avoidance at tetrahedral and octahedral sites yielded the best agreement with the spectroscopical results.     

Improving the efficiency of ant algorithms using adaptive refinement: Application to storm water network design

The ant algorithm is a new evolutionary optimization method proposed for the solution of discrete combinatorial optimization problems. Many engineering optimization problems involve decision variables of continuous nature. Application of the ant algorithm to the optimization of these continuous problems requires discretization of the continuous search space, thereby reducing the underlying continuous problem to a discrete optimization problem. The level of discretization of the continuous search space, however, could present some problems. Generally, coarse discretization of the continuous design variables could adversely affect the quality of the final solution while finer discretization would enlarge the scale of the problem leading to higher computation cost and, occasionally, to low quality solutions. An adaptive refinement procedure is introduced in this paper as a remedy for the problem just outlined. The method is based on the idea of limiting the originally wide search space to a smaller one once a locally converged solution is obtained. The smaller search space is designed to contain the locally optimum solution at its center. The resulting search space is discretized and a completely new search is conducted to find a better solution. The procedure is continued until no improvement can be made by further refinement. The method is applied to a benchmark problem in storm water network design discipline and the results are compared with those of existing methods. The method is shown to be very effective and efficient regarding the optimality of the solution, and the convergence characteristics of the resulting ant algorithm. Furthermore, the method proves itself capable of finding an optimal, or near-optimal solution, independent of the discretization level and the size of the colony used.