Geology ReportsSearch

USGS · 70202464

Isotopic and petrologic investigation, and a thermomechanical model of genesis of large-volume rhyolites in arc environments: Karymshina Volcanic Complex, Kamchatka, Russia

Abstract

The Kamchatka Peninsula of eastern Russia is currently one of the most volcanically active areas on Earth where a combination of >8 cm/yr subduction convergence rate and thick continental crust generates large silicic magma chambers, reflected by abundant large calderas and caldera complexes. This study examines the largest center of silicic 4-0.5 Ma Karymshina Volcanic Complex, which includes the 25 × 15 km Karymshina caldera, the largest in Kamchatka. A series of rhyolitic tuff eruptions at 4 Ma were followed by the main eruption at 1.78 Ma and produced an estimated 800 km 3 of rhyolitic ignimbrites followed by high-silica rhyolitic post-caldera extrusions. The postcaldera domes trace the 1.78 Ma right fracture and form a continuous compositional series with ignimbrites. We here present results of a geologic, petrologic, and isotopic study of the Karymshina eruptive complex, and present new Ar-Ar ages, and isotopic values of rocks for the oldest pre- 1.78 Ma caldera ignimbrites and intrusions, which include a diversity of compositions from basalts to rhyolites. Temporal trends in δ 18 O, 87 Sr/ 86 Sr, and 144 Nd/ 143 Nd indicate values comparable to neighboring volcanoes, increase in homogeneity, and temporal increase in mantle-derived Sr and Nd with increasing differentiation over the last 4 million years. Data are consistent with a batholithic scale magma chamber formed by primarily fractional crystallization of mantle derived composition and assimilation of Cretaceous and younger crust, driven by basaltic volcanism and mantle delaminations. All rocks have 35–45% quartz, plagioclase, biotite, and amphibole phenocrysts. Rhyolite-MELTS crystallization models favor shallow (2 kbar) differentiation conditions and varying quantities of assimilated amphibolite partial melt and hydrothermally-altered silicic rock. Thermomechanical modeling with a typical 0.001 km 3 /yr eruption rate of hydrous basalt into a 38 km Kamchatkan arc crust produces two magma bodies, one near the Moho and the other engulfing the entire section of upper crust. Rising basalts are trapped in the lower portion of an upper crustal magma body, which exists in a partially molten to solid state. Differentiation products of basalt periodically mix with the resident magma diluting its crustal isotopic signatures. At the end of the magmatism crust is thickened by 8 km. Thermomechanical modeling show that the most likely way to generate large spikes of rhyolitic magmatism is through delamination of cumulates and mantle lithosphere after many millions of years of crustal thickening. The paper also presents a chemical dataset for Pacific ashes from ODDP 882 and 883 and compares them to Karymshina ignimbrites and two other Pleistocene calderas studied by us in earlier works.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Ilya N. Bindeman, Vladimir L. Leonov, Dylan P. Colon, Aleksey N. Rogozin, Niccole Shipley, Brian Jicha, Matthew W. Loewen, Taras V. Gerya. 2019-01-24. Isotopic and petrologic investigation, and a thermomechanical model of genesis of large-volume rhyolites in arc environments: Karymshina Volcanic Complex, Kamchatka, Russia. https://doi.org/10.3389/feart.2018.00238

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Pliocene (Piacenzian) planktic foraminiferal biogeography: Insights into how planktic biogeography is linked to climate and socioeconomic impacts on marine ecosystems

The Piacenzian Age (3.60–2.58 Ma) of the Pliocene Epoch was characterized by globally warmer climates, higher sea levels, and atmospheric CO 2 similar to present. Utilizing a robust dataset of 2,101 samples and over 637,000 foraminifer specimens from 77 deep-sea core sites worldwide, we document planktic foraminifer biogeography and biodiversity during the mid-Piacenzian Warm Period (mPWP). Cluster analysis and multidimensional scaling reveal five major bioregions: tropical, warm subtropical, transitional, polar, and a distinct North Atlantic polar bioregion. Each bioregion was dominated by characteristic species with well-established temperature preferences. Analyses demonstrate higher species richness and evenness in low and mid-latitudes, with increased diversity associated with periods of climatic warming and poleward expansion of warm water assemblages. The long-term stability of biogeographic patterns underscores ecological conservatism but also highlights potential vulnerability to rapid anthropogenic climate change. Our findings emphasize the critical role of planktic foraminifers in reconstructing past ocean conditions and offer valuable insights into links between planktic biogeography, climate, and socioeconomic impacts on marine ecosystems. This work advances our understanding of marine ecosystem responses to climate extremes and provides a foundation for future regional and temporal analyses of planktic foraminifer biogeography under global change scenarios.

