Vertical Velocity Anomaly (VVA)

• Vertical Velocity Anomaly (VVA) is a deviation of vertical motion from a chosen reference state, exemplified in localized stellar kinematics, fluid turbulence, and atmospheric drafts.
• Different disciplines define VVA operationally—using local spatial backgrounds, model-predicted fields, or climatological means—to detect perturbers, diagnose intermittency, or quantify model deficiencies.
• Studies of VVA employ advanced statistical and analytical methods to reveal links between anomalous vertical motions and phenomena such as dark matter subhalo impacts, supergranular upflows, and turbulent intermittency.

Vertical velocity anomaly (VVA) is a context-dependent term for a physically significant departure in vertical motion from an adopted reference state. In Galactic dynamics, the term has been used for a localized perturbation in the mean stellar vertical velocity $\langle V_z\rangle$ and was reported as a coherent dip of $\simeq 7\,\mathrm{km\,s^{-1}}$ in the Milky Way disk near $(\ell,b,D)\simeq(15.8^\circ,-0.5^\circ,3.9$–$4.0\,\mathrm{kpc})$ (Udagawa et al., 27 Aug 2025). In other literatures, the same label is applied to heavy-tailed wall-normal fluctuations in turbulent wall flows, subsurface supergranular upflows in helioseismology, residuals between measured and quasigeostrophic predictions in laboratory vortices, departures from climatological oceanic vertical velocity, and statistically extreme drafts in the mesosphere and lower thermosphere (Banerjee et al., 19 Jun 2026, Svanda, 2012, Aulnette et al., 19 Dec 2025, Gao et al., 25 Sep 2025, Chau et al., 10 Dec 2025). The term therefore denotes a family of anomaly diagnostics rather than a single invariant observable.

1. Terminological scope and operational definitions

Across the cited literature, VVA is defined operationally rather than universally. The reference state may be a local spatial background, an axisymmetric Galactic model, a climatological mean, a model-predicted field, or Gaussian reference statistics. This variation is not incidental: it reflects different inferential goals, including detection of perturbers, diagnosis of intermittency, and quantification of model deficiency.

Domain	Operational definition	Representative paper
Galactic disk	Localized region where mean stellar vertical velocity differs significantly from the surrounding disk	(Udagawa et al., 27 Aug 2025)
Solar-neighbourhood Milky Way mapping	$$\delta V_z(R,\phi,z)\equiv\langle V_z\rangle(R,\phi,z)-\langle V_z\rangle_{\rm axisym}(R,z)$$	(Widmark et al., 2022)
Wall turbulence	$FF_w\equiv \langle w'^4\rangle/\sigma_w^4$, with $FF_w>3$ measuring heavy tails or intermittency	(Banerjee et al., 19 Jun 2026)
Solar supergranulation	Deviation of $v_z$ from zero mean, typically upflow at the cell center and downflow near cell boundaries	(Svanda, 2012)
Quasigeostrophic vortices	$$\mathrm{VVA}(x,y,z)\equiv w_{\rm meas}(x,y,z)-w_{\rm pred}(x,y,z)$$	(Aulnette et al., 19 Dec 2025)
Ocean reconstruction	$\mathrm{VVA}(x,y,z,t)=w(x,y,z,t)-\overline{w}(z)$	(Gao et al., 25 Sep 2025)
MLT rogue drafts	Extreme events with $\simeq 7\,\mathrm{km\,s^{-1}}$0 and anomaly parameter $\simeq 7\,\mathrm{km\,s^{-1}}$1	(Chau et al., 10 Dec 2025)

A common source of ambiguity is that the same acronym may refer either to a field variable itself, to a residual after subtraction of a reference model, or to a higher-order statistic of fluctuations. In the Galactic-disk discovery paper, VVA denotes a localized kinematic feature in stars (Udagawa et al., 27 Aug 2025). In Banerjee et al., by contrast, VVA is identified with intermittency in the flatness factor of wall-normal velocity fluctuations (Banerjee et al., 19 Jun 2026).

2. Localized stellar VVA in the Milky Way disk

