The cone-volume model projects 8.9 reference planets reached so far.

Central projection at the 2026 cutoff: 8.9 broader reference planets or 1.1 small habitable-zone planets. Already-emitted signals reach 314 reference planets by 2126.

Continuous modeled DSN uplinks and 70-m beam widths dominate the volume. Shading excludes these larger uncertainties. This is not a complete transmitter census.

The first curve combines all included emitted cones through time, counting overlapping directions once. The solid curve extrapolates Arecibo/Goldstone asteroid-radar history from annual detected-target activity; the dashed curve uses recovered radar geometry. Both retain the same five spacecraft tracks and selected intentional messages. No transmissions after the cutoff are assumed. Future growth is caused by light already in flight.

The broader population is calibrated at 1.6 planets per ordinary star, limited to 5 Earth masses–10 Jupiter masses and 0.5–10 AU. The source’s quoted 0.71–2.32 occurrence endpoints give 3.9–12.9 reference planets today, before beam-history and host-population transfer uncertainty. These are sensitivity endpoints for this extrapolation, not a confidence interval for local contacts. The population is not an exhaustive count of all planets. The HZ estimate uses a separate small-planet occurrence model; the populations overlap and must not be added. Fractional values are expectations over unknown planetary systems.

Already-emitted cones keep reaching more planets.

Central estimate: 2,305 ly³ reached so far. The next year's propagation adds 0.58 reference planets or 0.072 small HZ planets in expectation.

Year (Oct 7) Reference planets Small HZ planets Volume (ly³) Mean /2′ (ly³)
2026 8.88 1.09 2,305 0.00665
2036 16.02 1.97 4,156 0.012
2126 314.20 38.60 81,543 0.235
2326 4,542.79 558.11 1,178,965 3.4

These are population expectations; catalogue encounters occur intermittently. Numerical precision supports reproducibility and is not observational certainty.

“Average cone volume” is standardized as the combined reached volume divided by unique angular area expressed in full-width 2′ beam-disk equivalents. The current union is 346,457 such equivalents, averaging 0.0067 ly³ each. These are units of area, not a count of independent transmitting sessions. Wider beams correspond to many units; a moving beam creates a corridor; revisits add no new direction.

Because the beams have different ages, multiplying every direction by one average age would bias the result. Instead each direction keeps its earliest emission date. A short pulse has passed through the interior of its cone even though its photons now occupy a moving shell.

The radar model places first coverage uniformly within an empirical ±49.4° ecliptic band, covering 95% of the recovered Arecibo angular union. The annual coverage hazard uses 25.9 reference footprints per Arecibo detected target apparition and 81.4 per Goldstone apparition. Repeated radar coverage and overlap with spacecraft/message cones are subtracted. Unknown interior years are explicitly interpolated; activity before the first record is omitted. Other radar classes remain uncalibrated.

Historical hardware also changed: some Goldstone tracks used 13-cm radiation rather than 3.5 cm; Arecibo used 70 cm before its upgrade; the 1968 Icarus transmitter was 26 m. Modern area-per-target calibration cannot certify these earlier footprints. The 0.5×–2× scenarios test transfer sensitivity and do not guarantee bounds on the true total. Intentional Altair’s modeled 54′ width is an additional assumption (half/double-width alternatives), and its A7 star is outside the G/M HZ calibration.

Old signals add about 0.58 reference planets over the next year.

The small-HZ projection adds about 0.072 planets in the same year. These arrivals come from signals already in flight.

Propagation interval from cutoff New reference planets New small HZ planets
Next day 0.00157 0.000192
Next week 0.011 0.00135
Next year 0.584 0.0718

These smooth density-based increments are expected values, not scheduled discrete contacts. Catalogued stars arrive intermittently. No future transmissions are included.

New sky coverage can be extrapolated in a common beam-area unit.

The 2025 Goldstone activity benchmark extrapolates about 3,336 reference footprints/year, or 9.1/day, before overlap with all other source families. This is a historical activity benchmark; current operations require their own logs.

