Fluid‐fluxed melting and melt loss in a syntectonic contact metamorphic aureole from the Variscan eastern Pyrenees          下载免费PDF全文

C. Aguilar  M. Liesa  J. Reche  R. Powell 《Journal of Metamorphic Geology》2016,34(4):379-400

Open‐system behaviour through fluid influx and melt loss can produce a variety of migmatite morphologies and mineral assemblages from the same protolith composition. This is shown by different types of granulite facies migmatite from the contact aureole of the Ceret gabbro–diorite stock in the Roc de Frausa Massif (eastern Pyrenees). Patch, stromatic and schollen migmatites are identified in the inner contact aureole, whereas schollen migmatites and residual melanosomes are found as xenoliths inside the gabbro–diorite. Patch and schollen migmatites record D1 and D2 structures in folded melanosome and mostly preserve the high‐T D2 in granular or weakly foliated leucosome. Stromatic migmatites and residual melanosomes only preserve D2. The assemblage quartz–garnet–biotite–sillimanite–cordierite±K‐feldspar–plagioclase is present in patch and schollen migmatites, whereas stromatic migmatites and residual melanosomes contain a sub‐assemblage with no sillimanite and/or K‐feldspar. A decrease in X _Fe (molar Fe/(Fe + Mg)) in garnet, biotite and cordierite is observed from patch migmatites through schollen and stromatic migmatites to residual melanosomes. Whole‐rock compositions of patch, schollen and stromatic migmatites are similar to those of non‐migmatitic rocks from the surrounding area. These metasedimentary rocks are interpreted as the protoliths of the migmatites. A decrease in the silica content of migmatites from 63 to 40 wt% SiO₂ is accompanied by an increase in Al₂O₃ and MgO+FeO and by a depletion in alkalis. Thermodynamic modelling in the NCKFMASHTO system for the different types of migmatite provides peak metamorphic conditions ~7–8 kbar and 840 °C. A nearly isothermal decompression history down to 5.5 kbar was followed by isobaric cooling from 840 °C through 690 °C to lower temperatures. The preservation of granulite facies assemblages and the variation in mineral assemblages and chemical composition can be modelled by ongoing H₂O‐fluxed melting accompanied by melt loss. The fluids were probably released by the crystallizing gabbro–diorite, infiltrating the metasedimentary rocks and fluxing melting. Release of fluids and melt loss were probably favoured by coeval deformation (D2). The amount of melt remaining in the system varied considerably among the different types of migmatite. The whole‐rock compositions of the samples, the modelled compositions of melts at the solidus at 5.5 kbar and the residues show a good correlation.  相似文献   

Numerical modelling of orogenic processes and gold mineralisation in the southeastern part of the Yilgarn Craton,Western Australia     

P. Sorjonen-Ward  Y. Zhang  C. Zhao 《Australian Journal of Earth Sciences》2013,60(6):935-964

We use numerical modelling codes to simulate aspects of some current hypotheses for the origin of gold deposits and hydrothermal systems in the Yilgarn Craton of Western Australia. In particular, we investigate conceptual models advocating vertically continuous hydrothermal systems as well as those invoking extensive lateral flow and possible links with advection of heat by late orogenic granitic magmatism. Numerical models of part of the Eastern Goldfields Province and Southern Cross Province have been built with FLAC3D, to simulate crustal‐scale coupled interaction between deformation and fluid flow. These illustrate the potential for fluid focusing and mixing in shear zones, including downflow of meteoric water, lateral fluid flow driven by topographic elevation and upwards flow of fluids derived from melting and metamorphism in the deep crust. In some cases, downflow also occurs within the middle crust, at depths where fluid influx might trigger melting if the geothermal gradient were appropriate. The models indicate that tectonic wedging within a layered crust and diverging thrust systems that generate ‘pop‐up’ wedges may be important in facilitating efficient fluid upflow and downflow during uplift, while topographic elevation related to asymmetric thrust migration and loading tends to promote lateral fluid flow. However, the effect of topography appears more important than the precise depth or location of the site of fluid production in the deep crust. The effects of thermal convection and fluid‐fluid interaction have also been numerically modelled for a simplified section across the Kalgoorlie Terrane. Modelling under both hydrostatic and lithostatically overpressured pore‐pressure gradients has effectively delineated domains of convective fluid flow within the middle and upper crust, and has identified two generic sites that are favourable for fluid mixing, notably hangingwall and footwall environments in major shear zones, such as the Bardoc Shear, and in broad antiforms, such as the Goongarrie ‐ Mt Pleasant Antiform. The thermal effect of small plutons embedded in a regional metamorphic regime can cause significant lateral displacement of fluid convection patterns, over distances greater than pluton diameter, as well as more proximal effects on precipitation and dissolution of mineral species. However, these results are highly dependent on the pore‐pressure gradient and the permeability structure of the crust, and require magmatic and metamorphic fluid generation to be precisely timed with respect to deformation, thus reinforcing the dynamic feedback between deformation, magmatism and fluid production and migration.  相似文献   

