The Inﬂuence of Solar Activity on Snow Cover over the Qinghai–Tibet Plateau and Its Mechanism Analysis

: Using global ocean vertical temperature anomaly data, we identified that a significant response of the sea temperature anomaly (STA) to the solar radio flux (SRF) exists. We found that the STA exhibited a significant correlation with Asian summer and winter precipitation, among which the response from the Qinghai–Tibet Plateau (the QTP) was particularly noticeable. Based on NCEP/NCAR reanalysis data, the latent heat flux (LHF) anomaly, which plays a key role in winter precipitation in China, especially over the QTP, showed a significant response to the SRF in the Pacific. The results demonstrated the bottom-up mechanism of impact of solar activity (SA) on the plateau snow through sea–air interaction. Meanwhile, a top-down mechanism was also present. When the SRF was high, the stratospheric temperature in the low and mid-latitudes increased and the temperature gradient pointed to the pole to strengthen the westerly wind in the mid-latitudes. The EP flux showed that atmospheric long waves in the high altitudes propagated downward from the stratosphere to the troposphere. A westerly (easterly) wind anomaly occurred in the south (north) of the QTP at 500 hPa and the snowfall rate over the QTP tended to increase. When the SRF was low, the situation was the opposite, and the snowfall rate tended to decrease. The model results confirmed that when total solar irradiance (TSI) became stronger (weaker), both of the solar radiation fluxes at the top of the atmosphere and the surface temperature over the QTP increased (decreased), the vertical updraft intensified (weakened), and the snowfall rate tended to increase (decrease) accordingly. These conclusions are helpful to deepen the understanding of SA’s influence on the snow cover over the QTP.


Introduction
As the major energy driving the Earth's atmospheric circulation, solar radiation plays a key role in weather and climate [1][2][3][4][5][6][7][8][9][10][11].A well-known example of the significant impact of SA on global climate change is the cold event known as the Little Ice Age.Due to a decrease in the number of solar sunspots, the global climate became significantly colder in different regions of the world [12][13][14][15][16].
The Qinghai-Tibet Plateau snow cover is one of the physical factors that affects China's climate [17][18][19][20][21][22][23].The thermal conditions and snow cover on the QTP are sensitive to SA.When the SA intensifies (weakens), the plateau receives more (less) shortwave radiation and the temperature rises (drops); moreover, the snow cover in the winter and spring increases (decreases) [24][25][26].Although there have been some studies on the response of the snow albedo and melting rate to the different bands in the solar spectrum [27][28][29], there is still a large gap in terms of a deep understanding of the relationship between SA and snow cover variation on the QTP.Therefore, there is an urgent need for research on SA's impact on snow cover on the QTP.Based on the analysis and a significance test, Song Yan et al. [24] pointed out that on an interdecadal time scale, both the correlation of SA and winter snow depth on the QTP and the correlation of SA and the East Asian winter monsoon reached the most significant lag of 2-6 years, indicating that East Asia is the significant response area to SA. Plateau snow cover is closely related to the East Asian winter monsoon and summer monsoon.An observation showed that on the time scale of the solar 11-year cycle, SA may impact the plateau snow cover by affecting the East Asian winter monsoon [25,26].SA may also affect the correlation distribution of plateau snow cover and the following East Asian summer monsoon, resulting in abnormal summer precipitation in China.In this sense, the sun's signal is transmitted and amplified [25,26] by the plateau's snow cover.The study showed that in years of strong or weak solar activity, the correlation distribution between the plateau snow cover and precipitation in China in the following summer differed greatly and the precipitation belt in East China shifted northward (southward) in the strong (weak) year [25,27,30].This demonstrated that the modulation of SA affects the relationship between plateau snow cover and summer precipitation in China.
Solar activity also has a significant impact on the global STA and heat content [31,32], especially in the Pacific Ocean, Indian Ocean, and Atlantic Ocean.Through air-sea interaction, STA in these areas affects the global distribution of precipitation anomalies, especially over the QTP.When the STA index is high (low), the precipitation on the plateau is less (more) in summer and more (less) in winter [31].A calculation showed that there was a high correlation between winter precipitation and snow depth in the plateau (up to 99.9%).Hence, the amount of winter precipitation on the plateau reflects the amount of snow cover.The results of Zhou et al. (2021) confirmed the impact of how STA responds to SA in a precipitation anomaly, which transmits and amplifies the effect of SA.
The previous studies revealed two physical paths of solar activity's impact on climate, namely a top-down mechanism (ultraviolet irradiance) and a bottom-up mechanism (total solar irradiance) [5,9,[33][34][35][36][37][38].Although the ultraviolet irradiance can only accounts for 9% of the overall solar radiation energy, its 11-year-cycle variation range reaches as high as 6-8%, far surpassing that of the total solar irradiance energy, which can be as low as 0.1% [10], comprising 32% of the total radiation variation of solar irradiance [11].The energy variation caused by the ultraviolet irradiance affects the earth's atmosphere significantly.In particular, the ultraviolet irradiance can directly modulate the ozone distribution in the stratosphere, resulting in an increase (decrease) in the ozone amounts at low (high) latitude.The redistribution of heat energy at different latitudes leads to a change in the thermal balance in the stratosphere, and consequently to a change in atmospheric circulation.Moreover, the interaction between the stratosphere and troposphere as well as the vertical transfer anomaly of the long wave in the atmosphere affect the troposphere and cause the anomalous climate phenomena [10,11,32,39].This mechanism is called a top-down mechanism (also called a UV mechanism).Another mechanism for the impact of SA on the climate system is the bottom-up mechanism (also called the TSI mechanism), which refers to the direct influence of TSI variation on the heat distribution on the earth's surface.In particular, the ocean in the clear-sky region receives more solar radiation and causes redistribution of the sea surface temperature anomaly (SSTA), which affects oceanic evaporation and moisture content in the atmosphere through sea-air interaction and changes Hadley circulation and Walker circulation [34,38,[40][41][42][43][44].It also causes the anomalous vertical draft and determines cloud amounts in the subtropical area, resulting in the convergence of water vapor toward the precipitation in the ITCZ [9].The climate model verified that the effects of the two physical mechanisms are in phase [35,45], so the overlapping effect of the two paths amplifies SA's impact.According to numerical experiments, the contribution of the top-down mechanism is seemingly higher than that of the bottom-up mechanism [35].In recent years, more studies proved the two mechanisms of top-down and bottom-up paths in the impact of SA on East Asian climate and the tropical Pacific.For example, Zhao et al. [46] pointed out that the dynamical responses of the tropical summer monsoon in the low levels and the upper subtropical westerly jet to the 11-year solar cycle transmitted both bottom-up and top-down solar signals, respectively, and that the common actions of the monsoon and the jet likely amplified the solar signals at the northern boundary of the monsoon to some extent.Huo et al. (2021) also revealed that a distinct 11-year solar signal in the tropical Pacific SSTA was found through model simulation and observation.Moreover, it was independent of the ENSO cycle; however, both of the patterns were similar to some extent.This indicated that the atmospheric forcing caused by net shortwave radiation is the major contributor to the initial warming in the solar maximum years, but this warming response is amplified by the air-sea interaction during the lagged years (bottom-up) [47].Meanwhile, cloud cover patterns provide the positive feedback in the solar maximum years to amplify the initial solar signal but have a negative contribution during the lagged years, which also suggests a "bottom-up" influence of SA.
Previous studies indicated that SA affected the plateau snow cover and East Asian Monsoon significantly.Moreover, SA influenced the tropical sea-air system through the top-down and bottom-up mechanisms together.However, the physical paths through which SA affects precipitation and snow cover on the plateau are still not clear.Hence, based on the previous studies, this study aimed to reveal the physical phenomena and explore the possible paths through which SA affects the plateau snow cover.
This study was organized as follows.In Section 2, the data and methods used in this study are introduced briefly, including the model simulation.The main analysis of the mechanisms of SA's impact on the plateau snow cover is presented in Section 3. A description of numerical experiments on the response of plateau snow to SA is included in Section 3 as well.Finally, conclusions and discussions are provided in Section 4.

