长江中下游降水异常特征及其与全国降水和气温异常的关系

利用１９５１－１９９８年中国１６０站逐月降水和气温资料,分析了长江中下游降水异常的基本特征,并着重分析了它与全国其他地区降水和气温异常的关系。结果表明：长江中下游降水异常有明显的季节,年际和年代际变化特征、６－７月降水异常程度大,频数多,以７月为最;在近４８年中,该地区的降水异常表现为３个明显不同的气候时段,可选出１０个涝年和８个旱年。     

Paleomagnetic result of the Lower-Cretaceous from Luanping basin,Hebei Province and its tectonic implications*

A paleomagnetic study of about 95 samples from 16 sites sampled in the Early Cretaceous in Luanping basin in Hebei Province was reported. Stepwise thermal demagnetization was used to isolate magnetic components. Most samples have a characteristic direction with a high temperature component above 500°C. The tectonic-corrected data areD = 347.8°,I = 50.4°, α₉₅ = 7.l°, and the corresponding pole position is at 76.1°N, 346.3°E,with dp =6.4°,dm = 3.8°, paleolatitude λ = 31.1°N. This result indicates a counterclockwise post-Cretaceous rotation of 30.7° ±9.8° with respect to the stable Ordos basin in the west of North China Block, and a non-significant northward motion. This rotation could be related to local fault action or structural detachment, or regional NNW-NWWward motion and collision of Kula-Pacific plate with eastern China since the Early Cretaceous.     

Recent Researches on the Short-Term Climate Prediction at IAP－A Brief Review

回顾了近年来在中国科学院大气物理研究所开展的有关短期气候预测研究的进展。第一个短期气候数值预测是曾庆存等利用一个耦合了热带太平洋海洋环流模式的全球大气环流模式作出的。1997年，一个基于海气耦合模式的ENSO预测系统，包括一个海洋初始化方案被建立起来，同时也开展了基于海温异常的东亚气候可预测性研究。利用气候变动的准两年信号，王会军等提出了一个可以显著改进模式预测准确率的模式结果修正方案。为了考虑土壤湿度的初始异常对夏季气候的影响，一个利用大气资料如温度、降水等经验地反演土壤湿度的方法也被建立起来。还通过一系列的数值试验研究了 1998年夏季大水发生当中海温异常和大气环流初始异常的作用。     

山西省不同地貌形态类型区土地覆被变化的GIS分析

基于全国1:100万数字地貌数据库,全国2000、2005、2010年土地覆被数据和山西省行政区边界数据,采用土地覆被动态度、转移概率矩阵和地貌面积频度方法,分析了2000-2010年山西省基本地貌形态类型下的土地覆被变化。结果表明：（1）山西省的土地覆被类型以耕地、林地和草地为主,面积总和占省域面积的95%以上。面积变化的总体趋势是耕地和草地缩减,林地和建设用地持续增加。通过分析转移概率矩阵,表明2000-2005年和2005-2010年2个时段内,土地利用的主要转化趋势相似,即林地与草地之间大面积的相互转化,耕地、草地分别与建设用地之间的相互转化。（2）耕地和建设用地的分布,随着起伏度的增大而逐渐减小。林地分布随着起伏度的增加而逐步增大。从面积变化来看,耕地、林地和草地变化主要位于中起伏山地;水域变化主要集中在大起伏山地;建设用地变化主要分布于平原和台地;未利用地变化主要在小起伏山地。从动态度来看,耕地、草地和水域动态度,在大起伏山地地区最大;林地动态度在平原区最大;建设用地的动态度最大位于丘陵;未利用地在小起伏山地地区最大。（3）土地覆被类型之间的主要转化在地貌上也存着差异。平原地区林地的增加主要来自草地转化。台地和丘陵地区3个主要转化类型相似,2000-2005年主要是林地向草地退化,2005-2010年主要是林地和草地的相互转化。在小起伏山地地区,主要是草地和林地之间的相互转化,但草地转化为林地的面积较林地向草地的转化多。中起伏和大起伏山地地区主要以草地和耕地向林地的转化为主,也有较高比例的建设用地转化为林地。封山育林和退耕还林主要是在起伏较大的山地地区,毁林主要发生在地形平缓地区。     

青藏高原东南缘活动断裂的地震应变能释放系统研究

