Search

UNLABELLED: On 15 January 2022, Hunga volcano erupted, creating an extensive and high-reaching umbrella cloud over the open ocean, hindering traditional isopach mapping and fallout volume estimation. In MODIS satellite imagery, ocean surface water was discolored around Hunga following the eruption, which we attribute to ash fallout from the umbrella cloud. By relating intensity of ocean discoloration to fall deposit thicknesses in the Kingdom of Tonga, we develop a methodology for estimating airfall volume over the open ocean. Ash thickness measurements from 41 locations are used to fit a linear relationship between ash thickness and ocean reflectance. This produces a minimum airfall volume estimate of 1.8-0.4+0.3 km3. The whole eruption produced > 6.3 km3 of uncompacted pyroclastic material on the seafloor and a caldera volume change of 6 km3 DRE. Our fall estimates are consistent with the interpretation that most of the seafloor deposits were emplaced by gravity currents rather than fall deposits. Our proposed method does not account for the largest grain sizes, so is thus a minimum estimate. However, this new ocean-discoloration method provides an airfall volume estimate consistent with other independent measures of the plume and is thus effective for rapidly estimating fallout volumes in future volcanic eruptions over oceans. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s00445-024-01744-6.

Deformation of all materials necessitates the collective propagation of various microscopic defects. On Earth, fracturing gives way to crystal-plastic deformation with increasing depth resulting in a brittle-to-ductile transition (BDT) region that is key for estimating the integrated strength of tectonic plates, constraining the earthquake cycle, and utilizing deep geothermal resources. Here, we show that the crossing of a BDT in marble during deformation experiments in the laboratory is accompanied by systematic increase in the frequency of acoustic emissions suggesting a profound change in the mean size and propagation velocity of the active defects. We further identify dominant classes of emitted waveforms using unsupervised learning methods and show that their relative activity systematically changes as the rocks cross the brittle-ductile transition. As pressure increases, long-period signals are suppressed and short-period signals become dominant. At higher pressures, signals frequently come in avalanche-like patterns. We propose that these classes of waveforms correlate with individual dominant defect types. Complex mixed-mode events indicate that interactions between the defects are common over the whole pressure range, in agreement with postmortem microstructural observations. Our measurements provide unique, real-time data of microscale dynamics over a broad range of pressures (10 to 200 MPa) and can inform micromechanical models for semi-brittle deformation.

Kinetic calcium isotope effects can be used as growth-rate proxies for volcanic and subvolcanic minerals. Here, we analyze Ca isotopic compositions in experimental and natural samples and confirm that large kinetic effects (>2‰) can occur during magmatic plagioclase crystallization. Experiments confirm theoretical predictions that disequilibrium isotope effects depend mainly on the rates for crystal growth relative to liquid phase Ca diffusivity (R/D). Plagioclase phenocrysts from the 1915 Mount Lassen rhyodacite eruption, the ∼650-y-old Deadman Creek Dome eruption, and several mafic subvolcanic orbicules and plagioclase comb layers from Northern California have disequilibrium Ca isotopic compositions that suggest rapid crystal growth rates (>1 cm/y to 15 cm/y). The Ca isotope results, combined with complementary crystal-size distribution analyses, suggest that magmatic rejuvenation (and eruption) events, as reflected in crystal growth times, can be as short as ∼10-3 y. Although mafic systems are predicted to have shorter magmatic rejuvenation periods, we find similarly short timescales in both mafic and silicic systems. These results are consistent with a growing body of evidence suggesting that dominantly crystalline volcanic magma reservoirs can be rapidly reactivated by the injection of fresh magma prior to eruption. By focusing on a common mineral such as plagioclase, this approach can be applied across all major magmatic compositions, suggesting that Ca isotopes can be used as a tool for investigating the dynamics and timing of volcanic eruptions.

Late Cretaceous records of environmental change suggest that Deccan Traps (DT) volcanism contributed to the Cretaceous-Paleogene boundary (KPB) ecosystem crisis. However, testing this hypothesis requires identification of the KPB in the DT. We constrain the location of the KPB with high-precision argon-40/argon-39 data to be coincident with changes in the magmatic plumbing system. We also found that the DT did not erupt in three discrete large pulses and that >90% of DT volume erupted in