Data
The SRF data used in this paper were the 0.7 cm (2800 MHz) wavelength data for the sun from the data center of the National Oceanic and Atmospheric Administration (NOAA) from 1947 to 2018 (http://www.esrl.noaa.gov/psd/data/correlation/solar.data;unit: sfu; 1 sfu = 10 −22 W•m −2 •Hz −1 ).SRF data were used to depict the solar activity.The atmo- spheric circulation data were from the global reanalysis monthly data of NCEP/NCAR [48].The precipitation data from China were the rain gauge data from 160 stations of the National Meteorological Information Center, China Meteorological Administration.The snow cover data on the QTP were the daily snow depth data from National Meteorological Information Center, China Meteorological Administration; these data were sorted into the continuous time series from 51 stations .The study area of this paper (the QTP) is shown in Figure 1.For the detailed method of processing the plateau snow cover data, please refer to Song et al. (2011) [49].The monthly average water-equivalent snow depth grid data (1951-2017) were from the second edition of the 20th Century Reanalysis data set (20 CRv2) (the grid data represent the mass of melted water from the snow cover per unit area; unit: kg•m −2 ; resolution: 2 • × 2 • ).Compared with the observation data, the grid data were more applicable for depicting the spatial pattern of snow anomalies on the QTP and conducting an empirical orthogonal decomposition (EOF) analysis because the observation data for the QTP were not evenly distributed, with more data from the east, where it is densely populated, and less from the west, where it is sparsely populated.
The STA data were from the WOD09 data set (World Ocean Atlas, NCAR Climate Data Guide; 1955-2016) of the National Oceanic Data Center (NODC) at a grid resolution of 1° × 1°.The STA data for the 0~700 m vertical depth were considered to be 5-year running mean data to filter out the effects of ENSO.The global precipitation data were the CMAP monthly precipitation rate data (CPC merged analysis of precipitation) of the NOAA (January 1979 to September 2017) at a grid resolution of 2.5 × 2.5.

Power Spectrum Method
Power spectrum analysis is widely used as a method for frequency domain analysis based on a Fourier transform.In this way, the total energy of one time series is decomposed into components at different frequencies and then the variance contribution of different frequency waves is used to identify the main cycle of the sequence to determine the main frequency of the cycle; namely, the significant period implied in the sequence [31,50].In this paper, this method was used to detect the SA signal from STA caused by SA to determine whether the STA was of an 11-year periodic cycle.

Composite Mean Difference (CMD) Method
The CMD method is usually used to deduce the spatial distribution of climate element fields in response to SA [51].In this paper, the CMD method was used to distinguish the spatial distribution of the climate element difference field of high SA years (HSAYs) from low SA years (LSAYs).The normalized SRF data values larger (lower) than zero were ascribed to high (low) SA years (called HSAYs (LSAYs)) to obtain two groups of climate element data corresponding to HSAYs and LSAYs, respectively.The spatial distribution of climate elements in response to SRF was generated by calculating the difference between the two groups of data.

Eliassen-Palm (EP) Flux
EP flux, which was pointed out by Eliassen and Palm [52] in 1961, can represent the propagation of a planetary wave in the atmosphere.It can describe the propagation direction and size of atmospheric long waves along meridional and vertical directions.It can The STA data were from the WOD09 data set (World Ocean Atlas, NCAR Climate Data Guide; 1955-2016) of the National Oceanic Data Center (NODC) at a grid resolution of 1 • × 1 • .The STA data for the 0~700 m vertical depth were considered to be 5-year running mean data to filter out the effects of ENSO.The global precipitation data were the CMAP monthly precipitation rate data (CPC merged analysis of precipitation) of the NOAA (January 1979 to September 2017) at a grid resolution of 2.5 × 2.5.

Methodology 2.2.1. Power Spectrum Method
Power spectrum analysis is widely used as a method for frequency domain analysis based on a Fourier transform.In this way, the total energy of one time series is decomposed into components at different frequencies and then the variance contribution of different frequency waves is used to identify the main cycle of the sequence to determine the main frequency of the cycle; namely, the significant period implied in the sequence [31,50].In this paper, this method was used to detect the SA signal from STA caused by SA to determine whether the STA was of an 11-year periodic cycle.

Composite Mean Difference (CMD) Method
The CMD method is usually used to deduce the spatial distribution of climate element fields in response to SA [51].In this paper, the CMD method was used to distinguish the spatial distribution of the climate element difference field of high SA years (HSAYs) from low SA years (LSAYs).The normalized SRF data values larger (lower) than zero were ascribed to high (low) SA years (called HSAYs (LSAYs)) to obtain two groups of climate element data corresponding to HSAYs and LSAYs, respectively.The spatial distribution of climate elements in response to SRF was generated by calculating the difference between the two groups of data.

Eliassen-Palm (EP) Flux
EP flux, which was pointed out by Eliassen and Palm [52] in 1961, can represent the propagation of a planetary wave in the atmosphere.It can describe the propagation direction and size of atmospheric long waves along meridional and vertical directions.It can also express the size and direction of the eddy motion heat flux and momentum flux.The specific formula is as follows: The divergence of EP flux can also reflect the wave-flow interaction.When the divergence of the EP flux is larger (smaller) than zero, the average speed of the zonal wind will increase (decrease).In addition, EP flux divergence can also diagnose the meridional flux of potential vorticity and its distribution inhomogeneity with latitude can reflect the change of potential vorticity with time.

Model Simulation
Using the Earth system model from the National Center for Atmospheric Research, USA (NCAR CESM2.1.1),fully coupled numerical integration experiments, including the annual variation in TSI from 1979 to 2009, were carried out to analyze the variation in the climate system within the 30 years of time.The atmospheric horizontal resolution was 1.9 • × 2.5 • .The horizontal resolution of ocean mode and sea ice mode was a nominal 1 • × Based on the TSI variation between 1979 and 2009, the years with TSI values higher than a 1.0 standard deviation were regarded as high TSI years: 1979TSI years: , 1980TSI years: , 1981TSI years: , 1989TSI years: , 1990TSI years: , 2000TSI years: , 2001TSI years: , and 2002, altogether eight years; and the years with TSI values lower than a −1.0 standard deviation were categorized as low TSI years: 1985TSI years: , 1986TSI years: , 1995TSI years: , 1996TSI years: , 1997TSI years: , 2007TSI years: , 2008TSI years: , and 2009, also eight years in total.

Analysis of the Bottom-Up Mechanism of the Impact of SRF on Plateau Snow Cover
Song et al. (2016,2022) [24,30] pointed out that the impact of SRF on plateau snow cover is especially significant in peak and valley years of SRF.The SRF time series exhibits an obvious 11-year cycle with alternating peak and valley SA years (Figure 2).When considering the impact of SA on snow cover on the QTP, the effects of ENSO and volcano eruption must be excluded because they are strong interannual signals that cause climate anomalies and ENSO has a distinct interdecadal variation, which could induce anomalous snow cover on the QTP [54].Therefore, the peak years of SRF selected, excluding the impact of strong ENSO events (i.e., strong El Niño and La Niña events) and volcano eruption, were 1957/1958, 1967/1968, 1979/1980, 1990/1991, and 2001/2002; while the valley years were 1963/1964, 1975/1976, 1986/1987, 1995/1996, and 2008/2009, respectively.Figure 3a shows that the peak (valley) years of SRF were accompanied by higher (lower) winter precipitation over the plateau; the black dots represent a score of more than 0.05 in the significance test.The significant precipitation anomaly areas covered Northeast China, the Hetao area, and the QTP (the circled area), with the correlation coefficients distributed in "+-+" mode.According to our calculations, the correlation coefficients of winter precipitation and snow amount on the QTP were significantly correlated.This implied that the plateau winter precipitation could reflect the snow cover on the plateau.Figures 3b and 3c represent the differences in the average snow days at the gauge station in HSAYs and LSAYs in winter and in spring, respectively.It is obvious that the snow amount in HSAYs was more than in LSAYs over most area of the QTP whether in winter or spring.Based on our calculation results, we discovered that the average snow days in winters of HSAYs surpassed those of LSAYs by 50 days per year, whereas the average snow days in springs of HSAYs exceeded those of LSAYs by 57 days per year.

Analysis of the Bottom-Up Mechanism of the Impact of SRF on Plateau Snow Cover
Song et al. (2016,2022) [24,30] pointed out that the impact of SRF on plateau snow cover is especially significant in peak and valley years of SRF.The SRF time series exhibits an obvious 11-year cycle with alternating peak and valley SA years (Figure 2).When considering the impact of SA on snow cover on the QTP, the effects of ENSO and volcano eruption must be excluded because they are strong interannual signals that cause climate anomalies and ENSO has a distinct interdecadal variation, which could induce anomalous snow cover on the QTP [54].Therefore, the peak years of SRF selected, excluding the impact of strong ENSO events (i.e., strong El Niño and La Niña events) and volcano eruption, were 1957/1958, 1967/1968, 1979/1980, 1990/1991, and 2001/2002; while the valley years were 1963/1964, 1975/1976, 1986/1987, 1995/1996, and 2008/2009, respectively.Figure 3a shows that the peak (valley) years of SRF were accompanied by higher (lower) winter precipitation over the plateau; the black dots represent a score of more than 0.05 in the significance test.The significant precipitation anomaly areas covered Northeast China, the Hetao area, and the QTP (the circled area), with the correlation coefficients distributed in "+-+" mode.According to our calculations, the correlation coefficients of winter precipitation and snow amount on the QTP were significantly correlated.This implied that the plateau winter precipitation could reflect the snow cover on the plateau.Figure 3b and 3c represent the differences in the average snow days at the gauge station in HSAYs and LSAYs in winter and in spring, respectively.It is obvious that the snow amount in HSAYs was more than in LSAYs over most area of the QTP whether in winter or spring.Based on our calculation results, we discovered that the average snow days in winters of HSAYs surpassed those of LSAYs by 50 days per year, whereas the average snow days in springs of HSAYs exceeded those of LSAYs by 57 days per year.This study analyzed the spatial distribution of the synchronous and lag correlation between winter snow cover and the SRF (Figure 4).We found that there was a significant synchronous and 1-year lag correlation (it exceeded the 0.05 significant level) between the plateau snow and SA signals by calculating the correlation coefficients (omitted).A positive correlation was identified between snow over most area of the plateau and the SRF, while a local negative correlation only existed on the western edge and in the northeast and southeast of the QTP.The significantly correlated areas were mainly in the middle of the QTP along the Tanggula Mountains/Nianqing Tanggula Mountains/Hengduan Mountains and the eastern section of Bayankala Mountains and Animaqing Mountains.With the lag time increasing, the extent of the positively correlated regions gradually decreased and that of the negative correlation regions gradually increased.For the same amount of incident solar irradiance, there is increasing absorption of incident solar radiation on the cloudless underlying surface than that with clouds, which will cause an abnormal heat exchange between the underlying surface and the atmosphere.This can be reflected as the anomalous sensible heat flux (SHF) and latent heat flux (LHF) between the ocean and atmosphere.Figure 5 presents the difference field of the LHF between the peak and valley years of the SRF.It shows four areas with a significantly anomalous LHF in the Pacific.The detailed extents of these four areas (A, B, C, and D) are displayed in Table 2.The index Ilatent was defined as the sum of the values of anomalous LHFs in Zone A, Zone B, Zone C, and Zone D; i.e., Ilatent = LB + LC + LD − LA. Figure 6 presents the spatial distribution of the correlation between the LHF and winter precipitation in China.It shows that the positive plateau precipitation anomaly was very similar to that given in Figure 3a, indicating that the LHF in the central and western tropical Pacific responding to the SRF played a key role in winter precipitation anomalies over the QTP.The plateau snow cover anomaly affected precipitation in China in the following summer, transmitting the energy of SA from winter to summer.Comparatively speaking, the correlation between the SHF anomaly's response to the SRF and winter precipitation was relatively weaker than that of the LHF.This showed that SHF played less of a role in precipitation.The reason was most likely that precipitation is more closely related to the LHF than the SHF, which includes the water cycle.Jin et al. (2017) [55] pointed out that there was a significant correlation between the geomagnetic index and the global LHF that leads to global climate change by influencing the release of Joule energy from the inner earth and air-sea interaction.The connection of the geomagnetic index and SA is relatively For the same amount of incident solar irradiance, there is increasing absorption of incident solar radiation on the cloudless underlying surface than that with clouds, which will cause an abnormal heat exchange between the underlying surface and the atmosphere.This can be reflected as the anomalous sensible heat flux (SHF) and latent heat flux (LHF) between the ocean and atmosphere.Figure 5 presents the difference field of the LHF between the peak and valley years of the SRF.It shows four areas with a significantly anomalous LHF in the Pacific.The detailed extents of these four areas (A, B, C, and D) are displayed in Table 2.The index I latent was defined as the sum of the values of anomalous LHFs in Zone A, Zone B, Zone C, and Zone D; i.e., I latent = L B + L C + L D − L A . Figure 6 presents the spatial distribution of the correlation between the LHF and winter precipitation in China.It shows that the positive plateau precipitation anomaly was very similar to that given in Figure 3a, indicating that the LHF in the central and western tropical Pacific responding to the SRF played a key role in winter precipitation anomalies over the QTP.The plateau snow cover anomaly affected precipitation in China in the following summer, transmitting the energy of SA from winter to summer.Comparatively speaking, the correlation between the SHF anomaly's response to the SRF and winter precipitation was relatively weaker than that of the LHF.This showed that SHF played less of a role in precipitation.The reason was most likely that precipitation is more closely related to the LHF than the SHF, which includes the water cycle.Jin et al. (2017) [55] pointed out that there was a significant correlation between the geomagnetic index and the global LHF that leads to global climate change by influencing the release of Joule energy from the inner earth and air-sea interaction.The connection of the geomagnetic index and SA is relatively close.Therefore, it is also an example of the indirect influence of SA on global climate change by causing anomalous LHF [55].
cipitation.The reason was most likely that precipitation is more closely related to the LHF than the SHF, which includes the water cycle.Jin et al. (2017) [55] pointed out that there was a significant correlation between the geomagnetic index and the global LHF that leads to global climate change by influencing the release of Joule energy from the inner earth and air-sea interaction.The connection of the geomagnetic index and SA is relatively close.Therefore, it is also an example of the indirect influence of SA on global climate change by causing anomalous LHF [55].Table 2.The extent of anomalous latent heat flux regions.An LHF anomaly caused by SRF during the sun's 11-year cycle can also affect summer precipitation in the following summer season of China, hence extending and amplifying the impact of solar role in the climate from winter to the following summer.Meanwhile, by impacting the global vertical ocean temperature, SA could affect winter and summer precipitation over the QTP indirectly.This also amplifies the impact of SA.Here, the WOD09 data set (World Ocean Atlas, NCAR Climate Data Guide, 1955-2016) (grid point resolution: 1 • × 1 • ) of STA data from the National Oceanographic Data Center (NODC, USA) was utilized to determine the 5-year running mean of STA data at each level of the vertical range from 0 m to 700 m.In this way, ENSO signals were eliminated.The calculation results showed that the STA at 0-200 m responded to the SRF more significantly than at 300-700 m.Using the accumulated STA at 0-200 m, we analyzed the composite difference in STA response to the SRF between high and low years, as shown in Figure 7. Six regions with a significant response were generated based on the calculations; namely, the A-F regions, all of which responded significantly to SA.The detailed extent and significant period are exhibited in Table 3.According to the power spectrum analysis, six regions (A, B, C, D, E, and F) among the global significant response regions were detected to contain the similar significant period of the SRF; that is, a period of around 11.3 years.In addition, all six of the regions passed the 95% reliability test of red noise.The temporal-meridional profile of STA further verified that the six regions indeed exhibited an 11-year cycle of STA significantly, which we thought to be the solar signal [31].Comparatively, the six regions we obtained as shown in Figure 7 were in accordance with those shown in Figure 3 in Wang (2015) [56], including two high-responding areas over the Pacific Ocean.

Zone Range Mark
We obtained the average of the absolute values of STA in these six areas and defined the index of STA as Z = STB + STF-STA-STC-STD-STE.The correlation coefficient of the index and the SRF reached 0.473, thus passing the 0.001 significance test.This indicated that the index Z defined by the STA's response to the SRF could reflect the intensity of the SA.
The synchronous correlation between the STA index Z and winter (summer) precipitation in East Asia showed that there was a significant precipitation anomaly both in winter and in summer over the QTP (Figure 8).A higher (lower) index Z was accompanied by higher (lower) winter precipitation and lower (higher) summer precipitation over the  According to the power spectrum analysis, six regions (A, B, C, D, E, and F) among the global significant response regions were detected to contain the similar significant period of the SRF; that is, a period of around 11.3 years.In addition, all six of the regions passed the 95% reliability test of red noise.The temporal-meridional profile of STA further verified that the six regions indeed exhibited an 11-year cycle of STA significantly, which we thought to be the solar signal [31].Comparatively, the six regions we obtained as shown in Figure 7 were in accordance with those shown in Figure 3 in Wang (2015) [56], including two high-responding areas over the Pacific Ocean.
We obtained the average of the absolute values of STA in these six areas and defined the index of STA as Z = ST B + ST F −ST A −ST C −ST D −ST E .The correlation coefficient of the index and the SRF reached 0.473, thus passing the 0.001 significance test.This indicated that the index Z defined by the STA's response to the SRF could reflect the intensity of the SA.
The synchronous correlation between the STA index Z and winter (summer) precipitation in East Asia showed that there was a significant precipitation anomaly both in winter and in summer over the QTP (Figure 8).A higher (lower) index Z was accompanied by higher (lower) winter precipitation and lower (higher) summer precipitation over the QTP.This indicated that the STA could reflect the SA's signals and affected the plateau's precipitation, hence amplifying the impact of SA.
ther verified that the six regions indeed exhibited an 11-year cycle of STA significantly, which we thought to be the solar signal [31].Comparatively, the six regions we obtained as shown in Figure 7 were in accordance with those shown in Figure 3 in Wang (2015) [56], including two high-responding areas over the Pacific Ocean.
We obtained the average of the absolute values of STA in these six areas and defined the index of STA as Z = STB + STF-STA-STC-STD-STE.The correlation coefficient of the index and the SRF reached 0.473, thus passing the 0.001 significance test.This indicated that the index Z defined by the STA's response to the SRF could reflect the intensity of the SA.
The synchronous correlation between the STA index Z and winter (summer) precipitation in East Asia showed that there was a significant precipitation anomaly both in winter and in summer over the QTP (Figure 8).A higher (lower) index Z was accompanied by higher (lower) winter precipitation and lower (higher) summer precipitation over the QTP.This indicated that the STA could reflect the SA's signals and affected the plateau's precipitation, hence amplifying the impact of SA.Shading indicates the anomaly regions that reached a 0.1 significance level; the contour interval is 0.5%.(Units: latitude is ° N , and longitude is ° E.)

Analysis of the Top-Down Mechanism of the Impact of SRF on Plateau Snow Cover
The discussion in the previous subsection regarded the bottom-up mechanism of SA's impact on climate.SRF affects STA and causes ocean LHF anomalies, which further lead to winter and summer precipitation anomalies in East Asia.Consequently, winter precipitation increases significantly over the QTP in high years of SRF, resulting in higher

Analysis of the Top-Down Mechanism of the Impact of SRF on Plateau Snow Cover
The discussion in the previous subsection regarded the bottom-up mechanism of SA's impact on climate.SRF affects STA and causes ocean LHF anomalies, which further lead to winter and summer precipitation anomalies in East Asia.Consequently, winter precipitation increases significantly over the QTP in high years of SRF, resulting in higher amounts of plateau snow cover.This subsection begins to discuss the top-down mechanism.Figure 9 consists of two diagrams: the vertical profile of the winter temperature difference (Figure 9a) and the geopotential height difference (Figure 9b) with a zonal mean between the high and low years of SRF.The profiles show that during peak years of SA (from 400 hPa to 100 hPa of the northern hemisphere), the Arctic cools and the low and middle latitudes become warmer, which corresponds to the temperature gradient pointing to the pole and is conducive to the positive phase of the Arctic Oscillation and the strengthening of the westerly wind in the middle latitudes.Such a distribution of temperature anomalies is in accordance with Figure 1 in the paper by Frame and Gray (2010) [57].The mean temperature difference profile between the equator and the pole shows that the size of temperature gradient in the lower stratosphere and the upper troposphere is greater than zero, which is conducive to the strengthening of the westerly in the northern hemisphere (figure omitted).Such a temperature gradient can last until the spring, summer, or autumn of the next year and can also propagate downward to the lower troposphere in the spring and summer.The geopotential height anomaly profile demonstrated that in the troposphere and the middle and lower stratosphere of the Northern Hemisphere, negative anomalies were identified in the polar regions while positive anomalies were present in the middle and low latitudes.The pressure gradient force pointed from the equator to the pole, which was also conducive to the strengthening of the westerly.During the valley years of SA, the situation reversed.Figure 10 presents the profile of the zonal wind difference in winter between high and low SRF years.It shows that westerlies prevailed in the atmosphere from the stratosphere to the troposphere at 40-60° N of the mid-latitudes, easterlies prevailed at 25-40° N, and an east wind prevailed in the troposphere and upper middle stratosphere north of 65° N.These phenomena indicated that in the middle latitudes of the stratosphere, the zonal westerly wind can propagate into the troposphere, strengthening the westerlies in the mid-latitudes of the troposphere and enhancing the annular mode in the Northern Hemisphere.Thus, westerly winds prevailed in the tropical stratosphere during HSAYs and easterly winds prevailed during LSAYs.How does the westerly circulation anomaly of the stratosphere transfer down to the troposphere?One path is the downward vertical propagation of atmospheric long waves triggered by thermal anomalies in the stratosphere.Figure 11 presents the EP flux and its divergence in high and low years of SRF, respectively.It is clearly shown that long waves in the stratosphere were transferred to the troposphere through vertical propagation of  Figure 10 presents the profile of the zonal wind difference in winter between high and low SRF years.It shows that westerlies prevailed in the atmosphere from the stratosphere to the troposphere at 40-60° N of the mid-latitudes, easterlies prevailed at 25-40° N, and an east wind prevailed in the troposphere and upper middle stratosphere north of 65° N.These phenomena indicated that in the middle latitudes of the stratosphere, the zonal westerly wind can propagate into the troposphere, strengthening the westerlies in the mid-latitudes of the troposphere and enhancing the annular mode in the Northern Hemisphere.Thus, westerly winds prevailed in the tropical stratosphere during HSAYs and easterly winds prevailed during LSAYs.How does the westerly circulation anomaly of the stratosphere transfer down to the troposphere?One path is the downward vertical propagation of atmospheric long waves triggered by thermal anomalies in the stratosphere.Figure 11 presents the EP flux and its divergence in high and low years of SRF, respectively.It is clearly shown that long waves in the stratosphere were transferred to the troposphere through vertical propagation of How does the westerly circulation anomaly of the stratosphere transfer down to the troposphere?One path is the downward vertical propagation of atmospheric long waves triggered by thermal anomalies in the stratosphere.Figure 11 presents the EP flux and its divergence in high and low years of SRF, respectively.It is clearly shown that long waves in the stratosphere were transferred to the troposphere through vertical propagation of long waves during HSAYs.Between 60 • N and 80 • N at high latitudes, there were consistent downward flows of EP flux from the stratosphere to the troposphere, indicating the downward propagation of long waves.Consequently, the impact of SA was transferred downward into the troposphere.In the meridional direction, the sensible heat flux and momentum flux pointed from the pole to the lower latitude and the heat and momentum were transported to the middle latitudes, which was conducive to the warming and momentum accumulation in the middle latitudes.In the stratosphere, the EP flux arrows point from the pole to the equator, indicating a northward transport of potential vorticity.From 45 • N to 70 • N of the stratosphere, there was a strong divergence in the EP flux and the zonal westerly wind there was strengthened.On the contrary, in the southern and polar regions, there was a convergence of the EP flux, the zonal westerly wind weakened, and the easterly wind strengthened.Such a zonal wind distribution with latitude corresponded to the positive phase of the annular mode; that is, the positive phase of the Arctic Oscillation (AO).The above situation reversed during LSAYs.Composites of the EP flux vector in HSAYs and LSAYs were in accordance with Figure 1 in the paper by Matthes (2006) [39] between the solar max and min numerical experiments.
Atmosphere 2022, 13, x FOR PEER REVIEW 14 of 21 long waves during HSAYs.Between 60° N and 80° N at high latitudes, there were consistent downward flows of EP flux from the stratosphere to the troposphere, indicating the downward propagation of long waves.Consequently, the impact of SA was transferred downward into the troposphere.In the meridional direction, the sensible heat flux and momentum flux pointed from the pole to the lower latitude and the heat and momentum were transported to the middle latitudes, which was conducive to the warming and momentum accumulation in the middle latitudes.In the stratosphere, the EP flux arrows point from the pole to the equator, indicating a northward transport of potential vorticity.From 45° N to 70° N of the stratosphere, there was a strong divergence in the EP flux and the zonal westerly wind there was strengthened.On the contrary, in the southern and polar regions, there was a convergence of the EP flux, the zonal westerly wind weakened, and the easterly wind strengthened.Such a zonal wind distribution with latitude corresponded to the positive phase of the annular mode; that is, the positive phase of the Arctic Oscillation (AO).The above situation reversed during LSAYs.Composites of the EP flux vector in HSAYs and LSAYs were in accordance with Figure 1 in the paper by Matthes (2006) [39] between the solar max and min numerical experiments.The monthly EP flux showed that during high SRF years, SA signals were transmitted downward from November to December, causing a downward transfer of long waves.Accordingly, the westerly wind was strengthened in high-latitude area, accompanied by the barotropic structure and the positive phase of AO.The barotropic structure and the positive phase lasted until January of the following year.During low-SRF years, the wave vertical transmission signal weakened in high-latitude areas and long waves were transferred upward.The upward transmission reached its highest level in January, resulting in the accumulation of heat in the Arctic region.This situation lasted until February, and was conducive to the negative phase of the AO.The barotropic structure was identified from January to February.The figures above are omitted here.The monthly EP flux showed that during high SRF years, SA signals were transmitted downward from November to December, causing a downward transfer of long waves.Accordingly, the westerly wind was strengthened in high-latitude area, accompanied by the barotropic structure and the positive phase of AO.The barotropic structure and the positive phase lasted until January of the following year.During low-SRF years, the wave vertical transmission signal weakened in high-latitude areas and long waves were transferred upward.The upward transmission reached its highest level in January, resulting in the accumulation of heat in the Arctic region.This situation lasted until February, and was conducive to the negative phase of the AO.The barotropic structure was identified from January to February.The figures above are omitted here.
The interaction between the stratosphere and troposphere alters the distribution of atmospheric circulation field in the troposphere.It is the difference at 500 hPa zonal wind field between HSAYs and LSAYs (Figure 12a) that causes not only a significant easterly wind anomaly over the QTP, but also a strong westerly wind anomaly to the south of the QTP around the Indian Peninsula, the Bay of Bengal, and the Arabian Sea.Such an anomalous zonal wind shear with latitude can produce positive vorticity anomalies, which are favorable for the generation of anomalous cyclones there, and hence favorable for the rich precipitation over the QTP.In this way, the plateau snow cover increased accordingly, as shown in Figure 2b.Based on our calculations, we discovered that the average snow days in winters of HSAYs surpassed those of LSAYs by 50 days per year.The anomalous zonal wind caused by SA can last until spring continuously (Figure 12b), although the wind shear becomes weak.In particular, a significant westerly wind anomaly existed in the Bay of Bengal to the south of the QTP in spring.It was also favorable for the increase in spring snow on the QTP (Figure 2c).In fact, the average snow days in springs of HSAYs exceeded those of LSAYs by 57 days per year.
Atmosphere 2022, 13, x FOR PEER REVIEW 15 of 21 The interaction between the stratosphere and troposphere alters the distribution of atmospheric circulation field in the troposphere.It is the difference at 500 hPa zonal wind field between HSAYs and LSAYs (Figure 12a) that causes not only a significant easterly wind anomaly over the QTP, but also a strong westerly wind anomaly to the south of the QTP around the Indian Peninsula, the Bay of Bengal, and the Arabian Sea.Such an anomalous zonal wind shear with latitude can produce positive vorticity anomalies, which are favorable for the generation of anomalous cyclones there, and hence favorable for the rich precipitation over the QTP.In this way, the plateau snow cover increased accordingly, as shown in Figure 2b.Based on our calculations, we discovered that the average snow days in winters of HSAYs surpassed those of LSAYs by 50 days per year.The anomalous zonal wind caused by SA can last until spring continuously (Figure 12b), although the wind shear becomes weak.In particular, a significant westerly wind anomaly existed in the Bay of Bengal to the south of the QTP in spring.It was also favorable for the increase in spring snow on the QTP (Figure 2c).In fact, the average snow days in springs of HSAYs exceeded those of LSAYs by 57 days per year.The above analysis explained the top-down mechanism (ultraviolet irradiance mechanism) relevant to the impact of SA on plateau snow cover.Together with the analysis of the bottom-up mechanism (total solar irradiance mechanism), this paper proved that SA affects plateau snow cover through two physical paths (bottom-up and top-down mechanisms).Thus, the physical reason for the response of plateau snow to SA has been deeply explored.

Numerical Experiments on the Response of Plateau Snow to SA
By using the Climate Earth System Model (CESM) of the USA (NCAR CESM2.1.1),this paper studied the impact of TSI (total solar irradiance), including its annual variation, on plateau snow and the atmosphere circulation.
The simulation results showed a significant response of plateau snow to SA.Based on the high and low years of SA selected by using TSI standardized data in the model, the spatial distributions of the annual temperature at the reference height (TREHFT) (2 m air temperature) as well as the annual net solar flux at the top of the atmosphere (FSNT) in HSAY and LSAY are given in Figure 13.We found that on the QTP in HSAYs, the solar flux at the top of the atmosphere intensified and that the near-surface air temperature also The above analysis explained the top-down mechanism (ultraviolet irradiance mechanism) relevant to the impact of SA on plateau snow cover.Together with the analysis of the bottom-up mechanism (total solar irradiance mechanism), this paper proved that SA affects plateau snow cover through two physical paths (bottom-up and top-down mechanisms).Thus, the physical reason for the response of plateau snow to SA has been deeply explored.

Numerical Experiments on the Response of Plateau Snow to SA
By using the Climate Earth System Model (CESM) of the USA (NCAR CESM2.1.1),this paper studied the impact of TSI (total solar irradiance), including its annual variation, on plateau snow and the atmosphere circulation.
The simulation results showed a significant response of plateau snow to SA.Based on the high and low years of SA selected by using TSI standardized data in the model, the spatial distributions of the annual temperature at the reference height (TREHFT) (2 m air temperature) as well as the annual net solar flux at the top of the atmosphere (FSNT) in HSAY and LSAY are given in Figure 13.We found that on the QTP in HSAYs, the solar flux at the top of the atmosphere intensified and that the near-surface air temperature also increased significantly; both of these results passed the 0.05 confidence level significance test.The findings indicated that with a strong SA, both the heat flux at the top of the atmosphere and the near-surface air temperature over the QTP would increase, and vice versa.increased significantly; both of these results passed the 0.05 confidence level significance test.The findings indicated that with a strong SA, both the heat flux at the top of the atmosphere and the near-surface air temperature over the QTP would increase, and vice versa.The atmosphere circulation field showed that in HSAYs, there was a low-pressure anomaly (Figure 14a) and a positive vorticity shear on the plateau at 500 hPa (omitted), while a high-pressure anomaly (Figure 14b) and a negative vorticity shear were identified at 100 hPa (omitted), which was consistent with the observation given above.Furthermore, an anomalous ascending movement was found on the plateau area (Figure 14c), resulting in an excessive snowfall rate there (Figure 14d).The above circulation anomalies were conducive to increased snow on the QTP in HSAYs, while the opposite occurred in LSAYs.These results also were accordance with those of Song Yan et al. [24,25]; i.e., when SA became strong, the winter and spring snow increased over the QTP.The atmosphere circulation field showed that in HSAYs, there was a low-pressure anomaly (Figure 14a) and a positive vorticity shear on the plateau at 500 hPa (omitted), while a high-pressure anomaly (Figure 14b) and a negative vorticity shear were identified at 100 hPa (omitted), which was consistent with the observation given above.Furthermore, an anomalous ascending movement was found on the plateau area (Figure 14c), resulting in an excessive snowfall rate there (Figure 14d).The above circulation anomalies were conducive to increased snow on the QTP in HSAYs, while the opposite occurred in LSAYs.These results also were accordance with those of Song Yan et al. [24,25]; i.e., when SA became strong, the winter and spring snow increased over the QTP.

Conclusions
Using observation and model simulation methods, this study explained the possible mechanisms of SA's influence on the plateau snow cover using two aspects; namely, the bottom-up mechanism (the impact of ocean thermal anomalies on the atmosphere) and the top-down mechanism (the impact of the stratospheric atmosphere on tropospheric atmosphere).Moreover, with CESM, the response of the plateau snowfall rate to SA was

Conclusions
Using observation and model simulation methods, this study explained the possible mechanisms of SA's influence on the plateau snow cover using two aspects; namely, the bottom-up mechanism (the impact of ocean thermal anomalies on the atmosphere) and the top-down mechanism (the impact of the stratospheric atmosphere on tropospheric atmosphere).Moreover, with CESM, the response of the plateau snowfall rate to SA was verified.In fact, the two mechanisms worked together to affect snow cover on the QTP and the atmospheric circulation.However, the simulation results could not present such a combined effect.The reason was very likely that CESM was inappropriate, especially in the description of SA-triggered chemical and physical processes in the stratosphere.The model should be improved in the future to better describe the impact of SA on snow cover on the QTP.
It should be noted that compared with human activities and the internal variability of the climate, the impact of SA on the climate is relatively weak.Only when other influence factors are offset will SA begin to play an important role in climate anomalies.However, before this hypothesis is proved, we need to be prepared to learn more about the knowledge of SA's impact on the climate system.This requires the development and improvement of coupled climate models to better understand the physical mechanisms by which SA affects the climate.
Finally, the results obtained above are summarized as follows: (1) Solar activity could directly affect the amount and spatial distribution of snow cover on the QTP, which was consistent with previous results.We focused on analyzing the spatial pattern of correlations between the SRF and snow cover on the QTP.We found that when the SA was strong, the number of snow days and the amount of snow cover on the QTP in winter and in spring increased significantly.The contemporaneous and lag 1-year correlations were the most significant, and then the correlation gradually weakened with time.
(2) In addition to its direct impact, SA can also affect snow cover on the QTP indirectly.
SST anomalies at 0-200 m in the Pacific and the Atlantic respond to SA significantly, as does the LHF in the Pacific, resulting in abnormal sea-air interactions.The abnormal sea-air heat exchanges correlate significantly with the abnormal precipitation over the QTP.Generally speaking, in peak (valley) years of SA, winter precipitation over the QTP increases (decreases); in HSAYs (LSAYs), SST anomalies respond to SA more (less) significantly, while winter precipitation over the QTP increases (decreases) but summer precipitation decreases (increases).Therefore, we deduced that SA causes plateau snow cover anomalies indirectly by influencing the sea-air interaction.The above analysis regarded the bottom-up mechanism relevant to the impact of SA on plateau snow cover.Plateau snow cover anomalies caused by SA could further affect summer precipitation in the following flood season in China [30] through their nonlinear amplified effect on the atmosphere.In this sense, the effect of SA is amplified.(3) In addition to the bottom-up mechanism, the research results showed that there is also top-down mechanism relevant to the impact of SA on snow cover over the QTP.
In peak (valley) years of SA, the temperature in the stratosphere of the Arctic region drops (rises) while the low-and mid-latitude temperatures rise (drop) and the westerly wind strengthens around mid-latitude, which is conducive to AO presenting a positive (negative) phase.In HSAYs (LSAYs), long waves are transferred downward (upward) from the stratosphere (troposphere) to the troposphere (stratosphere) vertically.A strong divergence (convergence) of the EP flux exists between 45 • N and 70 • N of the stratosphere, and the zonal westerly wind strengthens (weakens).Such a distribution of the zonal wind along the latitude is favorable for the formation of the positive (negative) phase of the AO.(4) As the result of the interaction between the stratospheric atmosphere and tropospheric atmosphere, the zonal wind field at 500 hPa is featured with a significant easterly (westerly) wind anomaly in winter over the QTP in HSAYs (LSAYs).Meanwhile, a significant westerly (easterly) wind anomaly was identified to the south of the QTP over the Indian Peninsula, the Bay of Bengal, and the Arabian Sea.This anomaly was favorable for the generation of abnormal (anti)cyclones, and hence favorable for the increase (decrease) in plateau precipitation, leading to the increase (decrease) in plateau snow.Triggered by SA, such a meridional shear of the zonal wind can last until spring, so it continues to influence the spring snow anomaly over the QTP.(5) By using CESM of the USA (NCAR CESM2.1.1),this paper carried out a fully coupled numerical experiment that included TSI variation during 1979-2009 to analyze the impact of SA on snow cover on the QTP and the atmospheric circulation.The simulation results indicated that in HSAYs (LSAYs), the snowfall rate rose (fell) in the main part of the QTP.This finding verified the results of the previous observation analysis and proved the reliability of the conclusions.
The analysis of the mechanism of the impact of SA on plateau snow cover was based on the analysis of observation data.In the future, improved climate system models will be selected for numerical experiments to verify the above analysis results.The climate model needs to have a good simulation performance in simulating upper atmospheric chemical and physical processes, including those of the ionospheric atmosphere.It also needs to include an objective multilevel air-sea coupling system.Moreover, the spatial scope of the snow cover will not be limited to the QTP and will be extended to involve Eurasian snow, including the sea ice of the Arctic.The possible study findings will be more helpful in understanding the impact of SA on the entire cryosphere objectively.

Figure 1 .
Figure 1.The area of the Qinghai-Tibet Plateau.(Units:latitude is °N, and longitude is °E).

Figure 1 .
Figure 1.The area of the Qinghai-Tibet Plateau.(Units:latitude is • N, and longitude is • E).

Figure 2 .Figure 2 .
Figure 2. Time series of standardized SRF data with the solid (empty) circles representing peak (valley) years of SA (Song et al., 2016).(•represents the peak year of solar activity, and ○ represents the valley year of solar activity.).

Figure 3 .Figure 3 .
Figure 3. (a) Differences in winter precipitation in China between peak and valley years of SRF.The black dots represent stations that passed the 0.1 significance test.(b) Differences between high and Figure 3. (a) Differences in winter precipitation in China between peak and valley years of SRF.The black dots represent stations that passed the 0.1 significance test.(b) Differences between high and low SRF years of the plateau snow days in winter and (c) the plateau snow days in spring (• indicates the stations that passed the 0.1 significance test).(Units: latitude is • N, and longitude is • E).

Figure 4 .
Figure 4. Spatial distribution of the synchronous and lag correlation between winter snow over the QTP and SRF during 1951-2017: (a) in the same period; (b-k) in the lag of 1-10 years.The black dots represent the grid points that passed the 0.05 significance test.(Units:latitude is ° N , and longitude is ° E.)

Figure 4 .
Figure 4. Spatial distribution of the synchronous and lag correlation between winter snow over the QTP and SRF during 1951-2017: (a) in the same period; (b-k) in the lag of 1-10 years.The black dots represent the grid points that passed the 0.05 significance test.(Units:latitude is • N and longitude is • E).

Figure 5 .
Figure 5. Differences in latent heat flux between peak and valley years of SRF in winter.Stippling indicates the grid points that passed the 0.05 significance test W/m 2 ).(Units: latitude is °N, and longitude is °E.).

