Crustal Reservoirs Research Articles

Petrological and geochemical evidence for a hot crystallization path and a recharge filtering bypass at Antimilos, Milos volcanic field, Greece

Antimilos volcano in the South Aegean Volcanic Arc, Greece, comprises an andesite–dacite suite that follows a distinct evolutionary path than the main edifice of the Milos volcanic field, despite their proximity. Petrographic and geochemical analyses reveal that basaltic andesite to low-Si dacite lavas have similar phenocryst assemblages that indicate crystallization from hot, relatively dry magmas in an upper crustal storage region. Rare antecrystic high-Mg# clinopyroxene cores with low Y, low Dy, and high Sr contents record the cryptic involvement of amphibole, a phase nominally absent from the erupted products, in the deeper parts of the plumbing system. Low temperature antecrysts with textures recording various degrees of disequilibrium suggest a protracted history of interaction between the upper crustal reservoir and deeper mafic melts, forming mobile hybrid magmas that consequently erupt as highly mingled, crystal-rich lava domes. Antimilos magmas seem to have escaped recharge filtering in the upper crust and prolonged stalling, which is the process that is probably responsible for the paucity of mafic eruptions in the rest of the Milos volcanic system. Large extensional structures offshore of Antimilos promote rapid ascent of mafic melts, inhibiting prolonged stalling and interaction with the arc crust. This model highlights the dominant role of the regional stress field in generating petrologically distinct suites in the marginal parts of some volcanic fields.

Identifying crustal contributions in the Patagonian Chon Aike Silicic Large Igneous Province

The volcanic rocks of the Chon Aike Silicic Large Igneous Province (CASP) are recognized as magmas dominantly produced by crustal anatexis. Investigating the zircon of the CASP provides an opportunity to gain further insight into geochemical and isotopic differences of the potential magmatic sources (i.e., crust versus mantle), to identify crustal reservoirs that contributed to the felsic magmas during anatexis, and to quantify the contributions of the respective sources. We present a combined zircon oxygen and hafnium isotope and trace element dataset for 16 volcanic units of the two youngest volcanic phases in Patagonia, dated here with LA-ICP-MS U–Pb geochronology at ca. 148–153 Ma (El Quemado Complex, EQC) and ca. 159 Ma (western Chon Aike Formation, WCA). The EQC zircon have 18O-enriched values (δ18O from 7 to 9.5‰) with correspondingly negative initial εHf values (− 2.0 to − 8.0). The WCA zircon have δ18O values between 6 and 7‰ and εHf values ranging between − 4.0 and + 1.5. Binary δ18O-εHf mixing models require an average of 70 and 60% melt derived from partial melting of isotopically distinct metasedimentary basements for the EQC and WCA, respectively. Zircon trace element compositions are consistent with anatexis of sedimentary protoliths derived from LIL-depleted upper continental crustal sources. The overlap between a high heat flux environment (i.e., widespread extension and lithospheric thinning) during supercontinental breakup and a fertile metasedimentary crust was key in producing voluminous felsic volcanism via anatexis following the injection and emplacement of basaltic magmas into the lower crust.

Emergence of the 2.1 Ga Francevillian biota was preceded by unprecedented hydrothermally driven seawater eutrophication

Recently, two independent studies suggest that the emergence of putative fossilized macro-eukaryotes in the Paleoproterozoic Francevillian Basin, ~2.1 billion years ago, may be related to a rise in seawater Zn bioavailability. This explanation is reliant on their extraordinary high Zn content and association with light Zn isotopes characteristic of eukaryotic enrichment. However, the trigger and origin of rising seawater Zn supply to the basin remains unknown. This study unravels a transient episode of intense submarine hydrothermal activity that triggered the weathering of a nutrient-rich oceanic crust reservoir, related to the collision of the Congo-São Francisco cratons during the Eburnean-Transamazonian orogeny, as the source of abundant seawater dissolved Zn, together with a suite of essential trace metals and phosphate to the continental margin waters. Surprisingly, the initiation of hydrothermal weathering coincided with the rapid onset of a rare Paleoproterozoic seawater eutrophication event. This transition is marked by basin-wide redox stratification, high sediment loading with organic carbon (Corg) and nitrogen, elevated C/N ratios, a steep negative Corg and positive bulk N isotope excursion, positive Ce anomalies, and low Mn/Fe ratios. Importantly, the transient eutrophication event ended with a reversal to lower seawater phosphate levels that coincided with rapid seawater ventilation and the appearance of macrofossil bearing sediments in Franceville. We suggest that these unexpected, localized conditions, set the stage for the emergence of the Francevillian biota.

Magma plumbing architectures and timescales of magmatic processes during historical magmatism on the Reykjanes Peninsula, Iceland

