变参考慢度Born近似傅氏偏移

针对常规Born近似傅氏偏移方法对于剧烈横向变速介质不能精确成像的状况而提出了变参考慢度Born近似傅氏偏移，理论上解决了任意速度变化地质模型的偏移成像问题。此外，为进一步提高复杂地层的成像精度和波场延拓算子的稳定性，对散射波场的计算公式作了改进。将改进的方法运用于盐丘模型的正演和偏移试验，并与常规Born近似偏移方法相比较，可明显看出变参考慢度Born近似傅氏偏移方法在效果上要优于后者，其处理速度横向变化的能力大大增强。     

在数值预报中应用卫星反演资料的质量控制研究

资料变分同化作为提高数值预报准确率的有效手段之一，受到了越来越多的关注。如何将分辨率较高的卫星资料应用到数值预报中，弥补常规观测网时空分辨率的不足，成为数值天气预报的重要课题。以1996年6月28-29日发生在江淮流域的暴雨过程为例，依据其中尺度场和能量场特征，分别以低云量、急流、湿静力能为标准设计了三套质量控制方案。将质量控制后的反演湿度场变分同化并应用到中尺度MM5模式中，结果表明：经质量控制后，可以明显地改善降水区域及降水量的预报；对降水预报，特别是强降水天气过程的预报，该质量控制方案具有普遍意义。     

GPS折射角资料的变分同化试验

越来越多的新型观测资料为数值天气预报水平的进一步提高提供了许多新的机会。在各种新型的观测资料中,GPS(全球定位卫星系统)折射角资料无疑是非常重要的。GPS折射角资料具有分辨率高、全天候探测、覆盖全球等优点,实现对GPS折射角资料的变分同化,将具有非常重要的意义。文中介绍了如何获得及同化GPS折射角资料的原理。对GPS折射角资料的变分同化可以分为两种:间接同化和直接同化,文中对这两种方法都作了具体介绍。在变分同化的最小化过程中,计算效率无疑是最重要的,而优化步长的计算又直接关系到算法效率的成败。根据最小化算法的特点,通过数学推导,得出一种适合于各种最小化算法的计算优化步长的自适应方法。最后,还利用1995年10月11日的GPS折射角资料进行了数值试验,结果表明了变分同化方法和计算优化步长方法的有效性。     

Assimilating data into open ocean tidal models

The problem of deriving tidal fields from observations by reason of incompleteness and imperfectness of every data set practically available has an infinitely large number of allowable solutions fitting the data within measurement errors and hence can be treated as ill-posed. Therefore, interpolating the data always relies on some a priori assumptions concerning the tides, which provide a rule of sampling or, in other words, a regularization of the ill-posed problem. Data assimilation procedures used in large scale tide modeling are viewed in a common mathematical framework as such regularizations. It is shown that they all (basis functions expansion, parameter estimation, nudging, objective analysis, general inversion, and extended general inversion), including those (objective analysis and general inversion) originally formulated in stochastic terms, may be considered as utilizations of one of the three general methods suggested by the theory of ill-posed problems. The problem of grid refinement critical for inverse methods and nudging is discussed.     

扩展卡尔曼滤波数据同化试验--基于Lorenz(1960)系统的研究

本文采用Lorenz(1960)系统，在只考虑初始误差及观测误差而不考虑模式误差的情况下，利用扩展卡尔曼滤波(Extended Kalman Filter)数据同化方法进行了数值模拟试验。数值试验的结果表明：扩展卡尔曼滤波数据同化方法对系统状态的估计有较好的改善作用，能有效的抑制估计误差的增长；加大观测频率，可以进一步改善数据同化的效果，使估计误差进一步减小；由于模式误差的存在，系统的不稳定能量会不断的累积，出现了估计误差的异常增长和计算的不连续现象，在模式预报方程中的均值演变方程加人二阶偏差纠错项，对控制估计误差的异常增长，进一步改善数据同化的效果有较明显作用。     