Figure 5 .
Figure 5. Differences in latent heat flux between peak and valley years of SRF in winter.Stippling indicates the grid points that passed the 0.05 significance test (unit: W/m 2 ).(Units: latitude is • N, and longitude is • E).

Figure 6 .
Figure 6.Correlation between latent heat flux anomaly index and winter precipitation in China Stippling indicates stations that passed the 0.05 significance test.(Units: latitude is °N, and longitude is °E).An LHF anomaly caused by SRF during the sun's 11-year cycle can also affect sum mer precipitation in the following summer season of China, hence extending and ampli fying the impact of solar role in the climate from winter to the following summer.Mean while, by impacting the global vertical ocean temperature, SA could affect winter and summer precipitation over the QTP indirectly.This also amplifies the impact of SA.Here the WOD09 data set (World Ocean Atlas, NCAR Climate Data Guide, 1955-2016) (grid point resolution: 1° × 1°) of STA data from the National Oceanographic Data Cente (NODC, USA) was utilized to determine the 5-year running mean of STA data at each level of the vertical range from 0 m to 700 m.In this way, ENSO signals were eliminated The calculation results showed that the STA at 0-200 m responded to the SRF more sig nificantly than at 300-700 m.Using the accumulated STA at 0-200 m, we analyzed the composite difference in STA response to the SRF between high and low years, as shown

