skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Reconstructing cardiac electrical excitations from optical mapping recordings
The reconstruction of electrical excitation patterns through the unobserved depth of the tissue is essential to realizing the potential of computational models in cardiac medicine. We have utilized experimental optical-mapping recordings of cardiac electrical excitation on the epicardial and endocardial surfaces of a canine ventricle as observations directing a local ensemble transform Kalman filter data assimilation scheme. We demonstrate that the inclusion of explicit information about the stimulation protocol can marginally improve the confidence of the ensemble reconstruction and the reliability of the assimilation over time. Likewise, we consider the efficacy of stochastic modeling additions to the assimilation scheme in the context of experimentally derived observation sets. Approximation error is addressed at both the observation and modeling stages through the uncertainty of observations and the specification of the model used in the assimilation ensemble. We find that perturbative modifications to the observations have marginal to deleterious effects on the accuracy and robustness of the state reconstruction. Furthermore, we find that incorporating additional information from the observations into the model itself (in the case of stimulus and stochastic currents) has a marginal improvement on the reconstruction accuracy over a fully autonomous model, while complicating the model itself and thus introducing potential for new types of model errors. That the inclusion of explicit modeling information has negligible to negative effects on the reconstruction implies the need for new avenues for optimization of data assimilation schemes applied to cardiac electrical excitation.  more » « less
Award ID(s):
1806833
PAR ID:
10588630
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
American Institute of Physics
Date Published:
Journal Name:
Chaos: An Interdisciplinary Journal of Nonlinear Science
Volume:
33
Issue:
9
ISSN:
1054-1500
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; (, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences)
    null (Ed.)
    Modelling of cardiac electrical behaviour has led to important mechanistic insights, but important challenges, including uncertainty in model formulations and parameter values, make it difficult to obtain quantitatively accurate results. An alternative approach is combining models with observations from experiments to produce a data-informed reconstruction of system states over time. Here, we extend our earlier data-assimilation studies using an ensemble Kalman filter to reconstruct a three-dimensional time series of states with complex spatio-temporal dynamics using only surface observations of voltage. We consider the effects of several algorithmic and model parameters on the accuracy of reconstructions of known scroll-wave truth states using synthetic observations. In particular, we study the algorithm’s sensitivity to parameters governing different parts of the process and its robustness to several model-error conditions. We find that the algorithm can achieve an acceptable level of error in many cases, with the weakest performance occurring for model-error cases and more extreme parameter regimes with more complex dynamics. Analysis of the poorest-performing cases indicates an initial decrease in error followed by an increase when the ensemble spread is reduced. Our results suggest avenues for further improvement through increasing ensemble spread by incorporating additive inflation or using a parameter or multi-model ensemble. This article is part of the theme issue ‘Uncertainty quantification in cardiac and cardiovascular modelling and simulation’. 
    more » « less
  2. ; ; (, Monthly Weather Review)
    Abstract For data assimilation to provide faithful state estimates for dynamical models, specifications of observation uncertainty need to be as accurate as possible. Innovation-based methods based on Desroziers diagnostics, are commonly used to estimate observation uncertainty, but such methods can depend greatly on the prescribed background uncertainty. For ensemble data assimilation, this uncertainty comes from statistics calculated from ensemble forecasts, which require inflation and localization to address under sampling. In this work, we use an ensemble Kalman filter (EnKF) with a low-dimensional Lorenz model to investigate the interplay between the Desroziers method and inflation. Two inflation techniques are used for this purpose: 1) a rigorously tuned fixed multiplicative scheme and 2) an adaptive state-space scheme. We document how inaccuracies in observation uncertainty affect errors in EnKF posteriors and study the combined impacts of misspecified initial observation uncertainty, sampling error, and model error on Desroziers estimates. We find that whether observation uncertainty is over- or underestimated greatly affects the stability of data assimilation and the accuracy of Desroziers estimates and that preference should be given to initial overestimates. Inline estimates of Desroziers tend to remove the dependence between ensemble spread–skill and the initially prescribed observation error. In addition, we find that the inclusion of model error introduces spurious correlations in observation uncertainty estimates. Further, we note that the adaptive inflation scheme is less robust than fixed inflation at mitigating multiple sources of error. Last, sampling error strongly exacerbates existing sources of error and greatly degrades EnKF estimates, which translates into biased Desroziers estimates of observation error covariance. Significance StatementTo generate accurate predictions of various components of the Earth system, numerical models require an accurate specification of state variables at our current time. This step adopts a probabilistic consideration of our current state estimate versus information provided from environmental measurements of the true state. Various strategies exist for estimating uncertainty in observations within this framework, but are sensitive to a host of assumptions, which are investigated in this study. 
    more » « less
  3. ; ; ; ; (, La Matematica)
    Using a simplified two-stage ice sheet model, we explore the potential of statistical data assimilation methods to improve predictions of glacier melt, which has significant implications for reducing uncertainty in projections of sea level rise. Through twin experiments utilizing artificial data, we find that the ensemble Kalman filter improves simulations of glacier evolution initialized with incorrect initial conditions and parameters, providing us with better predictions of future glacier melt. We explore the number of observations necessary to produce an accurate model run. We also explore optimal observation assimilation schemes, and determine that deviations from the true glacier response that stem from having few data points in the pre-satellite era can be corrected with modern observation data. Our results show that statistical data assimilation methods have great potential to improve complex glacier models using real-world observations. 
    more » « less
  4. ; (, Remote Sensing)
    Ensemble Kalman filters are an efficient class of algorithms for large-scale ensemble data assimilation, but their performance is limited by their underlying Gaussian approximation. A two-step framework for ensemble data assimilation allows this approximation to be relaxed: The first step updates the ensemble in observation space, while the second step regresses the observation state update back to the state variables. This paper develops a new quantile-conserving ensemble filter based on kernel-density estimation and quadrature for the scalar first step of the two-step framework. It is shown to perform well in idealized non-Gaussian problems, as well as in an idealized model of assimilating observations of sea-ice concentration. 
    more » « less
  5. ; ; ; ; ; ; ; ; ; (, Science Advances)
    Multiparametric investigation of cardiac physiology is crucial for the diagnosis and therapy of heart disease. However, no method exists to simultaneously map multiple parameters that govern cardiac (patho)physiology from beating hearts in vivo. Here, we present a cardiac sensing platform that addresses this challenge, functioning with a wireless interface. Advanced fabrication and assembling strategies enable the heterogeneous integration of transparent microelectrodes, light-emitting diodes, photodiodes, and optical filters into a multilayer array structure on soft substrates. The microelectrodes exhibit superior electrochemical performance for measuring electrical potentials and excellent transparency for co-localized fluorescence measurement. The device shows excellent biocompatibility and records the fluorescence of calcium reporter with performance comparable to imaging cameras. Multiparametric in vivo mapping of electrical excitation, calcium dynamics, and their combined effects on cardiac excitation-contraction coupling is demonstrated during normal rhythm, arrhythmia, and treatment. This technology offers potential widespread use in cardiac research to support scientific discoveries and advance clinical life-saving diagnostics and therapies. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email