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.     

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     

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

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.     

Data assimilation (4D-VAR) to forecast flood in shallow-waters with sediment erosion

In this paper, the four-dimensional variational data assimilation technique (4D-VAR) is presented as a tool to forecast floods. Our study is limited to purely hydrological flows and supposes that the weather, here a big rain, has been already forecasted by meteorological services. The technique consists in minimizing, in the sense of Lagrange, the cost function: a measure of the difference between calculated data and available observations, here the water level. This is done under constraints that are the equations of the physical model. In our case, we modified the shallow-water equations to include a simplified sediment transport model. The steepest descent algorithm is then used to find the minimum. This is made possible because we can compute analytically the gradient of the cost function by using the adjoint equations of the model. As an application of the 4D-VAR technique, the overflowing of the Chicoutimi River at the Chute-Garneau dam, during the 1996 flood, is investigated. It is found that the 4D-VAR method reduces the error in the water height forecast even when the erosion model is not activated. In terms of Lyapunov exponents, we estimate the predictability horizon of such an event to be about half-an-hour after a big rain. However, this limit of predictability can be increased by using more observations or by using a finer computational grid.     

Model identification for hydrological forecasting under uncertainty

Methods for the identification of models for hydrological forecasting have to consider the specific nature of these models and the uncertainties present in the modeling process. Current approaches fail to fully incorporate these two aspects. In this paper we review the nature of hydrological models and the consequences of this nature for the task of model identification. We then continue to discuss the history (“The need for more POWER‘’), the current state (“Learning from other fields”) and the future (“Towards a general framework”) of model identification. The discussion closes with a list of desirable features for an identification framework under uncertainty and open research questions in need of answers before such a framework can be implemented.     

Assimilation and Fractional Crystallization Controlled by Transport Process of Crustal Melt: Implications from an Alkali Basalt-Dacite Suite from Rishiri Volcano, Japan

Mechanisms of fractional crystallization with simultaneous crustalassimilation (AFC) are examined for the Kutsugata and Tanetomilavas, an alkali basalt–dacite suite erupted sequentiallyfrom Rishiri Volcano, northern Japan. The major element variationswithin the suite can be explained by boundary layer fractionation;that is, mixing of a magma in the main part of the magma bodywith a fractionated interstitial melt transported from the mushyboundary layer at the floor. Systematic variations in SiO₂ correlatewith variations in the Pb, Sr and Nd isotopic compositions ofthe lavas. The geochemical variations of the lavas are explainedby a constant and relatively low ratio of assimilated mass tocrystallized mass (‘r value’). In the magma chamberin which the Kutsugata and Tanetomi magmas evolved, a strongthermal gradient was present and it is suggested that the marginalpart of the reservoir was completely solidified. The assimilantwas transported by crack flow from the partially fused floorcrust to the partially crystallized floor mush zone throughfractures in the solidified margin, formed mainly by thermalstresses resulting from cooling of the solidified margin andheating of the crust. The crustal melt was then mixed with thefractionated interstitial melt in the mushy zone, and the mixedmelt was further transported by compositional convection tothe main magma, causing its geochemical evolution to be characteristicof AFC. The volume flux of the assimilant from the crust tothe magma chamber is suggested to have decreased progressivelywith time (proportional to t^–1/2), and was about 3 x 10^–2m/year at t = 10 years and 1 x 10^–2 m/year at t = 100years. It has been commonly considered that the heat balancebetween magmas and the surrounding crust controls the couplingof assimilation and fractional crystallization processes (i.e.absolute value of r). However, it is inferred from this studythat the ratio of assimilated mass to crystallized mass canbe controlled by the transport process of the assimilant fromthe crust to magma chambers. KEY WORDS: assimilation and fractional crystallization; mass balance model; magma chamber; melt transport; Pb isotope     

THE EFFECTS OF BOGUS TYPHOON AND OBSERVED OCEANIC DATA ON THE ABILITY OF T_(213)L_(31) TO PREDICT TC TRACK

The abilities of typhoon (TC) track prediction by a medium-range forecast model T213L31 at National Meteorological Center are analyzed and its ability to improve its TC forecasts is discussed. The results show that about 57% of the TCs could be predicted by T213L31 but the initial position errors are large.The 43% area without the prediction of TC tracks is concentrated between 13°N and 20°N and east of 120°E and lack of conventional observation data is the main reason for the absence of TC prediction in this area.The adding of bogus TC could improve the ability of TC track prediction when there is no TC vortex in the analysis field, but could only have positive effects on the short-range TC track prediction when there is TC vortex in the T213L31 analysis field.     

“00.7”北京特大暴雨模拟中气象资料同化作用的评估

针对2000年7月4～5日北京地区的一次特大暴雨过程(24 h降水量达240 mm),文中利用MM5/WRF三维变分系统和MM5非静力模式,对此次特大暴雨过程中的各种气象监测资料(地基GPS大气柱水汽含量、常规探空、高空测风、地面常规观测和地面自动气象站)的同化作用通过观测系统数值试验进行了评估.结果表明与传统的客观分析方案相比较,MM5/WRF三维变分同化系统可直接引入非常规地基GPS大气柱水汽含量监测资料,提供更好的大气初始分析场.在三维变分同化方案下,各种大气监测资料均对改进此次特大暴雨模拟有不同程度的贡献,其中,常规探空和高空测风监测资料对改进预报结果的影响最大,地面常规观测和地面自动气象站观测资料作用次之,地基GPS大气柱水汽含量资料在与其他大气监测资料相互优势互补后,可很好地改善模式大气的分析质量,通过三维变分同化技术在区域数值天气预报模式初始场中引入地基GPS大气水汽监测网资料,使此次强降水个例的6 h和24 h测站降水预报的TS评分值在1,5,10和20 mm预报检验阈值下分别提高了1%～8%.研究结果对利用三维变分数值系统,评估气象监测网资料在改进高影响天气事件预报中的作用有借鉴意义.     

AMSR2辐射率资料同化对台风“山神”分析和预报的影响研究

在WRFDA-3DVar（Weather Research and Forecasting model's 3-dimensional variational data assimilation）的框架下,添加了新的探测器AMSR2（Advanced Microwave Scanning Radiometer 2）微波辐射率资料的同化模块,实现了AMSR2辐射率资料在中小尺度同化系统中的有效使用。台风"山神"（Son-Tinh）直接同化AMSR2资料的个例试验结果表明,AMSR2资料可以很好的探测出台风的形态,并且与没有同化该资料的控制试验相比,同化AMSR2辐射率资料可以有效提高模式分析场的质量,进一步提高了台风中心气压,最大风速和台风路径的预报。