Fast space-domain evaluation of geodetic surface integrals

Geodetic surface integrals play an important role in the numerical solution of geodetic boundary-value problems. In many cases they can be evaluated using fast methods in the frequency domain (FFT). However, this is not possible in general, because the domain of integration may be non-trivial (as is the surface of the Earth), the kernel function may not be of convolution type, or the data distribution may be heterogeneous. Therefore, fast evaluation strategies are also required in the space domain. They are more difficult to design because only one property is left where a more or less fast evaluation strategy can be built upon: the potential type of the kernel function. Consequently, the idea is not to replace well-established frequency domain techniques, but to supplement them. Our approach to this problem goes in two directions: (1) we use advanced cubature methods where the integration nodes automatically densify in the vicinity of the evaluation points; (2) we use powerful computer hardware, namely MIMD computers with distributed memory. This enables us to evaluate geodetic surface integrals of any practical complexity in reasonable time and accuracy. This is shown in a numerical example. Received: 7 May 1996 / Accepted:17 March 1997     

数值模式中长波辐射参数化方案的对比试验

应用一个简化的长波辐射参数方案，与郭晓岚－钱永甫模式原有长波辐射参数化方案进行了对比，分析了它们对１００ｈＰａ，３００ｈＰａ，５００ｈＰａ，７００ｈＰａ及地面的高度场、温度场、流场及气压场模拟结果的影响，发现不仅两者结果相当一致，而且与实测结果也符合很好。     

陆面特征非均一作用参数化及其对区域气候影响的数值模拟试验

本文首先改进了陆地下垫面特征非均一性的次网格尺度参数化方法，然后利用三维地气耦合的区域气候模式，设计一系列值试验，研究了下垫面特征改变对区域气候环境变化的影响，主要分析了陆地表面特征变化对我国苏南附近地区夏季温度化的影响，结果表明采用地下垫面特征非均一作用的次网格尺度参数化方法对于改进数值模拟结果的质量有一定效果。     

地形对大地电磁测深(MTS)资料的影响

地形影响是山区ＭＴＳ工作的主要问题之一，对此做了理论分析并用均匀半空间表面二维地形模型进行了数值模拟。结果表明，地形影响不仅与ＭＴＳ资料的种类、极化模式、场源电磁波周期及测点位置有关，而且还依赖于不同的地形形态和尺度。此外，还根据悬崖模型计算，给出了关于地形影响范围的估     

水汽方程积分计算格式的数值试验

该文把Easter(1993)的水汽方程积分方案应用到球坐标的有限区模式（LAFS）中，并作了两组比较性试验。Easter的方案是把水汽的通量方程和连续方程联立求解。在计算中不会出现水汽负值。而且稳定性好。本文所做的两组试验表明。Easter的方案比LAFS的原方案所产生的计算性扩散小。因而降水预报更接近实况。但CPU时间增加了16%。     

中国西北部“4.5”沙尘暴过程中尺度低压的数值模拟

利用改进型ＰＳＵ／ＮＣＡＲ中尺度数值模式（ＭＭ４标准版），取模式水平格距４０ｋｍ，４６×６１网格，垂直方向ａ取１５层，即从地面到模式顶（１００ｈＰａ），σ＝０．０、０．１、０．２、０．３、０．４、０．５、０．６、０．７、０．７８、０．８４、０．８９、０．９３、０．９６、０．９８、０．９９、１．００，采用ＮＣＡＲ的３０＇×３０＇地形资料，以常规观测资料作为初始场，较好地模拟了此次沙尘暴过程的海平面气压的演变和分布，特别是张掖、柴达本盆地以及敦煌附近的三个中低压。同时，模拟了张掖中低压与蒙古冷高压之间的甘肃河西沙尘暴东大风。敏感性试验表明，沙尘暴中低压的形成发展主要是受于物理过程制约；沙尘暴中尺度系统的研究与暴雨中尺度系统的研究是有区别的，积云对流参数化并不是特别重要，在设计研究沙尘暴的数值模式中，应当合理地处理其他的热力、动力过程及大气外强迫源的作用。模式水平格距、地形真实程度对模拟中低压的位置、中心强度有重要贡献；下垫面变化中低压强度有一定影响。张掖热低压的形成发展主要是在有利的环境形势下，特殊地势起了重要作用，表现为直接动力强迫和间接热力强迫。     

斜向坡变形破坏机制的数值模拟研究

本文讨论了斜向坡的旋转式变形和破坏机制，用有限元法计算了此种斜坡中的剩余下滑力的分布方式，并以此来解释在滑动过程中发生一定旋转的原因；论文还运用离散元法，模拟了此种滑坡的滑动过程。     

