Corotating Interaction Regions (CIRs)

• Corotating Interaction Regions (CIRs) are large-scale solar wind structures formed when fast winds overtake slower flows, compressing plasma and magnetic fields.
• They exhibit distinct internal regions with varied thermodynamic behaviors and polytropic indices that reveal differences in plasma heating mechanisms.
• CIRs drive recurrent geomagnetic storms, energetic particle acceleration, cosmic-ray modulation, and influence planetary magnetospheres and habitability.

Corotating Interaction Regions (CIRs) are large-scale, long-lived structures in the solar wind produced when a fast solar-wind stream, typically emanating from a coronal hole, overtakes slower wind emitted earlier from the Sun. The interaction compresses and heats the plasma, enhances density and magnetic field strength, and generates a recurrent heliospheric structure tied to solar rotation, so that comparable fast–slow stream patterns can sweep past $1\ \mathrm{au}$ on a timescale of about 27 days. In the inner and outer heliosphere, CIRs are a principal framework for describing stream interfaces, compression regions, shock formation, recurrent geomagnetic activity, energetic-particle acceleration, and the evolution of large-scale solar-wind structure into merged interaction regions (MIRs) (Dayeh et al., 19 Feb 2025).

1. Origin, morphology, and heliospheric setting

CIRs arise from the interaction of fast and slow streams that follow Parker-spiral trajectories. In the solar case, the fast stream is associated with long-lived coronal holes, whereas the slow wind is emitted earlier and is overtaken at larger heliocentric distance. Near $1\ \mathrm{au}$, the large-scale CIR structure is commonly described by a forward boundary on the slow-wind side, a stream interface separating compressed slow and compressed fast wind, a reverse boundary on the trailing fast-wind side, and the end of the high-speed stream farther downstream. Characteristic in-situ signatures include compression of density and magnetic field, plasma heating, a drop in density across the stream interface, an increase in temperature, and a small but clear jump in bulk flow speed (Dayeh et al., 19 Feb 2025).

The radial evolution of CIRs is central to their dynamics. They can be observed from about $0.17\ \mathrm{au}$ to beyond Jupiter’s orbit, and farther out multiple CIRs and other transient structures can merge into MIRs that shape the outer heliosphere (Dayeh et al., 19 Feb 2025). In the standard heliospheric picture, sufficiently strong compressions steepen with distance into a forward shock propagating into the slow wind and a reverse shock propagating into the fast wind; shocks are usually seen beyond $1\ \mathrm{au}$, although they can form as close as about $0.7\ \mathrm{au}$ (Wu et al., 2014). In a young-Sun analog, the fast–slow contrast is much larger, and the resulting SIR compression steepens into strong shocks with $v_{\rm shock} \gtrsim 1000\ \mathrm{km\ s^{-1}}$ at the orbital distances of Mercury, Venus, and Earth (Sur et al., 25 Sep 2025).

2. Internal regions and thermodynamic behavior

A widely used partition of CIRs distinguishes four physically distinct regions bounded by the forward boundary (FB), stream interface (SI), reverse boundary (RB), and the end of the high-speed stream (HSS). This segmentation is particularly useful for thermodynamic analysis because the compressed and uncompressed parts of the slow and fast wind behave differently.

Region	Interval	Qualitative state
Uncompressed slow solar wind	1 day before FB	Lower speed, higher density, lower temperature, weak Alfvénicity
Compressed slow solar wind	FB to SI	Elevated density, temperature, and magnetic-field magnitude
Compressed fast solar wind	SI to RB	Elevated density, temperature, and magnetic-field magnitude
Uncompressed fast solar wind	RB to HSS end	High speed, relatively lower density, elevated temperature, strong Alfvénic fluctuations

The thermodynamics of these regions have been analyzed using the polytropic relation

$$P\,n^{-\gamma} = T\,n^{1-\gamma} = \mathrm{const.},$$

or, equivalently,

$\ln T_p = (\gamma-1)\ln n_p + \mathrm{constant}.$

