the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 04 Jul 2025
Seismic surveys at the Per Geijer iron oxide-apatite mineralization in Kiruna, Sweden
Abstract. LKAB currently runs an intense exploration program focused on the Per Geijer (PG) iron oxide-apatite mineralization, close to the active Kiruna mine, in northern Sweden. The mineralization and the hosting volcanic units dip at an angle of 50–60 degrees and are overlain by quartzite. Greenstones are also present in the area. The iron mineralization is associated with high levels of phosphorus and REEs. As a major part of the mineralization occurs at depths >500 m, various geophysical surveys have been used to assist the exploration throughout the years, recently also including seismic measurements. An initial pilot reflection profile in 2021 showed promising results, and three additional profiles were acquired shortly afterwards, all with Vibroseis as a source. These were followed by downhole seismics (VSP) in two drill holes. During 2024 a 3D survey covering the Per Geijer mineralization was also acquired. Clear P-wave reflections are visible in the shot records from the surface seismics as well as in the downhole seismics, and the latter also show interpreted P-S reflections. The steeply dipping structures in the area complicate the data processing, and out-of-the-plane reflections are present in some of the 2D profiles, the latter confirmed via cross-dip analysis. The resulting migrated sections and volumes show prominent reflections, some of which correlate with the Per Geijer mineralization or units close to it and suggest some extensions of the known mineralization. Reflections are also observed from the northern end of the nearby Kiirunavaara ore. Sonic and density logs are available to support the interpretation. These indicate that the main cause of the reflections from the iron mineralization is the density contrast, rather than the velocity contrasts. The downhole seismics also allow a closer link to the drilling information. Further notable in both our 2D and 3D surface data is a high-amplitude reflection originating from a steeply dipping structure outside the survey area, which if projected to the surface correlates with older greenstones on the western side of the Kiirunavaara orebody.
- Preprint
(3750 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on egusphere-2025-3110', Anonymous Referee #1, 31 Jul 2025
- AC1: 'Reply on RC1', Niklas Juhojuntti, 16 Sep 2025
-
RC2: 'Comment on egusphere-2025-3110', Anonymous Referee #2, 06 Aug 2025
This paper presents a result of seismic studies of iron oxide-apatite mineralisation in Kiruna mine area of Northern Sweden. The paper is well structured and describes 2D and 3D surface seismic studies, downhole seismic data recorded with DAS and key rock physics properties defining the host and ore formations.
Limitations in seismic acquisition and poor seismic data quality, that is typical for crystalline environment, resulted in fragmented subsurface picture. Nevertheless, the integration and interpretation section of the studied geophysical data presents convincing analysis of iron mineralisation. Although, there is a great room for alternative interpretations.
Overall, I would recommend this paper for publication but would like to see more detailed comparison between the three types of seismic data at the pre-image level. One of the important decisions to make before any further seismic exploration here or in similar geological environments is the trade-off between the value and the cost of the presented seismic methods.
Assessment of the 2D, 3D and VSP gathers from close by locations, their respective frequency spectra and signal-to-noise ratio at raw and pre-processed stages will be helpful in evaluation of the significance of each dataset.
Some minor comments and suggestions:
Line 51: typo: accelerated wight-drop – should be - accelerated weight-drop
Line 120: typo: source have been used - source has been used
Line 141: However, as both holes were waterfilled the coupling should have been adequate.
DAS is sensitive to axial strain, so coupling via water is not adequate for shear contact of the cable with the media.
Figure 4: add arrows pointing at the reflections discussed in the text
Also, this figure is a good place to add illustrations of the 2D, 3D and VSP gathers from close by locations, their respective frequency spectra and signal-to-noise ratio at raw and pre-processed stages.
Line 181: more moderately dipping or sub-horizontal along Profiles 1 and 3 (Fig. 5) - there is no profile 3 in Fig 5
Figure 5: it will be more convenient for readers to see West on the left and East on the right side for profile 2 - to be consistent with the map.
Line 186: Note that the fold is varying much along Profile 2, which is clearly noticeable in the stacked time section (Fig. 5). - no CMP fold presented in the figure or elsewhere – need to add if refereeing to it in the text.
Also, as Profie2 is bent, the image deficiency might be due to the geometry and, hence, it is good to see CDP location at the bending point or schematic plot of the profile XY geometry on top pf the section to have a reference.
198: The reflection at about distance 4 km along the profile, and at a depth 1.5 km,
coincides with the Per Geijer Deep mineralization – indicate by arrows in the figure.
Line 217: All surface seismic data have been depth-converted using a constant velocity of 6000 m/s.
Why velocity used for surface seismic is 6000m/s here, but it was 5000m/s for the downhole processing (Line 243)?
Line 243: A velocity of 5000 m/s was chosen for the migration, and the inherent depth conversion.
Figure 7: East-West indications are missing; similar to profile 2 in Fig 5, the section should be flipped to be consistent with the map orientation: West - left, East – right.
Figure 8: West-East orientation is good here, should be the same in Fig 5 and Fig 7.
It would be good to have an example of inline from migrated 3D depth volume coinciding (or located as close as possible) with profile 2 for one-to-one comparison of the 2D and 3D products.
Figure 13A: indicate reflection from the ‘2.40-1.96 Ga basalt-andesite units’
Show amplitude spectra for all datasets 2D, 3D, VSP for raw and pre-processed, just before stack (or before migration)
Appendix:
Bandpass filter: 10-20-150-220 Hz
Bandpass filter: 6-12-130-200 Hz
Bandpass filter: 7-14-50-80 H - why frequency content is so different for 3D processing post-stack in comparison to the other two?
Citation: https://doi.org/10.5194/egusphere-2025-3110-RC2 -
AC2: 'Reply on RC2', Niklas Juhojuntti, 16 Sep 2025
I have addressed most of the comments from the reviewer, and revised the manuscript accordingly. The exception is the wish from the reviewer to include a comparison of the different survey types (2D, 3D and VSP), which I think is outside the scope of the paper.
-
AC2: 'Reply on RC2', Niklas Juhojuntti, 16 Sep 2025
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 760 | 52 | 21 | 833 | 40 | 46 |
- HTML: 760
- PDF: 52
- XML: 21
- Total: 833
- BibTeX: 40
- EndNote: 46
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
THe manuscript is almost ready for publication, however it would benefit from a minor revision, including improving the english.
Authors should take advantage of the figure captions by fully explaining what is in the figure and what is the reader supose to look for. In the current form most of figure captions are to short to fully follow the authors reasoning.
Some minor coments are also included in the text attached.