A review of the SCOSTEP’s 5-year scientific program VarSITI—Variability of the Sun and Its Terrestrial Impact

The Sun is a variable active-dynamo star, emitting radiation in all wavelengths and solar-wind plasma to the interplanetary space. The Earth is immersed in this radiation and solar wind, showing various responses in geospace and atmosphere. This Sun–Earth connection variates in time scales from milli-seconds to millennia and beyond. The solar activity, which has a ~11-year periodicity, is gradually declining in recent three solar cycles, suggesting a possibility of a grand minimum in near future. VarSITI—variability of the Sun and its terrestrial impact—was the 5-year program of the scientific committee on solar-terrestrial physics (SCOSTEP) in 2014–2018, focusing on this variability of the Sun and its consequences on the Earth. This paper reviews some background of SCOSTEP and its past programs, achievements of the 5-year VarSITI program, and remaining outstanding questions after VarSITI.


Introduction
The Sun, our home star, is a variable star. Its activity is due to the action of the solar dynamo, maintaining the solar magnetic field. The solar wind-the everexpanding solar atmosphere-drags the magnetic field throughout the whole heliosphere and carries transients like coronal mass ejections and high-speed solar wind streams. The Earth, together with the other planets in the solar system, is immersed in the solar wind and interacts with the plasma and magnetic fields originating from the Sun. The Sun also emits radiation in all wavelengths, with sporadically superposed contributions by solar flares, and modulates the flux of galactic cosmic rays. All these solar agents vary on time-scales from milli-seconds to millennia and beyond. Consequently, their impact on the terrestrial system also varies. The variations of the Sun and the geoeffective solar agents, and their varying impacts on the Earth, were studied by the SCOSTEP's program VarSITI-Variability of the Sun and Its Terrestrial Impact.

