首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
The product covariance model, the product–sum covariance model, and the integrated product and integrated product–sum models have the advantage of being easily fitted by the use of marginal variograms. These models and the use of the marginals are described in a series of papers by De Iaco, Myers, and Posa. Such models allow not only estimating values at nondata locations but also prediction in future times, hence, they are useful for analyzing air pollution data, meteorological data, or ground water data. These three kinds of data are nearly always multivariate and because the processes determining the deposition or dynamics will affect all variates, a multivariate approach is desirable. It is shown that the use of marginal variograms for space–time modeling can be extended to the multivariate case and in particular to the use of the Linear Coregionalization Model (LCM) for cokriging in space–time. An application to an environmental data set is given.  相似文献   

2.
Geostatistical space–time models are used increasingly for addressing environmental problems, such as monitoring acid deposition or global warming, and forecasting precipitation or stream flow. Each discipline approaches the problem of joint space–time modeling from its own perspective, a fact leading to a significant amount of overlapping models and, possibly, confusion. This paper attempts an annotated survey of models proposed in the literature, stating contributions and pinpointing shortcomings. Stochastic models that extend spatial statistics (geostatistics) to include the additional time dimension are presented with a common notation to facilitate comparison. Two conceptual viewpoints are distinguished: (1) approaches involving a single spatiotemporal random function model, and (2) approaches involving vectors of space random functions or vectors of time series. Links between these two viewpoints are then revealed; advantages and shortcomings are highlighted. Inference from space–time data is revisited, and assessment of joint space–time uncertainty via stochastic imaging is suggested.  相似文献   

3.
It is often observed that dynamic recrystallization results in a recrystallized grain size distribution with a mean grain size that is inversely related to the flow stress. However, it is still open to discussion if theoretical models that underpin recrystallized grain size–stress relations offer a satisfactorily microphysical basis. The temperature dependence of recrystallized grain size, predicted by most of these models, is rarely observed, possibly because it is usually not systematically investigated. In this study, samples of wet halite containing >10 ppm water (by weight) were deformed in axial compression at 50 MPa confining pressure. The evolution of the recrystallized grain size distribution with strain was investigated using experiments achieving natural strains of 0.07, 0.12 and 0.25 at a strain rate of 5×10−7 s−1 and a temperature of 125 °C. The stress and temperature dependence of recrystallized grain size was systematically investigated using experiments achieving fixed strains of 0.29–0.46 (and one to a strain of 0.68) at constant strain rates of 5×10−7–1×10−4 s−1 and temperatures of 75–240 °C, yielding stresses of 7–22 MPa. The microstructures and full grain size distributions of all samples were analyzed. The results showed that deformation occurred by a combination of dislocation creep and solution-precipitation creep. Dynamic recrystallization occurred in all samples and was dominated by fluid assisted grain boundary migration. During deformation, grain boundary migration results in a competition between grain growth due to the removal of grains with high internal strain energy and grain size reduction due to grain dissection (i.e. moving boundaries that crosscut or consume parts of neighbouring grains). At steady state, grain growth and grain size reduction processes balance, yielding constant flow stress and recrystallized grain size that is inversely related to stress and temperature. Evaluation of the recrystallized grain size data against the different models for the development of mean steady state recrystallized grain size revealed that the data are best described by a model based on the hypothesis that recrystallized grain size organizes itself in the boundary between the (grain size sensitive) solution-precipitation and (grain size insensitive) dislocation creep fields. Application of a piezometer, calibrated using the recrystallized grain size data, to natural halite rock revealed that paleostresses can vary significantly with temperature (up to a factor of 2.5 for T=50–200 °C) and that the existing temperature independent recrystallized grain size–stress piezometer may significantly underestimate flow stresses in natural halite rock.  相似文献   

4.
Comparison of kriging techniques in a space-time context   总被引：1，自引：0，他引：1  
Space-time processes constitute a particular class, requiring suitable tools in order to predict values in time and space, such as a space-time variogram or covariance function. The space-time co-variance function is defined and linked to the Linear Model of Coregionalization under second-order space-time stationarity. Simple and ordinary space-time kriging systems are compared to simple and ordinary cokriging and their differences for unbiasedness conditions are underlined. The ordinary space-time kriging estimation then is applied to simulated data. Prediction variances and prediction errors are compared with those for ordinary kriging and cokriging under different unbiasedness conditions using a cross-validation. The results show that space-time kriging tend to produce lower prediction variances and prediction errors that kriging and cokriging.  相似文献   