In the Galactic-disk usage of Udagawa et al., a vertical velocity anomaly is defined as a localized region in which the mean vertical velocity of stars differs significantly—by several $\simeq 7\,\mathrm{km\,s^{-1}}$2—from that in the immediately surrounding disk (Udagawa et al., 27 Aug 2025). The reported anomaly is a coherent dip in $\simeq 7\,\mathrm{km\,s^{-1}}$3 of $\simeq 7\,\mathrm{km\,s^{-1}}$4 over a region $\simeq 7\,\mathrm{km\,s^{-1}}$5 in the outer disk near $\simeq 7\,\mathrm{km\,s^{-1}}$6, $\simeq 7\,\mathrm{km\,s^{-1}}$7, and $\simeq 7\,\mathrm{km\,s^{-1}}$8–$\simeq 7\,\mathrm{km\,s^{-1}}$9 (Udagawa et al., 27 Aug 2025). Prior to this result, large-scale vertical kinematic disturbances in the Milky Way disk had chiefly been discussed in connection with massive perturbers such as the Sagittarius dwarf spheroidal galaxy, which produce global warps, corrugations, or breathing modes. The reported VVA differs by being highly localized both spatially and kinematically (Udagawa et al., 27 Aug 2025).

The measurement used Gaia DR2 positions, parallaxes, and proper motions. Radial velocities were available for a subset but were not used in the principal VVA measurement; the analysis relied on transverse motions (Udagawa et al., 27 Aug 2025). Starting from Galactic latitude $(\ell,b,D)\simeq(15.8^\circ,-0.5^\circ,3.9$0, longitude $(\ell,b,D)\simeq(15.8^\circ,-0.5^\circ,3.9$1, parallax $(\ell,b,D)\simeq(15.8^\circ,-0.5^\circ,3.9$2, and proper motion $(\ell,b,D)\simeq(15.8^\circ,-0.5^\circ,3.9$3, the vertical component was written as

$(\ell,b,D)\simeq(15.8^\circ,-0.5^\circ,3.9$4

with $(\ell,b,D)\simeq(15.8^\circ,-0.5^\circ,3.9$5 and $(\ell,b,D)\simeq(15.8^\circ,-0.5^\circ,3.9$6 converting $(\ell,b,D)\simeq(15.8^\circ,-0.5^\circ,3.9$7 at $(\ell,b,D)\simeq(15.8^\circ,-0.5^\circ,3.9$8 to $(\ell,b,D)\simeq(15.8^\circ,-0.5^\circ,3.9$9 (Udagawa et al., 27 Aug 2025). In this study, the $4.0\,\mathrm{kpc})$0 term dominated, and the $4.0\,\mathrm{kpc})$1 contribution was omitted because radial-velocity coverage was incomplete.

To isolate noncircular motions, the systemic LSR transverse velocity $4.0\,\mathrm{kpc})$2 was subtracted from the longitudinal velocity $4.0\,\mathrm{kpc})$3, with $4.0\,\mathrm{kpc})$4, $4.0\,\mathrm{kpc})$5, $4.0\,\mathrm{kpc})$6, and $4.0\,\mathrm{kpc})$7 for $4.0\,\mathrm{kpc})$8 (Udagawa et al., 27 Aug 2025). Stars in the region $4.0\,\mathrm{kpc})$9, $$\delta V_z(R,\phi,z)\equiv\langle V_z\rangle(R,\phi,z)-\langle V_z\rangle_{\rm axisym}(R,z)$$0, $$\delta V_z(R,\phi,z)\equiv\langle V_z\rangle(R,\phi,z)-\langle V_z\rangle_{\rm axisym}(R,z)$$1 were gridded into voxels of size $$\delta V_z(R,\phi,z)\equiv\langle V_z\rangle(R,\phi,z)-\langle V_z\rangle_{\rm axisym}(R,z)$$2; within each voxel the median $$\delta V_z(R,\phi,z)\equiv\langle V_z\rangle(R,\phi,z)-\langle V_z\rangle_{\rm axisym}(R,z)$$3 was computed, and voxels containing fewer than five stars were discarded (Udagawa et al., 27 Aug 2025). A one-sample z-test comparing the candidate region to an adjacent background yielded $$\delta V_z(R,\phi,z)\equiv\langle V_z\rangle(R,\phi,z)-\langle V_z\rangle_{\rm axisym}(R,z)$$4, corresponding to a $$\delta V_z(R,\phi,z)\equiv\langle V_z\rangle(R,\phi,z)-\langle V_z\rangle_{\rm axisym}(R,z)$$5 detection (Udagawa et al., 27 Aug 2025).

