Issue 
A&A
Volume 658, February 2022



Article Number  A186  
Number of page(s)  15  
Section  The Sun and the Heliosphere  
DOI  https://doi.org/10.1051/00046361/202140846  
Published online  24 February 2022 
Analytic modeling of recurrent Forbush decreases caused by corotating interaction regions
^{1}
Hvar Observatory, Faculty of Geodesy, University of Zagreb, Kačićeva 26, 10000 Zagreb, Croatia
email: bvrsnak@geof.hr
^{2}
Institut für Experimentelle und Angewandte Physik, ChristianAlbrechtsUniversität zu Kiel, ChristianAlbrechtsPlatz 4, 24118 Kiel, Germany
^{3}
Karlovac University of Applied Sciences, Trg J.J. Strossmayera 9, 47000 Karlovac, Croatia
Received:
22
March
2021
Accepted:
25
October
2021
Context. On scales of days, the galactic cosmic ray (GCR) flux is affected by coronal mass ejections and corotating interaction regions (CIRs), causing socalled Forbush decreases and recurrent Forbush decreases (RFDs), respectively.
Aims. We explain the properties and behavior of RFDs recorded at about 1 au that are caused by CIRs generated by solar wind highspeed streams (HSSs) that emanate from coronal holes.
Methods. We employed a convectiondiffusion GCR propagation model based on the FokkerPlanck equation and applied it to solar wind and interplanetary magnetic field properties at 1 au.
Results. Our analysis shows that the only two effects that are relevant for a plausible overall explanation of the observations are the enhanced convection effect caused by the increased velocity of the HSS and the reduced diffusion effect caused by the enhanced magnetic field and its fluctuations within the CIR and HSS structure. These two effects that we considered in the model are sufficient to explain not only the main signatures of RFDs, but also the sometimes observed “overrecovery” and secondary dips in RFD profiles. The explanation in terms of the convectiondiffusion GCR propagation hypothesis is tested by applying our model to the observations of a longlived CIR that recurred over 27 rotations in 2007–2008.
Conclusions. Our analysis demonstrates a very good match of the model results and observations.
Key words: Sun: heliosphere / solarterrestrial relations / solar wind / magnetohydrodynamics (MHD) / methods: analytical
© ESO 2022
1. Introduction
Recurrent high speed streams (HSSs) in the solar wind that originate from coronal holes cause various important effects in the heliosphere (e.g., Nolte et al. 1976; Schwenn 2006; Tsurutani et al. 2006; Vršnak et al. 2007; Rotter et al. 2012; Hofmeister et al. 2018, and references therein); for a review, see Cranmer (2009). The interaction of HSSs with the preceding slow solar wind creates the socalled corotating interaction regions (CIRs; e.g., Gosling & Pizzo 1999; Schwenn 2006; Tsurutani et al. 2006, and references therein); for a review, see Richardson (2018). These HSS and CIR structures cause recurrent moderate geomagnetic storms and substorms, where the Disturbance StormTime index (Dst) rarely goes below −100 nT (e.g., Tsurutani et al. 2006; Vršnak et al. 2007, 2017; Zhang et al. 2007; Verbanac et al. 2011a,b, 2013), as well as temporal depressions in the galactic cosmic ray (GCR) count rate, which are generally known as recurrent Forbush decreases (RFDs; for a review, see, e.g., Iucci et al. 1979; Kunow et al. 1995; Richardson 2004, 2018). Generally, Forbush decreases (FDs; Forbush 1937) are temporal, welldefined decreases in the GCR count rate that are caused by interplanetary coronal mass ejections (ICMEs) and CIR structures (e.g., Lockwood 1971; Cane 2000; Richardson 2004, 2018; Dumbović et al. 2011, 2012, and references therein); for reviews, see Cane (2000) and Richardson (2004). We note here that some authors exclusively use the term Forbush decrease for GCR depletions caused by ICMEs, and those that are caused by CIRs, they call “recurrent GCR depressions”. RFDs are usually relatively weak, and their amplitudes are considerably lower than those of “standard” FDs caused by ICMEs (e.g., Richardson 2004; Čalogović et al. 2009; Dumbović et al. 2011, 2012; Badruddin & Kumar 2016; Melkumyan et al. 2019, and references therein).
In recent years, a number of studies related to RFDs were published. For example, Kumar (2014), Badruddin & Kumar (2016), Ghanbari et al. (2019), and Guo et al. (2021) studied the behavior of RFDs by applying the superposed epoch analysis approach. We also refer to the observationbased studies by Modzelewska et al. (2020) and Shen et al. (2020), as well as to the numerical studies by Alania et al. (2011), Guo & Florinski (2016), Kopp et al. (2017), and Luo et al. (2020).
The interaction of the highspeed stream with the preceding slow solar wind creates the CIR structure by compressing the plasma in the region in which they interact (Gosling & Pizzo 1999; Jian et al. 2006; Richardson 2018). The compression region includes parts of the slow solar wind and parts of the HSS, divided by the socalled stream interface. The compressed interaction region is characterized by a continuously increasing flow speed, consistent with the continuity equation. The compression region is followed by the rarefaction region, in which the flow velocity gradually decreases, which is again consistent with the continuity equation. The whole CIR and HSS structure (especially the compression region) is characterized by significantly increased magnetic field fluctuations that are mainly related to the Alfvénicwave perturbations.
Transport of GCRs in the heliosphere is usually described in terms of the socalled Parker GCR transport equation (Parker 1965). It considers the GCR particle density at a given time and position, which is determined by the effect of the convection and diffusion (i.e., two terms from the standard FokkerPlanck equation) and by a term added by Parker, describing the effect of adiabatic cooling of particles within an expanding volume (or heating in the case of compression). Later on, the Parker GCR transport equation was employed in many studies concerning various aspects of the heliospheric propagation and behavior of GCRs (e.g., Gleeson & Axford 1968; Perko & Fisk 1983; Jokipii et al. 1993; Fisk et al. 1998; Potgieter 1998; Wibberenz et al. 1998; Wibberenz & Cane 2000; MoralesOlivares & CaballeroLopez 2010; Gieseler & Heber 2016; Zimbardo et al. 2017; Richardson 2018, and very many references therein).
The diffusion rate is determined by the diffusion coefficient, which depends on the meanfree path of particles, where the scattering is dominantly caused by the magnetic field fluctuations (Parker 1965). While the diffusion coefficients are calculated, the corresponding total fluctuating magnetic energy is usually computed using the root mean square variation in the vector of the interplanetary magnetic field (e.g., Zank et al. 1998). The particle drift effect can also be included in the transport equation (e.g., Fisk et al. 1998, and references therein), as well as the effect of mirroring on steep magnetic field gradients and shocks (e.g., Kirin et al. 2020, and references therein). The diffusion regulates the inward transport of particles, whereas the solar wind flow regulates the outward transport of GCRs. The adiabatic term describes changes in the GCR energy spectrum.
This paper aims to explain the behavior of recurrent Forbush decreases caused by CIR and HSS structures in the solar wind at heliospheric distances about 1 au, based on a convectiondiffusion approach. Global 3D timedependent simulations obtained by solving the GCR transport equation have shown that the smaller diffusion coefficient in the regions of stronger magnetic field related to CIRs results in a rapid local decrease in cosmic ray flux (Kota & Jokipii 1991). We also note that using the convectiondiffusion approach to explain the CIRrelated GCR modulation is not new (e.g., Chih & Lee 1986; Wibberenz et al. 1998). Richardson et al. (1996) considered the steadystate diffusionconvection model, which includes adiabatic deceleration and longitudinal variations in the solar wind speed. In this model, the effect of the convection due to the outward flow of the solar wind is locally balanced by the inward diffusion of GCRs. The model that we present in this paper is to a certain degree similar to the more sophisticated 3D numerical model by Alania et al. (2011). They regarded a steadystate heliosphere with CIR represented as the heliolongitudinal variation in solar wind speed and applied a full Parker transport equation. They found that the variation in GCR intensity is inversely correlated with the modulation parameter, which is proportional to the product of the solar wind velocity and the strength of the interplanetary magnetic field. Moreover, their calculations have shown that the heliospheric current sheet (HCS) structure in the model does not significantly affect the results. A similar result was also found by Guo & Florinski (2016) using the 3D coupled MHDCR transport model. They found that GCR variations depend on the ratio of diffusion coefficients in the fast and slow solar winds. Finally, we compare our results to the most recent dedicated 3D coupled MHDCR transport models by Kopp et al. (2017) and Luo et al. (2020). Their simulations show that enhanced convection, reduced diffusion, and drifts due the magnetic field enhancement may contribute to the recurrent FD. However, the simulations cannot distinguish the relative contribution of different mechanisms.
In this respect we note that an application of the modeling results to observations and the corresponding quantitative analysis is still lacking. Hereafter, the GCR count rate measured in situ by spacecraft (s/c) is denoted as CR, whereas for the analog in the model option, the symbol Φ is applied.
For the purpose of our analysis, we use some observational examples related to a longliving CIR structure that was recorded over 27 solar rotations in the period from June 2007 to May 2009. It was analyzed in detail by Dumbović et al. (2022, hereafter, Paper I; see also Kühl et al. 2013). In Sect. 2 the observational aspect is described in detail, and in Sect. 3 the analytical model is proposed. A comparison of model results with observations is given in Sects. 4 and 5. The outcome of the analysis is discussed in Sect. 6, and conclusions are drawn in Sect. 7.
2. CIR structure and the related GCR flux modulation
In Fig. 1, two examples of welldefined CIRs are presented. The events were recorded in the periods between the day of the year DOY = 184 and 191, 2007 (Carrington rotation 2059; Fig. 1a), and DOY = 58 and 68, 2008 (Carrington rotation 2067; Fig. 1b). Both CIR episodes were related to a longliving equatorial coronal hole that produced a persistent CIR and HSS that lasted for 27 rotations. The presented examples were recorded during its second and eleventh rotation (hereafter, rot2 and rot11, respectively.
Fig. 1. 1 au measurements of CIRs recorded from DOY 183 to 191, 2007 (rot2; left column) and from DOY 58 to 68, 2008 (rot11; right column). Panels a and e: flow speed, V. Panels b and f: magnetic field strength, B (red), and fluctuations (root mean square), δB (blue). Panels c and g: proton number density, N (red), and temperature, T (blue). Panels d and h: EPHIN cosmic ray countrate change, CR, expressed in percent of the preevent count rate. The dashed rectangle in both panels marks the period of the flow speed increase. 
In Figs. 1a–c and e–g, the 1 h resolution data from the OMNI database^{1} (King & Papitashvili 2005) are displayed, and Figs. 1d and h show the 1 h resolution data from the Electron Proton Helium Instrument (EPHIN; MüllerMellin et al. 1995) on board the Solar and Heliospheric Observatory (SOHO), measured by its Fdetector, which records hydrogen and helium nuclei above 53 MeV (MüllerMellin et al. 1995). However, it is important to note that although the lower cutoff of SOHOEPHINF is about 53 MeV, the response function of SOHOEPHIN is such that it increases with energy, that is, it is less sensitive to lowenergy particles (Dumbović et al. 2020).
The rot2 event began on DOY = 184.3, 2007, starting from a quietSun wind speed of V ≈ 340 km s^{−1} (Fig. 1a). The HSS speed attained a maximum of ≈610 km s^{−1} at DOY = 185.3, after which it gradually decreased to only ≈305 km s^{−1} on DOY = 190.8. The phase of the rising speed is marked in panels a–d of Fig. 1 by the dashed blue rectangle. The period of increasing flow speed corresponds to the enhanced magnetic field, B, and plasma density, N (red graphs in Figs. 1b and c, respectively) caused by the interaction of the HSS with the ambient slow solar wind. This period is also characterized by enhanced magnetic field fluctuations, δB, and rising temperature, T (blue graphs in Figs. 1b and c, respectively). The change of the EPHIN cosmic ray count rate, CR, expressed in percentiles of the preevent value, is presented in Fig. 1d. The CR(t) graph very closely corresponds to the inverted V(t) graph (for more examples, see Paper I). The CR data show a clear “overrecovery” phenomenon (e.g., Jämsén et al. 2007; Dumbović et al. 2012) after the main CR dip.
The rot11 CIR began on DOY = 58.5, 2008, starting from a speed of V ≈ 350 km s^{−1} (Fig. 1e). The HSS speed attained a maximum of ≈780 km s^{−1} at DOY = 61.2, after which it gradually decreased to ≈320 km s^{−1} on DOY = 68. The phase of the rising speed is marked in panels e–h of Fig. 1 by the dashed blue rectangle. The period of increasing flow speed again corresponds to the enhanced magnetic field, B, and plasma density, N (red graphs in Figs. 1f and g, respectively), as well as to the enhanced magnetic field fluctuations, δB, and rising temperature, T (blue graphs in Figs. 1f and g, respectively). The change in the EPHIN cosmic ray count rate is presented in Fig. 1h. Again, the CR(t) graph very closely corresponds to the inverted V(t) graph. The data presented in Figs. 1f–g show a distinct secondary peak in the B(t), δB(t), and N(t) graphs between DOY ≈ 65 and 66, corresponding to a weak increase in flow speed that is recognizable in the V(t) graph (Fig. 1e), and is also notable as a weak new rise in the T(t) graph (Fig. 1g). At the same time, CR(t) shows a distinct secondary decrease, with a minimum at DOY = 65.7 (Fig. 1h).
In Fig. 2a, a simplified generic temporal behavior of CIR structure and the associated RFD is sketched schematically (based on the detailed statistical study presented in Paper I). The abscissa shows the mean value of the relative time differences, Δt, of the onsets, maxima, and ends of V, N, B, and CR relative to the onset of the flowspeed increase. It is expressed in days (thus, Δt = 0 for the V onset). The yaxis in the y > 0 range represents the mean CIRdisturbance amplitudes, normalized with respect to the preCIR values, based on the mean values of the amplitudes of N, B, and V (to highlight the relative amplitude of V, it is multiplied by a factor of two). The [x_{i}, y_{i}] points obtained for the mentioned key parameters, that is, onsets, maxima, and so on, are finally connected by straight lines (this is just a schematic drawing). The graph includes the rarefaction region at the rear of HSS. In the y < 0 range, the mean RFD amplitude is shown as the percentage of the preCIR level.
Fig. 2. “Generic” temporal behavior of the flow speed, V (blue), magnetic field, B (red), density, N, (green), and the cosmic ray countrate change expressed in percent, CR (black), based on the study of a sample of CIRs presented in Paper I. Panel a: simplified, generic, presentation based on the mean values of the amplitudes and timing of the B, N, and CR onsets, peaks, and ends, relative to the onset of V (to see better the relative amplitude of V, it is multiplied by a factor of two). The peak times are marked by vertical dashed lines. Panel b: SEA results; the variations relative to the preCIR values are normalized to the peak values, ΔV/V_{max}, ΔB/B_{max}, ΔN/N_{max}, and ΔCR/CR_{max}. The time Δt = 0 is set at the onset of the V increase. 
The plasma density starts to increase before the onset of the flowspeed rise on average, whereas the magnetic field in this representation starts to increase approximately simultaneously with the flowspeed increase. After the compression period, a rarefaction phase occurs, characterized by weakening of the magnetic field and plasma density decrease (e.g., Hundhausen 1973; Gosling 1996; Gosling & Pizzo 1999; Richardson 2018, and references therein). The CR decrease occurs during the period of the speed increase, corresponding to the phase of the enhanced magnetic field. Before the CR decrease, as is quite often observed, a weak enhancement in the count rate is seen in the period when only the density is increased, like in the superepoch analysis (SEA) representation (Fig. 2b). The CR recovery phase closely follows the decrease in the flow speed. Regarding the SEA approach, it should be mentioned that Badruddin & Kumar (2016) analysed the differences between characteristics of CIRs that developed a frontal shock and those without the shock. We emphasize that in our sample, CIRs did not show a shock signatures, and our SEA results are quite similar to those presented by Badruddin & Kumar (2016).
Figure 2b shows an analogous, but more detailed, statisticsbased generic time profile of V, B, N, and CR, where the displayed results are based on the SEA, performed in Paper I, covering the complete set of 27 rotations. Unlike the schematic presentation shown in Fig. 2a, the graph displayed in Fig. 2b reveals that not only does N starts to increase before V, but so does B (but within the standard deviation of Δt). Moreover, the graph V(t) reveals a slight depletion before the main rise of V. The graph CR(t) shows the preCIR slight increase in CR in more detail (e.g., Lockwood 1971), which is hereafter denoted as “nose”. As we show below, the nose is most likely a consequence of the preevent phase of the decreased V and the steep magnetic field gradient related to the compression caused by the interaction of the fast and slow solar wind (Kirin et al. 2020).
The characteristics of the CIR structure and the associated RFD behavior presented in Figs. 1 and 2 are the basis for constructing a model for the CIRrelated RFDs. It is also used for a direct comparison of the observations and model results.
3. Model
3.1. Physical background
In the following, we present a simple model for CIR and HSS effects on the cosmic ray count rate in the heliocentric distance range of r = 1 au (the global situation is sketched and the coordinate system explained in Figs. 3a and b). The model describes a single rotation, that is, rotations are considered independent (differences from one to another rotation are due only to different quietSun values of Φ_{0}, V_{0}, and B_{0}, and the temporal behavior of Φ(t), V(t), and B(t) across the CIR and HSS at 1 au).
Fig. 3. Panel a: sketch of the global situation (upper panel: rest frame; lower panel: frame of reference rotating with the Sun); for details, see the main text. Panel b: schematic presentation of the geometry and CIR characteristics (fixed reference frame), as applied in the model (x_{0} and x_{e} represent the beginning and end of CIR and HSS at 1 au, respectively; for details, see the main text). Panel c: dependence of the radial CRgradient, G_{r}, for three specific s/c positions (x_{1}, x_{2}, x_{3}) as a function of heliocentric distance, r. The heliocentric distance of the CIR frontal edge for the three specific s/c positions, x_{i}, are denoted r_{1}, r_{2}, r_{3}. Panel d: corresponding radial distance dependence of the cosmic ray count rate, Φ_{xi}(r), between CIR and HSS outer boundary and s/c. The red line marks the behavior in a quietSun situation at the CIR outer boundary. Panel e: sketch of the resulting Forbush decrease profile, Φ(t), as measured by the s/c. 
The model is based on the FokkerPlanck diffusion equation,
which was considered by Parker (1965) as a staring point in developing a model for the heliospheric propagation of GCRs. In Eq. (1), Φ represents the GCR particle number density at a given time and position, V denotes the convective velocity, and κ_{ij} represents the diffusion coefficient tensor. Equation (1) does not consider the change in particle energy spectrum due to the adiabatic heating/cooling of particles (for a complete GCR transport equation, including the adiabatic effect, see Parker 1965). In other words, we consider in the following only the change in total particle number density, regardless of the shape of the energy spectrum.
We note that although the lower cutoff of SOHOEPHIN is about 53 MeV, the response function of SOHOEPHINF is such that it increases with energy, that is, it is less sensitive to lowenergy particles. On the other hand, the spectrum of GCRs decreases with energy. As a result, the spectrum of particles detected by EPHIN therefore looks like a belllike curve with a peak somewhere around 1 GeV (depending on IP conditions, see, e.g., Dumbović et al. 2020). Bearing in mind that at 1 au, the proton spectrum recorded by the EPHINF peaks somewhat below but close to 1 GeV, small changes in proton energy spectrum therefore do not significantly affect the spectral distribution in the energy range ≈1 GeV measured by s/c.
In the following, we consider a CIR structure in a stationary state, that is, we assume that its structure and characteristics do not change significantly over the time period required for its passage over the spacecraft. In this case, the timederivative ∂Φ/∂t = 0 at a given position within the CIR and HSS, so that Eq. (1) can be written as
that is,
which means that the effect of the convection due to the solar wind outward flow is locally balanced by the inward diffusion of GCRs. The CIRrelated decrease in Φ is a result of the suppressed inward diffusion due to the enhanced magnetic field and its fluctuations, and it is also due to the enhanced outward convection. The diffusion coefficient parallel to the field lines κ_{∥} is far larger at 1 au than the perpendicular one κ_{⊥}, that is, the diffusion perpendicular to the magnetic field is negligible compared to the diffusion throughout the field (e.g., Parker 1965; Fisk et al. 1998; Giacalone 1998; Richardson 2018). We emphasize that the value of κ_{∥} is quite ambiguous as it cannot be determined directly from observations. Its estimate therefore depends on a number of assumptions. Consequently, quite different values are found in various studies; they vary from 10^{17} to > 10^{19} m^{2} s^{−1}. For example, Parker (1965) estimated it to be 10^{17} − 10^{18} m^{2} s^{−1}, ≈3 × 10^{19} m^{2} s^{−1} was reported by Giacalone (1998), 10^{18} m^{2} s^{−1} was adopted by Gleeson & Axford (1968) and Schwenn & Marsch (1991), and Wibberenz et al. (1998) considered values 1.4 − 2.8 × 10^{18} m^{2} s^{−1}. In the following, we use as a starting point for the quietSun solar wind the radial diffusion coefficient of κ_{r} = 1.25 × 10^{18} m^{2} s^{−1}, which is obtained from the CR radial gradient at 1 au of G_{r} = 4.8%/au (Gleeson & Axford 1968) and for the solar wind speed of V = 400 km s^{−1}, and using κ_{r} = V/G_{r} (e.g., Gleeson & Axford 1968; Schwenn & Marsch 1991, see also Eq. (6)).
In constructing the model, we assumed that the diffusion coefficient is related to the magnetic field magnitude B, based on a number of previous studies (e.g., Potgieter 2013; Richardson 2018, and references therein). In particular, we considered the simplest option, κ ∝ 1/B (see, e.g., Potgieter 1998; Wibberenz et al. 1998; Wibberenz & Cane 2000). In addition, because the mean free path of particles depends on the magnetic field fluctuations (e.g., Wibberenz et al. 1998; Wibberenz & Cane 2000; Michałek & Ostrowski 2001; Shalchi & Schlickeiser 2004; Chhiber et al. 2017, and references therein), as does the diffusion coefficient, we also took the magnetic field fluctuations δB into account, considering the simplest possible option κ ∝ 1/δB. In this respect, we note that according to the model proposed by Qin (2002), which is based on standard quasilinear theory, the parallel mean free path should be proportional to B^{5/3}/δB^{2}. As we show below, this is not consistent with observations, which show that the diffusion coefficient is roughly inversely proportional to both B and δB. Under these two approximations, the perturbed value of the diffusion coefficient could therefore be estimated as B_{0}/B, or equivalently, κ_{δB}′=κ_{0}δB_{0}/δB, where κ_{0}, B_{0}, and δB_{0} represent the values in the unperturbed solar wind.
However, the relation is most likely oversimplified, that is, it is more complex than the simplest options we assumed (for other options, see, e.g., Wibberenz et al. 1998; Wibberenz & Cane 2000; Michałek & Ostrowski 2001; Shalchi & Schlickeiser 2004; Chhiber et al. 2017, and references therein). For example, the mean free path, thus also the diffusion coefficient, depends on the particle energy, meaning that from the observational point of view, it depends on the energy channel of measurements (which is not of concern for our study because we considered only one channel). The relation κ(δB) also depends on the spectrum of fluctuations. As a consequence, it probably depends on the heliocentric distance and on the solar cycle phase (in our case, this is not important as we considered a relatively short time period relative to the solar cycle duration, and we focused on heliocentric distances of about 1 au). To compensate for all these ambiguities, we introduced a free parameter k that is expected to be different from event to event due to differences in fluctuation timespatial characteristics, and it would also depend on the energy channel of the measurements.
Thus, we introduce the scaling coefficients k_{b} and k_{δb}, which we incorporated in b^{*} = k_{b}b and δb^{*} = k_{δb}δb, where we used the abbreviations b = B/B_{0} and δb = δB/δB_{0}. In other words, b^{*} and δb^{*} were used in the model calculations, where k_{b} and k_{δb} were determined by fitting the model results to the observations by minimizing the rms differences between the modeled and observed data. With this notation, the effective local radial gradient can be expressed as , where V is the local solar wind velocity, and κ^{*} is defined by the unperturbed value of κ_{0} and the effect of the perturbed magnetic field, as determined by the parameters b^{*} or δb^{*}, that is, κ^{*} = κ_{0}/b^{*} or κ^{*} = κ_{0}/δb^{*}, which means G^{*} = (V/κ_{0})b^{*} = G_{0}b^{*} and equivalently, G^{*} = (V/κ_{0})δb^{*} = G_{0}δb^{*}.
We now consider as the first step of the following considerations the situation in the reference frame that rotates with the Sun, that is, the frame whitin which the CIR and HSS are at rest, whereas the s/c passes through it. In this reference frame, the solar wind flow velocity, which is radial in the rest frame of reference, becomes aligned with the magnetic field lines because the velocity now has an additional component, corresponding to the solar rotational velocity (see Fig. 3a, and, e.g., Figs. 1, 4, and 5 in Hundhausen 1972). The advantage of the frame rotating with the Sun is that the diffusion and the wind flow are coaligned, both following the magnetic field lines, so that the diffusion is antiparallel with the convective flow. Because κ_{⊥} ≪ κ_{∥} means that the distribution of Φ along a given field line, Φ(ξ), where ξ represents the curvilinear coordinate along the field line, is independent of the distribution of Φ along neighboring field lines. Thus, considering the situation along a particular field line, Eq. (3) can be written as
that is,
where we have introduced the abbreviation for the GCR gradient along a given field line,
Equation (6) defines in firm physical terms a local gradient of Φ(ξ), G_{ξ}, along a given field line. If G_{ξ} is approximately constant over a distance range from ξ_{1} to ξ_{2} along a particular field line, or in other words, if we take the mean value of G_{ξ} over this segment of the field line, , the solution of Eq. (5) reads
For the vicinity of 1 au, we can write
where Φ_{1 au} represents the value of Φ at s/c located at 1 au, and we set ξ = 0 at the intersection of a given field line with the 1 au s/c trajectory.
This means that when we have the 1 au s/c measurements of Φ_{1au} (or analogously, at any other heliocentric distance at which a particular s/c is located), as well as the 1 au measurements of V and B (providing an estimate of G_{ξ}), Eq. (8) provides an estimate of Φ at any position [x, r] within or out of CIR and HSS. Moreover, if measurements of Φ(t) and V from two s/c located at two different heliocentric distances, r_{1} and r_{2}, were available, we could easily estimate the value of within the corresponding distance range ξ_{1} − ξ_{2} along a particular field line, assuming the nominal shape of the corresponding Parker spiral. This would also provide an estimate of the mean value within the distance range ξ_{1} − ξ_{2} by using Eq. (6). By also employing the measurements of B and δB, we can also find an empirical expression relating κ_{∥} to B or δB.
Finally, we emphasize that Eq. (8) clearly shows that after the s/c exits the HSS, in an ideal case in which the flow speed and the magnetic field attain preCIR values, that is, the value of G_{ξ} returns to the initial state, the cosmic ray particle density should also return to the preevent level. If the flow speed and the magnetic field values differ from those measured before the CIR, the postevent cosmic ray count rate level should be different from that recorded in the preevent state. If values of V and B are lower than in the preevent state, the socalled overrecovery would be observed, in which the cosmic ray count rate after the HSS is higher than prior to it. If values of V and B in the postHSS phase remain higher than in the preevent phase, the countrate will stay lower than it was before the CIR.
3.2. “Operative” approximation
Although the consideration presented above describes the basic behavior of Φ throughout the CIR and HSS structure, it is inconvenient for quantitative analysis and comparison with the measured values of CR at 1 au because the values of Φ_{1 au} at different 1 au s/c positions that are basically governed by Eq. (5) cannot be related directly to the preevent values Φ_{0}. This is primarily because in the absence of measurements by two (or more) s/c overtaken by the same CIR structure at different heliocentric distances, the value of Φ at a certain referent distance is not available (for uncertainties of Φ at large heliocentric distances, see, e.g., Fig. 3 in McKibben et al. 1982).
In the following, we therefore approach the problem in an alternative way that provides an approximative, but directly operative quantitative analysis of Φ_{1 au}(t). More precisely, the proposed procedure described below enables estimating the ratio of the value of Φ at 1 au along the xaxis (for the definition of the coordinate system we used, see Fig. 3; in Fig. 3b we show the situation in the fixed reference frame) relative to the preevent value Φ_{0}, that is, the values of Φ_{1 au}(x)/Φ_{0}, which can then be readily converted into the ratio Φ_{1 au}(t)/Φ_{0}, as measured by s/c during the CIR and HSS passage over it. To do this, we converted the dependence Φ_{x}(ξ) into an equivalent that defines the dependence of Φ on the radial coordinate, r, at a given xcoordinate, Φ_{x}(r). In other words, the gradient G_{ξ}(x) as estimated for 1 au was converted into the radial gradient, .
The procedure of obtaining the ratio Φ_{1 au}(x)/Φ_{0} was started by calculating the values of G_{ξ} along the xaxis G_{ξ}(x, 1 au). Then, taking into account the nominal orientation of the local Parker spiral (geometrically, the Archimedian spiral), we calculated the component of the gradient G_{ξ}(x, 1 au) in rdirection to obtain local values, G_{r}(x, 1 au). The Archimedian spiral is defined in polar coordinates [r, θ] as
where v_{r} is a constant radial speed and ω is a constant rotational speed. The winding angle, representing the tilt angle of the field lines relative to the radial direction, reads
The winding angle relates G_{ξ} and G_{r} as
This situation is sketched in Fig. 3b (fixed reference frame) for three different s/c positions, x_{i}, denoted as x_{1}, x_{2}, and x_{3}, where the corresponding values G_{ξ1}, G_{ξ2}, and G_{ξ3} are also sketched together with the related components G_{r1}, G_{r2}, and G_{r3}. The CIR and HSS structure is shaded in gray, and the frontal edge of the CIR is marked by the bold red line. The magnetic field lines are sketched by the thin blue lines. The heliocentric distance of the CIR boundary at coordinates x_{1}, x_{2}, and x_{3} are denoted as r_{1}, r_{2}, and r_{3}, respectively. Ahead of the CIR boundary lies a quietstate slow solar wind (shaded light gray), characterized by the GCR radial gradient G_{r0} and the fieldaligned gradient G_{ξ0}. At the bottom of Fig. 3b, the typical behavior of the 1 au solar wind flow speed, V, the magnetic field magnitude, B, and the cosmic rate count rate, CR, are sketched schematically.
In the next step, we estimated for each s/c position x within CIR and HSS the values of G_{r}(x, r) along the rdirection within the CIR and HSS structure, that is, from r = 1 au up to the outer CIR boundary located at r = r_{cir}(x), meaning that we converted the measurements obtained by s/c at 1 au into what is expected in the entire CIR and HSS structure within the range from its 1 au onset to its 1 au end. To illustrate this mapping procedure, we considered a moment t = t_{i} after t = 0, set at the CIR onset as measured by the s/c. At this moment, the plasma element that was at 1 au at t = 0 now already lies at r_{cir} = 1 au + V_{r0}t_{i}, where V_{r0} is the solar wind velocity measured at t = 0, so that r_{cir} represents the radial distance of the frontal CIR boundary at t = t_{i}. Analogously, we can track the trajectory of any plasma element that was at 1 au in the period from t = 0 to t = t_{i}. For example, the element that was at 1 au at t_{j} < t_{i}, showing the radial speed V_{j}, at t_{i} will be at the distance r_{j}(t_{i}) = 1 au + V_{j} (t_{i} − t_{j}). When we take into account that the solar wind speed around 1 au does not change with distance, and when we also assume that in the considered distance range the magnetic field does not change much, the physical conditions within a given plasma element remain mainly unchanged over its trajectory within the t_{i} − t_{j} interval. Thus, the values of G_{r} estimated from s/c measurements (as previously described) for any t_{j} in the period from t = 0 to t_{i} can be attributed to the corresponding distance r_{j} between r_{cir}(x_{i}) and 1 au. Based on this mapping of G_{r}(r) at each x_{i}, it is straightforward to estimate the mean value of G_{r} from 1 au to r_{cir} for a given s/c position x_{i}, that is, for the corresponding moment t_{i}.
As a digression, we note here that the increasing V_{r} with the increasing time t_{i}, that is, increasing distance x_{i} within the CIR and HSS structure, also means that faster plasma elements catch up the preceding slower elements. Thus, over a larger heliospheric distance range, the previous approximations cannot be applied because the plasma and magnetic field within a plasma element are compressed due to [∂V_{r}/∂r]_{x} < 0. However, the compression effect is partly compensated for by the 1/r^{2} lateral expansion of the structure.
The steepening of the frontal CIR structure leads to one more important effect, and that is the forwardshock formation (Gosling 1996; Gosling & Pizzo 1999). As the spatial V_{r}(r) profile causes a gradual steepening of the frontal part of the CIR, it will at a certain distance r_{s} inevitably transform into a shock. Denoting the radial speed of the frontal CIR element as V_{r0}, the maximum 1 au HSS radial speed as , their difference as , the time elapsed from the 1 au CIR onset to the moment at which the HSS attains as Δt_{max}, and the corresponding distance between the CIR outer boundary and 1 au as Δr, we can estimate the time needed for the shock formation as Δt = Δr/ΔV = V_{r0}Δt_{max}/ΔV.
We take as typical examples the measurements presented in Figs. 1a and e, where the frontal radial speed was V_{r0} ≈ 350 km s^{−1} and the maximum speeds within the HSS were and 800 km s^{−1}, achieved within ≈1 and 3 d, respectively. The corresponding times needed for the full completion of the shock after the frontal CIR edge passed over the 1 au s/c are then equal to 1.4 and 2.3 d, respectively, and the corresponding heliocentric distances are ≈1.5 and 2.1 au. This estimate is broadly consistent with observations as well as with numerical results (e.g., Gosling 1996; Gosling & Pizzo 1999).
We return to the previous considerations of the radial GCR gradient. The obtained values of G_{r}(x_{i}, 1 au), whose behavior is schematically sketched in Fig. 3c for three different s/c positions, x_{1}, x_{2}, and x_{3}, provide the estimate of the cosmic ray count rate variation along r for a given s/c position x_{i}. Beyond the outer CIR boundary, located at r_{cir}(x_{i}) the value of G_{r}(x_{i}, r_{cir}(x_{i})) is taken to be constant (e.g., McDonald et al. 1981), having the value G_{0r} (Fig. 3b), which can be chosen from some of the values reported in a number of previous studies, for instance, 4.8%/au, as reported by Gleeson & Axford (1968); 3.5%/au, as estimated by McDonald et al. (1981); 1−10%/au, and 3−7%/au over different phases of the solar cycle, as reported by Fisk et al. (1998) and Fujii & McDonald (1993), respectively; 2.7%/au, as found by de Simone et al. (2011) and Gieseler & Heber (2016); 2−6.6%/au, as determined for different energy ranges by Marquardt & Heber (2019). These values depend not only to the phase of the solar cycle and the heliospheric distances considered, but also on the proton energies taken into account, consistent with the predictions by Parker (1965).
This enables an estimate of the values of Φ(x, 1 au), starting from the quietstate value of Φ(x, r_{cir}) at the outer boundary of CIR and HSS as defined by the quietstate radial gradient G_{0r}. The behavior of Φ_{xi}(r) at a particular xcoordinate, x_{i}, is now described by the equation for the rdirection,
Equation (12) defines the behavior not only at 1 au, but rather, at any heliocentric distance r within, as well as beyond 1 au, as given by the solution of Eq. (12),
To obtain the values of Φ at 1 au for a given x_{i} coordinate, Eq. (13) has to be integrated step by step, starting from r_{cir}(x_{i}), where the cosmic ray count rate has the nominal quietstate value for that heliospheric distance, Φ_{0}(r_{cir}(x_{i})), down to r = 1 au to derive the needed value Φ_{1au}(x_{i}).
The heliocentric distance
is calculated by substituting into Eq. (14) the value Δr_{cir}(x_{i}) = V_{0r}Δt_{i}, where V_{0r} denotes the preevent radial solar wind speed, and Δt_{i} represents the time measured from the onset of CIR at 1 au, t_{0}, that is, Δt_{i} = t_{i} − t_{0}.
The obtained values of G_{r}(x, r) within the CIR and HSS structure define the behavior of G_{r}(r) for any coordinate x_{i}. This is schematically sketched in Fig. 3c for three different xcoordinates, x_{1}, x_{2}, and x_{3}, corresponding to segments at the beginning portions of the CIR, the place where G_{r}(1 au) is maximum, and a portion within the rear parts of HSS, respectively. Thus, Eq. (13) can now be evaluated to obtain Φ_{xi}(r), and thus Φ_{xi}(1 au), and consequently, Φ_{1 au}(t_{i}). This is schematically depicted in Figs. 3d and e, respectively.
On the one hand, bearing in mind all the approximations that we used, and on the other hand, the inherent noise in measurements and the process of smoothing, we simplify below the calculation of Φ_{xi}(1 au), i.e., Φ_{1 au}(t_{i}), where instead of integrating Eq. (13) step by step along the rdirection at a given xcoordinate, x_{i}, we estimate the mean value of G_{r}, along r at x_{i}, to get , and use it to obtain the value of Φ_{1 au}(x_{i}) as
where Δr(x_{i}) = r_{cir}(x_{i})−1au represents the radial extent of CIR and HSS beyond 1 au at a given s/c position x_{i}, as defined below Eq. (14).
Having the values of Φ_{1 au}(x_{i}), i.e., Φ_{1 au}(t_{i}), the values of Φ_{1 au}(t_{i}) can be straightforwardly converted into the change of Φ_{1 au} relative to the preCIR value Φ_{0}, expressed in percentiles,
which can be directly compared with the observations in the form as presented in Figs. 1d and h, or the statisticsbased generic form shown in Figs. 2a and b.
Finally, we note that the described model is to a certain degree similar to the propagating diffusive barrier approach (PDB, Wibberenz et al. 1998; Wibberenz & Cane 2000), which describes the effects of merged interaction regions (MIRs) at large heliocentric distances. However, we note that in PDB, which are valid for large heliocentric distances at which the magnetic field lines are almost concentric circles, the radially outward convection is balanced by the diffusion perpendicular to the magnetic field.
4. Comparison with observations
In the following, the model proposed in the previous section is applied to the observations, using the results of the study presented in Paper I, where the characteristics of a longliving CIR that was recorded over 27 solar rotations in the period from June 2007 to May 2009, are analyzed in detail. First, we consider in Sect. 4.1 the generic scheme shown in Fig. 2a, as well as the results of the SEA displayed in Fig. 2b. Then, in Sect. 4.2 we consider two specific examples, recorded from DOY 183 to 191, 2007 (rot2) and from DOY 58 to 68, 2008 (rot11). They are presented in Fig. 1.
As specified in Sect. 3, the main model input parameters are the solar wind radial flow speed, V_{r} and the magnetic field strength, B, or its fluctuations, δB. We also note that we employed the previously mentioned free parameter k when we relate the magnetic field strength (or its fluctuations) to the radial gradient of the GCR count rate, G_{r}, because the dependence of the diffusion coefficient on the magnetic field strength is not well known and the background heliospheric GCR radial gradient, G_{r0}, as inferred in a number of previous studies, spans quite a broad range of values. The GCR gradient should be higher in stronger fields, that is, roughly G_{r} ∝ B due to κ_{∥} ∝ 1/B (Potgieter 1998; Wibberenz et al. 1998; Wibberenz & Cane 2000). We emphasize that the free parameter k does not affect the overall behavior of the modeled cosmic ray count rate within CIR and HSS. It is only used to adjust the overall amplitude of the model outcome to the observed amplitudes of the cosmic ray count rate variations. In other words, when the value of k is determined from fitting the model results to the observations, we obtain a better insight into the relation between the diffusion coefficient and the magnetic field strength, or its fluctuations.
4.1. Generic SEA and “schematic” time profiles
In the left column of Fig. 4 we show a comparison of the outcome of SEA (see Fig. 2b), with the model result for scaling factors k_{B} = 0.8 and k_{δB} = 0.85, estimated by matching the model with the observations. The values of coefficients k were obtained by minimizing the rms differences between the modeled and observed data. The model output is based on a smoothed V_{r}(t) curve and smoothed curves b(t) = B(t)/B_{0} and δb(t) = δB(t)/δB_{0}. The smoothing is based on segmental linear fitting (hereafter, SLF), where linear fitting is applied to each segment between two neighboring knees in the measured time curves of a given parameter (for details, see Appendix A). We chose the SLF option instead of directly comparing the 1 h observed data with the model results simply to reduce the noise in the observedminusmodel data, especially in the cases of B and δB, which are characterized by quite high noise. In Fig. 4c the SEA SLFsmoothed V_{r}(t) curve (dashed black line) is compared with the inverted curve of the SEA SLFsmoothed cosmic ray count rate variation, CR(t) (blue lines), to illustrate a close similarity of the V_{r}(t) and −CR(t) behavior. In Fig. 4e the modeled SEAbased Φ(t) curves are directly compared to the SLFsmoothed SEA CR(t) curve. The red rectangle emphasizes the period of rising solar wind speed, ∂V_{r}/∂t > 0.
Fig. 4. Comparison of the model results and observations: SEA, based on the study in Paper I, shown in Fig. 2a (left) and the schematic generic behavior, based on Fig. 2b (right). Panel a: magnetic field, B (gray), and its fluctuations δB (green), normalized with respect to the preevent values, b = B/B_{0} and δb = δB/δB_{0}, multiplied by the free parameter k (written in the inset), based on SEA. Panel b: magnetic field, B, normalized with respect to the preevent values, b = B/B_{0}, multiplied by the free parameter k (written in the inset), based on the schematic behavior. Panels c and d: solarwind flow speed, V_{r} (dashed black line) and the inverted graph of the cosmic ray count rate, CR (blue line). Panels e and f: observed CR (blue line) compared to the model results (dashed gray line) based on the behavior of b; in panel e: model results for the δb option are presented as the dashed green line. The red rectangles mark the period of the rising solar wind speed. 
In Figs. 4b, d, and f, a comparison of the generic schematic time profiles (Fig. 2a) and the model B option results for k = 1.5 is presented, analogously to that concerning the generic SEAoption shown in Figs. 4a, c and e, respectively. In both options, model results and measured data match well. Moreover, the nose in the CR time profile corresponds to a slight preCIR depletion in the velocity curve.
4.2. Two specific examples: rot2 and rot11
In Fig. 5 we present a comparison of the model results and observations for rot2 (left column) and rot11 (right column). These two examples are chosen due to their welldefined and “clean” CIR and HSS structure, that is, there was no effect of other transient disturbances, such as ICMEs or another CIR that overlapped those we studied. Furthermore, they are different to a certain degree, not only by the amplitudes and perturbation duration, but also because rot2 shows the overrecovery of CR at the end of HSS (Fig. 5e, see also Fig. 1d), whereas rot11 shows a distinct secondary peak in B and δB at the rear of HSS with a corresponding secondary decrease of CR (Fig. 5f, see also Fig. 1h).
Fig. 5. Comparison of the model results and observations for two CIR events. Left: rot2 (DOY = 184 and 191, 2007). Right: rot11 (DOY = 58–68, 2008). Panels a and b: magnetic field and its fluctuations normalized with respect to the preevent values (b and δb, gray and green; respectively). Panels c and d: solarwind flow speed (dashed black) and the inverted CR(t) graph (blue). Panels e and f: observed CR(t) values (blue) compared to the model results based on the behavior of b (dashed gray) and δb (dashed green). 
Figures 5a and b show the SLFsmoothed behavior of the normalized magnetic field strength, b = B/B_{0}, and its fluctuations, δb = δB/δB_{0}, for rot2 and rot11, respectively. The displayed values are scaled by the factor k we used in the calculations (k_{b} = 0.65 and 0.9 for b, and k_{δb} = 0.3 and 0.75 for δb, respectively). In Figs. 5c and d, the SLFsmoothed behavior of V_{r} is drawn together with the SLFsmoothed inverted CR curve to emphasize the similarity of the two curves. In Figs. 5e and f, the observed behavior of CR (blue) is compared with the model results based on b(t) (dashed gray) and δb(t) (dashed green). The graphs are very similar to the modeled and observed behavior, including the overrecovery in the case of rot2, as well as the secondary dip in rot11.
5. Analysis of the model results
In Fig. 6 the model results (yaxes) are plotted versus the observed results (xaxes) for the rot2, rot11, SEA, and generic cases. The slopes of the linear leastsquares fits, y = c_{1}x + c_{0}, are found to be c_{1} ≈ 1 on average, and the corresponding yaxis intercept is c_{0} ≈ 0 (a perfect match corresponds to c_{1} = 1, c_{0} = 0). The correlation coefficients are high; they range from cc = 0.78 to 0.94. The results displayed in Fig. 6 show that the correlation coefficients are better for the boption than for the δboption, which is expected because of the various effects that have to be taken into account (high noise in the δbdata, a strong dependence on the time resolution, the unknown spatial scales and timescales of the fluctuations that affect CR, etc.). However, we note that the clear correlation between the model results and the observations for the δboption is consistent with the hypothesis that the magnetic filed fluctuations govern the GCR diffusion term of the GCR transport equation, as expected from the theoretical point of view (Parker 1965). We note that δB ∝ B (see Sect. 6), which basically indicates that the fluctuations are of comparable normalized values, meaning that the turbulence attains larger amplitudes in stronger fields. The relative values remain comparable, however, δB/B ≈ 0.1.
Fig. 6. Quantitative comparison of the model values (yaxes) with the observed values (xaxes). Panels a–c: data drawn in blue and red represent the model results based on the normalized magnetic field strength b and its normalized fluctuations δb, respectively. Panel d: outcome for generic schematic profile based on the b model option is displayed. In the insets, linear leastsquares fits are written, together with the corresponding correlation coefficients, cc. 
In Fig. 5e the CR overrecovery is evident in both the observed and modeled values for rot2, which is also reflected in Fig. 6a (dots in the top right quadrant; x > 0, y > 0, respectively). Moreover, the model reproduces the secondary dip of CR in the case of rot11 (Fig. 5f) very well.
In Fig. 7 we present the distributions of the observedminuscalculated CR values, O–C, for rot2, rot11, and SEA for model results based on the b and δboption. In all cases, the O–C values of the CR instantaneous amplitudes (expressed in percent; see Fig. 5) peak in the range from −0.5 to +0.5%. The best results are obtained for the rot2/boption and SEA/boption (Figs. 7b and e, respectively).
Fig. 7. Distributions of observedminuscalculated (O − C) values: Panel a: rot2 modeled using the normalized magnetic field values b. Panel b: rot2 modeled using the normalized magnetic field fluctuations δb. Panels c and d: the same as in panels a and b, but for rot11. Panels e and f: the same as in panels a and b, but now for the SEA results. 
In Table 1 we summarize the basic results of our comparison of the model results and observations. In Col. 1 all the considered options are specified. In Cols. 2 and 3 the average and median of the observedminuscalculated values are shown. In Cols. 4–6 the basic parameters of the linear leastsquares fit are presented (the slope, c_{1}, the yaxis intercept, c_{0}, and the correlation coefficient, cc, respectively). In the bottom row, the mean values for each column are displayed. The average and median values show a good overall match between the calculated and the observed values (the mean values presented in the last row are 0.05 ± 0.20% and 0.10 ± 0.22%, respectively). Generally, the slopes of the calculatedversusobserved linear leastsquares correlations are grouped around c_{1} = 1, with the mean value 0.99 ± 0.16. The yaxis intercept, c_{0}, ranges from −0.44 to +0.46, with the mean value of −0.04 ± 0.36, that is, c_{0} ≈ 0. The correlation coefficients range from cc = 0.78 to 0.94, with a mean value of 0.85 ± 0.06.
Analysis of the results.
6. Discussion
First of all, we emphasize that our results for the overall generic data, as well as for the two specific examples that were analyzed in more detail, are related to only one longlasting CIR. Thus, the results we obtained should not be a priori generalized to all CIR events. This can be seen already from the fact that the scaling factor k is quite different in the case of rot2 and rot11, as well as in comparison with the SEA and schematic generic profiles. Still, the value of k for the boptions is approximately the value of k ≈ 1 on average (from the four inferred values of k, see Table 1, the mean value is ). The values of k for the δboption seem to be somewhat lower (k_{δb} = 0.30, 0.75, and 0.85 for rot2, rot11, and SEA, respectively), with a mean value of , but they are still within the same order of magnitude.
The dependence δB(B) for rot2 and rot11, presented in Fig. 8, shows that δB and B are correlated (the correlation coefficients are cc = 0.59 and 0.78, respectively). Although the number of data points is relatively small (n = 15 and 27, respectively), the Ftest confidence is high (P = 98.5% and P ≫ 99%, respectively). When the data for both data samples are put together, the correlation coefficient adds up to cc = 0.70, with an Ftest confidence P ≫ 99%.
Fig. 8. Dependence of the level of the magnetic field fluctuations, δB, on the related magnetic field strength, B, for rot2 (red) and rot11 (blue), together with the corresponding linear leastsquares fit. The fit embracing all data points is shown by the dashed black line. 
The linear leastsquares fit for rot2 results in δB = (0.06 ± 0.02) B + (0.13 ± 0.12), for rot11 in δB = (0.13 ± 0.02) B − (0.21 ± 0.11), and for the joined data, the relation becomes δB = (0.098 ± 0.015) B − (0.06 ± 0.08). The slopes of c_{1} ≈ 0.1 and the yaxis intercepts of c_{0} ≈ 0 confirm the close correlation of these two magnetic field parameters, where the field fluctuations δB achieve about 10% of the instantaneous magnetic field strength B on average. The tight correlation of δB and B explains why the model boption and δboption both work well and why the scaling factor k in both options has a quite similar value.
At this point, it is appropriate to discuss the effect of the magnetic field fluctuations δB on the GCR particle propagation. The OMNIreported fluctuations δB represent the root mean square of the average magnetic field magnitude within the time interval defined by the time resolution^{2}, in this case, 1 h. The measurements show that during the increasedB phase, the level of the fluctuations most often varyies from 0.2 to 2 nT, for example. The typical duration of these episodes of increased fluctuations is δt = 1 − 5 h. The lengthscale equivalent, δx, can be calculated as δx = v_{rot}δt, where v_{rot} represents the velocity of the CIR rotation at 1 au. In this way, we obtain δx ≈ 0.01 − 0.05 au, that is, δx ≈ 1.5 − 7 × 10^{6} km, and this is somewhat larger than a typical 1 GeV proton gyroradius in the 5−10 nT environment. This means that these fluctuation clumps taken as an entity should not affect the propagation of these particles too much (Parker 1965), and thus should not affect the diffusion coefficient significantly. The considered fluctuation clumps are a result of many fluctuations on still smaller spatial scales that would also affect the propagation of 1 GeV protons, that is, the related diffusion coefficient and thus the GCR count rate, and this would need an additional, more careful comparison of the CR behavior with δB on different timescales.
Finally, we comment on the preCIR density increase and the associated V depletion (Figs. 1a and c, respectively; see also Fig. 2b). It can be speculated that a possible explanation might be related to the plasma pressure gradient within the frontal part of the CIR. This pressure gradient should cause a component of the flow speed along the field lines out of the frontal boundary of CIR that is compressed and thus carries additional plasma ahead of the CIR, that is, causes a weak preCIR density increase. As this flow component has a radial component that opposes the overall preCIR solar wind flow, this effect should also lead to the formation of the depletion in the radial speed time profile. In the situation shown in Fig. 1e, where there is no significant preCIR V depletion, there is no preCIR density increase either.
This scenario might be also related to the formation of the preCIR nose in the CR time profile. Namely, Fig. 2b shows that the nose occurs during the period of V depletion (see also Fig. 1d). According to the previously proposed model, a lower V, that is, a reduced convective effect, should lead to an increase in CR. The situation presented in Figs. 1e–h is somewhat different because the nose is still observed in CR, although there is no preCIR V depletion. On the other hand, inspecting the B and δB graphs in Fig. 1f, we find that in the period of the noserelated CR increase, there is a depletion in both B and δB, and according to the model, this should enable a more efficient diffusion and thus the CR increase.
Because the proposed explanation for this weak increase in GCR flux before the main FD decrease is quite speculative and cannot be fully confirmed by our analysis, we should also consider another possible explanation that is related to the mirror effect at a steep gradient of the magnetic field (including the situation in which a forward shock has already formed) in the region in which the fast solar wind compresses the plasma in the interaction with slow solar wind. In this respect, we note that in the situation in which cylindrical geometry cannot be applied and the adiabatic approximation is not valid, the mirror effect is quite complex and cannot easily be included in the hydrodynamical model such as employed here, that is, the numerical approach has to be applied (see Kirin et al. 2020). The results presented by Kirin et al. (2020) show that a certain population of particles tends to stay bound to the shock for a certain period of time, indicating that this effect may explain the weak increase in GCR flux before the main FD decrease. Related to the compression in the frontal parts of CIRs, we also emphasize the socalled snowplow effect (e.g., Guo et al. 2021), which is observed as the plasma density increase that starts even before the solar wind speed starts to rise, and which may also cause the GCR preCIR increase (Guo et al. 2021).
Finally, we note that although our model is based on 2D approach, it uses a realistic CIR structure in the equatorial plane as a structural background for a quite simple analytical model for the GCR modulation, which is based on the simplified Parker transport equation. Our model is similar to the steadystate diffusionconvection model by Richardson et al. (1996), which in addition considers the adiabatic deceleration and longitudinal variations in the solar wind speed. In both models, the effect of the convection due to the solar wind outward flow is locally balanced by the inward diffusion of GCRs. Richardson et al. (1996) concluded that the main cause of the CIRrelated GCR depression is the increase in solar wind speed, that is, the enhanced outwardconvection effects, whereas the enhanced magnetic field in the frontal parts of CIR as well as the turbulencerelated drifts play a less important role. In contrast, we have found that the magnetic field and turbulence enhancement play an important role in the GCR modulation. The results of our model are on the other hand very similar to the much more sophisticated 3D numerical model by Alania et al. (2011). They considered a steadystate heliosphere in which the CIR is represented as heliolongitudinal variation of solar wind speed and applied the full Parker transport equation. They found that the variation in GCR intensity is inversely correlated with the modulation parameter, which is proportional to the product of the solar wind velocity V and the strength of the interplanetary magnetic field B. This is exactly the case in our model. Moreover, their calculations have shown that changing the HCS structure in the model does not significantly affect the results (for a detailed study of the effect of HCSs on GCRs, see, e.g., Thomas et al. 2014).
Results similar to those obtained by Alania et al. (2011) were also found by Guo & Florinski (2016) using the 3D coupled MHDCR transport model. They also found that that GCR variations depend on the ratio of the diffusion coefficients in the fast and slow solar winds, as is the case in our model. Finally, we compared our results to the most recent dedicated 3D coupled MHDCR transport models by Kopp et al. (2017) and Luo et al. (2020). Their simulations show that enhanced convection, reduced diffusion, and drifts due the magnetic field enhancement may contribute to the recurrent FD. However, the simulations cannot distinguish the relative contribution of different mechanisms. Our results indicate that convection and diffusion are enough to describe the basic properties of recurrent FDs. On the other hand, we do not exclude the possibility that drifts are needed to explain the fine structures within the recurrent FD. We also note that in all of the models above, including our model, the transport parameters need to be tuned to successfully reproduce the observations. In the quantitative aspect, our simple analytical model therefore performs equally well as the much more sophisticated and physically more realistic models.
7. Conclusion
Our analysis explains the physical background of RFDs caused by CIR and HSS structures in terms of the convectiondiffusion approach. It is demonstrated that the enhanced convectioneffect caused by the increased velocity of HSS and the reduced diffusioneffect caused by the enhanced magnetic field fluctuations within the CIR and HSS structure are sufficient to explain the observations. This approach was verified by applying a relatively simple analytical model to the observations.
The main guiding observational fact that leads in this direction is the close similarity of the V(t) curve and the inverted CR(t) curve (Figs. 4 and 5; for more examples, see Paper I). However, the enhanced convection effect cannot explain the CR(t) behavior without an additional ingredient. After we inspected the in situ solar wind data, we realized that this has to be the increase in magnetic field strength B within the CIR structure, or a physically more appropriate option, the enhanced magnetic field fluctuations δB that reduce the CR diffusion.
The combined effect of the increased solar wind speed and the enhanced magnetic field fluctuations practically entirely explain the overall RFD behavior in the 1 au range, related to the CIR and HSS structure at these distances. The two effects taken together can also explain the phenomenon of the overrecovery that is sometimes observed after the RFD events. It is likely that the overrecovery occurs (based on the two employed examples) in the events in which the solar wind speed and the magnetic field fluctuations are both lower than in the preevent state.
Acknowledgments
This work has been supported by Croatian Science Foundation under the project 7549 “Millimeter and submillimeter observations of the solar chromosphere with ALMA”. B.V. and M.D. also acknowledge support by the Croatian Science Foundation under the project IP2020029893 (ICOHOSS). B.H. acknowledges the financial support via projects HE 3279/151 funded by the Deutsche Forschungsgemeinschaft (DFG).
References
 Alania, M. V., Modzelewska, R., & Wawrzynczak, A. 2011, Sol. Phys., 270, 629 [Google Scholar]
 Badruddin, B., & Kumar, A. 2016, Sol. Phys., 291, 559 [NASA ADS] [CrossRef] [Google Scholar]
 Čalogović, J., Vršnak, B., Temmer, M., & Veronig, A. M. 2009, in Universal Heliophysical Processes, eds. N. Gopalswamy, & D. F. Webb, IAU Symp., 257, 425 [Google Scholar]
 Cane, H. V. 2000, Space Sci. Rev., 93, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Chhiber, R., Subedi, P., Usmanov, A. V., et al. 2017, ApJS, 230, 21 [Google Scholar]
 Chih, P. P., & Lee, M. A. 1986, J. Geophys. Res., 91, 2903 [NASA ADS] [CrossRef] [Google Scholar]
 Cranmer, S. R. 2009, Liv. Rev. Sol. Phys., 6, 3 [Google Scholar]
 de Simone, N., di Felice, V., Gieseler, J., et al. 2011, Astrophys. Space Sci. Trans., 7, 425 [NASA ADS] [CrossRef] [Google Scholar]
 Dumbović, M., Vršnak, B., Čalogović, J., & Karlica, M. 2011, A&A, 531, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dumbović, M., Vršnak, B., Čalogović, J., & Župan, R. 2012, A&A, 538, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dumbović, M., Vršnak, B., Guo, J., et al. 2020, Sol. Phys., 295, 104 [CrossRef] [Google Scholar]
 Dumbović, M., Vršnak, B., Temmer, M., Heber, B., & Kühl, P. 2022, A&A, 658, A187 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fisk, L. A., Wenzel, K. P., Balogh, A., et al. 1998, Space Sci. Rev., 83, 179 [NASA ADS] [CrossRef] [Google Scholar]
 Forbush, S. E. 1937, Phys. Rev., 51, 1108 [NASA ADS] [CrossRef] [Google Scholar]
 Fujii, F., & McDonald, F. B. 1993, Int. Cosm. Ray Conf., 3, 477 [NASA ADS] [Google Scholar]
 Ghanbari, K., Florinski, V., Guo, X., Hu, Q., & Leske, R. 2019, ApJ, 882, 54 [NASA ADS] [CrossRef] [Google Scholar]
 Giacalone, J. 1998, Space Sci. Rev., 83, 351 [NASA ADS] [CrossRef] [Google Scholar]
 Gieseler, J., & Heber, B. 2016, A&A, 589, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gleeson, L. J., & Axford, W. I. 1968, Can. J. Phys. Suppl., 46, 937 [Google Scholar]
 Gosling, J. T. 1996, ARA&A, 34, 35 [NASA ADS] [CrossRef] [Google Scholar]
 Gosling, J. T., & Pizzo, V. J. 1999, Space Sci. Rev., 89, 21 [NASA ADS] [CrossRef] [Google Scholar]
 Guo, X., & Florinski, V. 2016, ApJ, 826, 65 [NASA ADS] [CrossRef] [Google Scholar]
 Guo, X., Florinski, V., Wang, C., & Ghanbari, K. 2021, ApJ, 910, 99 [NASA ADS] [CrossRef] [Google Scholar]
 Hofmeister, S. J., Veronig, A., Temmer, M., et al. 2018, J. Geophys. Res.: Space Phys., 123, 1738 [Google Scholar]
 Hundhausen, A. J. 1972, in Interplanetary Shock Waves and the Structure of Solar Wind Disturbances, eds. C. P. Sonett, P. J. Coleman, & J. M. Wilcox, 308, 393 [NASA ADS] [Google Scholar]
 Hundhausen, A. J. 1973, J. Geophys. Res., 78, 7996 [CrossRef] [Google Scholar]
 Iucci, N., Parisi, M., Storini, M., & Villoresi, G. 1979, Nuovo Cimento C Geophys. Space Phys. C, 2, 421 [NASA ADS] [CrossRef] [Google Scholar]
 Jämsén, T., Usoskin, I. G., Räihä, T., Sarkamo, J., & Kovaltsov, G. A. 2007, Adv. Space Res., 40, 342 [CrossRef] [Google Scholar]
 Jian, L., Russell, C. T., Luhmann, J. G., & Skoug, R. M. 2006, Sol. Phys., 239, 337 [NASA ADS] [CrossRef] [Google Scholar]
 Jokipii, J. R., Kota, J., & Merenyi, E. 1993, ApJ, 405, 782 [NASA ADS] [CrossRef] [Google Scholar]
 King, J. H., & Papitashvili, N. E. 2005, J. Geophys. Res.: Space Phys., 110, A02104 [CrossRef] [Google Scholar]
 Kirin, A., Vršnak, B., Dumbović, M., & Heber, B. 2020, Sol. Phys., 295, 28 [NASA ADS] [CrossRef] [Google Scholar]
 Kopp, A., Wiengarten, T., Fichtner, H., et al. 2017, ApJ, 837, 37 [NASA ADS] [CrossRef] [Google Scholar]
 Kota, J., & Jokipii, J. R. 1991, Geophys. Res. Lett., 18, 1797 [NASA ADS] [CrossRef] [Google Scholar]
 Kühl, P., Dresing, N., Dunzlaff, P., et al. 2013, Cent. Eur. Astrophys. Bull., 37, 643 [Google Scholar]
 Kumar, A., & Badruddin, B. 2014, Sol. Phys., 289, 4267 [NASA ADS] [CrossRef] [Google Scholar]
 Kunow, H., Dröge, W., Heber, B., et al. 1995, Space Sci. Rev., 72, 397 [CrossRef] [Google Scholar]
 Lockwood, J. A. 1971, Space Sci. Rev., 12, 658 [NASA ADS] [CrossRef] [Google Scholar]
 Luo, X., Zhang, M., Feng, X., et al. 2020, ApJ, 899, 90 [Google Scholar]
 Marquardt, J., & Heber, B. 2019, A&A, 625, A153 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 McDonald, F. B., Trainor, J. H., Lal, N., Van Hollebeke, M. A. I., & Webber, W. R. 1981, ApJ, 249, L71 [NASA ADS] [CrossRef] [Google Scholar]
 McKibben, R. B., Pyle, K. R., & Simpson, J. A. 1982, ApJ, 254, L23 [NASA ADS] [CrossRef] [Google Scholar]
 Melkumyan, A. A., Belov, A. V., Abunina, M. A., et al. 2019, Adv. Space Res., 63, 1100 [NASA ADS] [CrossRef] [Google Scholar]
 Michałek, G., & Ostrowski, M. 2001, Sol. Phys., 200, 177 [CrossRef] [Google Scholar]
 Modzelewska, R., Bazilevskaya, G. A., Boezio, M., et al. 2020, ApJ, 904, 3 [NASA ADS] [CrossRef] [Google Scholar]
 MoralesOlivares, O. G., & CaballeroLopez, R. A. 2010, Adv. Space Res., 46, 1313 [NASA ADS] [CrossRef] [Google Scholar]
 MüllerMellin, R., Kunow, H., Fleißner, V., et al. 1995, Sol. Phys., 162, 483 [Google Scholar]
 Nolte, J. T., Krieger, A. S., Timothy, A. F., et al. 1976, Sol. Phys., 46, 303 [NASA ADS] [CrossRef] [Google Scholar]
 Parker, E. N. 1965, Planet. Space Sci., 13, 9 [Google Scholar]
 Perko, J. S., & Fisk, L. A. 1983, J. Geophys. Res., 88, 9033 [NASA ADS] [CrossRef] [Google Scholar]
 Potgieter, M. S. 1998, Space Sci. Rev., 83, 147 [CrossRef] [Google Scholar]
 Potgieter, M. S. 2013, Liv. Rev. Sol. Phys., 10, 3 [Google Scholar]
 Qin, G. 2002, PhD Thesis, University of Delaware, USA [Google Scholar]
 Richardson, I. G. 2004, Space Sci. Rev., 111, 267 [Google Scholar]
 Richardson, I. G. 2018, Liv. Rev. Sol. Phys., 15, 1 [Google Scholar]
 Richardson, I. G., Wibberenz, G., & Cane, H. V. 1996, J. Geophys. Res., 101, 13483 [CrossRef] [Google Scholar]
 Rotter, T., Veronig, A. M., Temmer, M., & Vršnak, B. 2012, Sol. Phys., 281, 793 [CrossRef] [Google Scholar]
 Schwenn, R. 2006, Liv. Rev. Sol. Phys., 3, 2 [Google Scholar]
 Schwenn, R., & Marsch, E. 1991, Phys. Chem. Space, 20, 21 [Google Scholar]
 Shalchi, A., & Schlickeiser, R. 2004, ApJ, 604, 861 [NASA ADS] [CrossRef] [Google Scholar]
 Shen, Z., Qin, G., Zuo, P., Wei, F., & Xu, X. 2020, ApJ, 900, 143 [NASA ADS] [CrossRef] [Google Scholar]
 Thomas, S. R., Owens, M. J., Lockwood, M., & Scott, C. J. 2014, Sol. Phys., 289, 2653 [NASA ADS] [CrossRef] [Google Scholar]
 Tsurutani, B. T., McPherron, R. L., Gonzalez, W. D., et al. 2006, J. Geophys. Res., 111, 1 [Google Scholar]
 Verbanac, G., Vršnak, B., Veronig, A., & Temmer, M. 2011a, A&A, 526, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Verbanac, G., Vršnak, B., Živković, S., et al. 2011b, A&A, 533, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Verbanac, G., Živković, S., Vršnak, B., Bandić, M., & Hojsak, T. 2013, A&A, 558, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Vršnak, B., Temmer, M., & Veronig, A. 2007, Sol. Phys., 240, 315 [CrossRef] [Google Scholar]
 Vršnak, B., Dumbović, M., Čalogović, J., Verbanac, G., & Poljančić Beljan, I. 2017, Sol. Phys., 292, 140 [CrossRef] [Google Scholar]
 Wibberenz, G., & Cane, H. V. 2000, J. Geophys. Res., 105, 18,315 [Google Scholar]
 Wibberenz, G., Le Roux, J. A., Potgieter, M. S., & Bieber, J. W. 1998, Space Sci. Rev., 83, 309 [NASA ADS] [CrossRef] [Google Scholar]
 Zank, G. P., Matthaeus, W. H., Bieber, J. W., & Moraal, H. 1998, J. Geophys. Res., 103, 2085 [CrossRef] [Google Scholar]
 Zhang, J., Richardson, I. G., Webb, D. F., et al. 2007, J. Geophys. Res.: Space Phys., 112, A10102 [CrossRef] [Google Scholar]
 Zimbardo, G., Perri, S., Effenberger, F., & Fichtner, H. 2017, A&A, 607, A7 [Google Scholar]
Appendix A: Segmental linear fitting
In Fig A.1 the SLF procedure is illustrated. It is applied to the EPHIN GCR count rate, solar wind velocity time profile, and interplanetary magnetic field time profile. The procedure is as follows. First, the moments in a given curve are identified, where the slope of the curve (neglecting the noise) significantly changes (nodes). Then, the parts of the curve between nodes are fit by a linear fit (red lines).
Fig. A.1. Three examples of the SLF (rot11). Panel a): EPHIN GCR count rate. Panel b): solar wind velocity time profile. Panel c): interplanetary magnetic field time profile. 
All Tables
All Figures
Fig. 1. 1 au measurements of CIRs recorded from DOY 183 to 191, 2007 (rot2; left column) and from DOY 58 to 68, 2008 (rot11; right column). Panels a and e: flow speed, V. Panels b and f: magnetic field strength, B (red), and fluctuations (root mean square), δB (blue). Panels c and g: proton number density, N (red), and temperature, T (blue). Panels d and h: EPHIN cosmic ray countrate change, CR, expressed in percent of the preevent count rate. The dashed rectangle in both panels marks the period of the flow speed increase. 

