Melt Migration Research Articles

Zircon Coupled Dissolution–Precipitation Replacement During Melt–Rock Interaction Modifies Chemical Signatures Resulting in Misleading Ages

ABSTRACTMelt migration through Earth's crust drives well‐documented melt–rock reactions, locally changing rock assemblage and geochemistry. However, melt–zircon interaction remains understudied. We report on three zircon‐melt interaction events from the Pembroke Granulite, New Zealand. Primary zircon from gabbroic gneiss which was subject to minor post‐emplacement melt migration and primary zircon from younger dykes exhibit straightforward microstructures, microchemistry, and age data. In contrast, zircon from melt‐mediated reaction halos adjacent to the dykes and from melt‐fluxed high‐strain zones display dissolution modification of grains, micro‐porosity and blurred or truncated internal zoning typical of replacement by coupled dissolution‐precipitation. Replaced zircon domains show changed rare earth element patterns and redistributed or lost radiogenic Pb that generates ambiguous apparent spot date arrays, smeared over tens of millions of years. We conclude that the metamorphism and three melt–rock interaction events were brief, and the arrays misrepresent the true age and duration of the metamorphism. Pb‐loss persisted beyond the metamorphism, with porosity and inclusions formed during coupled dissolution‐precipitation making replaced zircon domains more susceptible to subsequent Pb‐loss compared to the structurally intact, primary magmatic zircon in the host gabbroic gneiss or dykes. We recommend conducting high‐resolution microstructural investigations upon recognition of spot date arrays observed in single samples to rule out the possibility of spurious arrays resulting from coupled dissolution‐precipitation.

Composite Melt–Rock Interactions in the Lowermost Continental Crust: Insights from a Dunite–Pyroxenite–Gabbronorite Association of the Mafic Complex from the Ivrea–Verbano Zone (Italian Alps)

Abstract The processes leading to the building of the continental crust through magmatic underplating are fundamentally unknown, mainly because of the rare accessibility to deep level sections of the continental crust. The Italian Alps expose the Permian Mafic Complex, an 8-km-thick gabbronorite–diorite batholith that intruded the lower continental crust during the post-Variscan transtensional tectonics. We present here a petrological and geochemical study of a concentric dunite–pyroxenite–gabbronorite association, called Monte Mazzucco sequence, enclosed at deep levels of the Mafic Complex, thereby allowing us to provide new insights into the magmatic processes driven by emplacement of mantle melts in deep crustal continental areas. The studied sequence includes a ~ 60-m-thick dunite lens, in which olivine (82 mol % forsterite) is associated with accessory Cr-spinel including blebs and lamellae made up of magnetite. The dunite lens is permeated by mm- to cm-scale thick magmatic veins, which range in composition from hornblende lherzolite to olivine hornblendite and hornblende websterite. The lens is mantled by a m-scale ring consisting of amphibole-bearing (≤1 vol %) websterite, and the websterite ring is in turn enclosed by amphibole-free gabbronorites. Both magmatic veins within the dunites and mantling websterites typically include an oxide association of Al-spinel and magnetite. Remarkably, the hornblende websterite veins and the mantling websterites are typically plagioclase-free and include clinopyroxene and amphibole with chondrite-normalized rare earth element patterns characterized by negative Eu anomaly. The mantling websterites display a subtle, gradual outward decrease of Mg# for orthopyroxene, clinopyroxene and accessory olivine, coupled with an increase of the negative Eu anomaly in clinopyroxene and amphibole. The enclosing gabbronorites are amphibole-free and have a chemically evolved signature. We propose a petrogenetic scenario including two major events of melt–dunite interaction. The first resulted from focused reactive melt infiltration and formed the magmatic veins within dunites. The hornblende lherzolite and the olivine hornblendite veins were produced by focused reactive melt migration through dunite grain boundaries, involving dissolution of olivine and recrystallization of Cr-spinel into Al-spinel and magnetite, whereas the hornblende websterite veins crystallized from melts penetrating through narrow fractures and recording earlier plagioclase fractionation. Most likely, the infiltrating melts were overall derived from an evolving H2O-rich magma emplaced below the dunite body. The second event of melt–dunite reactive interaction developed the websterite ring around dunites. We envision that the outermost domain of the dunite body was replaced by websterites in response to reaction with an invading H2O-poor melt, which had previously undergone plagioclase fractionation. The dunite replacement occurred under dynamic conditions, which promoted the reaction progress, thereby leading to total or almost total dissolution of precursor olivine, and started to form the lens shape of the Monte Mazzucco ultramafic association. The gabbronorites closely adjacent to the websterite ring represent the crystallization products of the invading melt.