Continental flood basalt (CFB) provinces are products of the largest known volumetric eruptions on Earth (~104 km3), with individual flow fields commonly covering well over 10 000 km2 with a mean lava thickness of over 5 m. Studies focusing on the emplacement style of such lava flows have relied extensively on morphological observations and comparisons with modern lava flows and experimental analogs. In the present study, we compare the textures of flood basalt lavas with those from different eruption settings all over the world using data collected from pre-existing literature to gain detailed insights into the style of eruption. Comparison of crystal size distribution data indicates that the eruption style of CFBs is similar to those of modern-day fissure eruptions (e.g. Iceland). This matches inferences based on observations of morphology. We also use a 1D thermal model to estimate the depth-dependent cooling rates within a single lava lobe and test the validity of assumptions built into the formulation of these models for large scale flood basalt lavas. The results reveal that, on average, flood basalt lavas need to conductively cool much faster than we would expect (up to order of ~102 times faster) to match the textural observations. The model is also frequently unable to replicate the observed depth-wise relative variations in length with depth for CFB lavas. Furthermore, the calculated cooling rates from crystal shapes also do not match those calculated from crystal lengths, indicating the assumptions in cooling flow models need to be modified for large CFB flow fields. Given the large areas of CFB flow fields and the relatively long eruption times inferred for the emplacement of individual flow fields, we hypothesize that inflation of lobes and formation of new lobes via breakouts combined with variable eruption rates are key processes that are missing when modeling the cooling of these flow fields. Accounting for these processes is essential to derive accurate cooling rates, which is important to better understand the environmental impact CFBs have at the time of emplacement. [ABSTRACT FROM AUTHOR]

High‐temperature–low‐pressure metamorphism is commonly associated with intermediate to felsic magmatism in continental orogenic belts. The heat budgets and transfer mechanisms responsible for such elevated temperatures and partial melting of the upper crust are uncertain. The Trois Seigneurs massif, French Pyrenees, preserves a structurally continuous record of Variscan high‐temperature–low‐pressure metamorphism through a sequence of upper‐to‐mid‐crustal Paleozoic metasedimentary rocks. Conventional thermobarometry and phase equilibria calculations show that metamorphic conditions span ~2.5 kbar, 575°C to suprasolidus conditions of ~6 kbar, 700°C. Peak temperatures depend strongly on depth: temperature gradients of 50–60°C/km are present through the uppermost 12 km of the section; deeper portions (12–20 km) define restricted temperature conditions of ~650–700°C. The lowest‐grade metamorphic rocks preserve the largest spread in monazite 206Pb*/238U dates, from c. 325–285 Ma, while the spread in dates is restricted to c. 305–290 Ma in the highest‐grade rocks. Within this spread, each sample yields a well‐defined population of monazite 206Pb*/238U dates with peaks at c. 305 Ma in the andalusite schists, 295 Ma in the sillimanite schists, and 300 Ma in the migmatite sample. Monazite trace‐element compositions capture a systematic change with decreasing date and increasing metamorphic grade, including a more negative Eu‐anomaly and decreasing Sr concentrations, consistent with co‐crystallizing feldspar; increasing HREE and Y contents, consistent with xenotime breakdown; and decreasing Th/U, reflecting increasing U content during breakdown of inherited zircon. Zircon rims from a granite unit that formed via partial melting of the Paleozoic sedimentary package yields a 206Pb/238U‐207Pb/235U concordia age of 304.1 ± 3.73 Ma. These rims have trace‐element compositions reflecting cogenetic apatite and zircon growth during granite formation. Zircon from a calc‐alkaline granodiorite intrusion preserves a 40 Ma record of melt‐related activity in the lower crust that preceded the regional thermal climax. We interpret these petrochronological data to show that the Trois Seigneurs field gradient including andalusite schist and biotite granite samples represents a genuine geotherm through Variscan orogenic crust during the regional thermal climax at 305 Ma. When combined with constraints from other Pyrenean massifs, the form of the geotherm is consistent with a thermal scenario in which heat is advected to the upper crust by intermediate‐composition magmas generated in the lower crust. A simple thermal model for this process indicates that anatexis in the upper crust may plausibly occur within 10 Ma of the initiation of the lower‐crustal melting. Such a thermal scenario, however, requires focusing of melt through a fertile lower crust and an elevated Moho heat flux. We suggest that this process may have controlled the attainment of high‐temperature–low‐pressure metamorphic conditions along the Variscan belt and may currently be operating in zones of post‐orogenic continental extension. [ABSTRACT FROM AUTHOR]

