Geochemical model of precambrian granitoid magmatic evolution in the Korosten pluton (ukrainian shield): petrogenetic aspects and genesis of complex ore mineralization in metasomatic zones
Preliminary geochemical dataset for the granitoids (rapakivi, granite-porphyries, veined granites) of Korosten anorthosite-rapakivi granite pluton (Ukrainian Shield) was studied. The data set includes the results of all major and selected trace elementsdeterminations (XRF, >300 samples), and is compared with melt crystallization and partial melting models of trace elements behavior. Only the Rayleigh model closely approximates the trace element data for granitoids. Typical incompatible behavior under approximately constant bulk distribution coefficient is determined for Rb (DRb = 0,5). Model f values (weight fraction of liquid in magma chamber) are calculated for each granitoid type (residual melt portion) from Rayleigh equation and Rb content in rocks (CRb) assuming minimum concentration in granitoids (169 ppm) as Rb content in parent magma (C0 Rb). C vs. f curves for trace (including P, Ti, S, Cl, F, Ca) and major elements are approximated by means of equations of C?C0?fD-1 form or polynomial onesrespectively. This set of equations is an idealized model of element behavior that demonstrates the depletion of Ba, Sr, Ti, Zr, P, S and enrichment of Th and Ga in the residual liquid as well as the inversion behavior of LREE, Y, F, Cl, Nb, Zn and Pb &during the melt fractional crystallization in magma chamber. Monotonous decrease of both Zr and P content indicates melt saturation in zircon and apatite. Therefore, model temperature (Tmodel) of the melt is estimated from equations of
zircon and ap&atite solubility. The temperature evolution in magma chamber was presented as Tmodel vs. f polynomial equation (Tmodel range: 900-720oC). Inversion in LREE content (f = 0,185) indicates the apatite/monazite replacement in the crystallizing assemblage&. Water content in melt for this f value and corresponding Tmodel was calculated from equation of monazite solubility which demonstrates its high H2Odependence. C0 H2O=2,36 wt% was estimated on this basis (assuming DH2O = 0,1) for the liquidus of ini&tial granite melt (Ptotal ~6,3 kbar). According to designed model water saturation limit was reached at f = 0,165 and H2O-fluid was extracted from the melt during its further evolution. Synchronous inversion of C/C0 vs. f curves shows the enrichment &of fluid with F, Cl, Nb, Zn, Pb etc. Such characteristics of fluid are in agreement with geochemical data for ore-bearing altered rocks and might testify for their genetic unity.
Було вивчено попередній набір геохімічних даних для гранітоїдів (ра&паківі, граніт-порфіри та жильні граніти) Коростенського анортозит-рапаківігранітного плутону (КП) Українського щита. Дані щодо мікроелементного складу гранітоїдів достатньо близько апроксимуються моделлю фракційної кристалізації Релея. Для Rb встано&влено типову несумісну поведінку з постійним комбінованим коефіцієнтом розподілу (DRb = 0,5), що дозволило провести розрахунок значення f (масової частки рідкої фази в магматичній камері) для кожного типу гранітоїдів (порції залишкового розплаву) при&ймаючи мінімальний вміст елемента в гранітоїдах (169 ррт) за його вміст у материнському розплаві (C0Rb). Залежність концентрацій від f для мікро- (включаючи P, Ti, S, Cl, F, Ca) і петрогенних елементів апроксимовані рівняннями виду C=C0 f D-1 або пол&іноміальними, відповідно. Одержаний набір рівнянь є ідеалізованою моделлю поведінки елементів, яка ілюструє характер еволюції розплаву в процесі фракційної кристалізації: виснаження Ba, Sr, Ti, Zr, P, S; збагачення Th, Ga; інверсію поведінки LREE, Y,& F, Cl, Nb, Zn и Pb. Монотонне зниження концентрації Zr і Р вказує на насичення розплаву цирконом та апатитом, що дозволило провести розрахунок модельних температур (Tmodel) користуючись експериментальними рівняннями розчинності циркону і апатиту в с&илікатних розплавах, а також отримати поліноміальне рівняння залежності Tmodel від f (діапазон Tmodel: 900-720oC). Вміст води у розплаві для моменту заміни апатиту монацитом у продуктах кристалізації (інверсія поведінки LREE; f = 0,185) розраховано з&