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     

An analytical solution to multi-component NAPL dissolution equations

This note presents a novel method for determining the changing composition of a multi-component NAPL body dissolving into moving groundwater, and the consequent changes in the aqueous phase solute concentrations in the surrounding pore water. A canonical system of coupled non-linear governing equations is derived which is suitable for representation of both pooled and residual configurations, and this is solved. Whereas previous authors have handled such problems numerically, it is shown that these governing equations succumb to analytical solution. By a suitable substitution, the equations become decoupled, and the problem collapses to a single first-order equation. The final result is expressed implicitly, with time as a function of the number of moles of the least soluble component, m₁. The number of moles of each other component is expressed explicitly in terms of m₁. It is shown that the time-m₁ relationship has a well behaved inverse. An example is given in which the analytic solution is verified against traditional finite difference analysis, and its computational efficiency is shown.     

Shear waves at a corrugated interface between anisotropic elastic and visco-elastic solid half-spaces

Rayleigh’s method of approximation is employed to find out the reflection and transmission coefficients due to an incident plane SH wave at a corrugated interface between a laterally and vertically inhomogeneous anisotropic elastic solid half-space and a laterally and vertically inhomogeneous isotropic visco-elastic solid half-space. The lateral and vertical inhomogeneities are described by the exponential variations of elastic parameters. The formulae of reflection and transmission coefficients are derived in closed form for the first-order approximation of the corrugation. The effects of the corrugation of the interface, the inhomogeneity, the anisotropy, the visco-elasticity and the frequency of the incident wave on these coefficients are studied analytically and numerically for a specific model containing a periodic interface. The results of earlier workers have been reduced as particular cases from the present formulation.     

The November 2002 eruption at Piton de la Fournaise volcano,La Réunion Island: ground deformation,seismicity, and pit crater collapse

An eruption on the eastern flank of Piton de la Fournaise volcano started on 16 November, 2002 after 10 months of quiescence. After a relatively constant level of activity during the first 13 days of the eruption, lava discharge, volcanic tremor and seismicity increased from 29 November to 3 December. Lava effusion suddenly ceased on 3 December while shallow earthquakes beneath the Dolomieu summit crater were still recorded at a rate of about one per minute. This unusual activity continued and increased in intensity over the next three weeks, ending with the formation of a pit crater within Dolomieu. Based on ground deformation, measured by rapid-static and continuous GPS and an extensometer, seismic data, and lava effusion patterns, the eruptive period is divided into five stages: 1) slow summit inflation and sporadic seismicity; 2) rapid summit inflation and a short seismic crisis; 3) rapid flank inflation, onset of summit deflation, sporadic seismicity, accompanied by stable effusion; 4) flank inflation, coupled with summit deflation, intense seismicity, and increased lava effusion; and finally 5) little deflation, intense shallow seismicity, and the end of lava effusion. We propose a model in which the pre-intrusive inflation of Stage 1 in the months preceding the eruption was caused by a magma body located near sea level. The magma reservoir was the source of an intrusion rising under the summit during Stage 2. In Stage 3, the magma ponded at a shallow level in the edifice while the lateral injection of a radial dike reached the surface on the eastern flank of the basaltic volcano, causing lava effusion. Pressure decrease in the magmatic plumbing system followed, resulting in upward migration of a collapse front, forming a subterranean column of debris by faulting and stoping. This caused intense shallow seismicity, increase in discharge of lava and volcanic tremor at the lateral vent in Stage 4 and, eventually the formation of a pit crater in Stage 5.     

长江口九段沙下段冲淤演变水动力机制分析

概述了长江口深水航道整治一期工程及工程对九段沙下段地形冲淤的影响.根据已建立的σ坐标系下三维非线性水流数学模型,用新测水文资料对模型进行验证.利用所建立的模型,结合由现场测量得到一期工程前后九段沙下段近期地形演变,数值模拟了北槽区域流通量、底层欧拉余流和北槽中下段平面水流特征,以综合分析北槽水域水动力变化对九段沙下段地形冲淤变化的影响.     

Precise determination of focal parameters for July 20, 1995 M_L=4.1 earthquake sequence in the Huailai basin

IntroductionSince the late 1970s, the quickly developed global digital seismograph network has been providing high quality recordings of large earthquakes in global scale, based on which digital seismology has made great progress. Compared with large earthquakes, moderate and small sized shocks have more frequent occurrence, and comprise clues to geological tectonics and tectonic stress field in a region. Preceding and following a large earthquake, usually occur numbers of small events that im…     

地下工程震波技术与应用

激发的震波在工程岩土体中传播，检测其波动力学参数，来测试岩土体的工程力学性质。原位震波波速测试技术可以对井巷的岩体模量，支护选型，围岩性质进行检测；利用地震勘探，面波勘探方法可以对矿层厚度，超前探测等方面进行有效检测；震波CT、导波勘探等探测方法对于采面内的地质异常探测有较好的针对性；通过声发射检测技术可以解决地下工程安全问题。实际应用说明，震波探测技术从地下工程的开拓掘进，采场形成到采后管理等过程中有着广泛的应用前景。一些技术已作为地下工程中的日常工作加以应用，但进一步的生产应用还需深入研究。     

浙江雷甸淡水无光珠中球文石的首次确认

球文石（Vaterite)是CaCO3很不稳定的同质多象变体，在自然界里很少出现，通过光学显微镜在淡水无光珠中发现球文石尚属首次，本文研究了淡水无光珠内球文石的显微特征、X射线粉晶衍射特征及显微红外光谱特征，研究结果表明淡水无光珍珠内球文石粉晶衍射由强至弱五个强峰的d值为0.3266,0.2723,0.3553,0.1811,0.2057nm；其红外光谱特征波数为1489，1450,1420,1087,1050,876,830,762,743cm^-1。球文石的存在是影响淡水珍珠品质的主要内在因素之一。     

Earth’s density flattening and hypothesis of latitudinal normal density

In this paper, the definition of latitudinal density and density flattening of the level ellipsoid is given, and integral formulas of latitudinal density for pole gravity and equator gravity are derived. According to the pole gravity condition and equator gravity condition for the level ellipsoid, latitudinal density distribution function of the level ellipsoid is obtained. It is proved mathematically that latitudinal density of the earth’s equator is larger than that of the pole, the earth’s density flattening calculated preliminarily is 1/322, and hypothesis of the earth’s latitudinal normal density is further proposed, so that theoretical preparation for studying the forming cause of the earth gravity in problems such as continent drift, mantle convection, and submarine extension is made well.     

The Green''''s function of the harmonic horizontal force applied at the interior of the saturated half-space soil

IntroductionThe wave propagation problems in saturated soil are very important for the civil engineering, geophysics and seismology. Biot (1956,1962) established the theory of wave propagation in saturated soil firstly, and hereafter many researchers have used Biot theory to study wave propagation problems in saturated soil. By using integral transform and potential function method, Philippacopoulos (1988) studied the Lamb(s problem of a vertical point force applied to the surface of saturate…