The anomaly is centered at $$\delta V_z(R,\phi,z)\equiv\langle V_z\rangle(R,\phi,z)-\langle V_z\rangle_{\rm axisym}(R,z)$$6–$$\delta V_z(R,\phi,z)\equiv\langle V_z\rangle(R,\phi,z)-\langle V_z\rangle_{\rm axisym}(R,z)$$7 and is therefore a small-scale perturbation embedded within the broader vertical-kinematic structure of the Galactic disk (Udagawa et al., 27 Aug 2025). This localization is central to its interpretation: it is not a manifestation of the known large-scale warp or solar-neighbourhood bending signatures, but a compact disturbance co-located with specific gaseous features.

3. Association with gaseous disturbances and dark-subhalo interpretation

The stellar VVA reported in the Milky Way disk is spatially coincident with a suite of gaseous disturbances previously interpreted as signatures of a dark matter subhalo collision (Udagawa et al., 27 Aug 2025). The anomaly lies at the root of an H I filament seen in the HI4PI survey: a narrow, vertical feature $$\delta V_z(R,\phi,z)\equiv\langle V_z\rangle(R,\phi,z)-\langle V_z\rangle_{\rm axisym}(R,z)$$8 long in $$\delta V_z(R,\phi,z)\equiv\langle V_z\rangle(R,\phi,z)-\langle V_z\rangle_{\rm axisym}(R,z)$$9 and spanning LSR velocities $FF_w\equiv \langle w'^4\rangle/\sigma_w^4$0–$FF_w\equiv \langle w'^4\rangle/\sigma_w^4$1, with a central void of diameter $FF_w\equiv \langle w'^4\rangle/\sigma_w^4$2 (Udagawa et al., 27 Aug 2025). A molecular shell, the CO $FF_w\equiv \langle w'^4\rangle/\sigma_w^4$3–$FF_w\equiv \langle w'^4\rangle/\sigma_w^4$4 feature “CO 16.134–0.553,” of radius $FF_w\equiv \langle w'^4\rangle/\sigma_w^4$5 and velocity width $FF_w\equiv \langle w'^4\rangle/\sigma_w^4$6 forms the eastern edge of the H I void (Udagawa et al., 27 Aug 2025). In the paper’s overlays, the stellar $FF_w\equiv \langle w'^4\rangle/\sigma_w^4$7 map shows near-perfect positional coincidence between the VVA, the H I void, the CO shell, and the H I filament (Udagawa et al., 27 Aug 2025).

The proposed dynamical origin is an ultra-compact dwarf (UCD)-sized dark matter subhalo plunging through the disk at $FF_w\equiv \langle w'^4\rangle/\sigma_w^4$8 (Udagawa et al., 27 Aug 2025). Under the impulse approximation for a point-mass perturber of mass $FF_w\equiv \langle w'^4\rangle/\sigma_w^4$9 and impact parameter $FF_w>3$0, the velocity kick is

$FF_w>3$1

Two mass constraints were given. A tidal-stability criterion,

$FF_w>3$2

implies $FF_w>3$3 for $FF_w>3$4 and $FF_w>3$5 (Udagawa et al., 27 Aug 2025). A second estimate based on the Hoyle–Lyttleton radius,

$FF_w>3$6

yields $FF_w>3$7 for $FF_w>3$8 and $FF_w>3$9 (Udagawa et al., 27 Aug 2025). The absence of a luminous source at the filament’s terminus was therefore interpreted as evidence for a dark or failed UCD (Udagawa et al., 27 Aug 2025).