本文以青藏高原东南缘为研究区域,利用G-R震级能量经验公式和Benioff地震应变能释放曲线,对该区域内1500年以来的历史地震应变能释放进行了系统性的研究。文中给出了各断裂带和断块区的地震应变能释放周期表,及相应的地震危险性。分析发现研究区域地震应变能的释放具有东强西弱,南强北弱的特征,整体上各断层断块区的历史地震应变能释放符合准周期模式,某些断层和断块区上的地震周期具有某种程度上的同步现象。青藏高原东南缘现今处于大释放期中,地震的危险性不能忽视。局部结果显示,安宁河-则木河断裂带、小江断裂带的危险性很高,对于这些危险区要重点跟踪研究。今后仍需结合不同研究方法来提高地震危险性评估的可靠性。     

方向观测后方交会的点位精度

以虚拟误差方程,经严密推证,给出了以方向作为独立观测值的后方交会点位精度估算公式。以大桥施工控制网为例,介绍了３个方向与４个方向后方交会的精度状况及其应用。     

中国数字化地震观测系统

介绍了中国地震局“九五”重点项目中国数字化地震观测系统，特别重点介绍了中国数字化地震台网、国家标准数字化地震台、中国地震局卫星数据通信网、中国数字地震台网中心的建设规模、设计思路、技术参数、设备结构、服务功能，以及中国数字地震台网的运行流程和运行状况。     

Precise revision of the garnet-muscovite geothermometer

The garnet-muscovite geothermometer was refined through empirical calibration by using natural rocks metamorphosed under the physical conditions of 238—1306 MPa and 490—700℃. Input temperatures and pressures were determined through simultaneously applying the garnet-biotite geothermometer and the garnet-biotite-plagioclase-quartz barometer, assuming that all FeO in muscovite and garnet be ferrous. Garnet was treated as the asymmetric quaternary solid solution, and muscovite as the symmetric binary solid solution. Input muscovite compositions include Fe atoms between 0.03—0.19 and Mg atoms between 0.04—0.16 on the basis of 11 oxygen atoms, and input garnet compositions include spessartine fractions between 0.01—0.289, grossular fractions between 0.028—0.273, and the Fe/Mg ratio between 3.387-18.986. The resulting garnet-muscovite geothermometer reproduces temperatures within (50℃ compared with the garnet-biotite thermometer. Total random error of ±37℃ of the new thermometer may stem from the pressure uncertainty of ±200 MPa, and uncertainties of ±5% of Fe and Mg components in muscovite, and ±5% of Fe, Mg, Mn and Ca components in garnet, altogether. When there exist 10%, 20%, 30%, 40% and 50% Fe3+ in muscovite, respectively, the computed garnet-muscovite temperatures will be 1—6℃, 2—12℃, 3—16℃, 5—24℃ and 7—29℃, respectively, lower than those obtained when assuming that all FeO be ferrous. The new garnet-muscovite geothermometer can efficiently reflect temperature change of typical prograde sequences and contact aureole rocks, and may be applied to low- to high-grade and low- to high-pressure metamorphic rocks.     

Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.     

利用流固耦合模型模拟地震和海啸全过程

海底地震引起的海啸过程在力学上是一个流固耦合问题。地震引起的海底变形会影响流体的运动,流体运动会影响地震引起的海底变形。海啸的数值模拟,通常采用浅水波控制方程,把地震引起的海底变形作为海啸波动的边界条件或初始条件,不考虑它们之间的相互作用。本文采用势流体的流固耦合有限元方法模拟了地震和海啸的全过程。地震过程的模拟与地震位错模型不同,在位错模型中,断层的位错是事先指定的;而在本文中,首先形成自重作用下的初始应力场,然后通过断层材料的突然软化引起的错动,模拟地震震源的动力学过程。模拟结果显示,在海面除了可以看到大振幅的海啸波外,还可以发现体波震相和面波震相。在600 km的海面震中距上,它们要比海啸波早到48分钟,在此处面波的最大平均振幅可达0.55 m,是相同震中距海底面波最大平均振幅的2倍。因此,海啸预警信息在海面可以比在地表更早地得到。海啸波的传播速度在水深3 km的开阔海面是175.8 m/s,它要比理想长波理论预测的大,其平均振幅为2 m,波长可达32 km. 到达大陆架后其速度、波长都减小,在岸边可以激起10 m高的巨浪,水平方向深入陆地达53 m。震中附近海面和地震断层上的最大垂直加速度分别为5.9 m/s2和16.3 m/s2,后者是前者的2.8倍。由此看来,海水是很好的减震器。海啸波的加速度到达岸边会衰减10倍。与加速度不同,海面震中处的振动速度为3.2 m/s, 是海底震源处的1.4倍。震源处的最大位移小于震中海面的最大位移, 其差就是海啸波源的振幅。值得注意的是,海底地震的最大位错在震后23 s达到,不是发生在断层滑动的开始。