The solar activity dependence of nonmigrating tides in electron density at low and middle latitudes observed by CHAMP and GRACE

In this paper we use more than a decade of in situ electron density observations from CHAMP and GRACE satellites to investigate the solar activity dependence of nonmigrating tides at both low and middle latitudes. The results indicate that the longitudinal patterns of F region electron density vary with season and latitude, which are exhibiting a wavenumber 4 (WN4) pattern around September equinox at low latitudes and WN1/WN2 patterns during local summer at the southern/northern middle latitudes. These wave patterns in the F region ionosphere can clearly be seen during both solar maximum and minimum years. At low latitudes the absolute amplitudes of DE3 (contributing to the WN4 pattern) are found to be highly related to the solar activity, showing larger amplitudes during solar maximum years. Similarly a solar activity dependence can also be found for the absolute amplitudes of D0, DW2 and DE1 (contributing to the WN1 and WN2 pattern) at middle latitudes. The relative amplitudes (normalized by the zonal mean) of these nonmigrating tides at both low and middle altitudes show little dependence on solar activity. We further found a clear modulation by the quasi-biennial oscillation (QBO) of the relative DE3 amplitudes in both satellite observations, which is consistent with the QBO dependence as reported for the E region temperatures and zonal wind. It also supports the strong coupling of the low-latitude nonmigrating tidal activity between the E and F regions. However, the QBO dependence cannot be found for the relative amplitudes of the nonmigrating tides at middle latitudes, which implies that these tides are generated in situ at F region altitudes.


