大地电磁二维对称各向异性介质的有限元数值模拟

假定垂直轴为二维对称各向异性介质主轴之一，构造走向与另一主轴方向成任一夹角，用伽勒金（Ｇａｌｅｒｋｉｎ）有限元法和矩形网格，优化地合成总体刚度矩阵，使二维各向异性介质的基本方程形成有限元代数方程组，求出各节点场值，并利用ＭＯＭ法求出辅助场，进而求出张量阻抗等响应函数．同时还对有关文献的计算模型进行了数值模拟检验     

Numerical Studies with a Regional Atmospheric Climate Model Based on Changes in the Roughness Length for Momentum and Heat Over Antarctica

A regional atmospheric climate model is used toexamine the effect of changes in the roughnesslengths of momentum (_z0m) and heat (_z0h)on the structure of the lower atmosphere and on thesurface energy fluxes over Antarctica. Fourexperiments were carried out in which _z0mand/or _z0h were altered with respect to acontrol experiment. The changes consisted of (1) alowering of _z0m from a field aggregated froma vegetation map with an orographic correction basedon the European Centre for Medium-Range WeatherForecasts _z0m field, to a constant value of10^-3 m; and (2) a lowering of _z0h from a valueequal to _z0m to a constant value of 10^-3 mor a value dependent on the wind speed via a surfacerenewal model. A reduction of _z0m results in theexpected increase in near-surface wind speed. It alsoresults in an increase in the depth of the layer in whichsouth-easterly near-surface winds prevail, and in adecrease in the strength of the large-scale flow overthe continent, in particular in summer. In theescarpment region a decrease of _z0m is foundto result in too high wind speeds. Surface temperatureson average decrease while atmospheric temperaturesincrease, resulting in an increase of near-surfacestatic stability. Changes in roughness lengths donot significantly change the temperature profiles.The surface fluxes, on average found reduced, aremodelled best by using the _z0h based on thesurface renewal method.     

宁南深层岩溶地下水系统三维数值模拟

宁南“南北古脊梁”岩溶裂隙水系统的储水空间以岩溶裂隙为主,岩溶水的空间分布明显受南北向大型断裂构造控制。根据水动力场和水化学场特征,划分了3个相对独立的岩溶水子系统。通过三维数值模拟计算,确定了马渠—洪河岩溶水子系统天然资源量和可开采资源量,为本区的地下水合理开发利用提供了科学依据。     

地基GPS气象学研究的主要问题及最新进展

目前大多数映射函数在低高度角时精度都比较差，发展适应低高度角的映射函数，提高低高度角观测值的利用率，研究大气的水平梯度是目前地基GPS气象学的主要研究任务，发展随时空变化的动态映射函数也是国际上的一个主要研究问题。由于PWV不能提供水汽分布的三维信息，求定GPS信号斜路径方向上的水汽SWV，以及利用SWV层析大气的垂直结构是目前国际上 GPS气象学的研究前沿。文中总结了目前国际上SWV的求定及层析水汽三维结构的主要方法，及存在的问题。把GPS探测的水汽同化到数值天气预报中，以提高数值预报的能力是GPS气象学研究的目的，作者介绍了目前水汽在数值预报模型中的同化情况。另外还回顾了利用GPS观测结果进行全球气候变化和大气折射环境研究的进展。     

The role of diapycnal mixing for the equilibrium response of thermohaline circulation

Using a zonally averaged, one-hemispheric numerical model of the thermohaline circulation, the dependence of the overturning strength on the surface equator-to-pole density difference is investigated. It is found that the qualitative behavior of the thermohaline circulation depends crucially on the nature of the small-scale vertical mixing in the interior of the ocean. Two different representations of this process are considered: constant vertical diffusivity and the case where the rate of mixing energy supply is taken to be a fixed quantity, implying that the vertical diffusivity decreases with increasing stability of the water column. When the stability-dependent diffusivity parameterization is applied, a weaker density difference is associated with a stronger circulation, contrary to the results for a fixed diffusivity. A counterintuitive consequence of the stability-dependent mixing is that the poleward atmospheric freshwater flux, which acts to reduce the thermally imposed density contrast, strengthens the thermally dominated circulation and its attendant poleward heat transport. However, for a critical value of the freshwater forcing, the thermally dominated branch of steady states becomes unstable, and is succeeded by strongly time-dependent states that oscillate between phases of forward and partly reversed circulation. When a constant vertical diffusivity is employed, on the other hand, the thermally dominated circulation is replaced by a steady salinity-dominated state with reversed flow. Thus in this model, the features of the vertical mixing are essential for the steady-state response to freshwater forcing as well as for the character of flow that is attained when the thermally dominated circulation becomes unstable.Responsible Editor: Jin-Song von Storch     

Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.     

