European iron ore production is primarily sourced from magnetite dominated iron oxide-apatite ore deposits in the northern Norrbotten ore province of northernmost Sweden. The Malmberget iron oxide-apatite deposit is at present the largest iron ore resource in Europe and is an amphibolite facies grade analogue of the world-famous Kiirunavaara iron oxide-apatite deposit. The Malmberget rock association is characterised by multiple phases of deformation, metamorphism, and alteration that resulted in a geometrically and petrologically complex deposit that is genetically ambiguous. Primary ore textures and emplacement structures of the Malmberget iron oxide-apatite deposit have largely been recrystallised during metamorphic overprinting and now comprise dominantly medium- to coarse-grained granoblastic magnetite. However, isotopic characteristics are preserved and when combined with trace element chemistry, these can be used to understand magmatic vs. hydrothermal origin of the deposit. To unravel the primary origin of the Malmberget magnetite ore, we combined magnetite trace element chemistry and Fe-O stable isotopes to investigate the massive magnetite in the Fabian-Kapten and ViRi ore bodies of the Malmberget iron oxide-apatite deposit. Trace element correlations indicate a high-temperature magmatic to a transitional high-temperature magmatic-hydrothermal origin of the Malmberget iron oxide-apatite ore deposit, with data plotting into fields of clear magmatic affinity in trace element discrimination diagrams. Fe-O data fall into established magmatic fields regardless of subsequent metamorphic modifications, underlining a dominantly (ortho-)magmatic origin of the investigated deposits. Despite an overall magmatic to magmatic-hydrothermal origin for the two ore bodies studied, Fe-O isotope equilibrium calculations of the magnetite suggest a possible temperature discrepancy between the Fabian-Kapten ore body and the ViRi ore body, the latter showing a more pronounced magmatic character. These variations in trace element contents and Fe-O isotopes can be explained by the proximity of the more magmatic signatures to the centre of the ore forming magmatic system.

The central Ossa Morena Zone (SW Iberia) hosts a regionally extensive ironstone level interbedded with bimodal volcanic rocks, limestone and shale of Lower-Middle Cambrian age. The stratabound ironstone includes dominant magnetite and hematite with locally abundant chert and barite. It is interpreted as being (sub-)exhalative at or near the seafloor and formed during a rifting event that postdated the Cadomian orogeny. In some places, such as in the Las Herrerías deposit, the ironstone is irregularly replaced by a chalcopyrite-rich ore; the Cu-rich mineralization is accompanied by the pervasive phyllic alteration of the hosting siliciclastic sediments. The highest copper grades are found when the ironstone is crosscut by WNW-ESE-trending late-Variscan extensional brittle-ductile structures that are interpreted as the feeder channels for deep hydrothermal fluids. A similar nearby copper-rich mineralization (Pallares) is is likely controlled by the tectonic contact between limestone and pyrite-rich black shale.

Sr-Nd whole-rock isotope geochemistry data suggests that the Sr in the ironstone (⁸⁷Sr/⁸⁶Sr_i ≈ 0.7088) is close to isotopic equilibrium with the local exhalative barite (0.7084–0.7086) and Cambrian seawater. The ironstone has a significantly more crustal εNd initial signature (<-1.8) than the coeval volcanic rocks (+5.2 to + 7.9). The younger sulfide mineralization inherited the Nd isotope composition of the ironstone but shows a significant enrichment in ⁸⁷Sr (⁸⁷Sr/⁸⁶Sr > 0.7091) that is interpreted as related with the input of genetically different and more crustally-derived hydrothermal fluids.

³⁹Ar-⁴⁰Ar dating of the phyllic alteration suggest that the copper mineralization was formed at ca. 332–330 My. These ages are coeval with those of small peraluminous granite intrusions that host Cu-Au vein-like mineralization and dated at 331.8 ± 1.6 Ma (LA ICPMS U-Pb zircon). Our interpretation is that the copper-rich mineralization at the Las Herrerías area is the distal expression of an intrusion-related hydrothermal system.

Numerical modelling shows that ironstone is an effective trap for copper precipitation due to the large changes in pH and fO2 that take place when copper-bearing acid and reduced fluids react with the brittle ironstone. The precipitation of chalcopyrite, however, is controlled by the amount of available reduced sulfur in the ore trap. The δ³⁴S values of the sulfides (+12.6 to + 21.6 ‰) suggest that the most likely source for the reduced sulfur is the thermogenic reduction of aqueous sulfate equilibrated with the exhalative barite (δ34S, +31.4 to + 35 ‰) with some minor input of reduced sulfur leached from the metasediments.This system could be considered as a variant of the IOCG clan. The formation of the ironstone and the copper mineralization, however, are separated by more than 200 My. Probably, many IOCG systems have a similar origin as Las Herrerías, with an ironstone being just a passive geochemical trap for the copper–gold mineralization.

