A 10.5 km² 3D seismic survey was acquired over the Kylylahti mine area (Outokumpu mineral district, eastern Finland) as a part of the COGITO-MIN (COst-effective Geophysical Imaging Techniques for supporting Ongoing MINeral exploration in Europe) project, which aimed at the development of cost-effective geophysical imaging methods for mineral exploration. The cost-effectiveness in our case was related to the fact that an active-source 3D seismic survey was accomplished by using the receiver spread originally designed for a 3D passive survey. The 3D array recorded Vibroseis and dynamite shots from an active-source 2D seismic survey, from a vertical seismic profiling experiment survey, as well as some additional “random” Vibroseis and dynamite shots made to complement the 3D source distribution. The resulting 3D survey was characterized by irregular shooting geometry and relatively large receiver intervals (50 m). Using this dataset, we evaluate the effectiveness of the standard time-imaging approach (post-stack and pre-stack time migration) compared to depth imaging (standard and specialized Kirchhoff pre-stack depth migration, KPreSDM). Standard time-domain processing and imaging failed to convincingly portray the first ~1500 m of the subsurface, which was the primary interest of the survey. With a standard KPreSDM, we managed to obtain a good image of the base of the Kylylahti formation bordering the extent of the mineralization-hosting Outokumpu assemblage rocks, but otherwise the image was very noisy in the shallower section. The specialized KPreSDM approach (i.e., coherency-based Fresnel volume migration) resulted in a much cleaner image of the shallow, steeply dipping events, as well as some additional deeper reflectors, possibly representing repetition of the contact between the Outokumpu assemblage and the surrounding Kalevian metasediments at depth.

A sparse 3D seismic survey was acquired over the Blötberget iron oxide deposits of the Ludvika Mines in south-central Sweden. The main aim of the survey was to delineate the deeper extension of the mineralisation and to better understand its 3D nature and associated fault systems for mine planning purposes. To obtain a high-quality seismic image in depth, we applied time-domain 3D acoustic full-waveform inversion (FWI) to build a high-resolution P-wave velocity model. This model was subsequently used for pre-stack depth imaging with reverse time migration (RTM) to produce the complementary reflectivity section. We developed a data preprocessing workflow and inversion strategy for the successful implementation of FWI in the hardrock environment. We obtained a high-fidelity velocity model using FWI and assessed its robustness. We extensively tested and optimised the parameters associated with the RTM method for subsequent depth imaging using different velocity models: a constant velocity model, a model built using first-arrival travel-time tomography and a velocity model derived by FWI. We compare our RTM results with a priori data available in the area. We conclude that, from all tested velocity models, the FWI velocity model in combination with the subsequent RTM step provided the most focussed image of the mineralisation and we successfully mapped its 3D geometrical nature. In particular, a major reflector interpreted as a cross-cutting fault, which is restricting the deeper extension of the mineralisation with depth, and several other fault structures which were earlier not imaged were also delineated. We believe that a thorough analysis of the depth images derived with the combined FWI–RTM approach that we present here can provide more details which will help with better estimation of areas with high mineralisation, better mine planning and safety measures.

Seismic imaging is now a well-established method in mineral exploration with many successful case studies. Seismic data are usually imaged in the time domain (post-stack or pre-stack time migration), but recently pre-stack depth imaging has shown clear advantages for irregular/sparse acquisitions and very complex targets. Here, we evaluate the effectiveness of both ray-based and wave-equation-based pre-stack depth imaging methodologies applied to crooked-line 2D seismic reflection profiles. Seismic data were acquired in the Kylylahti mining area in eastern Finland over severely folded, faulted and subvertical Kylylahti structure, and associated mineralization. We performed 3D ray-based imaging, i.e., industry-standard Kirchhoff migration and its improved version (coherency migration, CM), and wave-equation-based migration, i.e., reverse time migration (RTM) using a velocity model built from first-arrival traveltime tomography. Upon comparing the three different migrations against available geological data and models, it appeared that CM provided the least noisy and well-focused image, but failed to image the internal reflectivity of the Kylylahti formation. RTM was the only method that produced geologically plausible reflections inside the Kylylahti formation including a direct image of the previously known shallow massive sulfide mineralization.

