Major and trace element and Sr–Nd isotope signatures of lavas from the Central Lau Basin: Implications for the nature and influence of subduction components in the back-arc mantle

New major and trace element and Sr–Nd isotope data are presented for basaltic glasses from active spreading centers (Central Lau Spreading Center (CLSC), Relay Zone (RZ) and Eastern Lau Spreading Center (ELSC)) in the Central Lau Basin, SW Pacific. Basaltic lavas from the Central Lau Basin are mainly tholeiitic and are broadly similar in composition to mid-ocean ridge basalts (MORB). Their generally high ⁸⁷Sr/⁸⁶Sr ratios, combined with relatively low ¹⁴³Nd/¹⁴⁴Nd ratios are more akin to MORB from the Indian rather than Pacific Ocean. In detail, the CLSC, RZ and ELSC lavas are generally more enriched in large ion lithophile elements (Rb, Ba, Sr, and K) than average normal-MORB, which suggests that the mantle beneath the Central Lau Basin was modified by subducted slab-derived components. Fluid mobile/immobile trace element and Sr – Nd isotope ratios suggest that the subduction components were essentially transferred into the mantle via hydrous fluids derived from the subducted oceanic crust; contributions coming from the subducted sediments are minor. Compared to CLSC lavas, ELSC and RZ lavas show greater enrichment in fluid mobile elements and depletion in high field strength elements, especially Nb. Thus, with increasing distance away from the arc, the influence of subduction components in the mantle source of Lau Basin lavas diminishes. The amount of hydrous fluids also influences the degree of partial melting of the mantle beneath the Central Lau Basin, and hence the degree of melting also decreases with increasing distance from the arc.     

Nb/Ta,Zr/Hf and REE in the depleted mantle: implications for the differentiation history of the crust–mantle system

High-precision Nb, Ta, Zr, Hf, Sm, Nd and Lu concentration data of depleted mantle rocks from the Balmuccia peridotite complex (Ivrea Zone, Italian Alps) were determined by isotope dilution using multiple collector inductively coupled plasma mass spectrometry (MC-ICPMS) and thermal ionisation mass spectrometry (TIMS). The Zr/Hf ratios of all investigated samples from the Balmuccia peridotite complex are significantly lower than the chondritic value of 34.2, and the most depleted samples have Zr/Hf ratios as low as 10. Correlated Zr/Hf ratios and Zr abundances of the lherzolites preserve the trend of a mantle residue that has been depleted by fractional melting. This trend confirms experimental studies that predict Hf to behave more compatibly than Zr during mantle melting. Experimentally determined partition coefficients imply that the major Zr and Hf depletion most likely occurred in the spinel stability field, with (D_Zr/D_Hf)^cpx≈0.5, and not in the garnet stability field, where (D_Zr/D_Hf)^grt is probably close to one. However, minor amounts of melting must have also occurred in a garnet facies mantle, as indicated by low Sm/Lu ratios in the Balmuccia peridotites. The Nb/Ta ratios of most lherzolites are subchondritic and vary only from 7 to 10, with the exception of three samples that have higher Nb/Ta ratios (18–24). The overall low Nb/Ta ratios of most depleted mantle rocks confirm a higher compatibility of Ta in the mantle. The uniform Nb/Ta ratios in most samples imply that even in ‘depleted’ mantle domains the budget of the highly incompatible Nb and Ta is controlled by enrichment processes. Such a model is supported by the positive correlation of Zr/Nb with the Zr concentration. However, the overall enrichment was weak and did barely affect the moderately incompatible elements Zr and Hf. The new constraints from the partitioning behaviour of Zr–Hf and Nb–Ta provide important insights into processes that formed the Earth’s major silicate reservoirs. The correlation of Zr/Hf and Sm/Nd in depleted MORB can be assigned to previous melting events in the MORB source. However, such trends were unlikely produced during continental crust formation processes, where Sm/Nd and Zr/Hf are decoupled. The different fractionation behaviour of Zr/Hf and Sm/Nd in the depleted mantle (correlated) and the crust (decoupled) indicates that crustal growth by a simple partial melting process in the mantle has little effect on the mass budget of LREE and HFSE between crust and mantle. A more complex source composition, similar to that of modern subduction rocks, is needed to fractionate the LREE, but not Zr/Hf and the HREE.     