In the past 4000 years, three rifting episodes have occurred on the Reykjanes Peninsula (RP, SW Iceland) at time intervals of about 800-1000 years, with the activation of four of the five volcanic systems on the peninsula: Reykjanes, Svartsengi, Krýsuvík and Brennisteinsfjöll. In 2021 AD, 2022 AD and 2023 AD, three eruptions occurred in the Fagradalsfjall volcanic system, the only volcanic system that remained inactive during these past rifting episodes. Considering the eruptive history of the peninsula, this could mark the onset of a new eruptive episode. Although extensive petrochemical datasets have emerged from past rifting episodes, little is known about magma plumbing configuration and timescales of pre-eruptive magmatic processes. In this work, we focus on 16 basaltic lava units erupted across RP during a well-established rifting episode: the 800-1240 AD Fires. We analysed major and minor element contents of glasses, mineral phases and melt inclusions, and we reconstruct plumbing system architecture and constrain mixing-to-eruption timescales, with the goal of assessing likely eruption scenarios of a possible future eruptive episode. Mineral-melt equilibria suggests that most of crystal cores are antecrysts that did not crystallize from the magmas that carried them to the surface, but instead were scavenged from disaggregating mush layers genetically related to the erupted magmas. Magmas erupted in the easternmost part of the RP reveal larger contributions from geochemically enriched mantle-derived melts. However, MIs that record various extents of mixing of geochemically enriched and depleted components throughout the entire peninsula appear devoid of any spatial control of the underlying mantle. Independent clinopyroxene-melt and melt-based barometry, after correction for microlite crystallization, show consistent results and suggest that magmas from the western part of the peninsula were extracted from reservoirs located at about 7-10 km depth. In contrast, and similar to observations from recent eruptions at Fagradalsfjall, the easternmost system, Brennisteinsfjöll, was fed from deeper crustal reservoirs, at about 14-21 km depth. Fe-Mg diffusion chronometry applied to olivine rims suggests short mixing-to-eruption timescales on the order of days to weeks and likely faster ascent rates of magmas underneath Brennisteinsfjöll. Compared to 2021 AD Fagradalsfjall, mixing-to-eruption timescales were considerably shorter during the 800-1240 AD Fires. With inferred eruption durations on the order of months, the short olivine rim diffusion records are consistent with several discrete eruptive phases, which eventually assembled the individual magmatic units. In the light of possible future eruptive episodes on the RP, our work suggests that magmas are likely to accumulate at persistent depths across the RP, with key implication for interpreting pre-eruptive signals and eruption mitigation.

Hafnium isotope systematics of zircon in high-grade metamorphic rocks of the Anabar shield, Siberia: Radiogenic Hf without mantle input?

The timing for the onset of plate tectonics, along with the secular changes in the tectonic settings of continental crust formation, continue to be debated. Recent interpretations based on the increase in zircon 176Hf/177Hf ratios at the time of crystallisation (expressed as εHf(t) with respect to chondritic evolution) have been used to ascertain changes in geodynamic settings in the early Earth, specifically in the 3.8–3.6 Ga interval. This increase is widely interpreted as a change in magma generation, from source(s) dominated by ancient crust to source(s) dominated by juvenile inputs from the mantle. At issue, Hf isotope variations remain limited in the early Earth due to the long decay of the LuHf system. This feature, along with the scarcity of rocks and minerals of Eo/Mesoarchaean and Hadean ages, generate large uncertainties over the nature and the timing of the interactions between mantle and crustal reservoirs. The distinction between mantle and crustal sources becomes much clearer in the Palaeoproterozoic, and the study of ancient terranes with several billion years of protracted crustal evolution may hold the key to unlock complex mantle–crust interactions. Here we investigate high-grade metamorphic rocks from the Anabar shield, which contain zircons that crystallised between the Eoarchaean (oldest core at 3814 ± 16 Ma) and the Palaeoproterozoic (youngest core at 2251 ± 15 Ma and youngest rim at 1910 ± 21 Ma). The combination of in situ UPb and Hf isotope analyses in zircon indicates the formation of the continental crust in the Siberian Craton in the Eoarchaean, and a conspicuous metamorphic event at 2.0–1.9 Ga. We demonstrate that 2.0–1.9 Ga zircon ages reflect recrystallisation processes under subsolidus conditions, involving the breakdown of high-Lu/Hf minerals (i.e. garnet and pyroxene). The εHf(t) shift at 2.0–1.9 Ga towards more radiogenic values may not be related to a change in magmatic style and sources, but rather to resetting of the LuHf system in response to heating and metamorphic reactions on a mineral scale. Our findings challenge the widely-evoked mechanism of changes in tectonic style and magma sources to account for vertical arrays in the εHf(t) versus crystallisation age space. This calls for considering alternative options, including those based on petrographic data, when interpreting Hf isotope variations in the Hadean/Archaean detrital zircon record.

Magma process of the ore-forming gabbrodiorite porphyry of the Meishan iron deposit, eastern China: Insights from the phenocrysts chemistry