Vegetation response to rainfall seasonality and interannual variability in tropical dry forests          下载免费PDF全文

Rodolfo Souza  Xue Feng  Antonio Antonino  Suzana Montenegro  Eduardo Souza  Amilcare Porporato 《水文研究》2016,30(20):3583-3595

Projected changes in rainfall seasonality and interannual variability are expected to have severe impacts on arid and semi‐arid tropical vegetation, which is characterized by a fine‐tuned adaptation to extreme rainfall seasonality. To study the response of these ecosystems and the related changes in hydrological processes to changes in the amount and seasonality of rainfall, we focused on the caatinga biome, the typical seasonally dry forest in semi‐arid Northeast Brazil. We selected four sites across a gradient of rainfall amount and seasonality and analysed daily rainfall and biweekly Normalized Difference Vegetation Index (NDVI) data for hydrological years 2000 to 2014. Rainfall seasonal and interannual statistics were characterized by recently proposed metrics describing duration, timing and intensity of the wet season and compared to similar metrics of NDVI time series. The results show that the caatinga tends to have a more stable response with longer and less variable growing seasons (3.1 ± 0.1 months) compared to the duration wet seasons (2.0 ± 0.5 months). The ecosystem ability to buffer the interannual variability of rainfall is also evidenced by the stability in the timing of the growing season compared to the wet season, which results in variable delays (ranging from 0 to 2 months) between the peak of the rainfall season and the production of leaves by the ecosystem. The analyses show that the shape and size of the related hysteresis loops in the rainfall–NDVI relations are linked to the buffering effects of soil moisture and plant growth dynamics. Finally, model projections of vegetation response to different rainfall scenarios reveal the existence of a maximum in ecosystem productivity at intermediate levels of rainfall seasonality, suggesting a possible trade‐off in the effects of intensity (i.e. amount) and duration of the wet season on vegetation growth and related soil moisture dynamics and transpiration rates. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Accounting for soil parameter uncertainty in a physically based and distributed approach for rainfall‐triggered landslides          下载免费PDF全文

E. Arnone  Y. G. Dialynas  L. V. Noto  R. L. Bras 《水文研究》2016,30(6):927-944

In this study we propose a probabilistic approach for coupled distributed hydrological‐hillslope stability models that accounts for soil parameters uncertainty at basin scale. The geotechnical and soil retention curve parameters are treated as random variables across the basin and theoretical probability distributions of the Factor of Safety (FS) are estimated. The derived distributions are used to obtain the spatio‐temporal dynamics of probability of failure, in terms of parameters uncertainty, conditioned to soil moisture dynamics. The framework has been implemented in the tRIBS‐VEGGIE (Triangulated Irregular Network (TIN)‐based Real‐time Integrated Basin Simulator‐VEGetation Generator for Interactive Evolution)‐Landslide model and applied to a basin in the Luquillo Experimental Forest (Puerto Rico) where shallow landslides are common. In particular, the methodology was used to evaluate how the spatial and temporal patterns of precipitation, whose variability is significant over the basin, affect the distribution of probability of failure, through event scale analyses. Results indicate that hyetographs where heavy precipitation is near the end of the event lead to the most critical conditions in terms of probability of failure. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Timing of the younger dryas glacial maximum in western Norway     

Øystein S. Lohne  Jan Mangerud  John Inge Svendsen 《第四纪科学杂志》2012,27(1):81-88