Hazard mapping based on macroseismic data considering the influence of geological conditions

The object of this study is to consider directly the influence of regional geological conditions on the assessment of seismic hazard. It is assumed that macroseismic data at individual locations contain, in an average way, the influence of geological conditions.A Data Base referring to 199 historical (5) and instrumental (194, in the 1947–1993 period) events with macroseismic information in 1195 locations of Portugal was built. For any given seismic event, whenever macroseismic information was available at a location (town, village, etc.), an EMS-92 intensity value was estimated. To each one of those locations a geological unit, representing the most common type of soil, was assigned, based on the Geological Portuguese Map at a scale 1:500 000; the geological units were grouped into three categories: soft, intermediate and hard soils.The Data Base was used to determine the attenuation laws in terms of macroseismic intensity for the three different geological site conditions, using multiple linear regression analysis. The reasonability of the laws was tested by (i) checking residual distributions and (ii) comparing the map of isoseismals of important earthquakes with the isoseismals generated by the attenuation curves derived for each one of the three different soil classes, taking into consideration the soil class of each site. The main results of attenuation modeling are: high dispersion on macroseismic intensity data; all the models predict intensity values, for short hypocentral distances, lower than the ones observed; and for some important analyzed earthquakes and for the observed range of distances, the models confirm the expectancy that macroseismic intensity increases from hard to soft soil.The approach to obtain the hazard assessment at each location consisted in the use of the attenuation law specifically derived for the class of soil of that particular location. This method, which considers the influence of the regional geology, was illustrated with the mapping of hazard for the country for several return periods. Comparison with previous maps not taking into consideration the regional geological conditions emphasizes the importance of this new parameter. It can be concluded that (i) soil segmentation is clearly the cause for hazard increase in the region to the north of Lisbon, especially at sites with soft and intermediate soils as the ones in lower Tagus valley; the maximum increase on hazard is, in any case, less than one degree; (ii) when geological conditions are disregarded in the attenuation regression analysis, hazard pattern is similar to the one obtained for the case of hard soil everywhere.     

The mean and turbulence structure simulation of the monsoon trough boundary layer using a one-dimensional model withe-l ande-ε closures

An attempt has been made here to study the sensitivity of the mean and the turbulence structure of the monsoon trough boundary layer to the choice of the constants in the dissipation equation for two stations Delhi and Calcutta, using one-dimensional atmospheric boundary layer model withe-ε turbulence closure. An analytical discussion of the problems associated with the constants of the dissipation equation is presented. It is shown here that the choice of the constants in the dissipation equation is quite crucial and the turbulence structure is very sensitive to these constants. The modification of the dissipation equation adopted by earlier studies, that is, approximating the Tke generation (due to shear and buoyancy production) in theε-equation by max (shear production, shear + buoyancy production), can be avoided by a suitable choice of the constants suggested here. The observed turbulence structure is better simulated with these constants. The turbulence structure simulation with the constants recommended by Aupoixet al (1989) (which are interactive in time) for the monsoon region is shown to be qualitatively similar to the simulation obtained with the constants suggested here, thus implying that no universal constants exist to regulate dissipation rate. Simulations of the mean structure show little sensitivity to the type of the closure parameterization betweene-l ande-ε closures. However the turbulence structure simulation withe-ε. closure is far better compared to thee-l model simulations. The model simulations of temperature profiles compare quite well with the observations whenever the boundary layer is well mixed (neutral) or unstable. However the models are not able to simulate the nocturnal boundary layer (stable) temperature profiles. Moisture profiles are simulated reasonably better. With one-dimensional models, capturing observed wind variations is not up to the mark.     

Application of an optimum design technique for determining the coefficient of consolidation by using piezocone test data

For normally consolidated clay, several researchers have developed a number of theoretical time factors to determine the coefficient of consolidation from piezocone test results. However, depending on assumptions and analytical techniques, it could vary considerably, even for a specific degree of consolidation. In this paper a method is proposed to determine a consistent coefficient of consolidation by applying the concept of an optimum design technique over all ranges of the degree of consolidation. Initial excess pore pressure distribution is assumed to be capable of being obtained by the successive spherical cavity expansion theory. The dissipation of pore pressure is simulated by means of a two-dimensional linear-uncoupled axi-symmetric consolidation analysis. The minimization of differences between measured and predicted excess pore pressure was carried out by the BFGS unconstrained optimum design algorithm with a one-dimensional golden section search technique. By analyzing numerical examples and in-situ test results, it was found that the adopted optimum design technique gives consistent and convergent results.