An outline of the geology and geochemistry,and the possible petrogenetic evolution of the volcanic rocks of the Tonga-Kermadec-New Zealand island arc

The Pleistocene-Recent volcanism of this arc extends nearly linearly NNE from northern New Zealand for some 2800 km. Along its western margin lies an active marginal basin (Lau Basin and Havre Trough) which has its southern termination in the Taupo volcanic zone (TVZ, New Zealand). The New Zealand arc segment is developed within a continental crust, whereas the Tonga-Kermadec segments are developed on a ridge system within the oceanic basin. Submarine morphology suggests that the Kermadec volcanoes represent a less advanced stage of evolution relative to those of Tonga.Magmas erupted within the TVZ are dominantly rhyolitic (≈16,000 km³) with subordinate andesites and rare high-alumina tholeiites and dacites. The Kermadec Islands are dominated by tholeiites and basaltic andesites, with subordinate andesites and dacites. The Tongan Islands are dominated by basaltic andesites, with locally developed andesites and dacites. These Tonga-Kermadec lavas are characterised by subcalcic groundmass clinopyroxenes, whereas the younger group of TVZ andesites contain groundmass hypersthene and augite.Geochemically, the TVZ andesites are systematically enriched (relative to those of Tonga-Kermadec) in “incompatible” elements (e.g. K, Rb, Cs, Ba, light REE, U, Th, Zr, Pb), are less Fe-enriched, and contain more radiogenic Sr and Pb (excepting certain ²⁰⁷Pb/²⁰⁴Pb compositions). The evidence points to crustal equilibration of the TVZ andesites prior to eruption.A complete overlap of major and trace element chemistry (including TiO₂) is observed between the Kermadec-TVZ tholeiites and basaltic andesites, and the ocean floor tholeiites of the Lau Basin. Compared to the Tongan lavas, those of the Kermadecs exhibit a greater degree of chemical variability, also reflected in the greater heterogeneity in their Pb isotopic compositions. Moreover, many of the Tonga-Kermadec basaltic andesites exhibit more depleted “incompatible” trace element abundances than the Kermadec and TVZ tholeiites.The “primary” magmas of this arc are interpreted to be of basaltic andesite type, derived from Benioff zone melting (essentially anhydrous), but extensively modified by low-pressure crystal fractionation processes. The Kermadec tholeiites are explained as products of relatively shallow upper mantle partial fusion induced during the earlier stages of diapiric rise of Benioff zone-derived magmas, which are sufficiently hot to intersect the peridotite solidus. This should result in the production and intermixing of a series of magmas extending from olivine tholeiite to basaltic andesite composition. The voluminous rhyolites of TVZ are interpreted as the products of crustal fusion involving Mesozoic sediments.     

A magnetic model for deep plutonic bodies beneath the central Taupo Volcanic Zone, North Island, New Zealand

The central Taupo Volcanic Zone (TVZ) is a region of intense Quaternary rhyolitic volcanism and geothermal activity in the North Island of New Zealand from which about 14,000 km³ of pyroclastics and lavas have been erupted during the last 1.6 Ma. Analysis of aeromagnetic surveys over the TVZ showed the presence of long-wavelength (10 to 25 km) magnetic anomalies which roughly follow the trend of the currently active eastern TVZ, from the north of Lake Taupo to the east of Lake Rotorua. An interpretation of the long-wavelength magnetic anomalies using 3-D magnetic modelling suggests that these anomalies are caused by the magnetic effects of < 3 km thick sequence of volcanic rocks and deeper magnetised bodies within the non-magnetic upper crust (4–7 km depth) beneath the young (age < 0.7 Ma), currently active eastern TVZ. The deep magnetised bodies are interpreted as solidified rhyolitic sub-volcanic plutons that have cooled down to below their Curie temperature.Although the existence of plutonic bodies beneath the TVZ has been postulated prior to this study, this magnetic interpretation result appears to be the first geophysical model of such bodies.     

Geophysical evidence on the structure of the Taupo Volcanic Zone and its hydrothermal circulation