Mesozoic iron deposits in the Ningwu volcanic basin in eastern China are genetically linked to the gabbrodiorite porphyries, which have long been regarded as the subvolcanic intrusive of the andesitic volcanics. Compared with abundant research conducted on the mineralization process, the magmatic process of the ore-forming intrusion has been poorly understood. This study performed petrographic and in situ analyses of representative phenocrysts from the ore-forming gabbrodiorite porphyry of the Meishan iron deposit to gain new insights into the magmatic process. The phenocrysts are composed of clinopyroxene and plagioclase, which commonly exhibit disequilibrium features, such as reverse zoning and sieved textures, suggesting that the gabbrodiorite porphyry evolved in an open-system andesitic-magma reservoir into which mafic magma was repeatedly recharged. The zoned clinopyroxene phenocrysts show consistent trace-element distribution patterns and a narrow range of 87Sr/86Sr of 0.7037–0.7049, which suggest that the recharged mafic magma was the primitive basaltic magma and crystal fractionation dominated the evolution from basaltic magma to andesitic magma. Pressure estimations based on clinopyroxene phenocryst compositions indicate that the magma reservoirs were emplaced at the uppermost mantle. Comparison of clinopyroxene phenocryst compositions between the Meishan gabbrodiorite porphyry and published andesitic volcanics shows distinct rare earth element distribution patterns, which don't support a volcanic–subvolcanic relationship between the gabbrodiorite porphyry and andesitic volcanics. Estimated pressures indicate that the primitive magma of andesitic volcanics evolved in lower crust reservoir, which was much shallower than the reservoirs of the ore-forming gabbrodiorite porphyry. Therefore, the Meishan gabbrodiorite porphyry was not the subvolcanic intrusive of the andesitic volcanics and probably represented an independent magmatic pulse. This understanding suggests that the exploration target of iron deposits in the Ningwu basin is not necessarily limited to the contact zone between gabbrodiorite porphyry and andesitic volcanics.

Tracking timescales of magma reservoir recharge through caldera cycles at Santorini (Greece). Emphasis on an explosive eruption of Kameni Volcano

Pre-eruptive processes and their timescales are critical information for risk management at explosive volcanoes, and Santorini caldera (Greece) provides an excellent context in which to approach this subject. We ask two questions. First, are pre-eruptive processes the same for small and big eruptions? To investigate, we performed a multi-mineral diffusion timescale study of a small explosive eruption of Kameni Volcano and compared the results with those published for larger caldera-forming eruptions at Santorini. The Kameni dacite resembles products of larger eruptions in being crystal-poor, containing plagioclase with antecrystic cores and autocrystic rims, bearing orthopyroxene with sector zoning and phantom skeletal morphologies, and showing evidence for mixing of different silicic magmas prior to eruption. Diffusion timescales from Mg-Fe profiles in orthopyroxene and clinopyroxene phenocrysts are <1–23 years, and Mg diffusion modelling in plagioclase gives <10 years. Our physical model for the Kameni eruption is similar to those proposed for larger eruptions, where silicic melt produced in gabbroic to dioritic lower to middle crustal mush bodies is transferred (along with entrained mafic magma) to an upper crustal reservoir. Crystals grow in the hydrous silicic melts due to decompression, cooling, and magma mixing during ascent and injection into upper crust. We propose that large eruptions are preceded by similar processes as small ones, but on a larger scale. Our second question: do diffusion timescales relate to eruptive volume or position in a caldera cycle? For this, we obtained orthopyroxene Mg-Fe diffusion timescales for three additional eruptions, growing our orthopyroxene timescale database to seven eruptions of different sizes and cycle timings. No clear relationship exists between diffusion timescale and volume; however, timescales are systematically shorter (<0.01–10 years) early in a cycle and longer (1–5,000 years) late in a cycle. Thermal maturation and H2O-flushing of the crustal magma reservoir through the caldera cycle could explain this, as the reservoir would change from a rigid to more mushy state as the cycle progresses. This would change the mechanical response to melt input and allow accumulation of progressively larger melt layers in the upper crust, resulting in increasing crystal residence times.

Magma storage and transfer in the Villarrica volcanic chain, South Chile: MT insights into volcano-tectonic interactions

The active volcanic arc of the Andes is controlled by subduction of the Nazca and Antarctic Plates beneath the South American Plate. The Southern Volcanic Zone is the segment with the most active volcanoes in the Andes and the one where volcano-tectonic interactions are more evident. Magmatic fluids are there channeled through the brittle crust by faults and fractures, becoming an important region for understanding the mechanisms associated with melt migration and storage in subduction zones. The Liquiñe–Ofqui Fault System and Andean Transverse Faults, control the overall architecture of the volcanic arc and play an important role on magmatic transport and compositional partitioning. Evidence of this is the oblique chain composed of Villarrica, Quetrupillán and Lanín stratovolcanoes and a number of minor eruptive centers, where a wide range of magma compositions have been erupted with basalts and basaltic andesites dominating the suite. We used long–period and broad–band magnetotelluric stations deployed surrounding of Villarrica volcano to investigate the implications of crustal faulting and volcanism and their consequence on crustal reservoirs. Inversion of the data was used to generate three–dimensional electrical resistivity models. The resistivity distribution shows the upper crust as highly resistive, but below and east of Villarrica volcano, the model suggests the presence of a magmatic reservoir at shallow crustal levels (between 1.5 and 3 km b.s.l.), possibly corresponding to a transient magma storage. Meanwhile, the middle crust contains several intermediate to low conductive features interpreted as fluid pathways and/or melt storage regions, respectively, revealing the important role of fault systems. The lower crust also contains low resistivity zones indicating the presence of partial melt and/or fluids, associated with deep reservoirs (8–20 km), with a significant proportion of them likely non eruptible parts. This would suggest that melt is accumulated as highly crystallized mush or disconnected melt pockets and highlight the complex vertical extent of the structurally-controlled plumbing systems even in thin crust settings.

