Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging     

Nils-Otto Kitterr?d  Lars Gottschalk 《Stochastic Hydrology and Hydraulics》1997,11(6):459-482

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.  相似文献   

Dynamik der Schilfröhrichte am Bodensee unter dem Einfluss von Wasserstandsvariationen     

Michael Dienst  Klaus Schmieder  Wolfgang Ostendorp 《Limnologica》2004,34(1-2):29

On the basis of different sets of aerial photos the dynamics of the reed bed areas of Lake Constance were investigated in relation to the dynamics of the water levels. The objectives of the study were to quantify the changes of reed areas due to different flood events in the last decades and their recovery in the time periods between these events. The results should given information of the relevance of water level variations on reed bed dynamics and the regeneration times of reed beds after extreme disturbance events.Following the extreme flood at Lake Constance in 1999 the reed belts of Lake Constance lost approximately 30 ha (24%) of the lakeside reed beds. The loss is comparable to the situation in the late 1960s, when approximately 40 ha died back due to the extreme flood in 1965 and the high spring water levels in the subsequent years. In the time period between the extreme floods of 1965 and 1999, the reed areas expanded to nearly 85% of the area before 1965. The expansion rates increased with increasing distance to the flood event of 1965. Especially in periods with series of years of low spring water level the expansion rates were high.The damage degrees of the reed areas in the years 2000 and 2002 showed a clear relation to the elevation (i.e. average water level) of the stands. The damage degree increased with decreasing elevation. Furthermore the regeneration process of severely damaged stands was related to the elevation level of the stands. Whereas stands at high elevation regenerate fast, those at low elevation died off completely in the years after the extreme flood. This supports the hypothesis that the water level flutuations play a major role in the reed dynamics of Lake Constance.As a consequence of the climate change an increase in the frequency of high spring water levels is expected. Thus, it seems unlikely that reed stands will ever expand again to the same area as before 1965.  相似文献   

Structure and tectonics of convergent plate margins     

Aleš Špi?ák  Ond?ej ?adek  E.Robert Engdahl 《Physics of the Earth and Planetary Interiors》2004,141(4):241

  相似文献   

Influence of the zone of weakness on dip angle and shear heating of subducted slabs     

Martin Kuka?ka  Ctirad Matyska 《Physics of the Earth and Planetary Interiors》2004,141(4):243-252

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.  相似文献   

A Modeling Study of the Effects of Anomalous Snow Cover over the Tibetan Plateau upon the South Asian Summer Monsoon   总被引：1，自引：0，他引：1  

刘华强  孙照渤王举 闵锦忠 《大气科学进展》2004,21(6):964-975

The effect of anomalous snow cover over the Tibetan Plateau upon the South Asian summer monsoon is investigated by numerical simulations using the NCAR regional climate model (RegCM2) into which gravity wave drag has been introduced. The simulations adopt relatively realistic snow mass forcings based on Scanning Multi-channel Microwave Radiometer (SNINIR) pentad snow depth data. The physical mechanism and spatial structure of the sensitivity of the South Asian early summer monsoon to snow cover anomaly over the Tibetan Plateau are revealed. The main results are summarized as follows. The heavier than normal snow cover over the Plateau can obviously reduce the shortwave radiation absorbed by surface through the albedo effect, which is compensated by weaker upward sensible heat flux associated with colder surface temperature, whereas the effects of snow melting and evaporation are relatively smaller.The anomalies of surface heat fluxes can last until June and become unobvious in July. The decrease of the Plateau surface temperature caused by heavier snow cover reaches its maximum value from late April to early May. The atmospheric cooling in the mid-upper troposphere over the Plateau and its surrounding areas is most obvious in May and can keep a fairly strong intensity in June. In contrast, there is warming to the south of the Plateau in the mid-lower troposphere from April to June with a maximum value in May.The heavier snow cover over the Plateau can reduce the intensity of the South Asian summer monsoon and rainfall to some extent, but this influence is only obvious in early summer and almost disappears in later stages.  相似文献   

The measurement of sulfate mineral solubilities in the Na-K-Ca-Cl-SO₄-H₂O system at temperatures of 100, 150 and 200°C     

Daniela Freyer  Wolfgang Voigt 《Geochimica et cosmochimica acta》2004,68(2):307-318

