Earth-Shadow Deficit: Astro & DM Effects

• Earth-shadow deficit is a modulation of incident flux caused by Earth's geometric shadow and atmospheric interactions, affecting photons, cosmic rays, and dark matter signals.
• It combines geometric principles like umbra and penumbra with refraction, extinction, and scattering, leading to measurable brightness reduction in satellite observations.
• The effect induces diurnal modulation in direct-detection experiments, offering a tool to constrain dark matter properties and understand terrestrial particle attenuation.

The Earth-shadow deficit refers to the reduction or modulation in particle flux—whether photons, cosmic rays, or hypothetical particles such as dark matter (DM)—due to the opacity or modifying influence of the Earth along the propagation path to a detector or observer. In astrophysical and experimental contexts, this effect is central to understanding not only the illumination of artificial satellites during eclipses (via atmospheric refraction and extinction) but also the expected diurnal variations in underground DM direct-detection experiments caused by Earth-shadowing and scattering. The phenomenon quantifies the position-dependent suppression (or, in some regimes, enhancement) of an incident flux and encodes both geometric shadowing and physical interactions with the terrestrial environment.

1. Geometric Foundations of the Earth-Shadow Effect

The geometric determination of Earth's shadow is foundational for both satellite photometry and penetrating particle/DM flux calculations. The classic umbra and penumbra are defined by the projection of Earth's disk against a distant radiating source (e.g., the Sun), parameterized by:

• Earth radius $R_e = 6371$ km, solar radius $R_\odot = 6.96\times10^5$ km, and Earth–Sun distance $D \approx 1.496 \times 10^8$ km,
• Umbral and penumbral cone half-angles: $\theta_u = \arctan[(R_e-R_\odot)/D]$, $\theta_p = \arctan[(R_e+R_\odot)/D]$,
• A satellite at position $(s, h)$, where $s$ is the distance along the Earth–Sun line from the shadow limb tangent and $h$ is the perpendicular height above the solar tangent line.

For low-Earth-orbit satellites, the umbra and penumbra boundaries are well approximated as cylindrical, and the in-shadow ($|h| \leq s \tan\theta_u$), penumbral ($s\tan\theta_u < |h| \leq s\tan\theta_p$), and fully illuminated ($|h| > s\tan\theta_p$) regions are delineated directly from this geometry (Mallama, 2021).

For ambient cosmic or DM flux, the geometric shadow is the simple line-of-sight chord through the Earth, parameterized by zenith angle $\theta$ with chord length $L(\theta) = 2R_\oplus \cos\theta$ for $\theta \leq 90^\circ$ (neglecting atmospheric layers), relevant for attenuation modeling (Bernabei et al., 2015).

2. Atmospheric and Terrestrial Modifications: Refraction, Absorption, and Scattering

Beyond simple geometric shadowing, the Earth's gaseous and solid material mediates interactions that alter the apparent "shadow" deficit. In the context of satellite brightness:

Atmospheric refraction bends sunlight into the geometric shadow, producing a nonzero illumination within the eclipse boundary. The reference model defines an atmospheric scale height $H \approx 6.58$ km. The "half-intensity ray" is determined as that where the refracted intensity drops to half: $$\delta(z_0) \simeq \nu(z_0) \sqrt{2\pi a/H} = \pi/2 \implies \nu(z_0) = H/D_\mathrm{sat}$$ yielding specific $z_0$ (altitude) values for a given satellite distance $D_\mathrm{sat}$.

Atmospheric extinction (primarily Rayleigh scattering) introduces a wavelength-dependent reduction of intensity, encoded as: $T(\lambda, z) = \exp[-k(\lambda) m(z)]$ where $k(\lambda)$ is the extinction coefficient per unit air mass, and $m(z)$ is the integrated horizontal air-mass along refracted paths. The extinction is stronger at shorter wavelengths, with transmission factors ranging from $E(B) \approx 0.73$ (B-band, $0.44\,\mu$m) to $E(i) \approx 0.95$ (i-band, $0.75$–$0.90\,\mu$m). This leads to a spectral reddening of brightness during eclipse ingress/egress (Mallama, 2021).

For DM and high-penetration cosmic rays, the solid Earth induces attenuation via elastic nuclear scattering. A particle with scattering cross section $\sigma_n$ off nuclei and local average Earth density $\rho_E$ will have a surviving flux (for incident zenith angle $\theta$): $$\phi_\text{out}(\theta) = \phi_\text{in} \exp[-\sigma_n \rho_E L(\theta)]$$ (Bernabei et al., 2015). The chord length $L(\theta)$ and compositional stratification of the Earth (e.g., core, mantle) can be incorporated via layered models or effective column densities.

3. Analytic and Numerical Modeling of the Deficit

The integrated effect of these processes on flux or brightness is model dependent.

For satellites, the normalized eclipse brightness $F(s,\lambda)$ is computed via convolution of the refracted (and extinguished) solar intensity over the visible solar disk, incorporating limb-darkening: $$F(s, \lambda) = \frac{1}{\pi R_\odot^2}\iint_{|\rho|<R_\odot} \left[\frac{I_{\rm ref}(h(\rho,s))}{I_0}\right]\left[E(\lambda)^{m(z(\rho,s))}\right]\frac{I(\phi)}{I_0} dA$$ The dimmed apparent magnitude is then: $m(s, \lambda) = -2.5 \log_{10}[F(s, \lambda)]$ This model was tabulated in 1 km steps in $h$ for specific satellite distances and multiple photometric bands (Mallama, 2021).