Development characteristics and main controlling factors of Carboniferous volcanic reservoirs in the Shixi area, Junggar Basin

The Carboniferous volcanic reservoirs in the Shixi area of the Junggar Basin are complex and diverse. Identifying the characteristics and main factors controlling high-quality volcanic reservoirs is the key to increasing oil and gas reserves and production in this area. Through core observations, thin section identification, physical property and pore structure analyses, combined with production data, the main controlling factors and development modes of high-quality reservoirs were analysed. The results show that the Carboniferous strata in the Shixi area mainly contain andesite and dacite of overflow facies, followed by volcanic breccia and tuff of explosive facies. Volcanic reservoirs in the study area are high-porosity–low-permeability and medium-porosity–low-permeability reservoirs. Volcanic breccia of explosive facies has the best physical properties, showing the characteristics of high porosity and medium permeability. The reservoir space is mainly composed of gas cavities, corrosion pores and fractures, among which the corrosion pores are the most important reservoir spaces of the Carboniferous volcanic rocks. Lithology and lithofacies, weathering and corrosion, and fractures are the main factors controlling the development of high-quality volcanic reservoirs. Volcanic rocks that had experienced weathering and denudation for a long time developed a large number of secondary corrosion pores due to the corrosion of soluble minerals or volcanic ash. Fractures further improved the physical properties, causing volcanic rocks to eventually develop into weathering crust reservoirs. The physical properties of the volcanic rocks far away from the weathering crust were improved through primary gas cavities and structural fractures, and these volcanic rocks eventually developed into the inner reservoir.

Quantitative characterization of vertical zonation of Mesozoic granite weathering reservoirs in the coastal area of eastern Fujian Province, China

Weathering crust reservoirs have obvious vertical zonation, which is the focus of weathering crust reservoir research, but there is a lack of quantitative characterization indexes. To achieve the quantitative characterization of granite weathering crust reservoir and provide the basis for oil exploration of granite weathering crust buried hill reservoir, in this paper, the vertical zonation of granite weathering crust reservoir is quantitatively divided by testing and analyzing the uniaxial compressive strength (UCS), magnetic susceptibility (MS), permeability, and chemical index of alteration (CIA) of the Mesozoic granite weathering crust in the coastal area of eastern Fujian. The results show that the granite weathering crust reservoir can be divided into four zones vertically: a soil zone (SZ), weathered and dissolved zone (WDZ), fracture zone (FZ), and bedrock zone (BZ). A cataclastic area is developed in the FZ and BZ, in which structural fractures are well-developed, the fracture surface density is usually greater than 200 m/m2, and the contribution to the fractures in the rock mass is up to about 50%, making this the sweet spot of the reservoir. In the SZ, the rocks are loose, and the pores are well-developed. The UCS is less than 10 MPa, and the average rate of change of the UCS (Δσ) is 0.90. The average permeability is 2823.00 mD, and the average rate of change of the permeability (Δκ) is 5.13. The average CIA is 74.9%. The average clay mineral content is 7%. The rocks in the WDZ have been significantly weathered by physical and chemical processes, and the weathering fractures and dissolution pores are well-developed. The average UCS is 18.2 MPa, and the average Δσ is 0.70. The average permeability is 143.80 mD, and average Δκ is 4.17. The average CIA is 65.3%. The average clay mineral content is 4%. Under the influence of tectonic movement and physical weathering, the rocks in the FZ have developed structural fractures and a few weathered fractures. The average UCS is 57.9 MPa, and the average Δσ is 0.18. The average permeability is 5.50 mD, and the average Δκ is 2.55. The average CIA is 61.6%. The average clay mineral content is 2%. In the BZ, the rocks are intact and hard. The average UCS is 69.9 MPa, and the average Δσ is 0.13. The average permeability is 1.46 mD, and the average Δκ is 1.43. The average CIA is 57.8%. The average clay mineral content is less than 1%. The multi-parameter combination of the UCS, Δσ, permeability, Δκ, CIA, and clay mineral content achieved good results in the division of the zones of the weathering crust. The UCS increases gradually from top to bottom, while Δσ, permeability, Δκ, CIA, and clay mineral content all decrease gradually. In addition, based on the petrophysical parameters of the rocks, including the density, resistivity, and acoustic velocity, a good division effect was also achieved, which can provide a basis for the vertical zonation of the granite buried-hill weathering crust reservoir.

A dataset of microscopic images of carbonate rocks of Leikoupo Formation of Middle Triassic in the Central Sichuan