Textural and chemical variations in peridotite induced by hydrous melt migration in mantle under the back-arc spreading

This study presents the structural and petrological characteristics of mantle peridotites in relation to the rock-hydrous melt reaction found in the Hayachine ultramafic complex, NE Japan. We conducted sampling, microstructural observations, crystal-orientation analyses, and major element composition analyses of the major constituent minerals in the peridotites. We used 17 serpentinized peridotites that preserved better mantle textures. The peridotites are composed of lherzolite to harzburgite, consisting of olivine, orthopyroxene, clinopyroxene, amphibole, and spinel. The peridotites were classified into two textures according to olivine grain size: coarse-grained (ca. 2–3 mm) and fine-grained (ca. 0.3–0.5 mm). Fine-grained peridotites were newly found in this study and were characterized by aggregates of orthopyroxene and amphibole with a large number of spinel inclusions. Based on olivine crystal-preferred orientation (CPO), we found that the coarse-grained peridotites were further classified into Group 1 (A type CPO) and 2 (AG type CPO) and the fine-grained peridotites were accordingly classified into Group 3 (weak CPO). The systematic continuity in the chemical compositions of the minerals suggests that Group 1 peridotites partially melted to form Group 2 peridotites, followed by Group 3 peridotites, due to further reaction of Group 2 peridotites with hydrous melts. These textural and chemical variations in the peridotites could have resulted from rock-hydrous melt reactions under back-arc spreading and subsequent processes.

Equilibrium model approach to predict local chemical changes in recovery boiler deposits

A step-by-step equilibrium model approach was developed to estimate and describe the enrichment of K and Cl within recovery boiler deposits. The model predicts free melt movement towards colder temperatures within deposits. Due to a temperature gradient, the melt amount and composition differ across the deposit. The melt movement affects the local composition, which leads to changes in the local phase composition. The model predicts how the deposit profile changes due to melt migration into pores within the deposits. Changes in the deposit composition profile also affect the melting behavior locally, resulting in lower local melting temperatures, which can be detrimental to heat exchanger materials. The modeling results were compared to earlier published laboratory and full-scale boiler measurements, and there is a good agreement between the results. The model predicts a local decrease in the first melting temperature of recovery boiler deposits by ∼30 °C. These findings closely align with experimental results, shedding light on the intricate mechanisms of melt percolation and intra-deposit aging processes. The proposed step-by-step model offers a means to achieve more accurate estimations of locally prevalent first melting temperatures in recovery boiler deposits.

Evaluation of plume-induced continental crust growth rate in early Earth: Insight from integrated petrological-thermo-mechanical modeling

The origin of Earth’s felsic continental crust remains a mystery. The generation of felsic continental crust requires a two-stage partial melting from original mantle sources. There are two hypotheses for the continental crust generation in the early Earth. One is the subduction-related magmatism, e.g., island arcs, which produces intermediate to felsic magma that constitutes the early buoyant continental crust. The other is the magmatism induced by the mantle plume that creates thick basaltic crust and finally the continental crust. However, there is controversy about the origin of plate tectonics, which is an obstacle for simply applying the subduction-induced model in the early Earth. On the other hand, the efficiency of mantle plume-induced continental crust growth remains unknown. In this study, we develop a new numerical model, integrating the petrological-thermo-mechanical model with melt migration and crystallization, to evaluate the efficiency of continental crust production by mantle plumes in Earth’s history. Our results indicate that mantle plumes are considerably more effective for continental crust generation in the hot early Earth than that in the present Earth. The contribution of plume-induced continental crust growth may be greatly promoted by the possible high frequency of mantle plume generation in the early Earth than the present.

