The problem stems from a fundamental property of EEG. Electrodes do not measure absolute electrical potentials; they measure potential differences relative to a reference. Because no perfect, silent reference exists, researchers must choose one, and the choice shapes everything downstream. Previous work by the same group, based at the Universidad Autónoma del Estado de Morelos and the Universidad Nacional Autónoma de México, had already established that the median reference, in which each electrode’s signal is expressed relative to the median of all electrodes at every instant, produces functional brain networks that closely resemble those obtained from magnetoencephalography, a technique that is essentially reference-free. The median reference, it seemed, was the safest choice.

But EEG data also need to be filtered, typically into the classical frequency bands such as delta, theta, alpha, beta, and gamma, which are associated with different cognitive and physiological states. That raises a question that sounds almost trivial: should the signals be filtered first and then re-referenced to the median, a scheme the authors call FR, or re-referenced first and then filtered, a scheme they call RF? Intuition might favor RF, since re-referencing only once reduces the computational workload. The new study demonstrates that this shortcut comes at a steep and sometimes irreparable price.

The team, led by J. Daniel Arzate-Mena, Paola V. Olguín-Rodríguez, and Markus F. Müller, tested both schemes on three kinds of data: whole-night sleep EEG recordings from ten healthy young men, synthetic signals from controlled numerical models, and surrogate data that mimic the statistical properties of real EEG while containing no genuine correlations at all. The numerical models included Nf-Tori, which are linear superpositions of two thousand harmonic oscillations whose power spectra could be tuned to follow power laws of varying steepness, and IAAFT surrogates, which preserve the amplitude distribution and power spectrum of real recordings while destroying any true interdependence between channels. This combination allowed the researchers to know exactly what the correct answer should be, and to see precisely how each preprocessing scheme deviated from it.

The results were striking. When re-referencing was performed before filtering, the RF scheme, spurious positive cross-correlations appeared throughout the network, and they grew systematically stronger in the faster frequency bands. In the gamma band, during deep sleep, the average displacement of correlation values reached about 0.5, with individual estimates climbing as high as 0.75. The true spatial correlation structure was overshadowed by this flood of artificial correlations, simulating a kind of hypersynchronized brain dynamics that simply does not exist. In contrast, the FR scheme introduced only small, systematic negative shifts of roughly minus 0.025, which were not statistically significant.

Most alarming was the finding that the RF-induced artifacts are irreversible. The researchers applied a sophisticated statistical correction known as the Significant Average Correlation matrix, or SAC-matrix, which compares each pair of channels against surrogate data at a Bonferroni-corrected significance level to separate genuine correlations from random ones. For the FR scheme, this correction worked beautifully, essentially eliminating the tiny residual distortions. But for the RF scheme, even the SAC-matrix could not repair the damage. The artificial correlations survived the correction, meaning they become indistinguishable from genuine brain connectivity. Once the redundant information has been embedded into every channel by the flawed pipeline, no post-hoc statistical trick can tease it apart.

The mechanism behind the effect lies in the peculiar behavior of the median reference itself. At any given moment, the median of an odd number of channels equals the potential of one specific electrode, which is then subtracted from all the others. The researchers measured how long each electrode serves as the reference before another takes over. For steep power spectra, they found that high-amplitude, slow-frequency components can dominate the median signal for hundreds of consecutive sampling points. During these extended intervals, the same slow waveform is effectively stamped into every channel, embedding redundant information across the entire recording. When the data are subsequently filtered into fast frequency bands, the mismatch between the slow time scale of the reference signal and the fast time scale of the band creates the spurious correlations, with the effect most pronounced in the beta and gamma ranges.