Using Wind/SWE and MFI data at 96-second cadence, bi-Maxwellian proton moments, 5-minute sliding windows, a Bernoulli-integral stability filter with variance/mean below 10%, and weighted linear fits in $\ln T_p$–$\ln n_p$ space, a 117-event sample of CIRs yielded distinct regional polytropic indices. The mean values were approximately $1\ \mathrm{au}$0 for uncompressed slow wind, $1\ \mathrm{au}$1 for compressed slow wind, $1\ \mathrm{au}$2 for compressed fast wind, and $1\ \mathrm{au}$3 for uncompressed fast wind. These results imply nearly adiabatic behavior upstream in uncompressed slow wind, sub-adiabatic behavior and strong net heating in compressed slow wind, near-adiabatic behavior in compressed fast wind, and the largest super-adiabatic deviation in the trailing uncompressed fast wind (Dayeh et al., 19 Feb 2025).

The same study linked these thermodynamic differences to Alfvénic activity. For an individual CIR example, the correlation coefficient between $1\ \mathrm{au}$4 and $1\ \mathrm{au}$5 was about $1\ \mathrm{au}$6 in uncompressed slow wind, $1\ \mathrm{au}$7 in compressed slow wind, $1\ \mathrm{au}$8 in compressed fast wind, and $1\ \mathrm{au}$9 in uncompressed fast wind; over all 117 events, the downstream uncompressed fast wind showed correlation distributions strongly peaked near $0.17\ \mathrm{au}$0. This pattern places CIRs in a two-mechanism heating framework: compression dominates in compressed slow wind, whereas Alfvénic wave dissipation associated with fast streams strongly modifies the effective polytropic behavior downstream (Dayeh et al., 19 Feb 2025).

3. Energetic particles, composition, and cosmic-ray modulation

CIRs are major acceleration sites for energetic ions. At sufficiently large heliocentric distance, the forward and reverse shocks accelerate particles from solar-wind and suprathermal seed populations, and the resulting energetic-ion abundances mirror fast-solar-wind composition modified by a decreasing power law in mass-to-charge ratio, $0.17\ \mathrm{au}$1. In twelve Wind/LEMT events, the source or freezing-in temperatures inferred from abundance fits were $0.17\ \mathrm{au}$2 or $0.17\ \mathrm{au}$3, typical of the fast solar wind, and the fitted exponents were nearly always negative, ranging from about $0.17\ \mathrm{au}$4 to $0.17\ \mathrm{au}$5. One unusual 1982 event yielded $0.17\ \mathrm{au}$6, interpreted as a case where the reverse shock was at about $0.17\ \mathrm{au}$7, so the energetic-ion abundances directly matched fast-solar-wind abundances (Reames, 2018).

Transport strongly shapes the spectra observed at $0.17\ \mathrm{au}$8. A focused-transport analysis of a 2008 February 08 CIR event, observed by both ACE and STEREO-B, found a reverse shock at STEREO-B ($0.17\ \mathrm{au}$9) but not at ACE ($1\ \mathrm{au}$0). Assuming a stationary CIR in the corotating frame and transport strictly along Parker-spiral field lines without cross-field diffusion, the inferred shock-connection radii were about $1\ \mathrm{au}$1 for ACE and $1\ \mathrm{au}$2 for a later STEREO-B interval; the study concluded that modulation of sub-$1\ \mathrm{au}$3 particles is significant and that reasonable agreement with both spacecraft required the CIR shock to accelerate more particles at larger heliocentric distance than at smaller heliocentric distance (Zhao et al., 2017).

Recent modeling has clarified how perpendicular diffusion alters the expected $1\ \mathrm{au}$4 dependence. In a synthetic CIR built with EUHFORIA and followed with a focused-transport model, ion spectra without perpendicular diffusion showed a strong $1\ \mathrm{au}$5 dependence, whereas inclusion of perpendicular diffusion reduced the scaling to $1\ \mathrm{au}$6, and to $1\ \mathrm{au}$7 when $1\ \mathrm{au}$8 was taken to be $1\ \mathrm{au}$9-independent. This qualitatively agrees with the weak $0.7\ \mathrm{au}$0 dependence often inferred from CIR observations near $0.7\ \mathrm{au}$1 (Ding et al., 2024).

CIRs also modulate galactic cosmic rays. An analytic convection–diffusion model for recurrent Forbush decreases argued that two effects are sufficient to explain the observed profiles: enhanced convection caused by the increased velocity of the high-speed stream, and reduced diffusion caused by the enhanced magnetic field and its fluctuations within the CIR and trailing high-speed stream. Using a quiet radial diffusion coefficient $0.7\ \mathrm{au}$2 and a quiet radial gradient $0.7\ \mathrm{au}$3, the model reproduced the main depression, over-recovery, and secondary dips for a long-lived CIR that recurred over 27 rotations in 2007–2008 (Vrsnak et al., 2022).