The abundance of widely-developed carbonate rocks of the Leikoupo Formation in Central Sichuan Basin shows good exploration prospects and great exploration potential. Carbonate reservoirs in central Sichuan Basin are well formed, and various pores are well developed and preserved. Supergene dissolution, buried dissolution, and dolomitization have improved reservoir physical properties and is conducive to the reservoir permeability. Pore types including intragranular dissolved pores, intergranular dissolved pores, mold pores, etc. brings out connected dissolution fractures. At the top, a good paleoweathering crust reservoir is formed under karstification. By collecting the strata of Leikoupo Formation encountered in Well LG 19, Well LG 173, Well ZT 1, Well TN 112, Well Y 105, and Well PL 19 in the Central Sichuan Basin as the research object, we produced a dataset of rock micro thin section data. Leikoupo Formation in the Central Sichuan Basin is mainly composed of algal laminar dolomite, sandy dolomite, and breccia dolomite. The sedimentary facies are carbonate evaporation limited platform facies. The subfacies can be divided into tidal flat and intra platform beach subfacies. The tidal flat has microfacies such as cloud flat and gypsum cloud flat. And the intra platform beach can be subdivided into bioclastic beach and sandy beach. This dataset can be used to divide the sedimentary microfacies of the Leikoupo Formation in this area, study the control of the formation and distribution of micro relative reservoirs, determine the favorable reservoir microfacies, and provide effective guidance for oil and gas exploration and development. The image dataset can further provide basic data support for the production of machine learning image dataset.

Contraints on Petrogenesis and Fe Fertility of Intrusive Complexes in the Han–Xing Region, North China Craton from Apatite Geochemistry

The Han–Xing (Handan–Xingtai) region is famous for its endowment of skarn iron deposits in China. These deposits are mainly spatially and genetically associated with diorite rocks, but these rocks show different Fe ore potential. Major and trace elements composition of apatite from the Kuangshan and Fushan diorite complexes were investigated to explore the potential of apatite as a proxy of petrogenesis and Fe fertility of these rocks. All the investigated apatite grains are identified as fluorapatite, which is typical for magmatic apatite. The Sr, Y, Mn, and Heavy Rare Earth Elements (HREE) contents of apatite in the Kuangshan diorite complex are positively correlated with the increase of melt SiO2 content compared to that in the Fushan diorite complex. Apatite geochemistry indicates that magmas of the Fushan complex mainly experienced the fractional crystallization of hornblende in the deep crustal reservoirs, whereas the Kuangshan complex has experienced the fractional crystallization of hornblende in the deep and the shallow plagioclase fractional crystallization. The F, Cl and S content of the Kuangshan complex estimated by apatite volatile (F = 2632 ppm, Cl = 4100 ppm, SO3 = 140 ppm) is significantly higher than that of the Fushan complex (F = 2488 ppm, Cl = 3400 ppm, SO3 = 90 ppm). The Eu, Ce anomalies, Mn, and SO3 contents of apatite show that both of the two complexes have higher oxygen fugacity (Δ FMQ), but the oxygen fugacity of the Kuangshan complex calculated by Mn and SO3 content (Δ FMQ + 2.41) is higher than that of the Fushan complex (Δ FMQ + 1.77), which may also be one of the reasons for the great difference in ore-forming scale between the two complexes. Our results suggest that the high volatile contents and oxidation states of magma estimated by apatite, as well as the lower Sr/Y in apatite reflect favorable conditions for skarn iron mineralization. Therefore, our study shows that magmatic apatite geochemistry may be a useful tool to distinguish the Fe fertility of plutonic rocks related to skarn deposits.

Anatomy of a craton: Isotopic heterogeneity across an Archean crustal cross-section

The formation of granitoids, predominantly tonalite-trondhjemite-granodiorite (TTG) suites was a key process in the differentiation of Earth’s earliest crust. However, the relative importance of intracrustal reworking versus input from external sources remains unclear. U–Pb and Hf isotopic data collected from zircon in Archean granitoids frequently support reworking and partial melting of a common crustal reservoir over time to produce subsequent generations of melt. A less frequently observed isotopic signature, in a single region, implies the production of different generations of TTGs from unrelated sources. These end-member scenarios of granitoid generation have quite different implications for Archean crustal evolution. Here, we investigate TTGs from across a crustal cross-section known as the Kapuskasing Uplift in the Neoarchean Wawa terrane of the southern Superior Province. We report U–Pb and Hf isotopic data from zircon in TTGs from the lower, middle, and upper crust to evaluate their sources. Results indicate three groups of TTGs: a ca. 2885 Ma upper-crustal TTG with chondritic εHf, ca. 2835–2820 Ma lower-crustal TTGs with subchondritic εHf, and ca. 2750–2685 Ma TTGs with suprachondritic εHf emplaced at middle- to upper-crustal levels. Each group was derived from an isotopically distinct source. We propose that the Mesoarchean rocks exposed in the Kapuskasing Uplift represent a crustal fragment rifted from the margin of existing Mesoarchean terranes in the Superior Province. Subsequent Neoarchean juvenile magmatism dominates the upper and middle crust resulting in isotopic stratification of the Kapuskasing Uplift, as only the lower crust contains evidence of all three isotopic sources. Inferences about Archean crustal evolution in general rely on upper- to middle-crustal rocks (i.e., greenstone belts and grey gneiss), as these represent the majority of exposed Archean crust worldwide. The potential for the lower crust to be decoupled from overlying rocks implies that some parts of the Archean geological record may be preserved only at the deepest crustal levels—particularly rocks related to the most ancient sources.

The obscuring effect of magma recharge on the connection of volcanic-plutonic rocks