Introduction
The ionosphere-thermosphere (IT) system is a crucial part of the solar-terrestrial interaction region and is an important research objective for space weather.It is well known that the IT system varies noticeably with different spatial and temporal scales.This variation is caused not only by the effect of solar radiation, interplanetary, magnetospheric processes, but also by the waves propagating up from the lower atmosphere.The combined effects of the aforementioned various driving mechanisms determine the changes of composition, temperature and density of the ionized and neutral constituents in the IT system (Schunk and Sojka, 1996).
In the past decade much effort has been devoted to investigating the longitudinal variations of the ionospheric and thermospheric quantities, which are often related to lower atmospheric processes, especially to the influence of various waves.These waves are referred to as atmospheric migrating and nonmigrating tides with periods that are harmonics of a solar day (Forbes, 1995).The migrating tides are mainly generated in the troposphere and stratosphere through the absorption of solar energy by tropospheric water vapor and stratospheric ozone (Oberheide et al., 2002).From the perspective of a ground-based observer they propagate westward with the motion of the Sun and thus will not cause a longitude variability at a fixed local time around a constant latitude circle.The nonmigrating tides do not migrate with the Sun.Therefore they will cause longitudinal variations of the IT system if viewed at a fixed local time.According to previous studies nonmigrating tides are mainly excited by zonal asymmetries of driving mechanisms (e.g.topography, land-sea differences, longitude dependence of absorb-ing species) (Forbes et al., 2003), or by nonlinear interactions between the migrating diurnal tide and planetary waves (Hagan and Roble, 2001), or gravity waves (McLandress and Ward, 1994), or by the latent heat release in the troposphere (Hagan and Forbes, 2002).
In the universal time (UT) frame, these tides can be expressed as (Forbes et al., 2006) A n, s cos(n t + sλ − φ n, s ), (1 where t is the universal time, is the rotation rate of the Earth ( = 2π/24 h), λ is the longitude, n (n = 1, 2, 3, . ..) is the harmonic of a solar day, s (s = . . .−3, −2, −1, 0, 1, 2, 3 . ..) is the zonal wave number, A n, s is the amplitude, φ n, s is the phase of the tides and defined as the time when the maximum passes the zero degree longitude.
When viewed in the local time frame, Eq. ( 1) becomes where t LT (t LT = t + λ/ ) is local time.When s = n, the components denote migrating tides and are longitude independent.While when s = n, the components represent nonmigrating tides, with wave number |s − n| observed in longitude for a fixed local time.For labeling the tidal components we use the common convention.The first letter D, S, or T stands for diurnal, semidiurnal or terdiurnal; the second letter E or W for eastward or westward propagation and the final number quantifies the azimuthal wave number.D0 represents a wave that is increasing and decreasing simultaneously at all longitudes at a diurnal period.Stationary planetary waves are labeled SPWs and the number at the end quantifies the maxima around the globe.The phase of the SPW defines the longitude where the wave crest occurs.At low and equatorial latitudes the WN4/WN3 patterns of the ionospheric and thermospheric quantities and the related tidal components have been widely studied (Sagawa et al., 2005;Immel et al., 2006;England et al., 2006England et al., , 2010;;Lühr et al., 2008Lühr et al., , 2012;;Kil et al., 2007;Liu et al., 2009;Häusler and Lühr, 2009;Zhang et al., 2010;Xiong and Lühr, 2013;Xiong et al., 2014).These tides originate from the lower atmosphere and can propagate upward to the upper atmosphere, either by direct propagation or through the E region wind dynamo coupling mechanism (e.g.Wu et al., 2012).The results showed that at middle latitudes the most prominent tidal components are D0, DW2, SPW1 and DE1.From model simulations Jones Jr. et al. (2013) suggest that at middle latitudes the nonmigrating tides could be generated in situ through ion-neutral interactions due to the longitude dependence of the ionosphere, imposed by the geomagnetic field configuration.Some recent studies have reported the solar activity dependence of the nonmigrating tides at low latitudes.With 5 years (2007)(2008)(2009)(2010)(2011) observations of total electron content (TEC) from COSMIC Chang et al. (2013) reported that the absolute amplitudes of DE3 are directly related to the solar activity, while the relative amplitude (normalized by the zonal mean) is inversely related to the solar activity through 2010.The solar activity dependence of DE3 has also been checked by S. Wang et al. (2015) by using 14 years of TEC data from global ionosphere maps (GIMs) generated at JPL (Jet Propulsion Laboratory).They found that the absolute amplitudes of DE3 are highly related to the solar activity.However, the relative amplitudes show little dependence on solar activity.CHAMP and GRACE provide more than 10 years in situ electron density (Ne) observations, which offer us a good opportunity for investigating the solar activity dependence of the nonmigrating tides in the topside ionospheric F region.We will further compare the behavior of tides at low and middle latitudes.
In Sect. 2 we firstly introduce the data sets (ISDC, 2015) and processing approach.Then we will present the observational results from CHAMP and GRACE.The interpretation of the results and the comparison with earlier studies are given in Sect. 4.

Data sets and processing approach
The CHAMP spacecraft was launched on 15 July 2000 into a near-circular polar orbit (inclination: 87.3) with an initial altitude of 456 km.By the end of the mission, 19 September 2010, the orbit had decayed to 250 km.The orbital period is about 93 min, thus circling the Earth about 15 times per day.The orbital plane covers all local times within 131 days, at a rate of 1 h in local time per 11 days (Reigber et al., 2002).The Planar Langmuir Probe (PLP) on board the satellite measured the electron density and temperature every 15 s.The PLP electron density readings have been validated by comparison against digisonde observations at Jicamarca (McNamara et al., 2007).
GRACE consists of two spacecraft and was launched on 17 March 2002 into a near-circular, polar orbit (inclination: 89.5 • ) with an initial altitude of about 525 km (Tapley et al., 2004).At the end of 2014 the altitude of the two spacecraft were around 420 km.The local time of the orbital plane precesses by 4.5 min every day taking the mission 161 days to sample all local times.The two spacecraft follow each other at a distance of about 170-220 km.The total electron content (TEC) between the two spacecraft can be deduced from the K Band Ranging (KBR) data.When dividing the horizontal TEC by the distance between the spacecraft, the averaged electron density can be derived (Xiong et al., 2010;Lühr and Xiong, 2010).The GRACE electron density data are further validated by EISCAT and the Incoherent Scatter Radars (ISR) at Millstone Hill and Arecibo (Xiong et al., 2015a).
As already introduced, the amplitudes of the different tides at the lower atmosphere vary with seasons.The WN4 pattern at low latitudes is most prominent around September equinox, while the WN1/WN2 patterns at middle latitudes are most prominent during local summer in the Southern/Northern hemispheres.Therefore, observations from the late summer season have been selected for low-latitude research, and December/June solstice data have been selected for the southern/northern middle latitudes.To cover the 24 h local time, 131 and 161 days observations are needed for CHAMP and GRACE, respectively.It should be noticed that for the late summer season the 131/161 days centered on 15 August are selected because WN4 patterns at low latitude are most prominent during this period of time.Hereafter we will term this season September equinox.Days with planetary geomagnetic activity index (Kp) exceeding 3 are not included, to reduce the influence of magnetic disturbances.For each selected period we sort the electron density data into geomagnetic latitude (1 • ) × geographic longitude (15 • ) × local time (1 h) bins.For low latitudes, 5 to 20 • MLAT are selected to focus on the crest region of the Equatorial Ionization Anomaly (EIA).For middle latitudes, 40 to 60 • MLAT are selected.For the selected latitude bins, the longitudinal mean values have been subtracted hour by hour, in order to suppress the migrating tides.The longitudinal mean-free data is labeled Ne in the present study.For the next step these longitudinal mean-free data are further processed by the onedimensional Fourier transform which will bring forth the sum of observed tidal signatures for each wave number up to 4. A similar approach for obtaining the amplitudes and phases of various tidal signatures for the F region electron density can be found in Xiong and Lühr (2014).Based on the aforementioned approach the amplitudes of all interesting tidal components can be derived for each period, and their amplitudes will be further used for analyzing the solar activity dependence.

Results
Figure 1 shows the zonal (longitudinally averaged) mean electron density observed by CHAMP and GRACE satellites between ±70 • MLAT from 2000 to 2014.The blank areas are missing data.The orbital height at ascending node, as well as the solar radio flux index P10.7 are also presented in the upper panels.The P10.7 is defined as P10.7 = (F10.7 + F10.7A)/2, where F10.7A is the 81 days averaged value of F10.7.We can see that the solar maximum of the solar cycle occurred around 2000-2002, with P10.7 exceeding 250 sfu.Then the solar activity gradually decreases, and the solar minimum appears around 2008-2009, with P10.7 below 70 sfu.The recovery of solar activity appears around 2010, to the end of 2014 the P10.7 reaches occasionally 200 sfu.During their mission periods CHAMP has decayed from 450 km to about 300 km, and GRACE decayed from 500 km to about 400 km.The zonal mean electron densities from both satellites show a clear dependence on the solar activity.Take CHAMP for example, the maximum electron density exceeds 1.4 × 10 6 cm −3 around 2002, while the maximum value has decreased roughly to 6 × 10 5 cm −3 during 2008-2009 in spite of the orbital decay.The seasonal variation of the F region electron density can also be found.It maximizes around equinox seasons and minimizes around June solstice.For the latitude variation, the EIA crests gen- erally appear between ±8 and 20 • MLat, which is somewhat closer to the magnetic equator during solar minimum years.It is noteworthy that the obvious increase of GRACE electron density during 2014 in fact is additionally reinforced by its orbit decay.
For illustrating the longitudinal wave patterns, taking CHAMP as an example, Fig. 2   wave features the WN4 pattern at low latitudes is found to be prominent only during daytime.At low and equatorial latitudes most of the tides are generated in the lower atmosphere and propagate to the upper atmosphere through the E region wind dynamo coupling mechanism.Therefore the fountain effect is prominent only during daytime.At middle latitudes Ne is mostly positive in the southern hemisphere between 30 and 150 • W around 22:00 LT, while in the Northern Hemisphere two peaks of Ne can be found between 90-150 • W and 30-90 • E during 12:00-22:00 LT.These longitudes with higher electron density during summer night corresponds well with the mid-latitude summer night anomaly (MSNA), as reported by previous studies (e.g.Lin et al., 2010;Liu et al., 2010).The longitudinal patterns at low and middle latitudes from GRACE observations are similar to that of CHAMP, but are not shown here.
By the aforementioned method we obtained the amplitudes and phases of different tidal components that contribute to the longitudinal structures.It is found that the longitudinal WN4 structures of Ne at low latitude are dominated by the nonmigrating tide DE3, and the contributions from the other components (such as SPW4, DW5, SE2, SW6, TE1 and TW7) are somewhat less important.For the middle latitudes the longitudinal WN1 and WN2 patterns in the Southern and Northern Hemisphere are dominated by the nonmigrating tides D0, DW2, and DE1.Thus in the following we only focus on these four tidal components, to reveal their variation with the solar activity.
We first studied DE3 at low latitudes.Figure 3 1c and d).Therefore, the DE3 amplitudes are not obtained for these years.The relative amplitudes again show little dependence on the solar activity.However, except for the amplitude of DE3 in 2008, larger relative amplitudes during even years appear again.In fact, this 2-year oscillation in the DE3 tides has earlier been reported by Forbes et al. (2008) for migrating tides.They found a prominent quasi-biennial oscillation (QBO) related interannual variation in the temperature at 100 km.It is known that  QBO dominates the variability of the equatorial stratosphere, and it is easily seen as a downward propagating easterly or westerly wind regime, with a variable period of about 27 months (Baldwin et al., 2001).For comparison, the temporal variation of the observed equatorial stratospheric monthly mean zonal wind at level 20 hPa is added as the third panel in Fig. 3.The positive values of wind velocity mean westerly.It is shown that the relative amplitude of DE3 is larger when the zonal wind in the stratosphere is westerly and is smaller when the zonal wind is easterly.A clear QBO effect can be found in the relative amplitudes of DE3 from both satellites observation.We will come back to this point in the discussion part.
Figure 4 shows the temporal variation of the absolute amplitudes of D0 and DW2 components at southern middle latitudes during local summer.It is clear that the absolute amplitudes of the D0 and DW2 components, which contributed to the prominent WN1 patterns, are closely correlated with the variation of P10.7, especially in case of D0. Figure 5 shows the temporal variation of absolute amplitudes of D0 and DE1 components, which contribute to WN2 patterns at northern middle latitudes during local summer.It can be seen that maximum amplitudes of DE1 and D0 at CHAMP altitude emerge in 2002 with a value of about 0.4 × 10 5 cm −3 .Clearly, DE1 and D0 have the same trend of variation with the solar activity.As presented in the right panel of Figs. 4 and 5, the absolute amplitudes of the three tides from GRACE observations show similar solar activity dependence as that from CHAMP.
Figure 6 shows the relative amplitudes of D0 and DW2 at southern middle latitudes during local summer from both satellites observations.Figure 7 is the same as Fig. 6, but for DE1 and D0 at northern middle latitudes local summer.No clear dependence on solar activity can be found, and at Southern Hemisphere the relative amplitudes are varying around 48 and 18 % for the D0 and DW2 tides, respectively.For Northern Hemisphere, the amplitudes are about 10 and 16 % for the D0 and DE1 components, respectively.