多普勒雷达资料4DVAR同化反演的模拟研究

利用Sun等建立的同化模式和四维变分同化方法对多普勒雷达资料反演大气风场、热力场和微物理场进行了模拟试验研究.反演的基本思路是将4DVAR同化方法应用到三维云模式,定义价值函数表征雷达资料与模式预报结果之间的差别,通过极小化价值函数得到反演场,价值函数相对模式控制变量的梯度由伴随模式求取.试验结果表明,4DVAR同化技术能够从单(双)多普勒雷达资料反演大气三维风场、热力场和微物理场.各个变量反演精度高低与同化过程中变量受约束的大小程度呈正相关.速度场和雨水场反演精度较高,温度场、云水和水汽的反演精度次之,温度场的准确反演需要较长的同化时间.价值函数中加入背景场,哪怕是单点探空给出的平均场信息也有利于提高反演精度.在采用单部多普勒雷达资料进行反演时,速度场的反演误差较大.反演区相对雷达站的位置变化对速度场反演结果有一定的影响,而对其他变量的反演影响很小.两个时次的雷达观测资料基本足够提供反演所需的时间演变信息,同化更多时次的雷达资料,反演效果改进很小.雷达观测资料的缺值会显著降低同化效果,甚至可能导致同化失败,引入背景场可以改善这一状况.4DVAR同化技术对于雷达观测资料误差不太敏感.利用双多普勒雷达合成风场提供水平风场边界条件是比较准确可靠的.在反演主体离边界较远时,VAD风场也基本可用作水平风场边界条件.微物理场的反演对模式中的微物理参数化方案较敏感.     

利用ARGO资料改进海洋资料同化和海洋模式中的物理过程

国际ARGO计划的实施每年将可提供多达 10万个剖面 (0～ 2 0 0 0m水深 )的海水温度和盐度资料 ,这些资料的获取无疑将会大大促进海洋和大气科学的发展 ,使人们加深对海洋过程的了解 ,揭示海 气相互作用的机理 ,为长期天气预报和短期气候预测提供模式初始场 ,提高长期天气预报和短期气候预测的能力。如何利用这些资料开展研究工作以及在实际业务中应用这些资料是目前大气和海洋科学界的一个前沿课题。本研究将ARGO浮标资料引入了国家气候中心的NCC GODAS同化系统 ,结果分析表明 ,同化ARGO资料后所得到的海温场在三大洋中不仅在温度数值的大小 ,而且在分布形式方面都与观测场具有较好的一致性 ,可以很好地反映出观测到的冬季和夏季海温的分布形式以及海温的季节变化特征和异常特征。本研究还应用最新的ARGO海洋观测资料 ,通过建立新的热带西太平洋次表层海温参数化方案 ,改进了Zebiak Cane(1987)海洋模式 (ZC模式 ) ,克服了ZC模式几乎没有模拟赤道西太平洋表层和次表层海温变化能力的缺陷。在ZC模式中引用新的次表层海温参数化方案后 ,在赤道西太平洋不仅次表层海温的模拟得到了改善 ,对海面温度异常的模拟也有了较大的改进 ,不仅模拟出了赤道西太平洋表层和次表层海温异常的年际变化特征 ,也模拟出了与观     

Silica-rich Melts in Quartz Xenoliths from Vulcano Island and their Bearing on Processes of Crustal Anatexis and Crust-Magma Interaction beneath the Aeolian Arc, Southern Italy

