Lunar-Nodal Climate Hypothesis

MPRC-Derived Modulation of South Asian Hydroclimate

South Asian climate contains several unresolved variability patterns, including a seasonal temperature inversion in the Upper Indus Basin (UIB) and contrasting Pakistan hydroclimate outcomes during La Niña and solar ascending phases: near-identical forcing produced severe drought in 1998–2001 and catastrophic flooding in 2010 and 2022.

We test whether the 18.6-year lunar nodal cycle, represented by a zero-parameter MPRC-derived geometric index:

is a candidate modulator of these outcomes. The 1998–2001 drought occurred near nodal minimum (MF 1.00–1.04), whereas the catastrophic 2010 and 2022 floods occurred at elevated MF (1.28 and 1.30 respectively), motivating the hypothesis.

Arabian Sea JJAS SST anomaly is the strongest additional predictor by AIC. MF does not add statistically significant predictive skill for major flood occurrence in the 44-year record (phase-randomised null test p = 0.806), and the best-fit period scan identifies 21.5 years rather than 18.6 years. However, the MF coefficient is positive and stable across all leave-one-event-out iterations (range [1.46, 2.99], full = 2.06). The null test failure is consistent with expected low power: the record spans 2.4 nodal cycles, whereas reliable period detection requires approximately five complete cycles (~93 years).

The Pakistan Monsoon Behavioural Split

Pakistan and the Upper Indus Basin occupy the western edge of the South Asian monsoon. Multiple independent studies document a seasonal temperature inversion: winter maximum temperatures rising at +0.38°C per decade while summer maximum temperatures cool at −0.14°C per decade (Fowler and Archer 2006; Iqbal et al. 2021). Pre-monsoon temperatures have increased uniformly across all regions (Sheikh et al. 2009). These patterns are not reproduced by uniform greenhouse gas forcing scenarios.

The behavioural split in the Pakistan monsoon record is the primary motivation for this study. La Niña is a well-established precursor of enhanced Pakistan monsoon rainfall (Adnan et al. 2021). The ascending solar cycle phase is associated with intensified atmospheric circulation (Gray et al. 2010). Yet the coincidence of both conditions in 1998–2001 produced Pakistan's worst multi-year drought since independence (GCISC 2009), while the same joint condition produced catastrophic flooding in 2010 (1,985 deaths, $9.7B damage) and 2022 (1,739 deaths, $40B damage).

Existing attribution studies for 2022 emphasise anthropogenic warming, La Niña circulation, and subtropical high behaviour (Nanditha et al. 2023; Zhou et al. 2023) but do not identify a variable that explains why the same ENSO and solar configuration produced drought in 1998–2001.

We quantify this via a dimensionless index MF, derived from the MPRC geometric framework (Arshad 2026), that cycles between 1.00 (minor standstill) and 1.40 (major standstill) with zero free parameters. This paper reports both the motivating event-phase pattern and the results of formal statistical tests. We do not claim the mechanism is established; we claim it is testable, motivated by an unexplained event contrast, and falsifiable on a near-term observational timescale.

The Zero-Parameter Nodal Index

Climate tests in this paper do not require acceptance of the full MPRC framework; they require only that MF is a reproducible, zero-parameter time series with the 18.612-year nodal period.

Operational Definition

The Moon's orbital inclination to Earth's equator cycles between 18.3° at minor standstill and 28.6° at major standstill over the 18.612-year nodal precession cycle. Confirmed recent standstills: minor October 1997, major mid-2006, minor October 2015, major early 2025. The sum angle and MF are defined as:

The 40.3% range follows from the confirmed standstill inclination limits (18.3° and 28.6°) and Earth's fixed axial tilt (23.44°). No fitting. No free parameters. All MF values are independently reproducible from any standard lunar ephemeris.

Physical Motivation

The sin-squared form follows the same angular dependence used in tidal constituent nodal factor analysis (Haigh et al. 2011): the nodal modulation of tidal forcing is proportional to the squared sine of the inclination angle. MF is therefore consistent in functional form with, though not calibrated to, standard tidal harmonic nodal factors.