4. Geomagnetic, thermospheric, and planetary consequences

At Earth, CIRs are established drivers of recurrent space weather. They can produce strong enhancements in solar-wind dynamic pressure, sometimes forming shocks near $0.7\ \mathrm{au}$4, compress the dayside magnetopause, and sustain prolonged geomagnetic activity with $0.7\ \mathrm{au}$5 up to about $0.7\ \mathrm{au}$6–$0.7\ \mathrm{au}$7. Repeated CIRs contribute to long-duration disturbances in the magnetosphere, affecting aurorae, the ring current, the ionosphere–thermosphere system, and spacecraft charging (Dayeh et al., 19 Feb 2025).

The exceptionally quiet solar minimum of 2008 isolated these effects especially clearly. During that interval, recurrent high-speed streams from near-equatorial coronal holes were the primary contributor to geomagnetic activity, and 38 high-speed streams/CIRs were identified in 2008. Superposed-epoch analysis of CHAMP accelerometer densities normalized to $0.7\ \mathrm{au}$8 showed that thermospheric density changed by 75% on average, with a nighttime increase of 85% on average during the main phase; the relative changes were comparable at different latitudes, although the variability was largest at high latitudes. Lomb–Scargle spectra further showed strong 9- and 13.5-day periodicities in solar-wind parameters, $0.7\ \mathrm{au}$9, and thermospheric density, linked to the spatial distribution of low- to middle-latitude coronal holes (Lei et al., 2010).

For a young solar-type star, modeled CIRs become substantially more extreme. In a $v_{\rm shock} \gtrsim 1000\ \mathrm{km\ s^{-1}}$0 Ceti-based young-Sun simulation, the CIR compression region at $v_{\rm shock} \gtrsim 1000\ \mathrm{km\ s^{-1}}$1 had $v_{\rm shock} \gtrsim 1000\ \mathrm{km\ s^{-1}}$2 between $v_{\rm shock} \gtrsim 1000\ \mathrm{km\ s^{-1}}$3 and $v_{\rm shock} \gtrsim 1000\ \mathrm{km\ s^{-1}}$4, proton density between $v_{\rm shock} \gtrsim 1000\ \mathrm{km\ s^{-1}}$5 and $v_{\rm shock} \gtrsim 1000\ \mathrm{km\ s^{-1}}$6, and dynamic pressure peaking at $v_{\rm shock} \gtrsim 1000\ \mathrm{km\ s^{-1}}$7. The resulting early-Earth magnetospheric response included a minimum dayside magnetopause standoff distance of $v_{\rm shock} \gtrsim 1000\ \mathrm{km\ s^{-1}}$8, a simulated minimum SYM-H of $v_{\rm shock} \gtrsim 1000\ \mathrm{km\ s^{-1}}$9, peak cross-polar-cap potentials of $$P\,n^{-\gamma} = T\,n^{1-\gamma} = \mathrm{const.},$$0 in the Northern Hemisphere and $$P\,n^{-\gamma} = T\,n^{1-\gamma} = \mathrm{const.},$$1 in the Southern Hemisphere, and hemispheric Joule-heating maxima of $$P\,n^{-\gamma} = T\,n^{1-\gamma} = \mathrm{const.},$$2 and $$P\,n^{-\gamma} = T\,n^{1-\gamma} = \mathrm{const.},$$3, respectively. In that event, proton density, not speed, was the dominant contributor to the dynamic pressure and to the ensuing electron precipitation and Joule heating (Sur et al., 25 Sep 2025).

5. Radial evolution, stellar-age dependence, and habitability

CIRs are not static with respect to heliocentric distance or stellar age. In the present-day heliosphere, their compressions steepen with radial distance, shocks usually strengthen outward, and multiple interaction regions merge into MIRs beyond the inner heliosphere (Dayeh et al., 19 Feb 2025). Extending this idea over stellar evolution, a rotational-evolution study of a solar-mass star found that CIRs form closer to the star during early spin-up phases and migrate outward during spin down, with the minimum CIR radius inversely related to the stellar rotation rate (Waugh et al., 22 Jun 2026).

