<p>We investigate seismic velocity changes in the rock mass related to mining induced seismic events and ore exploitation by computing a one-month long 4D elastic model of Kiirunavaara mine (Sweden). We focus on a specific mine sector, where a single <InlineEquation ID="IEq1"> <EquationSource Format="TEX">\(\varvec{M_W}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <msub> <mi mathvariant="bold-italic">M</mi> <mi mathvariant="bold-italic">W</mi> </msub> </mrow> </math></EquationSource> </InlineEquation>=2.0 event occurred on May 22 (02:31 local time), damaging the infrastructure. We make use of P- and S-first-arrival times obtained from the permanent seismic system for computing the full 4D (continuous 3D volume in time) seismic velocity model of Kiruna mine using a trans-dimensional Monte Carlo sampling. The trans-dimensional approach guarantees that the resolution, both in space and in time, is strictly data-driven. Our results give the following insights into the velocity differences at the mining levels and at different time-length scales. (a) We observe a striking correlation between spatial variations of <InlineEquation ID="IEq2"> <EquationSource Format="TEX">\(\varvec{V_P}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <msub> <mi mathvariant="bold-italic">V</mi> <mi mathvariant="bold-italic">P</mi> </msub> </mrow> </math></EquationSource> </InlineEquation> and ore-body geometry, confirming the robustness of the velocity model. Clay zones appear as a low <InlineEquation ID="IEq3"> <EquationSource Format="TEX">\(\varvec{V_P/V_S}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <msub> <mi mathvariant="bold-italic">V</mi> <mi mathvariant="bold-italic">P</mi> </msub> <mo mathvariant="bold" stretchy="false">/</mo> <msub> <mi mathvariant="bold-italic">V</mi> <mi mathvariant="bold-italic">S</mi> </msub> </mrow> </math></EquationSource> </InlineEquation> ratio zones, as seen in previous tomographic studies. (b) High-frequency (hourly) fluctuations of the rock mass <InlineEquation ID="IEq4"> <EquationSource Format="TEX">\(\varvec{V_P}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <msub> <mi mathvariant="bold-italic">V</mi> <mi mathvariant="bold-italic">P</mi> </msub> </mrow> </math></EquationSource> </InlineEquation> around the ore-passes are highly correlated with seismic sequences in the same rock volumes. In particular, <InlineEquation ID="IEq5"> <EquationSource Format="TEX">\(\varvec{V_P}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <msub> <mi mathvariant="bold-italic">V</mi> <mi mathvariant="bold-italic">P</mi> </msub> </mrow> </math></EquationSource> </InlineEquation> increases rapidly when ore-passes are seismically active and <InlineEquation ID="IEq6"> <EquationSource Format="TEX">\(\varvec{V_P}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <msub> <mi mathvariant="bold-italic">V</mi> <mi mathvariant="bold-italic">P</mi> </msub> </mrow> </math></EquationSource> </InlineEquation> values keep a high value for few (1-4) hours after the end of the seismic sequence. (c) The smoothed velocity model, computed as averaged model over a 2-days moving window, suggests that low-frequency <InlineEquation ID="IEq7"> <EquationSource Format="TEX">\(\varvec{V_P}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <msub> <mi mathvariant="bold-italic">V</mi> <mi mathvariant="bold-italic">P</mi> </msub> </mrow> </math></EquationSource> </InlineEquation> fluctuations can be compared to stress cell measurements located closely.</p>