The MPRC framework (Arshad 2026, Chapter 7) provides a geometric derivation of the sin-squared partition from first principles, independently verified against Newton's law on eleven gravitational datasets. Whether or not the MPRC framework is accepted, MF is a well-defined, reproducible quantity.

Figure 1 (described): Physical basis of MF. (a) The MPRC gravitational partition from Chapter 7: T_drag = GMm/r² partitions via cos²(θ) + sin²(θ) = 1 (verified to ratio 1.00000000 on 11 datasets). (b) θ_sum variation over the nodal cycle with confirmed standstill dates. (c) MF showing the 40.3% range between standstills.

South Asian Climate Unresolved Patterns

Each anomaly below is cited from peer-reviewed literature or primary government sources. No findings from the present analysis are introduced in this section.

UIB Seasonal Temperature Inversion

Iqbal et al. (2021) document winter T_max +0.38°C/decade and spring T_max +0.35°C/decade, while summer T_max declined at −0.14°C/decade and autumn T_min at −0.22°C/decade across nine UIB stations for 1955–2016. Fowler and Archer (2006) confirm the pattern for 1961–2000. This seasonal inversion has no established causal mechanism.

Pre-Monsoon Thermal Intensification

Sheikh et al. (2009) found April–May temperatures increased in all regions across 52 stations for 1961–2000. The spatial uniformity suggests a common large-scale forcing.

The 1998–2001 La Niña Drought Paradox

Pakistan's 1998–2001 drought was the worst in recorded history (GCISC 2009). It occurred during confirmed La Niña conditions and the ascending phase of Solar Cycle 23 — the same joint configuration that produced major flooding in 1988 and catastrophic flooding in 2010 and 2022. No published study identifies why 2010 — also following a weak El Niño in 2009 — was catastrophic under the same configuration.

The 2010–2022 Catastrophe Clustering

The two worst Pakistan floods in recorded history are separated by exactly 12 years — one Schwabe solar cycle. Both occurred during solar ascending phase with active monsoon-season La Niña. Attribution studies for 2022 identify anthropogenic warming, La Niña circulation, and subtropical high behaviour as contributing factors. The present paper proposes the nodal cycle as an additional candidate modulating variable, complementary to these explanations.

SC22 versus SC25 Severity Asymmetry

Solar Cycles 22 and 25 are 33 years apart with near-identical peak sunspot numbers (158 vs 157). SC25 produced a catastrophic event (2022) while the SC22 epoch produced only major events (1988 major, 1992 severe). A variable beyond the solar cycle is required to explain this asymmetry.

Statistical Framework and Null Tests

Event Catalogue, 1982–2025

The primary analysis uses the 1982–2025 period (n = 44 years). NOAA OISST Arabian Sea SST data begin in September 1981 and the full regression panel requires complete covariate coverage. A supplementary descriptive analysis uses the 1971–2026 record for event-phase visualisation.

Flood severity is classified on a 0–3 ordinal scale: 0 = no significant event, 1 = major flood, 2 = severe flood, 3 = catastrophic. Classification sources: Federal Flood Commission Annual Reports (primary, 1971–2014), NDMA annual reports, UN OCHA situation reports (2010, 2022), and published academic literature. The binary outcome event_binary = 1 when severity ≥ 1. A separate catastrophic_binary identifies severity = 3 events (two cases: 2010 and 2022).

Covariates

Regression Models

Logistic regression (statsmodels, Python) with Firth penalisation for rare events. Model comparison by AIC. Thirteen model specifications tested. Standard errors from the observed information matrix; Wald p-values reported with the acknowledgment that the dataset is small and asymptotic approximations may be unreliable.

Null Tests

Four null tests were applied to the primary specification (event_binary ~ la_nina_jjas + solar_ascending + MF + jjas_sst_arabian_sea_anom), tested against the baseline (event_binary ~ la_nina_jjas + solar_ascending + jjas_sst_arabian_sea_anom):

Event-Phase Contrast and Statistical Tests

Event-Phase Contrast — Motivating Pattern