Numerical study on migration and solidification behavior of molten droplets on structure surface and in debris voids during debris bed melting process

The melt migration and solidification behavior on the surface of metal structure and in the debris voids is a key phenomenon of debris bed melting, which has an important impact on the assessment of reactor pressure vessel damage and melt leakage. In this study, the moving particle semi-implicit (MPS) method with high-order discretization scheme was coupled with a surface tension model that was established based on interface reconstruction. The migration and solidification behavior of UO2–ZrO2 and Fe–Zr melts on the metal structure surface and in the debris voids were investigated with the improved MPS method. The effects of temperature, velocity, contact angle and droplet diameter on the melt migration and solidification were analyzed. The results showed that the molten UO2–ZrO2 droplet on the metal structure surface presented three stages with solidification rate decreasing, while the solidification rate of molten Fe–Zr droplet had little change. The maximum spread factor and solidification rate increased with the increase of droplet falling velocity, but the spread factor was independent of the melt type at the same velocity. The high-temperature molten droplet penetrated the voids formed by the debris with diameter of 5 mm more easily compared to the voids formed by the debris with diameter of 3 mm. The influence of contact angle on the migration of molten droplet with initial velocity was small, and the maximum difference in droplet mass fraction was about 6%. Three groups of molten droplets with different diameters penetrated the voids, and the average increments of penetration mass fraction were 2.3%, 2.7% and 5.8%, respectively, with the increase of velocity. Near the bottom of the debris bed, the molten droplet without initial velocity solidified and blocked the flow channel, but the molten droplet with initial velocity of 0.5 m/s penetrated through the debris voids. The droplet with initial velocity of 0.5 m/s had a faster solidification rate in the voids compared to the droplet without initial velocity, and the droplet solidification mass fraction was 96.27% when the initial melt temperature was 1410 K.

Hydrous shear zones are sites of melt transfer in the lower arc crust: A case study from Fiordland, New Zealand

AbstractRecent studies of the Cretaceous lower arc crust exposed in Fiordland, New Zealand, conclude that shear zones are sites of melt migration and mass transfer through the deep crust. Here, we investigate the 4–10 km‐wide George Sound Shear Zone, which cuts the Western Fiordland Orthogneiss, comprising two main rock types: two‐pyroxene gneiss and hornblende gneiss. Previous studies infer a predominantly igneous origin for the two types of gneiss. However, this study finds that melt‐rock interaction within the George Sound Shear Zone formed the hornblende gneiss from the precursor two‐pyroxene gneiss. Petrographic analyses of samples collected in transects across the shear zone show hydration reaction textures ranging from rims of hornblende + quartz around pyroxene grains to complete replacement of pyroxene grains. Plagioclase is recrystallized and partially replaced by clinozoisite. Additionally, biotite mode increases towards the higher strain rocks in the shear zone. Backscatter images and polarized light microscopy show microstructures indicative of former melt‐present deformation, including (a) interconnected mineral films of quartz and K‐feldspar along grain boundaries, (b) grains that terminate with low apparent dihedral angles, (c) interstitial grains, with some (d) undulose extinction in plagioclase and (e) serrated grain boundaries. In addition, zircon microstructures are consistent with Zr mobility, further supporting the former presence of melt; geochemical data show enrichment of Zr in the hornblende gneiss as compared to the two‐pyroxene gneiss. From the above observations, it is inferred that a felsic to intermediate hydrous melt migrated through the George Sound Shear Zone reacting with the host two‐pyroxene gneiss of the Western Fiordland Orthogneiss. Melt migration along grain boundaries was deformation assisted, (i) causing hydration of pyroxene to hornblende + quartz, and plagioclase to clinozoisite, (ii) increasing proportions of biotite within the shear zone and (iii) causing depletion of Cr + Ni and Zr enrichment in the hydrated rock. Our interpretation is supported by published observations of pegmatite dyke swarms that intruded into the George Sound Shear Zone, the P‐T conditions of deformation and characterization of microstructures that contrast sharply with those typically found in mylonitic rocks formed under solid‐state metamorphic conditions. Our results confirm that hydrous shear zones within otherwise anhydrous country rock are retrogressive and may represent evidence of melt migration through zones of deformation.

