pith. machine review for the scientific record. sign in

arxiv: 2605.10562 · v1 · submitted 2026-05-11 · 🧮 math.NA · cs.NA

Recognition: no theorem link

Data-driven moving-window Bayesian inference for transient CO2-temperature network models of buildings

Clemens V. Verhoosel, Idoia Cortes Garcia, Stein K.F. Stoter, Zhijian Wang

Pith reviewed 2026-05-12 05:04 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Bayesian inferenceCO2-temperature modelRC networkbuilding thermal dynamicsmoving-window estimationsparse sensingairflow modelingoccupancy estimation
0
0 comments X

The pith

Coupling CO2-informed airflow with thermal dynamics in a network model allows robust temperature prediction in buildings using sparse sensor data.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a network model that connects multi-zone CO2 transport and thermal dynamics through shared latent variables for airflow and occupancy. Thermal parts use an RC network with added convective terms from airflow, while CO2 uses convective transport equations. A moving-window Bayesian procedure jointly infers the parameters, trajectories, and noise levels from limited measurements, providing forecasts with uncertainty bands. This setup is tested on synthetic data and a physical experiment, showing accurate reconstruction and quick adaptation after regime changes. The result supports practical building monitoring and energy assessment when full sensing is unavailable.

Core claim

The authors propose that linking CO2 and temperature dynamics via shared airflow and occupancy in a network model, calibrated by moving-window Bayesian inference, yields accurate posterior reconstructions and low-error forecasts while providing diagnostics for model mismatch at abrupt changes.

What carries the argument

The CO2-temperature network model with shared latent airflow and occupancy, calibrated via moving-window Bayesian inference that estimates thermal parameters, trajectories, and sensor noise.

Load-bearing premise

The RC network with airflow convective exchange and CO2 transport, linked by shared latent drivers, sufficiently captures the building's dynamics, and the moving-window inference tracks changes without major bias or delay.

What would settle it

Observing that forecast errors stay high or uncertainty bands do not adapt after a known regime shift in the physical experiment would indicate the claim fails.

Figures

Figures reproduced from arXiv: 2605.10562 by Clemens V. Verhoosel, Idoia Cortes Garcia, Stein K.F. Stoter, Zhijian Wang.