Within the paper’s cosmological framing, a dark subhalo of $v_z$0–$v_z$1 is fully consistent with $v_z$2CDM expectations for the subhalo mass function around Milky-Way-mass halos, $v_z$3 with $v_z$4 (Udagawa et al., 27 Aug 2025). The cited theoretical context includes high-resolution N-body simulations predicting $v_z$5–$v_z$6 subhalos in the mass range $v_z$7–$v_z$8 within a Milky-Way-type virial radius (Udagawa et al., 27 Aug 2025). In that sense, the reported VVA is presented as one of the very few direct dynamical detections at this low-mass scale.

4. Related astrophysical usages: Galactic mapping and helioseismology

A distinct Galactic usage appears in large-area stellar-disk mapping with Gaia DR3 and StarHorse distances, where the anomaly is defined as

$v_z$9

In the absence of bulk vertical streaming, $$\mathrm{VVA}(x,y,z)\equiv w_{\rm meas}(x,y,z)-w_{\rm pred}(x,y,z)$$0, so $$\mathrm{VVA}(x,y,z)\equiv w_{\rm meas}(x,y,z)-w_{\rm pred}(x,y,z)$$1 is the measured $$\mathrm{VVA}(x,y,z)\equiv w_{\rm meas}(x,y,z)-w_{\rm pred}(x,y,z)$$2 after subtraction of the Sun’s reflex motion, taken as $$\mathrm{VVA}(x,y,z)\equiv w_{\rm meas}(x,y,z)-w_{\rm pred}(x,y,z)$$3, and any smooth warp-induced gradient (Widmark et al., 2022). The analysis used stars with radial-velocity errors $$\mathrm{VVA}(x,y,z)\equiv w_{\rm meas}(x,y,z)-w_{\rm pred}(x,y,z)$$4, transformed them into Galactocentric $$\mathrm{VVA}(x,y,z)\equiv w_{\rm meas}(x,y,z)-w_{\rm pred}(x,y,z)$$5 with $$\mathrm{VVA}(x,y,z)\equiv w_{\rm meas}(x,y,z)-w_{\rm pred}(x,y,z)$$6 and $$\mathrm{VVA}(x,y,z)\equiv w_{\rm meas}(x,y,z)-w_{\rm pred}(x,y,z)$$7, and computed $$\mathrm{VVA}(x,y,z)\equiv w_{\rm meas}(x,y,z)-w_{\rm pred}(x,y,z)$$8 in spatial bins with $$\mathrm{VVA}(x,y,z)\equiv w_{\rm meas}(x,y,z)-w_{\rm pred}(x,y,z)$$9 and variable $\mathrm{VVA}(x,y,z,t)=w(x,y,z,t)-\overline{w}(z)$0, requiring $\mathrm{VVA}(x,y,z,t)=w(x,y,z,t)-\overline{w}(z)$1 effective stars per cell (Widmark et al., 2022). Two coherent anomalies emerged across all four magnitude bins: a Local-Spiral-Arm breathing mode centered near $\mathrm{VVA}(x,y,z,t)=w(x,y,z,t)-\overline{w}(z)$2 with peak $\mathrm{VVA}(x,y,z,t)=w(x,y,z,t)-\overline{w}(z)$3, peak $\mathrm{VVA}(x,y,z,t)=w(x,y,z,t)-\overline{w}(z)$4, wavelength $\mathrm{VVA}(x,y,z,t)=w(x,y,z,t)-\overline{w}(z)$5, and a $\mathrm{VVA}(x,y,z,t)=w(x,y,z,t)-\overline{w}(z)$6 phase shift between density and vertical motion; and a large-scale radial-gradient mode at $\mathrm{VVA}(x,y,z,t)=w(x,y,z,t)-\overline{w}(z)$7 with $\mathrm{VVA}(x,y,z,t)=w(x,y,z,t)-\overline{w}(z)$8 in $\mathrm{VVA}(x,y,z,t)=w(x,y,z,t)-\overline{w}(z)$9 and $\simeq 7\,\mathrm{km\,s^{-1}}$00 in density asymmetry (Widmark et al., 2022). Relative to the localized VVA of Udagawa et al., these are extended perturbative modes rather than a compact impact signature.

