Using rainfall-runoff modeling to interpret lake level data   总被引：2，自引：0，他引：2  

Lars Bengtsson  Joakim Malm 《Journal of Paleolimnology》1997,18(3):235-248

Using water balance computations, the behavior of different kinds of lakes is discussed. Simple analytical expressions relating water level to hydrological conditions and lake bathymetry are given. The importance of knowing the river basin area when analyzing lake levels is stressed. A conceptual rainfall-runoff model including lake routing is used to simulate runoff and lake levels and to compute quasi-steady state conditions and long-term transient situations. It is suggested that models can be used to construct curves relating lake levels to precipitation and lake evaporation. By comparing with paleo-lake levels, the annual precipitation related to these levels can be found, provided information is available about the seasonal distribution of the precipitation.  相似文献   

THE RETREAT OF THE FOREST IN SOUTHEAST ASIA: A CARTOGRAPHIC ASSESSMENT     

Stéphane Bernard  Rodolphe De  Koninck 《Singapore journal of tropical geography》1997,17(1):1-14

Rapid deforestation is a major problem throughout the tropical world. The conditions and the pace under which societies and economies of the Third World are currently evolving and growing, combined with the specificities of tropical forests, render the latter increasingly vulnerable. Among the major tropical areas of the world, Southeast Asia is perhaps the one where these conditions have had the most impact on the retreat of the forest cover over the last quarter of this century. This is illustrated through the presentation of two maps of the distribution of five basic forest formations in Southeast Asia circa 1970 and circa 1990. The maps are examined and compared, as well as confronted with statistical assessments of deforestation. Finally, the complex causes behind the retreat of the tropical forests as well as the implications of this retreat are briefly discussed.  相似文献   

410 km discontinuity sharpness and the form of the olivine α-β phase diagram: resolution of apparent seismic contradictions     

George R. Helffrich  Bernard J. Wood 《Geophysical Journal International》1996,126(2):F7-F12

  相似文献   

The replacement of anthophyllite by jimthompsonite: a model for hydration reactions in biopyriboles     

Bernard Henry Grobéty 《Contributions to Mineralogy and Petrology》1997,127(3):237-247

Anthophyllite crystals found in ultramafic lenses of the Lepontine Alps (Switzerland) contain coherent, submicroscopic intergrowths of ordered and disordered biopyribole polysomes. The chain width distributions of disordered polysomes were analyzed using high resolution transmission electron microscopy (HRTEM). Chains wider than triple were interpreted as intermediate products in the transformation of anthophyllite to the triple chain silicate jimthompsonite. The concentration of individual chain types is strongly correlated with the reaction progress. Based on observed zipper terminations and the transformation rules given by Veblen and Buseck (1980) a scheme of possible reaction paths leading from anthophyllite to jimthomp sonite is proposed. The reaction scheme and a simple kinetic model for elementary reactions allow modeling of the observed chain width distributions. The model suggests that the complex reaction paths involving steps with increasing and decreasing chain width are more important in the formation of jimthompsonite than the direct transformation from anthophyllite. The wide chains (>triple) occurring as intermediate products of the multi-step paths are structurally closer to talc than jimthompsonite. The back-transformation of these wide chains to triple chains is, therefore, a strong argument that jimthompsonite is a stable phase and not only a metastable intermediate product in the transformation of anthophyllite to talc. Received: 8 July 1996 / Accepted: 13 December 1996  相似文献   

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

Nils-Otto Kitterrød  Lars Gottschalk 《Stochastic Environmental Research and Risk Assessment (SERRA)》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.  相似文献   

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

Bathymetric and geochemical investigation of Kawah Ijen Crater Lake, East Java, Indonesia     

B. Takano  K. Suzuki  K. Sugimori  T. Ohba  S. M. Fazlullin  A. Bernard  S. Sumarti  R. Sukhyar  M. Hirabayashi 《Journal of Volcanology and Geothermal Research》2004,135(4):299-329

