Age of the magma chamber and its physicochemical state under Elbrus Greater Caucasus, Russia using zircon petrochronology and modeling insights.

Mount Elbrus, Europe's tallest and largely glaciated volcano, is made of silicic lavas and is known for Holocene eruptions, but the size and state of its magma chamber remain poorly constrained. We report high spatial resolution U-Th-Pb zircon ages, co-registered with oxygen and hafnium isotopic values, span ~ 0.6 Ma in each lava, documenting magmatic initiation that forms the current edifice. The best-fit thermochemical modeling constrains magmatic fluxes at 1.2 km3/1000 year by hot (900 °C), initially zircon-undersaturated dacite into a vertically extensive magma body since ~ 0.6 Ma, whereas a volcanic episode with eruptible magma only extends over the past 0.2 Ma, matching the age of oldest lavas. Simulations explain the total magma volume of ~ 180 km3, temporally oscillating δ18O and εHf values, and a wide range of zircon age distributions in each sample. These data provide insights into the current state (~ 200 km3 of melt in a vertically extensive system) and the potential for future activity of Elbrus calling for much-needed seismic imaging. Similar zircon records worldwide require continuous intrusive activity by magmatic accretion of silicic magmas generated at depths, and that zircon ages do not reflect eruption ages but predate them by ~ 103 to 105 years reflecting protracted dissolution-crystallization histories.

Pervasive Hydrothermal Events Associated with Large Igneous Provinces Documented by the Columbia River Basaltic Province.

The degree and extent of crustal hydrothermal alteration related to the eruption of large igneous provinces is poorly known and not easily recognizable in the field. We here report a new δ18O dataset for dikes and lavas from the Columbia River Basalt Group (16-15 Ma) in the western USA, and document that dikes on average are 1-2‰ more depleted in δ18O than basalt flows. We show that this observation is best explained with the involvement of heated meteoric  waters during their cooling in the crust. The largest 6-8‰ depletion is found around and inside a 10 m-thick feeder dike that intruded the 125 Ma Wallowa tonalitic batholith. This dike likely operated as a magma conduit for 4-7 years, based on the extent of heating and melting its host rocks. We show that this dike also created a hydrothermal system around its contacts extending up to 100 m into the surrounding bedrock. A model that considers (a) hydrothermal circulation around the dike, (b) magma flow and (c) oxygen isotope exchange rates, suggests that the hydrothermal system operated for ~150 years after the cessation of magma flow. In agreement with a previously published (U-Th)/He thermochronology profile, our model shows that rocks 100 m away from such a dike can be hydrothermally altered. Collectively, our sample set is the first documentation of the widespread hydrothermal alteration of the shallow crust caused by the intrusion of dikes and sills of the Columbia River Basalt Province. It is estimated that heating and hydrothermal alteration of sediments rich in organic matter and carbonates around the dikes and sills releases 18 Gt of greenhouse gases (CH4 and CO2). Furthermore, hydrothermal δ18O depletion of rocks around dikes covers 500-600 km3, which, when scaled to the total CRB province constitutes 31,000 km3 of low-δ18O rocks. These volumes of crust depleted in δ18O are sufficient to explain the abundant low-δ18O magmas in eastern Oregon and western Idaho. This work also demonstrates that the width and magnitude of δ18O depletion around dikes can identify them as feeders. Given this, we here interpret Paleoproterozoic dikes in Karelia with the world's lowest δ18O depletions (-27.8‰) as feeders to the coeval large igneous province aged 2.2-2.4 Ga that operated under the Snowball Earth glaciation conditions.

Hydrothermal alteration of kimberlite by convective flows of external water.

Kimberlite volcanism involves the emplacement of olivine-rich volcaniclastic deposits into volcanic vents or pipes. Kimberlite deposits are typically pervasively serpentinised as a result of the reaction of olivine and water within a temperature range of 130-400 °C or less. We present a model for the influx of ground water into hot kimberlite deposits coupled with progressive cooling and serpentisation. Large-pressure gradients cause influx and heating of water within the pipe with horizontal convergent flow in the host rock and along pipe margins, and upward flow within the pipe centre. Complete serpentisation is predicted for wide ranges of permeability of the host rocks and kimberlite deposits. For typical pipe dimensions, cooling times are centuries to a few millennia. Excess volume of serpentine results in filling of pore spaces, eventually inhibiting fluid flow. Fresh olivine is preserved in lithofacies with initial low porosity, and at the base of the pipe where deeper-level host rocks have low permeability, and the pipe is narrower leading to faster cooling. These predictions are consistent with fresh olivine and serpentine distribution in the Diavik A418 kimberlite pipe, (NWT, Canada) and with features of kimberlites of the Yakutian province in Russia affected by influx of ground water brines. Fast reactions and increases in the volume of solid products compared to the reactants result in self-sealing and low water-rock ratios (estimated at <0.2). Such low water-rock ratios result in only small changes in stable isotope compositions; for example, δO18 is predicted only to change slightly from mantle values. The model supports alteration of kimberlites predominantly by interactions with external non-magmatic fluids.