• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    A Study on Low Frequency Electromagnetic Cyclotron Waves in the Solar Wind

    2022-08-02 08:19:26HaiFengYangGuoQingZhaoHengQiangFengGilbertPiQiangLiuLiangXiangandQiuHuanLi

    Hai-Feng Yang, Guo-Qing Zhao , Heng-Qiang Feng, Gilbert Pi, Qiang Liu, Liang Xiang, and Qiu-Huan Li

    1 Institute of Space Physics, Luoyang Normal University, Luoyang 471934, China; zgqisp@163.com

    2 Henan Key Laboratory of Electromagnetic Transformation and Detection, Luoyang 471934, China

    3 Faculty of Mathematics and Physics, Charles University, Prague 121 16, Czech Republic

    Received 2022 March 17; revised 2022 April 9; accepted 2022 April 13; published 2022 May 20

    Abstract According to Wind observations between 2005 and 2015, this paper investigates the dependences of the occurrence of low frequency electromagnetic cyclotron waves (ECWs)on the plasma parameters, the solar cycle,and the orientations of alpha-proton drift velocity(Vd)and the ambient magnetic field(B).The occurrence rates of ECWs with respect to six plasma parameters are calculated. Results show that the preferential conditions for generation of left-handed(LH)ECWs are higher proton temperature(Tp),higher proton velocity(Vp),lower proton density (Np), stronger proton temperature anisotropy (T⊥/T‖), higher normalized alpha-proton drift velocity(Vd/VA), and higher normalized alpha particle density (Nα/Np), where T⊥and T‖ refer to proton temperatures perpendicular and parallel to B, and VA is the local Alfvén velocity. For right-handed (RH) ECWs, however, the dependences on these plasma parameters are not obvious.On the other hand,it is found that the occurrence rate of LH ECWs increases as the sunspot number decreases, and decreases as the sunspot number increases. Further investigation shows that the increased occurrence rate of LH ECWs is accompanied by an increase in the medians of Vp,Vd/VA,and Nα/Np.For RH ECWs,the occurrence rate appears to be nearly unrelated to the sunspot number,suggesting a negligible correlation with the solar cycle. In addition, a modified angle is introduced to include the factors of orientation of B(anti-sunward or sunward)and angle between Vd and B,simultaneously.It is found that the occurrence of LH ECWs has the strong preferential condition that Vd is anti-sunward,while a different situation arises for RH ECWs. These results are discussed in the context of the temperature-anisotropy-driven instabilities with the effect of alpha particles.

    Key words: (Sun:) solar wind – instabilities – waves

    1. Introduction

    The solar wind is a tenuous, magnetized plasma streaming outward from the Sun and permeating interplanetary space.It is mainly composed of electrons and protons, and carries a frozen-in magnetic field.Alpha particles and other ionic species are also present at low abundance levels.These particles in the solar wind are generally collisionless and far from thermodynamic equilibrium. They exhibit non-Maxwellian property corresponding to an anisotropy (Gary 1993). These nonthermal particles can serve as free-energy sources to induce plasma kinetic instabilities in the solar wind, resulting in the generation of kinetic waves(Hollweg 1975;Schwartz 1980;Tu& Marsch 1995; Marsch et al. 2004; Marsch 2006; Cranmer 2014; He et al. 2015; Hellinger & Trávní?ek 2016; Klein et al. 2018; Wilson et al. 2018; Xiang et al. 2021; Zhao et al.2022), such as electromagnetic cyclotron waves (ECWs) near the proton cyclotron frequency.

    ECWs have been researched extensively in various space environments,such as terrestrial foreshock regions(Smith et al.1985;Wong et al.1991),magnetosheath(Schwartz et al.1996;Soucek et al.2015),and planetary magnetosphere(Russell and Blancocano 2007;Rodríguez-Martínez et al.2010).It was only in recent years that ECWs in the solar wind were observed widely.A series of papers reported the observations of coherent low-frequency ECWs with a typical frequency of 0.1–0.5 Hz in the solar wind (Jian et al. 2009, 2010, 2014, 2016; Boardsen et al. 2015; Gary et al. 2016; Zhao et al.2017, 2018, 2019a, 2019b; Pi et al. 2022). The main conclusions from the above observations are as follows: (1)These waves are transverse waves and propagate mainly in the directions quasi-(anti)parallel to the ambient magnetic field(B);(2)they have narrow frequency band near the proton cyclotron frequency;(3)some wave characteristics vary moderately when heliocentric radial distance changes from 0.3 to 1 au; (4) lefthanded(LH)ECWs are usually the dominant waves compared with right-handed (RH) ECWs in the solar wind.

    Theoretically, plasma kinetic instabilities are often investigated using the linear Vlasov–Maxwell kinetic theory, and the plasma dispersion relations and anisotropy thresholds associated with different instabilities are calculated to explain the generation and propagation of kinetic waves (Gary et al.2000a, 2000b, 2006; Hellinger & Trávní?ek 2006; Hellinger et al. 2006; Podesta & Gary 2011; Maruca et al. 2012;Verscharen & Chandran 2013; Verscharen et al. 2013; Wicks et al.2016;Xiang et al.2018a,2018b).Two kinetic instabilities driven by proton temperature anisotropies are proposed to account for the generation of the LH and RH ECWs(Gary et al.1976; Gary 1993, 2015; Kasper et al. 2002; Marsch et al.2004;Hellinger&Trávní?ek 2006;Omidi et al.2014a,2014b;Yoon 2017). Specifically, for proton temperature anisotropy with T⊥>T‖, where T⊥and T‖refer to proton temperatures perpendicular and parallel to B, proton cyclotron instability is excited to produce LH cyclotron waves in the plasma frame.For a converse temperature anisotropy with T⊥<T‖, firehose instability may arise to generate RH magnetosonic waves in the plasma frame.

    Besides temperature anisotropies, differential flows between different ionic species or between the core and beam of an ionic species are also the free-energy sources for the plasma instabilities at ion-kinetic scales, resulting in the generation of kinetic waves.ECWs can be induced by the instability of the proton distribution function, where the positive gradient in the proton distribution between the core and beam region is a freeenergy source (Wicks et al. 2016). The presence of alphaproton differential flow may affect the properties of dispersion relations for temperature-anisotropy-driven instabilities and determine the propagation directions of unstable modes(Hellinger & Trávní?ek 2006; Podesta & Gary 2011; Seough& Nariyuki 2016). When alpha-proton drift velocity Vd=Vα-Vpis directed away from the Sun, where Vαand Vpare the velocities of alpha particles and protons respectively,the LH ECWs mainly propagate away from the Sun, and the RH ECWs mainly propagate toward the Sun in the plasma frame (Podesta & Gary 2011). It should be noted that the polarization of RH ECWs propagating toward the Sun in the plasma frame turns into LH ECWs in the spacecraft frame due to the presence of a large Doppler shift(Jian et al.2009,2010).

    Despite there having been some observations and theoretical works on the occurrence of ECWs in the solar wind, the generation mechanism of ECWs is inconclusive and there is still room for further study. To the best of our knowledge, no work has yet investigated the effect of the solar cycle on the occurrence of ECWs, although the time-dependent occurrence rate of ECWs has been shown recently(Zhao et al.2017).The sunspot number changes over the period of 11 yr, and solar wind parameters vary with changes in the sunspot number,which may affect the occurrence of ECWs in the solar wind.In this paper,we try to find out how and why the occurrence rates of ECWs change in a solar cycle. In addition, the preferential local plasma conditions and generation mechanisms of ECWs are studied further in this paper, although some preferential conditions with data from the STEREO mission were revealed(Zhao et al.2017).We also study the effects of the orientations of Vdas well as B on the occurrence of ECWs in the solar wind, because it is believed that Vdplays an important role in the generation and propagation of ECWs (Podesta &Gary 2011).

    The paper is organized as follows. Section 2 introduces the data and analysis methods.Section 3 shows the dependences of the occurrence of ECWs on the plasma parameters, the solar cycle, and the orientations of Vdand B. The summary with a brief discussion is presented in Section 4.

    2. Data and Analysis Methods

    The solar wind data used in the present paper are from the Solar Wind Experiment (SWE) and Magnetic Field Investigation(MFI)instruments onboard the Wind mission,which is a solar-terrestrial mission in a halo orbit around the L1 Lagrange point.SWE measures ions in the solar wind,and the solar wind velocity, density and temperature can be deduced.The 3D velocity distribution measurements of the ion component in the solar wind are made by a pair of Faraday cup analyzers. The velocity data used in this paper are produced via a nonlinear-least-squares bi-Maxwellian fit of ion spectrum from the Faraday cups (Kasper et al. 2006). MFI investigates the structures and fluctuations of the interplanetary magnetic field and measures the intensity and direction of the magnetic field vector.

    Using high-resolution magnetic field data, we conducted a survey of ECWs occurring in the years between 2005 and 2015, which cover the whole solar cycle. We employed an automatic wave detection procedure to identify the ECWs,which was developed by Zhao et al. (2017,2018),and used in previous works (Zhao et al. 2019a, 2019b). Here we only give the main steps of the procedure: (1) The reduced magnetic helicity is calculated, which is normalized and takes values in the range from -1 to 1, (2) the presence of the waves is identified, which is associated with a band of enhanced magnetic helicity spectrum as well as a power spectrum, and(3)the waves are recorded if they satisfy an amplitude criterion of 0.1 nT.Considering a small magnitude of Vdmay result in a large uncertainty and random orientation of Vdrelative to B(Kasper et al. 2006; Alterman et al. 2018), we discard the observations with Vd/Vp<1%.

    3. Statistical Results

    Based on the Wind data, we study the dependences of the occurrence of ECWs on the plasma parameters,the solar cycle,and the orientations of Vdand B in Sections 3.1–3.3,respectively. Note that the polarization is described in the spacecraft frame throughout the paper,except that we point out the plasma frame.

    Figure 1.Occurrence rates of LH(red square line),RH(blue circle line)ECWs,and PDDs of the ambient plasma(black triangle line)with respect to Tp,Vp,Np,T⊥/T‖, Vd/VA, and Nα/Np.

    3.1. Dependence of Occurrence of ECWs on the Plasma Parameters

    To identify the preferential local plasma conditions and generation mechanisms of ECWs, we calculate the occurrence rates of ECWs, which are defined as the ratios of the ECW number and the total sample number in a certain parameter cell(bin). It is worth noting that a previous article reported the preferential conditions for generation of LH ECWs with data from the STEREO mission (Zhao et al. 2017). Compared with the previous article,data from the Wind mission are used in this paper.The data cover a longer period,more plasma parameters,and larger value ranges of plasma parameters. Thus, more detailed and comprehensive preferential conditions are expected. In Figure 1 we show the occurrence rates of LH and RH ECWs with respect to six plasma parameters,including proton temperature (Tp), proton velocity (Vp), proton density(Np), proton temperature anisotropy (T⊥/T‖), normalized alpha-proton drift velocity (Vd/VA) and normalized alpha particle density(Nα/Np),where VAis the local Alfvén velocity.The probability density distributions (PDDs) of the ambient plasma are also given in each figure, by which the event number of LH or RH ECWs in a certain cell can be determined by N×P×O, where N is the total sample number in each figure, P and O are the PDD of the ambient plasma and the occurrence rate of LH or RH ECWs.

    Examining Figure 1 overall,three conclusions can be drawn.First, the occurrence rate of LH ECWs presents significant changes with the plasma parameters and is higher than that of RH ECWs in most parameter intervals. LH ECWs take place preferentially in plasma characterized by higher Tp, higher Vp,lower Np, larger difference between T⊥and T‖, higher Vd/VA,and higher Nα/Np. Second, the preferential conditions for generation of RH ECWs are ambiguous, because the occurrence rate of RH ECWs presents little changes with the plasma parameters compared to that of LH ECWs. Third, in some parameter intervals,the occurrence rate of RH ECWs can even exceed that of LH ECWs, such as at Tp<2.5×104K,Vp<350 km s-1, Np>11 cm-3, Nα/Np<0.02.

    The ECWs are supposed to be induced by the kinetic instabilities driven by temperature anisotropies and the maximum growth rates of the instabilities can be increased in the presence of alpha-proton drift velocity. Figure 1(d) shows the occurrence rate of LH ECWs has the maximum at around T⊥/T‖~0.14 and 9, which is far from temperature isotropy.Although the occurrence rate of RH ECWs does not vary as much as that of LH ECWs, we can still find that RH ECWs have a larger occurrence rate at stronger temperature anisotropy with T⊥/T‖=0.14 or 5. Figures 1(e) and (f) show the occurrence rates of ECWs versus Vd/VAand Nα/Np. The minimum of PDD at Vd/VA=0.05 is likely due to the discarding of the events with Vd/Vp<1%. It is found that the occurrence rate of LH ECWs increases with Vd/VAand Nα/Np,which is consistent with the theory of temperature-anisotropydriven instability with the effect of alpha particles (Podesta &Gary 2011).

    3.2. Dependence of Occurrence of ECWs on the Solar Cycle

    The Sun’s activity shifts through an approximately 11 yr cycle, with sunspot number, radiation level, and ejected material changing over time. Meanwhile, the solar wind parameters vary with the changes in the Sun’s activity, which may affect the generation of ECWs in the solar wind. The period of the measurements studied in this paper is 11 yr from 2005 to 2015,which covers a whole solar cycle.This provides an opportunity to study the effect of the sunspot number on the generation of ECWs.For this purpose,in Figure 2 we show the occurrence rates, the medians of Vp, Vd/VA, Nα/Np, and the sunspot number with respect to months in the years between 2005 and 2015. The monthly sunspot numbers used in Figure 2(e) are from the Solar Influences Data analysis Center(SIDC),and are smoothed over 13 months.To make it easier to compare the values in different months, we have included reference lines in the graph.We divide the solar cycle into four periods I, II, III, IV in Figure 2, corresponding to less and decreasing sunspot number, less and increasing sunspot number, more and increasing sunspot number, more and decreasing sunspot number, respectively.

    Overall, two conclusions can be drawn from Figure 2. First,the occurrence rate and the medians of Vp, Vd/VA, Nα/Npassociated with LH ECWs are clearly correlated with the solar cycle, and their medians are usually greater than the medians associated with RH ECWs and the ambient plasma.Second,the generation of LH ECWs mainly depends on whether the sunspot number increases or decreases, rather than whether they are at the maximum or minimum. As in Figure 2(a), the occurrence rate of LH ECWs increases (decreases) as the sunspot number decreases(increases)corresponding to periods I and IV (II and III). The increased occurrence rate of LH ECWs is accompanied by an increase in the medians of Vp,Vd/VA, and Nα/Np. Looking closely at Figures 2(c) and (d),Vd/VAand Nα/Npseem to play major roles in the increase in occurrence rates during period I and period IV, respectively.For RH ECWs, the occurrence rate does not show a clear correlation with the solar cycle,and the medians of Vp,Vd/VA,Nα/Npare almost the same as the medians associated with the ambient plasma during most of the solar cycle.

    Specifically,inspecting the changes in the occurrence rates in Figure 2(a), it can be found that, during the strong influence periods 2005 January to 2008 August and 2014 October to 2015 December (periods I and IV in Figure 2) the occurrence rate of LH ECWs fluctuates considerably above 1.3% and is significantly higher than that of RH ECWs. Through careful comparison, it is found that the strong influence periods correspond to the decreasing of the sunspot number,rather than its minimum. While for RH ECWs, the occurrence rate fluctuates in a range from 0.5% to 1.0% for most of the solar cycle and does not show a significant dependence on the sunspot number.

    Figure 2(b) shows the medians of Vpassociated with LH ECWs,RH ECWs,and the ambient plasma.One can find three phenomena. First, during the solar cycle with less and decreasing sunspot number (period I), the medians of Vpassociated with LH ECWs,RH ECWs and the ambient plasma are all significantly larger than in other times. This can be partly explained by the formation mechanism of the fast solar wind. The fast solar wind with speed of up to 800 km s-1mainly originates from the coronal holes, which occur at the poles and low latitudes,reaching their largest when activity on the Sun is at its minimum(Lotova et al. 2000;Barabash 2012;Pierrard et al. 2020). However, during period II with a small and increasing sunspot number, the medians of Vpassociated with LH ECWs, RH ECWs and the ambient plasma show significant decrease. Second, throughout the solar cycle, the median of Vpassociated with LH ECWs is generally larger than that associated with RH ECWs and the ambient plasma,which is consistent with the preferential condition of higher Vpshown in Figure 1(b). Third, the median of Vpassociated with RH ECWs is almost the same as that associated with the ambient plasma, which is consistent with the small change in the occurrence rate of RH ECWs against Vpshown in Figure 1(b).

    It is believed that the presence of alpha particles plays an important role in kinetic instabilities which can induce plasma kinetic waves. Thus, in Figures 2(c) and (d), we give the evolutions of the medians of Vd/VAand Nα/Np. One can see that the large (small) occurrence rate of LH ECWs mostly corresponds to large (small) medians of Vd/VAand Nα/Npduring periods I and IV (II and III). One should note that near the right boundary of the period I (2008 August), the occurrence rate of LH ECWs shows a sudden decrease. At the same time, the median of Vd/VAreaches the maximum,while the median of Nα/Npgreatly decreases. This tends to indicate that the sudden decrease in the median of Nα/Npleads to the decrease in the occurrence rate of LH ECWs, and the occurrence rate should be understood in the combination of Vd/VAand Nα/Np.

    Figure 2.Occurrence rates(a),medians of Vp(b),Vd/VA(c),Nα/Np(d),and the sunspot number(e)with respect to months in the years between 2005 and 2015.The red,blue,and black solid lines are for LH ECWs,RH ECWs,and the ambient plasma,respectively.The black horizontal dashed lines indicate the occurrence rate of 1.3%, Vp of 500 km s-1, Vd/VA of 0.35, Nα/Np of 3.5%, and sunspot number of 50.

    Using the data from the STEREO mission orbiting the Sun just inside of 1 au, the occurrence rates of LH and RH ECWs with respect to months in the years between 2007 and 2013 were given in Figure 1 of Zhao et al. (2017). Compared with the sunspot number in Figure 2(e) of the present paper, the occurrence rate of LH ECWs from the STEREO mission(Zhao et al. 2017) increases during the period of 2007 January–2008 June when the sunspot number decreases and decreases during the period of 2008 July–2010 January when the sunspot number increases. This confirms the main conclusion of Figure 2 from the Wind mission in the present paper. One may note that the occurrence rates of LH ECWs in Zhao et al.(2017) and in the present paper show moderate differences especially after 2010, which might be attributed to the increasing distance between the STEREO mission and the Wind mission.The separation angle of the STEREO-A mission with Earth increased from 0.2°in 2007 January to 64°in 2010 January and then to 150° in 2013 December, while the Wind mission was placed in a halo orbit around Sun–Earth gravitational equilibrium point L1. For RH ECWs from the STEREO mission (Zhao et al. 2017), the occurrence rate is nearly unrelated to months, which is also consistent with the present paper.

    Figure 3. Occurrence rates of LH (red square line), RH (blue circle line)ECWs, and PDDs of the ambient plasma (black triangle line) with respect to the angle θBR (a) and θVdB (b).

    3.3. Dependence of Occurrence of ECWs on the Orientations of Vd and B

    To find out the effects of orientations of vectors Vdand B on the occurrence of ECWs in the solar wind,Figure 3 presents the occurrence rates of LH and RH ECWs, and the PDDs of the ambient plasma against θBRand θVdB. Taking θBRas an example, the angle between B and R (the radial vector of the Sun) is defined as

    with range[0,π].An angle θBR<90°means that B is pointing outward from the Sun while an angle θBR>90° denotes sunward orientation.In the same way,one can define θVdBand θVdRas the angle between Vdand B and between Vdand R.

    Inspecting Figure 3(a) where the PDD and occurrence rates with respect to θBRare shown, several points can be found as follows. First, the occurrence rates of LH and RH ECWs have the maximum when B is parallel or anti-parallel to R. Second,the occurrence rate of LH ECWs is larger than that of RH ECWs at most angles. Third, the occurrence rates of LH and RH ECWs fall close to 0 when θBRapproaches 90°. We may attribute the decrease in the occurrence rate around 90° to the observation limit. When θBRis around 90°, the spacecraft probably crosses approximately perpendicularly to the wavevector of ECWs and fails to detect changes in the magnetic field.

    Figure 3(b) displays the PDD and occurrence rates with respect to θVdB. According to the PDD of the ambient plasma,one can see that most of the data are distributed at 0°–20° and 160°–180° bins, indicating that Vdis usually quasi-(anti)parallel to B. It is easy to understand by considering that a component of Vdperpendicular to B would lead to a net force that would result in a gyrotropic distribution within several cyclotron periods (Kasper et al. 2006). From Figure 3(b), one can find that the occurrence rate of LH ECWs is more than twice that of RH ECWs, and have the maximum when Vdis parallel or anti-parallel to B.One point worth noting is that the occurrence rate of LH or RH ECWs at 0°–20° bin is nearly as large as that at 160°–180° bin, indicating that whether Vdis parallel or anti-parallel to B has no significant effect on the generation of ECWs.

    To simultaneously study the effects of the orientations of the Vdand B on the occurrence of ECWs, we present the occurrence rates associated with LH and RH ECWs in the space of (θBR, θVdB) in Figures 4(a) and (b). For LH ECWs in Figure 4(a), the occurrence rate depends on θVdBand θBRsimultaneously. The occurrence rate reaches the maximum when Vdand R are both quasi-(anti)parallel to B and Vdis quasi-parallel to R.However,for RH ECWs in Figure 4(b),the dependency is not very clear. It seems that when B is antisunward, RH ECWs are more likely to be generated at θVdB>90°.

    The above discussions show that the occurrence rates of ECWs are related to both θVdBand θBR. To include the effects of θVdBand θBRon the occurrence rates of ECWs simultaneously, we introduce a modified angle θVdB-R. If θBR>90°,we define θVdB-R=180°-θVdB. If θBR≤90°, we define θVdB-R=θVdB. Considering most of the data are distributed at θVdBbins of 0°–20°and 160°–180°,in Figure 4(c)we give the PDDs of the ambient plasma with respect to θVdRat 0°<θVdB-R<20° and 160°<θVdB-R<180°. It is found that θVdB-Rless than 20°(greater than 160°)mainly corresponds to θVdRless (greater) than 90°, indicating the anti-sunward(sunward) orientation of Vd.

    Data distributions regulated by the modified angle θVdB-Rare shown in Figure 4(d).First,one can see that the PDD at 0°–20°bin is significantly greater than that at 160°–180°bin.This is reasonable by considering that alpha particles often flow faster than protons in the solar wind (McKenzie et al. 1978).Second,the occurrence rate of LH(RH)ECWs at 0°–20°bin is much larger (slightly smaller) than that at 160°–180°bin. This indicates that the generation of LH (RH) ECWs has the strong(weak) preferential condition that Vdis anti-sunward(sunward), which mainly corresponds to θVdR<90°(θVdR>90°),although the preferential condition for generation of RH ECWs is not as obvious as that of LH ECWs. This is similar to the conclusions from Figures 4(a)and(b),except that only one parameter θVdB-Ris used. Third, the occurrence rate of LH ECWs at 0°–20°(160°–180°)bin is larger(smaller)than that of RH ECWs,and the dominant polarization of the ECWs changes from LH to RH when θVdB-Ris greater than 145°.

    Figure 4.Color scale plots of occurrence rates associated with LH(a)and RH(b)ECWs in the space of(θBR,θVdB).(c)PDDs of the ambient plasma with respect to θVdR at 0°<θVdB-R <20°(black solid line)and 160°<θVdB-R <180°(gray dashed line).(d)Occurrence rates of LH(red square line),RH(blue circle line)ECWs,and PDD of the ambient plasma (black triangle line) with respect to the modified angle θVdB-R.

    The effect of Vdorientation on the generation of ECWs in Figure 4 could be interpreted by the theoretical investigations of the kinetic instabilities using the linear Vlasov–Maxwell kinetic theory(Podesta&Gary 2011).The LH and RH ECWs are preferentially generated by electromagnetic ion cyclotron instability and parallel firehose instability with a maximum growth rate occurring at k·Vd>0 and k·Vd<0,respectively,where k is the wavevector (Podesta & Gary 2011). The theoretical result predicts that the main propagation direction of LH(RH)ECWs is the same as(opposite to)the direction of Vd.With this in mind,the conclusions drawn from Figure 4 can be interpreted. At 0°–20° (160°–180°) bin where Vdis antisunward (sunward), LH (RH) ECWs propagate away from the Sun in the plasma frame and the polarization does not change in the spacecraft frame. However, the RH (LH) ECWs propagate toward the Sun in the plasma frame and the polarization turns into LH (RH) in the spacecraft frame due to the presence of a large Doppler shift. Thus, the LH (RH)ECWs are the dominant waves when Vdis anti-sunward(sunward).

    4. Summary and Discussion

    In-situ measurements by the STEREO mission (Zhao et al.2017)initially give the preferential conditions for generation of LH ECWs. In this paper, refined measurements are carried out by the Wind mission in a halo orbit around the L1 Lagrange point. We obtained the obvious preferential conditions for generation of LH ECWs as higher Tp, higher Vp, lower Np,which confirm the findings obtained by the STEREO mission(Zhao et al. 2017). Furthermore, we identified the preferential conditions for generation of LH ECWs including stronger proton temperature anisotropy, higher Vd/VA, and higher Nα/Np. For RH ECWs, the occurrence rate presents little changes with the plasma parameters and the preferential conditions do not be determined in this paper.

    Measurements spanning 11 yr were carried out to study the effects of the solar cycle on the occurrence of ECWs.The solar cycle is divided into four periods according to whether the sunspot number is more or less and increases or decreases.It is found that the occurrence rate of the LH ECWs increases(decreases)as the sunspot number decreases (increases),rather than when it is more or less. The increased occurrence rate of LH ECWs is accompanied by an increase in the medians of Vp,Vd/VA, and Nα/Np. Specifically, Vd/VA(Nα/Np) tends to play a major role in the increase of the occurrence rate during the period with less(more)and decreasing sunspot number.This is consistent with the preferential conditions for generation of LH ECWs and can be partially interpreted by the changes in properties of the solar wind during different periods of the solar cycle. However, for the RH ECWs, the occurrence rate fluctuates with a smaller amplitude and the dependence on the solar cycle is not clear. The medians of Vp, Vd/VA, Nα/Npfluctuate nearly in the same manner as those of the ambient plasma.

    The orientations of B and Vdin the solar wind are variable,and it will be interesting to study their influences on the generation of ECWs.With respect to θBR,the occurrence rate of LH ECWs is larger than that of RH ECWs in the entire range,and they both reach the maximum when B is parallel or antiparallel to R. The occurrence rates of both LH and RH ECWs approach zero when θBRis around 90°, which could be attributed to observation limitations. As for the study of θVdB,we limit our measurements to Vd/Vp>1%,and the discussions in this paper mainly focus on bins of 0°–20° and 160°–180°,where most of the ambient plasma is distributed. The occurrence rates of ECWs in (θBR, θVdB) space show the occurrence rate of LH ECWs reaches the maximum when Vdand R are both quasi-(anti)parallel to B and Vdis quasi-parallel to R. The dependence of the occurrence rate of RH ECWs on θVdBin (θBR, θVdB) space is not obvious. By introducing a modified angle θVdB-R, it is shown that the generation of LH(RH) ECWs has the strong (weak) preferential condition that Vdis anti-sunward (sunward), and LH (RH) ECWs dominate when Vdis anti-sunward(sunward),which mainly corresponds to θVdR<90° (θVdR>90°).

    In summary, with Wind observations between 2005 and 2015,we study the dependences of the occurrence of ECWs on the plasma parameters, the solar cycle, and the orientations of Vdand B.It is found that LH ECWs are more likely to occur at higher Tp, Vp, Vd/VA, Nα/Np, lower Np, stronger proton temperature anisotropy, on the other hand, as the sunspot number decreases and Vdis anti-sunward. The findings could be of potential importance for a better understanding of the generation of ECWs. Before concluding, two remarks on the present study may be appropriate. The preferential conditions for generation of RH ECWs on the plasma parameters or on the solar cycle are ambiguous, at least not as clear as LH ECWs,which suggests that RH ECWs have a relatively complex generation mechanism. Perhaps more forms of free-energy should be considered in the future,such as differential flow and temperature anisotropy of proton beam, or those of alpha particles. Second, the dependence of the occurrence of ECWs on the solar cycle is not fully understood. For instance, the occurrence rate, and the medians of Vp, Vd/VAand Nα/Npassociated with LH ECWs are quite different during periods with increasing and decreasing sunspot number, although they both correspond to large or small sunspot number.

    Acknowledgments

    This research was supported by NSFC under grant Nos.41874204, 42174202, 41974197, 12103018. Research by G.-Q. Zhao was supported partly by the Excellent Young Scientists Fund in Henan Province under grant No.222300420061. The authors acknowledge the Wind mission for the data, which can be obtained via the Coordinated Data Analysis Web (http://cdaweb.gsfc.nasa.gov/cdaweb/istp_public/). The authors also acknowledge the SIDC for the sunspot number, which can be obtained from Web (https://wwwbis.sidc.be/silso/datafiles).

    ORCID iDs

    Guo-Qing Zhao, https://orcid.org/0000-0002-1831-1451

    听说在线观看完整版免费高清| 又爽又黄a免费视频| 男人添女人高潮全过程视频| 日韩av不卡免费在线播放| 日韩欧美精品免费久久| 日韩人妻高清精品专区| 国产毛片在线视频| 丝袜美腿在线中文| 久久久久久九九精品二区国产| 欧美性感艳星| 高清午夜精品一区二区三区| 久久久精品免费免费高清| 国产精品一区www在线观看| av在线天堂中文字幕| 黄片wwwwww| 国产男女超爽视频在线观看| 国产精品人妻久久久久久| 91精品伊人久久大香线蕉| 九色成人免费人妻av| 婷婷色综合大香蕉| 熟女av电影| 尤物成人国产欧美一区二区三区| 国产视频内射| 日本欧美国产在线视频| 久久人人爽人人爽人人片va| 精品人妻偷拍中文字幕| 国产精品伦人一区二区| 18禁裸乳无遮挡免费网站照片| 国产一区二区亚洲精品在线观看| 婷婷色综合大香蕉| 永久网站在线| 日本熟妇午夜| tube8黄色片| 极品教师在线视频| 国产淫片久久久久久久久| 18禁在线无遮挡免费观看视频| 久久久精品欧美日韩精品| 亚洲国产av新网站| 777米奇影视久久| 狂野欧美激情性bbbbbb| 久久精品久久精品一区二区三区| 久久久久久久午夜电影| 久久久国产一区二区| 少妇丰满av| 人妻夜夜爽99麻豆av| 国产淫语在线视频| 国产精品久久久久久精品电影小说 | 中文精品一卡2卡3卡4更新| 黑人高潮一二区| 亚州av有码| 亚洲人成网站在线观看播放| 在现免费观看毛片| 色吧在线观看| 久久久亚洲精品成人影院| 精品久久久久久久人妻蜜臀av| 成人二区视频| 国产黄片视频在线免费观看| 日韩制服骚丝袜av| 在线天堂最新版资源| 国产探花在线观看一区二区| 国产精品99久久99久久久不卡 | 在线免费观看不下载黄p国产| 高清av免费在线| 亚洲一区二区三区欧美精品 | 1000部很黄的大片| 夜夜爽夜夜爽视频| 建设人人有责人人尽责人人享有的 | 九九在线视频观看精品| 国产69精品久久久久777片| 国产精品女同一区二区软件| 国产精品成人在线| 久久久久久久久久成人| 国产精品一二三区在线看| 又大又黄又爽视频免费| 99视频精品全部免费 在线| 一本久久精品| 婷婷色综合大香蕉| 伊人久久国产一区二区| 亚洲综合精品二区| 亚洲欧洲日产国产| 亚洲av日韩在线播放| 熟女av电影| 国产中年淑女户外野战色| 丰满人妻一区二区三区视频av| av在线亚洲专区| 国产精品嫩草影院av在线观看| 97在线视频观看| 国产精品无大码| 久久久久久久久久人人人人人人| 极品少妇高潮喷水抽搐| 性插视频无遮挡在线免费观看| 99久久精品热视频| 老司机影院毛片| 免费av观看视频| 久久久久国产精品人妻一区二区| 亚洲精品国产成人久久av| 亚洲真实伦在线观看| 国产成人精品一,二区| 一级毛片我不卡| 亚洲一级一片aⅴ在线观看| 久久女婷五月综合色啪小说 | 久久久色成人| 99热国产这里只有精品6| 久久鲁丝午夜福利片| 国产精品久久久久久精品古装| 少妇裸体淫交视频免费看高清| av国产免费在线观看| 人体艺术视频欧美日本| 老师上课跳d突然被开到最大视频| 国产亚洲av嫩草精品影院| 99久久精品热视频| 插逼视频在线观看| 国产色爽女视频免费观看| 国产亚洲av嫩草精品影院| 国产黄片视频在线免费观看| av免费在线看不卡| 91午夜精品亚洲一区二区三区| 免费看不卡的av| 亚洲av不卡在线观看| 中文乱码字字幕精品一区二区三区| 中文精品一卡2卡3卡4更新| 黄色一级大片看看| 国产在线一区二区三区精| 国产在视频线精品| 丰满人妻一区二区三区视频av| 少妇高潮的动态图| 久久亚洲国产成人精品v| 精品人妻一区二区三区麻豆| 偷拍熟女少妇极品色| 亚洲一级一片aⅴ在线观看| 久久国产乱子免费精品| 中文天堂在线官网| 成人高潮视频无遮挡免费网站| 男女无遮挡免费网站观看| 亚洲成人一二三区av| 丝袜脚勾引网站| 性色avwww在线观看| 亚洲精品国产成人久久av| av又黄又爽大尺度在线免费看| 搡女人真爽免费视频火全软件| 中文天堂在线官网| 久久女婷五月综合色啪小说 | 国产日韩欧美在线精品| 国产精品蜜桃在线观看| 亚洲欧美日韩另类电影网站 | 欧美极品一区二区三区四区| 日韩欧美精品免费久久| 丰满乱子伦码专区| 亚洲av男天堂| 亚洲天堂av无毛| 欧美日韩一区二区视频在线观看视频在线 | 国产男女内射视频| 亚洲一区二区三区欧美精品 | 黄片wwwwww| 永久免费av网站大全| 99热这里只有是精品50| 色视频www国产| 一级av片app| 色婷婷久久久亚洲欧美| 久久97久久精品| 国产成人福利小说| 我的女老师完整版在线观看| 久久精品国产自在天天线| 久久久久久久久久久免费av| 人妻 亚洲 视频| 久久人人爽人人片av| 高清av免费在线| 久久99蜜桃精品久久| 亚洲欧美日韩东京热| 小蜜桃在线观看免费完整版高清| 小蜜桃在线观看免费完整版高清| 日韩一区二区视频免费看| 免费黄频网站在线观看国产| 国产精品久久久久久精品古装| 乱系列少妇在线播放| 好男人视频免费观看在线| 激情 狠狠 欧美| av又黄又爽大尺度在线免费看| av播播在线观看一区| av专区在线播放| 欧美zozozo另类| 街头女战士在线观看网站| 夜夜爽夜夜爽视频| 视频中文字幕在线观看| 亚洲av欧美aⅴ国产| 老女人水多毛片| 嫩草影院入口| 又大又黄又爽视频免费| 欧美激情久久久久久爽电影| 国内精品宾馆在线| 久久久久久久精品精品| 搡女人真爽免费视频火全软件| 国产乱人视频| 美女脱内裤让男人舔精品视频| 日韩,欧美,国产一区二区三区| 中国美白少妇内射xxxbb| 男女无遮挡免费网站观看| 欧美97在线视频| 欧美一级a爱片免费观看看| 精品久久久噜噜| 精品酒店卫生间| 91午夜精品亚洲一区二区三区| 一区二区av电影网| 欧美性感艳星| 亚洲精品久久午夜乱码| 国产精品熟女久久久久浪| 国产亚洲精品久久久com| 高清欧美精品videossex| 欧美激情在线99| 亚洲精品久久午夜乱码| 免费观看无遮挡的男女| 69人妻影院| 亚洲国产欧美人成| 18禁裸乳无遮挡动漫免费视频 | 中文字幕av成人在线电影| 夜夜看夜夜爽夜夜摸| 国产精品人妻久久久久久| 白带黄色成豆腐渣| av在线蜜桃| 男人爽女人下面视频在线观看| 两个人的视频大全免费| 肉色欧美久久久久久久蜜桃 | 在线看a的网站| 国产老妇女一区| 国产高潮美女av| 成年女人在线观看亚洲视频 | 国产高清三级在线| 超碰97精品在线观看| 建设人人有责人人尽责人人享有的 | 精品午夜福利在线看| 国产av码专区亚洲av| 国产精品国产av在线观看| 成人综合一区亚洲| 中文在线观看免费www的网站| 香蕉精品网在线| 日韩欧美精品免费久久| 中文乱码字字幕精品一区二区三区| 精品久久久精品久久久| 一二三四中文在线观看免费高清| 国产精品女同一区二区软件| 一区二区三区乱码不卡18| 超碰97精品在线观看| 97在线人人人人妻| 免费黄频网站在线观看国产| 丝袜喷水一区| 免费人成在线观看视频色| 日韩电影二区| 成人欧美大片| 精品久久久久久久人妻蜜臀av| 人体艺术视频欧美日本| 久久久久久久久久人人人人人人| 久久精品国产亚洲av涩爱| 99久久精品一区二区三区| 高清视频免费观看一区二区| 久久久久国产网址| 精品少妇黑人巨大在线播放| 国产一区有黄有色的免费视频| 嘟嘟电影网在线观看| 超碰97精品在线观看| 尤物成人国产欧美一区二区三区| 涩涩av久久男人的天堂| 亚洲人成网站高清观看| 午夜精品一区二区三区免费看| 亚洲人成网站在线观看播放| 亚洲最大成人中文| 看免费成人av毛片| 99久久九九国产精品国产免费| 一级毛片黄色毛片免费观看视频| 三级国产精品欧美在线观看| 精品国产一区二区三区久久久樱花 | 亚洲av电影在线观看一区二区三区 | 69av精品久久久久久| 舔av片在线| 一级二级三级毛片免费看| 亚洲av免费在线观看| 亚洲最大成人av| 亚洲,一卡二卡三卡| 美女内射精品一级片tv| 久久韩国三级中文字幕| 国产免费一区二区三区四区乱码| 亚洲国产色片| 麻豆精品久久久久久蜜桃| 国产一级毛片在线| 亚洲不卡免费看| 国产亚洲av片在线观看秒播厂| 成人午夜精彩视频在线观看| 久久精品国产a三级三级三级| 国产成人精品一,二区| 成年女人在线观看亚洲视频 | 看十八女毛片水多多多| 久久女婷五月综合色啪小说 | 99热这里只有是精品50| 日韩 亚洲 欧美在线| 午夜激情福利司机影院| 18禁动态无遮挡网站| 亚洲天堂国产精品一区在线| 国产 一区 欧美 日韩| 亚洲人成网站在线观看播放| 青春草视频在线免费观看| 老女人水多毛片| 国产探花在线观看一区二区| 久久精品熟女亚洲av麻豆精品| 尾随美女入室| 十八禁网站网址无遮挡 | 少妇裸体淫交视频免费看高清| av在线播放精品| 国产在视频线精品| 毛片女人毛片| 免费黄色在线免费观看| 免费黄频网站在线观看国产| 亚洲无线观看免费| 亚洲一区二区三区欧美精品 | 婷婷色麻豆天堂久久| 国产精品成人在线| 国产女主播在线喷水免费视频网站| 在线观看av片永久免费下载| 亚洲成人中文字幕在线播放| 国产成年人精品一区二区| 哪个播放器可以免费观看大片| 亚州av有码| 精品久久久久久久久av| 亚洲精品国产色婷婷电影| 王馨瑶露胸无遮挡在线观看| 成年人午夜在线观看视频| 嘟嘟电影网在线观看| 夜夜爽夜夜爽视频| 日本熟妇午夜| 高清午夜精品一区二区三区| 美女脱内裤让男人舔精品视频| 日本免费在线观看一区| 我的老师免费观看完整版| 国产高潮美女av| 日产精品乱码卡一卡2卡三| 少妇人妻 视频| 中文精品一卡2卡3卡4更新| 久久综合国产亚洲精品| 精品少妇久久久久久888优播| 久久97久久精品| 午夜福利视频1000在线观看| 狂野欧美激情性bbbbbb| 最近2019中文字幕mv第一页| 亚洲精品国产av蜜桃| 天堂网av新在线| 亚洲精华国产精华液的使用体验| 交换朋友夫妻互换小说| 少妇被粗大猛烈的视频| 国产精品一及| 有码 亚洲区| 亚洲精品国产色婷婷电影| 97在线视频观看| 国产亚洲午夜精品一区二区久久 | 我要看日韩黄色一级片| 国产黄片视频在线免费观看| 欧美成人a在线观看| 国产精品久久久久久精品电影小说 | 中文精品一卡2卡3卡4更新| 在线观看美女被高潮喷水网站| 女人十人毛片免费观看3o分钟| 久久久久九九精品影院| 亚洲最大成人手机在线| 国产中年淑女户外野战色| 另类亚洲欧美激情| 一个人看的www免费观看视频| 99视频精品全部免费 在线| 91午夜精品亚洲一区二区三区| 我的女老师完整版在线观看| 春色校园在线视频观看| 超碰97精品在线观看| 日本三级黄在线观看| 欧美成人精品欧美一级黄| 精品人妻一区二区三区麻豆| 欧美日韩一区二区视频在线观看视频在线 | 免费不卡的大黄色大毛片视频在线观看| 国国产精品蜜臀av免费| 91精品伊人久久大香线蕉| 美女高潮的动态| 国产午夜精品一二区理论片| 成年版毛片免费区| 久久精品人妻少妇| 亚洲最大成人av| 日韩大片免费观看网站| 国产成人福利小说| 国产综合精华液| 卡戴珊不雅视频在线播放| 最近的中文字幕免费完整| 男女那种视频在线观看| 欧美日本视频| 成人亚洲欧美一区二区av| 黄片无遮挡物在线观看| 亚洲国产欧美人成| av国产精品久久久久影院| 五月天丁香电影| 成人鲁丝片一二三区免费| 内地一区二区视频在线| 国产伦在线观看视频一区| av在线蜜桃| 日韩成人av中文字幕在线观看| 丰满少妇做爰视频| av女优亚洲男人天堂| 老女人水多毛片| 免费av不卡在线播放| 街头女战士在线观看网站| 亚洲自偷自拍三级| 啦啦啦中文免费视频观看日本| 国产亚洲精品久久久com| 女的被弄到高潮叫床怎么办| 国产探花在线观看一区二区| 白带黄色成豆腐渣| 男的添女的下面高潮视频| 久久久成人免费电影| 久久韩国三级中文字幕| 青春草国产在线视频| av专区在线播放| 老司机影院毛片| 亚洲成人精品中文字幕电影| 男女边吃奶边做爰视频| 国产成人a∨麻豆精品| 成年人午夜在线观看视频| 别揉我奶头 嗯啊视频| 18+在线观看网站| 身体一侧抽搐| 亚洲精品日韩av片在线观看| 日韩精品有码人妻一区| 久久ye,这里只有精品| 日韩欧美精品免费久久| 亚洲国产成人一精品久久久| 精品视频人人做人人爽| 91精品一卡2卡3卡4卡| 五月开心婷婷网| 久久99精品国语久久久| a级毛片免费高清观看在线播放| 亚洲最大成人手机在线| 国产精品久久久久久精品古装| 大又大粗又爽又黄少妇毛片口| 女人久久www免费人成看片| 久久精品久久久久久噜噜老黄| 亚洲av电影在线观看一区二区三区 | 99久国产av精品国产电影| 欧美成人精品欧美一级黄| 91狼人影院| 一区二区三区精品91| 亚洲内射少妇av| 国产精品嫩草影院av在线观看| 午夜免费观看性视频| eeuss影院久久| 国产免费一区二区三区四区乱码| 九九爱精品视频在线观看| 美女cb高潮喷水在线观看| 人妻一区二区av| 最新中文字幕久久久久| 亚洲欧美一区二区三区国产| 成人亚洲精品一区在线观看 | 亚洲色图av天堂| 在线观看三级黄色| 国产精品久久久久久久久免| 1000部很黄的大片| 综合色丁香网| 在线观看一区二区三区激情| 国产极品天堂在线| 一区二区三区乱码不卡18| 午夜视频国产福利| 国产精品秋霞免费鲁丝片| 免费黄色在线免费观看| 久久久久久久久久久免费av| 亚洲欧洲日产国产| 极品教师在线视频| 国内少妇人妻偷人精品xxx网站| 嫩草影院入口| 麻豆成人av视频| 国产欧美另类精品又又久久亚洲欧美| 国产精品久久久久久精品古装| 91久久精品电影网| 国产黄片视频在线免费观看| 国产成人a区在线观看| 亚洲av成人精品一二三区| 久久人人爽人人片av| 国产av码专区亚洲av| 我要看日韩黄色一级片| 亚洲图色成人| 男人舔奶头视频| 国内精品宾馆在线| 国产精品.久久久| 一级a做视频免费观看| 亚洲av在线观看美女高潮| 国产乱人视频| 婷婷色综合www| 看免费成人av毛片| 一级二级三级毛片免费看| 日韩视频在线欧美| 精华霜和精华液先用哪个| 亚洲欧美一区二区三区国产| 夜夜看夜夜爽夜夜摸| 亚洲,一卡二卡三卡| 人人妻人人看人人澡| 国产成人免费观看mmmm| 91精品国产九色| 男的添女的下面高潮视频| 国语对白做爰xxxⅹ性视频网站| 久久鲁丝午夜福利片| 高清av免费在线| 精品人妻一区二区三区麻豆| 免费看a级黄色片| 中文乱码字字幕精品一区二区三区| 国产精品一区二区在线观看99| 97在线视频观看| 高清在线视频一区二区三区| 国产极品天堂在线| 五月天丁香电影| 你懂的网址亚洲精品在线观看| 女人十人毛片免费观看3o分钟| 国产精品不卡视频一区二区| 久久人人爽人人爽人人片va| 亚洲不卡免费看| 婷婷色综合大香蕉| 久久精品久久久久久久性| 黄色配什么色好看| 嫩草影院新地址| 亚洲欧美日韩东京热| 国产高清有码在线观看视频| 看免费成人av毛片| 在线观看免费高清a一片| 伦理电影大哥的女人| 亚洲国产精品成人综合色| 女人久久www免费人成看片| 男男h啪啪无遮挡| 亚洲成色77777| 欧美xxxx黑人xx丫x性爽| 大陆偷拍与自拍| 久久久久九九精品影院| 天堂中文最新版在线下载 | 国产69精品久久久久777片| av在线观看视频网站免费| 国产精品一二三区在线看| 自拍偷自拍亚洲精品老妇| 婷婷色av中文字幕| 少妇人妻一区二区三区视频| 欧美日韩亚洲高清精品| 视频区图区小说| 精品人妻一区二区三区麻豆| 日韩三级伦理在线观看| 男人添女人高潮全过程视频| 久久久国产一区二区| 久热这里只有精品99| 青春草亚洲视频在线观看| 大陆偷拍与自拍| 久久久久性生活片| 最近最新中文字幕免费大全7| 亚洲精华国产精华液的使用体验| 性色avwww在线观看| 自拍偷自拍亚洲精品老妇| 真实男女啪啪啪动态图| 视频中文字幕在线观看| 日本-黄色视频高清免费观看| 性色av一级| xxx大片免费视频| 看免费成人av毛片| 国产精品99久久99久久久不卡 | 波多野结衣巨乳人妻| 国产av国产精品国产| 亚洲av欧美aⅴ国产| 亚洲精品乱码久久久v下载方式| 男女国产视频网站| 99视频精品全部免费 在线| 少妇熟女欧美另类| 我的女老师完整版在线观看| 亚洲在久久综合| 99久国产av精品国产电影| 美女主播在线视频| 亚洲av中文字字幕乱码综合| 成人无遮挡网站| 国产精品福利在线免费观看| 真实男女啪啪啪动态图| 久久精品久久精品一区二区三区| 1000部很黄的大片| 国产久久久一区二区三区| 我要看日韩黄色一级片| 日本-黄色视频高清免费观看| 特大巨黑吊av在线直播| 欧美另类一区| 国产高清不卡午夜福利| 欧美丝袜亚洲另类| 久久久成人免费电影| 91狼人影院| 国产老妇伦熟女老妇高清| 成人国产麻豆网| 亚洲欧美成人精品一区二区| 国内精品宾馆在线| 免费播放大片免费观看视频在线观看| 午夜免费男女啪啪视频观看| 高清欧美精品videossex| 日日啪夜夜撸| 亚洲不卡免费看| 欧美性猛交╳xxx乱大交人| 身体一侧抽搐| 亚洲不卡免费看| 99九九线精品视频在线观看视频| 三级国产精品欧美在线观看| 亚洲久久久久久中文字幕| 舔av片在线| 身体一侧抽搐| 亚洲图色成人| 视频中文字幕在线观看| 国产乱人视频| 又黄又爽又刺激的免费视频.| 国产黄片视频在线免费观看| 欧美精品一区二区大全| 成年av动漫网址| 舔av片在线| 我的老师免费观看完整版| 免费黄网站久久成人精品| 搡老乐熟女国产| 我的老师免费观看完整版| 伊人久久国产一区二区| 久久99热这里只频精品6学生| 国产精品一区二区性色av|