The Taupo Volcanic Zone (TVZ) of New Zealand is characterised by extensive volcanism and by high rates of magma production. Associated with this volcanism are numerous high-temperature (> 250 °C) geothermal systems through which the natural heat output of 4200 ± 500 MW is channelled. Outside the geothermal fields the heat flow is negligible. The average heat flux from the central 6000 km² of the TVZ, which contains most of the geothermal fields, is 700 mW/m³. This heat flux appears to be more concentrated along the eastern margin of the TVZ.Schlumberger resistivity measurements (AB/2 of 500 m and 1000 m) have identified 17 distinct geothermal fields with natural heat outputs greater than 20 MW. An additional six, low-heat-output geothermal fields also occur, and may represent formerly more active systems now in decline. Two extinct fields have also been identified. The average spacing between fields is 10–15 km. The distribution of geothermal fields does not appear to be directly associated with individual volcanic features except for the geothermal system that occurs within Lake Taupo and which occupies the vent of the 1800 yr.B.P. Taupo eruption. The positions of the geothermal fields do not appear to have varied for at least the last 200,000 years. These data are consistent with a model of large-scale convection occurring throughout the TVZ, in which the geothermal fields represent the upper portion of the rising, high-temperature, convective plumes. The majority of the recharge to the convection system is provided by the downward movement of cold meteoric water between the fields which suppresses the heat flow in these regions.Gravity measurements indicate that to a depth of about 2.5 km the upper layers of the TVZ consist of low-density pyroclastic infill. A seismic refraction interface with velocity change from 3.2 km/s to 5.5 km/s occurs at a similar depth. The cross-sectional area of the convection plumes (identified electrically) appears to increase at depths of 1–2 km, consistent with a decrease in permeability at the depth at which the velocity and density increase.The seismicity is dominated by swarm activity which accounts for about half of all earthquakes and is highly variable in both space and time. The small number of seismic events (and swarms) that have well determined depths show a cut off of seismicity at depths of 7–9 km. The depth of the transition from brittle to ductile behaviour of the rocks is identified with the transition from a regime where heat is transported by (hydrothermal) convection and pore pressures are near-hydrostatic to a regime where heat transport is dominantly conductive and pore pressures are lithostatic. Within the convective region, temperatures are moderated by the circulation of water so that the depth of the transition from convective to conductive heat transfer can be linked to the bottom of the seismogenic zone. Rocks must become ductile within about 1 km of the bottom of the overlying convective zone.Seismic refraction studies suggest that the crust beneath the TVZ is highly thinned with a seismic velocity of about 7.5 km/ s, typical of the upper mantle, occurring at depth of 15 km. Seismological studies indicate the upper mantle is highly attenuating beneath the TVZ. Conductive heat transfer between the bottom of the convective system, at about 8 km, and the base of the material with crustal velocities, at 15 km, is not able to provide all the heat that is discharged at the surface. Repeated intrusion from the mantle may provide the additional heat transport required.     

The petrology, phase relations and tectonic setting of basalts from the taupo volcanic zone, New Zealand and the Kermadec Island arc - havre trough, SW Pacific

Volcanism in the Taupo Volcanic Zone (TVZ) and the Kermadec arc-Havre Trough (KAHT) is related to westward subduction of the Pacific Plate beneath the Indo-Australian Plate. The tectonic setting of the TVZ is continental whereas in KAHT it is oceanic and in these two settings the relative volumes of basalt differ markedly. In TVZ, basalts form a minor proportion (< 1%) of a dominant rhyolite (97%)-andesite association while in KAHT, basalts and basaltic andesites are the major rock types. Neither the convergence rate between the Pacific and Indo-Australian Plates nor the extension rates in the back-arc region or the dip of the Pacific Plate Wadati-Benioff zone differ appreciably between the oceanic and continental segments. The distance between the volcanic front and the axis of the back-arc basin decreases from the Kermadec arc to TVZ and the distance between trench and volcanic front increases from around 200 km in the Kermadec arc to 280 km in TVZ. These factors may prove significant in determining the extent to which arc and backarc volcanism in subduction settings are coupled.All basalts from the Kermadec arc are porphyritic (up to 60% phenocrysts) with assemblages generally dominated by plagioclase but with olivine, clinopyroxene and orthopyroxene. A single dredge sample from the Havre Trough back arc contains olivine and plagioclase microphenocrysts in glassy pillow rind and is mildly alkaline (< 1% normative nepheline) contrasting with the tholeiitic nature of the other basalts. Basalts from the TVZ contain phenocryst assemblages of olivine + plagioclase ± clinopyroxene; orthopyroxene phenocrysts occur only in the most evolved basalts and basaltic andesites from both TVZ and the Kermadec Arc.Sparsely porphyritic primitive compositions (Mg/(Mg+Fe²) > 70) are high in Al₂O₃ (>16.5%), and project in the olivine volume of the basalt tetrahedron. They contain olivine (Fo₈₇) phenocrysts and plagioclase (> An₆₀) microphenocrysts. These magmas have ratios of CaO/Al₂O₃, A1₂O₃/TiO₂ and CaO/TiO₂ in the range of MORB and MORB picrites and can evolve to the low-pressure MORB cotectic by crystallisation of olivine±plagiociase. Such rocks may be the parents of other magmas whose evolutionary pathways are complicated by interaction of crystal fractionation, crystal accumulation and mixing processes and the filtering action of crust of variable density and thickness. The interplay of these processes likely accounts for the scatter of data about the cotectic. More evolved rocks from both TVZ and KAHT contain clinopyroxene and orthopyroxene phenocrysts and their compositions merge with basaltic andesites and andesites. Stepwise least-squares modelling using phenocryst assemblages in proportions observed in the rocks suggest that crystal fractionation and accumulation processes can account for much of the diversity observed in the major-element compositions of all lavas.We conclude that the parental basaltic magmas for volcanism in the TVZ and KAHT segments are similar thereby implying grossly similar source mineralogy. We attribute the diversity to secondary processes influencing liquids as they ascended through complex plumbing systems in the sub arc mantle and cross.     