The contrast between La Niña + solar ascending years is visually clear: years at MF ≤ 1.09 produced drought or no event; years at MF ≥ 1.15 produced major to catastrophic events.

The 1989 anomaly (MF 1.39, La Niña weak by JJA, ASC, no major flood) and the 1988 major flood (not catastrophic despite MF 1.40) illustrate the limits of the three-factor framework. These are noted as open problems.

Seasonal Inversion Consistent with Nodal Phase

The UIB seasonal inversion gap acceleration after the mid-1990s is consistent with MF rising from the 1997 minor standstill toward the 2006 major standstill. Prediction: the gap should moderate toward 2034 as MF declines.

Regression and Null Test Results, 1982–2025

Arabian Sea JJAS SST anomaly is the strongest additional predictor by AIC. Adding MF to the baseline model does not improve AIC for major flood prediction. The MF coefficient is positive across all specifications but does not reach standard significance thresholds in the 44-year record.

Power Context

The null test failures must be interpreted in the context of statistical power. Reliable detection of an 18.6-year periodic signal via phase-randomisation tests requires approximately five complete cycles (~93 years). The 1982–2025 record spans 44 years = 2.4 nodal cycles. The estimated power of the phase-randomisation test at this record length is approximately 20–25% for a true signal of the observed effect size.

SC22 versus SC25 Asymmetry

The 33-year solar period asymmetry between SC22 and SC25 is consistent with the different nodal phases at the time of La Niña conjunctions. During SC22, La Niña in 1988–1989 occurred at MF = 1.40–1.39 (near the 1988 major standstill) but La Niña intensity was moderate (JJA ONI ~−0.8). During SC25, the sustained La Niña of 2020–2022 coincided with elevated MF throughout the 2006–2025 major standstill window, reaching MF = 1.30 in 2022 alongside strong La Niña (JJA ONI ~−1.4).

Near-Term Observational Tests

Claims, Limitations, and Required Next Steps

What This Paper Claims and Does Not Claim

Required Next Steps

Arabian Sea tidal analysis: Compute constituent-specific nodal factors for M2, K1, O1, and Mf in the Arabian Sea from TPXO9 tidal constituents or FES2014. Quantify whether the nodal modulation of internal tide generation at the Arabian Sea thermocline follows the MF proxy or a different functional form.

Mixed-layer depth verification: Extract an 18.6-year signal from Argo mixed-layer depth and ORAS5 ocean reanalysis in the Arabian Sea after removing ENSO, IOD, greenhouse trend, aerosol forcing, and internal variability. This is the critical mechanistic link.

Extended historical catalogue: Extend the Pakistan flood catalogue to 1950 using FFC historical records and British-era Indus discharge data (4.1 nodal cycles). This would approximately double the statistical power of the phase-randomisation test.

Regional model experiment: Run a regional Arabian Sea ocean model with and without the 18.6-year variation in tidal mixing amplitude, forced by observed ENSO boundary conditions, and compare Pakistan rainfall anomalies for 1996–2005 versus 2006–2015.

The MPRC Framework as Theoretical Motivation

The sin-squared functional form of MF is derived from the MPRC gravitational thread partition (Arshad 2026, Chapter 7), which recovers Newton's law exactly on eleven independent datasets. The framework provides a geometric reason why the perpendicular component of the gravitational thread scales as sin²(θ_sum). This derivation motivates the form of MF; it does not by itself establish that MF modulates Arabian Sea mixing or Pakistan flood severity. The climate mechanism must be established through the oceanographic and atmospheric tests described above.

An Exploratory Falsifiable Hypothesis

We present an exploratory hypothesis: that the 18.6-year lunar nodal cycle, quantified by a zero-parameter MPRC-derived index MF, is a candidate modulator of South Asian hydroclimate extremes. The motivating observation is that La Niña and solar ascending phase produced drought when MF was near its minimum (1998–2001, minor standstill) and catastrophic flooding when MF was elevated (2010, 2022, post-major standstill). This contrast is visually clear and sign-stable across leave-one-event-out samples.