At T > 100°C development of thermodynamic models suffers from missing experimental data, particularly for solubilities of sulfate minerals in mixed solutions. Solubilities in Na⁺-K⁺-Ca²⁺-Cl⁻-SO₄²⁻/H₂O subsystems were investigated at 150, 200°C and at selected compositions at 100°C. The apparatus used to examine solid-liquid phase equilibria under hydrothermal conditions has been described.In the system NaCl-CaSO₄-H₂O the missing anhydrite (CaSO₄) solubilities at high NaCl concentrations up to halite saturation have been determined. In the system Na₂SO₄-CaSO₄-H₂O the observed glauberite (Na₂SO₄ · CaSO₄) solubility is higher than that predicted by the high temperature model of Greenberg and Møller (1989), especially at 200°C. At high salt concentrations, solubilities of both anhydrite and glauberite increase with increasing temperature. Stability fields of the minerals syngenite (K₂SO₄ · CaSO₄ · H₂O) and goergeyite (K₂SO₄ · 5 CaSO₄ · H₂O) were determined, and a new phase was found at 200°C in the K₂SO₄-CaSO₄-H₂O system. Chemical and single crystal structure analysis give the formula K₂SO₄ · CaSO₄. The structure is isostructural with palmierite (K₂SO₄ · PbSO₄). The glaserite (“3 K₂SO₄ · Na₂SO₄”) appears as solid solution in the system Na₂SO₄-K₂SO₄-H₂O. Its solubility and stoichiometry was determined as a function of solution composition.  相似文献   

Tetrahedral Plots of the Phase Relations for Basalts     

Sven Maal?e  Richard N. Abbott Jr. 《Mathematical Geology》2005,37(8):869-893

The phase relations of quaternary systems are generally represented by projections onto ternary compositional planes. Such projections often obscure relationships that would only be evident in a three-dimensional tetrahedral plot. The tetrahedral plot requires that compositions of the minerals and melts be transformed into Cartesian coordinates. It is shown here how this transformation is carried out. The application is demonstrated by tetrahedral plots of experimental melt compositions of partially molten lherzolite. Furthermore, the plot can be used to evaluate whether or not a particular basaltic composition represents a primary melt. The methods are applicable to any four-component system.  相似文献   

On the Sr isotope and REE compositions of anhydrites from the TAG seafloor hydrothermal system   总被引：1，自引：0，他引：1  

Susan E. Humphris  Wolfgang Bach 《Geochimica et cosmochimica acta》2005,69(6):1511-1525

  相似文献   

Spatially controlled Fe and Si isotope variations: an alternative view on the formation of the Torres del Paine pluton     

Norbert?A.?GajosEmail author  Craig?C.?Lundstrom  Alexander?H.?Taylor 《Contributions to Mineralogy and Petrology》2016,171(11):93

We present new Fe and Si isotope ratio data for the Torres del Paine igneous complex in southern Chile. The multi-composition pluton consists of an approximately 1 km vertical exposure of homogenous granite overlying a contemporaneous 250-m-thick mafic gabbro suite. This first-of-its-kind spatially dependent Fe and Si isotope investigation of a convergent margin-related pluton aims to understand the nature of granite and silicic igneous rock formation. Results collected by MC-ICP-MS show a trend of increasing δ⁵⁶Fe and δ³⁰Si with increasing silica content as well as a systematic increase in δ⁵⁶Fe away from the mafic base of the pluton. The marginal Torres del Paine granites have heavier Fe isotope signatures (δ⁵⁶Fe = +0.25 ± 0.02 2se) compared to granites found in the interior pluton (δ⁵⁶Fe = +0.17 ± 0.02 2se). Cerro Toro country rock values are isotopically light in both Fe and Si isotopic systems (δ⁵⁶Fe = +0.05 ± 0.02 ‰; δ³⁰Si = ?0.38 ± 0.07 ‰). The variations in the Fe and Si isotopic data cannot be accounted for by local assimilation of the wall rocks, in situ fractional crystallization, late-stage fluid exsolution or some combination of these processes. Instead, we conclude that thermal diffusion or source magma variation is the most likely process producing Fe isotope ratio variations in the Torres del Paine pluton.  相似文献   

封面及目录     

本刊编辑部 《岩矿测试》2016,35(4):I-IV

  相似文献