The Seismic and Electromagnetics Methods for Deep mineral exploration (SEEMS DEEP) project is associated with the Koillismaa Layered Igneous Complex (KLIC) in north-eastern Finland. The KLIC is characterized by a Bouguer positive gravity and magnetic anomaly zone connecting the two exposed ends of the KLIC, i.e. the Koillismaa intrusion and the Näränkävaara intrusion. The KLIC has the potential to host several critical raw minerals like nickel and cobalt, which are in the European Union's critical raw material list. For this purpose, two reflection seismic profiles were acquired to map the regional reflectivity in the area with a focus on imaging the large-scale geological architecture of the KLIC and the associated mineralization. Seismic imaging delineated reflectivity up to a depth of approximately 5 to 6 km with several distinctive reflective packages at various depths representative of the presence of dykes, faults, and major lithological contacts in the area. Symmetrical orientation of several shallow reflectors indicates systematic fracturing in the area. The deep-seated magma conduit associated with the KLIC was successfully mapped, revealing a complex internal structure with more structural or compositional heterogeneity than was previously assumed to be formed from a single lithological unit. The magma conduit is characterised as a vertically extensive mafic–ultramafic package with the potential of hosting orthomagmatic critical raw material deposits of mostly Ni-Cu-PGE sulphide type. A system of overthrust – steep regional faults was imaged for the first time in the western part of the study area. It is interpreted that such steep fault systems might have played a key role in magma emplacement by acting as pathways for ascending melts and mineralising fluids. Results were compared against the available petrophysical data and a preliminary geological model based on the density model obtained from the gravity inversion with constraints from the drillhole data.

A sparse 3D seismic survey was acquired over the Blötberget iron-oxide deposits of the Ludvika Mines in south-central Sweden. The main aim of the survey was to delineate the deeper extension of the mineralisation and to better understand its 3D nature and associated fault systems for mine planning purposes. To obtain a high-quality seismic image in depth, we applied time-domain 3D acoustic full-waveform inversion (FWI) to build a high-resolution P-wave velocity model. This model was subsequently used for pre-stack depth imaging with reverse time migration (RTM) to produce the complementary reflectivity section. We developed a data preprocessing workflow and inversion strategy for the successful implementation of FWI in the hardrock environment. We obtained a high-fidelity velocity model using FWI and assessed its robustness. We extensively tested and optimised the parameters associated with the RTM method for subsequent depth imaging using different velocity models: a constant velocity model, a model built using first-arrival traveltime tomography and a velocity model derived by FWI. We compare our RTM results with a priori data available in the area. We conclude that, from all tested velocity models, the FWI velocity model in combination with the subsequent RTM step, provided the most focussed image of the mineralisation and we successfully mapped its 3D geometrical nature. In particular, a major reflector interpreted as a cross-cutting fault, which is restricting the deeper extension of the mineralisation with depth, and several other fault structures which were earlier not imaged were also delineated. We believe that a thorough analysis of the depth images derived with the combined FWIRTM approach that we presented here can provide more details which will help with better estimation of areas with high mineralization, better mine planning and safety measures.

A pioneering ~22 km² three-dimensional (3D) seismic survey was conducted in northeastern Finland within the SEEMS DEEP project to investigate the concealed Koillismaa Deep Intrusion (KDI), part of the 2.5–2.4 Ga Koillismaa–Näränkävaara Layered Intrusion Complex (KLIC). Utilizing low-cost piezoelectric nimble nodes and a single Vibroseis source, this survey represents one of the first 3D seismic investigations of a deep layered intrusion and its feeder system. The seismic data, supported by petrophysical measurements and synthetic modelling, reveal a complex internal architecture characterized by coherent reflectivity patterns indicative of modal layering. Interpretation suggests that the KDI is not a simple feeder but a more developed, chonolithic or funnel-tube-shaped intrusion. Integration of seismic data interpretation with gravity inversion and geological data enabled the construction of a Common Earth Model (CEM), which successfully reproduces the observed gravity anomaly. A major thrust fault, likely associated with the Oulujärvi Shear Zone, was imaged for the first time, offering new insights into the tectonic evolution and mineralization potential of the region. This study highlights the advances in 3D seismic methods in hardrock environment and establishes a methodological framework for future mineral exploration in similar geological settings.