=== versions ===
python 3.12.0 | numpy 2.5.4 | scipy 1.18.1 | rasterio 1.5.2 | matplotlib 3.11.2 | cmcrameri 1.10

=== inputs ===
read insar-volcano-deformation-python\data\svartsengi_los_2025.npz: grid 330 x 750, x -18.0..18.7 km, y -17.8..18.5 km from SENG
ascending  LOS: 123,486 pixels, -14.9 to +109.0 mm
descending LOS: 121,976 pixels, -27.4 to +104.0 mm

=== Okada (1985) check against Table 2 (x=2, y=3, d=4, dip=70, L=3, W=2, lambda=mu) ===
strike-slip  ux -8.689e-03 uy -4.298e-03 uz -2.747e-03 | published -8.689e-03 -4.298e-03 -2.747e-03 | max diff 4.2e-07
dip-slip     ux -4.682e-03 uy -3.527e-02 uz -3.564e-02 | published -4.682e-03 -3.527e-02 -3.564e-02 | max diff 2.7e-06
tensile      ux -2.660e-04 uy +1.056e-02 uz +3.214e-03 | published -2.660e-04 +1.056e-02 +3.214e-03 | max diff 4.1e-06
100 m x 100 m sill, 1 m opening, 4 km deep; uz in micrometres at r = 0, 2, 4, 8 km:
  Okada [298.338 170.812  52.758   5.339] | point crack formula [298.416 170.823  52.753   5.338] | Mogi, same volume [149.208 106.764  52.753  13.346]
  above the source a sill lifts the surface 2.00x more than a Mogi source of the same volume change

=== quadtree downsampling ===
ascending  123,486 pixels -> 706 quadtree points (block sizes 4-32 rows)
descending 121,976 pixels -> 643 quadtree points (block sizes 4-32 rows)

=== noise covariance from the Mogi residuals ===
ascending  residual std 4.8 mm; covariance C(r) = 27.4 mm2 * exp(-r / 2.1 km) (sigma 5.2 mm, e-folding distance 2.1 km)
descending residual std 5.2 mm; covariance C(r) = 28.8 mm2 * exp(-r / 4.1 km) (sigma 5.4 mm, e-folding distance 4.1 km)

=== weighted least-squares inversions ===
mogi : x0 (km) -1.19, y0 (km) -1.17, depth (km) 4.28, dV (1e6 m3) 9.33
       centre 63.8716 N, 22.4580 W; volume change 9.33 million m3 in 48 days = 2.25 m3/s; RMS misfit 6.7 mm at the quadtree points; chi2 436 for 1349 points and 10 parameters; 10 evaluations
sill : x0 (km) -1.24, y0 (km) -1.21, depth (km) 4.77, length (km) 5.87, width (km) 0.20, strike (deg) 58.88, opening (m) 6.38
       centre 63.8713 N, 22.4591 W; volume change 7.49 million m3 in 48 days = 1.81 m3/s; RMS misfit 6.0 mm at the quadtree points; chi2 380 for 1349 points and 13 parameters; 118 evaluations

=== Monte Carlo uncertainty (200 realisations of correlated noise) ===
mogi  x0 (km)          -1.19  95% interval    -1.37 ..    -1.00  (std 0.10)
mogi  y0 (km)          -1.17  95% interval    -1.31 ..    -1.00  (std 0.08)
mogi  depth (km)        4.28  95% interval     4.04 ..     4.51  (std 0.12)
mogi  dV (1e6 m3)       9.33  95% interval     8.39 ..    10.55  (std 0.61)
mogi  volume (1e6 m3)     9.33  95% interval     8.39 ..    10.55  -> 2.02-2.54 m3/s
sill  x0 (km)          -1.24  95% interval    -1.43 ..    -1.07  (std 0.09)
sill  y0 (km)          -1.21  95% interval    -1.37 ..    -1.07  (std 0.08)
sill  depth (km)        4.77  95% interval     4.29 ..     5.09  (std 0.20)
sill  length (km)       5.87  95% interval     5.18 ..     6.43  (std 0.33)
sill  width (km)        0.20  95% interval     0.20 ..     2.83  (std 0.78)
sill  strike (deg)     58.88  95% interval    51.70 ..    65.87  (std 3.72)
sill  opening (m)       6.38  95% interval     0.40 ..     7.84  (std 2.28)
sill  volume (1e6 m3)     7.49  95% interval     6.85 ..     8.37  -> 1.65-2.02 m3/s

=== comparison with Parks et al. (2025) ===
mogi : our centre is 0.35 km from the published Mogi location (63.869 N, 22.454 W); depth 4.28 km vs 4.0 km (Mogi) and 4.35-4.7 km (sill)
sill : our centre is 0.35 km from the published Mogi location (63.869 N, 22.454 W); depth 4.77 km vs 4.0 km (Mogi) and 4.35-4.7 km (sill)

figures written to output/