3-D BAROCLINIC NUMERICAL SIMULATION OF THE SOUTH CHINA SEA I. UPPER CIRCULATION

A three-dimensional baroclinic shelf sea model was employed to simulate the seasonal characteristics of the South China Sea (SCS) upper circulation. The results showed that: in summer, an anticyclonic eddy, after its formation between the Bashi Channel and Dongsha Islands in the northeastern SCS, moves southwestward until it disperses slowly. There exists a northward western boundary current along the east shore of the Indo-China Peninsula in the western SCS and an anticyclonic gyre in the southern SCS. But at the end of summer and beginning of autumn, a weak local cyclonic eddy forms in the Nansha Trough, then grows slowly and moves westward till it becomes a cyclonic gyre in the southern SCS in autumn. At the beginning of winter, there exists a cyclonic gyre in the northern and southern SCS, and there is a southward western boundary current along the east shore of the Indo-China Peninsula. But at the end of winter, an anticyclonic eddy grows and moves toward the western boundary after forming in the Nansha Trough. The eddy‘s movement induces a new opposite sign eddy on its eastern side, while the strength of the southward western boundary current gets weakened. This phenomenon continues till spring and causes eddies in the southern SCS.     

ODP208航次——早新生代极端气候：Walvis脊断面

DP2 0 8航次将在Walvis脊东北侧水深 2 5 0 0～ 4 80 0m处钻探 5个站位 ,目标是获取早新生代沉积的一个深部断面。本区以前的钻探 (DSDP74航次 )已回收古新世 /始新世和始新世 /渐新世界线的深海软泥和白垩 ,本次钻探的目的是获取含有深部端员端位的“关键”层位以及可能的白垩纪 /第三纪 (K/T)界线的完整沉积序列。DSDP74航次站位的已有资料和新的多道地震数据可确定具有早新生代连续沉积序列的站位。在每个站位运用双 /三管先进液压活塞和双管加长岩心管取心 ,高分辨率物性测量以及综合剖面的建立 ,可确保 10 0 %的回收沉积剖面。这…     

岩浆成矿体系的热演化和剥露史的数字模拟

通过对岩浆冷却过程的数字模拟研究,揭示出岩浆在冷却成矿过程中的温度分布和变化规律及影响因素．在此基础上,进一步应用高精度的温龄计组合来限定岩浆成矿体系的热演化和剥露历史,精确地计算出岩浆的初始侵位时间和深度、矿物结晶时间、冷却速率、冷却和暴露地表时间,以及剥露和剥蚀速率等重要参数,并将模拟结果应用于斑岩铜矿床的成矿研究中．研究表明,将精确的年龄测试手段与计算机模拟技术相结合,可为定量研究岩浆矿床的热演化和剥露史、深入了解矿床的成因机制提供一种有效方法．     

Application of Accelerometer Data in Precise Orbit Determination of GRACE -A and -B

We investigate how well the GRACE satellite orbits can be determined using the onboard GPS data combined with the accelerometer data.The preprocessing of the accelerometer data and the methods and models used in the orbit determination are presented.In order to assess the orbit accuracy,a number of tests are made,including external orbit comparison,and through Satellite Laser Ranging (SLR) residuals and K-band ranging (KBR) residuals.It is shown that the standard deviations of the position differences between the so-called precise science orbits (PSO) produced by GFZ,and the single-difference (SD) and zero-difference (ZD) dynamic orbits are about 7 cm and 6 cm,respectively.The independent SLR validation indicates that the overall root-mean-squared (RMS) errors of the SD solution for days 309-329 of 2002 are about 4.93cm and 5.22cm,for GRACE-A and B respectively; theoverall RMS errors of the ZD solution are about 4.25 cm and 4.71 cm,respectively.The relative accuracy between the two GRACE satellites is validated by the KBR data to be on a level of 1.29 cm for the SD,and 1.03 cm for the ZD solution.     

A study of rainfall-runoff response in a catchment using TOPMODEL

