Low-Cost High-Order Singular Value Decomposition for Tensor-Based Reconstruction from Sparse Sensor Measurements: Urban Flow and Air-Quality Applications
Reviewed by Pith T0 review T1 audit T2 compute T3 formal T4 reserved 2026-06-26 00:18 UTCgrok-4.3pith:3OQASMZPrecord.jsonopen to challenge →
The pith
lcHOSVD reconstructs urban velocity and pollutant fields from 1-4% sparse sensors by preserving tensor structure instead of flattening data into matrices.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
lcHOSVD is the first method to combine sparse sensing with HOSVD for field reconstruction. It preserves the natural tensor structure of high-dimensional velocity and concentration data, exploits correlations across spatial, temporal, and physical-variable dimensions, and reduces the computational cost of conventional HOSVD while delivering lower reconstruction errors than low-cost matrix SVD on urban flow and air-quality cases observed at only 1-4% of spatial locations; it is also more robust to anisotropic sensor distributions.
What carries the argument
lcHOSVD, the low-cost approximation to high-order singular value decomposition that operates on the full tensor to retain cross-dimensional correlations during sparse reconstruction.
If this is right
- lcHOSVD yields lower reconstruction errors than lcSVD on data with strong multidimensional coupling and heterogeneous dynamics.
- The tensor formulation remains accurate under uneven sensor placements that commonly occur in real monitoring networks.
- Computational cost is substantially lower than conventional HOSVD while retaining the ability to use tensor correlations.
- Reconstruction quality improves when spatial, temporal, and variable dimensions are kept coupled rather than flattened.
Where Pith is reading between the lines
- The same tensor-preserving approach could be tested on other sparse high-dimensional datasets such as climate or ocean models to check whether error reductions generalize.
- Robustness to anisotropic sampling suggests the method may improve data assimilation pipelines that ingest irregular sensor networks.
- Embedding lcHOSVD inside real-time digital-twin systems would allow direct comparison of reconstruction latency and accuracy against current matrix pipelines.
Load-bearing premise
The low-cost approximation to full HOSVD still captures enough cross-dimensional correlations that its error advantage over matrix methods is preserved on these urban datasets at 1-4% sensor coverage.
What would settle it
Apply both lcHOSVD and lcSVD to the urban flow dataset with 2% randomly chosen sensor locations and measure average reconstruction error; if the tensor method does not produce lower error, the performance claim is falsified.
Figures
read the original abstract
Urban flow and air-quality simulations generate high-dimensional datasets describing velocity and pollutant transport across multiple spatial, temporal, and physical-variable dimensions. Reconstructing these fields from sparse sensor measurements is a fundamental challenge in environmental monitoring, digital twins, forecasting, and data assimilation. Existing low-cost reconstruction approaches are commonly based on matrix decompositions, which require multidimensional datasets to be flattened into two-dimensional snapshot matrices, thereby discarding important structural information. This work introduces the low-cost High-Order Singular Value Decomposition (lcHOSVD), a novel tensor-based sparse-sensing reconstruction framework for high-dimensional environmental fields. To the authors' knowledge, this is the first methodology that combines sparse sensing and HOSVD for field reconstruction. Unlike matrix-based approaches, lcHOSVD preserves the natural tensor structure of the data, enabling the exploitation of correlations across spatial, temporal, and physical-variable dimensions while substantially reducing the computational requirements of conventional HOSVD. The methodology is applied to urban flow and air-quality datasets, where three-dimensional velocity and pollutant concentration fields are reconstructed using only 1-4% of the available spatial locations. While lcSVD provides larger computational speed-ups, lcHOSVD consistently achieves lower reconstruction errors in configurations characterized by strong multidimensional coupling and heterogeneous dynamics across dimensions. Additional sensor-anisotropy analyses demonstrate that the tensor formulation is significantly more robust to uneven sensor distributions, a common situation in practical environmental monitoring networks.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces lcHOSVD, a low-cost approximation to high-order singular value decomposition, as a tensor-based framework for reconstructing 3D velocity and pollutant fields from 1-4% sparse sensor measurements in urban flow and air-quality datasets. It claims to be the first combination of sparse sensing with HOSVD, preserving multi-way correlations across spatial, temporal, and variable modes better than matrix-based lcSVD, yielding lower reconstruction errors under heterogeneous dynamics and greater robustness to anisotropic sensor placement.
Significance. If the performance claims hold with rigorous validation, the work could meaningfully advance sparse reconstruction techniques in environmental fluid dynamics by demonstrating practical benefits of tensor structure preservation at reduced computational cost relative to full HOSVD, with direct relevance to digital twins and data assimilation.
major comments (2)
- [Abstract] Abstract: the central claim that lcHOSVD 'consistently achieves lower reconstruction errors' than lcSVD in strong multidimensional coupling cases lacks any supporting error metrics, cross-validation procedure, or dataset description in the provided text, preventing assessment of whether the low-cost approximation retains the necessary cross-mode correlations.
- [Abstract] Abstract: no derivation, error bound, or condition is given for when the low-cost HOSVD approximation begins to lose the cross-dimensional correlation exploitation that is asserted to drive the advantage over lcSVD; this is load-bearing for the claimed superiority at 1-4% coverage.
minor comments (1)
- [Abstract] Abstract: the phrase 'to the authors' knowledge' should be replaced by a precise literature search statement or removed.
Simulated Author's Rebuttal
We thank the referee for their comments on the abstract. The full manuscript contains the supporting details in Sections 2–4, but we agree the abstract would benefit from greater self-containment. We respond to each point below.
read point-by-point responses
-
Referee: [Abstract] Abstract: the central claim that lcHOSVD 'consistently achieves lower reconstruction errors' than lcSVD in strong multidimensional coupling cases lacks any supporting error metrics, cross-validation procedure, or dataset description in the provided text, preventing assessment of whether the low-cost approximation retains the necessary cross-mode correlations.
Authors: The error metrics, cross-validation procedure, and dataset descriptions appear in the full manuscript (Sections 2, 3, and 4). To address the concern, we will revise the abstract to include a concise reference to the quantitative improvements demonstrated in the experiments and the validation approach employed. revision: yes
-
Referee: [Abstract] Abstract: no derivation, error bound, or condition is given for when the low-cost HOSVD approximation begins to lose the cross-dimensional correlation exploitation that is asserted to drive the advantage over lcSVD; this is load-bearing for the claimed superiority at 1-4% coverage.
Authors: The derivation of the lcHOSVD low-cost approximation is given in Section 3.1. No theoretical error bound is provided; the conditions for retaining cross-mode correlations are instead characterized empirically in Section 4 across the 1–4% coverage range. We will revise the abstract to note that the reported advantage holds under the strong multidimensional coupling regimes examined in the numerical experiments. revision: partial
Circularity Check
No significant circularity detected
full rationale
The paper presents lcHOSVD as a novel tensor-based method for sparse reconstruction, with claims centered on its introduction, computational advantages over conventional HOSVD, and empirical performance on urban flow/air-quality datasets (lower errors than lcSVD at 1-4% sensor coverage). No derivation chain, equations, or fitted parameters are described that reduce outputs to inputs by construction, nor are there load-bearing self-citations or uniqueness theorems invoked. The work is methodological and data-driven rather than deriving predictions from self-referential fits, making the central claims self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
Sustainability of Europe’s mobility sys- tems,
European Environment Agency, “Sustainability of Europe’s mobility sys- tems,” tech. rep., European Environment Agency, Copenhagen, Denmark,
-
[2]
Available at:https://www.eea.europa.eu/en/analysis/publications/ sustainability-of-europes-mobility-systems
-
[3]
Europe’s Air Quality Status 2024,
European Environment Agency, “Europe’s Air Quality Status 2024,” tech. rep., European Environment Agency, Copenhagen, Denmark, 2024. Available at:https://www.eea.europa.eu/en/analysis/publications/ europes-air-quality-status-2024
2024
-
[4]
Computational fluid dynamics for urban physics: Importance, scales, possibilities, limitations and ten tips and tricks towards accurate and reliable simu- lations,
B. Blocken, “Computational fluid dynamics for urban physics: Importance, scales, possibilities, limitations and ten tips and tricks towards accurate and reliable simu- lations,”Building and environment, vol. 91, pp. 219–245, 2015
2015
-
[5]
Ten questions concerning modeling of near-field pollutant dispersion in the built environment,
Y. Tominaga and T. Stathopoulos, “Ten questions concerning modeling of near-field pollutant dispersion in the built environment,”Building and Environment, vol. 105, pp. 390–402, 2016
2016
-
[6]
Particle image velocimetry for combustion measurements: Applications and developments,
C. Fang and L. Hong, “Particle image velocimetry for combustion measurements: Applications and developments,”Chinese Journal of Aeronautics, vol. 31, no. 7, pp. 1407–1427, 2018
2018
-
[7]
Multiscale proper ortho- gonal decomposition (mpod) of tr-piv data—a case study on stationary and transient cylinder wake flows,
M. A. Mendez, D. Hess, B. B. Watz, and J.-M. Buchlin, “Multiscale proper ortho- gonal decomposition (mpod) of tr-piv data—a case study on stationary and transient cylinder wake flows,”Measurement Science and Technology, vol. 31, no. 9, p. 094014, 2020
2020
-
[8]
Importance of the nozzle-exit boundary-layer state in subsonic turbulent jets,
G. A. Brès, P. Jordan, V. Jaunet, M. Le Rallic, A. V. Cavalieri, A. Towne, S. K. Lele, T. Colonius, and O. T. Schmidt, “Importance of the nozzle-exit boundary-layer state in subsonic turbulent jets,”Journal of Fluid Mechanics, vol. 851, pp. 83–124, 2018
2018
-
[9]
Direct numerical simulation of non-premixed turbulent flames,
L. Vervisch and T. Poinsot, “Direct numerical simulation of non-premixed turbulent flames,”Annual review of fluid mechanics, vol. 30, no. 1, pp. 655–691, 1998
1998
-
[10]
Spatio-temporal flow structures in the three-dimensional wake of a circular cylinder,
S. L. Clainche, J. M. Pérez, and J. M. Vega, “Spatio-temporal flow structures in the three-dimensional wake of a circular cylinder,”Fluid Dynamics Research, vol. 50, no. 5, p. 051406, 2018
2018
-
[11]
Singular value decomposition and least squares solu- tions,
G. H. Golub and C. Reinsch, “Singular value decomposition and least squares solu- tions,” inLinear algebra, pp. 134–151, Springer, 1971
1971
-
[12]
The structure of inhomogeneous turbulent flows,
J. L. Lumley, “The structure of inhomogeneous turbulent flows,”Atmospheric turbu- lence and radio wave propagation, pp. 166–178, 1967. 32
1967
-
[13]
The proper orthogonal decomposition in the analysis of turbulent flows,
G. Berkooz, P. Holmes, and J. L. Lumley, “The proper orthogonal decomposition in the analysis of turbulent flows,”Annual review of fluid mechanics, vol. 25, no. 1, pp. 539–575, 1993
1993
-
[14]
Holmes,Turbulence, coherent structures, dynamical systems and symmetry
P. Holmes,Turbulence, coherent structures, dynamical systems and symmetry. Cam- bridge university press, 2012
2012
-
[15]
A predictive hybrid reduced order model based on proper ortho- gonal decomposition combined with deep learning architectures,
R. Abadía-Heredia, M. López-Martín, B. Carro, J. I. Arribas, J. M. Pérez, and S. Le Clainche, “A predictive hybrid reduced order model based on proper ortho- gonal decomposition combined with deep learning architectures,”Expert Systems with Applications, vol. 187, p. 115910, 2022
2022
-
[16]
Proper orthogonal decomposition of large-eddy simulation data over real urban morphology,
Y. Liu, C.-H. Liu, G. P. Brasseur, and C. Y. Chao, “Proper orthogonal decomposition of large-eddy simulation data over real urban morphology,”Sustainable Cities and Society, vol. 89, p. 104324, 2023
2023
-
[17]
Non-intrusive reduced order model of urban airflow with dynamic boundary conditions,
S. Xiang, X. Fu, J. Zhou, Y. Wang, Y. Zhang, X. Hu, J. Xu, H. Liu, J. Liu, J. Ma, et al., “Non-intrusive reduced order model of urban airflow with dynamic boundary conditions,”Building and Environment, vol. 187, p. 107397, 2021
2021
-
[18]
A reduced order model for turbulent flows in the urban environment using machine learning,
D. Xiao, C. Heaney, L. Mottet, F. Fang, W. Lin, I. Navon, Y. Guo, O. Matar, A. Robins, and C. Pain, “A reduced order model for turbulent flows in the urban environment using machine learning,”Building and Environment, vol. 148, pp. 323– 337, 2019
2019
-
[19]
Dynamic mode decomposition of numerical and experimental data,
P. J. Schmid, “Dynamic mode decomposition of numerical and experimental data,” Journal of fluid mechanics, vol. 656, pp. 5–28, 2010
2010
-
[20]
Reduced-order modelling of urban wind environment and gaseous pollutants dispersion in an urban-scale street canyon,
S. Ding and R. Yang, “Reduced-order modelling of urban wind environment and gaseous pollutants dispersion in an urban-scale street canyon,”Journal of Safety Science and Resilience, vol. 2, no. 4, pp. 238–245, 2021
2021
-
[21]
Stable reduced-order models for pollutant dispersion in the built environment,
L. Vervecken, J. Camps, and J. Meyers, “Stable reduced-order models for pollutant dispersion in the built environment,”Building and Environment, vol. 92, pp. 360–367, 2015
2015
-
[22]
Llm-rom: A novel framework for efficient spatiotem- poral prediction of urban pollutant dispersion,
P. Wu, Z. Qin, and Y. Yang, “Llm-rom: A novel framework for efficient spatiotem- poral prediction of urban pollutant dispersion,”AI, vol. 7, no. 3, p. 104, 2026
2026
-
[23]
A review of advances towards effi- cient reduced-order models (rom) for predicting urban airflow and pollutant disper- sion,
S. Masoumi-Verki, F. Haghighat, and U. Eicker, “A review of advances towards effi- cient reduced-order models (rom) for predicting urban airflow and pollutant disper- sion,”Building and Environment, vol. 216, p. 108966, 2022
2022
-
[24]
A multilinear singular value decomposition,
L. De Lathauwer, B. De Moor, and J. Vandewalle, “A multilinear singular value decomposition,”SIAM journal on Matrix Analysis and Applications, vol. 21, no. 4, pp. 1253–1278, 2000. 33
2000
-
[25]
On the best rank-1 and rank-(r 1, r 2,..., rn) approximation of higher-order tensors,
L. De Lathauwer, B. De Moor, and J. Vandewalle, “On the best rank-1 and rank-(r 1, r 2,..., rn) approximation of higher-order tensors,”SIAM journal on Matrix Analysis and Applications, vol. 21, no. 4, pp. 1324–1342, 2000
2000
-
[26]
Tensor decompositions and applications,
T. G. Kolda and B. W. Bader, “Tensor decompositions and applications,”SIAM review, vol. 51, no. 3, pp. 455–500, 2009
2009
-
[27]
Higher order dynamic mode decomposition,
S. Le Clainche and J. M. Vega, “Higher order dynamic mode decomposition,”SIAM Journal on Applied Dynamical Systems, vol. 16, no. 2, pp. 882–925, 2017
2017
-
[28]
Low-cost singular value decomposition with optimal sensor placement,
A. Hetherington and S. Le Clainche, “Low-cost singular value decomposition with optimal sensor placement,”Physics of Fluids, vol. 37, no. 8, 2025
2025
-
[29]
Data-driven sparse sensor placement for reconstruction: Demonstrating the benefits of exploiting known patterns,
K. Manohar, B. W. Brunton, J. N. Kutz, and S. L. Brunton, “Data-driven sparse sensor placement for reconstruction: Demonstrating the benefits of exploiting known patterns,”IEEE Control Systems Magazine, vol. 38, no. 3, pp. 63–86, 2018
2018
-
[30]
Pysensors: A python package for sparse sensor placement,
B. M. De Silva, K. Manohar, E. Clark, B. W. Brunton, S. L. Brunton, and J. N. Kutz, “Pysensors: A python package for sparse sensor placement,”arXiv preprint arXiv:2102.13476, 2021
-
[31]
A hierarchical machine learning framework for real-time reconstruction of urban wind fields from sparse sensor networks,
F. M. Nav, R. Snaiki, and F. Ding, “A hierarchical machine learning framework for real-time reconstruction of urban wind fields from sparse sensor networks,”Building and Environment, p. 114875, 2026
2026
-
[32]
A. Vishwasrao, S. B. C. Gutha, A. Cremades, K. Wijk, A. Patil, C. Gorle, B. J. McKeon, H. Azizpour, and R. Vinuesa, “Diff-sport: Diffusion-based sensor placement optimization and reconstruction of turbulent flows in urban environments,”arXiv preprint arXiv:2506.00214, 2025
-
[33]
A low-cost singular value decomposition-based data assimilation technique for analysis of heterogeneous combustion data,
P. Pillai, A. I. Hetherington, L. Saavedra, and S. Le Clainche, “A low-cost singular value decomposition-based data assimilation technique for analysis of heterogeneous combustion data,”Physics of Fluids, vol. 37, no. 8, 2025
2025
-
[34]
P. Jeanney, A. Hetherington, S. E. Ahmed, D. Lanceta, S. Saiz, J. M. Perez, and S. L. Clainche, “Ensemble kalman filter for data assimilation coupled with low-resolution computations techniques applied in fluid dynamics,”arXiv preprint arXiv:2507.00539, 2025
-
[35]
Modelflows-app: Data-drivenpost- processing and reduced order modelling tools,
A. Hetherington, A. Corrochano, R. Abadía-Heredia, E. Lazpita, E. Muñoz, P. Díaz, E.Maiora, M.López-Martín, andS.LeClainche, “Modelflows-app: Data-drivenpost- processing and reduced order modelling tools,”Computer Physics Communications, vol. 301, p. 109217, 2024
2024
-
[36]
Normalization: A Preprocessing Stage
S. Patro and K. K. Sahu, “Normalization: A preprocessing stage,”arXiv preprint arXiv:1503.06462, 2015. 34
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[37]
Eddies, streams, and convergence zones in turbulent flows,
J. C. Hunt, A. A. Wray, and P. Moin, “Eddies, streams, and convergence zones in turbulent flows,”Studying turbulence using numerical simulation databases, 2. Proceedings of the 1988 summer program, 1988
1988
-
[38]
Large-scale cfd analysis of urban airflow and pollutant transport over a 24-hour period with pod analysis: The vallecas district (madrid) case study,
P. Jeanney, C. García-Sánchez, S. Saiz, J. M. Pérez, and S. Le Clainche, “Large-scale cfd analysis of urban airflow and pollutant transport over a 24-hour period with pod analysis: The vallecas district (madrid) case study,”arXiv preprint, 2026
2026
-
[39]
E. Economista, “La gran transformación dentro de Vallecas: de zona degradada a nuevo barrio con 1.400 viviendas y una residencia de estudiantes.”https: //www.eleconomista.es/vivienda-inmobiliario/noticias/13060987/10/24/ la-gran-transformacion-dentro-de-vallecas-de-zona-degradada-a-nuevo-barrio-con-1400-viviendas-y-una-residencia-de-estudiantes. html, 10, 2024
-
[40]
Towards automatic reconstruction of 3d city models tailored for urban flow simulations,
I. Pađen, C. García-Sánchez, and H. Ledoux, “Towards automatic reconstruction of 3d city models tailored for urban flow simulations,”Frontiers in Built Environment, vol. 8, p. 899332, 2022
2022
-
[41]
Geoportal – red de vigilancia de la calidad del aire
Ayuntamiento de Madrid, “Geoportal – red de vigilancia de la calidad del aire.” https://geoportal.madrid.es/IDEAM_WBGEOPORTAL/index.iam, 2025
2025
-
[42]
A tensorial approach to computa- tional continuum mechanics using object-oriented techniques,
H. G. Weller, G. Tabor, H. Jasak, and C. Fureby, “A tensorial approach to computa- tional continuum mechanics using object-oriented techniques,”Computers in physics, vol. 12, no. 6, pp. 620–631, 1998
1998
-
[43]
A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows,
S. V. Patankar and D. B. Spalding, “A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows,” inNumerical prediction of flow, heat transfer, turbulence and combustion, pp. 54–73, Elsevier, 1983
1983
-
[44]
The numerical computation of turbulent flows,
B. E. Launder and D. B. Spalding, “The numerical computation of turbulent flows,” inNumerical prediction of flow, heat transfer, turbulence and combustion, pp. 96–116, Elsevier, 1983
1983
-
[45]
On the use of the k–εmodel in commercial cfd soft- ware to model the neutral atmospheric boundary layer,
D. Hargreaves and N. G. Wright, “On the use of the k–εmodel in commercial cfd soft- ware to model the neutral atmospheric boundary layer,”Journal of wind engineering and industrial aerodynamics, vol. 95, no. 5, pp. 355–369, 2007
2007
-
[46]
A comprehensive modelling approach for the neutral atmospheric boundary layer: consistent inflow conditions, wall function and turbulence model,
A. Parente, C. Gorlé, J. van Beeck, and C. Benocci, “A comprehensive modelling approach for the neutral atmospheric boundary layer: consistent inflow conditions, wall function and turbulence model,”Boundary-layer meteorology, vol. 140, no. 3, pp. 411–428, 2011
2011
-
[47]
EMEP/EEA air pollutant emission inventory guidebook 2023,
European Environment Agency, “EMEP/EEA air pollutant emission inventory guidebook 2023,” tech. rep., European Environment Agency, Copenhagen, Denmark, 2023. 35
2023
-
[48]
Scipy 1.0: fundamental al- gorithms for scientific computing in python,
P. Virtanen, R. Gommers, T. E. Oliphant, M. Haberland, T. Reddy, D. Cournapeau, E. Burovski, P. Peterson, W. Weckesser, J. Bright,et al., “Scipy 1.0: fundamental al- gorithms for scientific computing in python,”Nature methods, vol. 17, no. 3, pp. 261– 272, 2020
2020
-
[49]
Data-driven assessment of arch vortices in simplified urban flows,
Á. Martínez-Sánchez, E. Lazpita, A. Corrochano, S. Le Clainche, S. Hoyas, and R. Vinuesa, “Data-driven assessment of arch vortices in simplified urban flows,”In- ternational Journal of Heat and Fluid Flow, vol. 100, p. 109101, 2023. 36 Figure A1: From left to right: low-cost input, lcHOSVD reconstruction, and lcSVD reconstruction for the physical flow var...
2023
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.