In the text 
Fig. 2. “Generic” temporal behavior of the flow speed, V (blue), magnetic field, B (red), density, N, (green), and the cosmic ray countrate change expressed in percent, CR (black), based on the study of a sample of CIRs presented in Paper I. Panel a: simplified, generic, presentation based on the mean values of the amplitudes and timing of the B, N, and CR onsets, peaks, and ends, relative to the onset of V (to see better the relative amplitude of V, it is multiplied by a factor of two). The peak times are marked by vertical dashed lines. Panel b: SEA results; the variations relative to the preCIR values are normalized to the peak values, ΔV/V_{max}, ΔB/B_{max}, ΔN/N_{max}, and ΔCR/CR_{max}. The time Δt = 0 is set at the onset of the V increase. 

In the text 
Fig. 3. Panel a: sketch of the global situation (upper panel: rest frame; lower panel: frame of reference rotating with the Sun); for details, see the main text. Panel b: schematic presentation of the geometry and CIR characteristics (fixed reference frame), as applied in the model (x_{0} and x_{e} represent the beginning and end of CIR and HSS at 1 au, respectively; for details, see the main text). Panel c: dependence of the radial CRgradient, G_{r}, for three specific s/c positions (x_{1}, x_{2}, x_{3}) as a function of heliocentric distance, r. The heliocentric distance of the CIR frontal edge for the three specific s/c positions, x_{i}, are denoted r_{1}, r_{2}, r_{3}. Panel d: corresponding radial distance dependence of the cosmic ray count rate, Φ_{xi}(r), between CIR and HSS outer boundary and s/c. The red line marks the behavior in a quietSun situation at the CIR outer boundary. Panel e: sketch of the resulting Forbush decrease profile, Φ(t), as measured by the s/c. 

