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)数据同化方法进行了数值模拟试验。数值试验的结果表明：扩展卡尔曼滤波数据同化方法对系统状态的估计有较好的改善作用，能有效的抑制估计误差的增长；加大观测频率，可以进一步改善数据同化的效果，使估计误差进一步减小；由于模式误差的存在，系统的不稳定能量会不断的累积，出现了估计误差的异常增长和计算的不连续现象，在模式预报方程中的均值演变方程加人二阶偏差纠错项，对控制估计误差的异常增长，进一步改善数据同化的效果有较明显作用。     

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     

Indicator Simulation Accounting for Multiple-Point Statistics

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

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.     

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

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

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).     

The Application of ARGO Data to the Global Ocean Data Assimilation Operational System of NCC

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-…     

Establishment of apparent quantum yield and maximum ecosystem assimilation on Tibetan Plateau alpine meadow ecosystem

The alpine meadow is widely distributed on the Tibetan Plateau with an area of about 1.2×106kn2. Damxung County, located in the hinterland of the Tibetan Plateau, is the place covered with this typical vegetation. An open-path eddy covariance system was set up in Damxung rangeland station to measure the carbon flux of alpine meadow from July to October,2003. The continuous carbon flux data were used to analyze the relationship between net ecosystem carbon dioxide exchange (NEE) and photosynthetically active radiation (PAR), as well as the seasonal patterns of apparent quantum yield (α) and maximum ecosystem assimilation (Pmax).Results showed that the daytime NEE fitted fairly well with the PAR in a rectangular hyperbola function, with α declining in the order of peak growth period (0.0244 μmolCO2 · μmol-1pAR) ＞early growth period ＞ seed maturing period ＞ withering period (0.0098 μmolCO2 · μmol-1pAR).The Pmax did not change greatly during the first three periods, with an average of 0.433mgCO2· m-2· s-1, i.e. 9.829 μmolCO2· m-2· s-1. However, during the withering period, Pmax was only 0.35 mgCO2 · m-2 · s-1, i.e. 7.945 μmolCO2 · m-2 · s-1. Compared with other grassland ecosystems, the α of the Tibetan Plateau alpine meadow ecosystem was much lower.