This study precisely constrains the timing of the Younger Dryas (YD) glacial maximum in south‐western Norway by utilizing sediment records from lake basins. Two of the basins, located on the distal side of the mapped Herdla–Halsnøy Moraine, received meltwater directly from the ice sheet only when the ice margin reached its maximum extent during the YD. In the cores, the ice maximum is represented by well‐defined units with meltwater deposits, dominantly laminated silt. Plant macrofossils in the sediment sequences are common and we obtained 18 radiocarbon ages from one of the cores. By applying Bayesian age–depth modelling we obtained a precise date for this meltwater event and thereby also for the timing of the YD glacial maximum. We conclude that the ice‐sheet advance culminated at the Halsnøy Moraine at 11 760 ± 120 cal a BP, and that the ice margin stayed in this position for 170 ± 120 years. The subsequent retreat started at 11 590 ± 100 cal a BP, i.e. close to the YD/Holocene boundary. Withdrawal was probably triggered by abrupt climatic warming at this time. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Estimation of instantaneous peak flows from maximum mean daily flows using the HBV hydrological model          下载免费PDF全文

J. Ding  M. Wallner  H. Müller  U. Haberlandt 《水文研究》2016,30(9):1431-1448

The record length and quality of instantaneous peak flows (IPFs) have a great influence on flood design, but these high resolution flow data are not always available. The primary aim of this study is to compare different strategies to derive frequency distributions of IPFs using the Hydrologiska Byråns Vattenbalansavdelning (HBV) hydrologic model. The model is operated on a daily and an hourly time step for 18 catchments in the Aller‐Leine basin, Germany. Subsequently, general extreme value (GEV) distributions are fitted to the simulated annual series of daily and hourly extreme flows. The resulting maximum mean daily flow (MDF) quantiles from daily simulations are transferred into IPF quantiles using a multiple regression model, which enables a direct comparison with the simulated hourly quantiles. As long climate records with a high temporal resolution are not available, the hourly simulations require a disaggregation of the daily rainfall. Additionally, two calibrations strategies are applied: (1) a calibration on flow statistics; (2) a calibration on hydrographs. The results show that: (1) the multiple regression model is capable of predicting IPFs with the simulated MDFs; (2) both daily simulations with post‐correction of flows and hourly simulations with pre‐processing of precipitation enable a reasonable estimation of IPFs; (3) the best results are achieved using disaggregated rainfall for hourly modelling with calibration on flow statistics; and (4) if the IPF observations are not sufficient for model calibration on flow statistics, the transfer of MDFs via multiple regressions is a good alternative for estimating IPFs. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Identification of appropriate lags and temporal resolutions for low flow indicators in the River Rhine to forecast low flows with different lead times     

Mehmet C. Demirel  Martijn J. Booij  Arjen Y. Hoekstra 《水文研究》2013,27(19):2742-2758

The aim of this paper is to assess the relative importance of low flow indicators for the River Rhine and to identify their appropriate temporal lag and resolution. This is done in the context of low flow forecasting with lead times of 14 and 90 days. First, the Rhine basin is subdivided into seven sub‐basins. By considering the dominant processes in the sub‐basins, five low flow indicators were selected: precipitation, potential evapotranspiration, groundwater storage, snow storage and lake storage. Correlation analysis was then carried out to determine the relationship between observed low flows and preselected indicators with varying lags (days) and temporal resolutions (from 1 day to 7 months). The results show that the most important low flow indicators in the Alpine sub‐basins for forecasts with a lead time of 14 days are potential evapotranspiration with a large lag and temporal resolution, and lake levels with a small lag and temporal resolution. In the other sub‐basins groundwater levels with a small lag and temporal resolution are important in addition to potential evapotranspiration with a large lag and temporal resolution. The picture is slightly different for forecasts with a lead time of 90 days. The snow storage in the Alpine sub‐basins and the precipitation in the other sub‐basins also become relevant for low flows. Consequently, the most important low flow indicators in the Alpine sub‐basins for forecasts with a lead time of 90 days are potential evapotranspiration with a large lag and temporal resolution, lake levels with a small lag and temporal resolution and snow storage with a small lag and large temporal resolution. The resultant correlation maps provide appropriate lags and temporal resolutions for indicators to forecast low flows in the River Rhine with different lead times. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

Estimation of winter leaf area index and sky view fraction for snow modelling in boreal coniferous forests: consequences on snow mass and energy balance     