In helioseismology, the term denotes anomalous vertical plasma flow beneath an average solar supergranule rather than a stellar-kinematic perturbation in the Milky Way (Svanda, 2012). Using HMI Dopplergrams, finite-frequency kernels under the Born approximation, and SOLA inversion, Švanda constructed an average over $\simeq 7\,\mathrm{km\,s^{-1}}$01 independent supergranules (Svanda, 2012). The inversion targeted $\simeq 7\,\mathrm{km\,s^{-1}}$02 through a weighted combination of travel-time measurements,

$\simeq 7\,\mathrm{km\,s^{-1}}$03

with trade-off parameters controlling localization, noise, and cross-talk (Svanda, 2012). A good-localisation solution with single-supergranule noise of $\simeq 7\,\mathrm{km\,s^{-1}}$04 was combined with statistical averaging, reducing the random-noise level by $\simeq 7\,\mathrm{km\,s^{-1}}$05 and yielding an uncertainty $\simeq 7\,\mathrm{km\,s^{-1}}$06 in the final average (Svanda, 2012). The resulting VVA peaked at

$\simeq 7\,\mathrm{km\,s^{-1}}$07

with root-mean-square vertical velocity $\simeq 7\,\mathrm{km\,s^{-1}}$08 over a cell radius $\simeq 7\,\mathrm{km\,s^{-1}}$09, whereas averaged line-of-sight Dopplergrams gave a photospheric upflow of only $\simeq 7\,\mathrm{km\,s^{-1}}$10 (Svanda, 2012). The discrepancy was interpreted as evidence that supergranular upflows are concentrated in a narrow subsurface layer and are smeared out by inversions that favor low noise over localization (Svanda, 2012).

5. Fluid-dynamical, oceanic, and atmospheric formulations

In high-Reynolds-number wall turbulence, Banerjee et al. quantify VVA by the flatness factor of wall-normal velocity fluctuations,

$\simeq 7\,\mathrm{km\,s^{-1}}$11

with $\simeq 7\,\mathrm{km\,s^{-1}}$12 indicating heavy tails and large-scale intermittency (Banerjee et al., 19 Jun 2026). Starting from the wall-normal Navier–Stokes fluctuation equation, they argue that in the inertial sublayer the dominant balance is between the pure-inertial divergence of the fifth moment and the pressure–velocity interaction, while the mixed term $\simeq 7\,\mathrm{km\,s^{-1}}$13 is negligible (Banerjee et al., 19 Jun 2026). Using an extended Rotta closure, a quasi-Gaussian approximation, and a macro-scale choice $\simeq 7\,\mathrm{km\,s^{-1}}$14, they derive

$\simeq 7\,\mathrm{km\,s^{-1}}$15

where $\simeq 7\,\mathrm{km\,s^{-1}}$16 and $\simeq 7\,\mathrm{km\,s^{-1}}$17 (Banerjee et al., 19 Jun 2026). Empirically, thirteen laboratory experiments spanning $\simeq 7\,\mathrm{km\,s^{-1}}$18–$\simeq 7\,\mathrm{km\,s^{-1}}$19 and five near-neutral atmospheric surface-layer sites collapsed to a nearly constant plateau $\simeq 7\,\mathrm{km\,s^{-1}}$20 for $\simeq 7\,\mathrm{km\,s^{-1}}$21, with field and lab inversions clustering around $\simeq 7\,\mathrm{km\,s^{-1}}$22 (Banerjee et al., 19 Jun 2026). In this literature, the anomaly is a fourth-order intermittency diagnostic rather than a mean-velocity perturbation.

In quasigeostrophic laboratory vortices, the anomaly is defined pointwise as the discrepancy between measured and predicted vertical velocity:

$\simeq 7\,\mathrm{km\,s^{-1}}$23