Estimate and scope New 2′ area equivalents/year Average/day
Arecibo recovered historical mean, 2002–2019 1,382 3.8
Goldstone 2025 activity-transfer budget 3,336 9.1
Goldstone continuation scenario, new versus cutoff coverage 3,304 9.0
Continuing DSN scenario, Oct2026–Oct2027, new versus past DSN 265 0.73

The rows describe different periods and baselines and must not be added. Arecibo stopped in 2020; DSS-14 is offline for modernization during 2026–2028. No defensible observed global current “cones/day” total follows from these public records.

Each new 2′ area equivalent is about 0.00087 square degrees. A useful forecast of volume is obtained by retaining the dates of the new coverage, then propagating its front; it is not a fixed volume added per transmission. The table above reports the resulting average volume through time for the cones already emitted. A continuing-transmission scenario would require a future activity schedule, which is distinct from this primary already-emitted forecast.

Continuing radar activity could add about 3,304 new reference footprints per year.

Hypothetical Goldstone activity after a May 2028 restart, at its 2025 level, initially adds about 9.0 new 2′ area equivalents/day. After 100 years the new cones alone average 0.022 ly³ per unique reference footprint.

Years after assumed restart Cumulative new 2′ equivalents New equivalents/year then Mean reached volume per new footprint (ly³) Extra reference / HZ planets
1 3,304 3,304 0 0 / 0
10 33,023 3,301 1.76e-05 0.00223 / 0.000274
50 164,809 3,289 0.00277 1.76 / 0.216
100 328,855 3,273 0.0222 28.1 / 3.46

This is a separate future-activity scenario, not an observed operating forecast. It holds the central historical model fixed and varies future area transfer by 0.5×–2×. The worldwide current rate remains unknown.

Only sky unilluminated in the cutoff model can become a new direction here. Repeated future coverage is removed with the same uniform-band hazard, and the remaining new-direction rate slowly falls as available sky is used. The assumed restart date follows the scheduled modernization interval; actual activity and dates may differ. No other future transmitter is introduced.

Newly emitted cones shorter than4.24 ly cannot reach another known stellar system, so their planet projection remains zero initially. Their average cumulative volume grows with time because their pulses continue crossing new space. These additional volumes occupy new angular directions in the adopted model and may be combined with the old-emission projection only at matching dates and with the same conditional assumptions. They are not included in the first curve.

Reached means physical illumination under the declared model.

“Reached” means a signal’s leading front crosses a target; it does not imply detection, communication, life, or confirmed habitability. A small planet in a modeled habitable zone (HZ) receives a suitable range of stellar energy; its atmosphere and surface conditions are unknown. In the catalogue estimator the spacecraft model predicts 4 arrived objects, and the selected intentional records add distinct Altair. Two eligible red dwarfs imply 0.32 expected small HZ planets; this conditions on known hosts and cannot be added to the volume-density estimate.

The source inventory includes five spacecraft directions in a published Deep Space Network (DSN) study, public JPL Goldstone configuration records, recovered Arecibo asteroid transmit-run pages, and selected primary-documented intentional interstellar messages. The selected intentional records add one distinct named target, Altair, with nominal arrival around 2000. Ordinary leakage, military/aircraft radars, sidelobes, lunar/planetary radar, and other message campaigns require additional histories.

Spacecraft beams reach new objects intermittently.

Next nominal arrival: around 2028 via Voyager 2. The selected model adds 129 objects over the next 100 years.

Forecast interval End date (UTC) Cumulative objects Newly reached objects
Next day 2026-10-08 4 0
Next week 2026-10-14 4 0
Next calendar year 2027-10-07 4 0
Next 10 years 2036-10-07 7 3
Next 100 years 2126-10-07 133 129

These are conditional geometric-model forecasts, with uncertain dates. No new arrivals in the next year is specific to this inventory. A 100-year average is not a constant daily arrival rate.

The 100-year mean is 1.29 catalogue objects/year, 0.025/week, or 0.0035/day. Actual arrivals are discrete; all three immediate intervals above contain zero arrivals. Unique Gaia objects are not necessarily independent stellar systems.

