扩展卡尔曼滤波数据同化试验--基于Lorenz(1960)系统的研究   总被引：4，自引：0，他引：4       下载免费PDF全文

费剑锋  韩月琪  王云峰  侯志明 《气象科学》2004,24(4):413-423

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

多普勒雷达资料4DVAR同化反演的模拟研究   总被引：21，自引：5，他引：21  

许小永  郑国光  刘黎平 《气象学报》2004,62(4):410-422

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

利用ARGO资料改进海洋资料同化和海洋模式中的物理过程   总被引：8，自引：1，他引：8  

张人禾  刘益民  殷永红  史历 《气象学报》2004,62(5):613-622

国际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     

FREZZOTTI  MARIA-LUCE; PECCERILLO  ANGELO; ZANON  VITTORIO; NIKOGOSIAN  IGOR 《Journal of Petrology》2004,45(1):3-26

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   总被引：2，自引：2，他引：2  

HAASE  KARSTEN M.; GOLDSCHMIDT  BJORN; GARBE-SCHONBERG  C.-DIETER 《Journal of Petrology》2004,45(5):883-905

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  相似文献   

Indicator Simulation Accounting for Multiple-Point Statistics   总被引：7，自引：0，他引：7  

Julián M. Ortiz  Clayton V. Deutsch 《Mathematical Geology》2004,36(5):545-565

Geostatistical simulation aims at reproducing the variability of the real underlying phenomena. When nonlinear features or large-range connectivity is present, the traditional variogram-based simulation approaches do not provide good reproduction of those features. Connectivity of high and low values is often critical for grades in a mineral deposit. Multiple-point statistics can help to characterize these features. The use of multiple-point statistics in geostatistical simulation was proposed more than 10 years ago, on the basis of the use of training images to extract the statistics. This paper proposes the use of multiple-point statistics extracted from actual data. A method is developed to simulate continuous variables. The indicator kriging probabilities used in sequential indicator simulation are modified by probabilities extracted from multiple-point configurations. The correction is done under the assumption of conditional independence. The practical implementation of the method is illustrated with data from a porphyry copper mine.  相似文献   

Transformation of Residuals to Avoid Artifacts in Geostatistical Modelling with a Trend   总被引：1，自引：0，他引：1  

Oy Leuangthong  Clayton V. Deutsch 《Mathematical Geology》2004,36(3):287-305

Trend modelling is an important part of natural resource characterization. A common approach to account for a variable with a trend is to decompose it into a relatively smoothly varying trend and a more variable residual component. Then, the residuals are stochastically modelled independent of the trend. This decomposition can result in values outside the plausible range of variability, such as grades below zero or ratios that exceed 1.0. We transform the residuals conditional to the trend component to explicitly remove these complex features prior to geostatistical modelling. Back transformation of the modelled residual values allows the complex relations to be reproduced. A petroleum-related application shows the robustness of the proposed transformation. Furthermore, a mining application shows that when this conditional transformation is applied to the original variable, instead of the residual, simulated values are assured to be nonnegative.  相似文献   

鹅湖岩体岩石谱系单位的建立及其就位机制分析     

吴文革  谢卫红 《城市地质》2004,16(3):12-17,23

鹅湖岩体为S型花岗岩 ,具有结构和成分演化序列。运用岩石结构序列法 ,将其划分为 5个单元 ,归并为 1个超单元。按侵入体空间展布形式、岩体地质学、岩石学、构造等资料 ,探讨了岩体就位机制  相似文献   

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

Jean-Marie Dautria  Jean-Michel Liotard  Danièle Briot 《Comptes Rendus Geoscience》2004,336(11):971-981

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).  相似文献   

The Application of ARGO Data to the Global Ocean Data Assimilation Operational System of NCC   总被引：11，自引：0，他引：11       下载免费PDF全文

LIU Yimin  ZHANG Renhe  YIN Yonghong  NIU Tao 《Acta Meteorologica Sinica》2005,19(3):355-365

1. Introduction It is well-known that the state of ocean plays very important role in the climate change. But there is a paucity of the ocean observation data. The data distri- bution in the space, time and different components is very inhomogeneous, even in some areas, there are no any observation data. Hence, it brings some diffcul- ties to the scientists to study many problems relevant to ocean. This situation has been being changed since ARGO (Array for Real-time Geostrophic Oceanogra-…  相似文献