Formal regression and null tests on the 1982–2025 panel do not confirm statistically significant predictive skill for MF in major flood occurrence. Arabian Sea JJAS SST anomaly is the stronger observed predictor. The null test failures are consistent with expected low power from a 44-year record spanning only 2.4 nodal cycles; reliable period detection requires approximately 93 years.

The critical analytical next steps are TPXO Arabian Sea tidal constituent analysis, Argo mixed-layer depth verification, and extension of the flood catalogue to 1950.

Data and Code Availability

Analysis code, MPRC_CORE.py constants registry, event catalogue, null test scripts, and figure generation code are available from the author. MF time series are independently reproducible from any standard lunar ephemeris using:

OISST Arabian Sea SST data: NOAA OISST v2.1, ncei.noaa.gov. Solar cycle data: NOAA WDC-SILSO, Brussels. ENSO ONI: NOAA CPC. Flood catalogue: FFC Annual Report 2014 (primary 1971–2014), NDMA/UN OCHA (2015–2026).

References

Adnan S, Ullah K, Shuanglin L, Gao S, Khan AH, Mahmood R (2021) Relationship between ENSO and Pakistan rainfall and the influence of the QBO and geomagnetic activity. Int J Climatol 41:S1513–S1527.

Arshad M (2026) The MPRC Framework — A Discrete Geometric Theory of Nature. April 2026. muhammadarshad.github.io/pages-mprc/.

FFC — Government of Pakistan, Federal Flood Commission (2014) Annual Flood Report 2014. Federal Flood Commission, Islamabad.

Fowler HJ, Archer DR (2006) Conflicting signals of climatic change in the Upper Indus Basin. J Clim 19(17):4276–4293.

GCISC (2009) Climate of Pakistan — Trends and Variability. Global Change Impact Studies Centre, Islamabad.

Gray LJ, Beer J, Geller M et al. (2010) Solar influences on climate. Rev Geophys 48(4):RG4001.

Haigh ID, Eliot M, Pattiaratchi C (2011) Global influences of the 18.61 year nodal cycle and 8.85 year cycle of lunar perigee on high tidal levels. J Geophys Res 116:C06025.

Iqbal MF, Liu J, Li J, Saeed S, Bilal M, Shafiq MU (2021) Observed trends and variability of temperature and precipitation and their global teleconnections in the Upper Indus Basin, Hindukush-Karakoram-Himalaya. Atmosphere 12(8):973.

Nanditha JS, Kushwaha AP, Singh R et al. (2023) The Pakistan flood of August 2022: causes and implications. Earths Future 11(3):e2022EF003230.

Sheikh MM, Manzoor N, Adnan M et al. (2009) Climate profile and past climate changes in Pakistan. GCOS-Pakistan final report, GCISC-RR-08.

van Loon H, Meehl GA (2012) The Indian summer monsoon during peaks in the 11-year sunspot cycle. Geophys Res Lett 39:L13701.

Zhou S, Huang G, Huang P (2023) Changes in the East Asia-Pacific teleconnection and influences on summer precipitation over the Yangtze River basin. Clim Dyn 60:2151–2166.

MPRC Gravitational Partition — Derivation of MF

This appendix provides the MPRC derivation for readers who wish to evaluate the theoretical basis of MF. It is not a prerequisite for the climate results in Sections 3–6.

The Thread Force and Frame Conversion

Chapter 7 derives the gravitational thread force T_drag and proves its reduction to Newton's law under frame conversion:

The Sin-Squared Partition

The frame conversion identity cos² + sin² = 1 partitions T_drag into two orthogonal components:

This ratio is derived from the partition proof in Chapter 7. It varies between 1.000 (minor standstill, θ_sum = 41.74°) and 1.403 (major standstill, θ_sum = 52.04°) with zero free parameters.

Relationship to Standard Tidal Theory

MF is consistent in functional form with standard tidal harmonic nodal factors, which scale with sin or sin² of the lunar inclination for different constituents. MF is not calibrated to specific tidal constituents and should not be equated with constituent-specific nodal factors for the Arabian Sea. A constituent-specific calibration using TPXO9 is identified as required work in Section 7. MF serves as a first-order geometric proxy for the nodal modulation amplitude pending that calibration.