Selected model: Nominal 70-m model: S-band for four probes, X-band for New Horizons; Gaia light-arrival correction, moving-star interception, continuous modeled uplinks. It assumes continuous uplinks during the adopted mission windows and omits the known 2020 Voyager 2 uplink gap. It reconstructs nominal geometry, rather than proving an uplink occurred during each alignment. All forecast rows first filter emission epochs to the past. Fresh six-hour JPL ephemerides and 15-minute refinements replace the original coarse track; spherical geometry and stellar proper-motion conventions are checked. Remaining limitations include incoming geocentric spacecraft directions instead of ground-site outgoing pointings, missing radial velocities, beam power contours, and launch-time continuous-coverage assumptions. The stellar epoch accounts for the incoming catalogue light delay and the outgoing wavefront flight; the final version solves the moving-star interception rather than keeping distance fixed.

The earlier chat’s 1,296 figure counted star–mission pairs. The author catalogue has 1,278 unique Gaia identifiers, of which 1,277 have past modeled first-emission epochs at this cutoff. Its calendar dates used rounded arrival years. Those preserved paper-reproduction files are separate from the selected reanalysis. The selected reanalysis uses New Horizons’ X-band uplink width of 0.038° for a 70-m antenna; the original paper used the 0.128° S-band width for every mission.

The Altair transmission has a nominal arrival around 2000.

Selected primary records identify 17 nearby intentional stellar targets. Altair is the one whose nominal intended-target arrival predates the cutoff; none of these targets has a new arrival in the next day, week, or year.

Selected named target Earliest recorded date Nominal intended-target arrival
Altair / α Aql 1983-08-15 ~2000
HIP 74995 / Gliese 581 2008-10-09 ~2029
GJ 273 / Luyten’s Star 2017-10-16 ~2030
HIP 4872 2003-07-06 ~2036

These are intended-target travel-time estimates from date-only records and rounded distances, not reconstructed beam intersections or observed reception. The selected table is not an exhaustive intentional-message history.

The 1983 record describes a one-hour Stanford 46-m transmission aimed at Altair, approximately 16.7 ly away. Date-plus-distance arithmetic puts its nominal arrival around 2000. Altair is distinct from the four current DSN-model objects, giving five named/modelled objects in the combined current inventory under their respective assumptions. The reconstructed Arecibo archive adds no current nearby-catalogue object. This bounded count cannot establish the worldwide total.

The saved intentional table contains 21 date-target sessions, 17 distinct nearby stars and the additional M13 globular-cluster target. Source-rounded arithmetic predicts three additional selected stellar targets over ten years and 16 over 100 years. These should not be added to future DSN or radar catalogue totals without identifier cross-matching. Altair is an A7V star; the Sun-like G-star and M-dwarf HZ occurrence calibrations do not apply to it.

Sónar’s actual 2017 transmissions were on October 16–18, preceding the November announcement. The source’s 12.4-ly distance places the first nominal arrival around 2030. M13’s approximately 25,000-ly travel time places the intended-cluster arrival around year 26,974; foreground beam encounters were not reconstructed here.

Named stellar target Earliest transmission date Rounded source distance (ly) Nominal arrival year
alpha Aql 1983-08-15 16.7 ~2000
HD 186408 1999-05-24 70.5 ~2069
HD 190406 1999-06-30 57.6 ~2057
HD 178428 1999-06-30 68.3 ~2067
HD 190360 1999-07-01 51.8 ~2051
HD 197076 2001-08-29 68.5 ~2070
HD 95128 2001-09-03 45.9 ~2047
HD 50692 2001-09-03 56.3 ~2057
HD 126053 2001-09-03 57.4 ~2059
HD 193664 2001-09-04 57.4 ~2059
HD 76151 2001-09-04 55.7 ~2057
HIP 4872 2003-07-06 32.8 ~2036
HD 245409 2003-07-06 37.1 ~2040
HD 75732 2003-07-06 40.9 ~2044
HD 10307 2003-07-06 41.2 ~2044
HIP 74995 2008-10-09 20.3 ~2029
GJ 273 2017-10-16 12.4 ~2030

Host stars need different planet-occurrence estimates.

