Numerical simulations of katabatic jumps in coats land,Antartica

A non-hydrostatic numerical model, the Regional Atmospheric Modeling System (RAMS), has been used to investigate the development of katabatic jumps in Coats Land, Antarctica. In the control run with a 5 m s^-1downslope directed initial wind, a katabatic jump develops near the foot of the idealized slope. The jump is manifested as a rapid deceleration of the downslope flow and a change from supercritical to subcritical flow, in a hydraulic sense, i.e., the Froude number (Fr) of the flow changes from Fr > 1 to Fr> 1. Results from sensitivity experiments show that an increase in the upstream flow rate strengthens the jump, while an increase in the downstream inversion-layer depth results in a retreat of the jump. Hydraulic theory and Bernoulli's theorem have been used to explain the surface pressure change across the jump. It is found that hydraulic theory always underestimates the surface pressure change, while Bernoulli's theorem provides a satisfactory estimation. An analysis of the downs balance for the katabatic jump indicates that the important forces are those related to the pressure gradient, advection and, to a lesser extent, the turbulent momentum divergence. The development of katabatic jumps can be divided into two phases. In phase I, the t gradient force is nearly balanced by advection, while in phase II, the pressure gradient force is counterbalanced by turbulent momentum divergence. The upslope pressure gradient force associated with a pool of cold air over the ice shelf facilitates the formation of the katabatic jump.     

各向异性半空间上方垂直磁偶极子电磁场的球面波展开式

文中推导出单轴各向异性媒质半空间上方垂直磁偶极源电磁场的精确解析表达式. 应用圆柱波函数的球面波展开式和超几何函数理论，场分量中的索末菲尔德型积分被表示成快速、绝对收敛的球面波函数系展开式；展开系数是以物性参数为复宗量的勒让德多项式. 该展开式数学物理意义明显，并且不受场点和源点的位置、媒质的物性参数和频率等条件的限制. 利用本文的结果可十分方便地计算和分析任意场点处的电磁场分布.     

图们江流域中更新世火山泥石流分布及其特征

图们江流域中更新世火山泥石流仅分布在残留的Ⅲ阶地上，碎屑物的粒级范围很广，具正粒序层理，分选差，是在冰川作用或冰雪融化作用下快速搬运、沉积形成的，具有浊积岩的沉积模式。     

计算机图形软件编制地图的探讨

阐述了利用Coreldraw8．0图形软件来完成《连云港市市区图》的设计和具体编制的方法，利用Coreldraw8．0图形软件的特性来提高制图效率、成图质量。结合在地图编制工作中遇到的问题，对计算机制图和传统地图理论的结合、发展提出一些看法。     

High-resolution structures of the Landers fault zone inferred from aftershock waveform data

High-frequency body waves recorded by a temporary seismic array across the surface rupture trace of the 1992 Landers, California, earthquake were used to determine fault-zone structures down to the seismogenic depth. We first developed a technique to use generalized ray theory to compute synthetic seismograms for arbitrarily oriented tabular low-velocity fault-zone models. We then generated synthetic waveform record sections of a linear array across a vertical fault zone. They show that both arrival times and waveforms of P and S waves vary systematically across the fault due to transmissions and reflections from boundaries of the low-velocity fault zone. The waveform characteristics and arrival-time patterns in the record sections allow us to locate the boundaries of the fault zone and to determine its P - and S -wave velocities independently as well as its depth extent. Therefore, the trade-off between the fault-zone width and velocities can be avoided. Applying the method to the Landers waveform data reveals a low-velocity zone with a width of 270–360 m and a 35–60 per cent reduction in P and S velocities relative to the host rock. The analysis suggests that the low-velocity zone extends to a depth of ∼7 km. The western boundary of the low-velocity zone coincides with the observed main surface rupture trace.     

直线补给边界非稳定井流数值模拟研究

针对解析方法求解半无限含水层井流问题,对直线补给边界附近的井流采用数值方法进行模拟计算,并将两种方法进行比较,数值计算结果表明,存在补给边界时,抽水一定时间以后降深能达到稳定,但仅仅代表观测点处降深达到稳定,而整个渗流场并未达到稳定状态.     

SEAWAT 2000: modelling unstable flow and sensitivity to discretization levels and numerical schemes

A systematic analysis shows how results from the finite difference code SEAWAT are sensitive to choice of grid dimension, time step, and numerical scheme for unstable flow problems. Guidelines to assist in selecting appropriate combinations of these factors are suggested. While the SEAWAT code has been tested for a wide range of problems, the sensitivity of results to spatial and temporal discretization levels and numerical schemes has not been studied in detail for unstable flow problems. Here, the Elder-Voss-Souza benchmark problem has been used to systematically explore the sensitivity of SEAWAT output to spatio-temporal resolution and numerical solver choice. A grid size of 0.38 and 0.60% of the total domain length and depth respectively is found to be fine enough to deliver results with acceptable accuracy for most of the numerical schemes when Courant number (Cr) is 0.1. All numerical solvers produced similar results for extremely fine meshes; however, some schemes converged faster than others. For instance, the 3rd-order total variation-diminishing method (TVD3) scheme converged at a much coarser mesh than the standard finite difference methods (SFDM) upstream weighting (UW) scheme. The sensitivity of the results to Cr number depends on the numerical scheme as expected.     

Post-drilling destabilization of temperature profile in borehole Yaxcopoil-1, Mexico

As part of the International Continental Scientific Drilling Program (ICDP), the 1.5-km-deep borehole Yaxcopoil-1, located in the Chixculub meteor impact structure in Mexico, has undergone further study after drilling operations ceased. Temperature logs were repeated ten times at intervals 0.3–0.8, 15, 24 and 34 months after borehole shut-in. The logs bear a distinct signature of transient heat transfer by groundwater flow manifested by a gradual distortion of the linear temperature profile when a cold wave of 0.8–1.6°C amplitude was detected propagating downward from 145 to 312 m at a rate of 4–6 m/month. To understand the nature of this moving anomaly, a 20-day monitoring of the cold wave was carried out at a depth of 307 m that showed further cooling of 0.6°C during the first 16 days of the passage followed by temperature stabilisation. As an explanation of this unusual phenomenon, a theory is proposed, whereby the drilling mud has accumulated within the overlying and cooler highly porous and permeable karstic rocks during the drilling and migrates downward. The observed migration rate suggests a permeability higher than 10^?11 m². This indicates a high vulnerability to contamination of the only freshwater aquifer in the Yucatan region.     

Numerical simulation of groundwater flowing to horizontal seepage wells under a river

A horizontal seepage well, consisting of an interconnected vertical well, galleries, chambers and small-diameter radiating bores, is used to acquire relatively clean water that has been filtered through natural alluvial deposits in a riverbed. It has wide application, especially in arid and semi-arid areas. The lack of calculation formulae or models for horizontal seepage wells, up until now, has resulted in several false applications. Based on the analysis of groundwater flow characteristics, it has been concluded that several flow regimes coexist and hydraulic head loss exists in the horizontal seepage well. To avoid the difficulty of confirming the flux or head distribution in such a complex system, the model boundary of the whole horizontal seepage well has been moved to that of just the vertical well, and the well-aquifer system was treated as a heterogeneous medium, where the horizontal seepage well itself is a highly permeability medium. A mathematical model has been developed, based on the coupled seepage-pipe flow, by the introduction of equivalent hydraulic conductivity according to different flow regimes. Then a three-dimensional finite difference numerical model, based on the mathematical model, was developed and applied to a horizontal seepage well in China. The numerical model verified the groundwater flow characteristics of the horizontal seepage well. An erratum to this article can be found at