Migration and accumulation of hydrous mantle incipient melt in the Earth's asthenosphere: Constraints from in-situ falling sphere viscometry measurements

The Earth's oceanic lithosphere-asthenosphere boundary (LAB) is marked by a notable decrease in seismic wave velocities and an increase in attenuation. This phenomenon is likely attributed to the accumulation of partial melt at the top of the asthenosphere. Nevertheless, the process involving the upward migration and aggregation of low-degree partial melts, highly likely to be mantle incipient melts in the asthenosphere, remains underexplored. Viscosity is a key factor controlling the flow of melts, thus in this study, we used In-situ X-ray falling sphere viscometry experiments to determine the viscosity of the mantle incipient melt containing 1.8–4.0 wt.% water at 1.5–6 GPa and 2100 K. We found that water and pressure can effectively decrease the viscosity of the mantle incipient melt. Therefore, the mantle incipient melts are highly mobile with a high segregation velocity in the deeper part of the asthenosphere. However, during ascent, their mobility diminishes due to an increase in viscosity and a decrease in water content. These mechanisms potentially contribute to the accumulation of partial melt at the top of the asthenosphere.

Prolonged Multi‐Phase Magmatism Due To Plume‐Lithosphere Interaction as Applied to the High Arctic Large Igneous Province

AbstractThe widespread High Arctic Large Igneous Province (HALIP) exhibits prolonged melting over more than 50 Myr, an observation that is difficult to reconcile with the classic view that large igneous provinces (LIPs) originate from melting in plume heads. Hence, the suggested plume‐related origin and classification of HALIP as a LIP have been questioned. Here, we use numerical models that include melting and melt migration to investigate a rising plume interacting with lithosphere of variable thickness, that is, a basin‐to‐craton setting applicable to the Arctic. Models reveal that melt migration introduces significant spatial and temporal variations in melt volumes and pulses of melt production, including protracted melting for at least about 30–40 Myr, because of the dynamic feedback between migrating melt and local lithosphere thinning. For HALIP, plume material deflected from underneath the Greenland craton can re‐activate melting zones below the previously plume‐influenced Sverdrup Basin after a melt‐free period of about 10–15 Myr, even though the plume is already ∼500 km away. Hence, actively melting zones do not necessarily represent the location of the deeper plume stem at a given time, especially for secondary pulses. Additional processes such as (minor) plume flux variations or local lithospheric extension may alter the timing and volume of HALIP pulses, but are to first order not required to reproduce the long‐lived and multi‐pulse magmatism of HALIP. Since melting zones are always plume‐fed, we would expect HALIP magmatism to exhibit plume‐related trace element signatures throughout time, potentially shifting from mostly tholeiitic toward more alkalic compositions.

Deep versus shallow emplacement of sills and dykes: new insight from thermo-visco-elastic modelling

Rapid emplacement of a mafic dyke or sill at mid-crustal depth heats and possibly melts the felsic wall rock followed by solidification. Associated volume changes generate stresses, possibly enforcing brittle failure and melt migration. We model the evolution of melting, solidification, temperature, and stress including visco-elastic relaxation in 1D - dykes or -sills using realistic rock rheologies of the Weschnitz pluton (Odenwald). For deep emplacement (Case 1, 15.3 km) extensive contact melting of the wall rock occurs, for shallow emplacement (Case 2, 10 km) it is negligible. The stresses are zero at high melt fractions, but increase during solidification and cooling: The intrusion orthogonal stress is always zero. The intrusion parallel stress σ∥ within the intrusion is tensile (O(200 MPa)). It relaxes on a time scale between a few years (Case 1) and 0.6 m.y. (Case 2). Within the wall rock σ∥ is compressive during heating, but becomes tensile under solidification and cooling. Wall rock stresses relax on a time scale of months to 100 years. A Deborah number is defined based on viscous to thermal relaxation allowing generalization of our results. Adding lithostatic stresses, the total stresses of Case 1 remain below the brittle strength, while for Case 2 they may exceed it. Adding the lithostatic pressure to the melt pressure, the effective stresses exceed the brittle strength and intrusion orthogonal tensile fractures are predicted. Combined with the pressure gradient within the over-pressurized felsic melts generated in the wall rock, this explains the migration of felsic contact melt into shrinkage cracks of the mafic sill in the Weschnitz pluton.

