Upscaling permeability for three-dimensional fractured porous rocks with the multiple boundary method

Upscaling permeability of grid blocks is crucial for groundwater models. A novel upscaling method for three-dimensional fractured porous rocks is presented. The objective of the study was to compare this method with the commonly used Oda upscaling method and the volume averaging method. First, the multiple boundary method and its computational framework were defined for three-dimensional stochastic fracture networks. Then, the different upscaling methods were compared for a set of rotated fractures, for tortuous fractures, and for two discrete fracture networks. The results computed by the multiple boundary method are comparable with those of the other two methods and fit best the analytical solution for a set of rotated fractures. The errors in flow rate of the equivalent fracture model decrease when using the multiple boundary method. Furthermore, the errors of the equivalent fracture models increase from well-connected fracture networks to poorly connected ones. Finally, the diagonal components of the equivalent permeability tensors tend to follow a normal or log-normal distribution for the well-connected fracture network model with infinite fracture size. By contrast, they exhibit a power-law distribution for the poorly connected fracture network with multiple scale fractures. The study demonstrates the accuracy and the flexibility of the multiple boundary upscaling concept. This makes it attractive for being incorporated into any existing flow-based upscaling procedures, which helps in reducing the uncertainty of groundwater models.     

Contrasting sources of Late Paleozoic rhyolite magma in the Polish Lowlands: evidence from U–Pb ages and Hf and O isotope composition in zircon

International Journal of Earth Sciences - The Polish Lowlands, located southwest of the Teisseyre–Tornquist Zone, within Trans-European Suture Zone, were affected by bimodal, but dominantly...     

Shale basins,sulfur-deficient ore brines and the formation of exhalative base metal deposits

The massive sulfide deposits of the Iberian Pyrite Belt are interbedded with felsic volcanic rocks and shale, and underlain by several thousand meters of siliciclastic sedimentary rocks known as the PQ Group. Isotope geochemistry and regional geology are both consistent with equilibration of the ore-forming fluids with the PQ Group, prior to ore deposition near the former seafloor. The average Cu:Zn:Pb ratio of the PQ Group rocks (ca. 26:55:19) is similar to the weighted average of all the massive sulfide orebodies combined (ca. 25:52:23).The genetic relationship between massive sulfide deposits and a siliciclastic sedimentary metal source is explained here by a thermodynamic model, proposing that mildly reducing redox conditions imposed by equilibration with the sedimentary rocks are most critical for the formation of an effective ore-forming fluid. Relatively metal-rich but organic-poor pyrite-bearing shale undergoing dewatering of saline pore fluids is an effective source for the generation of sulfur-deficient but relatively iron and base metal-rich brines. Thus, we propose that the giant deposits of the Iberian Pyrite Belt owe their existence not to exceptionally metal-enriched (e.g., magmatic) fluids, but to the existence of a fairly ordinary but large metal source in reactive siliciclastic sediments, combined with an underlying igneous heat source and a particularly efficient mechanism of sulfide precipitation by mixing with H₂S-rich fluids at or near the seafloor.Essentially similar mineral equilibria are imposed when saline fluids are buffered by typical continental basement rocks. Leaching of retrograde minerals and possibly residual salts from their magmatic or metamorphic prehistory is expected to generate similar, variably metal-rich but relatively sulfide-deficient fluids. Thus, the existence of mildly reducing rocks can be the dominant chemical control in the source of fluids generating many volcanogenic, Irish-type or sedex deposits, many of which are known to precipitate their metal load in response to biogenic sulfide addition at the ore deposition site.     

Identifying low-temperature hydrothermal karst and palaeowaters using stable isotopes: a case study from an alpine cave,Entrische Kirche,Austria

The area south of the prominent east–west trending Salzach Valley at the northern rim of the Central Alps of Austria has long been known to host anomalously warm springs emerging from a highly deformed calcite marble (Klammkalk). This unit also hosts cavities whose shapes suggest a hydrothermal karst origin and which are lined by calcite spar. We report here petrographic and isotopic evidence suggesting that dissolution by ascending low-temperature thermal waters also played an important role in the origin of a large cave in this region, Entrische Kirche. A paleo cave wall, preserved behind a thick flowstone in the interior of this cave, revealed a brownish bleaching zone which contrasts to the medium grey colour of the unaltered marble beneath. Across this zone the C and O isotope values gradually decrease by 3 and 11‰, respectively. These compositions are very different from those of the speleothem above but are similar to phreatic calcite spar from hydrothermal karst cavities in other outcrops in the area, where the absence of two-phase fluid inclusions suggests a low-temperature (less than ca. 50°C) hydrothermal origin. U/Th dating of the flowstone capping the alteration zone yielded a minimum age of the thermal water invasion in Entrische Kirche of ca. 240 kyr. There is no evidence in Entrische Kirche that these palaeowaters reached the point of calcite precipitation, but it is physically conceivable that higher and as yet unexplored parts of this deep (ca. 900 m) cave contain cavities lined by phreatic cave spar.     

