Computing Solar Eclipses — Research

NASA SVS and Ernie Wright

workingupdated 2026-09-15svswrightrasterlunar-limblolashapefiles
  • The method is published. Wright, E. and Young, C. A. 2024, "A Raster-oriented Method for Creating Eclipse Maps", The Astronomical Journal 168:163, published 2024 September 19, DOI 10.3847/1538-3881/ad6b23, open access under CC BY 4.0 1.
  • The core idea is to test every map pixel, at every time step, against a lunar limb profile built from the SLDEM2015 and LOLA LDEM elevation models, with the observer lifted to the SRTM terrain height above the WGS 84 ellipsoid via the EGM96 geoid 1.
  • The datasets are ESRI shapefiles and KML in WGS 84 longitude and latitude: umbra polygons at 1 s and 10 s steps, path outlines, central line, duration contours at 30 s and obscuration contours at 1 per cent and 5 per cent 2.
  • The constants stated by SVS are DE421, Earth radius 6378.137 km, WGS 84 flattening, Moon datum radius 1737.4 km, Sun radius 696,000 km (959.645 arcseconds at 1 au) and, for 2017, ΔT\Delta T = 68.917 s 3.
  • The differences from other predictors come from the solar radius (959.63 versus 959.95 arcseconds), the limb, and terrain, not from the ephemeris. Espenak's own site says limb corrections move limits by 1 to 3 km 4.

What this topic covers

This topic covers the one modern eclipse-mapping method that is both published in full and backed by released data. It takes the 2024 Astronomical Journal paper formula by formula, inventories every SVS eclipse page and the files it releases down to the shapefile schema a developer will load, and sets the SVS path against the NASA eclipse site, EclipseWise, Jubier's maps and the Besselian Elements map to say which constant explains each disagreement.

Notes in this topic:

What this topic changes for the pipeline

The limb correction stops being a post-processing step on Besselian contacts and becomes a per-pixel test, so one code path produces the path limits, central line, duration contours and obscuration contours. The released umbra polygons and cities file become the pipeline's ground truth, fixing the target constants, resolutions and outputs. The limb profile, not the ephemeris, sets the resolution floor at about 600 m.

No public code exists, and no repository by Wright or SVS was found. Every SVS page and the paper PDF were read from Wayback Machine snapshots, because svs.gsfc.nasa.gov refused connections and iopscience.iop.org sat behind a bot wall at the time of writing (2026 September). The snapshot date is recorded in the note field of each row of refs.tsv.

References

  1. 1peer-reviewed Wright, E. and Young, C. A. 2024, A Raster-oriented Method for Creating Eclipse Maps, AJ 168, 163 The paper. Read in full from the Wayback capture of the IOP PDF dated 2024-11-19 (iopscience.iop.org/article/10.3847/1538-3881/ad6b23/pdf), extracted with pdftotext. CC BY 4.0. Sections 4 and 5 give the limb-profile construction and the limb test; Appendix A gives the SPICE calls.
  2. 2primary NASA SVS 5073, The 2023 and 2024 Solar Eclipses: Map and Data (2023) The data page: Making the Map statement (SRTM, LRO, DE421), layer descriptions for both zips, KML zips, cities JSON schema. Read from Wayback capture 2024-04-28.
  3. 3primary NASA SVS 4515, 2017 Path of Totality (2016) First limb- and terrain-corrected product. Constants table with WGS84, EGM96, 1737.4 km, 696,000 km (959.645 arcsec), DE421, EOP kernel, Delta T 68.917 s. Read from Wayback capture 2026-08-28.
  4. 4primary Espenak, F., Path of Total Solar Eclipse of 2024 Apr 08, NASA Eclipse Web Site VSOP87/ELP2000-85, Delta T 70.6 s, the 1 to 3 km and 1 to 3 s limb caveat, central line table, greatest eclipse 18:17:18.3 UT. Read via fetch tool.

In this section