Conservative estimates give 0.38–0.63 small HZ planets per Sun-like G star and about 0.16 per early M dwarf; these are mean planet counts with broad uncertainty.

Host / occurrence model Mean planets per star Published uncertainty Planet radius
G star, lower completeness model 0.38 68% credible interval 0.16–0.88 0.5–1.5 Earth radii
G star, upper completeness model 0.63 68% credible interval 0.25–1.57 0.5–1.5 Earth radii
Early M dwarf, conservative HZ 0.16 Approximate 68% interval 0.09–0.33 1.0–1.5 Earth radii

Different radius definitions prevent interpreting the mixed population as one uniform rocky-planet selection. An occurrence rate is not the fraction of stars that host a planet, nor the fraction with inhabited planets.

The four reached objects in the original model are GJ 1276 (white dwarf), GJ 1154 (M dwarf), 2MASSW J1507476-162738 (brown dwarf), and GJ 359 (M dwarf). Applying the M-dwarf calibration to the two eligible hosts gives 2 × 0.16 = 0.32 expected conservative-HZ planets; the optimistic-HZ model gives 0.48. White and brown dwarfs are excluded from these occurrence calibrations, rather than asserted incapable of planets. This remains conditional on the modeled exposure and on applying the early-M result to these hosts.

Dressing & Charbonneau’s selected cool-dwarf sample spans 2,661–3,999 K, with median 3,746 K and mostly early M dwarfs. Transferring its population mean to these M4/M4.5 hosts adds an unquantified population assumption; 0.32 is an expected planet count under that transfer, not a measurement of their planets.

Probe direction Gaia identifier Distance (pc) Nominal model arrival
Pioneer 10 2611561706216413696 8.54 2002-02-06
Voyager 2 3698534434669937024 8.09 2007-03-22
Voyager 2 6306068659857135232 7.41 2007-10-11
Pioneer 11 640862653425455360 12.34 2018-07-13

The local 10-pc census yields 0.0836 ordinary stars/pc³ and 0.0609 systems containing ordinary stars/pc³. M-dwarf and F/G/K counts give a central conservative small-HZ-planet density of 0.0164 planets/pc³, equivalent to about 0.20 per ordinary star. The conservative upper median gives 0.0199/pc³. An endpoint sensitivity spans 0.0085–0.0451/pc³; this combines marginal source endpoints and is not a joint credible interval. Binary suppression, local population transfer, and incomplete classifications are unquantified limitations.

Radar runs often revisit the same sky.

The retrieved June–September 2025 Goldstone subset contains 4,323 runs and approximately 1,413 beam-disk equivalents of unique nominal angular coverage.

This is new coverage within the retrieved subset, starting from an empty inventory. Historical beams may already have covered these directions. Reception/configuration windows give a nominal tracking envelope, rather than every actual pulse.

A “new direction” is defined here by the incremental union of circular half-power beam footprints. Its area divided by one 1.92-arcminute beam disk is an equivalent number of beam disks, which may be fractional. It is not a count of sessions, targets, or disjoint cones. Tracking can sweep many beam widths; repeated pointings can add no new area.

The 123 configuration windows cover 16 targets, 19 observing dates, and 30 target-days. Their total span is 71.1 hours of receive/configuration windows, not 71.1 transmitter-on hours. The angular union is 0.000346076 steradians, or 0.0028% of the sky. Monte Carlo integration standard error is 5.9e-07 sr; unresolved pulse gaps and historical pointing offsets produce additional source uncertainty.

Dividing by the 105-day calendar interval gives 13.5 equivalent new beam disks/day, or 4,914/year if that sparse sample were repeated at the same rate. This is an illustrative sample-rate scaling, not measured annual global activity. Nine windows with setup/pointing warnings remain in the nominal envelope; their exclusion is provided in the numerical sensitivity results.

Recovered Arecibo beams have no current catalogue arrivals.

The recovered archive contains 74,201 usable transmit intervals across 762 targets. Its nominal nearby-catalogue match has 0 arrivals today and 0 in the next year.

Forecast interval New nearby catalogue objects Expected new ordinary stars, volume model Expected new small HZ planets, volume model
Next day 0 9.7e-06 1.9e-06
Next week 0 6.8e-05 1.3e-05
Next calendar year 0 0.0038 0.00074