The classical QG $\simeq 7\,\mathrm{km\,s^{-1}}$24-equation predicted peak values $\simeq 7\,\mathrm{km\,s^{-1}}$25 ($\simeq 7\,\mathrm{km\,s^{-1}}$26) for a Gaussian vortex with $\simeq 7\,\mathrm{km\,s^{-1}}$27, $\simeq 7\,\mathrm{km\,s^{-1}}$28, and $\simeq 7\,\mathrm{km\,s^{-1}}$29, with multipolar lobes concentrated at the vortex periphery (Aulnette et al., 19 Dec 2025). However, independent measurements from horizontal-divergence integration and particle residence times both gave $\simeq 7\,\mathrm{km\,s^{-1}}$30, about five times larger (Aulnette et al., 19 Dec 2025). The extended $\simeq 7\,\mathrm{km\,s^{-1}}$31-equation,

$\simeq 7\,\mathrm{km\,s^{-1}}$32

was then used to test whether viscous and scalar diffusion could account for part of the missing vertical motion (Aulnette et al., 19 Dec 2025).

A different oceanographic usage appears in the VISION reconstruction framework, where the anomaly is defined relative to a climatological mean:

$\simeq 7\,\mathrm{km\,s^{-1}}$33

Here the target field is low-pass-filtered vertical velocity $\simeq 7\,\mathrm{km\,s^{-1}}$34 from the KD48 benchmark, a $\simeq 7\,\mathrm{km\,s^{-1}}$35 regional Kuroshio Extension subset of LLC4320 MITgcm with hourly snapshots over $\simeq 7\,\mathrm{km\,s^{-1}}$36 days and surface variables SSH, buoyancy, and zonal and meridional surface velocities (Gao et al., 25 Sep 2025). VISION uses a Dynamic Prompting paradigm with an availability mask, a Universal Observation Adapter, a State-conditioned Prompting module, and geometry- and scale-aware operators, and is trained with smooth-$\simeq 7\,\mathrm{km\,s^{-1}}$37 loss on a denoised target $\simeq 7\,\mathrm{km\,s^{-1}}$38 (Gao et al., 25 Sep 2025). On the reported benchmark, VISION achieved RMSE values of $\simeq 7\,\mathrm{km\,s^{-1}}$39, $\simeq 7\,\mathrm{km\,s^{-1}}$40, and $\simeq 7\,\mathrm{km\,s^{-1}}$41 (all $\simeq 7\,\mathrm{km\,s^{-1}}$42) at depths $\simeq 7\,\mathrm{km\,s^{-1}}$43, $\simeq 7\,\mathrm{km\,s^{-1}}$44, and $\simeq 7\,\mathrm{km\,s^{-1}}$45, with corresponding PCC values $\simeq 7\,\mathrm{km\,s^{-1}}$46, $\simeq 7\,\mathrm{km\,s^{-1}}$47, and $\simeq 7\,\mathrm{km\,s^{-1}}$48 (Gao et al., 25 Sep 2025). In this framework, VVA is the reconstructed instantaneous departure from a depth-dependent climatological mean.

In a steady linear $\simeq 7\,\mathrm{km\,s^{-1}}$49-layer atmospheric model on an $\simeq 7\,\mathrm{km\,s^{-1}}$50-plane, the relevant anomaly concept concerns the vertical structure of $\simeq 7\,\mathrm{km\,s^{-1}}$51 profiles. The model admits free modes satisfying

$\simeq 7\,\mathrm{km\,s^{-1}}$52

whose vertical structures resemble the first and second baroclinic modes (Ahmed et al., 18 Sep 2025). In the standard parameter regime, the first-baroclinic mode has a characteristic horizontal scale of $\simeq 7\,\mathrm{km\,s^{-1}}$53 and the second-baroclinic mode a smaller scale of $\simeq 7\,\mathrm{km\,s^{-1}}$54 (Ahmed et al., 18 Sep 2025). Strong-gradient surface-temperature forcing projects strongly onto the second mode and yields bottom-heavy $\simeq 7\,\mathrm{km\,s^{-1}}$55 profiles, whereas weak-gradient forcing projects strongly onto the first mode and yields top-heavy profiles (Ahmed et al., 18 Sep 2025). In that usage, “vertical velocity anomalies” refer to top- and bottom-heavy variants of tropical atmospheric ascent rather than local residuals with respect to a background field.

6. Extreme drafts in the mesosphere and lower thermosphere, and recurrent conceptual issues