Gravity profiles across the Cheyenne Belt, a precambrian crustal suture in southeastern Wyoming

Geologic discontinuities across the Cheyenne Belt of southeastern Wyoming have led to interpretations that this boundary is a major crustal suture separating the Archaean Wyoming Province to the north from accreted Proterozoic island arc terrains to the south. Gravity profiles across the Cheyenne Belt in three Precambrian-cored Laramide uplifts show a north to south decrease in gravity values of 50–100 mgal. These data indicate that the Proterozoic crust is more felsic (less dense) and/or thicker than Archaean crust. Seismic refraction data show thicker crust (48–54 km) in Colorado than in Wyoming (37–41 km). We model the gravity profiles in two ways: 1) thicker crust to the south and a south-dipping ramp in the Moho beneath and just south of the Cheyenne Belt; 2) thicker crust to the south combined with a mid-crustal density decrease of about 0.05 g/cm³. Differences in crustal thickness may have originated 1700 Ma ago because: 1) the gravity gradient is spatially related to the Cheyenne Belt which has been immobile since about 1650 Ma ago; 2) the N-S gradient is perpendicular to the trend of gravity gradients associated with local Laramide uplifs and sub-perpendicular to regional long-wavelength Laramide gradients and is therefore probably not a Laramide feature. Thus, gravity data support the interpretation that the Cheyenne Belt is a Proterozoic suture zone separating terrains of different crustal structure. The gravity “signature” of the Cheyenne Belt is different from “S”-shaped gravity anomalies associated with Proterozoic sutures of the Canadian Shield which suggests fundamental differences between continent-continent and island arc-continent collisional processes.     

Genesis of lavas of the Taupo Volcanic Zone, North Island, New Zealand

The Taupo Volcanic Zone forms part of the Taupo-Hikurangi subduction system, and comprises five volcanic centres: Tongariro, Taupo, Maroa, Okataina and Rotorua. Tongariro Volcanic Centre is formed almost entirely of andesite while the other four centres contain predominantly rhyolitic volcanics and later fissure eruptions of high-Al basalt. Estimated total volume of each lava type are as follows: 2 km³ of high-Al basalt (< 0.1%); 260 km³ of andesite (< 2.5%); 5 km³ of dacite (< 0.1%); > 10,000 km³ of rhyolite and ignimbrite (> 97.4%).The location of the andesites and vent alignments suggest a source from a subduction zone underlying the area. However, the lavas differ chemically from island-arc andesites such as those of Tonga; in particular by having higher contents of the alkali elements, light REE and Sr and Pb isotopes. This suggests some crustal contamination, and it is considered that this may occur beneath the wide accretionary prism of the subduction system. Amphibolite of the subduction zone will break down between 80 and 100 km and a partial melt will rise. A multi-stage process of magma genesis is then likely to occur. High-Al basalts are thought to be derived from partial melting of a garnet-free peridotite near the top of the mantle wedge overlying the subduction zone, locations of the vents controlled largely by faults within the crust. Rhyolites and ignimbrites were probably derived from partial melting of Mesozoic greywacke and argillite under the Taupo Volcanic Zone. Initial partial melting may have been due to hydration of the base of the crust; the “water” having come from dehydration of the downgoing slab. The partial melts would rise to form granodiorite plutons and final release of the magma to form rhyolites and ignimbrites was allowed because of extension within the Taupo graben.Dacites of the Bay of Plenty probably resulted from mixing of andesitic magma with small amounts of rhyolitic magma, but those on the eastern side of the Rotorua-Taupo area were more likely formed by a higher degree of partial melting of the Mesozoic greywacke-argillite basement. This may be due to intrusion of andesite magma on this side of the Taupo volcanic zone.     