5.
The interpretation of DSS (deep seismic soundings) profiles in Central and Eastern Alps is recalled in the paper and the models of the lower crust and Moho proposed several years ago are compared to the results of the TRANSALP seismic reflection profile. This evaluation highlights a good agreement as far as the geometry of the deep crustal structure is concerned. Therefore, the reliability of the interpretative models, previously exclusively based on DSS profiles, becomes improved. The deep structure beneath the whole Alpine range is examined reconsidering the map of the Moho boundary and the structural model already proposed for the central-eastern sector. Five main interpretative transects are put side by side, starting from the Western Alps and moving eastwards to the Swiss–Lombardian Central Alps (“European Geotraverse”), to the cross section from southern Bavaria to the Euganei Hills, to the TRANSALP profile, and finally to the easternmost profile available so far (southern Bavaria–Trieste). The comparison outlines lateral variations of the deep crustal structure as well as a sharp contrast between the Adria and the European lower crust and Moho. The transition from the Adria plate to the Dinaric domain remains, up to now, undefined.  相似文献   

6.
Barbolini  M.  Natale  L.  Savi  F. 《Natural Hazards》2002,25(3):225-244
Dynamical models for calculating snow avalanche motion have gained growingimportance in recent years for avalanche hazard assessment. Nevertheless, inherentuncertainties in their input-data specification, although well acknowledged, areusually not explicitly incorporated into the analysis and considered in the mappingresults. In particular, the estimate of avalanche release conditions is affected bystrong uncertainties when associated to a return period. These sources of error arenormally addressed through sensitivity analysis or conservative parameters estimate.However, each of these approaches has limitations in assessing the statistical implications of uncertainties.In the present paper the problem of release scenarios randomness is looked at following a Monte Carlo procedure. This statistical sampling-analysis method allows the evaluation of the probability distributions of relevant variables for avalanche hazard assessment – such as runout distance and impact pressure – once the release variables – essentially releasedepth and release length – are expressed in terms of probability distributions, accounting explicitly for inherent uncertainties in their definition. Both the theoretical framework of this procedure and its application to a real study case are presented. As initial step of this research in the present work the attention is mainly focused on flowing avalanches descending on open slopes. Therefore, the one-dimensional version of VARA dynamic models is usedfor avalanche simulations.  相似文献   

7.
Covariance models provide the basic measure of spatial continuity in geostatistics. Traditionally, a closed-form analytical model is fitted to allow for interpolation of sample Covariance values while ensuring the positive definiteness condition. For cokriging, the modeling task is made even more difficult because of the restriction imposed by the linear coregionalization model. Bochner's theorem maps the positive definite constraints into much simpler constraints on the Fourier transform of the covariance, that is the density spectrum. Accordingly, we propose to transform the experimental (cross) covariance tables into quasidensity spectrum tables using Fast Fourier Transform (FFT). These quasidensity spectrum tables are then smoothed under constraints of positivity and unit sum. A backtransform (FFT) yields permissible (jointly) positive definite (cross) covariance tables. At no point is any analytical modeling called for and the algorithm is not restricted by the linear coregionalization model. A case study shows the proposed covariance modeling to be easier and much faster than the traditional analytical covariance modeling, yet yields comparable kriging or simulation results.  相似文献   

8.
Using the Chiller computer program, we performed modeling of the mechanisms of the joint transport and deposition of Au and Sb from various ore-forming solutions during the formation of Au-Sb deposits. Three models are considered by the example of the Uderei Au-Sb deposit in the Yenisei Ridge: (1) simple cooling (cooling only), (2) iso-enthalpy boiling (P = f(T)), and (3) solution–rock interaction (rock titration model). The behavior of Sb(III) and Au(I) in the system Au–Sb–Fe–Cu–Pb–Zn–As–H2O–Cl–H2S–CO2 under hydrothermal conditions was studied. It is shown that both weakly alkaline (near-neutral) and reduced acidic Feaq2+-enriched low-chloride high-CO2 and high-chloride hydrothermal solutions play a crucial role in the formation of gold parageneses of Au-Sb ores.  相似文献   

9.
In this paper, the condition number of the stationary kriging matrix is studied for some well-known covariance models. Indeed, the robustness of the kriging weights is strongly affected by this measure. Such an analysis can justify the choice of a covariance function among other admissible models which could fit a given experimental covariance equally well.  相似文献   

