  Dynamic rupture simulations generate synthetic waveforms that account for
non-linear source, path, and site complexity. Here, we analyze millions of
spatially dense waveforms from 3D dynamic rupture simulations in a novel way to
illuminate the spectral fingerprints of earthquake physics. We define a
Brune-type equivalent near-field corner frequency ($f_c$) to analyze the
spatial variability of ground motion spectra and unravel their link to source
complexity. We first investigate a simple 3D strike-slip setup, including an
asperity and a barrier, and illustrate basic relations between source
properties and $f_c$ variations. Next, we analyze > 13,000,000 synthetic
near-field strong motion waveforms generated in three high-resolution dynamic
rupture simulations of real earthquakes, namely, the $M_w$ 7.1 2019 Ridgecrest
mainshock, the $M_w$ 6.4 Searles Valley foreshock, and the $M_w$ 7.3 1992
Landers earthquake. All scenarios consider 3D fault geometries, topography,
off-fault plasticity, viscoelastic attenuation, 3D velocity structure, and
resolve frequencies up to 1-2 Hz. Our analysis reveals pronounced and localized
patterns of elevated $f_c$, specifically in the vertical components. We
validate such $f_c$ variability in observed near-fault spectra. Using isochrone
analysis, we identify the complex dynamic mechanisms that explain the rays of
elevated $f_c$ and cause unexpectedly impulsive, localized, vertical ground
motions. While the vertical high frequencies are also associated with path
effects, rupture directivity, and coalescence of multiple rupture fronts, we
show that they are dominantly caused by rake-rotated surface-breaking rupture
fronts that decelerate due to fault heterogeneities or geometric complexity.
Our findings highlight the potential of spatially dense ground motion
observations for furthering our understanding of earthquake physics directly
from near-field data.