Tectonic controls on magma genesis and evolution in the northwestern United States

Field, chronologic, chemical, and isotopic data for late Cenozoic basaltic rocks from the northwestern United States illustrate the relationship between crustal structure and tectonic forces in controlling the genesis and evolution of continental volcanism. In the northwestern U.S., the first major episode of basaltic volcanism was triggered by crustal rifting in a “back-arc” environment, east of the westward-migrating volcanic arc created by the subduction of the Juan-de-Fuca plate beneath the North American plate. Rifting and volcanism were concentrated by pre-existing zones of crustal weakness associated with boundaries between the old Archean core of the continent and newly accreted terranes. Basalts erupted during this time (Columbia River, Steens Mountain) show evidence of significant fractionation histories including contamination by crust of varying age depending on the crustal structure at the eruption site. Presumably this reflects ponding and stagnation of primary magmas in the crust or at the crust-mantle interface due to their encounter with thick crust, not yet extended and still containing its low-density, easily fusible component. Continued rifting of this crust, and modification of its composition through extraction of rhyolitic partial melts and deposition of the fractionation products from primary basaltic melts, coupled with a shift in stress orientation roughly 10.5 Ma ago, allowed relatively unfractionated and uncontaminated magmas to begin reaching the surface. In the western part of the region (Oregon Plateau), these magmas tapped a mantle source similar to that which produced most of the ocean island basalts of the northern hemisphere. To the east (Snake River Plain), however, the mantle sampled by basaltic volcanism has isotopic characteristics suggesting it has preserved a record of incompatible element enrichment processes associated with the formation of the overlying Archean crustal section some 2.6 Ga ago.     

A Jurassic-Cretaceous island arc in the Ladakh-Himalayas

The Ladakh Mesozoic ophiolite belt (western Himalaya) contains a pile of volcanic thrust sheets (Dras unit) which differ significantly in structure and composition from the ophiolitic mélange zones. The Dras unit is composed of pillow lavas, doleritic sills, very irregular basaltic (?basaltic andesites) and dacitic flows intercalated with pyroclastics, volcanoclastic sediments and radiolarian cherts. According to fossil evidence, this volcanism must have been active between Upper Jurassic and Upper Cretaceous.The presence of relict primary minerals, such as magnesiochromite, clinopyroxene, hastingsitic hornblende and Ti-magnetite as well as distinctive bulk chemistries, suggests that the volcanics belong to island arc tholeiite and to calc-alkaline rock series, typical of present island arcs in the Caribbean and Pacific.Model calculations incorporating probed phenocryst phases indicate that in addition to olivine, clinopyroxene and plagioclase, amphibole and titanomagnetite are crucial fractionating phases in the development of the dacites from a primitive tholeiitic melt. The latter process must have taken place at about 1000°C and at moderate depth of 5–15 km within or underneath the island arc. Today, hornblende-bearing mafic cumulates appear in the vicinity of Kargil within and close to the Dras volcanics.In a Sr-evolution diagram, the Dras volcanics have yielded a “pseudo-isochron” with a low initial ratio of 0.7035 ± 0.0003, which is in the same range as the mean of modern island arc volcanics. However, a geologically unrealistic age of 263 m.y., is obtained from the slope of this isochron.The upper mantle is regarded as the source material for the island arc tholeiitic magmas. Enrichment in K, Ba, Sr and LREE supports the involvement of components derived from dehydration or incipient melting of subducted Tethyan oceanic crust in the mantle.     

Trace-element and isotopic constraints on the source of magmas in the active volcano and Mariana island arcs, Western Pacific