In the text 
Fig. 4. Comparison of the model results and observations: SEA, based on the study in Paper I, shown in Fig. 2a (left) and the schematic generic behavior, based on Fig. 2b (right). Panel a: magnetic field, B (gray), and its fluctuations δB (green), normalized with respect to the preevent values, b = B/B_{0} and δb = δB/δB_{0}, multiplied by the free parameter k (written in the inset), based on SEA. Panel b: magnetic field, B, normalized with respect to the preevent values, b = B/B_{0}, multiplied by the free parameter k (written in the inset), based on the schematic behavior. Panels c and d: solarwind flow speed, V_{r} (dashed black line) and the inverted graph of the cosmic ray count rate, CR (blue line). Panels e and f: observed CR (blue line) compared to the model results (dashed gray line) based on the behavior of b; in panel e: model results for the δb option are presented as the dashed green line. The red rectangles mark the period of the rising solar wind speed. 

In the text 
Fig. 5. Comparison of the model results and observations for two CIR events. Left: rot2 (DOY = 184 and 191, 2007). Right: rot11 (DOY = 58–68, 2008). Panels a and b: magnetic field and its fluctuations normalized with respect to the preevent values (b and δb, gray and green; respectively). Panels c and d: solarwind flow speed (dashed black) and the inverted CR(t) graph (blue). Panels e and f: observed CR(t) values (blue) compared to the model results based on the behavior of b (dashed gray) and δb (dashed green). 

In the text 
Fig. 6. Quantitative comparison of the model values (yaxes) with the observed values (xaxes). Panels a–c: data drawn in blue and red represent the model results based on the normalized magnetic field strength b and its normalized fluctuations δb, respectively. Panel d: outcome for generic schematic profile based on the b model option is displayed. In the insets, linear leastsquares fits are written, together with the corresponding correlation coefficients, cc. 

In the text 
Fig. 7. Distributions of observedminuscalculated (O − C) values: Panel a: rot2 modeled using the normalized magnetic field values b. Panel b: rot2 modeled using the normalized magnetic field fluctuations δb. Panels c and d: the same as in panels a and b, but for rot11. Panels e and f: the same as in panels a and b, but now for the SEA results. 

In the text 
Fig. 8. Dependence of the level of the magnetic field fluctuations, δB, on the related magnetic field strength, B, for rot2 (red) and rot11 (blue), together with the corresponding linear leastsquares fit. The fit embracing all data points is shown by the dashed black line. 

In the text 
Fig. A.1. Three examples of the SLF (rot11). Panel a): EPHIN GCR count rate. Panel b): solar wind velocity time profile. Panel c): interplanetary magnetic field time profile. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.