Catalogue objects and homogeneous-density expectations are different estimators. Zero catalogue matches is conditional on nominal pointing reconstruction, the chosen beam contour, and catalogue completeness.

Systematic retrieval recovered 75,219 raw rows from 834 advertised per-object run pages. Removing 486 identical 22-field measurement duplicates leaves 74,733 distinct measurements, of which 74,201 have verified transmitter start/end times. The other 532 rows remain in the audit. Usable intervals span August 2001–June 2020, all at 2.38 GHz. Their summed duration is 975.5 recorded window-hours; removing UTC overlaps gives 973.1 hours of window coverage. Overlapping archive representations are not counted as independent transmitter time. This retrieves 88.9% of the advertised run count; no contact estimate is scaled by that percentage.

Ground-site outgoing directions include light-time point-ahead and sample recorded transmit intervals at no more than 15″ spacing for the nominal 60″ radius. UCLA’s date is anchored at receiver start; 72 midnight-crossing transmit windows were shifted to the previous day using that documented convention. Current JPL orbit solutions reconstruct the intended tracking; original antenna offsets and power logs are unavailable. Residual simultaneous cross-target records involve 92 interval keys across 13 targets; the nominal archive envelope retains them and an exclusion sensitivity is saved separately. This cannot certify both incompatible pointings occurred. The union contains approximately 26,432 two-arcminute beam-disk equivalents. Its current cumulative volume is 8.56 ly³, giving 0.0206 expected ordinary stars and 0.00405 expected small HZ planets under the local homogeneous-density model. These are spatial expectations, rather than fractional catalogue observations.

The catalogue calculation uses 300,526 GCNS objects with median distance at most 100 pc. It predicts 0 new objects over ten years and 8 over 100 years. Eventually 168 nominal catalogue sources intersect the sampled beams; expanding the radius by 15″ gives 215 in a grazing sensitivity. That expansion diagnoses sampling at edges and is not a certified uncertainty interval. White/brown dwarfs and components are included in catalogue-object counts, so multiplying that total by an ordinary-star occurrence rate would be invalid. The density model independently gives 8.29 ordinary stars and 1.63 small HZ planets by 2126.

Historical new sky can be measured from dated radar runs.

Within the recovered Arecibo inventory, complete calendar years 2002–2019 averaged 1,382 new two-arcminute beam-disk equivalents/year, or 3.8/day.

This is a historical subset rate, starting from an empty 2001 inventory. It does not establish globally untouched sky, full-archive activity, or today's operational rate.

Calendar year New 2′ beam-disk equivalents Average/calendar day
2001 390 1.1
2002 822 2.3
2003 1,148 3.1
2004 980 2.7
2005 1,039 2.8
2006 550 1.5
2007 441 1.2
2008 1,038 2.8
2009 1,191 3.3
2010 610 1.7
2011 398 1.1
2012 1,737 4.7
2013 1,593 4.4
2014 1,981 5.4
2015 2,221 6.1
2016 1,326 3.6
2017 2,121 5.8
2018 2,188 6.0
2019 3,499 9.6
2020 1,159 3.2

“New” means sky area first covered within this recovered dated inventory. A continuous tracking sweep creates a ribbon of angular coverage, so the number of cones depends on the declared beam footprint convention. The annual calculation removes angular overlap across years. Individual run counts, target counts, and beam-disk equivalents answer different questions.

Recent radar fronts cover little interstellar volume.

The 2025 Goldstone sample reaches no other stellar system by today or the next year: its front is still inside the nearest-star distance. The two recovered Arecibo apparitions have zero current catalogue matches.

Recovered radar inventory Current catalogue arrivals Nominal future result
Goldstone, 2025 subset 0; front shorter than 1.4 ly Density-model expectation ≈0.29 stars by 2126
Arecibo, DP107 2008/2016 0 One nearby-catalogue star nominally around 2156

Arecibo is a selected asteroid case, not the complete archive. A density expectation can be fractional; the actual nearby catalogue is discrete and spatially clustered.