Analytical results of the relative and absolute abundance of LIL-incompatible trace elements (K, Rb, Cs, Sr, and Ba) and isotopic compositions ( , , and ) are summarized for fresh samples from active and dormant volcanoes of the Volcano and Mariana island arcs. The presence of thickened oceanic crust (T 15–20 km) beneath the arc indicates that while hybridization processes resulting in the modification of primitive magmas by anatectic mixing at shallow crustal levels cannot be neglected, the extent and effects of these processes on this arc's magmas are minimized. All components of the subducted plate disappear at the trench. This observation is used to reconstruct the composition of the crust in the Wadati-Benioff zone by estimating proportions of various lithologies in the crust of the subducted plate coupled with analyses from DSDP sites. Over 90% of the mass of the subducted crust consists of basaltic Layers II and III. Sediments and seamounts, containing the bulk of the incompatible elements, make up the rest. Bulk Western Pacific seafloor has , δ ¹⁸O +7.2, K/Rb 510, K/Ba 46, and K/Cs 13,500. Consideration of trace-element data and combined systematics limits the participation of sediments in magmagenesis to less than 1%, in accord with the earlier results of Pb-isotopic studies. Combined data indicate little, if any, involvement of altered basaltic seafloor in magmagenesis. Perhaps more important than mean isotopic and LIL-element ratios is the restricted range for lavas from along over 1000 km of this arc. Mixtures of mantle with either the subducted crust or derivative fluids should result in strong heterogeneities in the sources of individual volcanoes along the arc. Such heterogeneities would be due to: (1) gross variations of crustal materials supplied to the subduction zone; and (2) lesser efficiency of mixing processes accompanying induced convection between arc segments (parallel to the arc) as compared to that perpendicular to the arc. The absence of these heterogeneities indicates that either some process exists for the efficient mixing of mantle and subducted material parallel to the arc or that subducted materials play a negligible role in the generation of Mariana-Volcano arc melts.Consideration of plausible sources in the mantle indicates that (1) an unmodified MORB-like mantle cannot have generated the observed trace-element and isotopic composition of this arc's magmas, while (2) a mantle similar to that which has produced alkali-olivine basalts (AOB) of north Pacific “hot spot” chains is indistinguishable in many respects spects from the source of these arc lavas.     

Mean crustal velocity: a critical parameter for interpreting crustal structure and crustal growth

Mean crustal velocity is a critical parameter for genesis of continental crystalline crust because it is a function of mean crustal composition and therefore may be used to resolve continental crustal growth in space and time. Although the best values of mean crustal velocity are determined from wide-angle reflection measurements, most studied here necessarily come from vertical averages in crustal refraction determinations. The mode of 158 values of mean crustal velocity is 6.3 km/s, a velocity which corresponds to a mean crustal composition of granodiorite to felsic quartz diorite; Archean crust may be slightly more mafic. Mean crustal velocities range from 5.8 to 7.0 km/s. The lowest values invariably are found in thermally disturbed rift zones and the highest values correspond to velocities in gabbro. Velocities in island arcs may be as low as 6.0 km/s but are typically 6.5–6.9 km/s which corresponds to andesitic composition; estimates of island arc composition are andesitic. If values of mean crustal velocity are not biased, this observation suggests that continental crust did not grow simply by addition of island arc material. Possibilities are that crust formed from fusion of island arcs and was later changed to more felsic composition by addition of material from the mantle or that the late Archean episode of major crustal growth did not involve processes similar to younger island arcs. Some crustal blocks might be changed in composition and thickness by such processes as underplating, interthrusting, necking and sub-crustal erosion. Specially designed experiments are suggested to determine this parameter so critical for understanding genesis of continental crust.     

Andesites in island arcs and continental margins: Relationship to crustal evolution

Andesites of both island arc and continental margin environments contain petrologic evidence of mixing of mantle and crustal melts. Andesitic volcanism appears to involve addition of mantle-derived basaltic magma to the crust and fractionation of preexisting crustal material. Changes in andesitic volcanism with increasingly continental character of the crust reflect changes in a rhyolitic component derived from increasingly aged and fractionated crust. The initial stage in development of continental crust is partial melting of oceanic crust.     

Thermal state, rheology and seismicity in the pannonian basin, Hungary

On the basis of data on crustal structure and terrestrial heat flow, a 3-D geothermal model for the lithosphere in the Pannonian basin, Hungary, has been calculated. This model, together with information on crustal composition, laboratory data on rock friction, and certain assumptions about fluid conditions and strain-rate levels within the lithosphere, has been used to construct a rheological model of the area.The results obtained show a layered rheological structure where an aseismic part of the crust is “sandwiched” between an upper and a lower seismogenic crustal layers. According to the proposed rheological model, seismic activity in the upper crust may be expected down to depths of 10–12 km, which is confirmed well by the observed depth distribution of seismicity. The model also predicts a lower crustal seismogenic layer down to 20–22 km. Because of infrequent occurrences of deep earthquakes and/or a generally small number of reliable hypocenter depth determinations in the study area, this seismogenic zone is less constrained by observations.The depth of the different rheologic horizons within the crust is governed mainly by thermal conditions. The lower boundary of both seismogenic layers appears isothermal. Brittle-ductile transition in the upper crust coincides with the ˜200 °C isotherm, while in the lower crust it coincides with the ˜ 375 °C isotherm. The lowermost crust and the upper mantle beneath Hungary show ductile behavior, thus the possibility of siesmic activity at these horizons can be excluded.     