Figure 6 .
Figure 6.Correlation between latent heat flux anomaly index and winter precipitation in China.Stippling indicates stations that passed the 0.05 significance test.(Units: latitude is • N, and longitude is • E).

Figure 7 .
Figure 7.The composite mean difference (CMD) of 0-200 m accumulated annual STA response to SRF after 5-year running mean.The A-F regions had significant responses; a warm color indicates a positive response and a cold color indicates a negative response.All regions passed the 0.05 significance test of red noise (as verified by the difference value between peak and valley years of SRF).(Units: latitude is • N and longitude is • E).

Figure 8 .
Figure 8. Correlation of STA index Z with (a) winter precipitation rate and (b) summer precipitation rate in East Asia.

Figure 8 .
Figure 8. Correlation of STA index Z with (a) winter precipitation rate and (b) summer precipitation rate in East Asia.Shading indicates the anomaly regions that reached a 0.1 significance level; the contour interval is 0.5%.(Units: latitude is • N and longitude is • E).

Figure 9 .
Figure 9.The composite mean difference (CMD)-derived spatial patterns of (a) winter temperature and (b) winter geopotential height response to SRF.Shading indicates that the regions passed the 0.1 and 0.05 significance tests, respectively.(The solid blue line indicates the positive value, and the dot red line indicates the negative value.Units: latitude is °N, and longitude is °E).

