Natural Variability and Non-CO₂ Drivers

Underrepresentation of Solar Cycles, AMO, PDO, and ENSO

Solar cycles Solar variability (particularly the Schwabe (11-year) and Gleissberg (~80-90 year) cycles) modulates the total solar irradiance (TSI) reaching Earth, affecting stratospheric ozone, atmospheric circulation, and possibly cloud nucleation through solar magnetic activity and cosmic ray flux modulation.

Most GCMs treat solar forcing as quasi-static, using smoothed TSI averages or linearized inputs. This neglects shorter-term solar variability, including potential feedbacks via the Quasi-Biennial Oscillation (QBO) and stratosphere-troposphere coupling. As a result, models understate the potential role of solar minima (e.g., Maunder Minimum) in shaping decadal-to-centennial climate anomalies.

Atlantic Multidecadal Oscillation (AMO) and Pacific Decadal Oscillation (PDO) These large-scale ocean oscillations, governed by thermohaline circulation and internal ocean-atmosphere feedbacks, drive multi-decadal climate variability, influencing drought frequency, tropical cyclone activity, Arctic sea ice extent, and precipitation regimes.

Most GCMs fail to simulate the amplitude, phase, and teleconnections of these oscillations. This failure is compounded by:

• Lack of ocean data for model initialization
• Insufficient spatial resolution
• Ensemble averaging that smooths out internal variability

El Niño–Southern Oscillation (ENSO) ENSO is the dominant mode of interannual climate variability, altering global rainfall, temperature, and storm patterns.

CMIP-class models exhibit systematic biases:

• Failure to reproduce ENSO frequency and periodicity
• Underestimation of extreme El Niño/La Niña amplitudes
• Incorrect spatial patterning of SST anomalies (“cold tongue bias”)

These limitations distort the modeled climate signal and blur the boundary between internal variability and anthropogenic trends.

Misattribution of Decadal and Regional Variability to Anthropogenic CO₂

Decadal variability Empirical records show pronounced decadal-scale swings in temperature and climate anomalies, including:

• The 1940s-1970s global cooling trend, concurrent with post-war aerosol emissions and solar minima
• The 1998-2014 “global warming hiatus”, during which surface temperature trends flattened despite rising CO₂

GCMs consistently fail to replicate the timing or magnitude of these events without extensive tuning. This disconnect suggests model overfitting to monotonic forcing assumptions and underweighting of internal ocean-atmosphere modes.

Regional patterns Heatwaves, droughts, and cold anomalies (such as the 2010 Russian drought, 2015 California drought, or 2021 Texas freeze) are often retrospectively attributed to climate change. In many cases, however, regional circulation anomalies (e.g., jet stream displacement, blocking highs, ENSO teleconnections) offer more plausible explanations.

Policy missteps Flawed attribution can lead to:

• Overconfidence in long-term risk assessments
• Misallocated infrastructure resilience investments
• Premature regulatory mandates grounded in short-term anomalies

Parameterizing Sun-Earth Interactions as Static

Static assumptions Solar irradiance is often modeled as a linear or flat baseline trend, ignoring:

• Short-term TSI fluctuations
• Magnetic field modulations
• Solar spectral variability, especially in the UV band
• Solar wind and geomagnetic interactions with the upper atmosphere

Most models exclude cosmic ray-cloud interaction hypotheses, despite growing observational and theoretical support for such mechanisms in cloud condensation nuclei (CCN) modulation.

Physical limitations in models Indirect solar effects on ozone chemistry, meridional circulation, and stratospheric cooling are underrepresented. GCMs typically lack vertical resolution in the stratosphere to capture these dynamics, especially over polar regions where solar-terrestrial interactions are strongest.

Underestimation of role By flattening Sun-Earth interactions, models may misrepresent the attribution of both 20th-century warming and natural climate oscillations, especially over multi-decadal periods.

Climate Variability vs. Mean-State Bias in Models

Mean-state calibration GCMs are calibrated primarily to reproduce global mean surface temperature trends. This metric obscures the structure and dynamics of variability, leading to a systemic underappreciation of extremes and anomalies.

Underestimated variability Studies show that models:

• Underpredict the amplitude of ENSO and AMO cycles
• Miss low-frequency variability in polar and mid-latitude systems
• Fail to reproduce observed variance in surface temperature records, precipitation, and storm frequency

Blind spots in extreme event risk By focusing on average change, models may miss nonlinear threshold behavior and “black swan” events (rare but high-impact shifts like sudden polar vortex breakdowns, monsoon failures, or abrupt sea ice loss).

Extremes May Rise Without a Shift in the Mean

Empirical evidence Observed climate records show rising frequency, duration, and intensity of extreme weather events in some regions, even where the mean climate has changed little or not at all.

Examples:

• Increased heatwave duration without proportional rise in average temperature
• Rising flash flood frequency without consistent precipitation trends
• More frequent blocking patterns in mid-latitudes

Modeling gaps GCMs are not built to resolve:

• Subgrid-scale triggers of extremes (e.g., local convection, mesoscale circulation)
• Compound events driven by coincident drivers (e.g., drought + heatwave)
• Tail-risk behavior in atmospheric dynamics

This leaves extreme event modeling empirically unsupported, especially at the spatial and temporal scales relevant for risk management and adaptation policy.

Synthesis and Implications for Policy and Research

Natural variability is systematically underrepresented The persistent omission or simplification of solar variability, ocean-atmosphere cycles, and internal climate dynamics results in models that over-attribute observed changes to CO₂ forcing. This narrows the interpretation of past trends and artificially inflates future risk projections.

Attribution remains a core scientific challenge Disentangling the anthropogenic signal from natural variability requires:

• High-resolution, initialized decadal models
• Expanded use of observation-constrained hindcasts
• Acknowledgment of uncertainty ranges in attribution science

Research priorities:

• Improved model initialization with historical ocean-atmosphere states
• Expanded solar forcing representation, including indirect pathways
• Higher-resolution GCMs to capture mesoscale variability and extremes
• Methodological pluralism: including stochastic models, empirical reconstructions, and reduced-form frameworks