Figure 1
Figure 1. Figure 1: Schematic of the room layout. To formalize the model, we define the sets of nodal states as ϕ = {ϕi} and T = {Ti}. Denoting the “unknown” (to be inferred) parameters that determine the CO2 and temper￾atuer models as θ ϕ and θ T , and denoting the “known” model parameters as κ ϕ and κ T , the forward models are expressed as the residuals r ϕ [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Schematic of local-to-global correspondence in the network. Here, ei denotes the ith oriented edge of the graph, where qi flows in the specified direction. The CO2 evolution is modeled as purely convective, and hence the network is described by an oriented graph whose nodes represent zones that store CO2 mass and whose edges represent inter-zonal convective exchange. For each (globally counted) node ni , t… view at source ↗
Figure 3
Figure 3. Figure 3: , we use data D = Dϕ ∪ DT contained in the current inference window to infer the parameter vector θ = θ ϕ ∪ θ T of the coupled CO2-temperature model. The inference step produces not only a single parameter estimate, but a range of plausible parameter values consistent with the noisy measurements. We use the resulting parameter samples to predict how CO2 and temperature behave beyond the inference window, a… view at source ↗
Figure 4
Figure 4. Figure 4: Schematic of Metropolis–Hastings sampling in one dimension. [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Illustration of the benchmark problem room layout and network structure. 4.1.2. Parameter configuration and data generation Our benchmark problem considers two sets of parameters, parameters to be inferred (θ = θ ϕ ∪ θ T ), and constants or known parameters (κ = κ ϕ ∪ κ T ). The parameters to be inferred and the constants are listed as: θ = {nppl, q, R, C}, (12a) κ = {BCϕ , BCT , {Vi}, qexh, ϕexh, Qppl, cp… view at source ↗
Figure 6
Figure 6. Figure 6: Synthetic ground-truth benchmark results for rooms A, E, and F, with and without noise. 4.2.2. Moving-window prediction and parameter tracking [PITH_FULL_IMAGE:figures/full_fig_p018_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Posterior predictive CO2 trajectories in rooms A, E and F for the synthetic-data experiment. 19 [PITH_FULL_IMAGE:figures/full_fig_p019_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Posterior predictive temperature trajectories in rooms A, E and F for the synthetic-data exper￾iment. The corresponding evolution of the inferred parameters of the coupled CO2 and tem￾perature model are shown in [PITH_FULL_IMAGE:figures/full_fig_p020_8.png] view at source ↗
Figure 8
Figure 8. Figure 8: Tracking of model parameters under moving-window Bayesian inference, showing the posterior mean (solid lines) and 95% credible intervals (shaded bands). variability, which is expected because multiple resistance configurations can yield similar temperature trajectories over a finite window. Overall, it indicates that the moving-window Bayesian inference tracks RC parameters in a temporally coherent and phy… view at source ↗
Figure 9
Figure 9. Figure 9: Inferred CO2 noise levels in rooms A, E and F. Solid lines denote posterior means and shaded bands indicate 95% credible intervals as a function of the window start time. 4.2.4. Change window size and noise level To evaluate the sensitivity of the Bayesian inference to the moving-window size and measurement noise, we vary the inference window size from 10 to 60 data points under noise level of synthetic CO… view at source ↗
Figure 10
Figure 10. Figure 10: nRMSE of prediction and ground-truth data, evaluated under different inference window size 10–60 and paired measurement-noise settings (σCO2 , σT ). (a) shows CO2 and (b) shows temperature. 5. Experimental validation To assess the predictive performance of our modeling framework on real measurement data that correspond to time-varying operating conditions, we now implement the moving￾window Bayesian infer… view at source ↗
Figure 11
Figure 11. Figure 11: Experiment setup Atmospheric reference sensors are installed outside the enclosure to record background CO2 and temperature for boundary conditions. To simulate human-body CO2 and heat production in room F, a premixed CO2 gas cylinder and a heating module were used. The premixed gas had a nominal CO2 concentra￾tion of 5000 ppm, as supplied by Nippon Gases [31]. The dosing flow rate is regulated to 25 mL/m… view at source ↗
Figure 12
Figure 12. Figure 12: Posterior predictive CO2 concentration (left column) and temperature (right column) trajec￾tories in rooms A, D, E, and F, for the scaled physical experiment. 29 [PITH_FULL_IMAGE:figures/full_fig_p029_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Posterior mean (solid lines) and 95% credible intervals (shaded bands) of the inferred occupan￾cies and fa¸cade airflow, for the scaled physical experiment. 5.2.3. Prediction evaluation To assess how well each posterior generalizes beyond the inference interval, we compute normalized root-mean-square error (nRMSE) between the model prediction and experimen￾tal measurements over a horizon of 80 samples aft… view at source ↗
Figure 14
Figure 14. Figure 14: Mean nRMSE over an 80 samples post-window prediction horizon, averaged over 8 rooms, for CO2 and temperature. room-wise inspection (not shown in the figure) indicates that the largest forecast errors con￾centrate in the directly actuated zone (room F), which is expected because it experiences the strongest source terms and the sharpest regime change. The remaining zones show smaller errors driven primaril… view at source ↗
read the original abstract

In this work, we proposes a CO2-temperature network model that links multi-zone mass transport and thermal dynamics through shared latent drivers, airflow and occupancy. The thermal component is formulated as a resistance-capacitance (RC) network augmented with airflow-driven convective exchange, while the CO2 component is governed by inter-zonal convective transport. To calibrate the model and track time-varying operating conditions based on sparse sensing, we introduce a moving-window Bayesian inference procedure that jointly estimates thermal parameters, airflow and occupancy trajectories. The estimation also provides room-specific sensor noise levels, yielding posterior predictive forecasts with credible intervals. The framework is assessed using a controlled synthetic benchmark, and a scaled physical validation experiment using CO2 and temperature sensing. In both settings, the posterior accurately reconstructs trajectories within windows and delivers low forecast errors. When inference windows overlap abrupt regime transitions, the widened uncertainty bands and increased inferred noise levels provide an interpretable diagnostic of model-data mismatch, followed by rapid recovery once the new regime is observed. Overall, coupling CO2-informed airflow with thermal dynamics yields a robust approach for conductive and advective temperature prediction, supporting practical monitoring and energy-performance assessment under limited sensing.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper proposes a CO2-temperature network model for multi-zone buildings that couples thermal RC dynamics (augmented with airflow-driven convective exchange) to inter-zonal CO2 convective transport through shared latent drivers (airflow and occupancy). A moving-window Bayesian inference procedure jointly estimates thermal parameters, airflow trajectories, occupancy, and sensor noise levels from sparse measurements, yielding posterior predictive forecasts with credible intervals. Validation on a controlled synthetic benchmark and a scaled physical experiment reports accurate trajectory reconstruction, low forecast errors, and interpretable widening of uncertainty bands at abrupt regime transitions followed by rapid recovery.

Significance. If the moving-window procedure can be shown to attribute predictive performance specifically to the CO2-informed airflow coupling rather than to local parameter refitting, the approach would provide a practical, data-driven framework for transient building monitoring and energy assessment under limited sensing. The explicit handling of regime transitions via uncertainty diagnostics is a methodological strength that could support real-world deployment.

major comments (2)
  1. [Abstract; moving-window Bayesian inference procedure] Abstract and the description of the moving-window Bayesian inference procedure: the claim that coupling CO2-informed airflow with thermal dynamics 'yields a robust approach' for prediction depends on the joint estimation correctly resolving convective exchange and transport via the shared latent drivers. However, the procedure jointly estimates thermal parameters alongside airflow and occupancy trajectories within each window. If the RC parameters are permitted to vary (or be re-estimated) locally rather than held globally fixed, discrepancies in advective terms can be absorbed into parameter adjustments instead of being enforced by the CO2-temperature linkage. This directly affects the central claim of robustness at regime transitions, where only widened uncertainty and 'rapid recovery' are reported without quantifying whether recovery is due to correct prediction or refitting.
  2. [Validation on synthetic benchmark and physical experiment] Validation sections (synthetic benchmark and physical experiment): the reported low forecast errors and accurate reconstruction within windows do not include an ablation or sensitivity test that fixes thermal parameters globally across windows while allowing only airflow/occupancy to vary. Without such a control, it remains unclear whether the CO2 coupling is load-bearing for the performance or whether the joint estimation masks model mismatch.
minor comments (2)
  1. [Abstract] Abstract: grammatical error in 'we proposes' should be corrected to 'we propose'.
  2. [Abstract] The abstract mentions 'room-specific sensor noise levels' but does not clarify how these are parameterized or whether they are estimated per window or globally; this notation should be made explicit for reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

Thank you for the constructive feedback. We address the two major comments point-by-point below and commit to revisions that include the suggested ablation studies to better isolate the contribution of the CO2 coupling.

read point-by-point responses
  1. Referee: Abstract and the description of the moving-window Bayesian inference procedure: the claim that coupling CO2-informed airflow with thermal dynamics 'yields a robust approach' for prediction depends on the joint estimation correctly resolving convective exchange and transport via the shared latent drivers. However, the procedure jointly estimates thermal parameters alongside airflow and occupancy trajectories within each window. If the RC parameters are permitted to vary (or be re-estimated) locally rather than held globally fixed, discrepancies in advective terms can be absorbed into parameter adjustments instead of being enforced by the CO2-temperature linkage. This directly affects the central claim of robustness at regime transitions, where only widened uncertainty and 'rapid recovery' are reported without quantifying whether recovery is due to correct prediction or refitting.

    Authors: We appreciate this insightful comment highlighting a potential ambiguity in attributing performance to the model coupling versus local parameter estimation. The moving-window Bayesian inference is designed to adapt to changing conditions, but we recognize that without additional controls, the role of the CO2-temperature linkage in enforcing consistency may not be fully isolated. In the revised manuscript, we will add an explicit discussion of this aspect and perform a sensitivity analysis on the synthetic data to quantify the impact. This will support our claim of robustness by showing that the coupling aids in correct prediction rather than solely relying on refitting. revision: yes

  2. Referee: Validation sections (synthetic benchmark and physical experiment): the reported low forecast errors and accurate reconstruction within windows do not include an ablation or sensitivity test that fixes thermal parameters globally across windows while allowing only airflow/occupancy to vary. Without such a control, it remains unclear whether the CO2 coupling is load-bearing for the performance or whether the joint estimation masks model mismatch.

    Authors: We agree that the absence of such an ablation test leaves open the possibility that joint estimation is compensating for model discrepancies. To address this, we will implement and report an ablation study in the revised paper. Specifically, for the synthetic benchmark, we will re-run the inference with thermal RC parameters held fixed at their globally estimated values (from the full dataset or initial window) and compare the resulting forecast errors and reconstruction accuracy against the original results where parameters are allowed to vary per window. A similar analysis will be attempted for the physical experiment to the extent possible with available data and computation. revision: yes

Circularity Check

0 steps flagged

No significant circularity; data-driven inference is self-contained

full rationale

The paper describes a moving-window Bayesian procedure that jointly estimates thermal RC parameters, airflow, and occupancy trajectories directly from sparse CO2 and temperature measurements. Validation occurs on independent synthetic benchmarks and physical experiments, with performance assessed via posterior predictive forecasts and credible intervals. No derivation step reduces a claimed prediction to a fitted input by construction, no self-citation is invoked as a uniqueness theorem or load-bearing premise, and the model equations (RC network augmented with convective terms) are calibrated rather than tautologically defined from the outputs they produce. The approach therefore remains externally falsifiable against held-out data segments and does not exhibit any of the enumerated circularity patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 3 axioms · 1 invented entities

The central claim rests on standard building physics modeling assumptions and Bayesian statistical procedures; free parameters are the quantities estimated from data rather than ad hoc choices.

axioms (3)
  • domain assumption Building thermal dynamics can be represented as an RC network augmented with airflow-driven convective exchange.
    Core of the thermal component formulation in the abstract.
  • domain assumption CO2 dynamics are governed by inter-zonal convective transport linked to the same airflow.
    Core of the CO2 component and coupling in the abstract.
  • domain assumption Moving-window Bayesian inference can jointly estimate time-varying parameters, trajectories, and noise from sparse sensing.
    Core of the calibration procedure described.
invented entities (1)
  • shared latent drivers (airflow and occupancy) no independent evidence
    purpose: Link multi-zone mass transport and thermal dynamics
    Introduced as inferred variables connecting the CO2 and temperature components.

pith-pipeline@v0.9.0 · 5521 in / 1750 out tokens · 56481 ms · 2026-05-12T05:04:40.603231+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages

  1. [1]

    Energy Performance of Buildings Directive,

    European Commission, “Energy Performance of Buildings Directive,” 2024. Accessed: 2026-03-22

  2. [2]

    The state of district heating and cooling in europe: A literature-based assessment,

    V. Mun´ can, I. Mujan, D. Macura, and A. S. Andelkovi´ c, “The state of district heating and cooling in europe: A literature-based assessment,”Energy, vol. 304, p. 132191, 2024

  3. [3]

    Performance optimization of an office building with a dynamic fa¸ cade system,

    Z. Wang, M. Hajdukiewicz, H. Zheng, P.-J. Hoes, and M. G. L. C. Loomans, “Performance optimization of an office building with a dynamic fa¸ cade system,”International Journal of Ventilation, pp. 1–13, 2025

  4. [4]

    Occupancy-based hvac control systems in buildings: A state-of-the-art review,

    M. Esrafilian-Najafabadi and F. Haghighat, “Occupancy-based hvac control systems in buildings: A state-of-the-art review,”Building and Environment, vol. 197, p. 107810, 2021. 33

  5. [5]

    A workflow for evaluating occupant- centric controls using building simulation,

    B. W. Hobson, B. Huchuk, H. B. Gunay, and W. O’Brien, “A workflow for evaluating occupant- centric controls using building simulation,”Journal of Building Performance Simulation, vol. 14, no. 6, pp. 645–661, 2021

  6. [6]

    The information gap in occupant-centric building operations: Lessons learned from interviews with building operators in germany,

    J. Hahn, J.-K. Kazak, and A. Wagner, “The information gap in occupant-centric building operations: Lessons learned from interviews with building operators in germany,”Frontiers in Built Environment, vol. 8, p. 838859, 2022

  7. [7]

    Grey-box modeling and application for building energy simulations—a critical review,

    L. Zhang, Y. Li, Z. O’Neill, J. Chen, P. Im, and J. DeGraw, “Grey-box modeling and application for building energy simulations—a critical review,”Renewable and Sustainable Energy Reviews, vol. 146, p. 111174, 2021

  8. [8]

    Review on building energy model calibration by bayesian infer- ence,

    D. Hou, I. G. Hassan, and L. Wang, “Review on building energy model calibration by bayesian infer- ence,”Renewable and Sustainable Energy Reviews, vol. 143, p. 110930, 2021

  9. [9]

    Guidelines for the bayesian calibration of building energy models,

    A. Chong and K. Menberg, “Guidelines for the bayesian calibration of building energy models,”Energy and Buildings, vol. 174, pp. 527–547, 2018

  10. [10]

    Parameter identification approach to represent building thermal dynamics reducing tuning time of control system gains: A case study in a tropical climate,

    A. K. Rivera, J. S´ anchez, and M. Chen Austin, “Parameter identification approach to represent building thermal dynamics reducing tuning time of control system gains: A case study in a tropical climate,” Frontiers in Built Environment, vol. 8, p. 949426, 2022

  11. [11]

    Development of physically coherent grey-box models for residential buildings using a simplified adjustment method,

    E. M. Lecocq, J. Pascual, and J. Salom, “Development of physically coherent grey-box models for residential buildings using a simplified adjustment method,”Energy and Buildings, 2025. In press

  12. [12]

    Coupling the multizone airflow and contaminant transport software contam with energyplus using co-simulation,

    W. S. Dols, B. J. Polidoro, and S. J. Emmerich, “Coupling the multizone airflow and contaminant transport software contam with energyplus using co-simulation,”Building Simulation, vol. 9, no. 4, pp. 469–479, 2016

  13. [13]

    Carbon dioxide-based occupancy estimation using stochastic differential equations,

    S. Wolf, D. Cal` ı, J. Krogstie, and H. Madsen, “Carbon dioxide-based occupancy estimation using stochastic differential equations,”Applied Energy, vol. 236, pp. 32–41, 2019

  14. [14]

    Co 2-based grey- box model to estimate airflow rate and room occupancy,

    S. Wolf, M. Justo Alonso, D. Cal` ı, J. Krogstie, H. M. Mathisen, and H. Madsen, “Co 2-based grey- box model to estimate airflow rate and room occupancy,” inE3S Web of Conferences (CLIMA 2019), vol. 111, p. 04036, 2019

  15. [15]

    Simulation-based occupancy estimation in office buildings using co2 sensors,

    F. Jorissen, W. Boydens, and L. Helsen, “Simulation-based occupancy estimation in office buildings using co2 sensors,” inProceedings of Building Simulation 2017: 15th Conference of IBPSA, pp. 1073– 1082, 2017

  16. [16]

    Low-cost data-driven estimation of indoor occupancy based on carbon dioxide (co 2) concentration: A multi-scenario case study,

    X. Liang, J. Shim, O. Anderton, and D. Song, “Low-cost data-driven estimation of indoor occupancy based on carbon dioxide (co 2) concentration: A multi-scenario case study,”Journal of Building Engi- neering, vol. 96, p. 108180, 2024

  17. [17]

    A framework for the continuous calibration of building energy models with uncertainty,

    A. Chong and S. Chao, “A framework for the continuous calibration of building energy models with uncertainty,” inProceedings of Building Simulation 2019 (IBPSA), (Rome, Italy), 2019

  18. [18]

    Continuous-time bayesian calibration of energy models using bim and energy data,

    A. Chong, W. Xu, S. Chao, and N.-T. Ngo, “Continuous-time bayesian calibration of energy models using bim and energy data,”Energy and Buildings, vol. 194, pp. 177–190, 2019

  19. [19]

    Bayesian filtering for building occupancy estimation from carbon dioxide concentration,

    C. Jiang, Z. Chen, R. Su, M. K. Masood, and Y. C. Soh, “Bayesian filtering for building occupancy estimation from carbon dioxide concentration,”Energy and Buildings, vol. 206, p. 109566, 2020

  20. [20]

    Inference of building occupancy signals using moving horizon estimation and fourier regularization,

    V. M. Zavala, “Inference of building occupancy signals using moving horizon estimation and fourier regularization,” Tech. Rep. MCS/ANL-P4031-0213, Argonne National Laboratory, Argonne, IL, USA, 2013

  21. [21]

    SPICE2: a computer program to simulate semiconductor circuits,

    L. W. Nagel, “SPICE2: a computer program to simulate semiconductor circuits,” tech. rep., University of Berkeley, 1975

  22. [22]

    Development of moving window state and parameter estimators under maximum likelihood and bayesian frameworks,

    J. Valluru, R. Kumar,et al., “Development of moving window state and parameter estimators under maximum likelihood and bayesian frameworks,”Journal of Process Control, 2017

  23. [23]

    Gelman, J

    A. Gelman, J. B. Carlin, H. S. Stern, D. B. Dunson, A. Vehtari, and D. B. Rubin,Bayesian Data Analysis. Boca Raton, FL, USA: Chapman & Hall/CRC Texts in Statistical Science, 3 ed., 2013

  24. [25]

    Weak convergence and optimal scaling of random walk metropolis algorithms,

    G. O. Roberts, A. Gelman, and W. R. Gilks, “Weak convergence and optimal scaling of random walk metropolis algorithms,”The Annals of Applied Probability, vol. 7, no. 1, pp. 110–120, 1997. 34

  25. [26]

    A tutorial on adaptive mcmc,

    C. Andrieu and J. Thoms, “A tutorial on adaptive mcmc,”Statistics and Computing, vol. 18, no. 4, pp. 343–373, 2008

  26. [27]

    Robust adaptive metropolis algorithm with coerced acceptance rate,

    M. Vihola, “Robust adaptive metropolis algorithm with coerced acceptance rate,”Statistics and Com- puting, vol. 22, no. 5, pp. 997–1008, 2012

  27. [28]

    Improvement to en iso 52016-1:2017 hourly heat transfer through opaque envelope elements: A comparison with detailed simulation,

    L. Mazzarella, R. Scoccia, P. Colombo, and M. Motta, “Improvement to en iso 52016-1:2017 hourly heat transfer through opaque envelope elements: A comparison with detailed simulation,”Energy and Buildings, vol. 214, p. 109882, 2020

  28. [29]

    Vaisala carbocap carbon dioxide probe gmp252,

    Vaisala Oyj, “Vaisala carbocap carbon dioxide probe gmp252,” 2024

  29. [30]

    Vaisala humicap humidity and temperature probe hmp1,

    Vaisala Oyj, “Vaisala humicap humidity and temperature probe hmp1,” 2024

  30. [31]

    Certified gas standard for laboratories,

    Nippon Gases, “Certified gas standard for laboratories,” 2026

  31. [32]

    Air sampling pump gsp-300ft-2,

    GASTEC Corporation, “Air sampling pump gsp-300ft-2,” 2012

  32. [33]

    Bayesian Online Changepoint Detection

    R. P. Adams and D. J. C. MacKay, “Bayesian online changepoint detection,”arXiv preprint arXiv:0710.3742, 2007. 35