Numerical modelling of fluid-melt partitioning in porphyry systems: configuration, post-processing and datasets

This data publication supports the numerical study “Timescales of metal fluxes in porphyry copper systems from coupled physicochemical processes of magmas, rocks and fluids’’, which investigates the interaction of physicochemical processes driving porphyry copper systems. The coupled magmatic-hydrothermal model uses an implementation of the modelling platform CSMP++ and the algebraic multi-grid solver SAMG. A full description of the software architecture is provided in Gruzdeva et al. (2024a).

As a further development of the model of Gruzdeva et al. (2024b), the study incorporates fluid-melt partitioning of chemical components, using constant distribution coefficients D to analyse the partitioning of copper and salt during magma crystallization. The models are configured to represent a 2D cross-section of the upper crust, incorporating an initial magma reservoir at variable depths. To consider fluid evolution in 3D, we use radial coordinates in respect to the vertical centre axis of the model. This approach allows for the quantification of mass balance in 3D when representing the fluid dynamics of the magmatic-hydrothermal system. The location of the cupola of a magma reservoir is a critical parameter in the study of PCDs, in particular in relation to the dynamic of fluid release and mineralization processes. In the study, simulations are conducted with cupolas at a depth of 4 and 6 km, respectively.

The dataset is designed to facilitate the analysis of numerical simulations focusing on the partitioning behaviour of melts and hydrothermal fluids in ore-forming environments in porphyry copper systems. The dataset is stored in a structured format within the repository provided below, ensuring reproducibility and accessibility for further research.