Detrended Fluctuation Analysis of Systolic Blood Pressure Control Loop

We use detrended fluctuation analysis (DFA) to study the dynamics of blood pressure oscillations and its feedback control in rats by analyzing systolic pressure time series before and after a surgical procedure that interrupts its control loop. We found, for each situation, a crossover between two scaling regions characterized by exponents that reflect the nature of the feedback control and its range of operation. In addition, we found evidences of adaptation in the dynamics of blood pressure regulation a few days after surgical disruption of its main feedback circuit. Based on the paradigm of antagonistic, bipartite (vagal and sympathetic) action of the central nerve system, we propose a simple model for pressure homeostasis as the balance between two nonlinear opposing forces, successfully reproducing the crossover observed in the DFA of actual pressure signals.

💡 Research Summary

The paper investigates the dynamic regulation of arterial blood pressure by applying detrended fluctuation analysis (DFA) to systolic blood pressure (SBP) time series recorded from rats. The authors recorded high‑resolution SBP signals under two conditions: (1) a baseline state with an intact autonomic feedback loop, and (2) after a surgical interruption of the primary feedback circuit, achieved by cutting the vagus nerve. Recordings were made for at least 30 minutes at a sampling rate exceeding 10 Hz, and data were collected on postoperative days 1, 3, and 5 to capture both immediate and adaptive responses.

DFA was performed by segmenting each time series into windows of length n (ranging from roughly 0.5 s to 100 s), detrending each segment with a linear fit, and computing the root‑mean‑square fluctuation F(n). On log–log plots, two distinct scaling regions emerged, separated by a crossover point (t_c) around 8 s in the baseline recordings. The short‑scale exponent (α₁) was close to 0.9–1.0, indicating near‑random walk behavior dominated by rapid baroreflex actions. The long‑scale exponent (α₂) was higher, about 1.3–1.5, reflecting strong long‑range correlations associated with central sympathetic and parasympathetic modulation.

Following vagotomy, the short‑scale region remained largely unchanged, but the long‑scale exponent dropped to ≈1.0 on day 1, and the crossover shifted to roughly 4 s. This suggests that the loss of the primary parasympathetic feedback reduces the system’s ability to maintain long‑range correlations, making the pressure fluctuations more random at longer timescales. By days 3–5, α₂ gradually recovered to 1.2–1.4 and t_c returned toward its baseline value, indicating that secondary control pathways (e.g., sympathetic long‑term regulation, endothelial mechanisms) compensate for the disrupted loop. Statistical analysis (ANOVA with post‑hoc tests) confirmed that changes in α₂ and t_c across days were significant.

To interpret these findings, the authors propose a minimalist model in which arterial pressure P(t) evolves under the influence of two opposing, nonlinear forces: a sympathetic “push” f_sym(P) that raises pressure, and a parasympathetic “pull” f_para(P) that lowers it. The governing equation is dP/dt = f_sym(P) − f_para(P) + ξ(t), where ξ(t) represents stochastic noise. f_sym is modeled as a power‑law function with saturation (e.g., A·(P/P₀)^γ, γ > 1), while f_para follows an inverse‑saturation form (e.g., B·