Abstract The current debate on volcanic-plutonic connection is centered on whether efficient liquid-crystal segregation dominates the evolution of a mushy reservoir to produce evolved, crystal-poor rhyolite and cumulate leftover. However, magma recharge may remarkably influence the evolution of a mushy reservoir and obscure the evidence of liquid-crystal segregation. This complexity poses a challenge to exploring the connection of volcanic-plutonic rocks. This study investigates the Qinzhou Bay granitic complex (~250–248 Ma) from South China, which contains crystal-poor (<19 vol%) peraluminous rhyolites and subsequent crystal-rich (28–54 vol%) porphyries. Although the rhyolite and porphyry units have a close spatio-temporal link, they do not share a fractionation trend and similar whole-rock Sr-Nd-O isotopic compositions; thus, a direct connection is not evidenced. We further present textural analyses, mineral and melt inclusion compositions, thermobarometry (the combination of Ti-in-zircon thermometer and Ti-in-quartz thermobarometer), and thermodynamic modeling to examine the alternative interpretations, i.e., the two units may have intrinsically independent origins or the connection of the two units has been obscured. For the rhyolite unit, thermobarometric results reveal a polybaric storage system consisting of middle (>600 ± 80 MPa) and upper (~150 ± 40 and ~60 ± 20 MPa) crustal reservoirs. Variations in quartz Fe content and chlorine-rich, metaluminous melt inclusions suggest that magma hybridization with less-evolved metaluminous magmas occurred at both crustal levels. In particular, the elevated Fe contents in the quartz population that crystallized at the shallowest level (~60 ± 20 MPa) suggest that recharge magmas were directly injected into the shallowest reservoir. Deviation of the whole-rock composition from the liquid evolution trend recorded in melt inclusions suggests a combined effect of magma mixing and crystal-melt segregation processes in upper crustal reservoirs. Thermodynamic modeling and mass balance calculations suggest that the whole-rock composition of the rhyolite could be reproduced by mixing between regionally exposed dacites and segregated melts at crystallinities of 50–60% (using parental magma represented by the least-evolved melt inclusion). For the porphyry unit, thermobarometric results reveal magma storage at middle (more than 450 ± 40 to 550 ± 40 MPa) and upper (110 ± 20 to 140 ± 20 MPa) crustal levels. The small-scale oscillatory zonation of plagioclase, the pervasive resorption of quartz and alkali feldspar, and the presence of peraluminous microgranular enclaves in the porphyries suggest a recharge event of metasediment-sourced magmas, triggering reactivation and convection of the reservoir. Autoclastic and overgrowth textures of quartz, plagioclase, and alkali feldspar phenocrysts and development of columnar jointing suggest that the reactivated porphyritic magmas ascended and emplaced at ultrashallow levels (~30 ± 10 MPa). Because of the similar storage pressures, the porphyries may represent remobilized cumulates of rhyolitic magmas, whereas the texture and geochemistry of the cumulate-liquid pair were modified, a key factor rendering a cryptic connection between the rhyolite and porphyry. Alternatively, the plumbing systems feeding the rhyolite and porphyry units are horizontally independent or vertically discrete, but this circumstance is inconsistent with the same evolution trend of quartz Fe and Al contents of the rhyolite and porphyry. Our study highlights that whole-rock composition may record blended information of complex processes, and caution should be taken when whole-rock composition is used to extract information of a single process. Multi-method constraints are required to evaluate the influence of recharge processes on the modification of liquid-cumulate records, and big data analysis on the basis of geochemistry should be conducted with caution to avoid biased understanding.

Assessing crustal contamination in the 28 to 18 Ma Dulce dike swarm with Nd-Sr isotopic data, southwestern Colorado

Alkaline to subalkaline mafic dikes in the 28 to 18 Ma Dulce swarm were emplaced in a zone of incipient extension from southern Colorado into northern New Mexico on the northeastern boundary of the San Juan Basin. The 87Sr/86Sr ratios for the dikes are 0.70503 to 0.70584, akin to most post-28 Ma mafic rocks across the northern San Juan Basin. These data are consistent with melting of metasomatized subcontinental lithospheric mantle with little to no crustal contribution as revealed by the geochemical and Sr-Nd isotopic signatures of most 28–0.6 Ma mafic rocks in the region. Time-corrected εNd(t) values of −4.1 to −7.4 for rocks in the Dulce swarm, however, indicate that magma production involved the crust. A previous hypothesis for Dulce magmas was contamination of lithospheric mantle melts with up to 45% mafic lower crust ± 0.5% upper crust. In this investigation, six new whole-rock Sr and Nd isotopic analyses were combined with published data to further investigate the contamination of lithospheric mantle melts with different crustal reservoirs. The Nd isotope signatures of the Dulce swarm offer evidence for the long-term involvement (∼10 Ma) of lower crust in the production of rift-related mantle magmas. Isotopic mixing curves support previous hypotheses for the contamination of lithospheric mantle melts with 10 to 40 percent lower mafic crust. This provides further insight into regional variations in mantle magmas produced after 28 Ma in the Four Corners region that likely triggered crustal melting related to caldera complexes in the western San Juan Mountains.

