The Rock of Gibraltar

Most people have heard of the Rock of Gibraltar, yet few know anything of its geology. Recent studies now interpret the Rock as a partly overturned mass of Jurassic limestone, thrust into position during continent-continent collision by about 15–20 million years ago, and shaped largely by shoreline processes active during Pleistocene tectonic uplift. The Rock's geology has implications for currently developing interpretations of the geological evolution of the western Mediterranean region and for the fast-growing development of Gibraltar's urban environment.     

Constraining the chronology of Pleistocene glaciations on Svalbard: Kapp Ekholm re‐visited

The Kapp Ekholm site, in central Spitsbergen, shows alternating units of glaciomarine sandy silt and diamicton representing three glacial cycles and is key in reconstructing the Late Pleistocene glacial history of Svalbard. Part of the site is reinvestigated here by focusing on re‐dating two units (B and F) interpreted as interglacial/interstadial glaciomarine deposits, in order to constrain the controversial chronology. A combination of Optical Stimulated Luminescence (OSL) on quartz, infrared stimulated luminescence with a 50 °C readout temperature (IRSL₅₀) and post infrared‐infrared stimulated luminescence (pIR), both on feldspar, was applied. While Formation B was beyond the dateable range of OSL, IRSL₅₀ and pIR ages lead to the conclusion that this unit represents the Last Interglacial, Marine Isotope Stage (MIS) 5e, and the underlying diamicton the MIS 6 glacial. Formation F yielded ages implying that the formation represents the MIS 5a interstadial and the underlying diamicton is interpreted to represent the MIS 5b stadial. This agrees with conclusions drawn concerning the Pleistocene glaciations elsewhere on Svalbard.     

Titanite dates crystallization: Slow Pb diffusion during super‐solidus re‐equilibration

Titanite can be found in rocks of wide compositional range, is reactive, growing or regrowing during metamorphic and hydrothermal events, and is generally amenable to U–Pb geochronology. Experimental evidence suggest that titanite has a closure temperature for Pb ranging from 550 to 650°C, and thus titanite dates are commonly interpreted as cooling ages. However, this view has been challenged in recent years by evidence from natural titanite which suggests the closure temperature may be significantly higher (up to 800°C). Here, we investigate titanite in an enclave of migmatitic gneiss included within a granite intrusion. The titanite crystals exhibit textural features characteristic of fluid‐mediated mass transfer processes on length scales of <100 µm. These textural features are associated with variation in both Pb concentrations and distinct U–Pb isotopic compositions. Zr‐in‐titanite thermometry indicates that modification of the titanite occurred at temperatures in excess of 840°C, in the presence of a high‐T silicate melt. The Pb concentration gradients preserved in these titanite crystals are used to determine the diffusivity of Pb in titanite under high‐T conditions. We estimate diffusivities ranging from 2 × 10^?22 to 5 × 10^?25 m²/s. These results are significantly lower than experimental data predict yet are consistent with other empirical data on natural titanites, suggesting that Pb diffusivity is similar to that of Sr. Thus our data challenge the wide‐held assumption that U–Pb titanite dates only reflect cooling ages.     

Svalbard glaciers re‐advanced during the Pleistocene–Holocene transition

Despite warming regional conditions and our general understanding of the deglaciation, a variety of data suggest glaciers re‐advanced on Svalbard during the Lateglacial–early Holocene (LGEH). We present the first well‐dated end moraine formed during the LGEH in De Geerbukta, NE Spitsbergen. This landform was deposited by an outlet glacier re‐advancing into a fjord extending 4.4 km beyond the late Holocene (LH) maximum. Comparing the timing of the De Geerbukta glacier re‐advance to a synthesis of existing data including four palaeoclimate records and 15 other proposed glacier advances from Svalbard does not suggest any clear synchronicity in glacial and climatic events. Furthermore, we introduce six additional locations where glacier moraines have been wave‐washed or cut by postglacial raised marine shorelines, suggesting the landforms were deposited before or during high relative sea‐level stands, thus exhibiting a similar LGEH age. Contrary to current understanding, our new evidence suggests that the LGEH glaciers were more dynamic, exhibited re‐advances and extended well beyond the extensively studied LH glacial expansion. Given the widespread occurrence of the LGEH glacier deposits on Svalbard, we suggest that the culmination of the Neoglacial advances during the Little Ice Age does not mark the maximum extent of most Svalbard glaciers since deglaciation; it is just the most studied and most visible in the geological record.     