A bathymetric survey of Kawah Ijen crater lake was conducted by acoustic sounding in 1996 to compare the lake morphology with those measured in 1922, 1925 and 1938, and to calculate the present lake volume. Even though the lake experienced several hydrothermal eruptions, the maximum depth became shallower (182 m) than before (200 m), resulting in a reduced lake volume (3.0×10⁷ m³).Fifty-two major and minor constituents including rare earth elements and polythionates (PT) of the lake waters at various depths were determined by ICP-AES, ICP-MS and HPLC, respectively. These ions except for several volatile elements are taken up by lake fringe through congruent dissolution of pyroclastics of Kawah Ijen volcano. Most ions are homogeneously distributed throughout the lake, although PT showed a considerable vertical variation. Rare earth elements (REE) in the Kawah Ijen water as well as those from other hyper-acidic crater lakes show distribution patterns likely due to the three rock dissolution (preferential, congruent and residual) types, and their logarithmic concentrations linearly depend upon the pH values of the lake waters.Using the PT degradation kinetics data, production rates of PT, injection rates of SO₂ and H₂S into the lake were estimated to be 114, 86 and 30 tons/day, respectively. Also travel time of the spring water at the Banyupahit Riverhead from Kawah Ijen was estimated to be 600–1000 days through the consideration of decreasing rates of PT. Molten sulfur stocks containing Sn, Cu, Bi sulfides and Pb-barite exposed on the inner crater slope were presumed to be extinct molten sulfur pools at the former lake bottom. This was strongly supported by the barite precipitation temperature estimated through the consideration of the temperature dependence of Pb-chlorocomplex formation.  相似文献   

Multicomponent ionic diffusion in porewaters: Coulombic effects revisited     

Bernard P. Boudreau  Filip J.R. Meysman 《Earth and Planetary Science Letters》2004,222(2):653-666

The diffusion of an ion in porewaters cannot occur independently of the other ions in solution as a result of Coulombic coupling, as well as from other effects not considered here. Unfortunately, a longstanding disagreement exists about the correct form and meaning of the equations that describe Coulombic coupling in porewaters, i.e., Ben-Yaakov [Am. J. Sci. 281 (1981) 974] vs. Lasaga [Am. J. Sci. 281 (1981) 981]. This paper re-examines this controversy by reformulating the problem starting from fundamental concepts of mass and charge conservation. We show that these antagonistic formulations are both valid and, in fact, equivalent, when the different interpretations of charge balance are resolved. Most of the disagreements between Ben-Yaakov and Lasaga are then shown to result from differing methods of solution, not fundamental disparities in their models. We note, however, that the explanation for the concept of “stationary” gradients of nonreacting ions as given Ben-Yaakov is inaccurate, and such gradients do lead to diffusive fluxes that are counterbalanced by electrochemical migrational fluxes to produce no net flux (excluding advective flux). We further find that the bicarbonate diffusive flux will not balance the diffusional charge flux of sulfate during its reduction if advection is present. This latter imbalance generates compensating fluxes in the other nonreacting ions. We have applied our theory to a simplified case of sulfate reduction in a marine sediment. The results show that nonreacting ions do diffuse and that with normally expected values of porewater advection, the ratio of the bicarbonate to the sulfate flux can be far different than the ideal value of −2.  相似文献   

Fitting the Linear Model of Coregionalization by Generalized Least Squares   总被引：2，自引：0，他引：2  

Bernard Pelletier  Pierre Dutilleul  Guillaume Larocque  James W. Fyles 《Mathematical Geology》2004,36(3):323-343

In geostatistical studies, the fitting of the linear model of coregionalization (LMC) to direct and cross experimental semivariograms is usually performed with a weighted least-squares (WLS) procedure based on the number of pairs of observations at each lag. So far, no study has investigated the efficiency of other least-squares procedures, such as ordinary least squares (OLS), generalized least squares (GLS), and WLS with other weighing functions, in the context of the LMC. In this article, we compare the statistical properties of the sill estimators obtained with eight least-squares procedures for fitting the LMC: OLS, four WLS, and three GLS. The WLS procedures are based on approximations of the variance of semivariogram estimates at each distance lag. The GLS procedures use a variance–covariance matrix of semivariogram estimates that is (i) estimated using the fourth-order moments with sill estimates (GLS₁), (ii) calculated using the fourth-order moments with the theoretical sills (GLS₂), and (iii) based on an approximation using the correlation between semivariogram estimates in the case of spatial independence of the observations (GLS₃). The current algorithm for fitting the LMC by WLS while ensuring the positive semidefiniteness of sill matrix estimates is modified to include any least-squares procedure. A Monte Carlo study is performed for 16 scenarios corresponding to different combinations of the number of variables, number of spatial structures, values of ranges, and scale dependence of the correlations among variables. Simulation results show that the mean square error is accounted for mostly by the variance of the sill estimators instead of their squared bias. Overall, the estimated GLS₁ and theoretical GLS₂ are the most efficient, followed by the WLS procedure that is based on the number of pairs of observations and the average distance at each lag. On that basis, GLS₁ can be recommended for future studies using the LMC.  相似文献   

Réponse aux commentaires d'Azambre et al. et de Rapaille et al.: [C. R. Geoscience 334 (2002) 197]     

Christian Desreumaux  Bernard Clément  Richard Fabre  Bruno Martins-Campina 《Comptes Rendus Geoscience》2004,336(2):173-174

  相似文献