One of the very common in situ signatures of interplanetary coronal mass ejections (ICMEs), as well as other interplanetary transients, are Forbush decreases (FDs), i.e. short-term reductions in the galactic cosmic ray (GCR) flux. A two-step FD is often regarded as a textbook example, which presumably owes its specific morphology to the fact that the measuring instrument passed through the ICME head-on, encountering first the shock front (if developed), then the sheath and finally the CME magnetic structure. The interaction of GCRs and the shock/sheath region, as well as the CME magnetic structure, occurs all the way from Sun to Earth, therefore, FDs are expected to reflect the evolutionary properties of CMEs and their sheaths. We apply modelling to different ICME regions in order to obtain a generic two-step FD profile, which qualitatively agrees with our current observation-based understanding of FDs. We next adapt the models for energy dependence to enable comparison with different GCR measurement instruments (as they measure in different particle energy ranges). We test these modelling efforts against a set of multi-spacecraft observations of the same event, using the Forbush decrease model for the expanding flux rope (ForbMod). We find a reasonable agreement of the ForbMod model for the GCR depression in the CME magnetic structure with multi-spacecraft measurements, indicating that modelled FDs reflect well the CME evolution.

Forbush decreases (FDs) are depletions in the galactic cosmic ray (GCR) count rate that last typically for about a week and can be caused by coronal mass ejections (CMEs) or corotating interacting regions (CIRs). Fast CMEs that drive shocks cause large FDs that often show a two-step decrease where the first step is attributed to the shock/sheath region, while the second step is attributed to the closed magnetic structure. Since the difference in size of shock and sheath region is significant, and since there are observed effects that can be related to shocks and not necessarily to the sheath region we expect that the physical mechanisms governing the interaction with GCRs in these two regions are different. We therefore aim to analyse interaction of GCRs with heliospheric shocks only. We approximate the shock by a structure where the magnetic field linearly changes with position within this structure. We assume protons of different energy, different pitch angle and different incoming direction. We also vary the shock parameters such as the magnetic field strength and orientation, as well as the shock thickness. The results demonstrate that protons with higher energies are less likely to be reflected. Also, thicker shocks and shocks with stronger field reflect protons more efficiently.

We present an analytical diffusion-expansion Forbush decrease (FD) model ForbMod which is based on the widely used approach of the initially empty, closed magnetic structure (i.e. flux rope) which fills up slowly with particles by perpendicular diffusion. The model is restricted to explain only the depression caused by the magnetic structure of the interplanetary coronal mass ejection (ICME). We use remote CME observations and a 3D reconstruction method (the Graduated Cylindrical Shell method) to constrain initial boundary conditions of the FD model and take into account CME evolutionary properties by incorporating flux rope expansion. Several flux rope expansion modes are considered, which can lead to different FD characteristics. In general, the model is qualitatively in agreement with observations, whereas quantitative agreement depends on the diffusion coefficient and the expansion properties (interplay of the diffusion and the expansion). A case study was performed to explain the FD observed 2014 May 30. The observed FD was fitted quite well by ForbMod for all expansion modes using only the diffusion coefficient as a free parameter, where the diffusion parameter was found to correspond to expected range of values. Our study shows that in general the model is able to explain the global properties of FD caused by FR and can thus be used to help understand the underlying physics in case studies.