Quartz-rich xenoliths in lavas and pyroclastic rocks from VulcanoIsland, part of the Aeolian arc, Italy, contain silicic meltinclusions with high SiO₂ (73–80 wt %) and K₂O (3–6wt %) contents. Two types of inclusions can be distinguishedbased on their time of entrapment and incompatible trace element(ITE) concentrations. One type (late, ITE-enriched inclusions)has trace element characteristics that resemble those of themetamorphic rocks of the Calabro-Peloritano basement of theadjacent mainland. Other inclusions (early, ITE-depleted) havevariable Ba, Rb, Sr and Cs, and low Nb, Zr and rare earth element(REE) contents. Their REE patterns are unfractionated, witha marked positive Eu anomaly. Geochemical modelling suggeststhat the ITE-depleted inclusions cannot be derived from equilibriummelting of Calabro-Peloritano metamorphic rocks. ITE-enrichedinclusions can be modelled by large degrees (>80%) of meltingof basement gneisses and schists, leaving a quartz-rich residuerepresented by the quartz-rich xenoliths. Glass inclusions inquartz-rich xenoliths represent potential contaminants of Aeolianarc magmas. Interaction between calc-alkaline magmas and crustalanatectic melts with a composition similar to the analysed inclusionsmay generate significant enrichment in potassium in the magmas.However, ITE contents of the melt inclusions are comparablewith or lower than those of Vulcano calc-alkaline and potassicrocks. This precludes the possibility that potassic magmas inthe Aeolian arc may originate from calc-alkaline parents throughdifferent degrees of incorporation of crustal melts. KEY WORDS: melt inclusions; crustal anatexis; magma assimilation; xenoliths; Vulcano Island     

Petrogenesis of Tertiary Continental Intra-plate Lavas from the Westerwald Region, Germany

Tertiary volcanic rocks from the Westerwald region range frombasanites and alkali basalts to trachytes, whereas lavas fromthe margin of the Vogelsberg volcanic field consist of morealkaline basanites and alkali basalts. Heavy rare earth elementfractionation indicates that the primitive Westerwald magmasprobably represent melts of garnet peridotite. The Vogelsbergmelts formed in the spinel–garnet peridotite transitionregion with residual amphibole for some magmas suggesting meltingof relatively cold mantle. Assimilation of lower-crustal rocksand fractional crystallization altered the composition of lavasfrom the Westerwald and Vogelsberg region significantly. Thecontaminating lower crust beneath the Rhenish Massif has a differentisotopic composition from the lower continental crust beneaththe Hessian Depression and Vogelsberg, implying a compositionalboundary between the two crustal domains. The mantle sourceof the lavas from the Rhenish Massif has higher ²⁰⁶Pb/²⁰⁴Pband ⁸⁷Sr/⁸⁶Sr than the mantle source beneath the Vogelsbergand Hessian Depression. The 30–20 Ma volcanism of theWesterwald apparently had the same mantle source as the QuaternaryEifel lavas, suggesting that the magmas probably formed in apulsing mantle plume with a maximum excess temperature of 100°Cbeneath the Rhenish Massif. The relatively shallow melting ofamphibole-bearing peridotite beneath the Vogelsberg and HessianDepression may indicate an origin from a metasomatized portionof the thermal boundary layer. KEY WORDS: continental rift volcanism; basanites; trachytes; assimilation; fractional crystallization; partial melting     

Particularités de la contamination crustale des phonolites : exemple du Velay oriental (Massif central)

Sr and Nd isotopic compositions of one trachyte, eight phonolites and five basalts have been measured. The isotopic characteristics of the trachyte can be explained by a combined assimilation–fractional crystallization process within an upper crustal magmatic chamber. Some phonolites display isotopic signatures identical to basalts, suggesting that they have been protected against any crustal assimilation during their formation. Some others have low Sr contents, whereas they are enriched in radiogenic Sr (0.70451<⁸⁷Sr/⁸⁶Sr_i<0.71192), and display basaltic ¹⁴³Nd/¹⁴⁴Nd ratios. Both observations could be explained by very strong alkali feldspar fractionation and by subsequent very low assimilation of surrounding rocks (between 0.3 and 4%) during intrusion. To cite this article: J.-M. Dautria et al., C. R. Geoscience 336 (2004).