Fractional Empirical Orthogonal Functions for Geophysical Fields with Anomalous Transport: Theory and Validation

Abstract. Empirical Orthogonal Function (EOF) analysis and its rotated variant (REOF) are foundational tools in the geosciences for decomposing spatiotemporal variability. However, the standard methodology implicitly assumes Gaussian statistics and exponentially decaying correlations, assumptions that are violated in many geophysical systems exhibiting anomalous diffusion, heavy-tailed distributions, and long-range spatial correlations. We develop a theoretical framework for fractional EOF (fEOF) analysis that extends the standard methodology by incorporating the fractional Laplacian operator into the covariance structure. The governing dynamics are formulated using the Riemann–Liouville fractional time derivative of order μ > 0, which is not restricted to the interval (0,1] and thereby accommodates both subdiffusive and superdiffusive transport regimes within a single formalism. The resulting fractional covariance operator naturally captures power-law correlations characteristic of anomalous transport in geophysical flows. We prove that the eigenvalue spectrum of the fractional covariance operator exhibits enhanced power-law decay λ_m^(α) ∼ m^{−(1+α+β/d)}, where the spatial fractional order α ∈ (0,2) provides a tunable control parameter independent of the underlying spectral slope β. The temporal evolution of fractional principal components follows Mittag-Leffler relaxation, interpolating between stretched exponential and power-law regimes. We validate the theoretical predictions through three independent approaches: (i) exact analytical results for fractional Brownian surfaces across three Hurst exponents (H = 0.3, 0.5, 0.7), confirming eigenvalue steepening to within 6 % of theoretical predictions with finite-domain corrections identified; (ii) spectral analysis of fields generated by the space-time fractional diffusion equation across spectral slopes β = 2, 3, 4, recovering predicted exponents to within 2–6 %; and (iii) Monte Carlo experiments over 50 realizations demonstrating that the eigenvalue scaling is distribution-independent, holding identically for Gaussian and heavy-tailed Student-t fields while achieving 5- to 8-fold reductions in the number of modes required for 95 % variance capture. The sensitivity analysis across fractional orders α ∈ [0, 1.8] confirms the predicted linear steepening relation ν_α = ν₀ + α to within 3 % throughout. The framework applies to a broad class of geophysical fields exhibiting anomalous transport, including oceanic tracer dispersion, flood inundation dynamics, atmospheric constituent spreading, and soil moisture redistribution. Connections to Okubo's empirical oceanic diffusion scaling and the Forecasting Inundation Extents using REOF (FIER) framework are discussed as illustrative applications.