H2O-Controlled Eruptive Filtering on the Bimodality of Continental Volcanism Across Tectonic Settings

Abstract The contemporary basaltic and rhyolitic volcanism with a dearth of intermediate compositions creates a compositional gap (‘Daly gap’), which is known to be pronounced in continental rifts/hotspots but absent in continental arcs. The intrinsic mechanism(s) producing the Daly gap remain poorly understood. Similar bimodal spectra recorded in melt inclusions from continental rifts/hotspots and arcs suggest that similar processes modulate the bulk composition of erupted magmas. A compiled data set of melt inclusion composition, phenocryst fraction, pre-eruptive temperature and melt H2O content enables precise calculations of melt viscosities and bulk magma viscosities for >100 eruptions of global volcanoes in various tectonic settings. Analyses of the data set suggest lower viscosities of rhyolitic melts in arcs than in rifts/hotspots (~104.5 ± 0.2 versus ~105.5 ± 0.2 Pa·s, respectively), at comparable temperatures, because of ~1.5 wt % higher melt H2O contents in arcs. Combining numerical constraints on ascent dynamics and data set observations, our study defines an eruption threshold of ~106.0 ± 0.3 Pa·s at the base of a volcanic conduit rooted in the subvolcanic reservoir. The lower viscosities of rhyolitic melts in arcs facilitate stirring/re-homogenization of upper crustal reservoirs upon recharge and accommodation of more phenocrysts in erupted products, leading to the common appearance of crystal-rich (30–50%) andesites and dacites (with rhyolitic interstitial melts) with mixing signatures in magmatic provinces of continental arcs. In contrast, the lower mixing efficiency and eruptibility of the more viscous magmas present in continental rifts/hotspots lead to preferential eruption of crystal-poor (<~30%) rhyolites, enhancing compositional bimodality in the erupted records. This H2O-controlled eruptive filter modulates the presence versus absence of the Daly gap and is key to the distinct differentiation patterns across tectonic settings.

Insights into the Weathering Crust Reservoirs of Granitoids: A Case Study from Qinghai Oilfield of Qaidam Basin, Northwest China

With proven reserves of 9.836 × 1010 m3, the largest known natural gas reservoir among terrigenous basement rocks has been discovered within the granitoids of the northern Qaidam Basin. Due to their high heterogeneity, the genesis of basement reservoirs remains unknown. Herein, the structure of the weathering crust in granitoids and their potential controlling factors on the reservoir development mechanism are discussed using a multidisciplinary approach based on data from cores, thin sections, scanning electron microscopy (SEM), conventional and imaging logs, and physical property and major elements analyses. Moreover, the classification standard of the weathering crust structure is established. The dissolution belt holding diverse reservoir spaces accounts for more than 50% of the total porosity, while the disintegration belt is the main context for the development of cleavage fractures and crack fractures. The original pores exist mainly among the crystal grains of quartz and mica, while the secondary pores and fractures were generated by the alteration of aluminosilicate minerals as well as biotite or hornblende. The quality of these reservoirs is controlled by their mineral composition, tectonic uplift, faulting, and paleogeomorphology. The femic granitoid is the main reservoir-forming lithology in the case of dissolution, while the felsic granitoid is more likely to develop cracks. The formation of the disintegration belt is significantly linked to the presence of faulting. These belts were mostly induced by tectonic deformation along the Altyn fault belt from the late Oligocene to the early Miocene. The diversity in paleogeomorphology influences the extent of the weathering. The exhumation in the Altyn terrane from the late Jurassic to the Cenozoic corresponds to the weathering and hypergene leaching period of the weathering crust within granitoids. Three types of reservoirs are present in the rocks: fractured-porous (Type Ⅰ); porous (Type Ⅱ); and fractured (Type Ⅲ). The fractured-porous and fractured reservoirs were developed mainly in the granitic gneiss and granite, while the porous reservoir was formed in granitic diorite and granitic gneiss. The reservoirs that developed in the weathering crust of granitoids are dominated by Type Ⅰ and Type Ⅱ. The highest quality reservoir, which is the fractured-porous type, developed mainly in the dissolution belt of the weathering crust, and has a porosity ranging from 1.56% to 8.48% and a permeability ranging from 0.03 mD to 14.48 mD. The mechanisms of the development of weathering crust reservoirs provide further information for the hydrocarbon exploration of basement rocks worldwide.

Key factors controlling deep Carboniferous volcanic reservoirs in the east slope of Mahu Sag, Junggar Basin, NW China