Sirpa Rasmus  David Gustafsson  Harri Koivusalo  Ari Laurén  Achim Grelle  Olli‐Kalle Kauppinen  Ola Lagnvall  Anders Lindroth  Kai Rasmus  Magnus Svensson  Per Weslien 《水文研究》2013,27(20):2876-2891

Leaf area index (LAI) and canopy coverage are important parameters when modelling snow process in coniferous forests, controlling interception and transmitting radiation. Estimates of LAI and sky view factor show large variability depending on the estimation method used, and it is not clear how this is reflected in the calculated snow processes beneath the canopy. In this study, the winter LAI and sky view fraction were estimated using different optical and biomass‐based approximations in several boreal coniferous forest stands in Fennoscandia with different stand density, age and site latitude. The biomass‐based estimate of LAI derived from forest inventory data was close to the values derived from the optical measurements at most sites, suggesting that forest inventory data can be used as input to snow hydrological modelling. Heterogeneity of tree species and site fertility, as well as edge effects between different forest compartments, caused differences in the LAI estimates at some sites. A snow energy and mass balance model (SNOWPACK) was applied to detect how the differences in the estimated values of the winter LAI and sky view fraction were reflected in simulated snow processes. In the simulations, an increase in LAI and a decrease in sky view fraction changed the snow surface energy balance by decreasing shortwave radiation input and increasing longwave radiation input. Changes in LAI and sky view fraction affected directly snow accumulation through altered throughfall fraction and indirectly snowmelt through the changed surface energy balance. Changes in LAI and sky view fraction had a greater impact on mean incoming radiation beneath the canopy than on other energy fluxes. Snowmelt was affected more than snow accumulation. The effect of canopy parameters on evaporation loss from intercepted snow was comparable with the effect of variation in governing meteorological variables such as precipitation intensity and air temperature. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

Determination of synthetic wind functions for estimating open water evaporation with Computational Fluid Dynamics     

P. Vidal‐López  V. Martínez‐Alvarez  B. Gallego‐Elvira  B. Martín‐Górriz 《水文研究》2012,26(25):3945-3952

A computational fluid dynamics (CFD)‐based methodology is proposed to derive convective mass‐transfer coefficients (wind functions) that are required for estimating evaporation of water bodies with the mass‐transfer method. Three‐dimensional CFD was applied to model heat transfer in two water bodies: a Class‐A tank evaporimeter and an on‐farm artificial pond. The standard k–? model assuming isotropic turbulence was adopted to describe turbulent heat transport, whereas the heat and mass transfer analogy was assumed to derive the wind functions. The CFD‐derived wind functions were very similar to those empirically derived from the experimental water bodies. The evaporation rates calculated with the synthetic wind functions were in good agreement with hourly and daily evaporation measurements for the tank and pond, respectively. The proposed CFD‐approach is generalisable and cost effective, because it has low input data requirements. Besides, it provides additional capability of modelling the spatial distribution of the evaporation rate over the water surface. Although the application of CFD to water bodies evaporation modelling is still in development, it looks very promising. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

Fibre reinforced sands: from experiments to modelling and beyond     

A. Diambra  E. Ibraim  A. R. Russell  D. Muir Wood 《国际地质力学数值与分析法杂志》2013,37(15):2427-2455

Based on hypotheses derived directly from experimental observations of the triaxial behaviour, a constitutive model for fibre reinforced sands is built in this paper. Both the sand matrix and the fibres obey their own constitutive law, whereas their contributions are superimposed using a volumetric homogenization procedure. The Severn‐Trent sand model, which combines well‐known concepts such as critical state theory, Mohr‐Coulomb like strength criterion, bounding surface plasticity and kinematic hardening, is adopted for the sand matrix. Although the fibres are treated as discrete forces with defined orientation, an equivalent continuum stress for the fibre phase is derived to allow the superposition of effects of sand and fibres. The fibres are considered as purely tensile elements following a linear elastic constitutive rule. The strain in the fibres is expressed as a fraction of the strain in the reinforced sample so that imperfect bonding is assumed at the sand‐fibre interface. Only those fibres oriented within the tensile strain domain of the sample can mobilize tensile stress—the orientation of fibres is one of the key ingredients to capture the anisotropic behaviour of fibre reinforced soil that is observed for triaxial compression and extension loading. A further mechanism of partition of the volume of voids between the fibres and the sand matrix is introduced and shown to be fundamental for the simulation of the volumetric behaviour of fibre‐reinforced soils. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献