Accretion of “young” Phanerozoic subcontinental lithospheric mantle triggered by back-arc extension—the case of the Ivrea-Verbano Zone

The subcontinental lithospheric mantle (SCLM) beneath Phanerozoic regions is mostly constituted by fertile lherzolites, which sharply contrast with cratonic mantle made of highly-depleted peridotites. The question of whether this chemical difference results from lower degrees of melting associated with the formation of Phanerozoic SCLM or from the refertilization of ancient depleted SCLM remains a subject of debate. Additionally, the timing and geodynamic environment of accretion of the fertile SCLM in many Phanerozoic regions are poorly constrained. We here document new geochemical and Nd-Hf isotopic data for orogenic lherzolite massifs from the Ivrea-Verbano Zone (IVZ), Southern Alps. Even though a few Proterozoic Re depletion ages are locally preserved in these mantle bodies, our data reveal that the IVZ lherzolitic massifs were “recently” accreted to the SCLM in the Upper Devonian (ca. 370 Ma) during Pangea amalgamation, with a petrochemical evolution characterized by low-degree (~ 5–12%) depletion and nearly contemporaneous pervasive to focused melt migration. The lithospheric accretion putatively took place through asthenospheric upwelling triggered by Variscan intra-continental extension in a back-arc setting related to the subduction of the Rheic Ocean. We thus conclude that the fertile sections of Phanerozoic SCLM can be accreted during “recent” events of back-arc continental extension, even where Os isotopes preserve memories of melting events in much older times.

Simultaneous rift-scale inflation of a deep crustal sill network in Afar, East Africa

Decades of studies at divergent plate margins have revealed networks of magmatic sills at the crust-mantle boundary. However, a lack of direct observations of deep magma motion limits our understanding of magma inflow from the mantle into the lower crust and the mechanism of sill formation. Here, satellite geodesy reveals rift-scale deformation caused by magma inflow in the deep crust in the Afar rift (East Africa). Simultaneous inflation of four sills, laterally separated by 10s of km and at depths ranging 9–28 km, caused uplift across a ~ 100-km-wide zone, suggesting the sills are linked to a common mantle source. Our results show the supply of magma into the lower crust is temporally episodic, occurring across a network of sills. This process reflects inherent instability of melt migration through porous mantle flow and may be the fundamental process that builds the thick igneous crust beneath magmatic rifts and rifted margins globally.

Constraints on crustal compositional architecture across the North China-Altaids transition and implications for craton margin reworking

Abstract The phase of secular evolution of continents is manifested as the degree of compositional differentiation, modification, and maturation of continental crusts, which is vital in understanding the mechanism of continental evolution but is difficult to quantify. Here we use integrated passive- and active-source seismic profiling to conduct joint analysis and inversion and derive Vs and Vp/Vs section models across the North China craton (NCC) to southeastern Altaids boundary zone that bears a tectonic transition from a reworked ancient craton margin to a Phanerozoic accretionary orogen. We systematically exploited the imaged multiple physical properties as precise and delicate proxies to constrain the compositional architecture in the crust across this important tectonic transition subject to various crustal evolutional phases. Our Vs and Vp/Vs imaging, together with the existing isotopic data, characterizes the Yin Shan-Yan Shan belt as the northern NCC margin with layered homogeneous compositions that point to an evolved crust. However, the lower-crustal low-Vs/high-Vp/Vs signature that overlaps the shallowly dipping to horizontal reflective fabrics suggests the crust of the northern NCC margin has undergone considerable reworking through lower-crustal-stretching assisted melt migration and mixing since the late-Paleozoic to Mesozoic era. The process probably involved crust-mantle interaction and thus resulted in a compositionally modified ancient crustal basement. On the contrary, the southeastern Altaids domain manifests crustal complexity in compositions and structures inferred to be indicative of a juvenile crust of the Phanerozoic accretionary orogen. Our results provide deep physical-property constraints that shed new light on the crustal evolution of a complex craton margin.