Figure 10 .
Figure 10.Profile of zonal wind differences in winter between high and low SRF years.Light and heavy shading indicate the regions that passed the 0.1 and 0.05 significance tests, respectively (Units: latitude is ° N , and the vertical pressure is hPa).

Figure 9 .
Figure 9.The composite mean difference (CMD)-derived spatial patterns of (a) winter temperature and (b) winter geopotential height response to SRF.Shading indicates that the regions passed the 0.1 and 0.05 significance tests, respectively.(The solid blue line indicates the positive value, and the dot red line indicates the negative value.Units: latitude is • N, and longitude is • E).

Figure 10 Figure 9 .
Figure10presents the profile of the zonal wind difference in winter between high and low SRF years.It shows that westerlies prevailed in the atmosphere from the stratosphere to the troposphere at 40-60 • N of the mid-latitudes, easterlies prevailed at 25-40 • N, and an east wind prevailed in the troposphere and upper middle stratosphere north of 65 • N.These phenomena indicated that in the middle latitudes of the stratosphere, the zonal westerly wind can propagate into the troposphere, strengthening the westerlies in the mid-latitudes of the troposphere and enhancing the annular mode in the Northern Hemisphere.Thus, westerly winds prevailed in the tropical stratosphere during HSAYs and easterly winds prevailed during LSAYs.