10.
This paper argues, following Friedrich Nietzsche, that recent evaluations of Google Earth uncritically privilege the product’s Apollonian determinations at the expense of its Dionysian uncertainties. Specifically, when we understand Google Earth as a virtual globe composed of surveyed panoramas, sober rationalization, dystopic control, and transparent order – or, even, as a tool for participation and empowerment – we undersell its capacities as an alluring digital peep-box, an uncertain orb spangled with vertiginous paranoia, frenzied navigation, jubilatory dissolution, and intoxicating giddiness. We argue that the former interpretations not only risk foreclosing our theorizations about how Google Earth is actually used in various ways and different contexts, they also reproduce a one-dimensional and conservative reading of technology that can be traced back (at least) to the writings of Theodor Adorno. By drawing on the work of Walter Benjamin (Adorno’s critic and pen pal for more than a decade) we aim to ‘go beyond Apollo and Adorno’ by illustrating the extent to which Apollonian order and Dionysian love makes Google Earth go round. To do this, we examine Google Earth as a “digital peep-box” with an online collective that revels in its “Spot the Black Helicopter” competitions; illuminated profanities; alien and giant insect invaders; naked sunbathers; and crashed transport planes in Darfur.  相似文献   

11.
Mineral exploration in the Neoproterozoic Goiás Magmatic Arc, central Brazil, dates back to the beginning of the 1970s. The Goiás Magmatic Arc extends for more than 1000 km in the western and northern parts of Goiás, into Tocantins, and disappears under the Phanerozoic Parnaíba Basin. Two main areas of Neoproterozoic juvenile crust, the Arenópolis and the Mara Rosa arcs, are identified. They lie in the southern and northern sectors of the Goiás Arc, respectively, and are relatively well studied.The Goiás Magmatic Arc dominantly comprises tonalitic/dioritic orthogneisses and narrow NNE-striking volcano-sedimentary belts. Recent U–Pb zircon data indicate crystallization of the tonalite protoliths in two main episodes: the older between ca. 890 and 790 Ma and the younger at 670–600 Ma. Nd isotopic data indicate the very primitive nature of the original magmas, with TDM model ages mostly within the interval between 0.9 and 1.0 Ga and Nd(T) values between +3.0 and +4.6. In the Chapada–Mara Rosa area, the supracrustal rocks form three individual NNE belts, known as the eastern, central and western belts, separated from each other by metatonalites/metadiorites.Gold and Cu–Au deposits of the Mara Rosa area occur in four main associations: (i) Au–Ag–Ba (e.g., Zacarias), which are interpreted as stratiform, disseminated volcanogenic deposits, (ii) Cu–Au (e.g., Chapada) which has been interpreted either as volcanogenic or as a porphyry-type deposit, (iii) Au-only deposits (e.g., Posse), interpreted as an epigenetic disseminated deposit controlled by a mesozonal shear zone and (iv) Au–Cu–Bi (e.g., the Mundinho occurrence), which are considered as vein-type deposits controlled by magnetite-rich diorites.The gold and Cu–Au deposits located within the Goiás Magmatic Arc can be spatially and temporally related to the magmatic evolution of a collisional belt or, in other words, to an orogenic gold deposit model. These models are based on the continuous evolution of collisional plates, which can be subdivided into four stages with distinct magmatic characteristics: (i) subduction stage, (ii) syntectonic collisional magmatism stage, (iii) post-tectonic collisional magmatism stage and (iv) post-orogenic extension stage.  相似文献   

12.
A Pluri-Gaussian method is developed for facies variables in three dimensions to model vertical cyclicity related to facies ordering and rhythmicity. Cyclicity is generally characterised by shallowing-upward or deepening-upward sequences and rhythmicity by the repetition of facies at constant intervals along sequences. Both of these aspects are commonly observed in shallow-marine carbonate successions, especially in the vertical direction. A grid-free spectral simulation approach is developed, with a separable covariance allowing a dampened hole-effect to capture rhythmicity in the vertical direction and a different covariance in the lateral plane along strata, as in space-time models. In addition, facies ordering is created by using a spatial shift between two latent Gaussian functions in the Pluri-Gaussian approach. Rapid conditioning to data is performed via Gibbs sampling and kriging using the screening properties of separable covariances. The resulting facies transiograms can show complex patterns of cyclicity and rhythmicity. Finally, a three dimensional case study of shallow-marine carbonate deposits at outcrop shows the applicability of the modelling method.  相似文献   