The Zealandia Volcanic Complex: Further evidence of a lower crustal “hot zone” beneath the Mariana Intra‐oceanic Arc,Western Pacific

This paper addresses formation of felsic magmas in an intra‐oceanic magmatic arc. New bathymetric, petrologic, geochemical, and isotopic data for Zealandia Bank and two related volcanoes in the south‐central Mariana arc is presented and interpreted. These three volcanoes are remnants of an older andesitic volcano that evolved for some time and became dormant long enough for a carbonate platform to grow on its summit before reawakening as a rhyodacitic volcano. Zealandia lavas are transitional between low‐ and medium‐K and tholeiitic and calc‐alkaline suites. They define a bimodal suite with a gap of 56–58 wt% SiO₂; this suggests that mafic and felsic magmas have different origins. The magmatic system is powered by mantle‐derived basalts having low Zr/Y and flat rare earth element patterns. Two‐pyroxene thermometry yields equilibration temperatures of 1000–1100 °C for andesites and 900–1000 °C for dacites. Porphyritic basalts and andesites show textures expected for fractionating magmas but mostly fine‐grained felsic lavas do not. All lavas show trace element signatures expected for mantle and crustal sources that were strongly melt‐depleted and enriched by subduction‐related fluids and sediment melts. Sr and Nd isotopic compositions fall in the normal range of Mariana arc lavas. Felsic lavas show petrographic evidence of mixing with mafic magma. Zealandia Bank felsic magmatism supports the idea that a large mid‐ to lower‐crustal felsic magma body exists beneath the south‐central Mariana arc, indicating that MASH (mixing, assimilation, storage, and homogenization) zones can form beneath intra‐oceanic as well as continental arcs.     

Nd-isotope systematics of ∼2.7 Ga adakites, magnesian andesites, and arc basalts, Superior Province: evidence for shallow crustal recycling at Archean subduction zones

An association of adakite, magnesian andesite (MA), and Nb-enriched basalt (NEB) volcanic flows, which erupted within ‘normal’ intra-oceanic arc tholeiitic to calc-alkaline basalts, has recently been documented in ∼2.7 Ga Wawa greenstone belts. Large, positive initial ?_Nd values (+1.95 to +2.45) of the adakites signify that their basaltic precursors, with a short crustal residence, were derived from a long-term depleted mantle source. It is likely that the adakites represent the melts of subducted late Archean oceanic crust. Initial ?_Nd values in the MA (+0.14 to +1.68), Nb-enriched basalts and andesites (NEBA) (+1.11 to +2.05), and ‘normal’ intra-oceanic arc tholeiitic to calc-alkaline basalts and andesites (+1.44 to +2.44) overlap with, but extend to lower values than, the adakites. Large, tightly clustered ?_Nd values of the adakites, together with Th/Ce and Ce/Yb systematics of the arc basalts that rule out sediment melting, place the enriched source in the sub-arc mantle. Accordingly, isotopic data for the MA, NEBA, and ‘normal’ arc basalts can be explained by melting of an isotopically heterogeneous sub-arc mantle that had been variably enriched by recycling of continental material into the shallow mantle in late Archean subduction zones up to 200 Ma prior to the 2.7 Ga arc. If the late Archean Wawa adakites, MA, and basalts were generated by similar geodynamic processes as their counterparts in Cenozoic arcs, involving subduction of young and/or hot ocean lithosphere, then it is likely that late Archean oceanic crust, and arc crust, were also created and destroyed by modern plate tectonic-like geodynamic processes. This study suggests that crustal recycling through subduction zone processes played an important role for the generation of heterogeneity in the Archean upper mantle. In addition, the results of this study indicate that the Nd-isotope compositions of Archean arc- and plume-derived volcanic rocks are not very distinct, whereas Phanerozoic plumes and intra-oceanic arcs tend to have different Nd-isotopic compositions.     

Adakitic lavas in the Central Luzon back‐arc region,Philippines: lower crust partial melting products?

Abstract Volcanism in the back-arc side region of Central Luzon, Philippines, with respect to the Manila Trench is characterized by fewer and smaller volume volcanic centers compared to the adjacent forearc side-main volcanic arc igneous rocks. The back-arc side volcanic rocks which include basalts, basaltic andesites, andesites and dacites also contain more hydrous minerals (ie, hornblende and biotite). Adakite-like geochemical characteristics of these back-arc lavas, including elevated Sr, depleted heavy rare earth elements and high Sr/Y ratios, are unlikely to have formed by slab melting, be related to incipient subduction, slab window magmatism or plagioclase accumulation. Field and geochemical evidence show that these adakitic lavas were most probably formed by the partial melting of a garnet-bearing amphibolitic lower crust. Adakitic lavas are not necessarily arc–trench gap region slab melts.     