Paleoarchean sedimentation and Mesoarchean high-temperature metamorphism in northeastern Brazil: Tracing sources and depositional ages in polymetamorphic rocks

Supracrustal rocks offer a window into tectonic processes of the early Earth, since they are common in the Archean lithosphere. However, these rocks are usually affected by several episodes of metamorphism that can compromise their UPb systematics, leading to equivocal interpretations of depositional ages and sources. In northeast Brazil, supracrustal rocks are frequent within the Archean basement of the São José do Campestre Massif. These metapelitic migmatites show a high-temperature mineral assemblage, with garnet + sillimanite ± spinel and retrograde cordierite, with abundant anatectic melt migration at conditions of upper amphibolite to granulite facies. Zircon UPb dating coupled to trace elements analysis through LA-ICP-MS, as well as zircon internal zoning patterns suggest a maximum depositional age of 3305 ± 16 Ma followed by high-temperature metamorphism in the Mesoarchean, Paleoproterozoic and Neoproterozoic. Mesoarchean high temperature metamorphism occurred between 3084 ± 4 and 3006 ± 6 Ma and generated a wide range of textures that could be grouped in, at least, three stages of zircon growth. The first, during prograde heating, led to dissolution of detrital cores through the process of Ostwald Ripening and reprecipitation in oscillatory zoned rims. The second, probably at peak conditions, occurred above the stability of monazite, as evidenced by high Th/U ratios within zircon grains, rounded shape and sector-zoned cores. The third, during retrograde cooling, is mostly driven by garnet breakdown, and resulted in the crystallization of convoluted rims. Ti-in-zircon temperatures indicate minimum temperatures of 712 ± 21 °C for prograde/retrograde stages and 881 ± 50 °C for peak conditions. The Paleoarchean sedimentation and Mesoarchean metamorphism are coeval with similar events in Kaapvaal and Dharwar Cratons (South Africa and South India, respectively), but show no correlation to any Archean domains in South America to date. Solid-state recrystallization during the Paleoproterozoic (ca. 2.0 Ga) and the Neoproterozoic (ca. 0.6 Ga) correlates with orogenic events in both the Borborema Province and São Francisco Craton, suggesting a common evolution since the Rhyacian.

Thermal and Physical Properties of Barrovian Metamorphic Sequence Rocks in the Ailao Shan‐Red River Shear Zone, and Implications for Crustal Channel Flow

AbstractThe collisional history between Greater India and the Eurasian plate has been well constrained by the study of exhumed Barrovian metamorphic sequence (BMS) rocks in the Himalayan Range. However, in the southeastern Tibetan Plateau, the collisional records have been obscured by intense, regional‐scale strike‐slip overprinting and recrystallization. Here, in BMS rocks from the Ailao Shan–Red River shear zone (ARSZ), we report the first discovery of a >250 km long, high‐pressure (high‐P) granulite belt (>1.0 GPa), identified by the presence of relict kyanite and associated decompression reaction textures. Petrological phase equilibrium modeling showed that exposed micaschists in the region represent exhumed middle crust (20–25 km, 600–670°C), while the high‐P granulite rocks are remnants of thickened lower crust (45–55 km, 800–850°C). This indicates that the northeast edge of the ARSZ experienced an additional ∼25 km of uplift and exhumation compared to the southwest side, facilitated by brittle thrusting/imbrication along the Ailao Shan fault (micaschists) and ductile extrusion along the Red River fault (granulite). Geochronological study shows that the upper portion of the BMS preserves older metamorphic ages (52–34 Ma) than the lower portion (32–29 Ma), which was attributed to spatial variation in cooling rates. Using calculated P–T–t–d paths, we also examined variation in density and seismic wave speeds for BMS in the ARSZ. Our data correlate with fieldwork conducted elsewhere within the Himalayan Range indicating that the kyanite to sillimanite transition zone may serve as a “cap” for the horizontal migration of melt within the lower crust.