Multiple re‐advances of a Lake Vättern outlet glacier during Fennoscandian Ice Sheet retreat,south‐central Sweden

Lake Vättern represents a critical region geographically and dynamically in the deglaciation of the Fennoscandian Ice Sheet. The outlet glacier that occupied the basin and its behaviour during ice‐sheet retreat were key to the development and drainage of the Baltic Ice Lake, dammed just west of the basin, yet its geometry, extent, thickness, margin dynamics, timing and sensitivity to regional retreat forcing are rather poorly known. The submerged sediment archives of Lake Vättern represent a missing component of the regional Swedish deglaciation history. Newly collected geophysical data, including high‐resolution multibeam bathymetry of the lake floor and seismic reflection profiles of southern Lake Vättern, are used here together with a unique 74‐m sediment record recently acquired by drill coring, and with onshore LiDAR‐based geomorphological analysis, to investigate the deglacial environments and dynamics in the basin and its terrestrial environs. Five stratigraphical units comprise a thick subglacial package attributed to the last glacial period (and probably earlier), and an overlying >120‐m deglacial sequence. Three distinct retreat–re‐advance episodes occurred in southern Lake Vättern between the initial deglaciation and the Younger Dryas. In the most recent of these, ice overrode proglacial lake sediments and re‐advanced from north of Visingsö to the southern reaches of the lake, where ice up to 400 m thick encroached on land in a lobate fashion, moulding crag‐and‐tail lineations and depositing till above earlier glacifluvial sediments. This event precedes the Younger Dryas, which our data reveal was probably restricted to north‐central sectors of the basin. These dynamics, and their position within the regional retreat chronology, indicate a highly active ice margin during deglaciation, with retreat rates on average 175 m a^?1. The pronounced topography of the Vättern basin and its deep proglacial‐dammed lake are likely to have encouraged the dynamic behaviour of this major Fennoscandian outlet glacier.     

The Bølling‐age Blomvåg Beds,western Norway: implications for the Older Dryas glacial re‐advance and the age of the deglaciation

Blomvåg, on the western coast of Norway north of Bergen, is a classical site in Norwegian Quaternary science. Foreshore marine sediments, named the Blomvåg Beds and now dated to the Bølling‐Allerød from 14.8 to 13.3 cal. ka BP, contain the richest Lateglacial bone fauna in Norway, numerous mollusc shells, driftwood, and flint that some archaeologists consider as the oldest traces of humans in Norway. The main theme of this paper is that the Blomvåg Beds are overlain by a compact diamicton, named the Ulvøy Diamicton, which was interpreted previously as a basal till deposited during a glacial re‐advance into the ocean during the Older Dryas (c. 14 cal. ka BP). Sediment sections of the Blomvåg Beds and the Ulvøy Diamicton were exposed in ditches in a cemetery that was constructed in 1941–42 and have subsequently not been accessible. A number of radiocarbon and cosmogenic ¹⁰Be exposure ages demonstrate that the diamicton is not likely to be a till because minimum deglaciation ages (14.8–14.5 cal. ka BP) from the vicinity pre‐date the Ulvøy Diamicton. We now consider that sea ice and icebergs formed the Ulvøy Diamicton during the Younger Dryas. The Scandinavian Ice Sheet margin was located on the outermost coastal islands between at least c. 18.5 and 14.8 cal. ka BP; however, no ice‐marginal deposits have been found offshore from this long period. The Older Dryas ice margin in this area was located slightly inside the Younger Dryas margin, whereas farther south it was located slightly beyond the Younger Dryas margin.     

A re‐investigation of the low Reynolds number uniform flow past a porous spherical shell

This work reinvestigates the flow field of a uniform flow past a porous spherical shell based on Song and Huang's (2000) theory of laminar poroelastic media flow with proper boundary conditions. The analytical solution of this study not only indicates viscous effects inside the porous spherical shell, but it also preserves the continuities of tangential velocity and shear stress at the interfaces. The result reveals that as the porosity approaches unity, the flow is entirely comprised of the incident stream; as the porosity approaches zero, the well‐known Stokes' solution of a uniform flow past a sphere with low Reynolds number is obtained. Copyright © 2004 John Wiley & Sons, Ltd.     