That same work argued that the energetic-particle output of CIR shocks changes markedly with age. During the Hadean period, CIRs may have generated a number of energetic particles that is $$P\,n^{-\gamma} = T\,n^{1-\gamma} = \mathrm{const.},$$4 to $$P\,n^{-\gamma} = T\,n^{1-\gamma} = \mathrm{const.},$$5 times greater than for the present-day Sun, and the frequency and strength of CIR–planet interactions peak during early rapid-rotation phases. The study further concluded that, for stars with mass less than $$P\,n^{-\gamma} = T\,n^{1-\gamma} = \mathrm{const.},$$6, CIRs can form within the habitable zone whereas their shocks always form beyond it, implying that energetic-particle impacts may propagate inward from more distant regions (Waugh et al., 22 Jun 2026).

Young-Sun magnetospheric modeling adds temporal context to this picture. For the $$P\,n^{-\gamma} = T\,n^{1-\gamma} = \mathrm{const.},$$7 Ceti analog, 3D MHD simulations yielded about 3–4 CIR events per 9.2-day stellar rotation, corresponding to roughly 11–14% of the time spent in high-dynamic-pressure compressed states at $$P\,n^{-\gamma} = T\,n^{1-\gamma} = \mathrm{const.},$$8; the authors suggested that even younger, faster rotators could place early Earth analogs inside CIR structures more than 30% of the time (Sur et al., 25 Sep 2025). This suggests that CIRs can be treated not merely as recurrent solar-wind compressions but as long-lived agents of atmospheric erosion, upper-atmospheric heating, and magnetospheric forcing over stellar evolutionary timescales.

6. Terminology, forecasting, and broader usage

The terminology of interaction regions remains contested. A recent review traced the sequence from Belcher and Davis’s “interaction region (IR)” to Smith and Wolfe’s “corotating interaction region (CIR)” and Jian et al.’s “stream interaction region (SIR).” That review argued against using “SIR” as a broad label for transient and possibly localized stream interactions with poor recurrence, and instead recommended using specific names for distinct interplanetary structures such as CIRs, high-speed streams, ICME sheaths, and ICMEs. It also introduced “Super CIR (SCIR)” for a CIR associated with magnetic reconnection at the edge of a solar coronal hole with an embedded coronal jet; the SCIR of 6–7 April 2000 was reported to have exceptionally strong internal magnetic fields, both forward and reverse shocks, and to have caused a SYM-H $$P\,n^{-\gamma} = T\,n^{1-\gamma} = \mathrm{const.},$$9 superstorm (Tsurutani et al., 6 Jun 2026).

Operationally, CIR forecasting has become a distinct modeling problem because these structures are recurrent and often dominate ambient solar-wind variability. The “Helio1D” pipeline couples Multi-VP at $\ln T_p = (\gamma-1)\ln n_p + \mathrm{constant}.$0 to a 1D ideal-MHD propagation model and forecasts ambient solar-wind conditions at L1 with a lead time of 4 days, explicitly targeting CIRs and high-speed streams. The modeled CIR signature at $\ln T_p = (\gamma-1)\ln n_p + \mathrm{constant}.$1 is a transition from slow to fast wind in $\ln T_p = (\gamma-1)\ln n_p + \mathrm{constant}.$2, often accompanied by a polarity change in $\ln T_p = (\gamma-1)\ln n_p + \mathrm{constant}.$3 and enhancements in density and temperature near the stream interface. An ensemble of 21 virtual targets is used to estimate timing and magnitude uncertainty, and Fast Dynamic Time Warping is used to compare modeled and observed CIR/HSS structures. In the forecast statistics, fast-wind arrivals are typically early by 10–20 hours and fast-stream speeds are often underestimated by 50–100 $\ln T_p = (\gamma-1)\ln n_p + \mathrm{constant}.$4 (Kieokaew et al., 2023).

Outside heliophysics, the term CIR is also used in the winds of massive stars. In apparently single Galactic Wolf–Rayet stars, a systematic survey found that 23 of 68 stars showed large-scale variability, but only 12 of 54 stars were potentially of CIR-type, corresponding to about 22.1%; in that literature, CIRs denote large-scale spiral wind structures tied to stellar rotation and diagnosed through coherent spectroscopic or polarimetric variability (Chené et al., 2011). This broader usage does not alter the solar-wind meaning of CIRs, but it underscores that corotating stream interactions are a generic feature of structured rotating outflows.