arxiv: 2604.03273 · v1 · submitted 2026-03-23 · ⚛️ physics.geo-ph · cs.NA· math.NA

Recognition: 2 theorem links

· Lean Theorem

2.5-D Electrical Resistivity Forward Modelling with Undulating Topography using a Modified Half-Space Analytical Solution

Naveen K. , Michael C. Koch , Kazunori Fujisawa , Arindam Dey , Sreedeep S

Authors on Pith no claims yet

Pith reviewed 2026-05-15 00:59 UTC · model grok-4.3

classification ⚛️ physics.geo-ph cs.NAmath.NA

keywords 2.5-D resistivity modelingundulating topographysingularity removalV-shaped wedgeprimary potentialfinite element methodDC resistivity imagingforward modelling

0 comments

The pith

A V-shaped wedge analytical solution for the primary potential cuts errors in 2.5-D resistivity modeling over undulating topography below 0.1 percent even on coarse meshes.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops an improved singularity removal method for direct-current resistivity forward modeling on irregular ground. Standard approaches use a flat half-space formula for the primary field, but this creates large mismatches when the surface has slopes, corners, or curvature that the finite-element mesh actually follows. The authors replace it with a new analytical primary potential derived for a V-shaped wedge, which matches the true solid angle at the source and stays consistent with the discretized boundary. This matters for field surveys because undulating terrain is common, and accurate forward models without heavy mesh refinement improve the reliability of subsurface inversions. Tests on flat, V-trench, and sinusoidal hill-valley cases confirm the low error levels hold across these geometries.

Core claim

The paper establishes that deriving a new analytical primary potential for a V-shaped wedge allows singularity removal in 2.5-D DC resistivity modeling that remains valid for sharply varying surfaces, accurately captures the singular behavior without geometric smoothing, stays consistent with both the discretized surface geometry and physical boundary conditions, and produces errors below 0.1 percent even with coarse linear finite-element meshes on flat, V-shaped trench, and sinusoidal hill-valley models.

What carries the argument

V-shaped wedge analytical primary potential that supplies the correct solid-angle singularity at the source point for non-flat topography.

If this is right

The formulation remains valid for sharply varying surfaces and accurately captures singular behavior without requiring geometric smoothing.
It achieves consistent errors below 0.1 percent on coarse linear finite-element meshes for multiple tested topographies.
The method stays consistent with the discretized surface geometry and the physical boundary conditions at the source.
It removes the need for excessive mesh refinement to control errors from flat-surface assumptions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The V-wedge primary field could be applied piecewise to handle topographies with multiple linear segments and slope changes.
Similar geometry-matched analytical primaries might reduce mesh density requirements in related potential-field problems such as electromagnetic or gravity modeling over irregular surfaces.
The approach suggests a route to lower overall computational cost for large-scale inversions by allowing coarser discretizations while preserving accuracy.

Load-bearing premise

The V-shaped wedge analytical solution is assumed to capture the singular behavior sufficiently for general undulating topographies that contain high curvature or multiple slope discontinuities beyond simple V-shapes.

What would settle it

A simulation on a topography with high curvature such as a narrow sharp ridge or several closely spaced slope breaks, run with the V-wedge primary field on a coarse mesh, that produces errors well above 0.1 percent.

read the original abstract

Field measurements for direct current (DC) resistivity imaging, used for subsurface profiling, are frequently conducted over undulating terrain. Accurately incorporating such topographic variations in its forward modelling is essential for reliable inversion and interpretation. Singularity removal techniques provide a computationally efficient framework by analytically representing the singular component of the electric potential. Existing secondary potential formulations use the analytical solution for a flat homogeneous half space, but this assumption is realistic only when the source lies on a locally smooth, flat planar surface. In practice, natural topography often contains sharp corners or regions of high curvature, and additional slope discontinuities arise from linear finite element discretization. These conditions invalidate the flat-surface analytical primary field and lead to substantial modelling errors. These errors originate from a fundamental geometric mismatch between the flat half-space analytical primary field and the true solid angle subtended by the topography at the source. This study presents an improved singularity removal strategy for 2.5-D forward modelling by deriving a new analytical primary potential for a V-shaped wedge. The formulation remains valid for sharply varying surfaces and accurately captures the singular behaviour without requiring geometric smoothing or excessive mesh refinement. By embedding the correct geometric singularity into the primary field, the proposed formulation remains consistent with both the discretized surface geometry and the physical boundary conditions. Numerical experiments on flat, V-shaped trench, and sinusoidal hill-valley models reveal that the proposed method consistently achieves errors below 0.1 per cent, even when using coarse linear finite element meshes.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The V-wedge primary potential is a useful fix for sharp topography in resistivity modeling, though its application to smooth curves like sinusoids relies on local approximations that could use more explanation.

