黑龙江省未来30、50年气候预测

通过CCCma、CCSR、CSIRO、Gfdl和Hadley气候模式对黑龙江省其中包括齐齐哈尔、佳木斯、哈尔滨、牡丹江等4区未来50年(2005～2050年)，在GG、GS情景下数值模拟。结果表明(采用GS结果)，未来2030、2050年气温均有较大增高。其中2030年年平均气温可增高1．94℃；春季提高2．06℃：；夏季提高1．29℃；秋季提高1．79℃；冬季提高2．66℃，2050年又继续增加．年平均气温将提高2．42℃；春季提高2．13℃；夏季提高1．68℃；秋季提高2．56℃；冬季提高3．21℃。冬季是四季中增幅最大的季节，其次秋季、春季和夏季。如果按GG情景下，未来气温还要高出1℃。增温中心在西部齐齐哈尔，增温较小为牡丹江。从哈尔滨年蒸发量来看，2030年可增加11％，2050年可增加13％。     

中国西部退耕还林(草)和沙漠化土地绿化的区域性气候效应

中国正在实施的西部大开发必将导致较大规模的土地利用和地表覆盖(LULC)变化,其中最重要的两类LULC变化是退耕还林(草)和沙漠化土地绿化工程。利用最新的全球地表特征资料数据库(GLCCD)识别以上两项生态建设工程的实施规模与范围,并以此更新了目前比较流行的区域气候模式的下边界条件,利用修正了的区域气候模式研究上述大规模LULC变化可能对中国区域性气候造成的影响。使用的卫星资源数据库是由美国地质测绘局(USGS)和Nebraska-Lincoln大学牵头,根据NOAA卫星的AVHRR资料、USGS掌握的高分辨率数字化地图,以及生态区划和植被分布等资料组建的。     

偏心结构非线性地震反应分析的一种简化方法

本文建立了基于结构状态方程和偏心结构体系的层剪力-扭矩等效屈服EYST面概念的求解平扭耦联非线性地震反应的分离-予估-修正递推算法，编制了相应的计算程序。通过与三个单层偏心结构模型的破坏性地震模拟试验结果的对比，验证了小变形阶段平-扭耦联地震反应算法的实用性和可靠性。算例表明，该程序数据处理量小，运算快捷，便于应用。     

团树断层场地地震反应的数值模拟

团树断层场地在多次地震中出现了烈度异常。为了弄清烈度异常的原因，应用显式有限元方法，数值模拟了在脉冲地震波斜入射情形下断层场地的地震反应。结果表明：在基岩山脚附近土层场地地震动放大比基岩山区平均高4倍、比远离基岩山脚土层场地平均高2倍，与历史震害分布有较好的对应。     

Turbulence Statistics Inside and Over Forest: Influence on Footprint Prediction

Observations of wind statistics within and above a Scots pine forest are comparedwith those predicted from an analytical second-order closure model. The roughnesssublayer (RSL) effects, and the influence of stability on similarity functions, arestudied using observations. The commonly accepted forms of similarity functionsdescribe the influence of diabatic effects above the RSL well. According to earlierstudies they are expected also to apply within the RSL. As an exception, the averagewind speed normalised with friction velocity was found to be invariant with stabilityclose to the canopy top under unstable conditions. Lagrangian stochastic trajectorysimulations were used to evaluate the influence of canopy turbulence profiles onfootprint prediction. The main uncertainty was found to arise from parameterisationof the random forcing term in the Lagrangian velocity equation. The influence ofdiabatic conditions was studied, and it was found that thermal stability affectssignificantly the footprint function above the forest canopy, but significantuncertainty exists because of uncertainties in the formulation of stability functions.     

磁场重联中离子轨道的混合模拟研究