Sangay volcano, Ecuador: structural development, present activity and petrology

Sangay (5230 m), the southernmost active volcano of the Andean Northern Volcanic Zone (NVZ), sits 130 km above a >32-Ma-old slab, close to a major tear that separates two distinct subducting oceanic crusts. Southwards, Quaternary volcanism is absent along a 1600-km-long segment of the Andes. Three successive edifices of decreasing volume have formed the Sangay volcanic complex during the last 500 ka. Two former cones (Sangay I and II) have been largely destroyed by sector collapses that resulted in large debris avalanches that flowed out upon the Amazon plain. Sangay III, being constructed within the last avalanche amphitheater, has been active at least since 14 ka BP. Only the largest eruptions with unusually high Plinian columns are likely to represent a major hazard for the inhabited areas located 30 to 100 km west of the volcano. However, given the volcano's relief and unbuttressed eastern side, a future collapse must be considered, that would seriously affect an area of present-day colonization in the Amazon plain, 30 km east of the summit. Andesites greatly predominate at Sangay, there being few dacites and basalts. In order to explain the unusual characteristics of the Sangay suite—highest content of incompatible elements (except Y and HREE) of any NVZ suite, low Y and HREE values in the andesites and dacites, and high Nb/La of the only basalt found—a preliminary five-step model is proposed: (1) an enriched mantle (in comparison with an MORB source), or maybe a variably enriched mantle, at the site of the Sangay, prior to Quaternary volcanism; (2) metasomatism of this mantle by important volumes of slab-derived fluids enriched in soluble incompatible elements, due to the subduction of major oceanic fracture zones; (3) partial melting of this metasomatized mantle and generation of primitive basaltic melts with Nb/La values typical of the NVZ, which are parental to the entire Sangay suite but apparently never reach the surface and subordinate production of high Nb/La basaltic melts, maybe by lower degrees of melting at the periphery of the main site of magma formation, that only infrequently reach the surface; (4) AFC processes at the base of a 50-km-thick crust, where parental melts pond and fractionate while assimilating remelts of similar basaltic material previously underplated, producing andesites with low Y and HREE contents, due to garnet stability at this depth; (5) low-pressure fractionation and mixing processes higher in the crust. Both an enriched mantle under Sangay prior to volcanism and an important slab-derived input of fluids enriched in soluble incompatible elements, two parameters certainly related to the unique setting of the volcano at the southern termination of the NVZ, apparently account for the exceptionally high contents of incompatible elements of the Sangay suite. In addition, the low Cr/Ni values of the entire suite—another unique characteristic of the NVZ—also requires unusual fractionation processes involving Cr-spinel and/or clinopyroxene, either in the upper mantle or at the base of the crust.     

An evaluation of the global variations in the major element chemistry of arc basalts

Arc volcanoes occur at convergent margins with a wide range in subduction parameters, and variations in these parameters might be expected to lead to variations in the chemistry of magmas parental to arcs. Major element analyses from approximately 100 volcanic centers within 30 arcs, normalized to 6% MgO to minimize the effects of crystal fractionation, display wide variations. Na₂O and CaO at 6% MgO (Na_6.0 and Ca_6.0) correlate remarkably well with the thickness of the overlying crust. These systematics are consistent with two possible models. In the first model, the crust behaves as a chemical filter; where the crust is thick, magmas crystallize at higher pressure and interact more extensively with the arc crust. Modeling of high-pressure crystallization and assimilation, however, does not reproduce the associated variations in Na_6.0 and Ca_6.0 without calling upon complicated combinations of fractionating phases and assimilants. In the second model, crustal thickness determines the height of the mantle column available for melting beneath arc volcanoes. If melting begins beneath arcs at similar depths, then the column of mantle that undergoes decompression melting is much shorter beneath the thickest arc crust. The shorter mantle column for arcs built on thick crust will lead to smaller extents of melting in the mantle, and hence higher Na_6.0 and lower Ca_6.0 in the parental magmas. Modeling shows that variations in the extent of melting in the mantle can easily account for the associated variations in Ca_6.0 and Na_6.0. The abundances of the other major elements at 6% MgO do not correlate well with crustal thickness, or any other subduction parameter. Co-variation of some of these other major elements (e.g., Si_6.0 and Fe_6.0) within individual arcs suggests that they are strongly influenced by local crustal level processes that obscure partial melting systematics. Correction for the crustal processes improves the relationship between Na_6.0 and Ca_6.0 that is so readily explained by partial melting. The extents of melting in the mantle beneath arc volcanoes estimated from the ranges in Na_6.0 and Ca_6.0 are remarkably similar to those estimated beneath mid-ocean ridges. This observation provides further evidence that the mantle wedge, and not the slab, melts beneath arc volcanic fronts.     

