By using data from multiple instruments, we investigate ionospheric/thermospheric behavior during the period from 21 to 23 June 2015, when three interplanetary shocks (IS) of different intensities arrived at Earth. The first IS was registered at 16:45 UT on 21 June and caused ~50 nT increase in the SYM-H index. The second IS arrived at 5:45 UT on 22 June and induced an enhancement of the auroral/substorm activity that led to rapid increase of thermospheric neutral mass density and ionospheric vertical total electron content at high latitudes. Several hours later, topside electron content and electron density increased at low latitudes on the nightside. The third and much larger IS arrived at 18:30 UT on 22 June and initiated a major geomagnetic storm that lasted for many hours. The storm provoked significant effects in the thermosphere and ionosphere on both dayside and nightside. In the thermosphere, the dayside neutral mass density exceeded the quiet time levels by 300-500%, with stronger effects in the summer hemisphere. In the ionosphere, both positive and negative storm effects were observed on both dayside and nightside. We compared the ionospheric observations with simulations by the coupled Sami3 is Also a Model of the Ionosphere/Rice Convection Model (SAMI3/RCM) model. We find rather good agreement between the data and the model for the first phase of the storm, when the prompt penetration electric field (PPEF) was the principal driver. At the end of the storm main phase, when the ionospheric effects were, most likely, driven by a combination of PPEF and thermospheric winds, the modeling results agree less with the observations.

The most severe geomagnetic storm in solar cycle 24 started with a sudden storm commencement (SSC) at 04:45 UT on St. Patrick's day March 17, 2015. It occurred without any significant precursor X- or M-type solar flares and appeared as a two-stage geomagnetic storm with a minimum SYM-H value of -233 nT. In the response to the storm commencement in the first activation, a short-term positive effect in the ionospheric vertical electron content (VTEC) occurred at low- and mid-latitudes on the dayside. The second phase commencing around 12:30 UT lasted longer and caused significant and complex storm-time changes around the globe with hemispherical different ionospheric storm reactions in different longitudinal ranges. Swarm-C observations of the neutral mass density variation along the orbital path as well as Langmuir probe plasma measurements of all three Swarm satellites and global TEC records during the storm interval are used for physical interpretations and modelling of the positive/negative storm scenario. At mid-latitudes, positive storm signatures were observed in the Northern Hemisphere (NH) of the European sector, whereas a large positive storm occurred in the Southern Hemisphere (SH) of the American sector. The negative storm phase was found to be strongest in the Asian sector, in particular in the NH, but developed globally on March 18 at the beginning of the recovery phase. These observations pose a challenge for the global numerical modelling of thermosphere-ionosphere storm processes as the storm, which occurred around spring equinox, obviously signify the existence of other impact factors than seasonal dependence for hemispheric asymmetries to occur. First numerical simulation trials using the Potsdam version of the Upper Atmosphere Model (UAM-P) are presented to explain these peculiar ionospheric storm processes.