Discussion
The analysis presented above has provided the solar activity dependence of the nonmigrating tides in F region electron density at both low and middle latitudes.Within this section we will discuss our results in the context of previous results related to nonmigrating tides in the ionospherethermosphere system.

Nonmigrating tides at low latitudes
During the past decade lots of previous studies focused on the longitudinal patterns of various ionospheric and thermospheric quantities at the equatorial region.The related tidal components and the possible coupling mechanisms between the mesosphere and lower thermosphere (MLT) and the upper atmosphere have also been proposed (e. g.Sagawa et al., 2005;Immel et al. 2006;England et al., 2006England et al., , 2010;;Lühr et al., 2008Lühr et al., , 2012;;Xiong and Lühr, 2013).These tides are usually excited by thermal sources, stationary planetary wave-  tide (e.g., Angelats i Coll and Forbes, 2002) and tide-tide interactions (Hagan et al., 2009;Oberheide et al., 2011).They further propagate upward to the upper atmosphere either by direct propagation (Häusler and Lühr, 2009) or through the E region wind dynamo coupling mechanism (e.g.Wu et al., 2012).
At daytime the E region zonal electric field transports the plasma upward, thereby depleting the equatorial F region at the equator and accumulating plasma in the EIA crest regions.The modulation of the E layer electric field is mapped with a certain time delay into the F region EIA (Stolle et al., 2008).An interesting feature in our Fig.2a is that the WN4 longitudinal pattern in electron density at the crests of EIA is quite prominent during daytime hours but vanishes at post-midnight hours.During higher solar activity years (2001)(2002)(2003)(2004), the WN4 pattern at the EIA crest region can also be found between 18:00 and 21:00 LT.This could be explained by the pre-reversal enhancement (PRE), which is defined as a sharp increase in vertical ion drift after sunset via an enhancement of the eastward electric field strength in the F region (e.g., Woodman, 1970;Fesen et al., 2000).The increased E × B drift will cause an electron density trough above the dip equator.While during solar minimum, the PRE effect is low, therefore, the WN4 patterns in the EIA crest region are not so clear, especially during the years 2008-2009.It could be considered as an indication that when the E layer disappears during nighttime, other mechanisms for interaction between ions and neutrals may appear in the F region.