Mapping the hydraulic connection between a coalbed and adjacent aquifer: example of the coal-seam gas resource area,north Galilee Basin,Australia

Coal-seam gas production requires groundwater extraction from coal-bearing formations to reduce the hydraulic pressure and improve gas recovery. In layered sedimentary basins, the coalbeds are often separated from freshwater aquifers by low-permeability aquitards. However, hydraulic connection between the coalbed and aquifers is possible due to the heterogeneity in the aquitard such as the existence of conductive faults or sandy channel deposits. For coal-seam gas extraction operations, it is desirable to identify areas in a basin where the probability of hydraulic connection between the coalbed and aquifers is low in order to avoid unnecessary loss of groundwater from aquifers and gas production problems. A connection indicator, the groundwater age indictor (GAI), is proposed, to quantify the degree of hydraulic connection. The spatial distribution of GAI can indicate the optimum positions for gas/water extraction in the coalbed. Depressurizing the coalbed at locations with a low GAI would result in little or no interaction with the aquifer when compared to the other positions. The concept of GAI is validated on synthetic cases and is then applied to the north Galilee Basin, Australia, to assess the degree of hydraulic connection between the Aramac Coal Measure and the water-bearing formations in the Great Artesian Basin, which are separated by an aquitard, the Betts Creek Beds. It is found that the GAI is higher in the western part of the basin, indicating a higher risk to depressurization of the coalbed in this region due to the strong hydraulic connection between the coalbed and the overlying aquifer.     

Flood frequencies reveal Holocene rapid climate changes (Lower Moulouya River,northeastern Morocco)

Current high‐resolution palaeoenvironmental records reveal short‐term Holocene coolings. One of these major Holocene rapid climate changes occurred between 3.2 and 2.7 cal. ka BP. The sensitivity of river systems vis‐à‐vis slight and short‐term Holocene climatic variations is a subject of controversy in the scientific community. In this paper, we present a 4.0 to 1.4 cal ka BP palaeoflood record from the Lower Moulouya River (northeastern Morocco) to demonstrate the high sensitivity of semiarid rivers in the southwestern Mediterranean towards Holocene environmental changes. The Lower Moulouya flood deposits are characterised by thick, well‐stratified, predominantly clayey to silty overbank fine sediments. These cohesive sediments show evidence of excellent preservation conditions against fluvial erosion and contain a continuous record of mid to late Holocene flood sequences. The Moulouya palaeoflood record can be interpreted in the context of regional and global high‐resolution proxy data, revealing a strong coupling with Holocene rapid climate changes. The centennial‐scale Moulouya palaeohydrological history will be discussed with palaeoenvironmental data from the same record (palaeomagnetics, sedimentary charcoal record, anthracological analyses, snail analyses) to generate new ideas about the mid to late Holocene hydrological cycle in the southwestern Mediterranean. The deduced features of pronounced Lower Moulouya flooding and the decreased fire recurrence during Holocene cooling remain somewhat inconsistent with the interpretation of other palaeohydrological and paleaoecological records from the southwestern Mediterranean. However, enhanced Lower Moulouya flood frequencies between 3.2 and 2.7 cal. ka BP agree with increased floodplain aggradation in other major river systems of Mediterranean North Africa. Copyright © 2009 John Wiley & Sons, Ltd.     

Fluid and source magma evolution of the Questa porphyry Mo deposit,New Mexico,USA

Combined fluid inclusion microthermometry and microanalysis by laser ablation inductively coupled plasma mass spectrometry (LA-ICPMS) are used to constrain the hydrothermal processes forming a typical Climax-type porphyry Mo deposit. Molybdenum mineralisation at Questa occurred in two superimposed hydrothermal stages, a magmatic-hydrothermal breccia and later stockwork veining. In both stages, texturally earliest fluids were single-phase, of low salinity (~7 wt.% NaCl_equiv.) and intermediate-density. Upon decompression to ~300 bar, they boiled off a vapour phase, leaving behind a residual brine (up to 45 wt.% NaCl_equiv) at temperatures of ~420°C. The highest average Mo concentrations in this hot brine were ~500 μg/g, exceeding the Mo content of the intermediate-density input fluid by about an order of magnitude and reflecting pre-concentration of Mo by fluid phase separation prior to MoS₂ deposition from the brine. Molybdenum concentrations in brine inclusions, then, decrease down to 5 μg/g, recording Mo precipitation in response to cooling of the saline liquid to ~360°C. Molybdenite precipitation from a dense, residual and probably sulphide-depleted brine is proposed to explain the tabular shape of the ore body and the absence of Cu-Fe sulphides in contrast to the more common Cu-Mo deposits related to porphyry stocks. Cesium and Rb concentrations in the single-phase fluids of the breccia range from 2 to 8 and from 40 to 65 μg/g, respectively. In the stockwork veins, Cs and Rb concentrations are significantly higher (45–90 and 110–230 μg/g, respectively). Because Cs and Rb are incompatible and hydrothermally non-reactive elements, the systematic increase in their concentration requires two distinct pulses of fluid exsolution from a progressively more fractionated magma. By contrast, major element and ore metal concentrations of these two fluid pulses remain essentially constant. Mass balance calculations using fluid chemical data from LA-ICPMS suggest that at least 25 km³ of melt and 7 Gt of deep input fluid were necessary to provide the amount of Mo contained in the stockwork vein stage alone. While the absolute amounts of fluid and melt are uncertain, the well-constrained element ratios in the fluids together with empirical fluid/melt partition coefficients derived from the inclusion analyses suggest a high water content of the source melt of ~10%. In line with other circumstantial evidence, these results suggest that initial fluid exsolution may have occurred at a confining pressure exceeding 5 kbar. The source of the molybdenum-mineralising fluids probably was a particularly large magma chamber that crystallised and fractionated in the lower crust or at mid-crustal level, well below the shallow intrusions immediately underlying Questa and other porphyry molybdenum deposits. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Neoproterozoic calc‐alkaline peraluminous granitoids of the Deleihimmi pluton,Central Eastern Desert,Egypt: implications for transition from late‐ to post‐collisional tectonomagmatic evolution in the northern Arabian–Nubian Shield