What is SCOSTEP
SCOSTEP, the Scientific Committee on Solar-Terrestrial Physics, was originally established as the Inter-Union Commission on Solar-Terrestrial Physics (IUCSTP) by the XIth General Assembly of the International Council of Scientific Unions (ICSU; now the International Science Council, ISC) in January 1966. At its XIVth General Assembly in September 1972, ICSU reorganized IUCSTP as a special committee with responsibility for interdisciplinary solar-terrestrial physics programs of finite duration. The designation of SCOSTEP took effect in September 1973 after the ratification of a new Constitution by the ICSU Executive Board and General Committee. At the XVIIth ICSU General Assembly in September 1978, SCOSTEP became a Scientific Committee of ICSU charged with the long-term responsibility to promote international interdisciplinary programs in solar-terrestrial physics (https://scostep.org/). Under the International Science Council (ISC), SCOSTEP is the only organization to deal with the Sun-Earth connections.
Research in the framework of these programs had led to significant improvements in understanding the physical processes within each of the Sun-Earth domains: Sun, heliosphere, magnetosphere, thermosphere/ionosphere, middle atmosphere, lower atmosphere, and climate. Some cross-disciplinary progress was also achieved during STEP and S-RAMP in understanding the domain as a whole. (2004)(2005)(2006)(2007)(2008) In 2000, the SCOSTEP's Long Range Planning Committee was constructed with the purpose to define future research programs. This group came up with a new approach for solar-terrestrial physics that led to a new paradigm-end-to-end study of scientific "problems" within the system as a whole on both short and long time-scales. The program was named CAWSES-Climate and Weather of the Sun-Earth System. This new program was intended to enable researchers in various disciplines to work together to understand the effects of the Sun on the Earth's environment.

CAWSES
The fundamental problems for CAWSES as formulated by the Planning Committee were (Schmieder et al. 2004): -Can we link the end-to-end processes that produce geoeffective coronal mass ejections, facilitate their transfer through the heliosphere, their interaction with the magnetosphere, and the production of geomagnetic storms that affect the atmosphere?
-Can we identify evidence for long-term variations of solar luminosity related to solar activity and resultant impacts on global change, compared with other climate change mechanisms? -Can we reconcile apparent responses of the middle and lower atmosphere to solar activity, identify the physical mechanisms, in comparison with anthropogenic influences, and estimate future ozone changes? -To what extent are the magnetosphere and ionosphere-thermosphere systems modulated by solar activity on long time scales, including the solar cycle, and how do variations driven by different processes interact with dynamical and radiative forcing processes from below?
To answer these questions, four projects were defined by the CAWSES Science Steering Group: solar influences on climate, space weather: science and applications, atmospheric coupling processes, and space climatology.
1.2.2.1 Theme 1: Solar influence on climate The aim of this project was to investigate the effects of solar variability on the climate of the lower and middle atmosphere. Evidences were accumulated that solar activity indeed influences climate. Using cosmogenic radionuclides which are primarily produced by the galactic cosmic rays whose intensity is modulated by the open solar magnetic field, the strength of the interplanetary magnetic field (IMF) was reconstructed for the past 10,000 years. The close relationship between the strength of the IMF and the total solar irradiance (TSI) found during the modern space era allowed TSI in turn to be reconstructed in the past from the reconstructed IMF. TSI was compared with selected climate records, and a relationship was found between solar variability and climate (Beer and McCracken 2009). Considerable evidence was also found that solar variability on decadal timescales does influence climate; however, it was noted that the solar signal is not uniformly distributed over the globe. Some mechanisms were suggested for the amplification of the solar signal, as summarized by (Haigh 2009 and the references therein). Among then, the so called "Bottom-up" mechanism works in the tropical oceans where there is less cloud and variations of solar irradiance have more sensible effects. It includes changes in atmospheric circulation associated with anomalies in horizontal temperature gradient of the sea surface. The basis of the "Top-down" mechanism connecting the stratosphere to the underlying tropospheric climate is the spectral solar irradiance (SSI) variability which is much greater than the TSI variability. Other possible mechanisms are related to modulation of El Niño by solar activity through a change in ocean circulations induced by longitudinally asymmetric changes in sea surface temperatures, etc.
1.2.2.2 Theme 2: Space weather-science and applications "Space weather" is known as the shortterm variations of the Sun, propagation of energetic particles, and electro-magnetic emissions through interplanetary space, and their effects in the magnetosphere, ionosphere, and thermosphere that can influence the performance and reliability of space-borne and groundbased technological systems and can endanger human life or health. The goals of this project were to develop dependable, robust deterministic end-to-end models that predict conditions in geospace from a quantitative understanding of the observed phenomena.
In addition to Interball, Geotail, Cluster, SOHO, ACE missions, new spacecraft were launched during CAWS ES like Hinode, RHESSI, Solar-B, and STEREO. They allowed to study in unprecedented detail processes on the Sun responsible for solar flares, energetic particles, and geoeffective transients in the solar wind. The observed properties of CMEs and ICMEs were summarized, with an emphasis on those properties that determine the ability of CMEs in causing space weather. On the applied side, special attention was paid to the effects of solar events on satellite missions in near-Earth orbit.

Theme 3: Atmospheric coupling processes
Solar and magnetospheric inputs propagate downward through the atmosphere while tropospheric effects propagate upward to the thermosphere-ionosphere system. This project sought to gain greater understanding of atmospheric coupling processes by means of observations, theory, and modeling. Global observing campaigns were conducted to characterize tides and their influence from the troposphere to the thermosphere, gravity waves and turbulence, solar and galactic cosmic rays influence on minor constituents, and layers at the extra-tropical summer mesopause. In the framework of the new topic "Equatorial atmosphere coupling processes," low latitude ionospheric disturbances associated with geomagnetic activity and electrodynamic coupling of equatorial F region with conjugate E regions were studied. A project was dedicated to atmospheric wave interactions with the winter polar vortices.
1.2.2.4 Theme 4: Space climatology Climatology includes a description and understanding of the average properties and regular variations of the system and is concerned with the long-term effects (e.g., longer than the solar rotation period). The goal of this project was to provide the data necessary to study the climatology of the Sun-Earth system with an emphasis on both a critical review of the original observations and their interpretations. Two composites of total solar irradiance measurements (PMOD and ACRIM) were updated, and a third composite (IRMB) was additionally presented (Fröhlich 2009). New geomagnetic indices [e.g., the Inter-Hour Variability (IHV), the Inter-Diurnal Variability (IDV), and Polar Cap Potential (PCP)] were devised from which the long-term variations of the solar wind parameters can be derived (Svalgaard and Cliver 2007). Special attention was paid to the long-term changes and trends found in the upper atmosphere at heights above 50 km, i.e., in the mesosphere, thermosphere, and ionosphere (Laštovička, Akmaev and Emmert 2009).
A summarizing CAWSES symposium was held in Kyoto in 2007, and selected papers from the symposium were gathered in the monograph "Weather of the Sun-Earth System (CAWSES): Selected papers from the 2007 Kyoto Symposium" (Tsuda et al. 2009) which is freely available online at https://www.terrapub.co.jp/ onlineproceedings/ste/CAWSES2007/index.html.
In summary, the first phase of CAWSES led to a significant improvement in our understanding of the solar influence on climate changes, effects of sporadic solar events on the terrestrial system and satellite missions in near-Earth orbit, and the coupling of atmospheric layers through dynamical, electromagnetic, and photochemical processes. It was recognized, however, that there was still much that needs to be achieved, and therefore, the SCOSTEP endorsed the continuation of the CAWSES program as CAWSES II during 2009-2013.

CAWSES II (2009-2013)
CAWSES II had four task groups: solar influences on Earth's climate, geospace response to altered climate, short-term solar variability and geospace, and geospace response to variable inputs from the lower atmosphere. Additional activities were related to capacity building and E-science and informatics, including the creation of scientific databases and their effective utilization in solar-terrestrial physics research. Below is a brief summary of the overview of the achievements by the CAWS ES II working groups mostly based on Tsuda et al. (2015).
1.2.3.1 TG1: Solar influences on Earth's climate An important contribution of TG1 was the highlighting of the role of the spectral solar irradiance in the solar influences on climate. Before CAWSES II, the total solar irradiance (TSI) was the only factor considered to quantify the solar influences on climate change. In addition to the "bottom-up" mechanism relying on the influence of the TSI on the tropical sea surface temperature with resulting changes in the circulation (Meehl et al. 2008), the "top-down" mechanism was proposed highlighting the role of the spectral solar irradiance (SSI) forcing on chemical-dynamical coupling via its interactions with atmospheric ozone and connecting the stratosphere to the underlying tropospheric climate (Gray et al. 2010). A major development during the CAWSES-II period for studies of energetic particle precipitation effects on the middle and upper atmosphere was the shift of focus from the previously dominant topic of solar proton events (SPE) to investigations including energetic electron precipitation (EEP). These two topics together defined as energetic particle precipitation (EPP). In the high latitude atmosphere, the enhanced ionization from EPP leads to the production of HOx and NOx which have an important role in middle atmosphere ozone balance, thus providing a potential link to dynamics and regional climate (Rozanov et al. 2012). Significant advances were also made in the field of cosmic rays. An extensive review on the results of TG1 is given by Seppälä et al. (2014).
Solar influence on climate is now accepted as an important contribution to climate variability, particularly on regional scales. During CAWSES II, the main focus has moved not only from TSI towards understanding SSI variations and their impact, but as well shifting from the global responses to more regional responses. With better understanding of SSI, the importance of the topdown stratospheric UV mechanism has been widely accepted. Improved measurements of both TSI and SSI became available leading to more reliable solar cycle variation estimates, and a new lower value for TSI was recommended for the IPCC AR5 climate simulations (Kopp and Lean 2011).

TG2: Geospace response to altered climate
The main focus of TG2 was on the effect of anthropogenic emissions of greenhouse gases (GHG) on the geospace-the region of the atmosphere between the stratosphere and several thousand kilometers above ground where atmosphere-ionosphere-magnetosphere interactions occur. GHG cause warming in the troposphere but cooling in the upper atmosphere.
GHG are not the only driver of long-term changes and trends in the upper atmosphere and ionosphere. A regionally important driver is the secular change of the main magnetic field of the Earth. Model simulation showed that the secular variation of the Earth's magnetic field may be the dominant cause of trends in the Atlantic region ionosphere (Cnossen and Richmond 2013). In the mesosphere and lower thermosphere, important additional drivers appear to be stratospheric ozone, changes in CH 4 , and mesospheric water vapor concentration, though the modelled changes are much higher than the observed (Qian et al. 2011). The poorly known trends in atmospheric circulation and particularly in atmospheric wave activity with predominantly tropospheric sources (one of potential trend drivers) remained key open question in the upper atmosphere and ionosphere trends (Laštovička et al. 2014).
A sufficient amount of information had made it possible to create a first scenario of global change in the mesosphere, thermosphere, and ionosphere (Laštovička et al. 2006). During CAWSES-II period, this scenario had been continuously improved and updated (e.g., Laštovička et al. 2012;Danilov 2012;Laštovička 2013). The derived trends qualitatively agree with model simulations that show that they are the consequences of the enhanced greenhouse gases. However, some of the trends are still controversial or discrepant.
The key challenges now are to improve the accuracy of various parameters and to reduce the differences in trends between models and observations.
A summary of the CAWSES-II TG2 achievements is given by Laštovička et al. (2014).
1.2.3.3 TG3: Short-term solar variability and geospace Towards the understanding of the origin and emergency of solar magnetism, the development of improved mean field solar dynamo models that take into account physical ingredients such as mutli-cellular and/ or time varying meridional circulation, magnetic pumping, North-South asymmetry, irregularity in the cycles, or improved source function for the poloidal field, has been accomplished (e.g., Gopalswamy et al. 2015a, b, c). The first self-consistent models were developed for generation of magnetic wreaths at the base of the convection zone that have become unstable and started to rise to the top (surface) of the numerical domain (Nelson et al. 2011, Nelson 2013. Three-dimentional MHD numerical simulations successfully explained the dynamics of several CMEs and the associated shock waves (Shen et al. 2011).
The International Study for Earth-Affecting Solar Transients (ISEST) program was newly created within the frame of CAWSES Task Group 3. It continued during the VarSITI program, with the focused objective of improving the scientific understanding of the origin and propagation of solar transients and developing the prediction capacity of these transients' arrival and potential impact on the Earth.
1.2.3.4 TG4: Geospace response to variable inputs from the lower atmosphere The Task Group 4 focused on the geospace response to variable inputs from the lower atmosphere. This contains atmospheric waves generated by meteorological disturbances, their interaction with the background mean flow in the middle atmosphere and with the ionospheric plasma, and their relation to competing disturbances in the thermosphere and ionosphere generated by energy inputs from the magnetosphere. A comprehensive summary of Task Group 4 is given by Oberheide et al. (2015).
The major scientific findings during CAWSES-II Task Group 4 are for impacts from lower atmosphere to the middle and upper atmosphere and ionosphere through various atmospheric waves, i.e., gravity waves with periods from 5 min to a few hours; tides with periods of 6, 8, 12, and 24 h; and planetary waves with periods of days. Vadas and Crowley (2010) showed the features of the penetration of gravity waves generated from convective plume to the middle atmosphere and ionosphere. Tsugawa et al. (2011) reported clear concentric waves in ionospheric total electron content (TEC) map generated by catastrophic Tohoku earthquake and Tsunami on 11 March 2011. This event provided a strong visual evidence of gravity wave penetration to the ionosphere. The wavenumber 4 structure observed in the equatorial ionosphere by Immel et al. (2006) was recognized as a result of penetration of atmospheric tides to the E-and F-region ionosphere and electromagnetic coupling between E-and F-region ionosphere (e.g., England 2012). After the discovery of response of ionospheric TEC to the stratospheric sudden warming (SSW) induced by planetary wave collapse by Goncharenko and Zhang (2008), many evidences of planetary wave penetration to the thermosphere and ionosphere have been reported during the CAWSES-II period.
A summarizing CAWSES-II symposium was held in Nagoya in 2013, and selected papers from the symposium were gathered in the special issue of Earth Planets and Space at https://www.springeropen.com/collections/ ICS (overview: Yamamoto et al. 2016). Review papers of the four task groups and overall CAWAES-II (overview: Tsuda et al. 2015) were published in the special issue of Progress in Earth and Planetary Science which is freely available online at https://www.springeropen.com/ collections/cawses2reviewcollection.

VarSITI (2014-2018)
CAWSES coincided with the prolonged period of very low solar activity between sunspot cycles 23 and 24. CAWSES II started in 2009, the long-awaited beginning of cycle 24, and covered the rising phase of the cycle. It was even subtitled "Towards solar maximum." However, that solar maximum was a surprise to most of the community.
Both CAWSES and CAWSES II were mainly focused on the Earth's atmosphere: how it is affected by solar activity on short and long time scales, how anthropogenically induced changes in the atmosphere affect climate, how the solar effects in the atmosphere depend on a changing climate, and how upper atmospheric levels are influenced by lower levels.
Due to the observed decline in solar activity which puts an end to the Modern Grand Maximum, as shown in Fig. 1, by the time of the start of the next SCOSTEP's scientific program, more attention was directed to the Sun: are we facing the next grand minimum in solar activity or will this be just a regular secular minimum? In either case, what are the most extreme solar events that we can expect, and how will they affect the Earth's system and its technological infrastructure? In a period of low solar activity, what are the characteristics of geoeffective solar agents, and how does the terrestrial system respond to them?
Basically, all we know about the Sun and the geoeffective solar transients from space-borne and in situ measurements is for the period of the "space era," coinciding with the period of the grand maximum of solar activity in the second half of the twentieth century. The question arises whether the mechanisms and models based on these observations will hold true in a period of lower activity.
The definition of the VarSITI program was a community effort (Gopalswamy 2013). The primary input was a set of nine white papers received from the community spanning the interests of all the scientific unions and interdisciplinary bodies of ICSU represented in the SCO-STEP Bureau. Community input was also received during the International Symposium on Solar-Terrestrial Physics in Pune India in 2012 when a panel discussion was held on the next scientific program. After a thorough discussion among the Bureau members, a set of 30 experts from the SCOSTEP community was invited to a brainstorming forum organized by the International Fig. 1 Instrumental sunspot observations: group sunspot number (Hoyt and Schatten 1998)-blue line and international sunspot number (Clette et al. 2007)-red line. The shaded area indicates the "space era" which started with the launching of the first artificial satellites of the Earth for in situ measurements in the outer space Shiokawa and Georgieva Progress in Earth and Planetary Science (2021)

Project SEE: Solar Evolution and Extrema
During CAWSES II, physical models based on the solar dynamo theory, together with other approaches (statistics, precursors, correlations, neural networks, etc.), were employed for the first time in predicting the following sunspot cycle 24. A set of predictions of the amplitude of Solar Cycle 24 accumulated in 2008 ranged from very low to unprecedented high levels of solar activity (Fig. 2). The predictions formed an almost normal distribution, centered on the average amplitude of all preceding solar cycles (Pesnell 2012). The wide range of the predicted amplitudes, even derived using the same methods, demonstrated that at the time of the completion of CAWSES II, we were still far from understanding the physics of the solar cycle.
The goal of the project SEE was to answer the following questions: What is the largest solar eruption/flare possible? What is the expectation for periods with absence of activity? 3. Does our current best understanding of the evolution of solar irradiance and mass loss resolve the "Faint Young Sun" problem? What are the alternative solutions?

Forecasting solar activity
One of the main focuses of VarSITI program was the end of the modern grand maximum of solar activity and what we can expect after it. Different approaches are employed to forecast the future sunspot cycles: extrapolations and statistical relations between some characteristics of the average sunspot cycle and between consecutive cycles, precursor methods, data-driven fluxtransport simulations, machine learning, neural network, autoregressive models, etc. All prediction methods can be generically divided into three broad groups, including the majority of the above classifications: statistical or extrapolation methods, precursor methods, and modelbased methods. Extensive reviews of the prediction methods and some early forecasts of cycle 25 are given by Petrovay (2020). It should be noted here that a major event occurred during the VarSITI program though it was not part of it. After a series of "Sunspot Number Workshops" (http:// ssnworkshop.wikia.com/wiki/Home), on July 1, 2015, the Sunspot Index Data Center in Brussels terminated the more than 400 years long data series of the International relative sunspot number R Z and replaced it by a new entirely revised data series S N (Cliver et al. 2013;Clette et al. 2014). In this revised data series S N , the overall level of solar activity was significantly increased as compared to R Z . In the future solar activity forecasts, some authors still use the original R Z series, others use the new S N series, and the comparison of the predictions for cycle 25, as well as between cycle 24 and the predictions for cycle 25 is not straightforward. To avoid misunderstanding, for each forecast, we will indicate which series is used: V1 for the original R Z (Clette et al. 2007) or V2 for the recalibrated S N (Clette et al. 2014).
Extrapolation methods can theoretically yield very long-term predictions. An attempt to make a superlong-term forecast was made by Popova et al. (2018), where it is stated that two fundamental dependences were found for a solar dipole and a quadrupole, on the basis of which one can make a forecast for many cycles ahead. This forecast was an example of the methods based on spectral analysis in which it is assumed that the sunspot time series is a superposition of periodicities. This article was criticized by Usoskin (2018), which the authors disagreed with Zharkova et al. (2018).
Most of the predictions based on statistical and extrapolation methods are for a low cycle 25, similar (possibly somewhat higher of somewhat lower) to cycle 24. Some studies predict a very weak cycle (e.g., Covas et al. 2019), a cycle stronger than cycle 24 (Sarp et al. 2018), or stronger than both cycles 23 and 24 (Pala and Atici 2019).
In the precursor methods, the future sunspot maximum is predicted using its correlations with other heliophysical parameters observed before this maximum ("precursors"). Helal and Galal (2013) used a correlation between the number of spotless days and the upcoming solar maximum to estimate that Solar Cycle 25 will have R z = 118.2 (V1), peaking 4.0 years after the solar minimum. Yoshida (2014) used correlations between R z before minimum with the upcoming solar maximum, using the symmetries of the even/odd cycles to derive the prediction: R z = 115.4 ± 11.9 (V1). Based on the empirical correlation between the dipole moment during solar minimum and the strength of the subsequent cycle, Cameron et al. (2016) suggest that Cycle 25 will be of moderate amplitude, not much higher than cycle 24. Kirov et al. (2018) show that both the maximum sunspot number of a solar cycle and the minimum of geomagnetic activity preceding it can be predicted from the maximum value of the solar dipole magnetic field of the previous cycle which is observed well before the sunspot minimum. They estimated that the minimum geomagnetic activity between cycles 24 and 25 as quantified by the geomagnetic ap-index will be about 5, and the maximum sunspot number for cycle 25 about 50-55 (V1). Petrovay et al. (2018) focused on the rush-to-the-poles (RTTP) phenomenon in coronal green line emission. Considering various correlations between properties of the RTTP with the upcoming solar cycle, they found a correlation between the rate of the RTTP and the time delay until the maximum of the next solar cycle. On the basis of this correlation and the known internal regularities of the sunspot number series, they predict that, following a minimum in 2019, cycle 25 will peak in late 2024 at an amplitude of about 130 (V2). This slightly exceeds the amplitude of cycle 24, but it would still make cycle 25 a fairly weak cycle. Gopalswamy et al. (2018) used polar and low-latitude brightness temperatures from the Nobeyama Radioheliograph at 17 GHz as proxies to the polar and active region magnetic fields, respectively, and the locations of prominence eruptions as a proxy to the filament locations to study their time variation. They showed that the polar microwave brightness temperature is highly correlated with the polar magnetic field strength and the fast solar wind speed, and that the polar microwave brightness in one cycle is correlated with the lowlatitude brightness with a lag of about half a solar cycle. They used this correlation to predict the strength of the solar cycle 25: the smoothed sunspot numbers (V2) in the southern and northern hemispheres can be predicted as 89 and 59, respectively. These values indicate that cycle 25 will not be too different from cycle 24 in its strength. Gopalswamy et al. (2018) also combined the rush-to-thepole data from Nobeyama prominences with historical data going back to 1860 to study the north-south asymmetry of sign reversal at solar poles. They find that the reversal asymmetry has a quasi-periodicity of 3-5 cycles.
The model-based methods are based on the theory of solar dynamo. A very recent review of dynamo models of the solar cycle is presented by Charbonneau (2020). The flux-transport mechanism allows predictions by assimilating the observed magnetic field on the surface. Since the emergence of sunspot groups has random properties, making it impossible to accurately predict the solar cycle and strongly limiting the scope of cycle predictions, Jiang et al. (2018) developed a scheme to investigate the predictability of the solar activity over one cycle. When a cycle has been ongoing for more than 3 years, the sunspot group emergence can be predicted along with its uncertainty during the rest time of the cycle. The method gives maximum strength of cycle 25 in the range 125 ± 32 (2σ uncertainty range, V2), which is about 10% stronger than cycle 24 based on the mean value. Using an Advective Flux Transport (AFT) model, Hathaway and Upton (2016) predicted that cycle 25 will be similar in strength to the cycle 24, with an uncertainty of about 15%. AFT also predicted that the polar fields in the southern hemisphere would weaken in late 2016 and into 2017 before recovering. After two more years of observations, Upton and Hathaway (2018) examined the accuracy of the 2016 AFT prediction and found that the new observations track well with AFT's predictions for the last 2 years. The updated prediction for cycle 25 was that it will be slightly weaker than cycle 24, making it the weakest cycle on record in the last hundred years. As weak cycles are preceded by long extended minima, we may not reach the cycle24/ 25 minimum until 2021. Bhowmik and Nandy (2018) coupled 2D flux transport model to a 2D kinematic dynamo model and sequentially calculated the polar field of one cycle and the number of sunspots in the next cycle. They predicted the maximum of cycle 25 in 2023-2025 with the height of 109-139 (V2).
The general conclusion is that the beginning of the twentyfirst century will be characterized by one or two cycles with a fairly low or just low intensity. A more serious, Maunder-type decline of activity cannot be ruled out either, but in the next few decades, the level of solar activity is expected to be either average (more likely) or of the type of the Dalton minimum (less likely).
A review of the recent results on predicting the future solar activity in the framework of VarSITI program is given by Nandi et al. (this issue).

Extreme solar events
Extreme can be considered events on the tail of a distribution-for example, a CME with speed or magnetic field exceeding certain thresholds that very few observed so far events have exceeded . From the terrestrial point of view, extreme are solar events leading to disturbances in the terrestrial system that very few observed so far events have caused. It should be noted that there is no direct correspondence between the two types of extreme events: extreme geomagnetic storms are nearly always, but not always, produced by extreme solar events, and extreme solar events do not necessarily cause extreme geomagnetic storms (Schmieder 2018).
The ability of a solar active region to produce an extreme event depends of the maximum possible magnetic Shiokawa and Georgieva Progress in Earth and Planetary Science (2021)  potential energy that can be converted into magnetic kinetic energy. The visual indicators of active regions are sunspots.
The number and area of sunspots reflect the intensity of the solar magnetic fields in the solar active regions (Houtgast and van Sluiters 1948). Therefore, using historical data of active regions' areas, estimations can be made of the active regions' magnetic fields, respectively, of the magnetic potential energy, and of the ability of the solar active regions to produce extreme events. Figure 3 shows the cumulative distribution of the active regions' areas A for the period 1874-2016 . The cumulative number decreases slowly until the area reaches~1000 msh (millionths of solar hemisphere) and then decreases rapidly. The overall distribution can be fit to the Weibull's function, and its high area tail-also to a power law. An extreme active region (once in 100 years) has an area of~7000 msh according to the power law and6 000 msh according to the Weibull function. Figure 4 presents the cumulative distribution of CMEs' speed from SOHO/ LASCO (https://cdaw.gsfc.nasa.gov) with the average speeds of the various CME populations marked on the plot: metric type II radio bursts (m 2 ); magnetic clouds (MC); ejecta (EJ)-ICMEs lacking flux rope structure; interplanetary shocks (S) detected in the solar wind; geomagnetic storms (GM) caused by CME magnetic field or shock sheath; halo CMEs (Halo) that appear to surround the occulting disk of the coronagraph and propagating Earthward or anti-Earthward; decameter-hectometric (DH) type II bursts indicating electron acceleration by CME-driven shocks in the interplanetary medium; SEP events caused by CME-driven shocks; and ground level enhancement (GLE) in SEP events indicating the acceleration of GeV particles. The fastest CME observed by SOHO/LASCO occurred on November 10, 2004, with an average speed in the coronagraph field of view of 3387 km/s, and there are very few CMEs with speeds exceeding~3000 km/s. Figure 5a illustrates the occurrence frequency of CMEs with speeds exceeding certain values. From the Weibull distribution, once in 100 years an event may occur with V > 3800 km/s, and once in 1000 years-an event with V > 4700 km/s . Figure 5b gives the cumulative distribution and occurrence frequency of the CMEs' kinetic energy. The once in 100 years and once in 1000 years events have kinetic energies 4.4 × 10 33 and 9.8 × 10 33 erg, respectively, which are not much higher than the highest reported values. The CME with the highest kinetic energy of 4.2 × 10 33 erg was observed by SOHO/LASCO on September 9, 2005 (Gopalswamy 2018).

Flares
As in the case of CMEs, the extreme solar flares are defined as one in 100 years and one in 1000 years events. Figure 6 is the cumulative distribution of all flares recorded by various GOES satellites since 1969 in the 1-8 A°energy band.
For the flare intensity, the Weibull and the power law distributions give almost identical results: the 100-year and 1000-year event sizes are X43.9 and X101, respectively, according to the Weibull distribution, and X42 and X115 according to the power law distribution (Gopalswamy 2018). As mentioned above, the largest flare during the space age was registered by the GOES satellite on November 4, 2003, but as the GOES X-ray sensor saturated at a level of X17.4 for about 12 min during this event, it was suggested that the flare size should have been in the range X34-X48, with a mean value of X40 (Brodrick et al. 2005). Therefore, it is about the size of the 100-year maximum. The bolometric energy (the energy emitted across the whole electromagnetic spectrum) corresponding to an X100 flare is 10 33 erg (see e.g., Benz 2017). Flares with bolometric energies > 10 33 erg are considered super flares.
There is no agreement about how extreme a solar flare can be, and how often extreme flares can occur. The period of instrumental observations of the Sun is only 400 years which is far from good statistics. However, we can complement these observations with observations of solar-like stars. Superflare stars are generally characterized by larger chromospheric emissions than other stars, including the Sun, but superflare stars with activity levels lower than, or comparable to, the Sun do exist, suggesting that solar flares and superflares most likely share the same origin (Karoff et al. 2016). Therefore, if we observe 10,000 solar type stars (slowly rotating G-type main sequence stars with surface temperatures between 5600 and 6000 K) for 1 year, we can get the data similar to the data obtained from 10,000 years' observations of the Sun. Maehara et al. (2012) using Kepler satellite data, which include observations of 80,000 solar type stars, found 365 superflares with total bolometric energy from 10 33 to 10 36 erg on 148 solar type stars. Superflares of 10 34 erg on Sun-like stars were estimated to occur once every 800 year, and those of 10 35 erg-once every 5000 year. This is much more frequently than implied by Fig. 3 where the tail of the flare-size distribution suggests that super flares (> 10 33 erg) can occur on the Sun once in a millennium, and a 10 34 erg flare-only once in 125,000 years . Besides, historical surveys of bright or oddly colored Sun, short intense heat waves, or low-latitude auroras indicate that the Sun has likely not had any superflares in the last two millennia (Schaefer et al. 2000). Aulanier et al. (2010Aulanier et al. ( , 2013) used a dimensionless numerical 3D MHD simulation for solar eruptions. When the model is scaled to the strongest measured sunspot group's area and magnetic field (6000 msh and 3.5 kG, respectively), it can be used to calculate the size of the sunspot pair that is required to generate the solar flares of various energies. Further, observations show that large sunspots groups are always fragmented into several spots and are never involved in a given flare as a whole, only one bipole will be involved in the flare. Aulanier et al. (2013) estimated that only 30% of the area of a given sunspot group can be involved in a flare. Figure 7 from Schmieder (2018) illustrates the size of the sunspot group needed to host a bipole which can produce a flare of a given size. With these scalings, a superflare of 10 36 erg requires a very large single pair of spots, whose extent is 48 o in longitude/latitude, at the surface of a Sun-like star. Such spots have never been reported on the Sun (Aulanier et al. 2013).
These scalings resulted in a maximum solar flare energy of∼6 × 10 33 erg which is only about six times that of the strongest-ever directly observed flare of Nov. 4, 2003. Aulanier et al. (2013 argued that superflares of 10 34 erg are unrealistic for the Sun because of the observed fragmentation of the magnetic flux in the active regions. However, Toriumi and Takasao (2017) demonstrated that even one of the largest, highly fragmented sunspot groups, such as the one numbered 14585 by the Royal Greenwich Observatory (RGO), produced a flare eruption leaving AR-sized, gigantic flare ribbons (on July 25, 1946), which may point to the possibility that even larger ARs could occur and cause a superflare. For the largest sunspot group on 1947 April 8, RGO 14886 with a maximum spot area of 6132 msh, Toriumi and Wang (2019) estimated its total flux to be 2.0 × 10 23 Mx, leading to a flare with energy 1.5 × 10 34 erg. This is more than twice larger than the maximum possible energy estimated by Aulanier et al. (2013), and in the energy range of the weak superflares reported by Maehara et al. (2012) for numerous slowly rotating and isolated Sunlike stars, but is several orders of magnitude smaller than that of strong stellar superflares.
An open question is whether the current solar dynamo can produce much larger sunspot groups as required to produce stronger flares.

Terrestrial perspective
From the terrestrial perspective, a crucial question is how extreme geomagnetic storms the solar extreme events can produce. The Dstindex is available since 1957, and since then, there has been only one storm with Dst ≤ 500 nT. It had an intensity of Dst = − 589 nT and happened on March 13-14, 1989. Another very intense magnetic storm, reaching almost the level of a superstorm with Dst ≤ 490 nT, occurred on November 20, 2003. Only these two events can qualify as possible superstorms since the Dst-index is available.
Before the Dst-index became available, the geomagnetic storms can be estimated from magnetic field data from the Colaba and Alibag Observatories operating since 1841. The strongest geomagnetic storm since then is the famous Carrington event on September 1, 1859 (Lakhina and Tsurutani 2018). The depletion of the horizontal component of the magnetic field measured in the Colaba observatory in Mumbai, India, was 1600 nT. The transit time of the ICME from the Sun to the Earth was 17 h and 40 min, so the average shock transit speed is found to be 2380 km/s, and the estimated speed at Earth's orbit was estimated to have been 1850 km/s. Based on the empirical relation between the speed and magnetic field of the ejecta (Gonzalez et al. 1998), the magnetic field magnitude was calculated to have been 90 nT at 1 AU. Finally, from the empirical expression for the dependence of the Dst-index of the speed and magnetic field of the ejecta, the Dst-index of the Carrington event was estimated as − 1760 nT (Tsurutani et al. 2003). Tsurutani and Lakhina (2014) and Lakhina and Tsurutani (2018) estimated the maximum possible intensity of a superstorm. They assumed a maximum CME speed close to the Sun of 3000 km/s, and a 10% decrease on the way to the Earth, or a speed 2700 km/s at 1 AU. This would give shock transit time from the Sun to the Earth of~12.0 h (for comparison, the fastest ever recorded event on August 1972 took 14.6 h and the Carrington event 17.6 h). This would mean Alfvén Mach number of~63 and magnetosonic Mach number of~45. The largest magnetosonic Mach number of the shock observed so far is~28 and that is for the shock associated with the extreme ICME of July 23, 2012 . The ram pressure of the maximum intensity superstorm downstream of the ICME shock was calculated to be 244 nPa, which would cause the magnetopause to be pushed inward from its quiet time position of 11.9 RE to a new subsolar position at~5.0 RE from the center of the Earth, where RE is an Earth radius (6371 km). So far, the lowest magnetopause position detected is at 5.2 RE for the August 1972 storm. Using the empirical relationship between the speed and magnetic field strength of the ICME at 1 AU (Gonzalez et al. 1998), the magnetic cloud field strength of 127 nT was estimated by Tsurutani and Lakhina (2014). This yields a maximum strength of the interplanetary electric field of 340 mV/m, and if we accept a linear dependence of the intensity of all magnetic storms on the interplanetary electric field (Burton et al. 1975;Echer et al. 2008), the maximum possible intensity of a superstorm is expected to be twice the intensity of the Carrington storm, that is, Dst~3500 nT. Figure 8 shows the cumulative distribution of geomagnetic storms with Dst < 100 nT. According to the Weibull distribution which fits better the data points, a 100-year event has a size of 603 nT, and a 1000-year event has a size of 845 nT. The 100-year event is consistent with the March 1989 event, while the 1000-year event is half the estimated size of the Carrington storm. It should be noted, however, that some authors argue that the − 1760 nT of the Carrington event is an overestimation and suggest a Dst around − 900 nT as a nominal value for it . Yermolaev et al. (2018), assuming power law, calculated the tails of distributions in the area of moderate and strong storms (from − 150 to − 300 nT) where statistics of events are sufficiently high in the period 1963-2015 and then extrapolated the distributions to the area of extreme magnetic storms where statistics are poor. They found that the most probable waiting times for extreme magnetic storms with Dst ≤ 500, − 1000, and − 1700 nT are 24 years, 250 years, and 1500 years, respectively (with errors 10%, 40%, and 100%). If only the storms caused by magnetic clouds are considered, the waiting times reduce to 20, 120, and 500 years, respectively. A detailed review of the SEE results is given by Nandi et al. (2021).

Highlights of contribution from young scientists to the SEE project
As has been shown in the highlight of young scientists in VarSITI Newsletter, several active young scientists have contributed to the VarSITI/SEE project. Toriumi et al. (2014) studied solar flux emergence and formation of flare-productive active region on solar surface using both satellite observation and three-dimentional simulation. Jiang et al. (2014Jiang et al. ( , 2015 studied effects on sunspot group tilt angles on large-scale magnetic field on the solar surface, and relation of the solar surface field evolution to the weak solar cycle 24. Pal et al. (2017Pal et al. ( , 2018 reported Sun-to-earth analysis of magnetic helicity of the 2013 March 17-18 interplanetary coronal mass ejection and the dependence of coronal mass ejection properties on their solar source active region characteristics and associated flare reconnection flux. Holappa et al. (2014a, b) showed a long-term relationship between solar wind structures and geomagnetic activities using new principal component analysis. For solar effect on the Earth's atmosphere, Maliniemi et al. (2014Maliniemi et al. ( , 2016 reported spatial distribution of Northern Hemisphere winter temperatures during different phases of the solar cycle and the effect of geomagnetic activity on the Northern Annular Mode in the lower atmosphere.

Project ISEST/MiniMax24: International Study of Earthaffecting Solar Transients
The International Study for Earth-Affecting Solar Transients (ISEST) program started already within the framework of CAWSES II Task Group 3 ("Short-term Solar Variability and Geospace"). The summary of this activity is given by Gopalswamy et al. (2015a). The project continued during the VarSITI program, focused on the origin and propagation of solar transients, and their impact on the Earth. Solar transients are short time-scale (from minutes to days) energetic and/or eruptive events occurring on the Sun that directly cause transient disturbances in the Earth's space environment and possibly impose adverse effects on advanced technological systems of human society.
Solar transient events can be divided in the following four categories: (1) solar flares, (2) coronal mass ejections (CMEs) including their interplanetary counterparts ICMEs, (3) solar energetic particle (SEP) events, and (4) stream interaction regions (SIRs) including corotating interaction regions (CIRs). These four types of Earth-affecting transient events differ in their physical origin or processes, observational appearances, and their geoeffectiveness.
This overview is a brief summary of the detailed description of the working group activities by Zhang et al. (2021).

Data
In the last decade, a large fleet of spacecraft provided continuous measurements of the Sun and the interplanetary space. Unprecedented multiple-angle observations were provided by the Solar Terrestrial Relations Observatory Ahead/Behind (STEREO A/B). Thanks to these observations, our knowledge on CME morphological and kinematics properties has improved significantly.
Several studies focused on Sun-to-Earth analysis of CMEs by linking the low coronal behaviour and properties of the eruption with its observed in situ signature (D'Huys et al. 2017;Temmer et al. 2017). A number of studies also compared magnetic flux rope properties, such as the magnetic flux, the chirality, and its helicity sign determined from post-eruptive arcades, flare ribbons, and coronal dimmings measured close to the Sun with magnetic cloud properties at 1 AU (Qiu et al. 2007;Gopalswamy et al. 2017;Palmerio et al. 2017;James et al. 2017).
The weak solar cycle 24 was one of the main focuses of the VarSITI program. The properties of the transient events were extensively studied in the framework of ISEST/MiniMax. The overall rate of CMEs was found to have increased in solar cycle 24 relative to cycle 23, although the rate of fast and wide CMEs decreased. Marked changes are observed in the number of intense (Dst ≤ 100 nT) geomagnetic storms (Gopalswamy et al., 2015a, b, c) and high-energy (≥ 500 MeV) SEP events (Gopalswamy et al. 2014;Mewaldt et al. 2015). On the other hand, the number of halo CMEs in cycle 24 did not decrease significantly (Gopalswamy et al. 2015b). In fact, the number of halo CMEs as well as the total number of CMEs normalized to the sunspot number is larger in cycle 24 than in cycle 23.

Theory
Scolini et al. (2019) used proxies of magnetic flux estimates determined from post-flare arcades , flare ribbons (Kazachenko et al. 2017;Tschernitz et al. 2018), and coronal dimmings (Dissauer et al., 2018a, b) as initial input for the EUropean Heliospheric FORecasting Information Asset (EUHFORIA) model, to study the geoeffectiveness of the famous 2017 September events. Good agreement with the observed Dst profile was found for simulations using the optimized input and including CME-CME interactions.
Recently, there has been substantial development in the forward modelling of the magnetic structures using physics-based empirical models. The forward models typically assume a specific morphology of a magnetic structure, i.e., flux rope, and evolve it, assuming specific propagation and expansion (Wood et al. 2017). A step further is to include a specific flux rope magnetic field topology (e.g., Patsourakos et al. 2016;Isavnin 2016;Kay et al. 2017;Möstl et al. 2015).

Simulation
During the period of the VarSITI/ISEST, the main focus in the MHD space weather modeling was in the following areas: the background solar wind, CME initiation and propagation, interaction between CMEs and solar wind structures (CIR and HCS), and CME-CME interaction.
A detailed review of the recent progress in these areas is given in Zhang et al. (2021), and the main conclusions are as follows: There has been a significant increase in the number of 3-D MHD codes that has been successfully used to simulate the Sun-to-Earth propagation of CMEs; this has been the case most notably in Europe with EUHFORIA, Shiokawa and Georgieva Progress in Earth and Planetary Science (2021)  Heliospheric codes (starting typically at 0.1 AU) have been used with spheromak and/or flux rope CMEs, which bridges the gap between computationally intensive Sun-to-Earth simulations and heliospheric simulations with code models. These types of simulations may be used to investigate the magnetic field configuration inside CMEs as well as their arrival time, and are more physically consistent when investigating CME-CME interaction than simulations in which the CMEs do not have internal magnetic fields.
In parallel, there has been an effort to make the CME initialization quicker and easier to perform in coronal codes using out-of-equilibrium flux ropes. This paves the way for future real-time Sun-to-Earth simulations with magnetized CMEs initiated based on magnetograms, EUV images, and early coronagraphic images. It is well known that major changes in the CME properties, including its speed and orientation, may occur below 0.1 AU where heliospheric models are initiated. At this time, it is however unclear whether simulations with magnetized CMEs initiated at 0.1 AU using multi-viewpoints coronagraphic measurements will perform worse than simulations with magnetized CMEs initiated at the solar surface, in term of space weather forecasting capabilities. The number of Sun-to-Earth simulations of CMEs initiated at the solar surface with a realistic model is still relatively low, even though there has been effort in presenting the results near 1 AU of more complex initiation mechanisms.
There has been significant new physics included in the solar wind models, including more advanced thermodynamics treatment and the inclusion of Alfvén waves. There has not been significant work quantifying how these new additions affect the CME propagation and the resulting structure near 1 AU.
There has been progress towards coupling timedependent magnetic field models with coronal models and heliospheric models. Initiating CMEs by means of magnetofrictional or other self-consistent models based on solar observations or flux emergence may lead the way to a better physical understanding of CMEs and is also the only way space weather forecasting could provide information before the launch of a CME. Further improvements towards this coupling are expected in the next few years.

Campaign study
The task of ISEST Working Group 4 (campaign events) was to integrate theory, simulations, and observations for better understanding of the chain of activity from the Sun to the Earth for a few carefully selected events. One focus of WG 4 was on why forecasts fail and how we can improve our predictions. This included analyzing the complications in linking CMEs to ICMEs, usually observed only in situ at 1 AU.
After examination, WG 4 classified the studied events into three general categories (Webb and Nitta 2017) as follows: 1. "Textbook" cases in which the complete chain of well-observed events is relatively well understood from its solar source, through its heliospheric propagation, to its geo-effects. These cases involve forecasts that are successful in a general way. 2. Cases in which there were problems understanding the complete chain, but which we think we now understand it. Thus, something was missing in the chain of a well-observed event but, in retrospect after analysis, we now know why. These cases usually involve forecasts that failed because the events were not geoeffective or were underestimated. 3. Problem cases in which the chain is not complete, and we still do not understand why.
Fourteen campaign events were discussed and analyzed by WG 4 (Table 5-1 in Zhang et al. 2021). Figure 9 is an illustration of the first GLE event observed on September 10, 2017, on two planets: Earth and Mars. A detailed description of the event and the modeled geometry, kinematics, propagation, and interaction of the three CMEs launched around 2017-09-10 from their solar origin to their arrivals at Mars and Earth are given by Guo et al. (2018).

SEP events
3.2.5.1 SEP origin One of the more unusual solar phenomena is the long-duration gamma-ray flares (LDGR Fs). The prime characteristic of these events is delayed and prolonged γ-ray (> 100 MeV) emission after the impulsive phase. Recently, the Fermi Large Area Telescope (LAT) observed dozens of LDGRFs. Share et al. (2018) characterized and catalogued 30 solar eruptive events observed by Fermi/LAT from 2008 to 2016, referring to this emission as "late-phase gamma-ray emission" (LPGRE). In all the events studied, Share et al. (2018) found that the LPGRE is temporally and spectrally distinct from the impulsive phase emission and concluded that due to the broad range of the LPGRE characteristics it is unlikely that the mechanism of acceleration of particles by the CME-driven shock alone is able to account for the emission in all of the events under study.
3.2.5.2 SEP acceleration Some studies have provided new support for a significant contributory or dominant role for flare acceleration of high-energy protons in gradual SEP events, contrary to the generally accepted scenario favoring shock acceleration (e.g., Dierckxsens et al. 2015;Grechnev et al. 2015). Several other recent studies support the prevailing shock picture for gradual SEP events, such as the observation of a prompt SEP event at widespread locations in conjunction with the longitudinal propagation of a white-light shock (Lario et al. 2016), the SEP source temperatures (Reames 2015), and the hierarchical relationship between the fluence spectra of gradual SEP events and the kinematics of the CMEs.

SEP transport
The unprecedented orbits of the two-spacecraft STEREO mission provided well-separated observations at 1 AU and allowed to study the longitudinal distribution of SEPs and especially events with extraordinarily wide particle spreads in great detail. The STEREO mission has also enabled the identification and study of some extreme cases of widespread events with distributions up to 360°around the Sun.

SEP compositional results
Measurements of relative abundances of heavy elements and their isotopic and ionic charge state composition in SEP events have been used in a wide variety of ways to infer critical information about the origins of the seed populations and the physical conditions under which these populations are produced. These seed particles are accelerated by CME shocks or in solar jets and indicate the conditions of their transport through the solar corona, the interplanetary medium, and out into the heliosphere (e.g., see Reames 2015; Desai and Giacalone 2016).
3.2.5.5 Progress in SEP theory and modeling Valuable insight on particle acceleration at coronal shock waves Fig. 9 In situ observations of the 2017-09-10 events at Earth (left panels), Mars (f-h), and STA (i-k). a, i The energy-dependent particle flux (counts/cm 2 /sr/MeV/sec) measured by GOES and STA. b Count rate of several ground-based neutron monitors with different cut-off rigidities (each neutron monitor data set is normalized to the average value of the selected time range and shifted apart in y-axis). c, h, and j The solar wind speed (black, left Y-axis (yl)), proton density (blue, right Y-axis (year)), and temperature (red, year) for Earth (OMNI), Mars (MEX), and STA, respectively. The ENLIL modeled results at three locations are also plotted as dashed lines. d Plots of the Alfvén Mach number (black, yl), plasma (blue, year), and flow pressure (red, year) estimated at Earth. e, k The vector magnetic fields at Earth/STA in Geocentric Solar Ecliptic (GSE) or spacecraft Radial-Tangential-Normal (RTN) coordinate and negative ENLIL modeled magnetic field strength (dashed line). f The normalized count rate for downward particles stopping in RAD and penetrating RAD with the former/later approximating protons with energies larger than 175/ has been recently obtained by studying the evolution of CME and shocks in the corona, as well as their interaction with underlying magnetic fields and coronal plasma, both through data-driven or analytical modelling, as well as MHD simulations and combinations of different approaches. New information has been gathered about the relevant parameters for efficient shock acceleration (such as the Mach number, compression ratios, and geometry of shock waves), the primary acceleration regions along the shock, the role of coronal magnetic field configuration, and how these factors are related to the particle spectra observed in space.
The effect of large-scale streamer-like magnetic configuration on particle acceleration at coronal shocks has been investigated by considering a CME-driven shock propagating through a streamer-like magnetic field (Kong et al. 2017) or from its flank (Kong et al. 2019), i.e., when the streamer is rotated with respect to the CME propagation direction. By numerically solving the Parker transport equation with both parallel and perpendicular diffusion in such configuration, Kong et al. (2019) found that the primary sources for particle acceleration are located at different regions and vary significantly as the shock propagates and expands, depending on the particle energy and time.
The well-separated SEP observations with the STER EO mission have also lent themselves to study the transport of SEPs in more detail. The comparison of multispacecraft observations with results of 2D or 3D models solving the focused transport equation allowed to study not only transport along the mean magnetic field but also perpendicular to it. It was found that the role of efficient transport perpendicular to the mean magnetic field can eventually be much stronger than expected (e.g., Dresing et al. 2012) and might play an important role, among extended injection and acceleration regions, in the longitudinal spreading of SEPs (e.g., Dröge et al. 2014), but also in creating asymmetries in their longitudinal distribution at 1 AU (He and Wan 2015).
An important step taken in SEP transport modelling is the inclusion of a realistic solar wind background. Wijsen et al. (2019) used the data-driven EUHFORIA (Pomoell and Poedts 2018) to generate a background solar wind for their SEP transport code. This allows, on the one hand, the study of the effect of solar wind streams on the SEP propagation and also on adiabatic energy changes. On the other hand, case studies with a realistic background field will be possible accounting for the effects of transient structures leading to non-Parker field configurations.

Stream Interaction Regions/Co-rotating Interaction Regions (SIRs/CIRs)
The evolution of long-lived coronal holes, closely related to SIRs/CIRs, can be studied in detail using multiple views on the Sun from combined STEREO and Earth imagery. It was found that coronal holes undergo evolutionary patterns revealing a growing and declining phase where area increases and decreases again over several solar rotations (Heinemann et al. 2018(Heinemann et al. , 2019. The coronal hole area is found to be well related to the solar wind speed measured in-situ at 1 AU (e.g., Vršnak et al. 2007;Rotter et al. 2012). With this well-known areaspeed relation, empirical forecasting tools for the "pure" background solar wind on the basis of coronal hole area measurements are performed on a regular basis (Vršnak et al. 2007;Rotter et al. 2012). By understanding the photospheric and coronal evolutionary characteristics of coronal holes, one can aim to gain a better understanding and in turn improve the forecast of CIRs (Heinemann et al. 2018). The forecasts are found to be most successful for periods of low solar activity, as during increased solar activity transient events, such as CMEs, strongly disturb the rather stable solar wind outflow for several days (Temmer et al. 2017). As the high-speed solar wind streams emanate from coronal holes which are low dynamic structures, CIRs can be forecast with long lead times. Based on that, a variety of persistence models were developed. Under the assumption of persistence, in situ measurements of the solar wind plasma flow from L1 and varying STEREO vantage point provide a forecast for Earth position with lead times of up to 27 days (depending on the exact STEREO spacecraft position). Implementing the actual changes of coronal hole areas (from EUV data) into such simple forecasting tools can improve the forecast quality . Lead times with about 4.5 days could be achieved when using data from an instrument permanently located at the Lagrangian point L5 (60°behind Earth; ESA preparation for the future L5 mission "Lagrange").
Forecasting solar wind structures in interplanetary space serves also as important information for analytical CME propagation models and space weather models. Therefore, under the Space Situational Awareness Program of the European Space Agency, forecasting services using empirical, and numerical models for the solar wind are available (see http://swe.ssa.esa.int). As the forecasting/nowcasting quality is still not sufficient for producing reliable space weather alerts, we need to better understand and closely monitor SIRs/ CIRs. Main aims for the near future are to verify and evaluate background solar wind models, and with that improve the input for CME propagation models, to predict periodic and recurrent geomagnetic effects from CIRs more accurately. In general, single events are easier to forecast compared to multiple events covering CME-CME interactions or interactions of CMEs with CIRs. By better understanding the physics behind the coronal hole evolution, we may improve their forecasting capability. With that we will also gain more insight in the ejection and acceleration processes that define high-speed streams. For this we can exploit data from the NASA mission Parker Solar Probe (PSP) measuring the near-Sun space and with that regions where the solar wind actually gets accelerated.
Though harder to evaluate and forecast, interaction events between SIR/CIR and CMEs are of special interest. Interaction events can lead to significant increase in geomagnetic effects when compared to individual events of similar strength (e.g., He et al. 2018).
In a recent paper by Jian et al. (2019), physical properties of a large sample of slow-to-fast SIRs were investigated using STEREO-A and -B data. They identified 518 pristine SIRs, of which more than 50% are associated with crossings of the heliospheric current sheet (HCS) and are of slow speeds but higher densities, and of increased dynamic and total pressure compared to those without HCS. In that respect, HCS-related SIRs can be classified as more effective in terms of space weather.
The detection and extraction of reliable coronal hole areas from operational solar observations is extremely important, not only for solar wind forecasting using the area-speed relation, but also for investigating the magnetic open flux on the Sun. At present, most extraction methods focus on EUV observation taken by SDO/AIA, SOHO/EIT, GOES/SUVI, and/or STEREO/EUVI. A multi-wavelength approach was developed by Garton et al. (2018) in the form of the multi-thermal emission recognition algorithm CHIMERA. Recently, with the dawn of machine learning, new methods, utilizing the increased computational performance have also emerged to provide an additional tool to identify and extract coronal holes (e.g., Illarionov and Tlatov 2018).
To improve solar wind models and to ensure accurate space weather forecasting, the solar wind models have to be tested and validated. The validation of solar wind models is done by comparing the simulation results with in situ measurements.
The performance assessment of the EUHFORIA solarwind model was analyzed by Hinterreiter et al. (2019). The solar wind was modeled rather well for times of solar minimum but not during increased solar activity. In a study by Lee et al. (2009), the heliospheric models ENLIL/MAS and ENLIL/WSA were compared with in situ measurements from ACE andWind (time range: 2003-2006). They found that the model results give lower densities for faster solar wind fully agreeing with the solar wind momentum flux constancy.

Forecasting CMEs
The lack of critical observations (e.g., no routine observations of the magnetic field in the corona) and limitations in theory and models (e.g., idealized initial and boundary conditions) is currently not allowed to predict when a CME would occur. However, thanks to advances in our observational knowledge (e.g., new observations from STEREO, Hinode, SDO, IRIS), and in modeling (e.g., increase of realism in models, data-constrained and data-driven models) and in analysis and forecasting techniques (e.g., use of advanced statistical tools and machine learning methods) significant progress in our understanding and eventual prediction of CMEs has been achieved over the last decade. Tziotziou et al. (2012) found that helicity and magnetic free energy thresholds of~2 × 10 42 Mx 2 and~4 × 10 31 erg separate eruptive with non-eruptive ARs. Magnetic twist (Tw) (i.e., number of turns of field lines around a magnetic axis) of magnetic field lines is another parameter that is extensively used in CME onset studies. Magnetic twist is an integral part of magnetic helicity discussed above, and comparison of its properties/distribution in the pre-eruptive/eruptive configurations in the solar atmosphere and at 1 AU supplies important physical clues about CMEs .
The problem of predicting major CMEs could be mitigated to predicting major flares. Recent reviews on flare forecasting, including machine-learning schemes, can be found in Leka et al. (2019) and Park et al. (2020). Major conclusions from their extensive benchmarking of a large number of methods currently used in flare forecasting are that numerous such methods do better than climatology, no method clearly outperforms the others, and consideration of prior flare history improves the corresponding skill scores.
Machine learning has been used directly in CME predictions. Bobra and Ilonidis (2016) applied a Support Vector Machine classification scheme to 18 parameters derived from HMI vector magnetograms for more than 3000 ARs and found that only 6 amongst these parameters are sufficient to separate erupting and non-erupting ARs within 24 h from the corresponding measurements. These parameters (e.g., mean gradient of the horizontal magnetic field, mean current helicity, mean twist parameter, etc) are intensive (i.e., do not depend on the AR size but are spatial averages) and not extensive (i.e., dependent on the AR size and corresponding to spatial sums). Guennou et al. (2017) by means of MHD analyzed a set of eruptive and non-eruptive MHD simulations and found, in agreement with Bobra and Ilonidis (2016), that intensive parameters are more relevant to eruptivity. The availability of multi-viewpoint STEREO imaging observations allowed to monitor a number of modulations that CMEs experience in the corona and in the inner heliosphere in terms of deflections and rotations. CME rotations as well as deflections in the corona are simultaneously addressed by the ForeCAT model of Kay et al. (2015). The model has been applied to several CMEs, and it was able to reproduce the deflections and rotations that the CMEs undergone as derived from the analysis of the corresponding STEREO observations. The predicted deflections and rotations exhibit significant sensitivity on the input parameters and the employed background models.
Predicting the magnetic field distribution of CMEs/ ICMEs is of highest importance in heliophysics, as extended intervals of intense southward magnetic fields are a necessary condition for inducing geomagnetic storms. A thorough account of the state-of-the art in this important problem was recently given in Vourlidas et al. (2019).

Minimax24 non-flare target
The aim of the ISEST/MiniMax24 non-flare target service is to monitor potentially geoeffective phenomena which are not related to solar flares and to email alerts. This service, provided daily by the observer on duty, was first established in the scope of the SCOSTEP/CAWSES II "MiniMax24 Campaign" in 2013, which was declared as the year of "MiniMax24" to note that, even though the Sun is going through activity maximum conditions, the activity is rather low. The goal of the action was to understand and explain the current behavior of the Sun and its potential impact on human society and Earth's space environment through year-long scientific and outreach activities. The campaign team counted 37 institutions from 17 countries focused on the solar-terrestrial observations of solar eruptive events through the Max-Millenium program of solar flare research (http://solar. physics.montana.edu/max_millennium/) as well as coronal holes, filaments, and CIRs (i.e., non-flare related phenomena) through the newly established email alert service. By the end of the year of "MiniMax24," Mini-Max24 email list reached more than 140 participants from more than 30 countries. After the end of CAWSES II program, the action transcended from its original 1year-campaign scheme and was included in the new SCOSTEP program VarSITI as one of the working groups of the ISEST project.
3.2.9 Highlights of contribution from young scientists to the ISEST/Minimax24 project Many young scientists have been glowing to contribute to the ISEST/Minimax24 project. Thalmann et al. (2016) reported plasma and magnetic environment of a particular flare signature and coronal magnetic field based on combined SDO, RHESSI, and STEREO data. Cheng et al. (2014) made a tracking of magnetic flux rope from inner to outer corona. Sindhuja and Gopalswamy (2020) investigated properties of flux ropes in CME near the Sun.  also investigate evolution of 12 July 2012 CME from Sun to Earth. Kay et al. (2015Kay et al. ( , 2016Kay et al. ( , 2017 studied global trends of CME deflections, evolution of CME, and predictability of the magnetic field of earth-impacting CMEs, respectively. Shen et al. (2017Shen et al. ( , 2018 investigated geoeffectiveness of interplanetary coronal mass ejection (ICME) statistically from 1995 to 2014 and for a particular ICME event on September 8, 2017. Dumbović et al. (2018) developed an analytical model of Forbush cosmic-ray decreases caused by flux ropes in the solar wind associated with interplanetary coronal mass ejection (ICME). Dissauer et al. (2018aDissauer et al. ( , b, 2019 extensively investigated coronal dimmings associated with flares and CMEs.

Project SPeCIMEN: Specification and Prediction of the Coupled Inner-Magnetospheric Environment
The coupled inner-magnetospheric environment has been focused by the project "Specification and Prediction of the Coupled Inner-Magnetospheric Environment" (SPeCIMEN). The goal and objectives of this project were the quantitative prediction and specification of the Earth's inner magnetospheric environment based on Sun/solar wind driving inputs. Comprehensive review on the achievements of this project is provided by Kanekal and Miyoshi (2021). Here, we briefly review the highlights of the scientific results obtained on this topic. Figure 10a shows basic configuration of energetic electrons and ions and the waves that interact with these particles in the inner magnetosphere around the earth. The energetic electrons and ions rotate around the Earth due to curvature and gradient of ambient magnetic field with a time scale of minutes to hours and interact with extremely low frequency (ELF) and very low frequency (VLF) waves (electron cyclotron waves) at frequencies of kHz, ultra-low frequency (ULF) (Pc1) waves (ion cyclotron waves) at frequencies of~Hz, and ULF (Pc4-5) waves at frequencies of~mHz. These energetic particles Ground network instruments have been deployed extensively at subauroral latitudes by the PWING project since 2016 ). PWING stands for study of dynamical variation of Particles and Waves in the INner magnetosphere using Ground-based network observations. Modeling efforts to cover the inner magnetosphere has also been made, such as Block-Adaptive-Tree-Solar-wind-Roe-Upwind-Scheme with Comprehensive Ring Current Model (BATSRUS/CRCM, Glocer et al. 2013) and Ring current-Atmosphere interaction Model (RAM, Jordanova et al. 2010). These extensive observation networks and advanced modeling efforts make a golden era of research for the inner magnetosphere during the VarSITI program. Figure 11 shows an example of 1-year electron flux variation in 2013 at energies from 46 keV to 1.553 MeV in the inner magnetosphere observed by the Van Allen Probes (Reeves et al. 2016). The vertical axis is the L value, which is the radial distance from the Earth with a unit of Earth's radius. We can see dynamic variation of electron fluxes in the inner magnetosphere. The amplitude of the variation reaches a few orders of magnitudes. The enhancements of electron fluxes clearly correlated with solar wind speed enhancements and decreases of Dst index (geomagnetic storms). The high-flux regions are separated into two regions around L~4-6 and L < 2, forming outer and inner radiation belts, respectively. Clear differences can be also seen in the electron behavior at different energies and from event to event. Basically, no significant electron fluxes are seen at energies above 1 MeV in the inner belt at L value less than 3, as newly found by Fennel et al. (2015), while the inner belt clearly appears at energies below 459 keV. At further high energies, Van Allen Probes also discovered sharp inner boundary for the ultrarelativistic (energies higher than 5 MeV) electrons in the Earth's radiation belts . A comprehensive review of the results obtained by the Van Allen Probes mission has been provided by Baker et al. (2018).
The ERG satellite was launched in December 2016 and has provided several interesting observations particularly related to the wave-particle interactions in the inner magnetosphere. Figure 12 shows one-to-one correspondence between ELF/VLF chorus wave elements (panel c) that have a time scale of less than 1 s observed by ERG with electrons at energies of kilo-electronvolts that causes visible flush aurora observed by ground an all-sky camera (panels b and f) (Ozaki et al. 2019). ERG and the ground all-sky camera were nearly on the same geomagnetic field line. Clear one-to-one correspondence can be seen at four ELF/VLF chorus wave elements in panel (c) with the four auroral flashes in panel (f) at1 3:01:28 UT and~13:01:32 UT. This correspondence indicates that the observed ELF/VLF chorus waves scatter electrons along the magnetic field line to cause precipitation of the electrons and the flash aurora. Kasahara et al. (2018) did show one-to-one correspondence between ELF/VLF chorus waves and electron fluxes in the loss cone, which is the direct evidence of pitch-angle scattering of keV-energy electrons by ELF/VLF chorus waves. On the other hand, Teramoto et al. (2019) provided evidence of wave-particle interaction between ULF (Pc4-5) waves at frequencies of mHz with energetic electrons through drift resonance based on multisatellite coordinated observation by the Van Allen Probes and ERG satellites. More comprehensive results obtained from the ERG mission have been provided by two special issues in Earth Planets and Space in 2017 and in Geophysical Research Letters in 2018.
In the objectives of SPeCIMEN (specification and prediction), the modeling takes the part of quantitative prediction of the Earth's inner magnetospheric environment. BATSRUS/CRCM model (Glocer et al. 2013) provides global response of the outer and inner magnetosphere to the solar wind. RAM simulation (Jordanova et al. 2010(Jordanova et al. , 2012 provides distribution of whistler-mode chorus waves that contribute acceleration and loss of the radiation belt electrons. Geomagnetic Environment Modeling System for Integrated Studies (GEMSIS)-Ring Current and Radiation Belt models reproduce interaction of Pc5 magnetic pulsations with relativistic electrons (Kamiya et al. 2018). Not only by these physical models, artificial intelligence has been introduced for reconstruction of the inner magnetosphere dynamics. Figure 13 shows neural network reconstruction of the global plasma density as a function of L and MLT by Bortnik et al. (2016), using measurements by three THEMIS satellites between 2008 and 2014.

Highlights of contribution from young scientists to the SPeCIMEN project
More than ten active young scientists contributed to the researches of the SPeCIMEN project and wrote articles for highlight of young scientists in the VarSITI newsletter. For studies of whistler mode waves which possibly accelerate relativistic electrons, Li et al. (2014)   generated type of nonlinear wave structures. Píša et al. (2015Píša et al. ( , 2016 investigated Langmuir waves observed upstream of Saturn's foreshock and bow shock observed by Cassini, which also has some implications on highenergy plasma acceleration in the inner magnetosphere. For the magnetosphere-ionosphere coupling in the inner magnetosphere and ionospheric processes, Thomas et al. (2015Thomas et al. ( , 2016Thomas et al. ( , 2019 investigated magnetic field variations associated with Pi 2 and ionospheric current observed by the CHAMP and Swarm satellites just above the ionosphere and on the ground to show their characteristic in and near the ionosphere. McGranaghan et al. (2015aMcGranaghan et al. ( , b, 2016a extensively investigated high-latitude ionospheric conductivities based on satellite observations and developed a model of the high-latitude ionosphere to reconstruct its electrodynamic variations. Gokani et al. (2015) reported lightning-associated whistlers at very low latitudes (L = 1.08), suggesting ducted mode of propagation of these waves near the equator. McCormick et al. (2018) proposed spatial and temporal monitoring of the D-region ionosphere using lightning-induced broadband sferic measurements. These ionospheric topics also have strong relation with the ROSMIC project. The goals and objectives of the Project "Role of the Sun and the middle atmosphere/thermosphere/ionosphere in Climate (ROSMIC)" was to understand the impact of the Sun on the terrestrial middle atmosphere, lower thermosphere, and ionosphere and Earth's climate and its importance relative to anthropogenic forcing over various time scales from minutes to centuries. A comprehensive review on the achievements of this project is provided by Ward et al. (2021). Here, we briefly review the highlights of the scientific results obtained on this topic during the VarSITI program. The topic of ROSMIC can be divided to (1) coupling from above, (2) coupling from below, and (3) observation of long-term trends. The coupling from above comes from (1a) solar radiation in various wavelengths from X-ray, ultraviolet, and visible wavelengths, (1b) solar energetic particles (SEP) or solar proton events (SPE) and (1c) magnetospheric energetic particles (MEP). The effects of (1a) are often divided into total solar irradiance (TSI) and spectral solar irradiance (SSI). SEP and MEP are called as energetic particle precipitation (EPP). Figure 14 shows an example of the EPP effect on the middle atmosphere observed by the three different satellites for three different magnetic storms and in superposed epoch analysis based on satellite-based measurements by GOMOS, SABER, and MLS (Andersson et al. 2014). Both the particular storm events in Fig. 14a-c as well as superposed epoch analysis in Fig. 14d-i, there are significant ozone depletions at the middle atmosphere with the maximum loss of ozone from 10 to 30 % at altitudes between 70 and 78 km after the storms. This indicates a significant effect of ozone loss by EPP in the middle atmosphere. Funke et al. (2014) have shown intrusions of EPP-caused nitric oxide (NOy) into the stratosphere with a time scale of several months, based on measurements by Michelson Interferometer for Passive Atmospheric Sounding (MIPA S) on board the Envisat satellite. This gradual descent of NOy from the mesosphere to the middle and lower stratosphere contributes to the catalytic destruction of ozone. Efforts to quantify the total forcing of EPP to the atmosphere have also been made during the VarSITI interval. Isono et al (2014) reported NO column density enhancement 1-5 days after the commencement of geomagnetic storms based on ground-based millimeter-wave Together with these new measurements of EPP effects on the atmosphere, efforts to model the global effects have also been made. Andersson et al. (2018) combined the Whole Atmosphere Community Climate Model (WACCM) together with the EPP forcing model developed by van de Kamp et al. (2016) for the 6th phase of the Coupled Model Intercomparison Project (CMIP6) (Matthes et al. 2017). They concluded that mediumenergy (30 keV-1 MeV) electrons (MEE) can enhance the stratospheric ozone response to solar activity by a factor of 2. Thus, the EPP forcing can be significant for the long-term stratospheric ozone variation, implying their effects on the thermal structure of the middle and lower atmosphere.

Coupling from above
Not only the EPP effect but also understanding of the effects of irradiance variation of the Sun (total solar irradiance (TSI) and spectral solar irradiance (SSI)) on the atmosphere has also been progressed during the VarSITI interval. This part will be visited in the Section 3.3.3 below.

Coupling from below
During the CAWSES and CAWAES-II interval, penetration of sound waves, gravity waves, tides, and planetary waves beyond the mesopause region into the thermosphere and their interaction with the ionospheric plasma had been newly recognized (e.g., Vadas and Crowley 2010;Goncharenko and Zhang 2008;Oberheide et al. 2015). During VarSITI interval, further studies on the role of these waves in the thermosphere and ionosphere have been made. Yiğit et al. (2016) reviewed the effects of these waves on the structure of atmosphere, motivated by the IAGA/ICMA/SCOSTEP Workshop on Vertical Coupling in the Atmosphere-Ionosphere System. Figure 15 shows the height-latitude profiles of zonal winds without gravity wave forcing above 100 km (Fig. 15a) and with the forcing modeled by Yiğit et al. (2008) (Fig. 15b) and their differences (Fig. 15c). The gravity wave dissipation and their forcing to the mean zonal wind can affect the magnitude of the zonal-mean zonal wind and provide a deceleration mechanism in the mesopause region and lower thermosphere. New observations have been developed regarding connection from lower to upper atmosphere during the Var-SITI program. Gravity wave generation, propagation, and dissipation impact to the atmospheric dynamics have been reported in the special issue of "ACP (2017) sources, propagation, dissipation and impact of gravity waves" in the Atmospheric Chemistory and Phycis (https://acp. copernicus.org/articles/special_issue899.html) in relation to the German Role Of the Middle atmosphere In Climate (ROMIC) project (2013-2017). For example, Baumgarten et al. (2018) reported a 10-day continuous lidar sounding of temporal variability of tidal and gravity waves, indicating the importance of such continuous high-resolution measurements to detect interaction phenomena between waves for improving parametrization schemes of gravity waves in general circulation models. The ROMIC project has extended in 2018 for additional 3-4 years in Germany. New Program of the Antarctic Syowa MST/IS radar (PANSY) has been operational (Sato et al. 2014), and Interhemispheric Coupling Study by Observations and Modeling (ICSOM) campaigns (http://pansy.eps.s.u-tokyo.ac. jp/icsom/) have been carried out in association with Var-SITI to understand dynamical variation of the whole atmosphere during stratospheric sudden warming (SSW). Liu (2016) have shown possible influences of ENSO and QBO in the neutral density variations in the thermosphere using a 46-year-long dataset of the thermospheric density during 1967-2012.
Efforts of model extension to the thermosphere and ionosphere have also been done. Liu et al. (2018) have developed the NCAR Whole Atmosphere Community Climate Model with thermosphere and ionosphere extension (WACCM-X). This new model reproduces thermospheric composition, density, and temperature and ionospheric plasma density and ExB drifts including prereversal enhancement in the equatorial region. Pedatella and Liu (2018) noted using the WACCM-X simulations that the effects of lower atmosphere variability should be included to accurately capture smaller-scale features of the upper atmosphere response to geomagnetic storms. Shinagawa et al. (2017Shinagawa et al. ( , 2018 successfully reproduced global distribution of sporadic E layers and daily and seasonal variation of equatorial spread F, respectively, based on the Ground-to-topside model of Atmosphere and Ionosphere for Aeronomy (GAIA) model (Jin et al. 2011).

Long-term trend
The long-term trend has been extensively studied under the working group 3 of the ROSMIC project. Laštovička (2017) provided a review of recent progress in trends in the upper atmosphere, together with the special issue of the long-term changes and trends in the upper atmosphere (https://www.sciencedirect.com/journal/journalof-atmospheric-and-solar-terrestrial-physics/vol/163/ suppl/C). Figure 16 shows linear temperature trends from the 25-year data set measured by a sodium lidar at Fort Collins (41°N, 105°W) in the USA reported by She et al. (2015). By removing strong warming effect by Mt. Pinatubo eruption in 1991, they concluded a cooling trend starting from an insignificant value of 0.64 ± 0.99 K/decade at a 85-km altitude, increasing to a maximum of 2.8 ± 0.58 K/decade between 91 and 93 km, and then decreasing to a warming trend above 103 km. Berger and Lübken (2015) reported trends in the polar mesospheric clouds (PMCs) which respond to long-term changes in mesospheric temperatures at northern high latitudes for the summer in 1961-2013. They suggested that the thermal conditions near 83 km height, with a cooling of − 0.58 ± 0.32 K/decade, mainly determine trends in the PMCs, whereas cooling at lower heights, induced by stratospheric ozone, controls to a large extent the long-term behavior of PMC altitudes. Lübken et al. (2018) used model simulations over 138 years to study the impact of increasing CO 2 and H 2 O on the development of noctilucent clouds (NLCs) on centennial time scales. They concluded that NLC presumably existed centuries earlier, but the chance to observe them by the naked eye was extremely small before the twentieth century, whereas it is likely to see several NLC per season in the modern era.
The above studies of long-term trends are motivated to identify the anthropogenic effects, i.e., the effect of increasing greenhouse gases, on the mesosphere and thermosphere. On the other hand, solar cycle effects on the upper atmosphere have also been evaluated during the VarSITI interval.  Tsuchiya et al. (2019aTsuchiya et al. ( , 2019b reported clear anticorrelation between 11-year solar activity and nighttime medium-scale traveling ionospheric disturbances (MSTI Ds) at middle and subauroral latitudes. This anticorrelation can be explained by the linear growth rate of ionospheric Perkins instability, responding to the solarcycle variation of neutral-plasma collision frequency in the bottom-side ionosphere. This negative correlation of nighttime MSTIDs at middle latitudes to solar cycle is a clear contrast to the positive correlation of post-sunset plasma bubbles in the equatorial ionosphere (e.g., Nishioka et al. 2008).

Highlights of contribution from young scientists to the ROSMIC project
Many young scientists have been newly joined in the field of the upper atmosphere during the VarSITI interval and contributed to the ROSMIC project. From space to the atmosphere, Misios et al. (2015) investigate the troposphere response to the 11-year solar cycle based on the ensemble simulations of the twentieth century climate performed in the fifth phase of the Coupled Model Inter-Comparison Project (CMIP5). Wenzel et al. (2016) reported development of global ionospheric flare detection system (GIFDS) using very low frequency (VLF) radio transmissions in the northern hemisphere. Yadav et al. (2016) reported the impact of the St. Patrick's Day storm on March 15, 2017, on the evolutionary pattern of equatorial ionization anomaly over the Indian longitudes using high-resolution TEC maps. Karan et al. (2016), Karan and Pallamraju (2017) reported electrodynamic influence on the diurnal behavior of oxygen dayglow and longitudinal variations of equatorial thermospheric waves seen in the dayglow observed from the ground. rom low to upper atmosphere,  and  investigated responses of equatorial electrojet (EEJ) strength, total electron content (TEC), and oxygen dayglow in the Indian sector to the sudden stratospheric warming and solar activity. Trinh et al. (2015Trinh et al. ( , 2016Trinh et al. ( , 2018 reported a comprehensive global observation of gravity waves in the middle atmosphere and thermosphere using several low-altitude satellites, providing evidences of strong vertical coupling from lower to middle and upper atmosphere through gravity waves on global scale. Gao et al. (2015) reported double-layer structure in hydroxyl (OH) dayglow measured by the Thermosphere-Ionosphere-Mesosphere Energetics and Dynamics (TIMED) satellite. Stober et al. (2014Stober et al. ( , 2018, Stober and Chau (2015) investigated neutral density variation using meteor radar echoes in the mesosphere and developed novel technique to improve wind measurements in the mesosphere and lower thermosphere using a multi-static multi-frequency meteor radar. For the use of GNSS positioning, Macalalad et al. (2014) evaluated performance of different ionospheric models in single-frequency code-based differential GPS positioning. Seif et al. (2015Seif et al. ( , 2017Seif et al. ( , 2018 investigated daytime ionospheric scintillations at equatorial latitudes and their relationship with the sporadic E layer and the gradient-drift instability.

Implementation of the VarSITI program
The VarSITI program officially started on January 13, 2014. It was announced in a press release by the SCO-STEP's Bureau and in information for media in various countries provided by members of the SCOSTEP community. A video was prepared highlighting the scientific objectives of VarSITI (https://www.youtube.com/ watch?v=couR4MyxNPY).
A competition was announced for the VarSITI program's logo. The logo, selected by the voting during the Panel 3 of CAWSES-II Symposium in Nagoya, Japan (November 18-22, 2013), was made by Ms. Kadowaki of STEL, Nagoya University (Fig. 17).
In the following subsections, we describe various implementation activities of VarSITI. Table 1 summarizes some representative numbers of these activities.

Website
The program's website (http://varsiti.org/) contains information about VarSITI's projects and members, VarSITI-related meetings with lots of freely available presentations, journal articles and books, and links to numerous resources including real time data, data bases of space-borne and ground-based measurements, models, visualization, and forecasting tools. After the end of the VarSITI program, the website remains operational. It is maintained by Dimitar Danov from the Space Research and Technology Institute at the Bulgarian Academy of Sciences (Fig. 18).
For the duration of the program, the website had been visited more than 180 thousand times, as shown in Fig. 19.

Mailing lists
Five mailing lists were compiled and regularly updated: for each of the four projects and an integrated list for all projects. They contained the names and e-mail addresses of the scientists who have expressed their interest to participate in the program and to whom announcements were sent related to activities of VarSITI and its projects. The statistics by the end of the program is as follows: Since the start of the VarSITI program in January 2014, 60 articles, 49 highlights of young scientists, 85 meeting reports, and 24 short news from 41 countries were published in these 21 newsletter volumes.
The newsletter was distributed online and in printed form at international meetings. The full archive is available on VarSITI web site (http://varsiti.org/).

Financial support
An important tool for implementing VarSITI program was the SCOSTEP's financial support provided annually for the organization of meetings, creation of databases, and campaigns. The grants were distributed by the VarSITI co-chair following recommendations and endorsements by the Var-SITI projects' co-leaders. The grants were associated with the condition links to materials to be added to the VarSITI website and to be then freely available to all in the VarSITI community (e.g., presentations, virtual observatories, press releases, publications, etc.).
For the 5-year duration of the VarSITI program, we have organized or supported the following: 1. 64 meetings or sessions, including the three VarSITI General Symposia: VarSITI2016, VarSITI2017, and VarSITI2019 2. 16 databases 3. 1 campaign 4. 1 interdisciplinary project

VarSITI-related special issues in journals
Related to the above meetings organized for VarSITI and its projects, several topical/special issues have been published in international journals. The representative issues are listed as follows.

Database collection
The effort to collect VarSITI-related database was initiated after the SCOSTEP -World Data System (WDS) workshop "Global Data Activities for the Study of Solar-Terrestrial Variability" held during 28-30 September 2015 in Tokyo, Japan. VarSITI strives for the international collaboration in data analysis, modelling, and theory to understand how solar variability affects the Earth's environment. The provision and long-term preservation of quality-assessed data are common objectives for SCOSTEP and the WDS. The development of advanced data systems to enable scientists to perform multidisciplinary data-analysis are another common target. This workshop marked the beginning of collaborations between SCOSTEP and WDS, and a Letter of Agreement was signed by both organizations The principal objective of the workshop was to stimulate interaction among data providers (WDS members, data centers, data networks, etc.), data scientists, and data-oriented researchers of the SCOSTEP community. Data analysis of selected solar-terrestrial events was an important component of the workshop also, not only to develop the study of solar-terrestrial variability but also to establish a mutual feedback loop between "data users" and data providers.
As an outcome of the discussions during the workshop, it was recommended to have a unified metadata system that enables users to search for VarSITI data holdings from a single entry point, for example, a "Var-SITI Data Portal." After the VarSITI co-chairs asked via the VarSITI mailing list the members to provide links to their resources, a prototype of such a system has been realized in the VarSITI webpage under the title of "Var-SITI-related Database Resources" (http://varsiti.org/ → Resources → Databases). It comprises 124 different databases plus 10 data analysis and visualization tools. This database list can be used for students and young scientists for their easy access to the solar-terrestrial physics data, allowing interdisciplinary data analysis over the wide area from the Sun to the earth.

Capacity building of SCOSTEP/VarSITI
One of the major tasks of SCOSTEP is the capacity building of solar-terrestrial physics for students and young scientists. During the VarSITI interval in 2014-2018, SCOSTEP has made three large schools in collaboration with International Space Weather Initiative (ISWI) at Lima, Peru, on September 15-24, 2014, at Sangli Maharashtra, India, on November 7-17, 2016, and at Baku, Azerbaijan, in October 8-12, 2018 In addition to these SCOSTEP/ISWI schools, the Var-SITI co-chairs organized five schools in Indonesia, Nigeria, and Russia. The two schools in Nigeria were named as the International School on Equatorial and Low-Latitude Ionosphere (ISELLI) and ISELLI-2. ISELLI was held at Abuja, Nigeria, on September 14-18, 2015, and 65 students joined from 7 African countries. ISEL LI-2 was held at Ota, Nigeria, on September 11-15, 2017, and 38 students joined from 7 African countries. The two schools in Indonesia were named as the International School on Equatorial and Low-Latitude Ionosphere (ISELION) and ISELION2018. ISELION was held at Bandung, Indonesia, on March 16-20, 2015, and 39 students joined from 9 Asian countries. ISELION2018 was also held at Bandung, Indonesia, on March 5-9, 2018, and 40 students joined from 7 Asian countries. The school in Russia was held on July 9-12, 2017, in Irkutsk, Russia, in association with the VarSITI 2017 General Symposium (July 10-15, 2017) at the same venue. Thirty-five students joined in this school from 5 countries. The database list developed by VarSITI was introduced in these schools to make it available to access novel data from latest missions for young scientists and students in developing countries (Fig. 21).
6 Future topic for predictability of solar-terrestrial coupling As addressed above, there has been great progress during the 5-year VarSITI program on 2014-2018. These progresses of understanding solar-terrestrial coupling mechanisms make it possible to consider predictability of the variable coupling processes more quantitatively. The expansion of space use for better human life, such as GNSS positioning and broadcast satellites, requires more accurate prediction of space weather and hazardous events in space. In addition, threats to Earth's climate change also requests study of predictability on the effect of solar variability on Earth's climate. Under these circumstances the SCOSTEP decided to have the next program as PRESTO-Predictability of variable solarterrestrial coupling for 2020-2024 (http://www.issibj.ac. cn/Publications/Forum_Reports/201404/W0201906205 92906717714.pdf). The mission of PRESTO is to identify predictability of the variable solar-terrestrial coupling performance metrics through modeling, measurements, Shiokawa and Georgieva Progress in Earth and Planetary Science (2021)  and data analysis and to strengthen the communication between scientists and users. The details of the PRESTO definition are provided by Daglis et al. (2020). Figure 22 summarizes challenging scientific topics in the solar-terrestrial coupling processes that have been investigated during the VarSITI program, but still need further investigation during the PRESTO program. For short-term variabilities, which are shown in the leftside of Fig. 22, the following topics will stand even after the VarSITI interval.
6.1 Short-term variability from Sun to earth 6.1.1 Evolution of sunspot from solar interior to solar surface The sunspot evolution in the solar active region is one of the key processes to understand and predict the solar flare eruption (e.g., Toriumi and Takasao 2017 and references therein). It is still a difficult part to realistically model this procedure from the solar interior to solar surface and to use it for prediction of sunspot evolution.  6.1.2 Flare/CME prediction at the solar surface Many progresses and proposals of new models have been made for the flare and CME prediction during the VarSITI program, as summarized by Zhang et al. (2021) and as shown by e.g., Bamba and Kusano (2018) and Kusano et al. (2020). However, this topic is still a challenging key topic during the PRESTO interval.
6.1.3 Accurate measurements and prediction of flare UV and X-ray spectra to understand ionospheric consequences In order to understand the effect of solar flares on the ionosphere, the spectral information of the solar flare is essential, because the UV and X-ray spectral distribution provides height distribution of ionization in the ionosphere. In other words, UV and X-ray emissions at shorter (longer) wavelengths cause ionization at lower (higher) altitudes in the ionosphere. The flare spectra can be different for different flares, even if the flare class defined by the X-ray intensity is comparable. Recently, Flare Irradiance Spectral Model (FISM) (Chamberlin et al. 2008) has been widely used to model the flare effects on the ionosphere. The model provides flare UV spectra from 0.1 to 190 nm with 1 nm resolution with a time resolution of 60 s, based on the satellite measurements. However, this model does not provide real time data. Efforts to measure these spectra have been made during the VarSITI program (e.g., Woods et al. 2017). Prediction of the spectra is a further challenging issue.

Measurement and prediction of spectra of Solar Energetic Particles (SEP)
Similarly to the flare UV/X-ray spectra, the energy distribution of SEP determines the height variation of the ionization in the middle atmosphere (e.g., Turunen et al. 2009). Thus, measurement and prediction of the spectral information of SEP are an essential part to understand SEP effect on the atmosphere.

Accurate IMF-Bz modeling and prediction in the CME
During the VarSITI program, IMF-Bz modeling and prediction in the CME were the topic in a working group of ISEST/Minimax24 (see Zhang et al. 2021 in details). But it is still a challenging issue to make accurate prediction of IMF-Bz direction and intensity in the CME, which essentially controls intensity of geomagnetic disturbances in the Earth's magnetosphere. For example, Cho et al. (2017) inspected two abnormal CME-Storm events and suggested the importance of 3D magnetic field geometry of interplanetary Flux Ropes.
6.1.6 Modeling and prediction of IMF-Bz, speed, and dynamic pressure in the CIR Another solar wind structure that causes strong geomagnetic disturbances is the CIR, which comes mainly from solar coronal holes during solar quiet time (e.g., Zhang et al. 2007). Similarly to the CME, modeling and prediction of IMF-Bz, speed, and dynamic pressure in the CIR is necessary to predict geomagnetic disturbances in the Earth's magnetosphere and ionosphere.
6.1.7 Connection between the outer magnetosphere and the inner magnetosphere During the VarSITI program, the plasma dynamics in the inner magnetosphere and acceleration and dissipation of plasma through wave-particle interaction have been extensively studied by the SPeCIMEN project (Kanekal and Miyoshi 2021). The connection between the outer magnetosphere and the inner magnetosphere during storms and substorms provides boundary condition of global energy inputs from dayside magnetopause and nightside magnetotail to the inner magnetosphere. One issue is that the modeling using MHD approximation is applicable in the outer magnetosphere, while it is not in the inner magnetosphere due to significant magnetic gradient and curvature and strong coupling with the non-MHD ionosphere. This situation causes a gap of simulation codes used to describe the outer magnetosphere (MHD simulation) and inner magnetosphere (test particles and magnetosphere-ionosphere-thermosphere coupling codes). Several attempts have been made to connect these different models smoothly and self-consistently (e.g., Fok et al. 2014;Seki et al. 2018). However, the connection of these different models is still a challenging issue in the PRESTO program.
6.1.8 Connection between the magnetosphere and the ionosphere/thermosphere Another challenging issue of MHD-non-MHD modeling connection is the magnetosphere-ionosphere coupling, which is a key process for the energy input from geospace to the Earth's atmosphere (e.g., Wiltberger et al. 2017;Ebihara and Tanaka 2017). The self-consistent connection of the magnetosphere and ionosphere including neutral atmosphere is still a challenging issue in the PRESTO program.
6.1.9 Connection from the high-latitude energy input to the thermospheric dynamics and to the plasma bubble generation The plasma convection driven by MHD processes in the magnetosphere can cause associated neutral wind in the thermosphere through ion-neutral collision (e.g., Conde et al. 2001;Xu et al. 2019). The auroral particles and associated Joule heating also cause disturbances in the high-latitude thermosphere (e.g., Lu et al. 2016). The resultant neutral wind in the thermosphere further propagates to lower latitudes (e.g., Shiokawa et al. 2007) and causes westward neutral wind by the disturbance dynamo (Blanc and Richmond 1980), which may suppress ionospheric instabilities (plasma bubbles) in the equatorial evening terminator (e.g., Carter et al. 2014;Shinagawa et al. 2018). The prompt penetrating electric field (PPEF) from the magnetosphere to the equatorial ionosphere is another major force to control equatorial plasma bubbles during magnetic disturbances (storms and substorms) (e.g., Tulasi Ram et al. 2015). During the VarSITI program, these global processes have been extensively studied. However, further modeling and observation efforts should be made for accurate prediction of plasma bubbles, particularly for their day-to-day variability. The equatorial plasma bubbles are one of the most intense ionospheric disturbances on the Earth that significantly affect GNSS positioning and radiocommunication with satellites.
6.1.10 Connection from the high-latitude energy input to the composition and dynamics of the atmosphere During the VarSITI program, growing interest and significant studies have been made on the solar energetic particles (SEP) and magnetospheric energetic particles (MEP) (Ward et al. 2021 and references therein). These high-energy (keV-MeV) electrons and protons can cause significant ozone depletion in the middle atmosphere (e.g., Andersson et al. 2018) which propagates downward with a seasonal time scale. The challenging part is how these ozone depletion and associate temperature change in the middle atmosphere can change the dynamical coupling between the middle atmosphere and troposphere.

Short-term variability from the earth to geospace
From the Earth's lower atmosphere to the thermosphere and ionosphere, significant dynamical effect can be caused by penetration of acoustic sound waves (periods: less than minutes), gravity waves (several minutes to hours), tides (6 h to 1 day), and planetary waves (days). Below are some outstanding topics that will stand even after the VarSITI program.
6.2.1 Gravity wave penetration to the thermosphere and its consequences to the thermospheric dynamics and circulation During the CAWSES-II and the VarSITI program, a great progress has been made on the penetration of gravity waves into the thermosphere (e.g., Oberheide et al. 2015;Yiğit 2018). It has been well known that the gravity waves control global circulation of the middle atmosphere through their dissipation in the mesopause region. However, the role of the primary and secondary gravity waves that penetrate into the thermosphere on global circulation is still not fully understood. Several papers during the VarSITI program have suggested that the role can be significant (e.g., Miyoshi and Yigit 2019;. However, more comparison with global observation of the thermosphere and ionosphere will be needed to fully understand the gravity wave role in the thermospheric dynamics and circulation. 6.2.2 Gravity wave penetration to the thermosphere and its consequences to the ionospheric disturbances Not only the global dynamics, the gravity waves penetrated into the thermosphere can also generate ionospheric disturbances both directly and by seeding ionospheric plasma instabilities (e.g., Otsuka 2018; Moral et al. 2019;. Thus, the gravity wave in the thermosphere is one of the key parameters to understand and predict the ionospheric disturbances that can cause degradation of GNSS positioning accuracy (e.g., Tsugawa et al. 2018). Further investigation of both modeling and ground and satellite measurements should be done.

Variability of pre-reversal enhancement and generation of equatorial plasma bubbles
Pre-reversal enhancement in the evening solar terminator of the ionosphere is another key process that controls the generation of equatorial plasma bubbles. The equatorial plasma bubble is the strongest disturbance in the Earth's ionosphere and significantly affects the GNSS positioning. Recent global modeling such as by TIEGCM and GAIA successfully reproduces seasonal variability of pre-reversal enhancement, and thus, plasma bubbles (e.g., Carter et al. 2014;Shinagawa et al. 2018). However, prediction of day-to-day variability of plasma bubbles is still a challenging issue. The day-to-day variability of pre-reversal enhancement can be one of the causes of this plasma bubble variability (e.g., Ghosh et al. 2020).
6.2.4 GW/Tide/Planetary wave penetration to the thermosphere and ionosphere and its consequences in the plasmasphere As cited above, the penetration of these waves with different periods can cause significant effects on the thermosphere (neutral atmosphere) and the ionosphere (plasma). However, their effect to the magnetosphere has not well been understood. For example, storm enhanced density (SED) in the highlatitude ionosphere corresponds to a plume structure in the plasmasphere around the Earth (Foster et al. 2005). The clear conjugacy of nighttime traveling ionospheric disturbances between the northern and southern hemispheres (e.g., Otsuka et al. 2004;Narayanan et al. 2018) indicates that these ionospheric disturbances are mapped between the hemispheres along geomagnetic field line, including the equatorial plane of the near-Earth magnetosphere. The mapped density structures in the plasmasphere can contribute ducting and cross-field localization of ULF/ELF/ VLF waves in the magnetosphere that are responsible for acceleration and dissipation of high-energy (radiation belt) plasma particles. Conjugate measurements with satellites in the inner magnetosphere and ground instruments will be necessary to see these possible connections between the ionosphere and the magnetosphere.
6.3 Long-term variability 6.3.1 Solar dynamo evolution and axial dipole moment to predict variability of the 11-year solar cycle As cited in the above section on Solar Evolution and Extrema, the solar dynamo evolution to predict variability of the 11-year solar cycle is still a challenging issue even after the great efforts during the VarSITI program. The solar axial dipole moment during the solar cycle minimum can be an indicator of the activity of next solar cycle (e.g., Iijima et al. 2017;Jiang et al. 2018). The prediction of the 11-year solar cycle intensity is important for long-term prediction of the space weather condition that directly affects spacecraft operation in geospace (Bhowmik and Nandy 2018;Gopalswamy et al. 2018).
6.3.2 Long-term variation of solar total and spectral irradiance for atmospheric consequences The long-term variation of solar total and spectral irradiance (TSI and SSI) can be a key parameter to affect the Earth's climate change. Various modeling efforts have been made during the VarSITI program (e.g., Matthes et al. 2017). However, quantitative assessment of TSI and SSI effects on the atmosphere and long-term climate change is not fully obtained yet. These irradiances also cause significant variation in atmospheric density in the thermosphere where most of the low-earth orbiting satellites are flying (e.g., Liu et al. 2018). The thermospheric density variation controls the lifetime of these satellites through atmospheric drag.
6.3.3 High-energy plasma effect on the atmospheric variability, i.e., composition change and its dynamical effect via ozone concentration The energetic particle precipitation (EPP) on the atmosphere and subsequent ozone depletion have been observed and modeled in a seasonal time scale during the VarSITI program (e.g., Funke et al. 2014;Turunen et al. 2016;Andersson et al. 2018). However, their consequences on climate change are the difficult part to approach. As cited before, the key issue is the effect of the ozone depletion on the temperature variation in the middle atmosphere and their connection to the troposphere through various types of atmospheric waves.
Long-term data analysis and modeling efforts should be done further during the PRESTO program.

Solar dependence of tropospheric parameters and their connecting mechanisms
It is well known that the atmospheric density and temperature in the thermosphere significantly vary during the 11-year solar cycle. This is clearly due to the variation of the solar UV emission. We also see solar variation in the tropospheric parameters, such as a 27day solar rotation in lightning activity (Miyahara et al. 2018) and 11-year solar cycle in Pacific Walker Circulation (Misios et al. 2019). The physical mechanisms that create these solar variations in the troposphere have not been quantitatively explained. The variation of solar irradiance in the visible wavelength that reaches the troposphere seems to be too small (~0.1%) to explain the observed variations. However, the amplitude of variation in solar UV emission is much larger than that in visible wavelength, though the UV emission cannot reach the troposphere due to absorption in the thermosphere and in the ozone layer in the middle atmosphere. This reads the idea of atmospheric wave coupling to cause global circulation and tropospheric climate change (e.g., Kodera and Kuroda 2002) as below.
6.3.5 Atmospheric wave coupling to cause global circulation and tropospheric climate change The atmospheric waves, i.e., gravity waves, tides, and planetary waves, can connect the middle atmosphere and the troposphere through dynamical coupling with global circulation, resulting in climate change in the troposphere driven by middle atmosphere variation (e.g., Kodera and Kuroda 2002;Lu et al. 2017). This process is one of the key issues to understand possible connection between the solar variabilities and the Earth's climate, because the solar UV emission and energetic particles can cause significant ozone depletion and associated temperature variation in the middle atmosphere.

Separation of anthropogenic effect and solar variability effect
The global warming of the troposphere due to anthropogenic effect, such as CO 2 and methane, is more quantitatively understood than that of the solar forcing effect. However, the actual observation of climate change contains both effects in several different time scales. Modeling efforts have been made to separate these two mechanisms (e.g., Egorova et al. 2018;Solomon et al. 2019). Comparison of these modeling efforts with longterm observation becomes increasingly important to identify the mechanisms on solar forcing to climate change. Shiokawa and Georgieva Progress in Earth and Planetary Science (2021)

Consequences on human life
The red characters in Fig. 22 indicate the consequences on human life, which are (a) plasma damage to space/air vehicles and astronauts by energetic particles from the Sun and the magnetosphere; (b) degradation of radio communication and GNSS positioning by equatorial plasma instabilities (plasma bubbles), auroras, and traveling ionospheric disturbances; (c) geomagnetically induced currents (GIC) caused by auroral electrojet currents and flowing through electric supply network and pipelines; (d) satellite orbital change due to atmospheric drag that is controlled by solar UV emission and auroral plasma; and (e) climate change possibly caused by energetic particles and solar UV radiation that changes dynamics of the middle atmosphere through ozone depletion. When we consider these consequences, the parameters used in scientific research, such as the sunspot number and the Kp and Dst indices, may not be useful, and more appropriate parameters, such as temporal and spatial scales and amplitudes of ionospheric plasma variation, radiation dose and its temporal variation, and atmospheric temperature variation, will be necessary for space users and climate scientists. Close communication with space users and climate scientists will be necessary for solar-terrestrial physics for space weather and space climate applications.

Summary
To understand variability of the Sun and its consequence on Earth, communication of different fields and interdisciplinary researches are essential. During the 5-year program of VarSITI, major efforts have been made to encourage more communication between solar and interplanetary scientists and magnetosphere, ionosphere, and atmosphere scientists. The efforts were made by developing mailing lists and websites; by distributing VarSITI newsletters; by supporting interdisciplinary meetings, databases, and campaigns; and by creating database list. More than 1000 scientists in solarterrestrial physics have joined in this program.
Although it is still challenging to predict the intensity of the next solar maximum, several key procedures of solar dynamo have been understood during the VarSTII interval under the SEE project. The interplanetary connection from the solar surface to the Earth's magnetopause has been investigated in detail by combining several in situ and remote-sensing satellites and ground instruments under the ISEST/Minimax24 project. Study of the inner magnetosphere was in the golden era by the THEMIS, Van Allen Probes, and ERG satellites as well as multi-point ground measurements and the global modeling efforts that connect outer and inner magnetosphere and ionosphere. Coordination efforts have been done for these inner-magnetosphere data under the SPe-CIMEN project. The high-energy plasma input to the atmosphere and its connection to the climate change in the middle and lower atmosphere have been extensively studied, as well as short-term effects from the lower atmosphere to the thermosphere and the ionosphere, under the ROSMIC project. All these scientific achievements obtained during the VarSITI program allow direction to the quantitative assessment of predictability of short-and long-term variation of solar-terrestrial coupling system. However, there are still many challenging issues in all regions from the solar interior to the Earth's atmosphere to make prediction of their variability. Increasing the use of space by human beings, such as GNSS positioning and broadcast satellites, and increasing threat to the Earth's climate change, are resulting in increasing requirements for reliable prediction of the variability of the Sun-Earth system. Under these circumstances, it is quite natural to have the next SCOSTEP program, PRESTO, to focus on the predictability of variable solar-terrestrial coupling.