Nonmigrating tides at middle latitudes
Different from the WN4 at low latitudes, WN1/WN2 patterns at middle latitudes show little local time difference between solar maximum and minimum years, except for a reduced magnitude.It implies that at middle latitudes the wave patterns in the F region may be less related to the E region.Jones Jr. et al. (2013) suggested that at middle latitudes the nonmigrating tides including DE1, D0, and DW2 could be generated in situ through ion-neutral interactions due to the longitude-dependent ionospheric properties imposed by the geomagnetic field configuration.Similar WN1 patterns in thermospheric mass density and zonal wind are also found at the southern middle latitudes, and the most prominent tidal components are D0 and DW2 (e.g.Xiong et al., 2015b).One interesting feature of D0 is that it shows an anti-phase diurnal variation between hemispheres both in the F region electron density as well as in the thermospheric mass density and zonal wind (Xiong et al., 2015b).By employing the global ionosphere-thermosphere model (GITM) simulation, H. Wang et al. (2015) found that the thermospheric zonal wind plays an important role in causing the WN1/WN2 patterns of the F region electron density at middle latitudes, but additional contributions are needed to explain the full extension of observed tides.

Solar activity dependence of the nonmigrating tides
The F region background electron density is largely controlled by the solar activity, especially the extreme ultraviolet (EUV) radiation.That is why the absolute amplitudes of the F region nonmigrating tides are directly related to P10.7 at both low and middle latitudes, as shown in Figs.3-5.This is in accord with the result of Chang et al. (2013).By analyzing the TEC observations from the FORMOSAT-3/COSMIC during 2007-2011 they reported that the absolute amplitudes of all the tides contributing to WN4 (DE3, SE2 and SPW4) are directly related to the level of solar activity.Similar results have been reported by S. Wang et al. (2015), based on 14 years (1999-2013) of TEC observations from the global ionospheric maps of JPL.They also found that the absolute amplitudes of DE3 and SPW4 exhibit a saturation effect on the solar activity.This apparent saturation of nonmigrating tidal amplitudes could be explained by the saturation of the TEC values at high solar activity (Liu et al., 2009).In our in situ N e observations no saturation effect has been found for the absolute DE3 amplitudes.When looking at the relative amplitudes of these tides at low latitudes, S. Wang et al. (2015) found little dependence on solar activity, while Chang et al. (2013) claimed that relative amplitudes of WN4 are inversely related to solar activity through 2010.From their Fig.5a we can see that at 15 • N the relative amplitudes of DE3 show little variations from 2007 to 2011, and the peak values are around 7 % during the 5-year period.Compared to the solar maximum from 2001 to 2003, the solar activity during 2007-2011 in fact is rather low, therefore it may not be sufficient from their data to conclude that relative amplitudes of DE3 are inversely related to solar activity.The relatively larger amplitude of DE3 in 2008 (as shown in their Fig.4a) in our eyes is more reasonably explained by the modulation of the propagation conditions from the troposphere and stratosphere by QBO.Additionally, the relatively small amplitude of DE3 in 2010 may also be related to QBO.It can be seen from our Fig. 3 that the transition time of wind from easterly to westerly during 2009 to 2010 is longer than during previous cycles.Forbes et al. (2008) found a QBO variation of order ±10-15 % in the migrating diurnal and semidiurnal tidal amplitudes in temperature measurements obtained by the SABER instrument on board TIMED between 100 and 116 km.Similar QBO variations in the DE3 tide amplitude have also been reported by Oberheide et al. (2009), from the zonal wind observations at 105 km.The QBO effect described above is mainly found in the E region observations.As shown in our Fig. 3, when normalized by the zonal mean values the QBO effect in the DE3 amplitudes can also be found in the F region electron density, especially during 2001 to 2007.The relatively larger amplitudes of the relative DE3, during the periods of westerly stratospheric mean zonal winds are consistent with the results as reported by Forbes et al. (2008) and Oberheide et al. (2009).The similarity of QBO effects found in both E and F regions also supports the inference that the low latitude DE3 nonmigrating tides at all levels are driven from below.
At middle latitudes, the absolute amplitudes of the nonmigrating tides are also highly related to the solar activity.However, the QBO effect cannot be found in their relative amplitudes, which implies that the nonmigrating tides in the F region are not strongly coupled to E region dynamics at these latitudes.