The simplicity of Topography-based hydrological model (TOPMODEL), as a way of reflecting the topographic controls on soil water storage and runoff generation, has become more attractive and more popular for land surface process study since digital elevation models (DEMs) have become widely available. In this paper, the effect of the topography index on soil water storage distribution, which is the key to TOPMODEL, is explained. Then a simple water cycle model for estimating other components of the surface water cycle is developed, which is implemented into the TOPMODEL to integrate the water cycle of the catchment. Using the output of a DEM from 100 m×100 m resolution data and a single flow direction algorithm, the index distribution function is calculated for a catchment (around 2500 km2 )in the upper reaches of the Yangtze River under different channel initiation thresholds. Finally, the daily and monthly rainfall-runoff response from 1960 to 1987 for the catchment is simulated with the TOPMODEL coup     

流域降雨径流路径的数字模拟技术

流域内降雨径流流动过程的水文分析具有重要意义。此文基于地理信息系统的栅格系统,利用数字高程模型(DEM)提供地形特征的能力,通过一系列相连的数据集,用计算机对流域降雨径流路径进行了数字模拟,其可以实现以正确的水文顺序跟踪地表径流在流域空间范围内的流动路径,并确定地表径流在流域空间内流动路径的水文计算等级,为基于流域降雨径流关系的建模研究提供了空间分析基础。最后把该项研究成果应用于黄家二岔小流域,证实该模拟技术具有一定的可靠性和实用性。     

Consistency Problem with Tracer Advection in the Atmospheric Model GAMIL

The radon transport test, which is a widely used test case for atmospheric transport models, is carried out to evaluate the tracer advection schemes in the Grid-Point Atmospheric Model of IAP-LASG （GAMIL）. Two of the three available schemes in the model are found to be associated with significant biases in the polar regions and in the upper part of the atmosphere, which implies potentially large errors in the simulation of ozone-like tracers. Theoretical analyses show that inconsistency exists between the advection schemes and the discrete continuity equation in the dynamical core of GAMIL and consequently leads to spurious sources and sinks in the tracer transport equation. The impact of this type of inconsistency is demonstrated by idealized tests and identified as the cause of the aforementioned biases. Other potential effects of this inconsistency are also discussed. Results of this study provide some hints for choosing suitable advection schemes in the GAMIL model. At least for the polax-region-concentrated atmospheric components and the closely correlated chemical species, the Flux-Form Semi-Lagrangian advection scheme produces more reasonable simulations of the large-scale transport processes without significantly increasing the computational expense.     

Numerical Physical Mechanism and Model of Turbulent Cascades in a Barotropic Atmosphere

In a barotropic atmosphere, new Reynolds mean momentum equauons including turbulent viscosit,dispersion, and instability are used not only to derive the KdV-Burgers-Kuramoto equation but also to analyze the physical mechanism of the cascades of energy and enstrophy. It shows that it is the effects of dispersion and instability that result in the inverse cascade. Then based on the conservation laws of the energy and enstrophy, a cascade model is put forward and the processes of the cascades are described.     

Numerical Simulation on Coseismic Effect of the November 14,2001 Great Kunlun Earthquake,Northern Tibet,China

The November 14, 2001 Ms8.1 Kunlun Mountains earthquake in northern Tibet is the largest earthquake occurring on the Chinese mainland since 1950. We apply a three-dimensional （3-D） finite element numerical procedure to model the coseismic displacement and stress fields of the earthquake based on field investigations. We then further investigate the stress interaction between the Ms8.1 earthquake and the intensive aftershocks. Our primary calculation shows that the coseismic displacement field is centralized around the east Kunlun fault zone. And the attenuation of coseismic displacements on the south side of Kunlun fault zone is larger than that on the north side. The calculated coseismic stress field also indicates that the calculated maximal shear stress field is centralized around the east Kunlun fault zone; the directions of the coseismic major principal stress are opposite to that of the background crustal stress field of the Qinghai-Xizang （Tibet） Plateau. It indicates that the earthquake relaxes the crustal stress state in the Qinghai-Xizang （Tibet） Plateau. Finally, we study the stress interaction between Ms8.1 earthquake and its intensive aftershocks. The calculated Coulomb stress changes of the Ms8.1 great earthquake are in favor of triggering 4 aftershocks.