在使用2.5维混合模拟方法研究了Petschek模型稳态驱动磁场 重联演化的基础上， 本文考察了计算域内各典型区域中粒子分布函数的变化，描绘了重联区不同位置几种类型的 非Maxwell分布函数. 结果表明，磁场重联会将重联区少部分粒子加速到很高的能量，不同 加速程度的粒子将形成球壳状的速度分布. 粒子的轨道特征表明，在重联区中出流的粒子， 有一部分被磁镜捕获，其回旋半径大于重联区宽度，并构成整个流体速度的低速部分. 另外 ，在X中性点附近进入重联区的粒子沿磁力线向出流区以三种形式漂移，分别为：沿磁力线 逃逸、捕获在磁镜中随流体运动、横越磁力线漂移，其比例分别约为70%，20%和10%。     

黄海、东海海域出海气旋发展过程中尺度数值模拟

利用MM5中尺度模式对1999年6月两个出海气旋发展过程进行 数值模拟. 数值模拟的气旋出海后移动路径与实际情况基本一致. 在数值模拟基础上重点讨 论了出海气旋发展过程潜热通量和感热通量的分布及其演变情况. 气旋出海后在气旋中心区 南方和东方存在负潜热通量和感热通量区. 出海气旋的东移和发展，其前方强大正热通量区 的存在可能是重要原因之一.     

1966年以来华北地区一系列七级大震破裂过程的数值模拟

用非连续变形分析方法(DDA+FEM)数值模拟，在华北地区各地块相互制约的块体系统环境中，地块边界断层上发生1966年邢台地震、 1969年渤海地震、 1975年海城地震、 1976年唐山地震等不同类型破裂模式大震的破裂过程. 数值模拟结果给出每个大震释放的主应力场，最大剪应力变化等值线图，地震前后位移变化矢量图，发震断层滑移随时间变化以及走滑错距和应力降等震源参数. 这些结果与地震的震源机制，用地震波资料研究得到的震源参数,宏观等震线,地表观测的水平位移矢量图基本一致. 其中1969年渤海地震正交破裂模式的结果与宏观等震线及小震分布图像更接近. 1976年唐山地震复杂震源模式与该震早期余震分布图像更相符. 表明用DDA+FEM方法进行数值模拟研究，能较好地模拟地震破裂过程.     

热层风场的理论模拟与分析

通过求解中性大气Navier Stokes动量方程建立了一个时变的三维风场理论模式，利用目前新版的中性大气模式NRLMSISE 00及国际电离层参考模式IRI2000作为输入参数给出热层风场. 基于该模式，计算得到中等太阳活动年磁静日风场的变化形态及其受电场和离子曳力的影响. 同时，将Navier Stokes动量方程作不同形式的简化，并利用简化模式与本文的模式计算结果的对比，分析中性大气Navier Stokes动量方程中黏性项以及非线性项(U·Δ)U的作用. 结果表明，本文所建立伪三维风场模式给出的结果更为合理，而简化模式在某些地区尤其在低纬和赤道区不适用，黏性项及非线性项的作用不可忽略. 本文所建立的风场模式将对研究电离层动力学过程、电离层与热层的耦合过程以及空间天气学研究都有着重要意义.     

Numerical Simulation of Drag Partition over Rough Surfaces

We present a numerical simulation of drag partition over rough surfaces. A computational fluid dynamics model is applied with high resolution to simulatingturbulent flows over arrays of roughness elements positioned on asmooth surface. The skin drag on the surface and the pressure drag on the roughnesselements are computed. The simulated drag partition compares well with wind-tunnelmeasurements and theoretical estimates for similar rough surfaces. This confirms that the computational approach offers an alternative to wind-tunnel and field experiments in studying drag and drag partition. The model is then applied to studying drag partition over rough surfaces with various roughness configurations. It is shown that drag partition depends not only on the magnitude of the roughness frontal area but also on the sizes and arrangement of roughness elements, because (1) the pressure drag coefficient is sensitive to roughness-element dimensions and (2) the arragement of roughness elements lead to different interferences of turbulent wakes. The impact ofthe latter factor is not insignificant.