The severity of the problem is governed by the power spectral exponent beta, the slope of the signal’s power spectrum. When the team tuned their synthetic signals to follow power laws from white noise, with beta equal to zero, up to steep red noise with beta equal to three, the artificial correlations grew with every increase in steepness. This matters because beta values of two or greater are not exotic. They characterize slow-wave sleep, the post-ictal period following focal epileptic seizures, anesthesia, neonatal recordings, and patients in vegetative states or coma. In the sleep data, the spectral slope correlated strongly with the amount of spurious gamma-band correlation under the RF scheme, with a correlation coefficient of roughly 0.85, while the FR scheme showed no such dependence. Deep sleep stages, with slopes approaching minus three, were hit hardest.

The implications reach well beyond sleep research. Any study that constructs functional brain networks from multivariate EEG recordings, whether investigating epilepsy, disorders of consciousness, infant development, or insomnia, could be producing distorted connectivity maps if the preprocessing order is wrong. The authors also showed that the global average reference, an alternative proposed back in 1950, suffers from its own flaw: high-amplitude channels dominate the average and induce spurious coupling in low-amplitude channels, a pattern that even the SAC-matrix cannot fully correct. The median reference, by contrast, is far more robust to amplitude differences caused by muscle or eye artifacts.

The team’s recommendation is clear: filter first, then re-reference, and always apply surrogate-corrected cross-correlations when analyzing signals with steep power spectra. The FR scheme demands more computation, since each filtered signal must be re-referenced individually, but the resulting distortions are an order of magnitude smaller and can be almost completely eliminated. Although the study was conducted on EEG, the authors emphasize that the principle applies to any multivariate recording in which spatial re-referencing and band-pass filtering are applied sequentially. In an era when brain connectivity studies inform everything from consciousness research to clinical diagnosis, the message is a sobering reminder that in signal processing, as in life, the order of operations is not a technical detail. It can be the difference between discovering how the brain truly communicates and chasing ghosts that the pipeline itself created.

Subject of Research: The effect of preprocessing step order on spurious correlations in EEG functional network analysis

Article Title: Irreversible functional network deformation induced by preprocessing scheme order in multivariate EEG recordings

Article References: Arzate-Mena, J. D., Olguín-Rodríguez, P. V., Corsi-Cabrera, M., López-Alvarado, E. A., & Müller, M. F. (2026). Irreversible functional network deformation induced by preprocessing scheme order in multivariate EEG recordings. Frontiers in Neuroinformatics, 20, Article 1907920. https://doi.org/10.3389/fninf.2026.1907920

Image Credits: AI Generated

DOI: 10.3389/fninf.2026.1907920

Keywords: EEG, preprocessing, functional connectivity, median reference, spurious correlations, spectral slope, sleep, surrogate data, cross-correlation, neuroinformatics, band-pass filtering, brain networks

The research, published in Physiological Reports, focused on high-frequency heart rate variability, or HF-HRV, a measure derived from electrocardiogram recordings that reflects beat-to-beat fluctuations in heart rhythm within a specific frequency band. HF-HRV is closely associated with respiratory sinus arrhythmia, the natural speeding and slowing of the heart with breathing, which in turn is linked to parasympathetic, or vagal, nervous system activity. Because this branch of the autonomic nervous system is thought to support calmness, social engagement, and emotional regulation, HF-HRV has become a popular window onto the physiological states of caregivers and children during interaction.

Previous studies of physiological synchrony have typically averaged heart rate variability over fixed time windows and compared those averages between partners. Such approaches can detect gross similarity but may miss a crucial feature: the possibility that one person’s physiological rhythm precedes or follows the other’s by a meaningful interval. Eye contact, vocal exchanges, and physical touch between parents and children can change within seconds, and the physiological consequences of those moments may ripple forward or backward in time. The Japanese team, led by Aoi Sakurai, set out to build and test a method that explicitly accounts for these temporal shifts.

The study recruited 15 mother-infant dyads with healthy infants aged five to seven months, recruited between December 2022 and March 2023 through local childcare newspapers and social media in Sendai City. Each pair was scheduled for two recording sessions roughly two weeks apart. Mothers were asked to feed and change their infants beforehand to ensure calm, and then both wore small wireless ECG sensors, positioned at the upper sternum for mothers and near the clavicular region for infants, while spending five minutes in unrestricted free play on a yoga mat, talking, touching, and playing as they naturally would.