13.
The maximum expected ground motion in Greece is estimated for shallow earthquakes using a deterministic seismic hazard analysis (DSHA). In order to accomplish this analysis the input data include an homogeneous catalogue of earthquakes for the period 426 BC–2003, a seismogenic source model with representative focal mechanisms and a set of velocity models. Because of the discrete character of the earthquake catalogue and of errors in location of single seismic events, a smoothing algorithm is applied to the catalogue of the main shocks to get a spatially smoothed distribution of magnitude. Based on the selected input parameters synthetic seismograms for an upper frequency content of 1 Hz are computed on a grid of 0.2° × 0.2°. The resultant horizontal components for displacement, velocity, acceleration and DGA (Design Ground Acceleration) are mapped. The maps which depict these results cannot be compared with previously published maps based on probabilistic methodologies as the latter were compiled for a mean return period of 476 years. Therefore, in order to validate our deterministic analysis, the final results are compared with PGA estimated from the maximum observed macroseismic intensity in Greece during the period 426 BC–2003.Since the results are obtained for point sources, with the frequency content scaled with moment magnitude, some sensitivity tests are performed to assess the influence of the finite extent of fault related to large events. Sensitivity tests are also performed to investigate the changes in the peak ground motion quantities when varying the crustal velocity models in some seismogenic areas. The ratios and the relative differences between the results obtained using different models are mapped and their mean value computed. The results highlight the importance in the deterministic approach of using good and reliable velocity models.  相似文献   

14.
Three conflicting models are currently proposed for the location and tectonic setting of the Eurasian continental margin and adjacent Tethys ocean in the Balkan region during Mesozoic–Early Tertiary time. Model 1 places the Eurasian margin within the Rhodope zone relatively close to the Moesian platform. A Tethyan oceanic basin was located to the south bordering a large “Serbo-Pelagonian” microcontinent. Model 2 correlates an integral “Serbo-Pelagonian” continental unit with the Eurasian margin and locates the Tethys further southwest. Model 3 envisages the Pelagonian zone and the Serbo-Macedonian zone as conjugate continental units separated by a Tethyan ocean that was sutured in Early Tertiary time to create the Vardar zone of northern Greece and former Yugoslavia. These published alternatives are tested in this paper based on a study of the tectono-stratigraphy of a completely exposed transect located in the Voras Mountains of northernmost Greece. The outcrop extends across the Vardar zone, from the Pelagonian zone in the west to the Serbo-Macedonian zone in the east.Within the Voras Massif, six east-dipping imbricate thrust sheets are recognised. Of these, Units 1–4 correlate with the regional Pelagonian zone in the west (and related Almopias sub-zone). By contrast, Units 5–6 show a contrasting tectono-stratigraphy and correlate with the Paikon Massif and the Serbo-Macedonian zone to the east. These units form a stack of thrust sheets, with Unit 1 at the base and Unit 6 at the top. Unstacking these thrust sheets places ophiolitic units between the Pelagonian zone and the Serbo-Macedonian zone, as in Model 3. Additional implications are, first, that the Paikon Massif cannot be seen as a window of Pelagonian basement, as in Model 1, and, secondly, Jurassic andesitic volcanics of the Paikon Massif locally preserve a gneissose continental basement, ruling out a recently suggested origin as an intra-oceanic arc.We envisage that the Almopias (Vardar) ocean rifted in Triassic time, followed by seafloor spreading. The Almopias ocean was consumed beneath the Serbo-Macedonian margin in Jurassic time, generating subduction-related arc volcanism in the Paikon Massif and related units. Ophiolites were emplaced onto the Pelagonian margin in the west and covered by Late Jurassic (pre-Kimmeridgian) conglomerates. Other ophiolitic rocks formed within the Vardar zone (Ano Garefi ophiolite, Unit 4) in latest Jurassic–Early Cretaceous time and were not deformed until Early Tertiary time. The Vardar zone finally sutured in the Early Tertiary creating the present imbricate thrust structure of the Voras Mountains.  相似文献   