Formation of coronae topography and fractures via plume buoyancy and melting

Coronae are volcanic-tectonic features seemingly unique to the surface of Venus. As such they offer potential insights into differences in geodynamic processes on Venus and Earth. Largely due to limitations in processing power and computational complexity, models of corona formation generally ignore the effect of melt and magmatism despite the presence of volcanic features like domes and lava flows. Here we present a new model of corona formation integrating visco-plastic rheology with two-phase flow melt migration. Additionally, we model brittle failure of surface material, viscous flow of partially molten material at depth, and cracking of solid rock from fluid pressure in one model. The model consists of a two-dimensional half space with a round plume body ascending through an asthenosphere and interacting with a Venusian lithosphere and crust. The size and temperature of the plume body, lithosphere thickness, crustal thickness, and regional strain rate are all varied over 20 model runs. Most models show deflection of the plume as it rises due to the increased viscosity of the base of the lithosphere. As the plume spreads and rises against the base of the lithosphere, melting occurs. Melt rises and accumulates in the lithosphere which supports surface topography. Brittle failure occurs at the surface due to the induced stress of the migrating plume body and rising melt. Increased temperatures and sizes of initial plumes result in larger topographies occurring over smaller timescales. Decreased temperatures and sizes of initial plumes result in smaller topographies over larger timescales. Increased lithosphere thickness generates smaller topographies. Thinner crusts result in more viscous lithosphere and more strain accumulation, creating surface depressions seen in approximately half of real corona. Both thicker crusts and increased regional strain result in lower viscosities at the base of the lithosphere, preventing the deflection of the plume body, which instead rises through the lithosphere relatively unimpeded. The varied model parameters generally result in an evolutionary sequence of topographies and fracture patterns which resemble actual corona, without requiring a collapsing dome to create outer rims, as is often invoked for corona formation via upwelling plumes. Corona rim topography is generally created within 105 and 106 year timescales, while multiple rims or central depressions are generally observed over 106 year timescales. Observed topography is actively supported by partially molten regions in the lithosphere, supporting the notion of an active squishy-lid Venus. These partially molten regions are not centrally located and may provide a source for observed lava flows away from corona centers.

Peraluminous Rare Metal Granites in Iberia: Geochemical, Mineralogical, Geothermobarometric, and Petrogenetic Constraints

The intensive variables, geochemical, mineralogical, and petrogenetic constraints of the Iberian peraluminous rare metal granites (RMGs), many of them unknown, are presented. The mineral chemistry of ore and gangue minerals, whole rock analyses, geothermobarometry, melt water and phosphorus contents, mass balance, and Rayleigh modeling were performed to achieve these objectives. These procedures allow us to distinguish two main contrasting granitic types: Nb-Ta-rich and Nb-Ta-poor granites. The former have lower crystallization temperatures, higher water contents, and lower emplacement pressures than Nb-Ta-poor granites. Nb-Ta-rich granites also have higher fluoride contents, strong fractionation into geochemical twins, higher Na contents, and different evolutionary trends. At the deposit scale, the fractional crystallization of micas properly explains the variation in the Ta/Nb ratio in both Nb-Ta-poor and Nb-Ta-rich RMGs, although in higher-grade granites, the variation is not as clear due to the action of fluids. Fluid phase separation processes especially occurred in the Nb-Ta rich granites, thus transporting halogens and metals that increased the grades in the top and sometimes in the core of granites. Gas-driven filter pressing processes facilitated the migration of fluid and melt near solidus melt in Nb-Ta-rich granites. The geochemical signature of the Iberian rare metal granites mainly follows the trends of two-mica granites and P-rich cordierite granites, but also of granodiorites.

Evidence of Topographic Change Recorded by Lava Flows at Atete and Aruru Coronae on Venus

