Total: 1
Spatial coupling and heterogeneity can generate complex power-grid frequency dynamics. We investigate the contributions of spatial structure, intrinsic nonlinear dynamics, and external driving using the Kuramoto model (KM), with and without uncorrelated noise, and compare its predictions with real-grid measurements. We analyse one month of 100 ms frequency measurements from 16 Hungarian transmission-grid sites and simulate the KM on a 387-node transmission network and two-dimensional lattices. The measured power spectra exhibit two scaling regimes, with exponents of approximately 1.63 and 2.23. One-second frequency increments decorrelate rapidly, while correlations of their local inverse variance decay slowly, initially as a power law, $τ^{-0.26(2)}$, up to 15 minutes. A superstatistical Gaussian mixture based on measured local variances captures much of the non-Gaussian increment distribution, but not the rare extreme events. We also estimate the non-stationary long-time frequency autocorrelation exponent. KM simulations on low-dimensional networks exhibit algebraically decaying autocorrelations and power spectra, consistent with the measurements and the Wiener-Khinchin relation in the steady state. On the heterogeneous transmission network, partial synchronization and frequency locking produce persistent spatial patterns and multi-peaked frequency distributions. However, neither deterministic KM dynamics nor independent Gaussian driving reproduces the measured spectra and increment statistics simultaneously. Moreover, empirical frequency fluctuations are largely shared across substations, whereas simulated fluctuations mainly reflect relative node motion. These findings suggest that spatially coupled nonlinear dynamics contribute to observed frequency fluctuations alongside external driving, including 15-minute power dispatches.