Figure 10 .
Figure 10.Profile of zonal wind differences in winter between high and low SRF years.Light and heavy shading indicate the regions that passed the 0.1 and 0.05 significance tests, respectively.(Units: latitude is ° N , and the vertical pressure is hPa).

Figure 10 .
Figure 10.Profile of zonal wind differences in winter between high and low SRF years.Light and heavy shading indicate the regions that passed the 0.1 and 0.05 significance tests, respectively.(Units: latitude is • N and the vertical pressure is hPa).

Figure 11 .
Figure 11.EP flux (arrow) and divergence (shading) in winter of (a) HSAYs and (b) LSAYs.Stippling indicates the grid points that passed the 0.1 significance test).Units: latitude is °N, and the vertical pressure is hPa).

Figure 11 .
Figure 11.EP flux (arrow) and divergence (shading) in winter of (a) HSAYs and (b) LSAYs.Stippling indicates the grid points that passed the 0.1 significance test).Units: latitude is • N, and the vertical pressure is hPa).

Figure 12 .
Figure 12.Differences between high and low SRF years of 500 hPa zonal wind field (a) in winter and (b) in spring.Stippling indicates the grid points that passed the 0.1 significance test.(Units: latitude is °N , and longitude is °E).

Figure 12 .
Figure 12.Differences between high and low SRF years of 500 hPa zonal wind field (a) in winter and (b) in spring.Stippling indicates the grid points that passed the 0.1 significance test.(Units: latitude is • N and longitude is • E).

