A high-order discretization consisting of a tensor product of the Fourier collocation and discontinuous Galerkin methods is presented for numerical modeling of magma dynamics.The physical model is an advection-reactio...A high-order discretization consisting of a tensor product of the Fourier collocation and discontinuous Galerkin methods is presented for numerical modeling of magma dynamics.The physical model is an advection-reaction type system consisting of two hyperbolic equations and one elliptic equation.The high-order solution basis allows for accurate and efficient representation of compaction-dissolution waves that are predicted from linear theory.The discontinuous Galerkin method provides a robust and efficient solution to the eigenvalue problem formed by linear stability analysis of the physical system.New insights into the processes of melt generation and segregation,such as melt channel bifurcation,are revealed from two-dimensional time-dependent simulations.展开更多
Discontinuous Galerkin(DG)and matrix-free finite element methods with a novel projective pressure estimation are combined to enable the numerical modeling of magma dynamics in 2D and 3D using the library deal.II.The p...Discontinuous Galerkin(DG)and matrix-free finite element methods with a novel projective pressure estimation are combined to enable the numerical modeling of magma dynamics in 2D and 3D using the library deal.II.The physical model is an advection-reaction type system consisting of two hyperbolic equations to evolve porosity and soluble mineral abundance at local chemical equilibrium and one elliptic equation to recover global pressure.A combination of a discontinuous Galerkin method for the advection equations and a finite element method for the elliptic equation provide a robust and efficient solution to the channel regime problems of the physical system in 3D.A projective and adaptively applied pressure estimation is employed to significantly reduce the computational wall time without impacting the overall physical reliability in the modeling of important features of melt segregation,such as melt channel bifurcation in 2D and 3D time dependent simulations.展开更多
The Taiwan mountain belt, one of the youngest orogenies in the world, is caused by the collision of the Luzon arc with the Eurasian margin, which leads to post-collisional extension and magmatism in the Northern Taiwa...The Taiwan mountain belt, one of the youngest orogenies in the world, is caused by the collision of the Luzon arc with the Eurasian margin, which leads to post-collisional extension and magmatism in the Northern Taiwan Volcanic Zone(NTVZ). The magma chamber process in this region has not previously been elucidated in detail. In this paper, the textural and compositional features of plagioclase phenocrysts in basalt from the Tatun Volcanic Group(TTVG) were studied to restrict the dynamics of magma system. Results show that the magma melts in TTVG are mainly sourced from the underlying MORB-like mantle wedge but influenced by incorporation of subduction components, causing the elevated Sr/Y and Ba/Y ratios in magma melts. The subduction components are mainly transported in the form of sediment melt. The plagioclase phenocrysts in the TTVG volcanic rocks are generally coarsely core-sieved with a clear rim. The An contents in the rims of plagioclase are much lower than those of cores, and elevated FeO concentrations are detected in the plagioclase rims. We propose there exists a double-layer magma chamber in this region. The core of the plagioclase was crystalized in the deeper quiescent magma chamber(~21 km), which was subsequently partially dissolved during the ascent of magma melt under H_(2)O-undersaturated condition, forming the typical coarsely sieved texture and synneusis. When this crystal-rich melt migrates into the shallower chamber, water saturation is reached and more sodic plagioclase formed as the rim of phenocryst. Due to the considerably higher fO_(2) in the shallow chamber than in the deeper one, and the distribution of Fe between plagioclase and melt positively correlates with fO_(2), the FeO content in the plagioclase rim elevates in conjunction with increasing fO_(2).展开更多
The Jinchuan Ni-Cu sulfide deposits,NW China,are hosted in small ultramafic intrusions that were emplaced into Paleoproterozoic metamorphic rocks. The ultramafic intrusions were previously thought to be the segments o...The Jinchuan Ni-Cu sulfide deposits,NW China,are hosted in small ultramafic intrusions that were emplaced into Paleoproterozoic metamorphic rocks. The ultramafic intrusions were previously thought to be the segments of a single elongate intrusion that was dismembered by late faults into eastern and western portions,each of which have distinct stratigraphic sequences.展开更多
We review recent progress in studying silicate, carbonate, and metallic liquids of geological and geophysical importance at high pressure and temperature, using the large-volume high-pressure devices at the third-gene...We review recent progress in studying silicate, carbonate, and metallic liquids of geological and geophysical importance at high pressure and temperature, using the large-volume high-pressure devices at the third-generation synchrotron facility of the Advanced Photon Source, Argonne National Laboratory. These integrated high-pressure facilities now offer a unique combination of experimental techniques that allow researchers to investigate structure, density, elasticity, viscosity, and interfacial tension of geo-liquids under high pressure, in a coordinated and systematic fashion. Experimental techniques are described, along with scientific highlights. Future developments are also discussed.展开更多
文摘A high-order discretization consisting of a tensor product of the Fourier collocation and discontinuous Galerkin methods is presented for numerical modeling of magma dynamics.The physical model is an advection-reaction type system consisting of two hyperbolic equations and one elliptic equation.The high-order solution basis allows for accurate and efficient representation of compaction-dissolution waves that are predicted from linear theory.The discontinuous Galerkin method provides a robust and efficient solution to the eigenvalue problem formed by linear stability analysis of the physical system.New insights into the processes of melt generation and segregation,such as melt channel bifurcation,are revealed from two-dimensional time-dependent simulations.
基金This work was supported by NSF grant OCE-1156706.
文摘Discontinuous Galerkin(DG)and matrix-free finite element methods with a novel projective pressure estimation are combined to enable the numerical modeling of magma dynamics in 2D and 3D using the library deal.II.The physical model is an advection-reaction type system consisting of two hyperbolic equations to evolve porosity and soluble mineral abundance at local chemical equilibrium and one elliptic equation to recover global pressure.A combination of a discontinuous Galerkin method for the advection equations and a finite element method for the elliptic equation provide a robust and efficient solution to the channel regime problems of the physical system in 3D.A projective and adaptively applied pressure estimation is employed to significantly reduce the computational wall time without impacting the overall physical reliability in the modeling of important features of melt segregation,such as melt channel bifurcation in 2D and 3D time dependent simulations.
基金funded by the Shandong Provincial Natural Science Foundation, China (No. ZR2019BD010)the National Natural Science Foundation of China (No. 42276085)+2 种基金Open Fund of the Key Laboratory of Marine Geology and Environment, Chinese Academy of Sciences (No. MGE2021KG01)Pilot National Laboratory for Marine Science and Technology (Qingdao) (No. JCZX202014)the Key Laboratory Submarine Geoscience and Prospecting Techniques, Ministry of Education, Ocean university of China。
文摘The Taiwan mountain belt, one of the youngest orogenies in the world, is caused by the collision of the Luzon arc with the Eurasian margin, which leads to post-collisional extension and magmatism in the Northern Taiwan Volcanic Zone(NTVZ). The magma chamber process in this region has not previously been elucidated in detail. In this paper, the textural and compositional features of plagioclase phenocrysts in basalt from the Tatun Volcanic Group(TTVG) were studied to restrict the dynamics of magma system. Results show that the magma melts in TTVG are mainly sourced from the underlying MORB-like mantle wedge but influenced by incorporation of subduction components, causing the elevated Sr/Y and Ba/Y ratios in magma melts. The subduction components are mainly transported in the form of sediment melt. The plagioclase phenocrysts in the TTVG volcanic rocks are generally coarsely core-sieved with a clear rim. The An contents in the rims of plagioclase are much lower than those of cores, and elevated FeO concentrations are detected in the plagioclase rims. We propose there exists a double-layer magma chamber in this region. The core of the plagioclase was crystalized in the deeper quiescent magma chamber(~21 km), which was subsequently partially dissolved during the ascent of magma melt under H_(2)O-undersaturated condition, forming the typical coarsely sieved texture and synneusis. When this crystal-rich melt migrates into the shallower chamber, water saturation is reached and more sodic plagioclase formed as the rim of phenocryst. Due to the considerably higher fO_(2) in the shallow chamber than in the deeper one, and the distribution of Fe between plagioclase and melt positively correlates with fO_(2), the FeO content in the plagioclase rim elevates in conjunction with increasing fO_(2).
文摘The Jinchuan Ni-Cu sulfide deposits,NW China,are hosted in small ultramafic intrusions that were emplaced into Paleoproterozoic metamorphic rocks. The ultramafic intrusions were previously thought to be the segments of a single elongate intrusion that was dismembered by late faults into eastern and western portions,each of which have distinct stratigraphic sequences.
基金support from the National Science Foundation (Nos. EAR-0001088, 0711057, and 1214376)Guoyin Shen acknowledges support from the DOE (Nos. DE-NA0001974 and DE-FG02-99ER45775)+5 种基金COMPRES for the support in developing the PEP system. Portions of this work were performed at GeoS oilE nviroC ARS (Sector 13), Advanced Photon Source (APS), Argonne National LaboratoryGeo Soil Enviro CARS is supported by the National Science Foundation-Earth Sciences (No. EAR-1128799)Department of Energy-Geo Sciences (No. DE-FG02-94ER14466)HPCAT operations are supported by DOE-NNSA under Award (Nos. DE-NA0001974)DOE-BES under Award (No. DE-FG02-99ER45775), with partial instrumentation funding by NSFUse of the Advanced Photon Source was supported by the U. S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract (No. DE-AC02-06CH11357)
文摘We review recent progress in studying silicate, carbonate, and metallic liquids of geological and geophysical importance at high pressure and temperature, using the large-volume high-pressure devices at the third-generation synchrotron facility of the Advanced Photon Source, Argonne National Laboratory. These integrated high-pressure facilities now offer a unique combination of experimental techniques that allow researchers to investigate structure, density, elasticity, viscosity, and interfacial tension of geo-liquids under high pressure, in a coordinated and systematic fashion. Experimental techniques are described, along with scientific highlights. Future developments are also discussed.