15.
Denghai Bai  Maxwell A. Meju   《Tectonophysics》2003,364(3-4):135-146
Magnetotelluric (MT) geophysical profiling has been applied to the determination of the deep structure of the Longling–Ruili fault (LRF), part of a convergent strike-slip fault system, underneath thick Caenozoic cover in Ruili basin in southwestern Yunnan, China. The recorded MT data have been inverted using a two-dimensional (2-D) nonlinear conjugate gradients scheme with a variety of smooth starting models, and the resulting models show common subsurface conductivity structures that are deemed geological significant. The models show the presence of a conductive (5–60 Ω m) cover sequence that is thickest (1–1.5 km) in the centre of the basin and rapidly pinches out towards the margins. A half-graben structure is interpreted for the Ruili basin. This is underlain by about 7–10 km thick upper crustal layer of high resistivity (>200–4000 Ω m) that is dissected by steep faults, which we interpret to flatten at depth and root into an underlying mid-crustal conductive layer at about 10 km depth. The mid-crustal layer does not appear to have been severely affected by faulting; we interpret it as a zone of partial melt or intracrustal detachment. The MT models suggest SE directed thrusting of basement rocks in the area. The Longling–Ruili fault is interpreted as a NW-dipping feature bounding one of the identified upper crustal fragments underneath Ruili city. We suggest that MT imaging is a potent tool for deep subsurface mapping in this terrain.  相似文献   

16.
The Central India Tectonic Zone (CITZ) is a prominent divide and a major suture zone between the North Indian and South Indian crustal blocks. The resistive upper crust as modeled in the magnetotelluric data from CITZ suggests a dominant tonalite–trdondhjemite–granodiorite composition associated with an accretionary complex characterized by mainly felsic rock components. The highly conductive bodies in this zone might represent mafic/ultramafic-layered intrusives derived from a deeper reservoir of underplated basaltic magma related to the formation of the Cretaceous Deccan flood basalts. The uniformly thick mafic lower crust below the cratons on both sides of the suture is interpreted as the accreted remnants of Archaean and Paleoproterozoic subducted slabs. We redefine the nature of deep faults traversing the CITZ, which were described as steep and penetrating the Moho by previous workers, and classify them as listric faults with gentle dips at depth.Seismic reflection data from the eastern side of the suture suggest a northwestward subduction of the Bhandara Craton. Reflection data from the central part of the CITZ show northerly dip in the southern part suggesting northward subduction of the Dharwar Craton. However, an opposite trend is observed in the northern part of the suture with a southward dip of the Bundelkhand craton. Based on these features, and in conjunction with existing magnetotelluric models, we propose a double-sided subduction history along the CITZ. This would be similar to the ongoing subduction–accretion process in the western Pacific region, which possibly led to the development of paired collision-type and Pacific-type orogens. One important feature is the domal structure along the central part of the suture with a thick felsic crust occurring between mafic and intermediate crust. The high resistivity felsic domain suggests underplated sediments/felsic crust that would have caused the doming. Our model also accounts for the extrusion of regional metamorphic belts at the orogenic core, and the occurrence of high pressure–ultrahigh-temperature paired metamorphic belts within the suture.  相似文献   

17.
This article illustrates the use of linear and nonlinear regression models to obtain quadratic estimates of covariance parameters. These models lead to new insights into the motivation behind estimation methods, the relationships between different methods, and the relationship of covariance estimation to prediction. In particular, we derive the standard estimating equations for minimum norm quadratic unbiased translation invariant estimates (MINQUEs) from an appropriate linear model. Connections between the linear model, minimum variance quadratic unbiased translation invariant estimates (MIVQUEs), and MINQUEs are examined and we provide a minimum norm justification for the use of one-step normal theory maximum likelihood estimates. A nonlinear regression model is used to define MINQUEs for nonlinear covariance structures and obtain REML estimates. Finally, the equivalence of predictions under various models is examined when covariance parameters are estimated. In particular, we establish that when using MINQUE, iterative MINQUE, or restricted maximum likelihood (REML) estimates, the choice between a stationary covariance function and an intrinsically stationary semivariogram is irrelevant to predictions and estimated prediction variances.  相似文献   