Frictional properties affect the propagation of high-amplitude seismic waves across rock fractures and faults. Laboratory evidence suggests that these properties can be measured in active seismic surveys, potentially offering a route to characterizing friction in situ. We present experimental results from a subresonance torsional modulus and attenuation apparatus that utilizes micron-scale sinusoidal oscillations to probe the nonlinear stress-strain relation at a range of strain amplitudes and rates. Nonlinear effects are further quantified using harmonic distortion; however, time series data best illuminate underlying physical processes. The low-frequency stress-strain hysteretic loops show stiffening at the sinusoid's static ends, but stiffening is reduced above a threshold frequency. This shape is determined by harmonic generation in the strain; the stress signal has no harmonics, confirming that the fractured sample is the source of the nonlinearity. These qualitative observations suggest the presence of rate-dependent friction and are consistent between fractures in three different rock types. We propose that static friction at the low strain rate part of the cycle, when given sufficient “healing” time at low oscillation frequencies, causes this stiffening cusp shape in the hysteresis loop. While rate-and-state friction is commonly used to represent dynamic friction, it cannot capture static friction or negative slip velocities. So we implement another dynamic friction model, based on the work of Dahl, which describes this process and produces similar results. Since the two models have a similar form, parameterizations of field data could constraint fault model inputs, such as specific location velocity strengthening or weakening properties.

We present the first results from the polarimetry mode of the Gemini Planet Imager (GPI), which uses anew integral field polarimetry architecture to provide high contrast linear polarimetry with minimal systematic biases between the orthogonal polarizations. We describe the design, data reduction methods, and performance of polarimetry with GPI. Point-spread function (PSF) subtraction via differential polarimetry suppresses unpolarized starlight by a factor of over 100, and provides sensitivity to circumstellar dust reaching the photon noise limit for these observations. in the case of the circumstellar disk around HR 4796A,GPI's advanced adaptive optics system reveals the disk clearly even prior to PSF subtraction. In polarized light, the disk is seen all the way in to its semi-minor axis for the first time. The disk exhibits surprisingly strong asymmetry in polarized intensity, with the west side ≳9 times brighter than the east side despite the fact that the east side is slightly brighter in total intensity. Based on a synthesis of the total and polarized intensities, we now believe that the west side is closer to us, contrary to most prior interpretations. Forward scattering by relatively large silicate dust particles leads to the strong polarized intensity on thewest side, and the ring must be slightly optically thick in order to explain the lower brightness in total intensity there. These findings suggest that the ring is geometrically narrow and dynamically cold, perhaps shepherded by larger bodies in the same manner as Saturn's Fring.

Deformation of all materials necessitates the collective propagation of various microscopic defects. On Earth, fracturing gives way to crystal-plastic deformation with increasing depth resulting in a "brittle-to-ductile" transition (BDT) region that is key for estimating the integrated strength of tectonic plates, constraining the earthquake cycle, and utilizing deep geothermal resources. Here, we show that the crossing of a BDT in marble during deformation experiments in the laboratory is accompanied by systematic increase in the frequency of acoustic emissions suggesting a profound change in the mean size and propagation velocity of the active defects. We further identify dominant classes of emitted waveforms using unsupervised learning methods and show that their relative activity systematically changes as the rocks cross the brittle-ductile transition. As pressure increases, long-period signals are suppressed and short-period signals become dominant. At higher pressures, signals frequently come in avalanche-like patterns. We propose that these classes of waveforms correlate with individual dominant defect types. Complex mixed-mode events indicate that interactions between the defects are common over the whole pressure range, in agreement with postmortem microstructural observations. Our measurements provide unique, real-time data of microscale dynamics over a broad range of pressures (10 to 200 MPa) and can inform micromechanical models for semi-brittle deformation. [ABSTRACT FROM AUTHOR]