AbstractCoronae are quasi‐circular volcano‐tectonic features on Venus. Four critical observations have been identified within the coronae population, including fracture annuli, a wide range of diameters, complex and varied topography, and various types of associated volcanic features. Geophysical models have attempted to replicate their formation from a variety of lithospheric processes but struggle to recreate all four of the critical observations. Volcanism is an often overlooked characteristic in corona formation models. Paleotopographic techniques are applied to lava flows associated with Atete Corona and Aruru Corona in the Beta‐Atla‐Themis region to investigate post‐emplacement changes in topography. Our results indicate marked divergence between lava flow orientation and the modern slope within the fracture annuli of the coronae. Intra‐annular flows at Atete Corona were emplaced on a surface that was reoriented by up to ∼180°. Intra‐annular flows at Aruru Corona were emplaced on a surface that was reoriented between 90° and 145°. Lava flows on the flanks of both coronae diverge less than those within the fractures. This finding suggests a progression of volcanism that starts at the interior of the corona and migrates outward to the fracture annuli. The role of volcanism, both intrusive and extrusive, is likely to play a more substantial role in the corona formation. Incorporating melt migration into geophysical models could significantly enhance our comprehension of the processes underlying the formation of coronae.

The critical role of deformation-assisted melt migration in the formation of oceanic core complexes

Oceanic core complexes provide an accessible window into deep processes occurring at slow and ultra-slow-spreading mid-ocean ridges. We analyse samples from IODP ocean drilling of core complexes at the Atlantis Bank, Atlantis Massif, and near the Kane Transform at the South West Indian and Mid-Atlantic ridges. We correlate secondary minerals, including oxides, with sites of melt migration. We interpret changes to mineral assemblage and microchemistry, reaction textures and melt-pseudomorph microstructures as fingerprints of open-system melt-mediated processes. This micro-scale information is combined with a macro-scale review of legacy mineral chemistry data to show that melt-fluxed rocks share remarkably similar characteristics across the three core complexes investigated. These are rich in oxides and have olivine, orthopyroxene and clinopyroxene grains that are chemically distinct from oxide-poor gabbros. We propose that oceanic crust fluxed with external melt can be recognised by the following key features: (1) high modes of secondary minerals, such as oxides and olivine, (2) microstructural evidence for the former presence of melt, and (3) mineral chemistry differences between primary and secondary olivine, orthopyroxene, clinopyroxene and plagioclase. Importantly, olivine has previously only been reported as primary, that is, having crystallised from magma. However, our results show that gabbros with secondary olivine are reliable indictors of melt–rock interaction during deformation-assisted diffuse melt migration through the gabbroic oceanic crust. Finally, we propose a new model for the formation of oceanic core complexes where deformation-assisted melt migration plays a critical role in strain localisation, exhumation and evolution of the core complex.

A Scaling for the Permeability of Loose Magma Mush Validated Using X‐Ray Computed Tomography of Packed Confectionary in 3D and Estimation Methods From 2D Crystal Shapes

AbstractMelt percolation through partially molten “mushy” regions of the crust underpins models for magma migration, accumulation, and processes that prime systems for eruption. Knowledge of the hydraulic properties of magma mush, specifically permeability, is required for accurate predictions of melt migration rates and accumulation timescales. Previous studies, validated for cuboidal crystal analogs, show that crystal shape exerts a first‐order control on the permeability, and is tested here for anisometric natural crystal shapes using X‐ray CT 3D data sets of magma mush analogs made from packed confectionary particles arranged randomly. We use a lattice‐Boltzmann fluid flow simulation tool to determine the permeability of the analogue melt phase network between the packed particles. We find excellent agreement with our data sets to within ∼0.1 log units, when the specific surface area is measured. To extend this to more typical cases where the specific surface area is unknown, we use the shape and size of the objects determined in both 3D and 2D to estimate the specific surface area assuming a cuboid approximation. These approximate solutions give good results to within ∼0.5 log units of the measured permeability and offer a method by which permeability could be estimated from a thin section of a cumulate or pluton sample. Our shape‐sensitive approach is more accurate than existing models for permeability of magma mush, most approximating natural crystal shapes to spheres. We therefore propose that these could be implemented in dynamic magma mush models for melt movement in the crust to produce more accurate flux predictions.