Figure 13 .
Figure 13.Differences in high and low SA years of (a) annual net solar flux (unit: W/m 2 ) at the top of the atmosphere and (b) annual temperature (unit: °C) at the reference height.Stippling indicates the grid points that passed the 0.05 significance test, shading indicates plateau areas, solid and red lines represent the isolines greater than zero, and dashed and blue lines represent the isolines less than zero.(Units: latitude is °N, and longitude is °E).

Figure 13 .
Figure 13.Differences in high and low SA years of (a) annual net solar flux (unit: W/m 2 ) at the top of the atmosphere and (b) annual temperature (unit: • C) at the reference height.Stippling indicates the grid points that passed the 0.05 significance test, shading indicates plateau areas, solid and red lines represent the isolines greater than zero, and dashed and blue lines represent the isolines less than zero.(Units: latitude is • N, and longitude is • E).

Atmosphere 2022 , 21 Figure 14 .
Figure 14.Spatial distribution of annual average differences between high and low TSI years of (a) geopotential height field at 500 hPa (unit: m), (b) geopotential height field at 100 hPa (unit: m), and (c) vertical velocity ω (unit: Pa/s), and (d) snowfall rate (unit: mm/d).Stippling indicates the grid points that passed the 0.05 significance test; shading indicates plateau areas.(The red line indicates the positive value, and the blue line indicates the negative value).

Figure 14 .
Figure 14.Spatial distribution of annual average differences between high and low TSI years of (a) geopotential height field at 500 hPa (unit: m), (b) geopotential height field at 100 hPa (unit: m), and (c) vertical velocity ω (unit: Pa/s), and (d) snowfall rate (unit: mm/d).Stippling indicates the grid points that passed the 0.05 significance test; shading indicates plateau areas.(The red line indicates the positive value, and the blue line indicates the negative value).

Table 1 .
1 • .The detailed model parameters are shown in Table1.The atmospheric component of the CESM was the Community Atmosphere Model, version 6 (CAM6).It used the finite volume (FV) dynamical core with a hybrid sigma coordinate.The Community Land Model (CLM5) was chosen for the land component mode, which updated the parameter settings for new snow density to provide a more realistic description of the effect of temperature and wind field on snow density.The ocean component was the Parallel Ocean Model parameters.

Table 2 .
The extent of anomalous latent heat flux regions.

Table 3 .
Extents and significant period of 0-200 m accumulated annual STA response to SRF.