Frontiers in Earth Science

Estimation of dynamic geologic CO2 storage resources in the Illinois Basin, including effects of brine extraction, anisotropy, and hydrogeologic heterogeneity

Since the vast majority of carbon dioxide (CO 2 ) storage resources in the United States are in deep saline aquifers, optimizing the use of these saline storage resources could be crucial for efficient development of geologic CO 2 storage (GCS) resources and basin- or larger-scale deployment of GCS in the country. Maximum CO 2 injection rates can be enhanced by extracting brine from the CO 2 storage unit. However, disposal of the extracted brine is both a technological and economic challenge. The lowest-cost option would likely be reinjection of the extracted brine into another formation above or below the CO 2 storage unit. Therefore, it is important to estimate brine injectivity as it will constrain the potential to increase CO 2 injectivity at an injection site that has access to multiple geologic storage units where either CO 2 or brine can be injected. Using a simulation-optimization framework, coupled with a non-isothermal, multiphase CO 2 -water-salt equation-of-state module, we developed a computationally efficient method for evaluating optimization of simultaneous CO 2 injection, brine extraction, and brine (re)injection at hypothetical injection sites deployed across a geologic basin. The Illinois basin is ideal for testing our methodology because it contains multiple geologic storage units with seals in between them to isolate injection of CO 2 in one unit from interfering with the injection of either brine or CO 2 in another unit above or below it. In addition, we investigated the relative effects of variation in key geologic parameters as well as two reservoir structures (hydrogeologic heterogeneity/anisotropy and homogeneity/isotropy) on CO 2 injectivities and enhancement of CO 2 injectivity through extracting brine. Results suggest that permeability, depth, and especially thickness of the storage unit could be the most influential parameters determining CO 2 injectivity. They also suggest that only injecting CO 2 into the storage unit with the greatest injectivity, enhancing that unit’s injectivity by extracting brine, and disposing of the produced brine in other suitable units could maximize total CO 2 injectivity in limited regions of the basin. At the majority of simulated injection sites, however, we found that injecting CO 2 into all of the accessible and suitable storage units was more likely to maximize the CO 2 storage resource.

Illinois, Kentucky, Indiana

An exploration of the relative influence of physical models for Omori’s law

Omori’s law states that the rate of aftershocks decays as a function of inverse time. There are multiple physical explanations that we reduce into a nonlinear mixed effects relation of three terms: (1) a Rate/State expression that can account for static/dynamic and viscoelastic triggering caused directly by the mainshock, (2) a fluid diffusion triggering term, and (3) a randomized secondary triggering (cascade) term. We fit free physical-model parameters to an observed aftershock sequence through two nonlinear regression methods to find the relative contributions of physics-based models in an observed aftershock sequence. Results from both methods show that Rate/State models overpredict aftershock rates by ∼0–30%. Secondary aftershocks cause a net negative contribution (seismicity rate reduction that corrects overprediction by other terms) ranging between ∼0 and 30%. All regression solutions yield negative secondary triggering contributions without being guided to do so. A physical explanation for this is that aftershock occurrence relieves stress from the crust, ultimately causing the sequence to extinguish itself. Fluid diffusion triggering contributions range from ∼0 to 20%. Diffusion processes are observed to be shorter in time than the full duration of an aftershock sequence and they are also spatially limited, diminishing their influence. Our results apply to an aftershock decay curve from the 2016 Central Apennines earthquake sequence, meaning that our specific results may not be general. Our primary conclusion is that any one physical model cannot alone fit the observed sequence as well as the combination of three we investigated.

Frontiers in Earth Science