Sr-Nd-Pb isotopes of the Early Paleozoic maficultramafic dykes and basalts from South Qinling belt and their implications for mantle composition

Late Early Paleozoic mafic-ultramafic dykes and volcanic rocks from the South Qinling belt are char- acterized by εNd( t ) = 3.28― 5.02, (87Sr/86Sr)i= 0.70341― 0.70555, (206Pb/204Pb)i = 17.256― 18.993, (207Pb/204Pb)i= 15.505―15.642, (208Pb/204Pb)i=37.125―38.968, ?8/4=21.18―774.43, ?7/4=8.11―18.82. These charac- teristics suggest that they derived from a Middle Neoproterozoic mantle with isotopic compositions of mixed HIMU, EMII and minor EMI components. We interpret that these rocks were melting products of depleted mantle modified by subducted ancient oceanic crust and continental margin sediments along the northern margin of Yangtze block during Early Neoproterozoic.     

Adakitic rocks associated with the Shilu copper–molybdenum deposit in the Yangchun Basin,South China,and their tectonic implications

South China is famous for the extensive magmatism and polymetallic mineralization that took place there in the Mesozoic. Shilu is a large porphyry–skarn Cu–Mo deposit in the Yangchun Basin, South China. The lithology of the Shilu intrusion is granodiorite and quartz diorite, both of which are high-K calc-alkaline series, with high Sr([400 ppm) content along with low Y and Yb contents. Most of the samples have characteristics of adakite except for a few samples that have slightly higher Y and Yb contents, which may be plausibly explained by crustal contamination. Laser Ablation Inductively Coupled Plasma Mass Spectrometry zircon U–Pb dating revealed ages between 106.6 ± 1.3 and 103.9 ± 0.5 Ma, with multiple magmatic pulses. Molybdenite Re–Os isochron age of 102.2 ± 2.9 Ma(MSWD = 9.4) was determined, which is identical to the youngest zircon U–Pb age(103.9 ± 0.5 Ma) within error.The Shilu intrusion has high oxygen fugacity as indicated by high zircon Ce~(4+)/Ce~(3+) and Eu_N/Eu_N* ratios. Considering the geochemical characteristics(high Sr, and low Y and Yb contents), high oxygen fugacity, and copper mineralization of the Shilu intrusion, it was most likely formed by partial melting of a subducted young oceanic slab. Whole-rock Sr–Nd isotope-, zircon Hf isotope-, and whole-rock trace element analyses show that Shilu adakitic magmas may have interacted with type II enriched mantle and/or crustal materials during ascent. South China was affected by the Pacific tectonic regime to the east and the Neo-Tethys tectonic regime to the south in the Cretaceous. Based on the Pacific Plate drifting and rotation history, it is hard to explain how the Pacific Plate would have subducted and melted, forming adakitic rocks in the Shilu region. Considering the tectonic history of Southeast Asia and the South China Sea, the Neo-Tethys trench should have been much closer to the South China Block in the Cretaceous, and thus have had a greater impact on the South China Block. Based on the subduction direction, time of subduction,and distance between the Neo-Tethys subduction zone and the Shilu deposit, subduction of the Neo-Tethys ridge is the best mechanism for explaining the Shilu adakitic rocks and Cu–Mo mineralization.     

On the Vp/Vs–Mg# correlation in mantle peridotites: Implications for the identification of thermal and compositional anomalies in the upper mantle

We use thermodynamically self-consistent and hybrid methods to analyze the correlation of important physical parameters (e.g. bulk density, elastic moduli) with bulk Mg# and modal composition in mantle peridotites at upper mantle conditions. Temperature (anharmonic and anelastic), pressure and compositional derivatives for all these parameters are evaluated. The results show that the widely used correlations between Vp/Vs and Mg# in peridotites are strictly valid only for garnet-bearing assemblages at temperatures < 900 °C. The correlation breaks down when: i) spinel is the stable Al-rich phase in the assemblage and ii) when anelastic attenuation of seismic velocities becomes important (T ? 900 °C). This implies that the range of applicability of published Vp/Vs–Mg# correlations for the upper mantle is limited to a depth interval between the spinel–garnet phase transition and the 900 °C isotherm. We use numerical simulations to show that this depth interval is virtually nonexistent in lithospheres thinner than ～ 140 km and can comprise up to ～ 50% of the lithospheric mantle in thick (> 220 km) lithospheric domains. In addition, we show that for most of the upper mantle the expected Δ(Vp/Vs) values associated with compositional variations are smaller than the resolution limit of current seismological methods. All these considerations suggest that the Vp/Vs ratio is not a reliable measure of compositional variations and that for large parts of the upper mantle compositional anomalies cannot be separated from thermal anomalies on the basis of seismological studies only. We further confirm that the only reliable indicator of compositional anomalies in a peridotitic mantle is the ratio of density to shear wave velocities (ρ/Vs). Our results demonstrate that geophysical–petrological models (forward or inverse) that model these two fields (i.e. density and Vs) self-consistently within a robust thermodynamic framework are necessary for characterizing the small-scale thermal and compositional structure of the lithosphere and sublithospheric upper mantle.     

Plume capture by divergent plate motions: implications for the distribution of hotspots,geochemistry of mid-ocean ridge basalts,and estimates of the heat flux at the core–mantle boundary

The coexistence of stationary mantle plumes with plate-scale flow is problematic in geodynamics. We present results from laboratory experiments aimed at understanding the effects of an imposed large-scale circulation on thermal convection at high Rayleigh number (10⁶≤Ra≤10⁹) in a fluid with a temperature-dependent viscosity. In a large tank, a layer of corn syrup is heated from below while being stirred by large-scale flow due to the opposing motions of a pair of conveyor belts immersed in the syrup at the top of the tank. Three regimes are observed, depending on the ratio V of the imposed horizontal flow velocity to the rise velocity of plumes ascending from the hot boundary, and on the ratio λ of the viscosity of the interior fluid to the viscosity of the hottest fluid in contact with the bottom boundary. When V≪1 and λ≥1, large-scale circulation has a negligible effect on convection and the heat flux is due to the formation and rise of randomly spaced plumes. When V>10 and λ>100, plume formation is suppressed entirely, and the heat flux is carried by a sheet-like upwelling located in the center of the tank. At intermediate V, and depending on λ, established plume conduits are advected along the bottom boundary and ascending plumes are focused towards the central upwelling. Heat transfer across the layer occurs through a combination of ascending plumes and large-scale flow. Scaling analyses show that the bottom boundary layer thickness and, in turn, the basal heat flux q depend on the Peclet number, Pe, and λ. When λ>10, q∝Pe^1/2 and when λ→1, q∝(Peλ)^1/3, consistent with classical scalings. When applied to the Earth, our results suggest that plate-driven mantle flow focuses ascending plumes towards upwellings in the central Pacific and Africa as well as into mid-ocean ridges. Furthermore, plumes may be captured by strong upwelling flow beneath fast-spreading ridges. This behavior may explain why hotspots are more abundant near slow-spreading ridges than fast-spreading ridges and may also explain some observed variations of mid-ocean ridge basalt (MORB) geochemistry with spreading rate. Moreover, our results suggest that a potentially significant fraction of the core heat flux is due to plumes that are drawn into upwelling flows beneath ridges and not observed as hotspots.     

Refining locations of the 2005–2006 recurrent earthquakes in Mukacheve,West Ukraine,and implications for their source mechanism and the local tectonics

Three groups of recurrent events have been identified within the 2005–2006 Mukacheve series of about forty M _SH 0.3–1.8 near-surface earthquakes. After a cross-correlation analysis of their waveforms and subsequent relocation with application of the adjusted arrivals and station terms it turned out that they most likely belong to at least two distinct crustal layers separated by an interface at about 4 km; the seismic activation started deeper and progressed upwards. As concerns polarities of first arrivals, the same source mechanism was inferred within individual groups, indicating a normal faulting along the Az ≈ 120° with a slight component of right lateral slip; this finding is consistent with data available on some other local sources, as well as with an ambient field of extensional stresses in the local crust and a significant horizontal gradient of vertical movements across the epicenter area of the Mukacheve series.     

The nature and formation of cristobalite at the Soufrière Hills volcano, Montserrat: implications for the petrology and stability of silicic lava domes

Cristobalite is commonly found in the dome lava of silicic volcanoes but is not a primary magmatic phase; its presence indicates that the composition and micro-structure of dome lavas evolve during, and after, emplacement. Nine temporally and mineralogically diverse dome samples from the Soufrière Hills volcano (SHV), Montserrat, are analysed to provide the first detailed assessment of the nature and mode of cristobalite formation in a volcanic dome. The dome rocks contain up to 11 wt.% cristobalite, as defined by X-ray diffraction. Prismatic and platy forms of cristobalite, identified by scanning electron microscopy (SEM), are commonly found in pores and fractures, suggesting that they have precipitated from a vapour phase. Feathery crystallites and micro-crystals of cristobalite and quartz associated with volcanic glass, identified using SEM-Raman, are interpreted to have formed by varying amounts of devitrification. We discuss mechanisms of silica transport and cristobalite formation, and their implications for petrological interpretations and dome stability. We conclude: (1) that silica may be transported in the vapour phase locally, or from one part of the magmatic system to another; (2) that the potential for transport of silica into the dome should not be neglected in petrological and geochemical studies because the addition of non-magmatic phases may affect whole rock composition; and (3) that the extent of cristobalite mineralisation in the dome at SHV is sufficient to reduce porosity—hence, permeability—and may impact on the mechanical strength of the dome rock, thereby potentially affecting dome stability.     

U–Pb ages of zircons from metamorphic rocks in the upper sequence of the Hidaka Metamorphic Belt,Hokkaido, Japan: Identification of two metamorphic events and implications for regional tectonics

The Hidaka Metamorphic Belt is a well-known example of island-arc crustal section, in which metamorphic grade increases westwards from unmetamorphosed sediment up to granulite facies. It is divided into lower (granulite to amphibolite facies) and upper (amphibolite to greenschist facies) metamorphic sequences. The metamorphic age of the belt was considered to be ~55 Ma, based on Rb – Sr whole-rock isochron ages for granulites and related S-type tonalities. However, zircons from the granulites in the lower sequence yield U – Pb ages of ~21 – 19 Ma, and a preliminary report on zircons from pelitic gneiss in the upper sequence gives a U – Pb age of ~40 Ma. In this paper we provide new zircon U – Pb ages from two pelitic gneisses in the upper sequence to assess the metamorphic age and also the maximum depositional age of the sedimentary protolith. The weighted mean ²⁰⁶Pb/²³⁸U ages from a biotite gneiss in the central area of the belt yield 39.6 ± 0.9 Ma for newly grown metamorphic rims and 53.1 ± 0.9 Ma for the youngest detrital cores. The ages of zircons from a cordierite–biotite gneiss in the southern area are 35.9 ± 0.7 Ma for metamorphic rims and 46.5 ± 2.8 Ma for the youngest detrital cores. These results indicate that metamorphism of the upper sequence took place at ~40 – 36 Ma, and that the sedimentary protolith was deposited after ~53 – 47 Ma. These metamorphic ages are consistent with the reported ages of ~37–36 Ma plutonic rocks in the upper sequence, but contrast with the ~21–19 Ma ages of metamorphic and plutonic rocks in the lower sequence. Therefore, we conclude that the upper and lower metamorphic sequences developed independently but coupled with each other before ~19 Ma as a result of dextral reverse tectonic movement.     

Modelling the biological half-life and seasonality of ¹⁴C in Fucus vesiculosus from the east coast of Ireland: implications for the estimation of future trends

Radiocarbon levels were recorded in Fucus vesiculosus samples collected on a monthly basis over a three-year period at a site on the east coast of Ireland. The resulting data was analysed using a numerical model which estimates the transit times from the Sellafield plant to the sampling location, and the mean availability time of ¹⁴C in seaweed. With the inclusion of a model parameter allowing for seasonal variability in uptake by the Fucus, good correlation was observed between the predicted and measured concentrations. Future temporal trends of ¹⁴C Fucus concentrations along the eastern Irish coastline were modelled with the application of three possible prospective discharge scenarios, predicting ¹⁴C Fucus concentrations to reduce to ambient background levels within 2.5-years of discharges being set to zero. Such projections may prove helpful in assessing the consequences of discharge management and policy making in the context of the OSPAR convention.     

Rare earth element and yttrium compositions of the Paleoproterozoic Yuanjiacun BIF in the Lüliang area and their implications for the Great Oxidation Event(GOE)

In China, most Precambrian banded iron formations(BIFs) are situated in the North China Craton. The Yuanjiacun iron deposit, located in the Lüliang area, is arguably the most representative Superior-type BIF. This iron deposit is coherent with the sedimentary rock succession of the Yuanjiacun Formation in the lower Lüliang Group, and was interpreted to be deposited at 2.3–2.1 Ga, based on ages of overlying and underlying volcanic strata. This age overlaps with the time range of the Great Oxidation Event(GOE, 2.4–2.2 Ga). The Yuanjiacun BIF consists mainly of subhedral-xenomorphic magnetite and quartz and rarely other minerals with a lower degree of metamorphism, from greenschist to lower amphibolite facies. The geochemical characteristics of this BIF are similar to those of Superior-type BIFs. Prominent positive La, Y, and Eu anomalies normalized by the Post Archean Australian Shale(PAAS) indicate that the primary chemical precipitate is a result of solutions that represent mixtures of seawater and high-T hydrothermal fluids. The contamination from crustal detritus found is negligible based on low abundances of Al2O3 and TiO2(0.5%) and of trace elements such as Th, Hf, Zr, and Sc(1.5 ppm), as well as the lack of co-variations between Al2O3 and TiO2. In particular, the Yuanjiacun BIF samples do not display significant negative Ce anomalies like those of the Archean iron formations, but rather, the Yuanjiacun BIF samples exhibit prominent positive Ce anomalies, low Y/Ho ratios, and high light to heavy REE((Pr/Yb)SN) ratios, which are essentially consistent with the late Paleoproterozoic(2.0 Ga) BIFs around the world. These characteristics of the Yuanjiacun BIF samples imply that the ancient ocean(2.3–2.1 Ga) was redox-stratified from oxic shallow water to deeper anoxic water. The specific redox conditions of the ancient ocean may be related to the GOE, which gave rise to the oxidation of Ce and Mn in the upper water, and to the presence of a Mn oxide shuttle in the ocean, resulting in varying REE patterns due to the precipitation and dissolution of this Mn oxide shuttle under different redox states. Therefore, the Yuanjiacun BIF appears to have formed near the redoxcline and lower-level reduced marine water.     

Th–U–total Pb geochronology of authigenic monazite in the Adelaide rift complex,South Australia,and implications for the age of the type Sturtian and Marinoan glacial deposits

The Adelaide rift complex in South Australia contains the type sections for Sturtian and Marinoan glacial deposits. The litho- and chemo-stratigraphy of these deposits play a central role in evaluating global Neoproterozoic ice age hypotheses and Rodinia supercontinent reconstructions, but reliable depositional age constraints have been extremely limited. We report results of in situ Th–U–total Pb (electron microprobe) dating of detrital and authigenic monazite in two samples from the Umberatana Group (Sturtian Holowilena Ironstone and pre-Marinoan Enorama Shale) in the Central Flinders Ranges. Several texturally and chemically distinct detrital and authigenic populations are recognized. Detrital dates range from 1600 Ma to 760 Ma and most relate to well-known orogenic or igneous events in surrounding cratonic regions. Authigenic monazite grew in three or more pulses ranging from 680 Ma to 500 Ma. The date of 680 ± 23 Ma (2σ) for the earliest generation of authigenic monazite in sandstone from the Enorama Shale (1) provides an estimate for the age of the base of the Trezona carbon isotopic anomaly just beneath the Marinoan glacial deposits, (2) provides an absolute minimum age constraint on the underlying Sturtian glacial deposits, and (3) supports proposed correlations between type Marinoan deposits and precisely dated glacial deposits in Namibia and China, which bracket the presumed Marinoan equivalents between 655 and 635 Ma. This age is inconsistent with a Re–Os isochron age of 643 ± 2.4 Ma (2σ) on shales near the bottom of the Sturtian–Marinoan interglacial succession, stratigraphically > 3000 m below the Enorama Shale sample, and militate against the hypothesis that the type Marinoan is correlative with the 580 Ma Gaskiers glaciation. Monazite growth near 600 Ma and again at about 500 Ma probably represent hydrothermal fluid-flow events, the latter of which also corresponds to the well-known Delamerian Orogeny during which the Adelaide sediments were folded into their present structural pattern.     

The geochronology and origin of mantle sources for late cenozoic intraplate volcanism in the frontal part of the Arabian plate in the Karacadağ neovolcanic area of Turkey. Part 1. The results of isotope-geochronological studies

This paper considers results from isotope-geochronological (K-Ar) studies of the products of Neogene-Quaternary volcanism in the Karacada? area, which is situated within the northern frontal part of the Arabian plate. It was found that magmatic activity has been evolving at this location for at least the last 11–10 Myr and was distinctly discrete in character. Three stages of volcanism have been identified: (I) Early or Miocene, ～11–6.7 Ma; (II) Middle or Pliocene-Early Quaternary, 4–1 Ma; and (III) Late or Late Quaternary, 0.4–0.1 Ma. The most recent manifestations of magmatic activity in the region date back to about 100000 years ago. An analysis of the spatial distribution of volcanic centers of different ages in the Karacada? neovolcanic area shows that the magmatism of that region involved a lateral migration of activity from northwest to southeast along a major regional tectonic fault. The migration was caused by the movement of local tension zones where the lithosphere was thinner and deep-seated mantle magmas were ascending.     

Quartz OSL dating of archaeological sites in Xiao Qaidam Lake of the NE Qinghai–Tibetan Plateau and its implications for palaeoenvironmental changes

In the past decades, archaeologists have found evidences for prehistorical human activity in the Qinghai–Tibetan Plateau (QTP). In 1982, some Paleolithic stone tools were found in a section from a terrace of the Xiao Qaidam Lake in the Qaidam Basin, NE of the QTP. The age of this Paleolithic site has remained unknown by far. Some believed that the age of human inhabitation in this Paleolithic site was about 30 ka. In this study, quartz optically stimulated luminescence was used to date 10 samples collected from four sections in the Xiao Qaidam Lake, using the single-aliquot regeneration-dose protocol. The two samples from section XCDH2, which is from a lake terrace about ～7–8 m above the present lake level and in which the top gravel layer contains stone tools, were not well-bleached before deposition. Their ages (>101 and >159 ka) determined by SAR should be considered minimums. OSL dating results of six samples from two sections (XCDH1 and XCDH3) of an adjacent lake terrace, which is ～12 m above the present lake level, suggest two possibilities for the age of the tool-bearing gravel layer: (1) younger than ～3 ka if the lake terrace of XCDH2 is younger than the terrace represented by XCDH1 and XCDH3; or (2) between ～3 and 11 ka if these two terraces are part of the deposit of the same time period. In either case, the age of the archaeological layer should be much younger than the previously proposed ～30 ka. As the climate in the early Holocene after 11 ka was increasingly warm and the Xiao Qaidam Lake area could be suitable for human inhabitation then, we deduce that the age range of ～3–11 ka is more likely the time frame for this archaeological site. The age of 3.1 ± 0.3 ka for the surface of terrace XCDH1/XCDH3 suggests a significant lake level decrease after this time and a corresponding arid event at ～3 ka; the lake level did not reach this level again after that time. Section XCDH4 is more than 40 m above the present lake level, and two samples gave ages of 37 ± 4 and 51 ± 4 ka. These two dates and the dates from the other sections demonstrate that two lake levels higher than present existed for Xiao Qaidam Lake, one at ～12 m and dated ～3–11 ka and the another at >40 m and dated ～37–51 ka.     

The geochronology and origin of mantle sources for late cenozoic intraplate volcanism in the frontal part of the Arabian plate in the Karacadağ neovolcanic area of Turkey. Part 2. The results of geochemical and isotope (Sr-Nd-Pb) studies

A geochemical and isotope-geochemical (Sr-Nd-Pb) study has been carried out for the Karacada? neovolcanic area, which is situated within the frontal part of the Arabian plate. The obtained data and the results of petrological modeling show that the petrogenesis of parental magmas in the Karacada? neovolcanic area involved two compositionally different mantle sources; one consisted of garnet-bearing peridotites of the asthenosphere mantle and the other was spinel-bearing peridotites of the enriched subcontinental lithosphere mantle. During early stages in the evolution of the magmatic system, deep-seated asthenospheric magmas were ascending to the surface while intensively interacting with the melts that had been generated at upper mantle depths. The interaction gradually diminished, so that the later effusive rocks mostly have compositions that are similar to those of the primitive asthenospheric magmas. It is shown that a significant (up to 17–18 wt % of the mantle melt) assimilation of crustal material could take place only during the initial phases of the magmatism. Periodic replenishment of the magma chambers by primitive magmas, which resulted in an observable high degree of homogeneity in the composition of young effusive rocks, was also of importance in the petrogenesis of lavas during the evolution of volcanic activity.