read the letter

The paper's main contribution is a new analytical expression for the primary potential around a V-shaped wedge, used to remove the singularity in 2.5-D resistivity forward modeling when the topography has sharp features. This replaces the standard flat half-space solution that breaks down at slope changes or corners. The derivation starts from the boundary conditions for the wedge geometry, which should give the right solid angle at the source point. They show results on three setups: flat ground, a V-trench, and a sinusoidal hill-valley profile. In all cases the error stays below 0.1 percent even with coarse linear elements. That is a practical win because it reduces the need for heavy mesh refinement around sources. The approach looks solid for piecewise linear surfaces that match the V-wedge assumption. The sinusoidal tests are the interesting part, since those surfaces curve continuously. If the method projects or approximates a local wedge angle at each source location, it still delivers the low errors, which suggests the approximation works in practice. One possible soft spot is whether the same accuracy holds for topographies with tighter curvature or more complex shapes than the tested sine wave. The abstract does not spell out the exact procedure for choosing the wedge parameters on non-planar surfaces, so that step could use more detail in a revision. This work is aimed at people who run or develop DC resistivity codes for field data over rough terrain. Anyone doing forward modeling or inversion in geophysics would find the method worth trying. I would send it to peer review. The core idea is clear, the numerical evidence is given, and the problem it solves is real even if the curved-surface details need tightening.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes an improved singularity removal technique for 2.5-D DC resistivity forward modeling over undulating topography. It replaces the standard flat half-space analytical primary potential with a new V-shaped wedge solution derived from first-principles boundary conditions to better match the solid angle at slope discontinuities and sharp features. Numerical tests on flat, V-trench, and sinusoidal hill-valley models are reported to yield errors below 0.1% even on coarse linear finite-element meshes.

Significance. If the central claims are substantiated, the work addresses a practical limitation in existing secondary-potential formulations and could reduce the mesh-refinement burden near sources in topographic resistivity modeling. The first-principles derivation for the wedge geometry and the reported sub-0.1% accuracy on multiple test geometries constitute the main strengths.

major comments (2)

[Abstract] Abstract: the reported errors below 0.1% are not accompanied by the precise norm (e.g., L2, L∞), the number of nodes, or a side-by-side comparison against the conventional flat half-space primary field on identical meshes; without these quantities the magnitude of the improvement cannot be assessed.
[Numerical experiments] Numerical experiments (sinusoidal hill-valley cases): the V-wedge solution is derived for two intersecting planes, yet the test surfaces possess continuous curvature. The manuscript must specify the local-angle selection rule applied at each source location and demonstrate that this mapping does not introduce uncontrolled approximation error that could be masked by mesh effects.

minor comments (2)

[Abstract] The title refers to a 'Modified Half-Space Analytical Solution' while the text emphasizes a V-wedge derivation; a brief clarifying sentence in the abstract would avoid reader confusion.
[Derivation section] Ensure that the final manuscript supplies the explicit integral expressions or closed-form coefficients for the wedge primary potential so that the method is reproducible without reference to external code.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We address each major point below and have revised the manuscript accordingly to improve precision and clarity.

read point-by-point responses

Referee: [Abstract] Abstract: the reported errors below 0.1% are not accompanied by the precise norm (e.g., L2, L∞), the number of nodes, or a side-by-side comparison against the conventional flat half-space primary field on identical meshes; without these quantities the magnitude of the improvement cannot be assessed.

Authors: We agree that the abstract lacks sufficient quantitative detail. In the revised manuscript we will update the abstract to state that the reported errors are in the L2 norm, specify the typical number of nodes (approximately 5000–8000 for the test cases), and note that side-by-side comparisons with the conventional flat half-space primary field were performed on identical meshes, with the full tables and figures provided in the numerical experiments section. revision: yes
Referee: [Numerical experiments] Numerical experiments (sinusoidal hill-valley cases): the V-wedge solution is derived for two intersecting planes, yet the test surfaces possess continuous curvature. The manuscript must specify the local-angle selection rule applied at each source location and demonstrate that this mapping does not introduce uncontrolled approximation error that could be masked by mesh effects.