Tectonic evolution of the Lachlan Fold Belt, southeastern Australia: constraints from coupled numerical models of crustal deformation and surface erosion driven by subduction of the underlying mantle

We have used a coupled thermo-mechanical finite-element (FE) model of crustal deformation driven by mantle/oceanic subduction to demonstrate that the tectonic evolution of the Lachlan Fold Belt (LFB) during the Mid-Palaeozoic (Late Ordovician to Early Carboniferous) can be linked to continuous subduction along a single subduction zone. This contrasts with most models proposed to date which assume that separate subduction zones were active beneath the western, central and eastern sections of the Lachlan Orogen. We demonstrate how the existing data on the structural, volcanic and erosional evolution of the Lachlan Fold Belt can be accounted for by our model. We focus particularly on the timing of fault movement in the various sectors of the orogen. We demonstrate that the presence of the weak basal decollement on which most of the Lachlan Fold Belt is constructed effectively decouples crustal structures from those in the underlying mantle. The patterns of faulting in the upper crust appears therefore to be controlled by lateral strength contrasts inherited from previous orogenic events rather than the location of one or several subduction zones. The model also predicts that the uplift and deep exhumation of the Wagga-Omeo Metamorphic Belt (WOMB) is associated with the advection of this terrane above the subduction point and is the only tectonic event that gives us direct constraints on the location of the subduction zone. We also discuss the implications of our model for the nature of the basement underlying the present-day orogen.     

MODELING OF THE DIVERSION CHANNEL CLOSURE FOR THE THIRD STAGE OF THE THREE GORGES PROJECT

The closure of third stage diversion channel for the Three Gorges Dam is characterized by large closure discharge, large drop in water surface, and high gap velocity. 1D and 2D flow mathematical models were used in this paper to simulate the hydraulic conditions during the closure. The variation of discharge in the diversion channel and the drops in water surface shared by the upstream and downstream cofferdams were computed using the ID model, and the detailed hydraulic patterns in the diversion channel were simulated using the 2D model. The computed results indicate that the designed closure scheme for discharge of 9,010 m^3/s was feasible for construction, while the designed closure scheme for discharge of 12,200 m^3/s was inapplicable.     

Bearing capacity of rock over mined cavities in Nottingham

A significant geohazard is created in Nottingham, UK, by hundreds of man-made caves cut in the weak sandstone beneath the city centre. Stability of the caves has been assessed by a single full-scale loading test, by numerical modelling with FLAC and by physical modelling in plaster. For typical caves 4 m wide, bearing capacity of the rock roof rises from 2 MPa where it is 1 m thick to 8 MPa where 3 m thick. Stability decreases over wider caves and where the loading pad edge is over the edge of the cave. Numerical modelling of a very wide cave revealed the failure mechanisms and also showed that an internal support wall increased roof bearing capacity by 50%. Local building regulations that require 3-5 m of rock cover over the sandstone caves appear to be conservative. In stronger rocks, including karstic limestone, a guideline that cover thickness exceeds 70% of the cave width appears to be appropriate.     

Proposed support design, Kaletepe tunnel, Turkey

In this paper, preliminary support design of Kaletepe tunnel, located on Bilecik-Istanbul highway, Turkey, was analyzed by empirical and numerical methods. The rock mass rating (RMR) and rock mass quality (Q) systems were employed for empirical rock mass quality determination. Numerical analysis for the stress–strain distribution of the tunnel excavation and support systems was also carried out. The applied support performance was investigated at different sections of the tunnel route. It was seen that empirical and numerical approaches showed similar results. This indicates that when the empirical method is supported by numerical method, the preliminary support design will be more reliable.