Tectonics and paleomagnetism of structural arcs of the Pamir-Punjab syntaxis

The Pamir-Punjab syntaxis consists of two structural arcs with the crests facing northward—the Pamir and Hindu-Kush-Karakorum arcs. These arcs are mutually disharmonious, and the exterior (Pamir) arc is more tight, as compared to the inner one. Paleomagnetic study of the Pamir arc has shown that the structures of the future Pamirs had a northeastern strike in the Paleogene and Cretaceous, and they occurred on the eastern limb of the Darvaz-Kopetdag structural arc, whose crest faces south. The Pamir arc originated after the Paleogene during the process of the formation of the Pamir-Punjab syntaxis. Knowledge of the kinematics of the Pamir arc, combined with data on the geometry of the syntaxis and the character of its boundaries, enable us to choose a model of the development of the syntaxis. The process of “plastic flow” of crustal masses around the underthrusting segment of the Indian plate was likely of paramount importance.     

Relative depletion of niobium in some arc magmas and the continental crust: partitioning of K, Nb, La and Ce during melt/rock reaction in the upper mantle

Depletion of Nb relative to K and La is characteristic of lavas in subduction-related magmatic arcs, as distinct from mid-ocean ridge basalts. Nb depletion is also characteristic of the continental crust. This and other geochemical similarities between the continental crust and high-Mg# andesite magmas found in arcs suggests that the continental crust may have formed by accretion of andesites. Previous studies have shown that the major element characteristics of high-Mg# andesites may be produced by melt/rock reaction in the upper mantle. In this paper, new data on partitioning of K, Nb, La and Ce between garnet, orthopyroxene and clinopyroxene in mantle xenoliths, and on partitioning of Nb and La between orthopyroxene and liquid, show that garnet and orthopyroxene have Nb crystal/liquid distribution coefficients which are much larger than those of K and La. Similar fractionations of Nb from K and La are expected in spinel and olivine. For this reason, reactions between migrating melt and large masses of mantle peridotite can produce substantial depletion of Nb in derivative liquids. Modeling shows that reaction between ascending, mantle-derived melts and mantle peridotite is a viable mechanism for producing the trace element characteristics of high-Mg# andesite magmas and the continental crust.

Alternatively, small-degree melts of metabasalt and/or metasediment in the subducting slab may leave rutile in their residue, and will thus have large Nb depletions relative to K and La [1]. Slab melts are too rich in light rare earth elements and other incompatible elements, and too poor in compatible elements, to be parental to arc magmas. However, ascending slab melts may be modified by reaction with the mantle. Our new data permit modeling of the trace element effects of reaction between small-degree melts of the slab and mantle peridotite. Modeling shows that this type of reaction is also a viable mechanism for producing the trace element characteristics of high-Mg# andesites and the continental crust. These findings, in combination with previous results, suggest that melt/rock reaction in the upper mantle has been an important process in forming the continental crust and mantle lithosphere.     

Rhyolite volcanism in the Krafla central volcano,north-east Iceland

At the Krafla central volcano in north-east Iceland, two main phases of rhyolite volcanism are identified. The earlier phase (last interglacial) is related to the formation of a caldera, whereas the second phase (last glacial) is related to the emplacement of a ring dike. Subsequently, only minor amounts of rhyolite have been erupted. The volcanic products of Krafla are volumetrically bimodal. Geochemically, there is a series of basaltic to basalto-andesitic rocks and a cluster of rhyolitic rocks. Rocks of intermediate to silicic composition (icelandites and dacites) show clear signs of mixing. The rhyolites are Fe-rich (tholeiitic), and aphyric to slightly porphyritic (plagioclase, augite, pigeonite, fayalitic olivine and magnetite). They are minimum melts on the quartz-plagioclase cotectic plane in the granite system (Qz-Or-Ab-An). The rhyolites at Krafla were produced by near-solidus, rather than nearliquidus fractionation. They are interpreted as silicic minimum melts of hydrothermally altered crust, mainly of basaltic composition. They were primarily generated on the peripheries of an active basaltic magma chamber or intrusive domain, where sufficient volumes of crust were subjected to temperatures favorable for rhyolite genesis (850–950° C). The silicic melts were extracted crystal-free from their source in response to crustal deformation.