Flux‐Free Fusion of Silicate Rock Preceding Acid Digestion for ICP‐MS Bulk Analysis

We present a new method for the decomposition of silicate rocks by flux‐free fusion in preparation for whole‐rock trace element determination (Sc, Rb, Sr, Y, Zr, Nb, Cs, Ba, rare earth elements and Hf) that is especially applicable to zircon‐bearing felsic rocks. The method was verified by analyses of RMs of mafic (JB‐1a, JB‐2, JGb‐1) and felsic rocks (JG‐3, JR‐3, JSd‐1, GSP‐2, G‐2). Pellets of powdered sample (up to 500 mg) without flux were weighed and placed in a clean platinum crucible. The samples were then fused in a Siliconit® tube furnace and quenched to room temperature. The optimum condition for the fusion of granitic rock was determined to be heating for 2–3 min at 1600 °C. The fused glass in the platinum crucible after heating was decomposed using HF and HClO₄ in a Teflon® beaker. Decomposed and diluted sample solutions were analysed using a quadrupole inductively coupled plasma‐mass spectrometer. Replicate analyses (n = 4 or 5) of the RMs revealed that analytical uncertainties were generally < 3% for all elements except Zr and Hf (～ 6%) in JG‐3. These higher uncertainties may be attributed to sample heterogeneity. Our analytical results for the RMs agreed well with recommended concentrations and recently published concentrations, indicating complete decomposition of our rock samples during fusion.     

Precise Determination of Silicon Isotopes in Silicate Rock Reference Materials by MC‐ICP‐MS

A HF‐free sample preparation method was used to purify silicon in twelve geological RMs. Silicon isotope compositions were determined using a Neptune instrument multi‐collector‐ICP‐MS in high‐resolution mode, which allowed separation of the silicon isotope plateaus from their interferences. A 1 μg g^‐1 Mg spike was added to each sample and standard solution for online mass bias drift correction. δ³⁰Si and δ²⁹Si values are expressed in per mil (‰), relative to the NIST SRM 8546 (NBS‐28) international isotopic RM. The total variation of δ³⁰Si in the geological reference samples analysed in this study ranged from ‐0.13‰ to ‐0.29‰. Comparison with δ²⁹Si values shows that these isotopic fractionations were mass dependent. IRMM‐17 yielded a δ³⁰Si value of ‐1.41 ± 0.07‰ (2s, n = 12) in agreement with previous data. The long‐term reproducibility for natural samples obtained on BHVO‐2 yielded δ³⁰Si = ‐0.27 ± 0.08‰ (2s, n = 42) on a 12 month time scale. An in‐house Si reference sample was produced to check for the long‐term reproducibility of a mono‐elemental sample solution; this yielded a comparable uncertainty of ± 0.07‰ (2s, n = 24) over 5 months.     

The initial garnet‐in reaction involving siderite–rhodochrosite,garnet re‐equilibration and P–T–t paths of graphitic schists in the Black Hills orogen,South Dakota,USA

It is generally thought that garnet in metapelites is produced by continuous reactions involving chlorite or chloritoid. Recent publications have suggested that the equilibrium temperatures of garnet‐in reactions may be significantly overstepped in regionally metamorphosed terranes. The growth of small spessartine–almandine garnet crystals on Mn‐siderite at the garnet isograd in graphitic metapelites in the Proterozoic Black Hills orogen, South Dakota, demonstrates that Mn‐siderite was the principal reactant that produced the initial garnet in the schists. Moreover, the positions of garnet compositions in isobaric, T–(C/H) pseudosections for the schists show that the temperature of the garnet‐in reaction from Mn‐siderite was overstepped minimally at the most. In the Black Hills, garnet was initially produced during regional metamorphism beginning at c. 1755 Ma due to the collision of Wyoming and Superior cratons, and was subsequently partially or fully re‐equilibrated at more elevated temperatures and pressures during intrusion of the Harney Peak Granite (HPG) at c. 1715 Ma. Garnet occurs in graphitic schists in garnet, staurolite and sillimanite zones, the latter being a product of contact metamorphism by HPG. During metamorphism, coexisting fluid contained both CO₂ and CH₄. In the garnet zone, garnet crystals contain petrographically distinct cores with inclusions of quartz, graphite and other minerals. Centres of the cores have distinctly elevated Y concentrations that mark the positions of garnet nucleation. The elevated Y is thought to have come from the Mn‐siderite onto which Y was probably absorbed during precipitation in an ocean. In the upper garnet and staurolite zones, the cores were overgrown by inclusion‐poor mantles. Mantles are highly zoned and have more elevated Fe and Mg and lower Mn and Ca than cores. The growth of mantles is attributed to late‐orogenic heating by leucogranite magmas and attendant influx of H₂O that caused consumption of graphite in rock matrices. A portion of the Proterozoic terrane that includes the HPG is surrounded by four large faults. In this ‘HPG block’, garnet is inclusion‐poor and its composition does not preserve its early growth history. This garnet appears to have re‐equilibrated by internal diffusion of its major components and/or recrystallization of an earlier inclusion‐rich garnet. It has equilibrated within the kyanite stability range, and together with remnant kyanite in the high‐strain aureole of the HPG, indicates that the HPG block had a ≥6 kbar history. The HPG block has undergone decompression during emplacement of the HPG. The decompression is evident in occurrences of retrograde andalusite and cordierite in the thermal aureole of the HPG. The data support a polybaric metamorphic history of the Black Hills orogen with different segments of the orogen having their own clockwise P–T–t paths.     