Turning raw infant ECG into usable data proved to be the study’s greatest technical hurdle. Of the 30 recorded sessions, 14 were excluded because of device failures, session interruptions, or electrode detachment and movement artifacts, leaving 16 analyzable sessions from 10 unique dyads. The researchers processed the surviving signals with a custom MATLAB pipeline that combined automated QRS detection with careful visual inspection and manual correction, flagging suspicious beats as missing values rather than interpolating them. Artifact rates in the retained data were remarkably low, at roughly 0.29 percent for mothers and 0.02 percent for infants, a testament to the labor-intensive cleaning process.

The core of the analysis was a lag-sensitive cross-correlation technique. The team computed HF-HRV power in 30-second windows advanced in half-second steps, using frequency bands of 0.15 to 0.40 hertz for mothers and 0.24 to 1.04 hertz for infants, reflecting developmental differences in breathing rate. For each maternal reference window, the infant’s signal was shifted across a search range spanning minus 30 to plus 30 seconds, and the correlation between the two traces was calculated at every offset. When a positive correlation coefficient of 0.8 or higher emerged as a local peak, that moment was marked as a synchronization point, classified as Lag-plus if the infant’s pattern aligned after a positive offset, Lag-minus if the mother’s pattern aligned after a negative offset, or Lag-zero for exact alignment.

The results were striking. Highly correlated intervals accounted for a mean of 74.5 percent of the total analyzed interaction time, a figure well above the variable or lower synchrony levels often reported in previous mother-infant studies. More intriguing still was the distribution of lags: within the highly correlated time, Lag-minus intervals made up about 54.4 percent and Lag-plus intervals about 44.6 percent, while exact zero-lag alignment accounted for a mere 1.1 percent. In other words, under this analytical definition, mother-infant heart rhythm synchrony was almost always offset in time, and the direction of the offset varied both across sessions and within them, with some dyads showing a bias toward one direction over limited periods.

Because high correlations between two slowly fluctuating biological signals can arise by chance from autocorrelation and shared measurement conditions, the team also compared their observations against a null model. For each session, they generated 100 surrogate datasets by randomly permuting the order of 30-second segments within the maternal and infant time series, destroying the original correspondence while preserving local structure. Observed synchrony exceeded the session-specific surrogate distribution in 13 of the 16 sessions, with a mean observed-minus-surrogate difference of about 0.205, and the positive difference held across all 10 dyads in a sensitivity analysis. This suggests that the alignment was more than coincidental temporal structure, though the authors are careful to frame the surrogate test as exploratory, since it relied on a single null model and a modest number of iterations.

The theoretical backdrop is Polyvagal Theory, which positions the vagal system as central to physiological co-regulation in close relationships. The findings are consistent with the idea that caregivers and infants mutually regulate one another’s autonomic states, and they echo earlier reports of heart rhythm synchrony between pregnant women and their fetuses, between marital partners, and even between dogs and their owners. Yet the authors caution against overinterpretation. The lag categories should not be read as evidence that one partner behaviorally leads or regulates the other, and the study measured HF-HRV without recording respiration, so the results reflect temporal physiological coordination rather than respiration-controlled estimates of vagal activity.

The limitations are candidly acknowledged. The sample was small and limited to mothers and five- to seven-month-old infants, with no fathers, no baseline or still-face conditions, and no non-partner dyad pairings to rule out shared environment or common movement. No synchronized behavioral coding was performed, so the physiological alignment cannot yet be tied to specific events like eye contact, vocalization, or touch. The predefined plus-or-minus 30-second lag range may itself have shaped the findings. Still, the study demonstrates the feasibility of a lag-sensitive approach that could, if replicated in larger and more diverse samples, eventually help researchers probe whether physiological synchrony differs in populations facing psychosocial risk, such as mothers with postpartum depression, offering an objective, minimally invasive complement to the behavioral observations that have long dominated attachment research.

