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

    Temperature Analysis of Flaring (AR 11283) and Non-flaring (AR 12194)Coronal Loops

    2022-05-24 14:21:08NargesFathalianSeyedehSomayehHosseiniRadNasibehAlipourandHosseinSafari

    Narges Fathalian ,Seyedeh Somayeh Hosseini Rad ,Nasibeh Alipour ,and Hossein Safari

    1 Department of Physics,Payame Noor University (PNU),P.O.Box 19395-3697,Tehran,Iran;narges.fathalian@gmail.com

    2 Department of Physics,Faculty of Science,University of Zanjan,P.O.Box 45195-313,Zanjan,Iran

    Abstract Here,we study the temperature structure of flaring and non-flaring coronal loops,using extracted loops from images taken in six extreme ultraviolet channels recorded by Atmospheric Imaging Assembly/Solar Dynamics Observatory.We use data for loops of an X2.1-class-flaring active region(AR 11283)during 22:10 UT until 23:00 UT,on 2011 September 6;and a non-flaring active region(AR 12194)during 08:00:00 UT until 09:00:00 UT on 2014 October 26.By using the spatially synthesized Gaussian differential emission measure(DEM)forward-fitting method,we calculate the peak temperatures for each strip of the loops.We apply the Lomb–Scargle method to compute the oscillation periods for the temperature series of each strip.The periods of the temperature oscillations for the flaring loops ranged from 7 to 28.4 minutes.These temperature oscillations show very close behavior to the slow-mode oscillation.We observe that the temperature oscillations in the flaring loops started at least around 10 minutes before the transverse oscillations and continue for a long time duration even after the transverse oscillations ended.The temperature amplitudes increased during the flaring time(20 minutes)in the flaring loops.The periods of the temperatures obtained for the non-flaring loops ranged from 8.5 to 30 minutes,but their significances are less(below 0.5)in comparison with the flaring ones(near to one).Hence the detected temperature periods for the non-flaring loops’strips are less probable in comparison with the flaring ones,and maybe they are just fluctuations.Based on our confined observations,it seems that the flaring loops’periods show more diversity and their temperatures have wider ranges of variation than the non-flaring ones.More accurate commentary in this respect requires more extensive statistical research and broader observations.

    Key words:Sun:flares–Sun:corona–Sun:oscillations–stars:imaging–stars:solar-type

    1.Introduction

    Analyzing the thermal structure of coronal loops is of considerable interest,especially as these magnetic loops have an essential role in heating the solar chromosphere and corona.Such analysis can help to describe how the process of solar flaring is correlated with the loop’s thermal structure.

    Detections of coronal waves have a historical preview and have been reported several times(e.g.,Aschwanden et al.1999;Berghmans&Clette 1999;Nakariakov et al.1999;De Moortel et al.2000;Wang et al.2003;Verwichte et al.2004;Wang &Solanki 2004;De Moortel &Brady 2007;Ballai et al.2011).Coronal seismology and MHD waves have been reviewed widely by De Moortel (2005),Nakariakov &Verwichte(2005),Aschwanden (2006),Banerjee et al.(2007) and De Moortel&Nakariakov(2012).Along with the development of the observations,transverse and longitudinal oscillations have also been studied theoretically (e.g.,Gruszecki et al.2006;Pascoe et al.2007;Fathalian &Safari 2010;Fathalian et al.2010;Luna et al.2010).Coronal seismology techniques help to elicit the information from observations of oscillatory phenomena and the results can be interpreted by using theoretical models(see for e.g.,Roberts et al.1984;Goossens et al.1992).Oscillatory patterns and processes which happen during solar flares have been interesting and the subject of investigations from different approaches (e.g.,Nakariakov et al.2010;Anfinogentov et al.2013;Nisticò et al.2013;Hindman &Jain 2014;Russell et al.2015).As we know,transverse loop oscillations usually occur in response to a close filament or flare(Wills-Davey &Thompson 1999).

    Rapidly decaying long-period oscillations are mostly interpreted as global (or fundamental mode) standing slow magnetoacoustic waves (reviewed by Liu &Ofman 2014,and Wang 2011,also see Ofman &Wang 2002,and for slowmode observed in fan-loops see Pant et al.2017).They often occur in hot coronal loops of active regions (ARs),associated with tiny (or micro) flares.Increasing evidence has suggested that the harmonic type of decaying pulsations detected in intensity plots of solar and stellar flares are possibly caused by standing slow-mode waves (see reviews by Van Doorsselaere et al.2016 and McLaughlin et al.2018).Excitation,propagation and damping mechanisms of slow-mode waves have been studied theoretically (e.g.,Wang et al.2007,2015;Jess et al.2016;Nakariakov et al.2017;Nisticò et al.2017;Kolotkov et al.2019;Krishna Prasad et al.2019;Reale et al.2019;Wang&Ofman 2019).To have a complete overview of slow-mode magnetoacoustic waves in coronal loops see the review by Wang et al.(2021).

    Investigating and comparing the thermal structures and oscillations of coronal loops in loops of flaring and non-flaring ARs could help us in better understanding the loops’ material oscillations and the flare impact on them.Several different methods have been developed to investigate the thermal structure of the coronal loops and loop strands.The thermal stability of the coronal loops was the subject of research done by Habbal &Rosner (1979) (and references cited therein).McClymont &Craig (1985) stated that a pressure fluctuation must assist asymmetric coronal temperature perturbation.They concluded that coronal loops are fairly stable in the case of uniform heating.Van Doorsselaere et al.(2011) used spectroscopic line ratios to obtain the required temperature (via CHIANTI code)and estimated the adiabatic index of the corona.The dependence of coronal loop temperature on loop length and magnetic field strength is also a favorite topic.For instance,Dahlburg et al.(2018)probed the temperature properties of solar coronal loops over a wide range of lengths and magnetic field strengths via numerical simulations and observed a very high correlation between magnetic field strength and a maximum of the temperature.The effect of temperature inhomogeneity on the periods and the damping times of the standing slow-modes in stratified solar coronal loops was also studied(e.g.,Abedini et al.2012).Fathalian (2019) estimated the loop temperature using the intensity ratios and the Atmospheric Imaging Assembly(AIA) response functions in different wavelengths.Different emission measure computations and methods have been developed to estimate the temperature in the corona,which led to various discussions.Schmelz et al.(2010)analyzed a coronal loop,which was observed on 2010 August 3 by AIA.They took some differential emission measure (DEM) curves,claiming a multithermal rather than an isothermal DEM distribution(for the cross-sectional temperature of the loop).After that,Aschwanden&Boerner (2011) criticized the method of background subtraction which Schmelz et al.had applied.They claimed that the background subtraction method caused their inferred result of a multithermal loop.Aschwanden &Boerner (2011)analyzed a set of a hundred loops and concluded that 66%of the loops could be fitted with a narrowband single-Gaussian DEM model.In this regard,some attention was paid to the instrumental limitations and ability of AIA and Guennou et al.(2012a,2012b) discussed the accuracy of the DEM diagnostics of solar plasmas with respect to the AIA instrument of Solar Dynamics Observatory (SDO).The abovementioned controversy of whether the cross-field temperatures of coronal loops are multithermal or isothermal continued with Schmelz et al.(2013) (similar to Schmelz et al.2011).They analyzed twelve loops to investigate their cross-field temperature distributions and revealed the loops’ substructure.Based on their achievements,the warmer loops entail broader DEMs.Thereafter,Schmelz et al.(2014)found indications of a relationship between the DEM weighted-temperature and the cross-field DEM width for coronal loops.They argued that cooler loops tend to have narrower DEM widths.This could imply that fewer strands are seen emitting in the later cooling phase,which they claim could potentially resolve the abovementioned controversy.In this subject,Aschwanden et al.(2015)(as well as Aschwanden et al.2013)developed a method to extract the loop temperature which is based on a Gaussian fit for the DEM,named spatially synthesized Gaussian DEM forward-fitting method.

    This paper aims to analyze and compare thermal oscillations of coronal loops in flaring and non-flaring ARs,11283 and 12194,respectively.The contents of this paper are as follows:In Section 2,we introduce the considered flaring and nonflaring ARs and describe the data employed and the time and properties of the flare which occurred in the AR.In Section 3,we explain the method we use to analyze the time-series of temperatures in different strips of the loops.Section 4 reports our results obtained related to flaring and non-flaring regions.In Section 5,we briefly state a summary of this work.

    2.Data

    We investigate the thermal structure and treatment of loops in a flaring region to see if it follows the transverse oscillations of the loops,and we examine the thermal fluctuations at the flare time.For this purpose,we select a high energy flare x2.1 for which the transverse oscillations of two loops have been analyzed by Jain et al.(2015).They analyzed intensity variations in the wavelength 171 ? in two coronal loops of this region and detected obvious transverse oscillation with periods of roughly 2 minutes and decay times of 5 minutes for these loops at the flare time.To see the specific thermal properties of the flaring loops,as a blind test,we select a nonflaring AR,extract its loops and analyze their thermal treatment.Then we compare the temperature treatment of the loops at the flaring region with the loops of the non-flaring region to see the differences.

    The temperature analysis done here uses extreme ultraviolet(EUV) images from the AIA onboard the SDO.AIA has ten different wavelength channels,three in white light and ultraviolet(UV),and the other seven in EUV channels.Among these seven,the 304 filter,which is mostly sensitive to chromospheric temperatures (in order of T=104.7K),not the corona,is not taken into account (Aschwanden et al.2015).Therefore,we consider the images of the events in the six wavelengths (94,131,171,193,211,335 ?).These are covering the coronal temperature range from T ≈0.6 to T ≥16 MK.

    The two below data sets are finally selected to study thermal variations and coronal loop oscillations in flaring or non-flaring ARs.A few distinct loops are visible in the regions.Finally,these loops are chosen:

    1.Three loops of the x-flaring AR 11283:Observationally,X-class flares rarely happen around the loops with the specification we are looking for.So,this selected line-ofsight(LOS)X-flare,which occurs near the loops,is a rare case.We consider EUV images of NOAA AR 11283,in the time period of 22:10 UT until 23:00 UT on 2011 September 6 with the cadence of 12 s.This period of time is selected since no other flare happens then.A few distinct loops are visible and follow-able here during this period.Loop shapes in our AR change permanently;therefore,it is difficult or impossible to follow a loop over a very long time.Hence,it is not useful to extend the time interval of this region to the time before the flare.The transverse oscillations of two loops in this region were analyzed before by Jain et al.(2015).We mark these loops by A and B in Figure 1(b).They detected fundamental mode oscillation with periods of roughly 2 minutes and decay time of 5 minutes for these loops.We are curious to see the loops’ thermal oscillations (if any)or thermal fluctuations in this condition.Figure 1(a)(left) displays AR 11283 and the area,indicated by the white box,is featured in a zoom-in view in Figure 1(b)(right)and the five selected parts of the center of the three chosen loops are shown by red lines (a movie of the region is available in this link).As is clear in the movie,these three loops oscillate together and their oscillations decay simultaneously.The center of Figure 1(a) has coordinates at (230,165) arcsec and its width and height are 450″×456″/ 750×775 pixels.The flare occurring in this AR is an X2.1 class flare located close to the disk center at latitude 14° north and longitude 18° westThis flare initiates at 22:12 UT,ends about 22:24 UT with the peak at 22:20 UT and is associated with a coronal mass ejection (CME) which occurs from 2011 September 6,21:36:05T to 2011 September 7,02:24:05T,with the radial velocity of 469 km s?1,angular width of 252°and position angle of 275°(for more details look at the LASCO CME catalog).3Based on data on these websites:https://solarflare.njit.edu/webapp.html and https://www.swpc.noaa.gov/.

    Figure 1.(a)AIA image of the AR 11283 on 2011 September 6,22:10 UT as seen in the 171 ? filter.(b)Zoom-in view of the area marked by a box is on the left.The selected loops are distinguished in red.The loops A and B are the same loops studied by Jain et al.(2015) (see Figure 3(a) in Jain et al.2015).

    2.Three loops of non-flaring AR 12194:As a blind test,we select three loops of the non-flaring (nonf hereafter) AR 12194 in the smooth time period of 08:00:00 UT until 09:00:00 UT of 2014 October 26.The center of Figure 2(a) has coordinates at (0,?264) arcsec and its width and height are 615″×615″/1025×1025 pixels.We consider the images of the selected area with the cadence of 12 s in the same six wavelengths mentioned above.These loops are relatively motionless and do not show any transverse oscillation(see the region’s movie in the link).We select the loops in such a way that they do not have any crossing over the neighboring loops (from our perspective)during this time.In Figure 2 the selected loops are distinguished in red in the mentioned AR.The size of the final cut of the nonf region (presented on the right) is 351×401 pixels.

    Figure 2.(a)The NOAA AR 12194 on 2014 October 26,at 08:00:00 UT in 171 ? recorded by AIA/SDO.(b)Zoom-in view of the area,marked by a box in the left,the loops are distinguished in red.

    Figure 3.From top to bottom:The time-series of the temperature oscillations for the first three strips of Loop A(strips 1 to 3 from top to bottom,respectively),and the first two strips of Loop B1.The horizontal axis is time and the vertical axis is the logarithm of temperature.The red lines mark the initial and final time of the flare x2.1.

    The data set is primarily downloaded at level 1 with a pixel resolution of 0.6″.We use the standard aia_prep.pro subroutine available in the SDO package SolarSoftWare(SSW)library to adjust the screen scale between the four arms of the AIA.This pre-processing step increases the data level from 1 to 1.5,so that finally no jump or sudden movement is observed in the image series.We also utilize the drot_map.pro subroutine to correct the differential rotation effect.According to the movie made by pre-processed images,the most obvious loops (marked in the abovementioned figures) are selected in each region(with obvious transverse oscillations in the case of the flaring AR).

    3.Temperature Analysis Method

    We extract the selected loop segment pixels,for each loop,and calculate the normal vectors to each point of the loop’s direction.Then by using these data,we straighten each loop in a considered box with the thickness of 15–40 pixels (macropixels,depending on the available empty area around each loop and the distance to the neighboring loop).The area around the loop is needed for calculations of background subtraction.The selected loop segment is cut in all wavelengths and at the same considered box from the set of images.These loop images are necessary prerequisites for our thermal analysis process.Then the loop is divided into different strips and its best division in terms of pixel intervals is considered.To do thermal analysis,we use the spatially synthesized Gaussian DEM forward-fitting method proposed by Aschwanden et al.(2015).

    The images in the above six wavelength filters are considered to calculate the temperature in each strip of the loop.The DEM function is considered a single-Gaussian function relative to the temperature determined by the forward fitting method.To obtain the temperature for each loop,we divided the loop into narrow strips,and then the intensity flux was averaged over each strip.The number of each strip is displayed with the index i.One of the usual methods to subtract the background from observed data is fitting a single-Gaussian cospatial function with a linear function on the flux profile.The DEM for each strip is considered to be a single-Gaussian DEM in terms of the logarithm of the temperature,which has three free parameters (Aschwanden &Boerner 2011)

    In this expression,Tp,iis the DEM peak temperature,EMp,iis the peak EM function and σT,iis the logarithmic width of the temperature for that strip.To calculate the backgroundsubtracted fluxes (for each strip) we use Equation (6) of Aschwanden &Boerner (2011) (below)

    Here,Rλ(T) is the instrumental temperature response function of each wavelength filter λ,which is obtained by the code aia_get_response.pro in the SSW package.As time has passed,the calibration of the AIA response functions has partly changed.Here,we use the updated calibration of the temperature response functions,for each of the AIA temperature filters,according to the CHIANTI Version 2019 code available in SSW.After forward-fitting the Gaussian DEM to the background-subtracted observed fluxes in multiple wavelengths,the three-fitting parameters,temperature width (σT,i),peak of temperature (Tp,i) and peak emission measure (EMp,i),are found by minimizing

    Our data sample is uneven because of omitting some damaged images in between.Therefore to analyze the temperature oscillations,we use the Lomb–Scargle method.This method implements the periodogram technique,in thecase where the observation times are unevenly spaced(Scargle 1982).The Lomb–Scargle periodogram method is useful in cases where the periodicity of data treatment is not immediately apparent.This method allows efficient computation of a Fourier-like power spectrum estimator from unevenly sampled data,resulting in an intuitive means of determining the period of oscillation (VanderPlas 2018).Therefore we use the Lomb–Scargle Periodogram to efficiently evaluate and estimate the periods of temperature oscillations in our loops.We select the first period related to the highest power frequency,which isobtained by this method.We consider the achieved periods with the highest significances and amplitudes.The most significant (highest) periods observed in temperature (minute)for flaring and nonf loops are listed in Tables 1 and 2,respectively.To estimate the significance of the periods,we computed the probability values (p-values).In the Lomb–Scargle method,the significance returned here is the false alarm probability of the null hypothesis,i.e.,as the data are composed of independent Gaussian random variables.Accordingly,low probability values (p-value less than 0.05) indicate a high degree of significance in the associated periodic signal.

    Table 1 The Properties Observed for the Loop Segments of the Flaring AR

    Table 2 The Properties Observed for the Loop Segments of the Non Flaring AR

    4.Results

    4.1.Temperature Analysis of Flaring Active Region Loops

    Thenceforth,the temperature time-series of different strips of the selected loops are calculated using the method described in Section 3.In the following figures,the vertical axis shows the logarithm of the temperature and the horizontal axis displays the time duration.To be comparable by eyes,all the forthcoming figures (which feature temperature oscillations of the loops)have been co-scaled in the(log)temperature range of 5.7–6.9.The color maps are shown for each temperature map.Loops A,B1,B2,C1 and C2 are subdivided into 25,11,8,12 and 6 strips,respectively.Each strip’s length is equal to 4 pixels (macro-pixel) for all loops in this paper.For brevity,a few strips’temperature oscillations are presented here.Figure 3 displays the time-series of temperature oscillations for the first three strips of Loop A,and the first two strips of Loop B1.We calculated the errors for each point (temperature) but removed them in the presentation to avoid overcrowding the figures.As we observe in Figures 3 and 4,the temperature oscillations started and increase around 22:12 before the flare peak time(22:20)and are mostly continuing after the flare ended(22:24).These temperature oscillations follow the transverse loop oscillations observed by Jain et al.(2015).As Jain et al.reported,Loops A and B have transverse oscillations with periods of roughly 2 minutes and decay times of 5 minutes,starting at 22:18 around the flare peak time (23:20) and decaying after the flare ended (22:24).So as we observe,the temperature oscillations in these flaring loops happen before the start of their transverse oscillations and are continuing even in the time interval after the transverse oscillations decay,although the temperature oscillations do not decay as rapidly as the transverse oscillations do,and conversely,the loop temperature increases at the end of the oscillating mode (see Figure 4,the temperature map of Loop A,for instance).

    Figure 4.Temperature map of the flaring Loops A,B1,B2,C1 and C2(from top to bottom respectively)as a time series.The vertical axis is distance along the loop in Mm,and the horizontal axis is time.The colorbar on the left shows the colors considered for the temperature range.

    We calculate the temperature oscillation periods using the Lomb–Scargle method.We consider the thermal oscillation periods with the highest significances.As this method shows,the most powerful period in the range of data time-series(listed in Table 1) is from 7 to 28.4 minutes observed in the strips of the marked loops of this flaring region.These loops of the flaring region also show some short periods in temperature oscillations in which some are less than 10 minutes(as listed in Table 1).These short periods are more frequently observed in the loops of the flaring AR.Such short periods are very scarce for the loops of the nonf AR (compare Tables 1 and 2).

    The first column in Table 1 is the number of every strip along the loop.The second column is the period of the most powerful frequency observed for the loop strips,calculated by the Lomb–Scargle method.The third column shows the maximum of log(T) minus its minimum in each strip.The columns of Table 2 are exactly the same as those in Table 1;the only difference is that Table 2 is for the nonf loops.

    Loop A has a length of 42.3(Mm)which is the length of the selected part of the loop marked in Figure 1(b).The mean of the parameter (Max(logT)–Min(logT))for the strips of Loop A is 1.21.The mean of the temperature (log) of this loop over time is 6.15±0.25.Loop B1,divided into 11 strips,has a length of 20.24 (Mm).The mean of (Max(logT)–Min(logT))and the mean of temperature for this loop are 1.10,and 6.28±0.22 respectively.Loop B2,which has eight strips,with a length of 15.61 (Mm),has a mean temperature (log) of 6.21±0.21.The mean of (Max(logT)–Min(logT)) is 0.81 through this loop segment.Loops C1 and C2,divided into 12 and 6 strips,have lengths of 22.08 and 11.06 (Mm),mean temperatures of 6.25±0.22 and 6.14±0.25 (log),and means(Max(logT)–Min(logT)) of 1.48 and 0.88,respectively.

    We observe that despite the temperature oscillations,the flaring loops show a temperature rise at the end of the considered time interval (Figure 3).As their temperature maps also affirm,the oscillations follow with a relatively sensible rise in the final temperature of the loop segments (Figure 4),although in the case of transverse oscillations,the loops oscillate as the flare occurs and then the oscillations decay and stop.In the case of temperature oscillations,the temperatures of the various strips of the loops oscillate and,at the end of the flare occurrence,they arrive at a relatively higher value of temperature on average.

    Figure 4 depicts the temperature maps of flaring loops A,B1,B2,C1 and C2,from top to bottom respectively,as a time series.In each plot,the vertical axis is the distance along the loop segment in Mm and the horizontal axis represents time.The color bar(on the left)signifies the temperature range.Each separate grid part on the map corresponds to one strip.Figure 4 demonstrates that the temperature for most of the strips increased,bypassing a few oscillations.Before the end of the time duration,some strips become hotter (yellow ones) and some cooler(blue ones).Loop B1 is colder at the early times of the duration and becomes hotter at the middle and end times with a swing to lower temperatures again(see Figure 4).There are some temperature fluctuations at the middle times (the red and green stripes) while at the end the strip temperatures are smoother with less fluctuations.The temperature map of the loop segment B2 (Figure 4) shows that at the beginning of the time duration,the first strips of the loop are hotter,and the last ones are colder,but at the end times this pattern is reversed in this loop segment.In loop segment C1 (Figure 4),the temperature fluctuations are mainly observed to start after the end of the flare (22:24),and at the end time (23:00) the temperature is much higher than the beginning.The temperature is also increasing after the flare time(22:24)for Loop C2(see Figure 4).This happens with some oscillations in the strips’ temperatures.So as Figure 4 demonstrates,the temperature increases with some fluctuation in most of the flaring loops’ strips after the flare time.According to these temperature maps,the temperature fluctuations in the flaring loops are increasing at the flaring time and around 20 minutes after that.

    We expect the flaring loops to cool down as a result of heat conduction and radiative cooling.Hence this relative temperature increase should be scrutinized.As we investigated,this temperature rise is also followed by an intensity time-series.As the intensity time-series show,the related intensity in Loop A of the flaring AR increases at the end of the time duration.To be assured,the authors also checked the wavelength of Fe XVIII which has a peak formation temperature of 7×106K (Ugarte-Urra &Warren 2014).By applying the method developed by Warren et al.(2012) the contribution of the Fe XVIII emission line can be isolated from AIA 94 ? to analyze the evolution of hot plasma in the loops.We do it to omit contamination from the cooler plasma (mostly around 1 MK) which also contributes to this AIA channel (Boerner et al.2012).This is done by subtracting the contaminating warm (i.e.,around 1 MK) component to the bandpass.This warm contribution is calculated from a weighted combination of the emission from the AIA 171 ? and 193 ? channels dominated by Fe X and Fe XII emission,respectively.This intensity analysis is done directly and it has not gone through any other process like the thermal analysis.For this purpose,we applied the formulation(1) used by Li et al.(2015).Plots in Figure 5 display the intensity map and the mean intensity variation of the wavelength Fe XVIII,for Loop A of the flaring region,in the top and bottom panels respectively.As these plots show,this intensity is also higher at the end of the time duration with respect to the flare time.It seems to us that the expected cooling has not occurred in these flaring loops yet,even after the flare occurrence in the probed duration due to some plausible reasons.We consider that the mentioned simultaneous CME(see Section 2) which this flare is associated with could cause this increase in temperature.We can be sure that the source of this CME is AR 11283 (Romano et al.2015).This CME is in our flare region,hence the loops receive energy even after the flare occurrence and it is probably the reason why the expected cooling does not occur.

    Figure 5.Normalized intensity map of flaring Loop A for the wavelength Fe XVIII,and mean intensity of Fe XVIII (from top to bottom).The vertical axis is the distance along the loop in Mm for the first plot,and normalized intensity is in the second plot.The horizontal axis is time.The colorbar on the left shows the colors considered for the intensity range.

    The thermal oscillation periods,obtained with the Lomb–Scargle method,do not have the same significance in all strips of the loops,but for most strips of the flaring loops,the significances are very near to one.To be assured about these oscillations,we examined the intensity time-series for each strip of the loops and we observed that this loop’s intensities also manifest intensity oscillations (i.e.,alongside the loop).The most probable dominant periods observed in intensity for wavelength of 171 ? is 18.22,and 16.7 minutes for strips of F-Loop A,16.7,and 18.22 minutes for strips of F-Loop B1,16.70,and 12.52 for F-Loop B2,and 16.7 for F-Loop C1 and F-Loop C2.These periods are on the same order as the observed thermal oscillation periods.The intensity in this time series has not passed any thermal process but still shows oscillation periods close to thermal ones,so we think these results confirm the observation of thermal oscillations.

    4.2.Temperature Analysis of Non-flaring Active Region Loops

    The temperature time-series for different strips of the selected loops of the nonf AR 12194 are calculated using the Lomb–Scargle method.In the following plots (Figure 6),the vertical axis represents the logarithm of the temperature and the horizontal axis corresponds to the time duration.Figure 6 displays the time-series of temperature variations for the first two strips of the nonf Loops A and B.These figures are all coscaled in the range of 5.7–6.9 for the logarithm of temperature(like the flaring loops range).The most powerful periods,observed in most of these nonf loops’strips(listed in Table 2),are from 8.5 to 30 minutes.Comparing the periods of the loops in the flaring region(Table 1)with the nonf ones(Table 2),we see that the temperature periods of the flaring loops have lower values on average and have more diversity than the nonf ones.As Tables 1 and 2 affirm,the mean temperatures of nonf-loops are lower in comparison with the f-loops,a fact we also expected from common sense.The parameter(Max(logT)–Min(logT)) in nonf-loops’ strips is less than that for the flaring loops’ strips.

    Nonf Loop A,divided into 11 strips,has the length of 19.91(Mm) which is the length of the selected part of the loop marked in Figure 2(b).The mean of (Max(logT)–Min(logT))for the strips of nonf Loop A is 0.81.The mean of temperature(log)of this loop segment over time is 5.93±0.10.Nonf Loop B,divided into six strips,has the length 11.11 (Mm),and the mean temperature (log) and the mean of (Max(logT)–Min(logT)) for this loop are 5.99±0.13 and 0.62 respectively.Nonf Loop C,which has five strips,with a length of 10.13(Mm),has the mean temperature (log) of 5.82±0.12 and the mean (Max(logT)–Min(logT)) of 0.56.

    The highest period observed for the temperature oscillations of these nonf loops strips is reported in Table 2.As we observe,the temperature periods in these nonf loops are mostly longer than those of the flaring loops (compare the values listed in Tables 1 and 2).Therefore the temperature oscillations of these loops are a little slower than the flaring ones.

    Figure 7 features the temperature maps of the nonf Loops A,B and C,from top to bottom respectively,as a time series.In each plot,the vertical axis is the distance along the loop in Mm and the horizontal axis is time.The color bar in the left shows the colors considered for the temperature range.Each separated colored part in the map is one strip.These color maps are all plotted with the same color range for the loops of the flaring region.

    As Figure 7 demonstrates,the temperatures of these strips associated with nonf loops are lower and smoother in comparison with the flaring ones (Figure 4).Furthermore,much increase in the temperatures of the strips,which was obvious in the loops of the flaring region toward the end times,is not observed here.The temperatures are also totally lower in the nonf loops in comparison with the flaring loops.Conversely,it seems that different strips of the nonf loops have relatively more similar temperature fluctuations.

    As Figure 8 shows,the peaks of the observed temperature periods for the loops’ strips of the flaring AR (blue ones) and nonf AR (red ones),are around 18 minutes and 30 minutes,respectively.The temperature periods’diversity is higher in the loops’strips of the flaring AR,and shorter temperature periods(less than 10 minutes,nearer to the transverse oscillations periods) are observed in the case of the flaring loops’ strips in comparison with the nonf ones.In addition,Figure 9 shows that the increase and decrease of temperature range,or the difference between maximum and minimum of the temperature value (max(log(T))–min(log(T))),are much higher on average for the loops’ strips of the flaring AR in comparison with the loops’ strips of the nonf one.

    Figure 7.From top to bottom:Temperature map of the nonf Loops A,B and C as a time-series.The vertical axis is distance along the loop in Mm,and the horizontal axis is time.The color-bar in the left displays the colors considered for the temperature range.

    Figure 8.Histogram of the temperature period percentages for the loops’ strips of the flaring (blue bars) and nonf (red bars) ARs.The horizontal axis displays the temperature periods in minutes.

    Figure 9.Histogram of the parameter of (max(log(T))–min(log(T))) for each strip of the loops of the flaring (blue bars) and nonf (red bars) ARs.

    5.Summary

    We reported the temperature oscillations of coronal loops of a flaring AR.We selected the flaring AR 11283 to investigate the thermal structure and treatment of its loops.This region includes a high energy flare x2.1 and the transverse oscillations of two loops associated with it were analyzed before by Jain et al.(2015).They analyzed intensity variations in the wavelength 171 ? in two coronal loops of this region and detected obvious transverse oscillation with periods of roughly 2 minutes and decay times of 5 minutes for these loops(Loops A and B in Figure 1(b)) at the flare time.We were curious to know if the temperature variations follow the transverse oscillations of the loops,or if there is any relation or correlation between them.We also wanted to investigate the thermal fluctuations at the flare time.As a blind test to see the specific thermal properties of the flaring loops,we selected an LOS nonf AR (12194),extracted three segments of its loops and analyzed their thermal treatment.Then we compared the temperature treatment of the loops at the flaring region with the loops of the nonf region to see the differences.We were eager to observe the probable discrepancies between flaring and nonf loops in this respect.

    Here we used data of three loops of the flaring AR (AR 11283)around the time of the Flare X2.1,from 22:10 UT until 23:00 UT on 2011 September 6,plus three loops of the nonf AR(AR 12194),from 08:00:00 UT until 09:00:00 UT of 2014 October 26 (marked in Figures 1 and 2).To calculate the time series of the loop temperature values,we first extracted the loop pixels in each image and then displayed the loop straightly for all the images in the time series of different wavelengths.To do thermal analysis,we used the spatially synthesized Gaussian DEM forward-fitting method proposed by Aschwanden et al.(2015).We calculated the peak temperatures for each strip of the loops.Then we applied the Lomb–Scargle method to analyze temperature oscillations of the time-series for each strip of the loops.

    We observed temperature oscillations which are following the transverse loop oscillations observed by Jain et al.(2015)for the flaring loops.Furthermore,the temperature oscillations in these flaring loops happen before the transverse oscillations start and continue even in the time duration after the transverse oscillations decay.As observed,the temperature oscillations do not decay as rapidly as the transverse oscillations do.Conversely,the strips’ temperatures increase at the end of the oscillating mode and a rather sensible rise is observed in the final temperatures of the f-loops’ segments.The ranges of the obtained periods are from 7 to 28.4 minutes for the flaring loops,and from 8.5 to 30 minutes for the nonf loops.With the onset of an X-flare in F-Loop A,which has a distinct transverse oscillation in the flaring time with period of roughly 2 minutes and decay time of 5 minutes,a temperature oscillation is observed with periods of roughly 10–28.5 minutes in different segments of this loop,and as the transverse oscillation decays in this interval,no special definite decay is observed in its temperature oscillations.

    The temperature periods of the flaring loops are rather shorter than the temperature periods of the nonf loops.The loops of the flaring region display some short temperature oscillation periods in which some are less than 10 minutes(Table 1).Such short periods are more frequently observed for the loops of the flaring AR and,in the case of the nonf ones,are very scarce.We observed that the periods of the flaring loops have more diversity than those of the nonf ones.Based on our confined observations,the nonf loops’ periods are longer and their temperature values are all lower.So,our research showed that thermal structures of the flaring loops differ from the nonf ones in the ways described above.As temperature maps demonstrate,the temperature fluctuations are increasing at the flaring time and around 20 minutes after,in the flaring loops.This happens with some oscillations in strips’ temperature.Conversely,it seems that different strips of the nonf loops have relatively more similar temperature fluctuations.The temperatures are both higher on average in the flaring loops’segments as expected.The significances of the periods,obtained by the Lomb–Scargle method,are calculated for each strip of each loop and the results affirm that these significances for the loops’strips of the flaring region are high and close to one,while for the loops’strips of the nonf region are less than 0.5.Hence the detected periods in the flaring loops’ strips have high significances (near to one) and are oscillations,whereas the detected periods in the nonf loops’strips have less significances in comparison with the flaring ones,and maybe they are just fluctuations.

    Using this method for the coronal loops showed that the oscillation modes obtained for the temperatures of the flaring loops are very close to those of the spatial slow-mode oscillations of the coronal loops,so the origin of temperature oscillation is probably slow-mode waves.Such kind of oscillations often occur in hot coronal loops (l og(T) >6) of ARs,especially the ones associated with small(or micro)flares(Wang et al.2021).The loops of our flaring AR are also hot loops with the mean temperature above this range.They also display intensity oscillations.Hence we think the above evidence confirms the slow-mode oscillations for flaring loops.The temperature of the nonf loops is lower (log(T) <6) and,as discussed above,we believe that the observed oscillationlike periods in nonf loops should be more probably related to the high amplitude fluctuations.

    Comparing the loops of the flaring and nonf regions,we observed that the amplitudes of the fluctuations show a discrepancy.The averages of the parameter (Max(logT)–Min(logT)) in F-Loop A,F-Loop B1,F-Loop B2,F-Loop C1 and F-Loop C2 are 1.21,1.10,0.81,1.48 and 0.88,respectively,and for the nonf region,means of (Max(logT)–Min(logT)) are 0.81,0.62 and 0.56,for nonfloops A,B and C respectively.Therefore the values of the quantity averages of (Max(logT)–Min(logT))for these nonf loops manifest a difference from the flaring ones and are lower.

    Loops of the nonf AR 12194 have a relatively uniform temperature at the beginning of the time interval,which rises slightly at its end.As the Solar Monitor reports in the neighborhood of this region,the flaring AR 12192 happens between multiple flares,and there is a c4.6 class flare occurring at 9:44 UT.Therefore,a possible explanation could be that the abovementioned slight temperature rise in the loops of AR 12194 (in the time interval 8:00–9:00) originated from the influence of an increase in the energy during the pre-flare conditions which exist in AR 12192.

    Hence as our study shows,the temperature of coronal loops of flaring AR changes in an oscillatory manner.Compared with these nonf loops,the flaring loops show higher temperatures on average and higher oscillation periods with higher peaks and deeper valleys.More accurate commentary in this respect requires more extensive statistical research and broader observations.

    Acknowledgments

    The author Narges Fathalian wishes to also express her thanks for the technical support and comments which were received from Dr.Farhad Daii and Dr.Mohsen Javaherian regarding this work.

    成人毛片a级毛片在线播放| 日日摸夜夜添夜夜爱| √禁漫天堂资源中文www| av又黄又爽大尺度在线免费看| 少妇人妻 视频| 精品国产一区二区久久| 亚洲综合色惰| 免费不卡的大黄色大毛片视频在线观看| 国产97色在线日韩免费| 男女国产视频网站| 久久精品夜色国产| 日韩在线高清观看一区二区三区| 男女下面插进去视频免费观看| 国产97色在线日韩免费| 国产精品熟女久久久久浪| 男女午夜视频在线观看| 午夜福利在线免费观看网站| 国产精品秋霞免费鲁丝片| 久久97久久精品| 亚洲精品成人av观看孕妇| 精品人妻在线不人妻| 国产成人精品久久久久久| 菩萨蛮人人尽说江南好唐韦庄| 亚洲,一卡二卡三卡| 日韩一区二区三区影片| 国产伦理片在线播放av一区| 制服诱惑二区| 日本爱情动作片www.在线观看| 国产精品99久久99久久久不卡 | 黄色 视频免费看| 国产亚洲精品第一综合不卡| 涩涩av久久男人的天堂| 国产色婷婷99| 午夜福利乱码中文字幕| 亚洲视频免费观看视频| 国产伦理片在线播放av一区| 欧美 亚洲 国产 日韩一| 中文字幕精品免费在线观看视频| 999精品在线视频| 精品亚洲成国产av| 国产精品不卡视频一区二区| 国产免费现黄频在线看| 免费女性裸体啪啪无遮挡网站| 91午夜精品亚洲一区二区三区| 亚洲精品自拍成人| 男人操女人黄网站| 国产精品一国产av| 蜜桃国产av成人99| 如何舔出高潮| av视频免费观看在线观看| 韩国精品一区二区三区| 十八禁高潮呻吟视频| 亚洲精品久久成人aⅴ小说| 欧美老熟妇乱子伦牲交| 久久久久久久久久人人人人人人| 一级毛片黄色毛片免费观看视频| 免费观看av网站的网址| 国产黄色免费在线视频| 国产成人精品福利久久| 欧美激情极品国产一区二区三区| 伦理电影大哥的女人| 看十八女毛片水多多多| 大陆偷拍与自拍| 欧美日韩一区二区视频在线观看视频在线| 1024视频免费在线观看| 国产一区二区三区av在线| 中文字幕色久视频| 国产成人aa在线观看| 午夜福利乱码中文字幕| 成人黄色视频免费在线看| 亚洲情色 制服丝袜| 韩国高清视频一区二区三区| 欧美激情极品国产一区二区三区| 七月丁香在线播放| 又黄又粗又硬又大视频| 国产精品香港三级国产av潘金莲 | 久热久热在线精品观看| 久久久国产精品麻豆| 视频区图区小说| 一二三四中文在线观看免费高清| videos熟女内射| 精品国产乱码久久久久久男人| 女人精品久久久久毛片| 日本av手机在线免费观看| 亚洲国产最新在线播放| 91在线精品国自产拍蜜月| videossex国产| 极品少妇高潮喷水抽搐| 制服丝袜香蕉在线| 久久ye,这里只有精品| 日韩 亚洲 欧美在线| 精品国产一区二区三区四区第35| 成人毛片a级毛片在线播放| 亚洲欧洲国产日韩| 男的添女的下面高潮视频| av网站在线播放免费| 亚洲国产精品成人久久小说| 久久久精品免费免费高清| 精品福利永久在线观看| 菩萨蛮人人尽说江南好唐韦庄| 国产成人精品婷婷| 免费播放大片免费观看视频在线观看| 美女脱内裤让男人舔精品视频| 亚洲精品美女久久久久99蜜臀 | 看免费av毛片| 性色avwww在线观看| 亚洲三区欧美一区| 少妇的逼水好多| 久久av网站| 国产成人精品在线电影| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日韩精品免费视频一区二区三区| 亚洲精品一二三| 性高湖久久久久久久久免费观看| 国产熟女欧美一区二区| 久久久久久久久久久久大奶| 亚洲熟女精品中文字幕| 久久久久精品人妻al黑| 综合色丁香网| 亚洲国产欧美网| 亚洲五月色婷婷综合| 18禁观看日本| 人妻一区二区av| 亚洲欧美一区二区三区黑人 | 久久精品久久精品一区二区三区| 2021少妇久久久久久久久久久| 少妇熟女欧美另类| 亚洲第一av免费看| 免费大片黄手机在线观看| 久久精品人人爽人人爽视色| 成人黄色视频免费在线看| 人人澡人人妻人| 春色校园在线视频观看| 精品少妇一区二区三区视频日本电影 | 久久精品国产亚洲av天美| 97在线人人人人妻| 少妇被粗大的猛进出69影院| 黑丝袜美女国产一区| 新久久久久国产一级毛片| 国产精品久久久久久久久免| 国产亚洲av片在线观看秒播厂| 久久久欧美国产精品| 叶爱在线成人免费视频播放| 母亲3免费完整高清在线观看 | 日日摸夜夜添夜夜爱| 国产片特级美女逼逼视频| 自线自在国产av| 宅男免费午夜| 国产一区二区激情短视频 | 性少妇av在线| 日本爱情动作片www.在线观看| 亚洲国产精品成人久久小说| 丁香六月天网| 捣出白浆h1v1| 最近中文字幕2019免费版| 国产极品天堂在线| 免费播放大片免费观看视频在线观看| av网站在线播放免费| 91精品国产国语对白视频| 老司机亚洲免费影院| 亚洲伊人色综图| 久久久久精品人妻al黑| 2022亚洲国产成人精品| 久久久久视频综合| 伦理电影免费视频| videossex国产| 如何舔出高潮| 黄片小视频在线播放| 久久人妻熟女aⅴ| 亚洲av成人精品一二三区| 男女边吃奶边做爰视频| 亚洲精品第二区| 午夜福利一区二区在线看| 亚洲一码二码三码区别大吗| 各种免费的搞黄视频| 又黄又粗又硬又大视频| 精品亚洲乱码少妇综合久久| 成人免费观看视频高清| 成人毛片a级毛片在线播放| 国产成人91sexporn| 午夜福利影视在线免费观看| av在线app专区| 亚洲av电影在线进入| 中文字幕人妻丝袜制服| 国产国语露脸激情在线看| 青春草亚洲视频在线观看| 国产在线一区二区三区精| 男女啪啪激烈高潮av片| 免费观看性生交大片5| 精品亚洲成国产av| 中文字幕人妻熟女乱码| 九色亚洲精品在线播放| 99精国产麻豆久久婷婷| 9色porny在线观看| 国产淫语在线视频| 五月伊人婷婷丁香| 日韩 亚洲 欧美在线| 五月开心婷婷网| 国产激情久久老熟女| 老汉色∧v一级毛片| 人体艺术视频欧美日本| 免费观看性生交大片5| 久久国产精品男人的天堂亚洲| 只有这里有精品99| 久久久久精品性色| 熟女电影av网| 久久97久久精品| 日本欧美视频一区| 成人国产麻豆网| 水蜜桃什么品种好| av.在线天堂| 亚洲av电影在线观看一区二区三区| 乱人伦中国视频| 日韩av不卡免费在线播放| 一区二区三区四区激情视频| 午夜福利一区二区在线看| 搡女人真爽免费视频火全软件| 人妻少妇偷人精品九色| 国产在线免费精品| 久久ye,这里只有精品| 极品少妇高潮喷水抽搐| 欧美日韩一级在线毛片| 丰满饥渴人妻一区二区三| 叶爱在线成人免费视频播放| 色视频在线一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 91国产中文字幕| 我要看黄色一级片免费的| 高清视频免费观看一区二区| 2018国产大陆天天弄谢| 精品亚洲成a人片在线观看| 日产精品乱码卡一卡2卡三| 99九九在线精品视频| 亚洲,欧美精品.| 一级毛片 在线播放| 热re99久久精品国产66热6| 香蕉丝袜av| 9色porny在线观看| 少妇被粗大的猛进出69影院| 自拍欧美九色日韩亚洲蝌蚪91| 成人二区视频| 18在线观看网站| 高清欧美精品videossex| 亚洲久久久国产精品| 午夜激情久久久久久久| 国产亚洲av片在线观看秒播厂| 国产一区二区三区av在线| 亚洲一区中文字幕在线| 最近中文字幕2019免费版| 久久久久久久久免费视频了| 午夜免费观看性视频| 国产一区二区在线观看av| 国产亚洲一区二区精品| 久久久久久久久久久久大奶| 老司机亚洲免费影院| 国产白丝娇喘喷水9色精品| 搡老乐熟女国产| 天天影视国产精品| 婷婷色综合大香蕉| 少妇 在线观看| 香蕉丝袜av| 啦啦啦在线免费观看视频4| 国产午夜精品一二区理论片| 国产激情久久老熟女| 亚洲少妇的诱惑av| 精品国产乱码久久久久久小说| a 毛片基地| 久久国产精品大桥未久av| 午夜免费鲁丝| 久热这里只有精品99| 精品久久久久久电影网| 成年av动漫网址| 青草久久国产| 美女国产高潮福利片在线看| 久久久久久久久久人人人人人人| av在线播放精品| 成人午夜精彩视频在线观看| 国产97色在线日韩免费| 亚洲少妇的诱惑av| 亚洲国产av影院在线观看| 另类精品久久| 男女啪啪激烈高潮av片| 人妻少妇偷人精品九色| 欧美激情极品国产一区二区三区| 久久久久视频综合| 老女人水多毛片| 99热国产这里只有精品6| 一区二区三区激情视频| 日韩av在线免费看完整版不卡| 美国免费a级毛片| 99re6热这里在线精品视频| 夫妻午夜视频| 青春草视频在线免费观看| 日日撸夜夜添| 乱人伦中国视频| 成人黄色视频免费在线看| 少妇 在线观看| 高清黄色对白视频在线免费看| 久久精品国产综合久久久| 宅男免费午夜| 中国国产av一级| 成年动漫av网址| 亚洲精品乱久久久久久| 丝袜喷水一区| 可以免费在线观看a视频的电影网站 | 亚洲av综合色区一区| 欧美精品人与动牲交sv欧美| av女优亚洲男人天堂| 中文字幕人妻丝袜制服| 春色校园在线视频观看| 夫妻性生交免费视频一级片| av卡一久久| 在线天堂中文资源库| 久久ye,这里只有精品| 成年动漫av网址| 可以免费在线观看a视频的电影网站 | 人妻 亚洲 视频| 国产精品亚洲av一区麻豆 | 欧美最新免费一区二区三区| 免费观看无遮挡的男女| 在线天堂中文资源库| 国产精品女同一区二区软件| 欧美+日韩+精品| 欧美少妇被猛烈插入视频| 亚洲三级黄色毛片| 日本猛色少妇xxxxx猛交久久| 91成人精品电影| 国产一区二区在线观看av| 三上悠亚av全集在线观看| 国产精品一国产av| 国产精品久久久久久精品古装| 丝袜人妻中文字幕| 侵犯人妻中文字幕一二三四区| 男的添女的下面高潮视频| 热99国产精品久久久久久7| 人体艺术视频欧美日本| 中文字幕精品免费在线观看视频| 欧美xxⅹ黑人| 热99国产精品久久久久久7| 精品国产乱码久久久久久男人| 国产免费一区二区三区四区乱码| 五月开心婷婷网| 啦啦啦中文免费视频观看日本| 深夜精品福利| 天堂8中文在线网| 亚洲视频免费观看视频| 精品第一国产精品| 欧美+日韩+精品| 精品第一国产精品| 亚洲精品美女久久久久99蜜臀 | 少妇的丰满在线观看| 久久久久久伊人网av| 街头女战士在线观看网站| 欧美 亚洲 国产 日韩一| 国产亚洲精品久久久久5区| 国产熟女xx| 国产成人精品在线电影| 女人爽到高潮嗷嗷叫在线视频| 国产精品自产拍在线观看55亚洲| 久久精品亚洲精品国产色婷小说| 免费观看人在逋| 韩国av一区二区三区四区| 嫩草影院精品99| 在线观看舔阴道视频| 大陆偷拍与自拍| 丝袜美足系列| 99精国产麻豆久久婷婷| 国产亚洲精品一区二区www| 久久中文字幕人妻熟女| 99国产精品一区二区三区| 18禁美女被吸乳视频| 亚洲va日本ⅴa欧美va伊人久久| 久久人妻熟女aⅴ| 午夜日韩欧美国产| 日韩欧美免费精品| 美女国产高潮福利片在线看| 亚洲一区二区三区欧美精品| 妹子高潮喷水视频| 一区在线观看完整版| 国产一区二区三区视频了| 国产精华一区二区三区| 国产精品综合久久久久久久免费 | 欧美丝袜亚洲另类 | 精品久久久久久久毛片微露脸| 色播在线永久视频| 丝袜美腿诱惑在线| 欧美日韩瑟瑟在线播放| 性色av乱码一区二区三区2| 亚洲一区二区三区不卡视频| 大型av网站在线播放| 欧美日韩精品网址| 国产精华一区二区三区| 欧美一区二区精品小视频在线| 欧美色视频一区免费| 亚洲五月色婷婷综合| 免费人成视频x8x8入口观看| 国产视频一区二区在线看| 国产精华一区二区三区| 日韩免费av在线播放| 国产成人av教育| av天堂在线播放| e午夜精品久久久久久久| 男女午夜视频在线观看| 成人三级做爰电影| 久久精品国产清高在天天线| 亚洲精品久久午夜乱码| 国产免费现黄频在线看| 久久精品亚洲熟妇少妇任你| 日韩人妻精品一区2区三区| 黑人巨大精品欧美一区二区mp4| 高清毛片免费观看视频网站 | 国产精品 国内视频| 黑人巨大精品欧美一区二区mp4| 国产精华一区二区三区| 大型黄色视频在线免费观看| 欧美丝袜亚洲另类 | 亚洲一区二区三区欧美精品| 亚洲一区二区三区色噜噜 | 丰满饥渴人妻一区二区三| 国产欧美日韩一区二区精品| 成人国产一区最新在线观看| 夜夜夜夜夜久久久久| 色在线成人网| av网站在线播放免费| 黄色女人牲交| 人妻丰满熟妇av一区二区三区| 女性生殖器流出的白浆| 国产精品久久久久成人av| 国产成人精品久久二区二区91| 老熟妇乱子伦视频在线观看| 午夜免费激情av| 国产高清视频在线播放一区| 首页视频小说图片口味搜索| 黄色女人牲交| 窝窝影院91人妻| av在线天堂中文字幕 | 大陆偷拍与自拍| 国产精品野战在线观看 | 麻豆成人av在线观看| 欧美成狂野欧美在线观看| 大型黄色视频在线免费观看| 老司机午夜福利在线观看视频| 两个人免费观看高清视频| 两个人看的免费小视频| 国产精品爽爽va在线观看网站 | 午夜免费观看网址| 法律面前人人平等表现在哪些方面| 在线永久观看黄色视频| 亚洲色图综合在线观看| 99精国产麻豆久久婷婷| 97人妻天天添夜夜摸| 大香蕉久久成人网| 韩国精品一区二区三区| 欧美另类亚洲清纯唯美| 亚洲男人天堂网一区| 成人精品一区二区免费| 国产成人精品久久二区二区91| 一边摸一边做爽爽视频免费| 视频区图区小说| 美女大奶头视频| 级片在线观看| 亚洲精品国产精品久久久不卡| 一本综合久久免费| 亚洲一区二区三区欧美精品| 亚洲专区中文字幕在线| 首页视频小说图片口味搜索| 久久天堂一区二区三区四区| 啦啦啦在线免费观看视频4| 欧美激情高清一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 一区在线观看完整版| 最好的美女福利视频网| 欧美中文综合在线视频| 亚洲熟妇熟女久久| 热99国产精品久久久久久7| 日韩大码丰满熟妇| 乱人伦中国视频| 免费久久久久久久精品成人欧美视频| 99热只有精品国产| 窝窝影院91人妻| 欧美激情高清一区二区三区| 国产精品亚洲一级av第二区| 久久久精品国产亚洲av高清涩受| 亚洲片人在线观看| 可以在线观看毛片的网站| 黄频高清免费视频| 亚洲欧美激情在线| 日本撒尿小便嘘嘘汇集6| 欧美黑人精品巨大| 日日摸夜夜添夜夜添小说| 国产伦人伦偷精品视频| 在线观看www视频免费| 日本vs欧美在线观看视频| 在线观看免费视频网站a站| 人妻丰满熟妇av一区二区三区| 俄罗斯特黄特色一大片| 一边摸一边做爽爽视频免费| 亚洲人成伊人成综合网2020| 国产免费av片在线观看野外av| 啦啦啦在线免费观看视频4| 精品免费久久久久久久清纯| 国产99白浆流出| 十分钟在线观看高清视频www| 国产精品国产av在线观看| 一级毛片精品| 亚洲黑人精品在线| 久久久久久人人人人人| 校园春色视频在线观看| 国产黄a三级三级三级人| 中文字幕精品免费在线观看视频| 亚洲五月婷婷丁香| 色在线成人网| 久久人人爽av亚洲精品天堂| 极品教师在线免费播放| 亚洲熟女毛片儿| 日本wwww免费看| 一区二区日韩欧美中文字幕| 丁香六月欧美| 亚洲中文av在线| 日日摸夜夜添夜夜添小说| 国产精品九九99| 日韩大尺度精品在线看网址 | 日韩三级视频一区二区三区| 久久人妻av系列| 精品国产美女av久久久久小说| 国产高清国产精品国产三级| 亚洲成人国产一区在线观看| 国产成人av激情在线播放| 女生性感内裤真人,穿戴方法视频| 美女福利国产在线| 51午夜福利影视在线观看| 亚洲欧美日韩另类电影网站| 国产欧美日韩一区二区精品| 午夜免费成人在线视频| 熟女少妇亚洲综合色aaa.| 亚洲精品一二三| 国产精品野战在线观看 | 欧美不卡视频在线免费观看 | 国产成人精品久久二区二区免费| 久久中文字幕人妻熟女| 亚洲色图av天堂| 欧美精品一区二区免费开放| 久久久久久人人人人人| 精品无人区乱码1区二区| 免费看a级黄色片| 99久久综合精品五月天人人| 国产不卡一卡二| 最新美女视频免费是黄的| 在线天堂中文资源库| 日日爽夜夜爽网站| 亚洲片人在线观看| 中国美女看黄片| 国产一区二区三区视频了| 精品无人区乱码1区二区| 国产一区在线观看成人免费| 免费久久久久久久精品成人欧美视频| 亚洲男人天堂网一区| 精品日产1卡2卡| 亚洲成人免费av在线播放| 熟女少妇亚洲综合色aaa.| 午夜影院日韩av| 日本三级黄在线观看| 亚洲精品一二三| 精品一区二区三区av网在线观看| 日韩欧美三级三区| 亚洲国产中文字幕在线视频| 黄色女人牲交| 久久国产精品人妻蜜桃| 国产精品久久久人人做人人爽| 亚洲欧美日韩高清在线视频| 午夜视频精品福利| 久久久国产成人精品二区 | 久9热在线精品视频| 黑人巨大精品欧美一区二区蜜桃| 国产97色在线日韩免费| 国产又爽黄色视频| 国产成人啪精品午夜网站| 亚洲精品一二三| 久久 成人 亚洲| 亚洲五月天丁香| 亚洲av成人不卡在线观看播放网| 亚洲精品国产色婷婷电影| 久久中文字幕一级| 999精品在线视频| 人妻久久中文字幕网| 久久久久久久午夜电影 | bbb黄色大片| 亚洲视频免费观看视频| netflix在线观看网站| 午夜精品在线福利| 一区福利在线观看| 91在线观看av| 男人操女人黄网站| 两个人免费观看高清视频| 美国免费a级毛片| 别揉我奶头~嗯~啊~动态视频| 午夜福利一区二区在线看| 日韩免费高清中文字幕av| 久9热在线精品视频| 国产免费av片在线观看野外av| 香蕉久久夜色| 91国产中文字幕| 一a级毛片在线观看| 精品国产国语对白av| 波多野结衣一区麻豆| av欧美777| 亚洲av成人一区二区三| 欧美日韩视频精品一区| 在线观看免费视频网站a站| 女人爽到高潮嗷嗷叫在线视频| 久9热在线精品视频| 另类亚洲欧美激情| 亚洲自偷自拍图片 自拍| 国产欧美日韩综合在线一区二区| 日韩欧美一区二区三区在线观看|