In the mesosphere and lower thermosphere (MLT), VVA is identified with “Rogue Vertical Drafts” (RVDs), defined by two simultaneous criteria: $\simeq 7\,\mathrm{km\,s^{-1}}$56 and exceedance of five standard deviations of the local vertical-velocity distribution (Chau et al., 10 Dec 2025). The anomaly parameter is

$\simeq 7\,\mathrm{km\,s^{-1}}$57

so rogue events satisfy $\simeq 7\,\mathrm{km\,s^{-1}}$58 (Chau et al., 10 Dec 2025). In multi-year summer observations over Northern Norway, the vertical-velocity distribution was approximated by a Gaussian with $\simeq 7\,\mathrm{km\,s^{-1}}$59 and $\simeq 7\,\mathrm{km\,s^{-1}}$60, giving a five-sigma threshold of about $\simeq 7\,\mathrm{km\,s^{-1}}$61 (Chau et al., 10 Dec 2025). Assuming independent samples every $\simeq 7\,\mathrm{km\,s^{-1}}$62, the expected waiting time between extreme events is modified by the fact that RVDs typically span $\simeq 7\,\mathrm{km\,s^{-1}}$63 of consecutive samples, yielding an average recurrence of about one event every $\simeq 7\,\mathrm{km\,s^{-1}}$64 days in summer over Northern Norway (Chau et al., 10 Dec 2025). Volumetric radar imaging further indicates varicose updraft–downdraft pairs with vertical extent $\simeq 7\,\mathrm{km\,s^{-1}}$65–$\simeq 7\,\mathrm{km\,s^{-1}}$66, horizontal widths $\simeq 7\,\mathrm{km\,s^{-1}}$67–$\simeq 7\,\mathrm{km\,s^{-1}}$68, occasional elongation to $\simeq 7\,\mathrm{km\,s^{-1}}$69, and durations from a few minutes to tens of minutes (Chau et al., 10 Dec 2025).

The broader literature reveals several recurring issues. First, the anomaly baseline is field-specific: local background in the Galactic-disk impact candidate (Udagawa et al., 27 Aug 2025), axisymmetric symmetry assumptions in Gaia DR3 disk mapping (Widmark et al., 2022), zero mean in supergranular inversions (Svanda, 2012), model-predicted $\simeq 7\,\mathrm{km\,s^{-1}}$70 in quasigeostrophic vortices (Aulnette et al., 19 Dec 2025), climatology in data-driven ocean reconstruction (Gao et al., 25 Sep 2025), and Gaussian reference statistics in MLT extremes (Chau et al., 10 Dec 2025). Second, multiple papers frame VVA through model–data mismatch. In wall turbulence, a down-gradient closure predicts $\simeq 7\,\mathrm{km\,s^{-1}}$71 when skewness is constant and therefore fails to reproduce the observed inertial-sublayer plateau, whereas the extended Rotta–QGA approach recovers $\simeq 7\,\mathrm{km\,s^{-1}}$72 and the observed weak decline and outer-layer upturn (Banerjee et al., 19 Jun 2026). In quasigeostrophic vortices, the classical $\simeq 7\,\mathrm{km\,s^{-1}}$73-equation underestimates measured vertical velocity by a factor of about five, motivating inclusion of dissipative terms (Aulnette et al., 19 Dec 2025). In helioseismology, the tension is between a subsurface updraft of $\simeq 7\,\mathrm{km\,s^{-1}}$74 and a photospheric signal of only $\simeq 7\,\mathrm{km\,s^{-1}}$75, attributed to localization–noise trade-offs in inversion kernels (Svanda, 2012).

A plausible implication is that “vertical velocity anomaly” is best understood as a methodological category for statistically or dynamically exceptional vertical motion, with the exact quantity determined by the governing equations, measurement geometry, and null model of the field under study. In Galactic dynamics, the term currently has particular prominence because the localized Milky Way VVA has been interpreted as rare observational evidence for a low-mass dark matter subhalo or dark, ultra-compact dwarf-sized intruder, linking stellar kinematics, gas morphology, and $\simeq 7\,\mathrm{km\,s^{-1}}$76CDM substructure directly (Udagawa et al., 27 Aug 2025).