The Arecibo case has 449 actual transmit windows, 9.04 transmitter-on hours, and nominal outgoing directions reconstructed with JPL Horizons. Four 2008 rows lacking transmitter times remain in the raw data and are explicitly excluded from the geometric reconstruction. The combined nominal sky union is approximately 131 two-arcminute beam disks. The current homogeneous-density expectation is only 0.00015 ordinary stars and 2.9e-05 small HZ planets; the catalogue finds zero. These values are an expectation under spatial randomization, rather than a claim of fractional observed stars.

Forecast tables provide cumulative ordinary stars, stellar systems, small-HZ planets, G-star-only HZ planets, and increments over the next day, week and calendar year. For example, at 2126 the 2025 Goldstone envelope gives 0.0086 expected new ordinary stars in its following year, 0.00016/week, or 2.3e-05/day; expected HZ planets scale by the defined population mixture. This smooth rate describes a density model, while the DSN table forecasts discrete catalogue arrivals.

Beam age and angular coverage determine reach.

Even 100,000 distinct fixed two-arcminute beams aged 60 years reach only about 4.6 ordinary stars and 0.9 small HZ planets in a homogeneous local-density model.

The number and ages of real independent directions are unknown. These are explicit sensitivity scenarios, not a fitted global estimate or an upper bound on moving beams.

Assumed distinct fixed beams Signal age (years) Ordinary stars expected Small HZ planets expected New ordinary stars in next year
1,000 20 0.00169 0.000332 0.000269
1,000 40 0.0136 0.00268 0.00105
1,000 60 0.0461 0.00906 0.00234
10,000 20 0.0169 0.00332 0.00269
10,000 40 0.136 0.0268 0.0105
10,000 60 0.461 0.0906 0.0234
100,000 20 0.169 0.0332 0.0269
100,000 40 1.36 0.268 0.105
100,000 60 4.61 0.906 0.234

The full 84,054-run Arecibo count cannot be substituted for the first column. Each run has a different time, duration and swept footprint; runs overlap and one run may traverse many beam widths. An all-source estimate needs the dated angular union, rather than a tally of radar records. Nor can one use a 12,000-ly cone as the current frontier: a 1960s signal has traveled only about 60 ly.

Today’s global transmission rate remains unresolved.

Arecibo sends no new beams. JPL lists Goldstone DSS-14 offline for modernization from March 2026 to May 2028. Other antennas and facilities prevent inferring a zero global rate.

The available history counts detections and scheduled observations, not all transmissions or previously untouched directions.

The JPL detection history contains 2,082 station–apparition rows across 1968–2026. Goldstone-family distinct detected targets number 55 in 2024, 41 in 2025, and one in the retrieved 2026 snapshot. These are not independent cones, and transmissions producing no echo are absent. Scheduled observations may be cancelled. The worldwide new-direction/day/year quantity remains unresolved for DSN uplinks to all spacecraft, other planetary-radar facilities, and additional radar types.

Even all 84,054 advertised UCLA asteroid runs would omit lunar/planetary targets, failed observations, and other facilities. Required fields are target, UTC transmitter-on/off times, site, frequency/power, original pointing or orbit-predicts, and beam offsets. Actual DSN uplink duty logs are needed to test its continuous-transmission model.

For a conditional estimate of continuing spacecraft beams, the next table assumes continuous tracking from October 7, 2026 to October 7, 2027. New coverage is measured against all five modeled past spacecraft tracks, with each mission’s own beam width. It is not compared with every historical human transmission.

Mission Full beam width New native beam-disk equivalents/year Average/day
Voyager 1 7.68′ 2.7 0.007
Voyager 2 7.68′ 14.0 0.038
New Horizons 2.28′ 14.3 0.039

The combined increment is approximately 0.23 square degrees/year, or 0.0006 square degrees/day. The mission footprints are widely separated. Native beam-disk equivalents have different widths and cannot simply be summed; sky area provides the common unit. This scenario assumes a 70-m antenna, continuous uplink duty, and geocentric incoming spacecraft directions. It estimates geometry under those assumptions, not a measured current operations rate.

Methods and the meaning of a light cone