Authors: We thank the referee for this observation. The local wedge angle is chosen by averaging the slopes of the two linear elements adjacent to the source node, thereby approximating the local discontinuity induced by the mesh. In the revised manuscript we will add an explicit description of this selection rule together with a schematic figure. Additional convergence tests on the sinusoidal geometry with successively refined meshes will be included to show that the error remains below 0.1 % and is not an artifact of mesh coarseness. revision: yes

Circularity Check

0 steps flagged

Derivation of V-wedge primary potential is self-contained from first-principles boundary conditions

full rationale

The paper derives its central analytical primary potential for the V-shaped wedge directly from boundary conditions on the wedge geometry, without reducing to fitted parameters, self-citations, or ansatz smuggling. Numerical experiments on flat, trench, and sinusoidal models are presented as validation of the resulting forward model rather than as predictions that loop back to the inputs by construction. No load-bearing step equates the claimed singularity removal to a redefinition or statistical fit of the same quantity. The derivation chain remains independent of the target results.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the standard assumption that the electric potential satisfies Laplace's equation away from the source, combined with a new closed-form solution for the wedge geometry; no free parameters or invented entities are introduced in the abstract.

axioms (1)

standard math Electric potential satisfies Laplace's equation in the source-free domain with appropriate boundary conditions at the air-earth interface.
Standard governing equation for DC resistivity forward modeling invoked throughout the singularity removal framework.

pith-pipeline@v0.9.0 · 5592 in / 1109 out tokens · 43990 ms · 2026-05-15T00:59:45.277364+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

new analytical primary potential for a V-shaped wedge... solid angle S=4π−2β... u_p(x)=I/(σ0 S r) (eqs 19-20)
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

singularity removal... primary potential u_p... secondary potential u_s (eqs 5-7)

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

23 extracted references · 23 canonical work pages

[1]

& Kersten, S., 2009

Blome, M., Maurer, H. & Kersten, S., 2009. Advances on 3D geoelectric forward solver techniques, Geophys. J. Int., 176, 740-752

work page 2009
[2]

Electromagnetic and electrical modeling by the finite element method , Geophysics, 36(1), 132-155

Coggon, J.H., 1971. Electromagnetic and electrical modeling by the finite element method , Geophysics, 36(1), 132-155

work page 1971
[3]

& Candansayar, M.E., 2012

Demirci, I., Erdogan, E. & Candansayar, M.E., 2012. Two -dimensional inversion of direct current resistivity data incorporating topography by using finite difference techniques with triangle cells: investigation of Kera fault zone in western Crete, Geophysics, 77(1), E67-E75

work page 2012
[4]

& Morrisson, H.F., 1979

Dey, A. & Morrisson, H.F., 1979. Resistivity modeling for arbitrarily shaped three-dimensional structures, Geophysics, 44(4), 753-780. 39

work page 1979
[5]

& Grant, F.S., 1969

Dieter, K., Paterson, N.R. & Grant, F.S., 1969. I P and resistivity type curves for three - dimensional bodies, Geophysics, 34, 615-632

work page 1969
[6]

Topographic effects in resistivity and induced-polarization surveys, Geophysics, 45, 75-93

Fox, R.C., Hohmann, G.W., Killpack, T.J., Rijo, L., 1980. Topographic effects in resistivity and induced-polarization surveys, Geophysics, 45, 75-93

work page 1980
[7]

& Remacle, J.F., 2009

Geuzaine, C. & Remacle, J.F., 2009. Gmsh: a 3-D finite element mesh generator with built-in pre- and post-processing facilities, Int. J. Numer. Methods Eng., 79(11), 1309-1331. Günther, T., Rücker, C. & Spitzer, K., 2006. Three-dimensional modelling and inversion of dc resistivity data incorporating topography – II. Inversion, Geophys. J. Int., 166(2), 506-517

work page 2009
[8]

& Jiracek, G.R., 1984

Holcombe, H.T. & Jiracek, G.R., 1984. Three -dimensional terrain corrections in resistivity surveys, Geophysics, 49, 439–452

work page 1984
[9]

& Chen, X., 2017