18.
Horizontal extension of a previously thickened crust could be the principal mechanism that caused the development of widespread extensional basins throughout the North China block (Hua-Bei region) during the Mesozoic. We develop here a regional tectonic model for the evolution of the lithosphere in the North China block, based on thin sheet models of lithospheric deformation, with numerical solutions obtained using the finite element method. The tectonic evolution of this region is defined conceptually by two stages in our simplified tectonic model: the first stage is dominated by N–S shortening, and the second by E–W extension. We associate the N–S shortening with the Triassic continental collision between the North and South China blocks, assuming that the Tan-Lu Fault system defines the eastern boundary of the North China block. The late Mesozoic E–W extension that created the Mesozoic basin systems requires a change in the regional stress state that could have been triggered by either or both of the following factors: First, gravitational instability of the lithosphere triggered by crustal convergence might have removed the lower layers of the thickened mantle lithosphere and thus caused a rapid increase in the local gravitational potential energy of the lithosphere. Secondly, a change to the constraining stress on the eastern boundary of the North China block, that might have been caused by roll-back of the subducting Pacific slab, could have reduced the E–W horizontal stress enough to activate extension. Our simulations show that widespread thickening of the North China block by as much as 50% can be explained by the collision with South China in the Triassic and Jurassic. If convergence then ceases, E–W extension can occur in the model if the eastern boundary of the region can move outwards. We find that such extension may occur, restoring crustal thickness of order 30 km within a period of 50 Myr or less, if the depth-averaged constitutive relation of the lithosphere is Newtonian, and if the Argand number (the ratio of buoyancy-derived stress to viscous stress) is greater than about 4. Widespread convective thinning of the lithosphere is not required in order to drive the extension with these parameters. If, however, the lithospheric viscosity is non-Newtonian (with strain-rate proportional to the third power of stress) the extensional phase would not occur in a geologically plausible time unless the Argand number were significantly increased by a lithospheric thinning event that was triggered by crustal thickening ratios as low as 1.5.  相似文献   

19.
In studies that involve a finite sample size of spatial data it is often of interest to test (statistically) the assumption that the marginal (or univariate) distribution of the data is Gaussian (normal). This may be important per se because, for example, a data transformation may be desired if the normality hypothesis is rejected, or it may provide a way of testing other hypotheses, such as lognormality, by testing the normality of the logarithms of the observations. The most commonly used tests, such as the Kolmogorov–Smirnov (K–S), chi-square (2), and Shapiro–Wilks (S–W) tests, are designed on the assumption that the observations are independent and identically distributed (iid). In geostatistical applications, however, this is not usually the case unless the spatial covariance (semivariogram) function is a pure nugget variance. If the covariance structure has a (practical) range greater than the minimum distance between observations, the data are correlated and the standard tests cannot be applied to the probability density function (pdf) or cumulative probability function (cdf) estimated directly from the data. The problem with correlated data arises not from the correlation per se but from cases in which correlated data are clustered rather than being located on a regular grid. In these cases inferences requiring iid assumptions may be seriously biased because of the spatial correlation among the observations. If unbiased (i.e., de-clustered) estimates of the pdf or cdf are obtained, then normality tests, such as K-S, 2, or S–W, can be applied using the unbiased estimates and an effective number of samples equivalent to the iid case. There are three questions to be addressed in these cases: Is the distribution ergodic?  相似文献   

20.
Phase equilibria in the system CaO–MgO–SiO2–CO2–H2O–NaCl are calculated to illustrate phase relations in metacarbonates over a wide-range of P–T–X[H2O–CO2–NaCl] conditions. Calculations are performed using the equation of state of Duan et al. (Geochim Cosmochim Acta 59:2869–2882, 1995) for H2O–CO2–NaCl fluids and the internally consistent data set of Gottschalk (Eur J Mineral 9:175–223, 1997) for thermodynamic properties of solids. Results are presented in isothermal-isobarical plots showing stable mineral assemblages as a function of fluid composition. It is shown that in contact-metamorphic P–T regimes the presence of very small concentrations of NaCl in the fluid causes almost all decarbonation reactions to proceed within the two fluid solvus of the H2O–CO2–NaCl system. Substantial flow of magma-derived fluids into marbles has been documented for many contact aureoles by shifts in stable isotope geochemistry of the host rocks and by the progress of volatile-producing mineral reactions controlled by fluid compositions. Time-integrated fluid fluxes have been estimated by combining fluid advection/dispersion models with the spatial arrangement of mineral reactions and isotopic resetting. All existing models assume that minerals react in the presence of a single phase H2O–CO2 fluid and do not allow for the effect that fluid immiscibility has on the flow patterns. It is shown that fluids emanating from calc-alkaline melts that crystallize at shallow depths are brines. Their salinity may vary depending mainly on pressure and fraction of crystallized melt. Infiltration-driven decarbonation reactions in the host rocks inevitably proceed at the boundaries of the two fluid solvus where the produced CO2 is immiscible and may separate from the brine as a low salinity, low density H2O–CO2 fluid. Most parameters of fluid–rock interaction in contact aureoles that are derived from progress of mineral reactions and stable isotope resetting are probably incorrect because fluid phase separation is disregarded.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号