Considering the complex factors controlling volcanic reservoirs, the Carboniferous strata in the eastern slope area of the Mahu Sag (ESMS) in the northwestern Junggar Basin (NJB) were investigated using rock cores, thin sections, scanning electron microscope (SEM), physical properties, major elements, X-ray fluorescence (XRF), well logging, and seismic data. The volcanic rocks revealed by drilling are mostly weathering crust reservoirs (WCRs), the formation of which in and around the study area is significantly controlled by weathering and leaching (WL). Most types of volcanic rock can be improved by long-term weathering. Favorable reservoirs in the ESMS are often developed within 150 m below the tectonic unconformity boundary at the top of the Carboniferous. The longer the weathering duration, the better are the overall quality of the WCRs. Weathering duration of about 40 Ma is probably an important threshold in the NJB. Ultra-long leaching of atmospheric water and strong late dissolution of acidic fluids before oil and gas accumulations are important for reservoir development and petroleum accumulation in volcanic strata filled with authigenic minerals, especially calcite. The early regional tectonic movement affected the volcanic eruption and controlled the lithofacies distribution. The linear density of fractures was negatively correlated with the distance from the main controlling fault. Owing to the relatively weak filling, high-angle fractures contribute significantly to the reservoir permeability. The physical properties of volcanic breccia are better than those of tuff, and the porosity, permeability, and fracture density of andesite are higher than those of basalt. The physical properties of near-source facies belts of a volcanic edifice are better than those of far-source facies ones. Favorable exploration areas are the structural highs and fault zones where the duration of WL is more than 40 Ma, explosive facies and effusive facies near the crater are developed, or the inherited ancient buried hills transformed by faults and fractures near excellent source rocks, where the dissolution of atmospheric water and organic acidic fluids are strong.

A rapid change in magma plumbing taps porphyry copper deposit-forming magmas

Porphyry-type deposits are a vital source of green technology metals such as copper and molybdenum. They typically form in subduction-related settings from large, long-lived magmatic systems. The most widely accepted model for their formation requires that mantle-derived magmas undergo an increase in volatiles and ore-forming constituents in mid- to lower crustal reservoirs over millions of years, however, this is mostly based on observations from shallow, sporadically exposed parts of porphyry systems. To examine this paradigm, we have evaluated the timeframe and geochemical signatures of magmatism in a ~ 8 km palaeodepth cross-section through plutonic and volcanic rocks of the classic Yerington magmatic system, Nevada. We show that the magmas in the upper parts of the system (< 8 km) underwent a major and rapid change in chemistry over a period of < 200 kyrs that is coincident with the initiation of ore formation. We attribute this change to a shift from extraction of quartz monzodiorite and quartz monzonite magmas evolving in mid-crustal reservoirs, and that had relatively poor ore-forming potential, to extraction of volatile-rich granitic magmas from greater (~ 30 km) depths. As the granites crystallised, late stage melts were intruded through the carapace as aplite dykes which contain traceable expressions of the porphyry deposit-forming fluids. The rapid nature of the shift in ore-forming potential narrows the temporal-geochemical footprint of magmas associated with porphyry mineralisation and provides new constraints for exploration models.

Petrogenesis of submarine volcanic rocks dredged from Dokdo, Ulleungdo, and the neighboring seamounts in the East Sea: Constraints from mineral chemistry, geochemistry, and 40Ar/39Ar ages

Several seamounts, including the Dokdo (DD), Ulleungdo (UD) Anyongbok (AYB), Simheungtack (SHT), and Isabu, were formed from intraplate volcanism in the southwest of the East Sea back-arc basin during the Quaternary. In the present study, mineralogy and geochemistry data generated using submarine volcanic rocks collected from the DD, UD, AYB, and SHT locations were utilized to investigate characteristics of the source, evolution, and storage conditions of the associated magmas. Analyses of representative samples from the DD, UD, and AYB using the 40 Ar/ 39 Ar method yielded plateau ages of 1.566–2.658, 0.533–0.714, and 1.736 Ma, respectively. The dredged samples consist of trachybasalt, phonotephrite, trachyandesite, and trachyte and have oceanic island basalt (OIB)-like characteristics. According to numerical modeling using trace element data, the primary magma linked to the formation of these samples can originate from ≤5% partial melting of an enriched garnet–peridotite facies metasomatized by melt derived from a subducted slab. The composition of major oxides, from the trachybasalts to trachytes, and mass balance calculations suggest that the trachyandesites and trachytes were formed via fractional crystallization from a basaltic melt. These diverse rock compositions probably highlight pulses of magma that separated from the main reservoir in the shallow crust. Results from mineral–melt thermometry indicate that the olivine and plagioclase in these basaltic rocks crystallized at temperatures of 1213 ± 56 °C and 1163 ± 23 °C, respectively. Thermobarometric calculations using clinopyroxene–melt pairs indicate that the Quaternary volcanism was associated with two magma storage regions that correspondingly involved pressures of 791 ± 253 MPa (1098 ± 34 °C) and 187 ± 123 MPa (938 ± 77 °C). The primary magma was likely present in a reservoir beneath the Moho (17–32 km), and it evolved by fractional crystallization to produce rocks, such as trachybasalts and phonotephrites. Crystallization and differentiation of the basaltic magma that was transported to crustal reservoirs explain the evolved compositions. The trachytic magma generated were present at depths of 2–11 km, and crystallization produced a paragenetic assemblage comprising clinopyroxene, amphibole, plagioclase, biotite, apatite, and titanomagnetite. • Quaternary intraplate volcanisms occurred in the East Sea back-arc basin. • New 40 Ar/ 39 Ar ages of submarine volcanic rocks were obtained. • Basaltic rocks were originated from an enriched garnet-bearing mantle source. • Trachytes were formed by fractional crystallization from basaltic magma. • The Quaternary volcanism had at least two magma storage regions.