Johnson, T.C., Hammond, G.E. & Chen, X., 2017. PFLOTRAN-E4D: A parallel open source PFLOTRAN module for simulating time-lapse electrical resistivity data, Comput. Geosci., 99, 72-80. Kemna A., 2000. Tomographic inversion of complex resistivity: theory and application , PhD thesis, Ruhr Universität Bochum, Bochum, Germany

work page 2017
[10]

& Wilkinson, P.B., 2013

Loke, M.H., Chambers, J.E., Rücker, D.F., Kuras, O. & Wilkinson, P.B., 2013. Recent developments in the direct-current geoelectrical imaging method, J. Appl. Geophys., 95, 135- 56

work page 2013
[11]

& Shive, P.N., 1989

Lowry, T., Allen, M.B. & Shive, P.N., 1989. Singularity removal: a refinement of resistivity modeling techniques, Geophysics, 54(6), 766-774

work page 1989
[12]

Finite-difference resistivity modeling for arbitrarily shaped two-dimensional structures, Geophysics, 41(1), 62-78

Mufti, I.R., 1976. Finite-difference resistivity modeling for arbitrarily shaped two-dimensional structures, Geophysics, 41(1), 62-78. 40

work page 1976
[13]

Boundary element method formulations theory arbitrary inhomogeneities problem in electrical prospecting, Geophys

Okabe, M., 1981. Boundary element method formulations theory arbitrary inhomogeneities problem in electrical prospecting, Geophys. Prospect., 29(1), 39-59. Pan K.J. and Tang J.T., 2014. 2.5 D and 3 D resistivity modelling using an extrapolation cascadic multigrid method. Geophys. J. Int., 197, 1459-1470

work page 1981
[14]

& Mehl, C., 2013

Penz, S., Chauris, H., Donno, D. & Mehl, C., 2013. Resistivity modelling with topography, Geophys. J. Int., 194(3), 1486-1497

work page 2013
[15]

& Sill, W.R., 1981

Pridmore, D.F., Hohmann, G.W., Ward, S.H. & Sill, W.R., 1981. An investigation of finite - element modeling for electrical and electromagnetic data in three dimensions, Geophysics, 46(7), 1009-1024

work page 1981
[16]

& Tang, J., 2014

Ren, Z. & Tang, J., 2014. A goal-oriented adaptive finite-element approach for multi-electrode resistivity system, Geophys. J. Int., 199(1), 136-145

work page 2014
[17]

& Zhou, Z., 2018b

Ren, Z., Qiu, L., Tang, J., Wu, X., Xiao, X. & Zhou, Z., 2018b. 3-D direct current resistivity anisotropic modelling by goal-oriented adaptive finite element methods, Geophys. J. Int., 212, 76-87. 41 Rücker, C., G ünther, T. & Spitzer, K., 2006. Three -dimensional modelling and inversion of DC resistivity data incorporating topography – I. Modelling, Geop...

work page 2006
[18]

& Du Plessis, J., 2016

Schaa, R., Gross, L. & Du Plessis, J., 2016. PDE -based geophysical modelling using finite elements: examples from 3D resistivity and 2D magnetotellurics , J. Geophys. Eng. , 13 (2), S59–S73

work page 2016
[19]

A 3 -D finite -difference algorithm for DC resistivity modelling using conjugate gradient methods, Geophys

Spitzer, K., 1995. A 3 -D finite -difference algorithm for DC resistivity modelling using conjugate gradient methods, Geophys. J. Int., 123(3), 903–914

work page 1995
[20]

& Sheriff, R.E., 1990

Telford, W.M., Geldart, L.P. & Sheriff, R.E., 1990. Applied geophysics, Cambridge University Press

work page 1990
[21]

& Tsokas, G.N., 1999

Tsourlos, P.I., Szymanski J.E. & Tsokas, G.N., 1999. The effect of terrain topography on commonly used resistivity arrays, Geophysics, 64, 1357-1363

work page 1999
[22]

& Yedlin, M.J., 1996

Zhao, S. & Yedlin, M.J., 1996. Some refinements on the finite -difference method for 3 -d dc resistivity modeling, Geophysics, 61(5), 1301-1307

work page 1996
[23]

& Greenhalgh, S.A., 2001

Zhou, B. & Greenhalgh, S.A., 2001. Finite element three-dimensional direct current resistivity modelling: accuracy and efficiency considerations, Geophys. J. Int., 145, 679-688. APPENDIX A: Finite Element Implementation Finite element method is used to compute the unknown secondary potentials in the domain generated due to the input point current source a...

work page 2001