Subject of Research: Lag-sensitive analysis of mother-infant heart rate variability synchrony during free interaction

Article Title: Lag‐sensitive cross‐correlation analysis of mother–infant HF‐HRV synchrony during free interaction: An exploratory feasibility study

Article References: Sakurai, A., Koide, K., Kasahara, Y., Nakamura, Y., & Yoshizawa, T. (2026). Lag‐sensitive cross‐correlation analysis of mother–infant HF ‐ HRV synchrony during free interaction: An exploratory feasibility study. Physiological Reports, 14(18), Article e71116. https://doi.org/10.14814/phy2.71116

Image Credits: AI Generated

DOI: 10.14814/phy2.71116

Keywords: heart rate variability, mother-infant synchrony, HF-HRV, cross-correlation, autonomic nervous system, respiratory sinus arrhythmia, polyvagal theory, attachment, physiological co-regulation, surrogate data, free play, infant development

The green line at 530.3 nanometers is the brightest emission line in the solar corona, produced by highly ionized iron atoms, Fe XIV, stripped of thirteen of their electrons by the extreme coronal temperatures of roughly two million kelvin. Unlike the photospheric surface we see with the naked eye, the corona traces magnetic structures that arch high above active regions, and its emission intensity varies with the strength and distribution of the underlying magnetic field. Because ground-based observatories, coordinated through networks that produced the green corona database and the coronal index of solar activity, have monitored this emission consistently since the mid-twentieth century, the green line provides an exceptionally long and homogeneous record of large-scale coronal behavior that modern space-based instruments simply cannot match in duration.

On the sunspot side, the study departs from tradition in a crucial way. Conventional analyses of the solar cycle often rely on the sunspot number or the total sunspot area, both of which collapse the entire two-dimensional distribution of spottedness into a single scalar per day. Takalo instead constructed a Gaussian representation of the sunspot activity field, built from the individual areas and latitudes of every recorded sunspot group. This activity-field model preserves where on the solar disk the magnetic flux emerges, capturing the characteristic equatorward migration of sunspot belts as each cycle progresses and the poleward drift of its remnants. The result is a mathematical description of the large-scale spatial organization of solar magnetic activity that is far more realistic than any single-number index.

The central discovery of the study is a remarkably stable, double-peaked correlation structure in latitude. When coronal green-line emission at each latitude band is correlated against the modeled sunspot activity field at the same latitude, the correlation profile shows two strong maxima that track the well-known active-region belts on either side of the solar equator. These peaks, typically centered in the low and mid-latitudes where sunspot groups preferentially emerge, remain strikingly similar in shape and position from one cycle to the next. Even between weak cycles and strong cycles, and between even-numbered and odd-numbered cycles with their differing amplitudes and parity, the latitude profiles of the corona–sunspot relationship barely change, indicating that the large-scale coupling between photospheric activity and coronal emission is largely independent of how vigorous any particular cycle happens to be.

This finding carries weight because solar cycles are far from identical twins. Cycle 19, which peaked in the late 1950s, was the most active in the observational record, while Cycle 24 was notably feeble, and the even–odd alternation in cycle strength, sometimes called the Gnevyshev–Ohl rule, hints at deep-seated differences in the solar dynamo from one cycle to the next. Yet the corona appears to respond to the sunspot engine in a nearly invariant geometric pattern. The double-peaked structure is a direct signature of the fact that green coronal emission is strongest not randomly across the disk but precisely where active-region magnetic fields cluster, confirming quantitatively what decades of qualitative observation have suggested: the corona is, in a very real sense, a magnetic portrait of the activity belts painted in glowing iron light.