For galactic dark matter, a single-scatter analytic approach models the perturbed velocity distribution at the detector: $$\tilde{f}(\mathbf{v}, \gamma) = f_0(\mathbf{v})\, e^{ -\sum_i \frac{d_{\mathrm{eff},i}(\cos\theta)}{\bar\lambda_i(v)} } + \sum_i \int d\Omega' \frac{d_{\mathrm{eff},i}(\cos\theta)}{\bar\lambda_i(v')} \frac{v'}{v} f_0(\mathbf{v}') P_i(\mathbf{v}' \to \mathbf{v})$$ where $f_0$ is the initial Maxwellian distribution, the first term ($f_A$) represents attenuation, and the second ($f_D$) describes the gain from deflected particles (Kavanagh et al., 2016). The formulation respects detailed balance as required for elastic scattering.

4. Diurnal Modulation and Experimental Observations

A critical implication of the Earth-shadow deficit is the generation of time-dependent modulation in observed rates or apparent brightness as the Earth rotates:

• For satellites: The eclipse deficit evolves on orbital timescales (tens of minutes), and was validated via lightcurve measurements with the MMT-9 robotic telescope, finding observed dimming slopes consistent with model predictions and an RMS residual of 0.3–0.5 mag in the V band. The light curves lie well below the geometric (no atmosphere) prediction, confirming the necessity of including atmospheric refraction and extinction (Mallama, 2021).
• For DM detection: The daily modulation arises from the changing path length through Earth as the laboratory rotates with respect to the incoming DM wind. The time-dependent event rate is: $$R(t) = R_0 [1+\mathcal{A}(\lambda)\cos(\omega t-\phi(\lambda))]+\ldots$$ with $\mathcal{A}$ and $\phi$ encoding amplitude and phase, both strongly dependent on latitude and DM–nucleon interaction operator. Typical amplitudes range from $\sim2\%$ enhancement (Gran Sasso, $\mathcal{O}_1$), to $10$–$30\%$ suppression (Southern hemisphere, velocity-suppressed operators), with distinctive phase shifts correlated to the form of the interaction (Kavanagh et al., 2016).

In direct detection data (e.g., DAMA/LIBRA phase1), systematic searches for a diurnal rate modulation (amplitude $A_d$) yield upper limits: $A_d \lesssim 2\times10^{-3}$ cpd/kg/keV, excluding substantial regions of $(m_\mathrm{DM}, \sigma_n, \xi)$ parameter space for high cross section, low halo-fraction scenarios (Bernabei et al., 2015).

5. Dependences on Interaction Physics and Geospecific Parameters

The magnitude and qualitative nature of the Earth-shadow deficit are sensitive to:

• Wavelength (satellite photometry) via the extinction coefficient $k(\lambda)$, resulting in stronger deficit in blue bands.
• DM–Nucleon operator structure: Spin-independent and isotropic operators ($\mathcal{O}_1$, $\mathcal{O}_4$) yield moderate attenuation and near-constant enhancement; velocity-suppressed and momentum-suppressed operators ($\mathcal{O}_8$, $\mathcal{O}_{12}$) produce pronounced daily modulations and distinctive phase behavior due to their angular scattering properties.
• Detector latitude and geography: Mid-latitude detectors observe modest deficit or enhancement; Southern hemisphere sites may see large suppression. The phase of modulation shifts with latitude and interaction type (Kavanagh et al., 2016).
• Satellite orbital parameters and local terrain/cloud effects: The effective height in the eclipse tables can be shifted by $Δh ≈ H·\ln(D_\mathrm{sat}'/D_\mathrm{sat})$ and a terrain/cloud cutoff term $\varepsilon$.

6. Experimental Constraints and Modeling Tools

The empirical detection or constraints on the Earth-shadow deficit yield nontrivial bounds on astrophysical and particle models:

• In satellite observations, the deficit tables provided by Mallama (Mallama, 2021) enable prediction of dimming for arbitrary low-Earth-orbit trajectories and photometric bands.
• For DM, the lack of an observed diurnal Earth-shadow deficit in DAMA/LIBRA-phase1 data imposes stringent upper limits on the halo fraction $\xi$ for large cross-section DM candidates. A statistically significant diurnal modulation ($\Delta R/R_0 \gtrsim 5\%$) is within reach of current low-threshold experiments (CRESST-III, SuperCDMS SNOLAB, DAMIC) with upcoming exposures (Bernabei et al., 2015, Kavanagh et al., 2016).
• The EarthShadow code (publicly released) computes the post-Earth-scattering DM velocity distribution for arbitrary detector location, allowing for realistic Earth compositional models and interaction operators (Kavanagh et al., 2016).
• Linear interpolation between tabulated values in $(h,\lambda)$ or over filter bands is straightforward for both satellite and DM scenarios (Mallama, 2021).

7. Outlook and Limitations

The analysis and interpretation of Earth-shadow deficits are increasingly relevant for precision astrometric catalogs, LEO satellite tracking, and the next generation of dark matter direct detection at low energy thresholds. The main limitations stem from uncertainties in high-altitude atmospheric structure, variability in terrain/cloud cutoff for photometry, and inhomogeneities in Earth's density for particle astroparticle attenuation. Models continue to be refined to address these uncertainties with improved observational validation and more granular mapping of both the physical and parameteric dependencies of the deficit.

The Earth-shadow deficit remains a crucial diagnostic for identifying and constraining astrophysical signals subject to terrestrial mediation, as well as a testing ground for the interplay between radiative transfer physics, celestial geometry, and underground particle physics (Mallama, 2021, Bernabei et al., 2015, Kavanagh et al., 2016).