To support economic decisions and exploration targeting, as well as to understand processes controlling the mineralization, three-dimensional structural and lithological boundary models of the Kiruna mining district have been built using surface (outcrop observations and measurements) and subsurface (drill hole data and mine wall mapping) data. Rule-based hybrid implicit-explicit modeling techniques were used to create district-scale models of areas with high disproportion in data resolution characterized by dense, clustered, and distant data spacing. Densely sampled areas were integrated with established conceptual studies using geologic conditions and the addition of synthetic data, leading to variably constrained surfaces that facilitate the visualization, interpretation, and further integration of the geologic models. This modeling approach proved to be efficient in integrating local, frequently sampled areas with district-scale, sparsely sampled regions. Dominantly S-plunging lineation on N-S–trending fracture planes, characteristic fracture mineral fill, and weak rock mass at the ore contact indicated by poor core orientation quality and rock quality description suggest that ore-parallel fractures in the Kiirunavaara area were more commonly reactivated. Slight variation in the angular relationship of fracture sets situated in different fault-bounded blocks suggests that strain accommodation across the orebodies was uneven. The location of brittle faults identified in drill core, deposit-scale structural analysis, and aeromagnetic geophysical maps indicate a close relationship between fault locations and the iron oxide-apatite mineralization, suggesting that uneven stress accommodation and proximity of conjugate fault sets played an important role in juxtaposing blocks from different crustal depths and control the location of the iron oxide-apatite orebodies.

The northern Norrbotten ore province in Sweden is one of the most mineralized areas in Europe. Iron, gold and/or copper deposits occur as iron oxide-apatite-style mineralization (IOA, Kiruna-type) as well as iron oxide-copper–gold (IOCG) style. Regardless of mineralization style, most deposits appear to be spatially controlled by a set of crustal-scale Palaeoproterozoic shear zones which share similar structural characteristics and deformation histories.

Reappraisal of regional geological and geophysical data, coupled with structural mapping, suggests crustal-scale shear zones form continuous c. N-S-trending zones extending from the Skellefte district in the south into the northern Norrbotten ore province. One example from Norrbotten is a zone that extends SSW from Karesuando in the north towards Svappavaara. While this structure has traditionally been inferred to continue SW towards Arjeplog (i.e. the Karesuando – Arjeplog Deformation Zone; KADZ), we favour its deflection SSE into the Nautanen-Aitik trend, making it a continuous, IOCG-bearing, crustal-scale deformation zone. Similar shear zone geometries can be observed in analogous zones to the west. Most of these crustal scale structures record at least two time-separated deformation events of regional significance. IOA and IOCG deposits form in different tectonic environments, separated in time and overprinting each other.

As part of the larger mineral systems approach to Cu-bearing mineralization in northern Norrbotten, this study utilizes structural geology to set the classic Pahtohavare Cu ± Au deposits into an up-to-date tectonic framework. The Pahtohavare Cu ± Au deposits, situated only 5 km southwest of the Kiirunavaara world-class iron oxide–apatite (IOA) deposit, have a dubious timing, and their link to IOA formation is not constrained. The study area contains both epigenic Cu ± Au (Pahtohavare) and iron oxide–copper–gold (IOCG; Rakkurijärvi) mineral occurrences which are hosted in bedrock that has been folded and bound by two shear zones trending northeast to southwest and northwest to southeast to the east and southwest, respectively. Structural mapping and petrographic investigation of the area reveal a noncylindrical, SE-plunging anticline. The cleavage measurements mirror the fold geometry, which characterizes the fold as F₂ associated with the late phase of the Svecokarelian orogeny. Porphyroclasts with pressure shadows, mylonitic fabrics, and foliation trails in porphyroblasts indicate S₀/S₁  is a tectonic fabric. The epigenetic Pahtohavare Cu ± Au mineralization sits in brittle–ductile structures that cross-cut an earlier foliation and the F₂ fold, indicating that the timing of the deposits occurred syn- to post-F₂ folding, at least ca. 80 Myr after the Kiirunavaara IOA formation. A 3D model and cross-sections of the Pahtohavare–Rakkurijärvi area and a new structural framework of the district are presented and used to suggest that the shear zones bounding the area are likely reactivated early structures that have played a critical role in ore formation in the Kiruna mining district.