Summary
In this study we have used more than 10 years of F region electron density observations from CHAMP and GRACE to analyze the solar activity dependence of nonmigrating tides at low and middle latitudes.The main results are summarized as follows.
The longitudinal WN4 pattern during September equinox at low latitudes and the WN1/WN2 patterns during local summer at the southern/northern middle latitudes in the F region electron density are robust features throughout the 23/24 solar cycles.
The absolute amplitudes of DE3 at low latitudes as well as DE1, D0 and DW2 at middle latitudes are highly related to the solar activity, all show larger amplitudes during solar maximum years.While their relative amplitudes (normalized by the zonal mean electron density) show little dependence on solar activity.
At low latitudes, a clear modulation by the QBO can be found in the relative amplitudes of DE3 both from CHAMP and GRACE observations, which is consistent with the QBO effect as reported for E region quantities.It also supports the strong coupling of the nonmigrating tides between the E and F region.However, no QBO effect can be found in the tidal components DE1, D0 and DW2 at middle altitudes, which implies that these tides are generated in situ at F region altitudes.

Figure 1 .
Figure 1.(a) The solar radio flux index P10.7,(b) the ascending node orbital heights of CHAMP and GRACE satellites, as well as the zonal (longitudinally averaged) mean electron density observed by (c) CHAMP and (d) GRACE satellites between ±70 • MLAT from 2000 through 2014.