Perhaps the most subtle part of the analysis concerns timing. Using lag-correlation techniques, the study finds that the strongest correlations between the coronal emission and the sunspot activity field generally occur at positive lags, which would naively suggest that the corona follows the sunspots with a specific, measurable delay. But here the author applied a rigorous statistical safeguard. Surrogate-data tests, in which the temporal ordering of one record is scrambled while its statistical properties are preserved, show that these formal lag maxima are usually not statistically distinguishable from the correlations at neighboring lag values. The lag-correlation profiles instead display broad, flat positive-lag plateaus, meaning the apparent delay is not a sharply defined physical delay at all.

What the plateaus do reveal is something arguably more interesting: temporal persistence and memory in the large-scale coronal magnetic structures. Rather than the corona slavishly tracking each new burst of sunspot emergence after a fixed interval, the green-line emission behaves like a slowly evolving reservoir that integrates and smooths the photospheric activity underneath it. Magnetic fields introduced into the corona by emerging active regions persist, relax, and redistribute over extended timescales, so the coronal response is broad and smeared rather than impulsive. This interpretation aligns with a growing body of work emphasizing that the large-scale distribution of magnetic flux, spanning scales far beyond individual sunspots, governs the evolution of the corona as a global system.

The comparison between modeling approaches proved decisive in this regard. When the Gaussian activity-field representation was used as the predictor of coronal emission, it produced substantially stronger and more coherent correlations than conventional sunspot number or sunspot area measures ever achieved. In other words, knowing not just how many sunspots exist but where they sit on the disk at any moment allows a far better reconstruction of the green-line corona’s behavior. This supports the central physical conclusion of the paper: that it is the large-scale spatial distribution of magnetic activity, rather than the presence of individual sunspots alone, that controls the evolution of the large-scale corona. For modelers of the solar atmosphere, this is an invitation to move beyond scalar activity indices toward genuinely two-dimensional descriptions of the emerging flux.

The practical implications extend beyond pure solar physics. The green corona has long served as a proxy for solar activity in space-weather and solar-terrestrial research, and the coronal index built from it has been used to study everything from geomagnetic disturbance to long-term solar variability. Establishing that the corona’s response to sunspot activity is stable in latitude, amplitude-independent, and temporally diffuse strengthens confidence in using coronal records to extend our knowledge of solar behavior backward in time and, potentially, to refine reconstructions of past solar cycles. It also cautions against overinterpreting any single lag measurement as a physical delay mechanism, a lesson in statistical humility that applies across solar-terrestrial science.

Spanning seven decades and seven solar cycles, this analysis demonstrates that beneath the Sun’s famously variable surface behavior lies an unexpectedly disciplined architecture. The green-line corona and the sunspot belts are locked together in a latitude-dependent partnership whose geometry survives even the wildest swings in cycle strength, while their timing relationship reflects the lingering memory of magnetic structures rather than any simple cause-and-effect delay. As the Sun progresses through Cycle 25 and observers accumulate fresh coronal data, the stable patterns documented here provide a robust baseline against which any genuine, long-term change in solar behavior can be recognized. In the quiet green glow of iron atoms two million degrees hot, the Sun has written down, cycle after cycle, a remarkably consistent record of its own magnetic heartbeat.

Subject of Research: Latitude-dependent relationships between green-line coronal emission and sunspot activity across Solar Cycles 18–24

Article Title: Large-Scale Latitude-Time Relationships Between the Green-Line Corona and Sunspot Activity During Solar Cycles 18 – 24

Article References: Large-Scale Latitude-Time Relationships Between the Green-Line Corona and Sunspot Activity During Solar Cycles 18 – 24. (n.d.). https://doi.org/10.1007/s11207-026-02729-8

Image Credits: AI Generated

DOI: 10.1007/s11207-026-02729-8

Keywords: solar corona, green coronal line, sunspots, solar cycle, solar magnetic fields, activity belts, lag correlation, surrogate data, solar physics, coronal index, solar dynamo, space weather