High‐Precision Cd Isotope Measurements of Soil and Rock Reference Materials by MC‐ICP‐MS with Double Spike Correction

This study presents a high‐precision Cd isotope measurement method for soil and rock reference materials using MC‐ICP‐MS with double spike correction. The effects of molecular interferences (e.g., ¹⁰⁹Ag¹H⁺, ⁹⁴Zr¹⁶O⁺, ⁹⁴Mo¹⁶O⁺ and ⁷⁰Zn⁴⁰Ar⁺) and isobaric interferences (e.g., Pd, In and Sn) to Cd isotope measurements were quantitatively evaluated. When the measured solution has Ag/Cd ≤ 5, Zn/Cd ≤ 0.02, Mo/Cd ≤ 0.4, Zr/Cd ≤ 0.001, Pd/Cd ≤ 5 × 10^?5 and In/Cd ≤ 10^?3, the measured Cd isotope data were not significantly affected. The intermediate measurement precision of pure Cd solutions (BAM I012 Cd, Münster Cd and AAS Cd) was better than ± 0.05‰ (2s) for δ^114/110Cd. The δ^114/110Cd values of soil reference materials (NIST SRM 2709, 2709a, 2710, 2710a, 2711, 2711a and GSS‐1) relative to NIST SRM 3108 were in the range of ?0.251 to 0.632‰, the δ^114/110Cd values of rock reference materials (BCR‐2, BIR‐1, BHVO‐2, W‐2, AGV‐2, GSP‐2 and COQ‐1) varied from ?0.196‰ to 0.098‰, and that of the manganese nodule (NOD‐P‐1) was 0.163 ± 0.040‰ (2s, n = 8). The large variation in Cd isotopes in soils and igneous rocks indicates that they can be more widely used to study magmatic and supergene processes.     

Zinc Isotopic Compositions of NIST SRM 683 and Whole‐Rock Reference Materials

In this study the homogeneity of the zinc isotopic composition in the NIST SRM 683 reference material was examined by measuring the Zn isotopic signature in microdrilled sample powders from two metal nuggets. Zinc was purified using AG MP‐1M resin and then measured by MC‐ICP‐MS. Instrumental mass bias was corrected using the “sample‐standard bracketing” method and empirical external normalisation with Cu doping. After evaluating the potential effects of varying acid mass fractions and different matrices, high‐precision Zn isotope data were obtained with an intermediate measurement precision better than ± 0.05‰ (δ⁶⁶Zn, 2s) over a period of 5 months. The δ⁶⁶Zn_JMC‐Lyon mean values of eighty‐four and fourteen drilled powders from two nuggets were 0.11 ± 0.02‰ and 0.12 ± 0.02‰, respectively, indicating that NIST SRM 683 is a good isotopic reference material with homogeneous Zn isotopes. The Zn isotopic compositions of seventeen rock reference materials were also determined, and their δ⁶⁶Zn values were in agreement with most previously published data within 2s. The δ⁶⁶Zn values of most of the rock reference materials analysed were in the range 0.22–0.36‰, except for GSP‐2 (1.07 ± 0.06‰, n = 12), NOD‐A‐1 (0.96 ± 0.03‰, n = 6) and NOD‐P‐1 (0.78 ± 0.03‰, n = 6). These comprehensive data should serve as reference values for quality assurance and interlaboratory calibration exercises.     