The widely distributed late‐collisional calc‐alkaline granitoids in the northern Arabian–Nubian Shield (ANS) have a geodynamic interest as they represent significant addition of material into the ANS juvenile crust in a short time interval (∼630–590 Ma). The Deleihimmi granitoids in the Egyptian Central Eastern Desert are, therefore, particularly interesting since they form a multiphase pluton composed largely of late‐collisional biotite granitoids enclosing granodiorite microgranular enclaves and intruded by leuco‐ and muscovite granites. Geochemically, different granitoid phases share some features and distinctly vary in others. They display slightly peraluminous (ASI = 1–1.16), non‐alkaline (calc‐alkaline and highly fractionated calc‐alkaline), I‐type affinities. Both biotite granitoids and leucogranites show similar rare earth element (REE) patterns [(La/Lu)_N = 3.04–2.92 and 1.9–1.14; Eu/Eu* = 0.26–0.19 and 0.11–0.08, respectively) and related most likely by closed system crystal fractionation of a common parent. On the other hand, the late phase muscovite granites have distinctive geochemical features typical of rare‐metal granites. They are remarkably depleted in Sr and Ba (4–35 and 13–18 ppm, respectively), and enriched in Rb (381–473 ppm) and many rare metals. Moreover, their REE patterns show a tetrad effect (TE_1,3 = 1.13 and 1.29) and pronounced negative Eu anomalies (Eu/Eu* = 0.07 and 0.08), implying extensive open system fractionation via fluid–rock interaction during the magmatic stage. Origin of the calc‐alkaline granitoids by high degree of partial melting of mafic lower crust with subsequent crystal fractionation is advocated. The broad distribution of late‐collisional calc‐alkaline granitoids in the northern ANS is related most likely to large areal and intensive lithospheric delamination subsequent to slab break‐off and crustal/mantle thickening. Such delamination caused both crustal uplift and partial melting of the remaining mantle lithosphere in response to asthenospheric uprise. The melts produced underplate the lower crust to promote its melting. The presence of microgranular enclaves, resulting from mingling of mantle‐derived mafic magma with felsic crustal‐derived liquid, favours this process. The derivation of the late‐phase rare‐metal granites by open system fractionation via fluid interaction is almost related to the onset of extension above the rising asthenosphere that results in mantle degassing during the switch to post‐collisional stage. Consequently, the switch from late‐ to post‐collisional stage of crustal evolution in the northern ANS could be potentially significant not only geodynamically but also economically. Copyright © 2011 John Wiley & Sons, Ltd.     

Probability Density Functions as Botanical-Climatological Transfer Functions for Climate Reconstruction

We present a new procedure, the pdf method (pdf=probability density function), for reconstructing Quaternary climate utilizing botanical data. The procedure includes the advantages of the indicator species method by considering the fossil and modern presence and absence of taxa rather than their frequencies, thus avoiding the need for modern analog plant communities. Overcoming the problematic use of absolute limits to describe climate response ranges is the main progress of the pdf method in comparison to the indicator species method. This advantage results from estimating probability density functions (pdfs) for monthly mean January and July temperature conditional on the present day occurrence of single taxa. Gaussian distributions sufficiently approximate pdfs of many, although not all, studied taxa. On the assumption of statistical independence, the procedure calculates a joint pdf as the product of the pdfs of the individual taxa. This algorithm weights each taxon according to the extent of its climate response range expressed by its covariance structure. We interpret the maximum of the resulting pdf as the most likely climate and its confidence interval as the uncertainty range. To avoid an artificial reduction of uncertainty arising from the use of numerous similar pdfs, a preselection method is proposed based on the Mahalanobis distance between pdfs. The pdf method was applied to the Carpinus phase of a profile from Gröbern, Germany, that spans the last interglaciation (Eemian). The reconstructed most probable January and July temperatures of about 0.0°C and 18.4°C barely differ from the modern values of −0.5°C and 18.3°C.