The Kiruna mining district, Sweden, known for the type locality of Kiruna-type iron oxide-apatite (IOA) deposits, also hosts several Cu-mineralized deposits including iron oxide-copper-gold (IOCG), exhalative stratiform Cu-(Fe-Zn), and structurally controlled to stratabound Cu +/- Au. However the relationship between the IOA and Cu-systems has not been contextualized within the regional tectonic evolution. A broader mineral systems approach is taken to assess the timing of energy drive(s) within a regional tectonic framework by conducting U-Pb zircon geochronology on intrusions from areas where Cu-mineralization is spatially proximal. Results unanimously yield U-Pb ages from the early Svecokarelian orogeny (ca. 1923-1867 Ma including age uncertainties), except one sample from the Archean basement (2698 +/- 3 Ma), indicating that a distinct thermal drive from magmatic activity was prominent for the early orogenic phase. A weighted average Pb-207/Pb-206 age of 1877 +/- 10 Ma of an iron-oxide-enriched gabbroic pluton overlaps in age with the Kiirunavaara IOA deposit and is suggested as a candidate for contributing mafic signatures to the IOA ore. The results leave the role of a late energy drive (and subsequent late Cu-mineralization and/or remobilization) ambiguous, despite evidence showing a late regional magmatic-style hydrothermal alteration is present in the district.

The Nautanen deformation zone in the Gällivare area of northern Sweden is a highly Cu-mineralized, magnetite-rich, large-scale shear zone with a long-lived (~100 m.y.) deformation, hydrothermal alteration, and mineralization history. This composite structure hosts the Aitik porphyry Cu-Au-Ag ± Mo deposit and several Cu-Au ± Fe ± Ag ± Mo occurrences assigned to the iron oxide copper-gold (IOCG) deposit class. The Nautanen deformation zone was a locus for polyphase deformation and intermittent metasomatic-hydrothermal activity that overprinted middle Orosirian (ca. 1.90–1.88 Ga) continental arc-related volcanic-plutonic rocks. The deformation zone is characterized by intense shearing fabrics that form a series of subvertical to moderately W-dipping, NNW-SSE–trending, first-order shear zones with oblique reverse kinematics and related NNE-SSW–oriented second-order shear zones that control hydrothermal alteration patterns and Cu-Au mineralization.

Hydrothermal alteration in the study area formed during several phases. Volcanic-volcaniclastic rocks to the east and west of the Nautanen deformation zone display low to moderately intense, pervasive to selectively pervasive (i.e., patchy zones or bands, disseminations) sericite ± feldspar, amphibole + biotite + magnetite ± tourmaline, and K-feldspar + hematite alteration. Both the amphibole + biotite and K-feldspar + hematite associations occur adjacent to NNW- and NE-oriented deformation zones and are locally associated with minor sulfide. Within the deformation zone, a moderate to intense biotite + amphibole + garnet + magnetite + tourmaline + sericite alteration assemblage is typically associated with chalcopyrite + pyrrhotite + pyrite and forms linear and subparallel, mainly NNW-oriented seams, bands, and zones that locally appear to overprint possibly earlier scapolite + sericite ± feldspar alteration. Late-stage epidote ± quartz ± feldspar alteration (retrograde saussuritization) forms selectively pervasive zones and epidote veinlets across the area and is partly related to brittle faulting.

A magnetite-amphibole-biotite–rich, penetrative S₁ foliation records shortening during early Svecokarelian-related deformation (D₁) and can be related to ca. 1.88 to 1.87 Ga arc accretion processes and basin inversion that overlaps with regional peak metamorphism to near mid-amphibolite facies conditions and a potential initial Cu mineralization event. Folding and repeated shearing along the Nautanen deformation zone can be assigned to a second, late-Svecokarelian deformation event (D₂ stage, ca. 1.82–1.79 Ga) taking place at a higher crustal level. This D₂ deformation phase is related to late-stage accretionary processes active during a transition to a stage of postorogenic collapse, and it was accompanied by abundant, syntectonic intrusions. D₂-related magmatism produced high-temperature and low-pressure conditions and represents a regional magmatic-hydrothermal event that controlled the recrystallization/remobilization of magnetite, biotite, and amphibole. Associated shear zone reactivation during D₂ favors the utilization of the Nautanen deformation zone as a fluid conduit, which preferentially controlled the siting and formation of epigenetic Cu-Au mineralization with distinctive IOCG characteristics within second-order shear zones.