Investigating the Influence of Non‐Spectral Matrix Effects in the Determination of Twenty‐Two Trace Elements in Rock Samples by ICP‐QMS

The influence of non‐spectral matrix effects on the determination of twenty‐two trace elements (Rb, Sr, Y, Cs, Ba, lanthanides, Pb, Th and U) in rock samples using ICP‐MS was investigated. Three types of multi‐element solutions were synthesised containing the twenty‐two trace elements, In, Tl and ten major rock‐forming elements with varying mass fractions mimicking the compositions of basalt, peridotite and dolomite. The synthetic solutions were conditioned to have dilution factors (DF) of 1000–10000. The extent of sensitivity suppression relative to the DF = 10000 solution became more significant for smaller DF solutions, which was not constant across different elements in a single solution but displayed general dependence on m/z. This indicates that at least two internal standards (e.g., In and Tl) are required for the correction of sensitivity variation. On the basis of the results, a new isotope dilution‐internal standardisation method for the determination of twenty‐two trace elements with ICP‐MS was developed, in which the sensitivity variation was corrected by monitoring two enriched isotopes, ¹¹³In and ²⁰³Tl. This method, coupled with the quantitative correction of interference from oxides and hydroxides, achieved precise determination of twenty‐two trace elements in some rock reference materials with reproducibilities of ±2% for basaltic to andesitic samples.     

Garnet re‐equilibration by coupled dissolution–reprecipitation: evidence from textural,major element and oxygen isotope zoning of ‘cloudy’ garnet

The analysis of texture, major element and oxygen isotope compositions of cloudy garnet crystals from a metapelite sampled on Ikaria Island (Greece) is used to assess the model of growth and re‐equilibration of these garnet crystals and to reconstruct the pressure–temperature–fluid history of the sample. Garnet crystals show complex textural and chemical zoning. Garnet cores (100–200 μm) are devoid of fluid inclusions. They are characterized by growth zoning demonstrated by a bell‐shaped profile of spessartine component (7–3 mol.%), an increase in grossular from 14 to 22 mol.% and δ¹⁸O values between 9.5 ± 0.3‰ and 10.4 ± 0.2‰. Garnet inner rims (90–130 μm) are fluid inclusion‐rich and show a decreasing grossular component from 22 to 5 mol.%. The trend of the spessartine component observed in the inner rim allows two domains to be distinguished. In contrast to domain I, where the spessartine content shows the same trend as in the core, the spessartine content of domain II increases outwards from 2 to 14 mol.%. The δ¹⁸O values decrease towards the margins of the crystals to a lowest value of 7.4 ± 0.2‰. The outer rims (<10 μm) are devoid of fluid inclusions and have the same chemical composition as the outermost part of domain II of the inner rim. Garnet crystals underwent a four‐stage history. Stage 1: garnet growth during the prograde path in a closed system for oxygen. Garnet cores are remnants of this growth stage. Stage 2: garnet re‐equilibration by coupled dissolution–reprecipitation at the temperature peak (630 < T < 650 °C). This causes the creation of porosity as the coupled dissolution–reprecipitation process allows chemical (Ca) and isotopic (O) exchange between garnet inner rims and the matrix. The formation of the outer rim is related to the closure of porosity. Stage 3: garnet mode decreases during the early retrograde path, but garnet is still a stable phase. The resulting garnet composition is characterized by an increasing Mn content in the inner rim’s domain II caused by intracrystalline diffusion. Stage 4: dissolution of garnet during the late retrograde path as garnet is not a stable phase anymore. This last stage forms corroded garnet. This study shows that coupled dissolution–reprecipitation is a possible re‐equilibration process for garnet in metamorphic rocks and that intra‐mineral porosity is an efficient pathway for chemical and isotopic exchange between garnet and the matrix, even for otherwise slow diffusing elements.     

An Eocene/Oligocene blueschist‐/greenschist facies P–T loop from the Cycladic Blueschist Unit on Naxos Island,Greece: Deformation‐related re‐equilibration vs. thermal relaxation