For a fixed beam with half-angle α and signal age τ in years:

Ω = 2π (1 − cos α)                 solid angle (steradians)
R = τ light-years                 leading-front radius
V_ever = Ω R³ / 3                 cumulative swept volume (ly³)
E[stars] = n_stars × V_union       compatible volume units required
E[HZ planets] = V_union × (n_M η_M + n_FGK η_FGK)
dE[stars]/dt = n_stars × Ω R² c    one nonoverlapping fixed beam
t_arrival ≈ t_emission + distance/c

A pulse emitted over an interval occupies a thin radial shell at a given time. “Ever touched” includes the interior through which that shell has already passed. It does not mean those stars are still receiving photons. The volume calculations remove the known empty region within 4.24 ly of the Sun; they are local-population expectations beyond it.

For moving beams, the study samples actual transmit intervals (Arecibo), nominal configuration envelopes (Goldstone), or continuous-tracking models (DSN). Dense points are no more than one quarter of the beam radius apart. A cap-mixture sampler chooses disks in proportion to their individual solid angles and weights each ray by the sum of all disk areas divided by its overlap multiplicity. The earliest included emission on a ray determines its leading front. This removes repeated angular coverage and earlier physical illumination in the defined disk inventory. Numerical integration uses 120,000 samples with analytical identical/disjoint/mixed-width controls, independent-seed checks of the radar union, and beam-width sensitivity. Discretized tracks still differ slightly from continuous tracks.

The radar-history extrapolation uses a first-coverage hazard g(t)/B in an ecliptic band of area B. Its first-emission density is [g(t)/B]exp(−∫g/B). Annual 16-node Gauss–Legendre quadrature integrates this density against the reached volume per ray. The retained DSN/message union is added; the expected minimum of the two fronts is subtracted on shared directions. Thus neither repeated radar coverage nor overlap with the backbone is counted twice. Arecibo’s area-per-target calibration already measures new coverage within its recovered history; applying the band hazard is a further explicit cross-family/older-history overlap approximation.

Density-based forecasts describe expectations under a homogeneous local population. Catalogue-based forecasts identify particular modeled targets; they are separate estimators and should not be summed as a global census. A directional Galactic population model is required for kiloparsec projections. A beam’s half-power contour is a declared geometric convention, not a hard electromagnetic edge.

The stellar centre is used as a proxy for planets orbiting it. At interstellar distances the nominal footprints are much larger than typical HZ orbits, but a grazing stellar-centre alignment does not guarantee that every planet is inside the chosen contour. Orbital phases and individual planetary positions were not modeled.

Geometric illumination also differs from detectability. Sheikh et al.’s conditional maxima use a receiver in the beam, matched bandwidth, sufficient integration and optimistic near-term telescope sensitivity: approximately 12,000 ly for former Arecibo planetary radar and 65 ly for DSN uplinks. They cannot be assigned to all waveforms/facilities, or mistaken for distances already traveled.

Editable source, data and verification

Analysis cutoff: 2026-10-07 00:00 UTC. All source archives were retrieved for this study. Numerical results are traceable to saved JSON/CSV exports and documented Python functions managed with uv. The selected DSN source is data/dsn/dsn_corrected_unique_objects_mixed_xband_linear_intersection_extended.csv. Independent review covers scientific assumptions, arithmetic, source alignment, rendered layout and controls; unresolved inventory limitations remain explicit.

Download the source Markdown, build manifest, combined recovered timeline, historical extrapolation timeline, day/week/year arrivals, global model manifest, future new-cone scenario, future-scenario manifest, annual activity calibration, explicit year imputations, population parameters, DSN arrival forecast, intentional target sessions, Goldstone volume forecasts, Arecibo archive volume forecasts, Arecibo annual new directions, Arecibo catalogue matches, Arecibo case forecasts, future DSN direction scenario, historical scenarios, and assumptions and source limitations.

The complete local project also retains raw archives, cached JPL requests/responses, primary papers, preserved original DSN results, alternative beam-width/epoch models, and the independent review record. No data points were inferred from screenshots. Plotly is embedded so this single HTML works offline; its runtime adds several megabytes.