Figure 2 .
Figure 2. Local time vs. longitudinal distribution of Ne variations observed by CHAMP at low latitudes (a, September equinox), northern middle latitudes (b, June solstice), and southern middle latitudes (c, December solstice) from 2001 to 2009.
presents the local time vs. longitudinal distribution of Ne at low latitudes (top panel, September equinox), northern middle latitudes (mid-dle panel, June solstice), and southern middle latitudes (bottom, December solstice) from 2001 to 2009.Here Ne means the longitudinal mean value of electron density has been subtracted hour by hour.We can clearly see that Ne shows longitudinal WN4, WN2/WN1 patterns at low and middle latitudes, respectively.These robust wave features can be found throughout the CHAMP mission period with somewhat larger amplitudes during higher solar activity years.When looking at the local time variations of these

Figure 3 .
Figure 3. Temporal variation of absolute (top) and relative (middle) amplitudes of the DE3 tidal component observed by CHAMP (left) and GRACE (right) for September equinox at the northern EIA crest region, as well as averaged values of P10.7 solar radio flux for the corresponding same time period.(bottom) Monthly mean of stratospheric (20 hPa) zonal wind at equatorial latitudes.
presents the temporal variation of the absolute (top) and relative (middle) amplitudes of DE3 components observed by CHAMP (left) and GRACE (right) for September equinox at the northern EIA crest region.For deriving the relative DE3 amplitudes, Ne is first normalized by the zonal mean value hour by hour, then decomposed with the same approach described above.The mean values of P10.7 during the corresponding study period are also presented.The absolute amplitude of DE3 observed by CHAMP is gradually decreasing towards lower solar activity, with smallest amplitudes appearing in 2009 when the level of solar activity reaches a minimum.The absence of a point for 2008 is due to the partly lacking observations.Different from the absolute amplitudes the relative amplitudes of DE3 show little dependence on solar activity.However, the amplitudes in 2002, 2004 and 2006 are somewhat larger than the previous years.At the GRACE altitude the absolute amplitudes of DE3 from GRACE observations are also directly correlated with the solar activity.It gradually decreases with the P10.7 from 2002 to 2009 and then increases again.Similar to the CHAMP observation in 2008, the interrupted observations of GRACE from 2012 to 2014 cannot provide the 24 h coverage during September equinox (see the blanks in Fig.

Figure 4 .
Figure 4. Temporal variation of absolute amplitudes of D0 and DW2 components observed by CHAMP (left) and GRACE (right) for local summer at Southern Hemisphere middle latitudes, as well as averaged values of P10.7 solar radio flux for the corresponding period of time.

Figure 5 .
Figure 5.The same as Fig. 4, but for D0 and DE1 tidal components at local summer in the Northern Hemisphere.

Figure 6 .
Figure 6.The same as Fig. 4, but for relative amplitudes of D0 and DW2.

Figure 7 .
Figure 7.The same as Fig. 5, but for relative amplitudes of D0 and DE1.