Geothermobarometric and geochronological work indicates a complete Eocene/early Oligocene blueschist/greenschist facies metamorphic cycle of the Cycladic Blueschist Unit on Naxos Island in the Aegean Sea region. Using the average pressure–temperature (P–T) method of thermocalc coupled with detailed textural work, we separate an early blueschist facies event at 576 ± 16 to 619 ± 32°C and 15.5 ± 0.5 to 16.3 ± 0.9 kbar from a subsequent greenschist facies overprint at 384 ± 30°C and 3.8 ± 1.1 kbar. Multi‐mineral Rb–Sr isochron dating yields crystallization ages for near peak‐pressure blueschist facies assemblages between 40.5 ± 1.0 and 38.3 ± 0.5 Ma. The greenschist facies overprint commonly did not result in complete resetting of age signatures. Maximum ages for the end of greenschist facies reworking, obtained from disequilibrium patterns, cluster near c. 32 Ma, with one sample showing rejuvenation at c. 27 Ma. We conclude that the high‐P rocks from south Naxos were exhumed to upper mid‐crustal levels in the late Eocene and early Oligocene at rates of 7.4 ± 4.6 km/Ma, completing a full blueschist‐/greenschist facies metamorphic cycle soon after subduction within c. 8 Ma. The greenschist facies overprint of the blueschist facies rocks from south Naxos resulted from rapid exhumation and associated deformation/fluid‐controlled metamorphic re‐equilibration, and is unrelated to the strong high‐T metamorphism associated with the Miocene formation of the Naxos migmatite dome. It follows that the Miocene thermal overprint had no impact on rock textures or Sr isotopic signatures, and that the rocks of south Naxos underwent three metamorphic events, one more than hitherto envisaged.     

The opening of the Plio-Quaternary Gibraltar Strait: assessing the size of a cataclysm

Abstract

The widespread Lago-Mare facies show that the Mediterranean was fully isolated from the World Ocean at the end of the Messinian salinity crisis (MSC): the Plio-Quaternary Gibraltar Strait has no geographical feature in common with the Miocene Portals, which were continentalised before the end of the MSC. A copy of a budget model has been modified to study the infilling of the Mediterranean. The yearly Atlantic water supply follows an exponential increase law. The level of the basins does not change significantly during the 26 early years, but the refill of the Mediterranean basins is completed within the next 10–11 years. The high velocity of the flow of Atlantic water can be reconfirmed by hydraulic calculations. The morphological study of the area suggests that the Gibraltar Strait originates from an eastwards flowing stream, which drained the eastern slope of an emergent Gibraltar Isthmus to the Mediterranean Basins at the end of the Messinian. Similar canyons are still found in the marine topography towards the Alboran Sea. The present Spartel and Camarinal Sills both originate from the coeval mass sliding of the northern bank of the strait. The strait of Gibraltar results from the work of a stream, which managed water-piracy from an Ocean. © 2002 Éditions scientifiques et médicales Elsevier SAS. All rights reserved.     

Evaluation of Major to Ultra Trace Element Bulk Rock Chemical Analysis of Nanoparticulate Pressed Powder Pellets by LA‐ICP‐MS

An efficient, clean procedure for the measurement of element mass fractions in bulk rock nanoparticulate pressed powder pellets (PPPs) by 193 nm laser ablation ICP‐MS is presented. Samples were pulverised by wet milling and pelletised with microcrystalline cellulose as a binder, allowing non‐cohesive materials such as quartz or ceramics to be processed. The LA‐ICP‐MS PPP analytical procedure was optimised and evaluated using six different geological reference materials (JP‐1, UB‐N, BCR‐2, GSP‐2, OKUM and MUH‐1), with rigorous procedural blank quantification employing synthetic quartz. Measurement trueness of the procedure was equivalent to that achieved by solution ICP‐MS and LA‐ICP‐MS analysis of glass. The measurement repeatability was as low as 0.5–2% (1s, n = 6) and, accordingly, PPP homogeneity could be demonstrated. Calibration based on the reference glasses NIST SRM 610, NIST SRM 612, BCR‐2G and GSD‐1G revealed matrix effects for glass and PPP measurement with NIST SRM 61×; using basalt glasses eliminated this problem. Most significantly, trace elements not commonly measured (flux elements Li, B; chalcophile elements As, Sb, Tl, In, Bi) could be quantified. The PPP‐LA‐ICP‐MS method overcomes common problems and limitations in analytical geochemistry and thus represents an efficient and accurate alternative for bulk rock analysis.