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

    Impact of Updraft on Neutralized Charge Rate by Lightning in Thunderstorms:A Simulation Case Study

    2015-12-12 12:57:03WANGFei王飛ZHANGYijun張義軍andZHENGDong鄭棟
    Journal of Meteorological Research 2015年6期
    關鍵詞:義軍王飛

    WANG Fei(王飛),ZHANG Yijun(張義軍),and ZHENG Dong(鄭棟)

    1 State Key Laboratory of Severe Weather,Chinese Academy of Meteorological Sciences (CAMS),Beijing 100081

    2 Laboratory of Lightning Physics and Protection Engineering,CAMS,Beijing 100081

    Impact of Updraft on Neutralized Charge Rate by Lightning in Thunderstorms:A Simulation Case Study

    WANG Fei1,2*(王飛),ZHANG Yijun1,2(張義軍),and ZHENG Dong1,2(鄭棟)

    1 State Key Laboratory of Severe Weather,Chinese Academy of Meteorological Sciences (CAMS),Beijing 100081

    2 Laboratory of Lightning Physics and Protection Engineering,CAMS,Beijing 100081

    The rate of neutralized charge by lightning(RNCL)is an important parameter indicating the intensity of lightning activity.The total charging rate(CR),the CR of one kind of polarity(e.g.,negative)charge (CROP),and the outflow rate of charge on precipitation(ORCP)are proposed as key factors impacting RNCL,based on the principle of conservation of one kind of polarity charge in a thunderstorm.In this paper,the impacts of updraft on CR and CROP are analyzed by using a 3D cloud resolution model for a strong storm that occurred in Beijing on 6 september 2008.The results show that updraft both promotes and inhibits RNCL at the same time.(1)Updraft always has a positive influence on CR.The correlation coefficient between the updraft volume and CR can reach 0.96.Strengthening of the updraft facilitates strengthening of RNCL through this positive influence.(2)Strengthening of the updraft also promotes reinforcement of CROP.The correlation coefficient between the updraft volume and CROP is high(about 0.9),but this promotion restrains the strengthening of RNCL because the strengthening of CROP will,most of the time,inhibit the increasing of RNCL.(3)Additionally,increasing of ORCP depresses the strengthening of RNCL.In terms of magnitude,the peak of ORCP is equal to the peak of CR.Because precipitation mainly appears after the lightning activity finishes,the depression effect of ORCP on RNCL can be ignored during the active lightning period.

    lightning activity,updraft,neutralized charge,charging rate,charge amount

    1.Introduction

    Lightning activity is an important phenomenon of strong convective systems,and has a close relationship with observed atmospheric convection.For example,many researchers have found that the radar echo top of a thunderstorm is related to the lightning rate(Shackford,1960;Jacobson and Krider,1976; Williams,1985).Based on this correlation,Price and Rind(1992)established a parameterized scheme between the cloud top and lightning rate in a global numerical model.Also,the lifting condensation level (LCL)is found to have a very high correlation with the log value of the lightning rate(Williams et al.,2005; Fuchs et al.,2014),and the combination of dry bulb temperature and convective available potential energy (CAPE)has an evident correlation with the lightning rate,although the correlation coefficient is not very high(Fuchs et al.,2014).Williams and Stanfill(2002) found that high frequency lightning activity tended to be observed in thunderstorms with high LCL and moderate CAPE,rather than in thunderstorms with high CAPE and low LCL.

    The above studies indicate that many parameters characterizing strong convection bear a relationship with lightning activity.In fact,most of these param-

    eters directly reflect the characteristics of updraft.Observations have shown that the existence of a strong updraft in the mixed-phase region(–40–0℃)is necessary to produce lightning(Workman and Reynolds, 1949;Williams and Lhermitte,1983;Dye et al.,1988; Rutledge et al.,1992;Carey and Rutledge,1996;Petersen et al.,1996,1999;Pickering et al.,1998;Baker et al.,1999;Black and Hallett,1999;Wang et al., 2009;Heymsfield et al.,2010;Palucki et al.,2011;Yao et al.,2013;Zhang et al.,2013;Reinhart et al.,2014; Wang et al.,2015).Moreover,quantitative observational studies on the role of convective updrafts have been carried out.Lang and Rutledge(2002)reported that thunderstorms that predominantly produced positive cloud-to-ground lightning were accompanied by a large area of strong updraft(with speed>20 m s?1).Tessendorf et al.(2005)and Wiens et al.(2005) found a good correlation between the updraft volume with speed exceeding 10 m s?1and lightning activity.Deierling and Petersen(2008)found that updraft volumes above the–5℃isotherm,with vertical velocities greater than 5 and 10 m s?1,were well correlated with mean total lightning activity;the correlation coefficients were 0.93 and 0.92,respectively.They also indicated that the relationship between updraft volume and total lightning activity varied as a function of temperature and,depending on the temperature,also with region.Barthe et al.(2010)also found that the maximum updraft speed correlated well with lightning frequency only in thunderstorms with a very strong updraft.

    Although the close relationship between updraft and lightning activity has been confirmed by many studies,it is still not clear as to the mechanism through which updraft affects lightning activity. Much research has focused on the impact of updraft on the charging process,while other studies have also indicated that updraft should have an evident influence on the discharging process(Zhang et al.,1999,2003; Miller et al.,2001).However,most of the work carried out on these two problems amounts to logical speculation based on observations,because of a lack of direct observations of updrafts and the charging process.Carey and Rutledge(1996)found that the volume and mass of graupel appearing in the upper part of a strong updraft correlated closely with the lightning rate,suggesting that updraft may affect the lightning frequency through its impact on the charging process occurring on graupel.Lang and Rutledge (2002)indicated that a strong updraft produced more super-cooled water and caused more collision and intensive charging between ice particles,resulting in the close relationship between updraft and lightning activity.Deierling and Petersen(2008)proposed that strengthening of the updraft produces more ice particles in the mixed-phase region,and thus a higher number of collisions with subsequent charge separation, which relate the updraft to the lightning frequency via positive correlation.They believed that the relationship between the updraft volume and lightning frequency is less regime-variant as the cold-cloud active precipitation-based charging zone is approached. Meanwhile,theoretical analysis has indicated that the updraft area—in particular,the area in the subfreezing layer of cloud—should be related to the generator current density,which in turn is linked to the electrical activity of the storm(Boccippio,2002).

    Lightning frequency is,in many studies,taken to be the main parameter representative of the intensity of lightning activity.Usually,a thunderstorm with deeper convection will produce more frequent lightning activity(McCaul et al.,2009).However,recent observations have also revealed that the lightning frequency in a particular period is relative to the location of the deep convective updraft:lightning flashes tend to be smaller and more frequent in and near a strong updraft,and the flashes will show the opposite tendency when they are far away from the strong updraft(Carey et al.,2005;Kuhlman et al.,2009;Weiss et al.,2012; Bruning and MacGorman,2013).Recent research indicates that the flash length scales and flash areas are generally anti-correlated to the lightning flash rate (Mecikalski et al.,2015),reminding us that the lightning frequency cannot represent the characteristics of lightning activity intensity comprehensively.In this paper,we define RNCL,which is the charge neutralized by lightning flashes on the whole-thunderstorm scale in a specific period of time,to represent the in-

    tensity of lightning activity.This parameter considers the effects of lightning frequency and flash scales comprehensively,and thus should be more appropriate for the lightning activity intensity.

    Additionally,most previous studies support the positive influence of updraft on the charging process, and subsequently,the lightning activity. However, some researchers have also found that the relationship between updraft and lightning activity could sometimes vary,or be evident only under certain conditions (Deierling and Petersen,2008).This suggests that the impact of updraft on lightning activity is complex, and that present understanding is insufficient.Accordingly,the present study first aims to identify the main factors affecting RNCL,and then examine the relationships between updraft and these factors,such that the impact of updraft on RNCL can be analyzed comprehensively.This study is important not only for understanding the relationship between updraft and electrical activity in a thunderstorm,but also for furthering research progress on lightning activity.

    2.Model and definitions

    2.1Model

    The model used here is a 3D cloud resolution model constructed firstly by Hu and He(1988),and then coupled with charging and discharging processes by Tan et al.(2006).The nonhydrostatic balance equation is adapted to construct the dynamic framework.The air is treated as a tack free and compressible fluid.Hydrometeors in clouds include five categories:cloud drop,ice crystal,rain drop,graupel,and hail.The sedimentation and the dragging force of all hydrometeors,except for cloud drop,are considered. Twenty seven kinds of cloud microphysical processes, such as evaporation-condensation,coagulation,automatic transformation,nucleation-multiplication,and melting-freezing,are simulated in the model.The initiation field of the model can be established based on sounding data.The convection is initiated by a damp and hot bubble,which is a restricted perturbation field with higher humidity and temperature than the environment added in the low level artificially,and developed according to the calculation result of buoyancy in the vertical motion equation.

    Three charging schemes,including the noninductive charging process in the process of collision and rebound between graupel/hail and ice crystal(Gardiner et al.,1985;Pereyra et al.,2000;Tan et al.,2006),the inductive charging process between graupel/hail and cloud drop in the dry growth process(Ziegler et al,1991),and the charging process due to the surface temperature difference in the ice crystal multiplication process(Hallet and Saunders,1979), are used here to simulate the charging process in the thunderstorm.

    In the discharging scheme,the electrical field threshold of lightning initiation was defined by Marshall et al.(1995)as a function of air density,which is related to height;the main points of the random discharging model proposed by Mansell et al.(2002)are adopted in the discharging scheme;and the propagation of positive and negative leaders and the inductive charge in the lightning channel are all calculated based on the concept of the bidirectional leader(Kasemir, 1960,2012).The lightning channel is treated as a good conductor with resistance when the electrical potential of a new grid of the lightning channel is calculated by using the formula:

    where φiniis the reference potential,usually replaced by the environmental potential value;p is the polarity of charge in the channel;n is the number of path segments that the channel has propagated;is the inner electrical field in the channel,set to a constant of 500 V m?1;and diis the length of every path segment.After a new step of the lightning channel is finished,the distribution of the electric field outside the channel is recalculated.The inductive charge density within the channel,is calculated as follows:

    where ε is the dielectric constant of air,and n indicates the normal direction of the conductor surface.When there is no probability for development of both ends

    of the lightning channel,or either end of the channel reaches the border of the simulation domain,this simulated discharging process will stop.After a discharging process stops,all charges in the channel will be redistributed according to the proportion of surface area of hydrometeors within and near the channel(the original charge on the hydrometeors will be ignored). The redistribution formula(Ziegler and MacGorman, 1994;Mansell et al.,2002)is as follows:

    The horizontal resolution and vertical resolution of the model are set to 1 and 0.5 km,respectively.The time step is 1 s,and the simulation time is set to 60 min.

    2.2 Definitions

    In the analysis,RNCL is the total charge neutralized by all lightning flashes in the thunderstorm, occurring within the specified time interval.CR is the total charge of one kind of polarity(e.g.,negative) of two polarities obtained by all kinds of hydrometeors in the thunderstorm,in the same time interval, through the non-inductive charging process,the inductive charging process,and the charging process in the ice crystal multiplication process.The total amount of one kind of polarity charge(TAOC)in the thunderstorm is the rest of TAOC in the thunderstorm after all lightning flashes in the time interval have finished, and CROP in the thunderstorm is the change between TAOC in the present time interval and TAOC in the previous interval.ORCP is the total charge of one kind of polarity(e.g.,negative)on the precipitation falling out of the thunderstorm in the time interval,and the precipitation is all hydrometeors that fall down to the height at or below 1 km.The updraft region is defined as the region where the vertical speed is upward and greater than 1 m s?1.

    3.General case description

    A series of strong convective systems occurred over Beijing from the night of 6 September to the morning of 7 September 2008.Based on the sounding at 1200 UTC 6 September in Beijing,a thunderstorm simulation was carried out by using the 3D charging-discharging cloud resolution model.The simulation results are shown in Fig.1.The first simulated lightning occurs at the 16th minute,and ends at the 31st minute. In total,the simulated lightning activity lasts for about 15 min.The simulated lightning frequency increases at first,before the 19th minute,when the frequency reaches its peak of 12 flash min?1,and then decreases(Fig.1a).The observation of the first thunderstorm cell of the system,based on on SAFIR3000,indicates the lightning frequency of

    the thunderstorm at the first 15 minutes to be similar to the simulated lightning frequency.The observed lightning frequency also increases at first,and reaches its peak of 14 flash min?1at the 6th minute after the start of lighting activity.The observed peak value is very close to that value of the simulated lightning frequency.

    Fig.1.(a)Comparison between the observed and simulated lightning flash frequency and(b)charge structure(color scale)along the vertical profile of the positive charge center at the 24th minute,and the distribution(line)of the heights of lightning radiation.Sources from SAFIR3000.

    According to the height distribution of the radiation sources of the lightning observed by SAFIR3000, the peak distribution range is from 9 to 11 km,which indicates the height of the main positive charge layer should also be located near the height of 9–11 km. Along the vertical profile of the simulated thunderstorm,together with the position of the maximum positive charge density when the simulated lighting frequency reaches its peak,the charge structure in the thunderstorm illustrates a dipolar structure,with the main positive charge layer above the main negative charge layer.The height range of the main positive charge layer is also from 9 to 11 km,which is the same as the height range of the main positive charge layer derived from the observation.

    The comparisons above suggest that,in terms of the lighting frequency and charge structure,the simulated thunderstorm is very similar to the observed thunderstorm in the system.Therefore,the simulated results are considered credible and appropriate for further analysis.

    In the simulated results,RNCL keeps increasing after lightning activity starts at the 13th minute. From the 23rd to 24th minutes,RNCL reaches a peak of about 408 C min?1;then,it decreases gradually until lightning activity ends at the 31st minute.

    4.Simulation result analysis

    From the perspective of microphysical processes, lightning is not only a simple neutralized charge process,but also a complex process of interaction with charge structure.In fact,lightning is a non-transient process.During the process,charge that is widely distributed in the thunderstorm(in spatial terms)mixes with charge deposited in the lightning channel through the processes of turbulent exchange,advection,gravity setting and charging.Part of the charge is neutralized,and part of the residual charge will participate in the reconstruction of the charge structure before the next lightning process begins(Tan et al.,2007).However,on larger spatial and longer temporal scales,the complexity of an individual lightning process can be neglected if the balance between the source and remit of the charge in the whole thunderstorm is considered; only the overall result of all lightning in the thunderstorm within a particular time interval can be shown. The source and remit of the charge with one kind of polarity in a thunderstorm should be balanced,as in the following formula:

    4.1CR

    The model used in this study includes three kinds of charging processes:the non-inductive charging process during the collision and rebound between graupel/hail and ice crystal(Gardiner et al.,1985);the inductive process between graupel/hail and cloud drop (Ziegler et al,1991);and the charging process during ice crystal multiplication(Hallett and Saunders, 1979).Although these three charging processes are only part of all charging processes in the thunderstorm,many studies have demonstrated that at least the first two charging processes are the main charging processes generating most of the charge in the thunderstorm.In particular,the non-inductive charging process considered here should be the most important charging process in the thunderstorm(Takahashi, 1978,1984).

    The simulation verifies that the CR of the noninductive charging process during the lightning activity is considerably greater than the CRs of the inductive charging process and the ice crystal multiplication charging process(Fig.2a).The mean CR of the non-inductive charging process is about 6.54×109nC s?1,which is almost three times the mean CR of the inductive charging process.Even from the 17th to 20th minutes,the CR of the non-inductive charging process is one order of magnitude higher than the CR of the inductive charging process.The mean CR of the ice crystal multiplication charging process is the smallest;the magnitude is only 107nC s?1,which is much less than the first two CRs.This indicates that the CR of the non-inductive charging process is dominant,at least during the time of main lightning activity.Therefore,using the CR of the non-inductive charging process to represent the overall CR during lightning activity should be reasonable.

    A comparison of the updraft volume and the CR of the non-inductive charging process is illustrated in Fig.2b.It reveals that they have a similar trend: before the 19th minute,they keep increasing.Then, after the 19th minute,both of them begin to decrease. Statistically,there is good correlation between them; the correlation coefficient is 0.87.This suggests that the updraft volume affects the intensity of the CR of the non-inductive charging process,which process is representative of the total CR in the thunderstorm during lightning activity.The positive correlation between the CR of the non-inductive charging process and the updraft volume reflects the fact that the CR as a whole is also positively correlated with the updraft volume.This result is coincident with many observations,in which researchers have found good correlation between the updraft volume and lightning frequency (Tessendorf et al.,2005;Wiens et al.,2005;Deierling and Petersen,2008;Barthe and Barth,2008;Barthe et al.,2010).

    The spatial relationship between the region with the non-inductive charging process and the updraft (Fig.3)indicates that:(1)The main non-inductive charging process occurs in the top part of the updraft when the updraft is strong(Fig.3a);(2)The main non-inductive charging process still occurs near the top part of the updraft when the updraft begins to weaken and the vertical channel of the updraft becomes gradually inclined(Fig.3b).However,the region with the non-inductive charging process becomes wider and its center starts to deviate to the center of the updraft in the horizontal direction,although an obvious charging process still occurs around the core of the updraft;(3)The main non-inductive charging process still occurs over the region with weak updraft when the updraft weakens considerably(Fig.3c),and the region with the non-inductive charging process and the region with the updraft show evident separation in the horizontal direction;(4)When the updraft continues to weaken and the vertical channel of the updraft almost disappears,the region with the evident noninductive charging process also lessens until it disappears,but it is still distributed near the region with

    a certain updraft speed(Figs.3d–f).

    Fig.2.(a)Efficiencies of three types of charging mechanism and(b)evolution of the updraft volume and non-inductive charging rate.

    Fig.3.Charging rate(contours;absolute charging rate>0.01 nC m?3s?1)and wind field(vectors;only winds with speed greater than 1 m s?1are shown)in the profile of the non-inductive charging rate center at the(a)13th,(b)16th, (c)19th,(d)22nd,(e)25th,and(f)28th minute.

    The region with the non-inductive charging process should bear a close relationship with the updraft region. Indeed,when a strong updraft exists,the main region with the non-inductive charging process is always located in and near the updraft region.Despite the evident updraft disappearing and the region with the non-inductive charging process not coinciding with the updraft region spatially,the non-inductive charging process nevertheless occurs near the region with weak updraft speed.This should be the reason that causes updraft to correlate well with the CR of the non-inductive charging process.This correlation should lead to a positive influence of updraft on the non-inductive charging process.

    4.2CROP

    The charging process on hydrometeors primarily occurs within and near the region of updraft. Charged hydrometeors are carried to a wider region

    by air transport.This makes the distribution region of charge become wider and wider,and the process should bring about the relationship between updraft and CROP.

    In the case of negative charge in the thunderstorm,the total amount of negative charge reaches a peak of–1.3×1012nC at the 21st minute,then begins to decrease,and the total amount keeps below 5 ×1011nC near the end of lightning activity(Fig.4a). The peak CROP for negative charge appears at the 16th minute and reaches about 5.5 C s?1.

    As we know,the total amount of charge is decided by the volume and density.Here,0.1 and 1 nC m?3are chosen as criteria of charge density to discuss the relationship of charge density with TAOC and the volume.The comparisons between TAOC for negative charge and the volumes with the two kinds of negative charge density criteria show that the trend of TAOC is obviously consistent with the trend of the volume with high charge density(Fig.4):they reach a peak at the same time and correlate well(correlation coefficient: 0.93).The volume with low charge density keeps increasing,which is different to the change of TAOC.This suggests that CROP should have a closer relationship with the change of the volume with high charge density.

    Fig.4.Evolution of the(a)total mass of negative charge on hydrometeors and(b)volumes with negative charge density on hydrometeors greater than 0.1 and 1 nC m?3.

    Meanwhile,the spatial relationship between updraft and the distribution of hydrometeors with negative surface charge(Fig.5)shows that the first region with negative charge appears in the upper part of the updraft.The distribution region of negative charge firstly extends in the vertical direction,and then spreads horizontally along the outflow near the top of the updraft.When the updraft is strong enough (Figs.5a–d),the negative charge distribution region also changes evidently;when the evident updraft channel disappears(Figs.5e–j),the change in the negative charge distribution region is inconspicuous,and the main change is expressed via a change in shape and height.

    A comparison between the volume with negative charge on hydrometeors and the volumes with updraft speeds greater than 1 and 5 m s?1shows that three curves have similar change if time is not taken into account(Fig.6).The volumes of both updrafts reach their peak earlier than the volume with negative charge on hydrometeors.This is reasonable because when the updraft begins to decrease,the region of hydrometeors with negative charge will keep expanding under the influence of the weakening updraft.Consequently,the peak time for the volume with negative charge on hydrometeors will arrive later than for the updraft volume.

    These results suggest that CROP should have a close relationship with the change of updraft volume.Although they might not be synchronous,their changes are similar.This makes updraft have a positive impact on CROP.According to Eq.(5),the positive influence of updraft on CROP,(Qt+?t?Qt),will be offset by that on CR.This causes updraft to inhibit the strengthening of RNCL through its positive impact on CROP,although this inhibition should,most of the time,be weaker than the strengthening effect of updraft on RNCL through its impact on CR.

    4.3ORCP

    Besides the neutralization of lightning,the charge in a thunderstorm could be reduced through precipitation(referred to as cold cloud precipitation),discharge in the upper atmosphere,and so on.Most of these phenomena have been studied,but quantitative observation is rare,and theoretical studies are primarily in the exploration stage.Therefore,in this study,only

    the reduction of charge through precipitation is analyzed.

    Fig.5.Profiles of negative charge on hydrometeors and updrafts at the(a)14th,(b)16th,(c)18th,(d)20th,(e)22nd, (f)24th,(g)26th,(h)28th,(i)30th,and(j)32nd minute.Lines indicate the contours of absolute charge density greater than 0.1 nC m?3.Vectors indicate the wind.

    Fig.6.Evolutions of the updraft volumes for w>1 and w>5 m s?1,and the updraft volume with negative charge on hydrometeors.

    In the case of negative charge in the thunderstorm,the magnitude of CR is the same as the magnitude of ORCP,and the maximum of ORCP is even slightly higher than that of CR(Fig.7).This indicates that the value of ORCP is significant.ORCP will severely weaken the promoting effect of CR on RNCL,and could even lead to a stop in lightning activity if they occur synchronously.However,in the thunderstorm simulated here,they are not synchronous. The main precipitation occurs after when the lightning activity finishes.During the lightning activity, CR has a significant advantage compared to ORCP.

    Fig.7.Evolutions of the charging rate by charging mechanisms in clouds and the charge fallout rate via precipitation.

    The observation of the first thunderstorm shows that the precipitation was weak over a very long period during this thunderstorm.The automated rain gauge in Beijing recorded precipitation of several millimeters only during about 2 h after the first lightning was detected.After the lightning activity had lasted for nearly 2 h,the recorded precipitation began to increase evidently,and the maximum reached 168 mm h?1.The tendency of the most evident precipitation lagging behind the main lightning activity is consistent between the observation and simulation.Indeed,such non-synchronization between precipitation and lightning activity has also been found in many other observations(e.g.,Zheng et al.,2010).Liu and Zipser (2008)analyzed the diurnal variation of precipitation and lightning activity in different regions of the tropics using satellite data.They found that the peak time of lightning activity was asynchronous with the peak time of precipitation in four of five regions,based on simultaneous precipitation and lightning data.These observations suggest that,on the whole,the time with frequent lightning flashes may not overlap with the time of heavy precipitation.Therefore,the simulation result with respect to the temporal relationship between lightning activity and the main precipitation should be typical for most thunderstorms. Consequently,ORCP(Qout)could be negligible compared with CR(Qchar)during most of the lightning activity period.

    From another point of view,it is conceivable that the non-synchronization of frequent lightning activity and heavy precipitation appears simply because the heavy ORCP seriously impairs the lightning activity.

    4.4 Discussion of the overall effect

    According to the above analysis,Eq.(5)could be simplified as follows for most of the lightning activity period:

    The polynomial fitting of the updraft volume with CR(Qchar)and CROP(Qt+?t?Qt)is given separately in Fig.8.Note that in order to express the different contributions of CR and CROP intuitively, the negative value of CROP for negative charge shows the increasing of TAOC.As can be seen,they both

    correlate well with the updraft volume;the correlation coefficients are 0.96 and 0.90,respectively.The trends of the fitting curves indicate that:(1)CR and CROP both have accelerated increases with increasing updraft volume;(2)CR keeps increasing as the updraft volume grows;(3)CROP increases with the growth of updraft volume only when the updraft volume is greater than 133 km3,and the value of CROP remains positive when the updraft volume is below 280 km3.This is because the volume and the mass of negative charge keep decreasing,although there is still a small region with weak updraft,while the updraft weakens considerably.Therefore,CROP is positive. However,when the updraft volume is greater than the threshold,CROP also increasingly accelerates.

    A comparison between both of the accelerations shows that CR has higher acceleration than CROP with the growth in updraft volume.According to the fitting equations(Fig.8)and Eq.(6),RNCL(y1+y2), synthesizing the results of CR and CROP,also illustrates an increasing acceleration with the growth in updraft volume.The inhibition of updraft on lightning activity becomes relatively evident only when the updraft is sufficiently strong,e.g.,the updraft volume is greater than 280 km3.As the updraft volume becomes larger in Fig.8,so does the difference between y1and y1+y2.During the final period of lightning activity,CROP becomes positive,which reduces the decreasing speed of RNCL.

    5.Conclusions

    To improve understanding of the factors impacting upon lightning activity,and according to the principle of the balance between the source and remit of one kind of polarity charge on the whole-thunderstorm scale,this study proposes that there should be three factors affecting RNCL,which represents the intensity of lightning activity,from the macro perspective of the thunderstorm:CR,CROP,and ORCP.Based on a simulation by using a 3D cloud resolution model coupled with charging and discharging processes,the impact of updraft on CR and CROP was analyzed. The results show that:

    1)Updraft has a positive influence on CR.The size of the region with the charging process is the key factor for CR.The main charging process always occurs within and near the updraft region.This causes a positive correlation between the updraft volume and the volume with the charging process.The correlation

    coefficient can reach approximately 0.96.Strengthening of the updraft promotes an increasing of CR through this correlation,and further,an increasing of RNCL.

    Fig.8.Combined effect of the charging rate and changing rate of total negative charge in the thunderstorm.

    2)Updraft has a positive influence on CROP.The delivery and divergence of charge in both the vertical and horizontal directions caused by the updraft lead to increasing TAOC,while the charging process keeps producing new charge.Therefore,strengthening of the updraft also promotes an increasing of TAOC over most of the time.This results in a positive correlation between the updraft volume and CROP,with a correlation coefficient of about 0.9.However,this positive influence of updraft on CROP will inhibit the increasing of RNCL.

    Additionally,the magnitude of ORCP is comparable to the magnitude of CR.The peak of ORCP(23 C s?1)in the simulated thunderstorm is even slightly greater than the peak of CR(20 C s?1).Because the main precipitation occurs after lightning activity, the inhibitory effect of ORCP on RNCL is negligible here.According to above analysis,strengthening of the updraft can not only promote lightning activity, but also inhibit it.Indeed,updraft has opposite influences on lightning activity simultaneously.On the whole,however,promotion is the dominant effect.

    A more comprehensive understanding of updraft on lightning activity can not only promote understanding of the relationship between dynamics and electrical activity in thunderstorms,but also provide a theoretical basis for explaining the relationship between lightning activity and other observed weather phenomena.It also supports further research on the construction of a logical and reliable lightning warning method,as well as a warning method for other severe weather using lightning data.

    Baker,M.B.,A.M.Blyth,H.J.Christian,et al.,1999: Relationships between lightning activity and various thundercloud parameters:Satellite and modelling studies.Atmos.Res.,51,221–236.

    Barthe,C.,and M.C.Barth,2008:Evaluation of a new lightning-produced NOx parameterization for cloud resolving models and its associated uncertainties. Atmos.Chem.Phys.,8,4691–4710.

    Barthe,C.,W.Deierling,and M.C.Barth,2010:Estimation of total lightning from various storm parameters:A cloud-resolving model study.J.Geophys. Res.,115,doi:10.1029/2010JD014405.

    Black,R.A.,and J.Hallett,1999:Electrification of the hurricane.J.Atmos.Sci.,56,2004–2028.

    Boccippio,D.J.,2002:Lightning scaling relations revisited.J.Atmos.Sci.,59,1086–1104.

    Bruning,E.C.,and D.R.MacGorman,2013:Theory and observations of controls on lightning flash size spectra.J.Atmos.Sci.,70,4012–4029.

    Carey,L.D.,and S.A.Rutledge,1996:A multiparameter radar case study of the microphysical and kinematic evolution of a lightning producing storm.Meteor. Atmos.Phys.,59,33–64.

    Carey,L.D.,M.J.Murphy,T.L.McCormick,et al., 2005:Lightning location relative to storm structure in a leading-line trailing-stratiform mesoscale convective system.J.Gophys.Res.,110,D03105,doi: 10.1029/2003JD004371.

    Deierling,W.,and W.A.Petersen,2008:Total lightning activity as an indicator of updraft characteristics. J.Geophys. Res.,113,D16210,doi: 10.1029/2007JD009598.

    Dye,J.E.,J.J.Jones,A.J.Weinheimer,et al.,1988: Observations within two regions of charge during initial thunderstorm electrification.Quart.J.Roy. Meteor.Soc.,114,1271–1290.

    Fuchs,B.,S.Rutledge,T.Lang,et al.,2014:Factors affecting lightning behavior in various regions of the United States.5th International Lightning Meteorological Conference.20–21 March,Tucson,Arizona, USA.

    Gardiner,B.,D.Lamb,R.L.Pitter,et al.,1985:Measurements of initial potential gradient and particle charges in a Montana summer thunderstorm.J. Geophys.Res.,90,6079–6086.

    Hallett,J.,and C.P.R.Saunders,1979:Charge separation associated with secondary ice crystal production.J.Atmos.Sci.,36,2230–2235.

    Heymsfield,G.M.,L.Tian,A.J.Heymsfield,et al., 2010:Characteristics of deep tropical and subtropical convection from nadir-viewing high-altitude airborne Doppler radar.J.Atmos.Sci.,67,285–308.

    Hu Zhijin and He Guanfang,1988:Numerical simulation of microphysical processes in cumulonimbus—Part I:Microphysical model.Acta Meteor.Sinica,2, 471–489.(in Chinese)

    Jacobson,E.A.,and E.P.Krider,1976:Electrostatic field changes produced by Florida lightning.J.Atmos.Sci.,33,103–117.

    Kasemir,H.W.,1960:A contribution to the electrostatic theory of a lightning discharge.J.Geophys.Res., 65,1873–1878.

    Kasemir,H.W.,2012:Theoretical and Experimental Determination of Field,Charge and Current on an Aircraft hit by Natural or Triggered Lightning. American Geophysical Union Washington,D.C., 885–693.

    Kuhlman,K.M.,D.R.MacGorman,M.I.Biggerstaff, et al.,2009:Lightning initiation in the anvils of two supercell storms.Geophys.Res.Lett.,36,L07802, doi:10.1029/2008GL036650.

    Lang,T.J.,and S.A.Rutledge,2002:Relationships between convective storm kinematics,precipitation, and lightning.Mon.Wea.Rev.,130,2492–2506.

    Liu,C.T.,and E.J.Zipser,2008:Diurnal cycles of precipitation,clouds,and lightning in the tropics from 9 years of TRMM observations.Geophys.Res. Lett.,35,L04819,doi:10.1029/2007GL032437.

    Marshall,T.C.,M.P.McCarthy,and W.D.Rust,1995: Electric field magnitudes and lightning initiation in thunderstorms.J.Geophys.Res.,100,7097–7103.

    Mansell,E.R.,D.R.MacGorman,C.L.Ziegler,et al.,2002:Simulated three-dimensional branched lightning in a numericalthunderstorm model. J.Geophy. Res.,107,ACL2-1-ACL2-12,doi: 10.1029/2000JD000244.

    McCaul,E.W.Jr.,S.J.Goodman,K.M.LaCasse,et al.,2009:Forecasting lightning threat using cloudresolving model simulations.Wea.Forecasting,24, 709–729.

    Mecikalski,R.M.,A.L.Bain,and L.D.Carey,2015: Radar and lightning observations of deep moist convection across northern Alabama during DC3:21 May 2012.Mon.Wea.Rev.,143,2774–2794.

    Miller,K.,A.Gadian,C.Saunders,et al.,2001:Modelling and observations of thundercloud electrification and lightning.Atmos.Res.,58,89–115.

    Palucki,J.L.,M.I.Biggerstaff,D.R.MacGorman, et al.,2011:Comparison between low-flash and non-lightning-producing convective areas within a mature mesoscale convective system.Wea.Forecasting,26,468–486.

    Pereyra,R.G.,E.E.Avila,N.E.Catellano,et al.,2000: A laboratory study of graupel charging.J.Geophys. Res.,105,20803–20812.

    Petersen,W.A.,S.A.Rutledge,and R.W.Orville, 1996: Cloud-to-ground lightning observations to TOGA COARE:Selected results and lightning location algorithm.Mon.Wea.Rev.,124,602–620.

    Petersen,W.A.,S.A.Rutledge,R.C.Cifelli,et al., 1999:Shipborne Dual-Doppler operations during TOGA COARE:Integrated observations of storm kinematics and electrification.Bull.Amer.Meteor. Soc.,80,81–97.

    Pickering,K.E.,Y.S.Wang,W.-K.Tao,et al.,1998: Vertical distributions of lightning NOx for use in regional and global chemical transport models.J. Geophys.Res.,103,31203–31216.

    Price,C.,and D.Rind,1992:A simple lightning parameterization for calculating global lightning distributions.J.Geophys.Res.,97,9919–9933.

    Reinhart,B.,H.Fuelberg,R.Blakeslee,et al.,2014: Understanding the relationships between lightning, cloud microphysics,and airborneradar-derived storm structure during Hurricane Karl(2010).Mon. Wea.Rev.,142,590–605.

    Rutledge,S.A.,E.R.Williams,and T.D.Keenan, 1992:The down upper Doppler and electricity experiment(DUNDEE):Overview and preliminary results.Bull.Amer.Meteor.Soc.,73,3–16.

    Shackford,C.R.,1960: Radarindications ofa precipitation-lightning relationship in New England thunderstorms.J.Meteor.,17,15–19.

    Takahashi,T.,1978:Riming electrification as a charge generation mechanism in thunderstorms.J.Atmos. Sci.,35,1536–1548.

    Takahashi,T.,1984:Thunderstorm electrification—A numerical study.J.Atmos.Sci.,41,2541–2558.

    Tan Yongbo,Tao Shanchang,Zhu Baoyou,et al.,2006: Thesimulation ofdouble-layerand branchingstructure of lightning in thunderstorm.Sci.China (Ser.D.),36,486–496.(in Chinese)

    Tan Yongbo,Tao Shanchang,Zhu Baoyou,et al.,2007: A simulation of the effects of intra-cloud lightning discharges on the charges and electrostatic potential distributions in a thundercloud.Chinese J.Geophys.,50,1053–1065.(in Chinese)

    Tessendorf,S.A.,L.J.Miller,K.C.Wiens,et al., 2005:The 29 June 2000 supercell observed during STEPS.Part I:Kinematics and microphysics.J. Atmos.Sci.,62,4127–4150.

    Wang Fang,Xiao Wenan,Lei Hengchi,et al.,2009:Numerical simulation of electricity characteristic of a thunderstorm case in summer of Jilin.Plateau Meteor.,28,385–394.(in Chinese)

    Wang Fei,Zhang Yijun,Zheng Dong,et al.,2015:Impact of the vertical velocity field on charging processes and charge separation in a simulated thunderstorm. J.Meteor.Res.,29,328–343.

    Weiss,S.A.,D.R.MacGorman,and K.M.Calhoun, 2012:Lightning in the anvils of supercell thunderstorms.Mon.Wea.Rev.,140,2064–2079.

    Wiens,K.C.,S.A.Rutledge,and S.A.Tessendorf,2005: The 29 June 2000 supercell observed during STEPS. Part II:Lightning and charge structure.J.Atmos. Sci.,62,4151–4177.

    Williams,E.R.,and R.M.Lhermitte,1983:Radar tests of the precipitation hypothesis for thunderstorm electrification.J.Geophys.Res.,88,10984–10992.

    Williams,E.R.,1985:Large-scale charge separation in thunderclouds.J.Geophys.Res.,90,6013–6025.

    Williams,E.R.,and S.Stanfill,2002:The physical origin of the land-ocean contrast in lightning activity. Comp.Rend.Phys.,3,1277–1292.

    Williams,E.R.,V.Mushtak,D.Rosenfeld,et al.,2005: Thermodynamic conditions favorable to superlative thunderstorm updrafts,mixed phase microphysics and lightning flash rate.Atmos.Res.,76,288–306.

    Workman,E.J.,and S.E.Reynolds,1949:Electrical activity as related to thunderstorm cell growth.Bull. Amer.Meteor.Soc.,30,142–149.

    Yao Wen,Zhang Yijun,Meng Qing,et al.,2013:A comparison of the characteristics of total and cloud-toground lightning activities in hailstorms.J.Meteor. Res.,27,282–293,doi:10.1007/s13351-013-0212-x.

    Zhang Wenjuan,Zhang Yijun,and Zhou Xiuji,2013: Lightning activity and precipitation characteristics of Typhoon Molave(2009)around its landfall.J. Meteor.Res.,27,742–757,doi:10.1007/s13351-013-0510-3.

    Zhang Yijun,Yan Muhong,and Liu Xinsheng,1999: Simulation study of discharge processes in thunderstorm.Chin.Sci.Bull.,44,2098–2102.(in Chinese)

    Zhang Yijun,P.R.Krehbiel,Liu Xinsheng,et al.,2003: Three dimensions structure of lightning discharge channel.Plateau Meteor.,22,217–220.(in Chinese)

    Zheng Dong,Zhang Yijun,Meng Qing,et al.,2010:Relationship between lightning activities and surface precipitation in thunderstorm weather in Beijing.J. Appl.Meteor.Sci.,21,287–297.

    Ziegler,C.L.,D.R.MacGorman,J.E.Dye,et al., 1991:A model evaluation of noninductive graupelice charging in the early electrification of a mountain thunderstorm.J.Geophys.Res.,96,12833–12855.

    Ziegler,C.L.,and D.R.MacGorman,1994:Observed lightning morphology relative to modeled space charge and electric field distributions in a tornadic storm.J.Atmos.Sci.,51,833–851.

    :Wang Fei,Zhang Yijun,and Zheng Dong,2015:Impact of updraft on neutralized charge rate by lightning in thunderstorms:A simulation case study.J.Meteor.Res.,29(6),997–1010,

    10.1007/s13351-015-5023-9.

    Supported by the National Natural Science Foundation of China(41205001),National(Key)Basic Research and Development (973)Program of China(2014CB441406),National Major Scientific Instruments and Equipment Development Special(61327810), and Basic Research Funds of CAMS(2012Y005 and 2013Z006).

    ?Corresponding author:feiwang@camscma.cn.

    ?The Chinese Meteorological Society and Springer-Verlag Berlin Heidelberg 2015

    (Received March 11,2015;in final form August 27,2015)

    猜你喜歡
    義軍王飛
    滬指失守3000點
    家人春節(jié)照
    當代作家(2023年3期)2023-04-23 01:58:26
    坐公交車鬧出的笑話
    A novel similarity measure for mining missing links in long-path networks
    以微粒觀構建思維模型用于突破高三復習離子反應的難點
    市場觀察(1)
    市場觀察(2)
    A Similarity Transformation and the Decay Mode Solutions for Three-Dimensional Cylindrical Kadomtsev-Petviashvili Equation
    應用數學(2019年2期)2019-03-30 08:21:34
    市場觀察(1)
    新出唐代張淮澄墓志所見歸義軍史事考
    敦煌學輯刊(2017年1期)2017-11-10 02:32:16
    乱人视频在线观看| 欧美中文综合在线视频| 国产亚洲精品久久久com| 亚洲精品久久国产高清桃花| 亚洲性夜色夜夜综合| 一二三四社区在线视频社区8| 五月伊人婷婷丁香| 欧美又色又爽又黄视频| 日韩中文字幕欧美一区二区| 国产亚洲精品一区二区www| av专区在线播放| 亚洲精品456在线播放app | 真实男女啪啪啪动态图| 99热6这里只有精品| 在线观看免费午夜福利视频| 黄片大片在线免费观看| 免费搜索国产男女视频| 69人妻影院| 91久久精品国产一区二区成人 | 好男人在线观看高清免费视频| 午夜免费男女啪啪视频观看 | 在线观看日韩欧美| 一级作爱视频免费观看| a在线观看视频网站| 亚洲av第一区精品v没综合| 亚洲av一区综合| 俄罗斯特黄特色一大片| 在线观看一区二区三区| 又粗又爽又猛毛片免费看| 亚洲精品亚洲一区二区| 最好的美女福利视频网| 99国产极品粉嫩在线观看| 波多野结衣高清无吗| 日韩人妻高清精品专区| 日韩人妻高清精品专区| 亚洲真实伦在线观看| 国产真实伦视频高清在线观看 | 别揉我奶头~嗯~啊~动态视频| 精品一区二区三区视频在线 | 欧美最新免费一区二区三区 | 又爽又黄无遮挡网站| 国产麻豆成人av免费视频| 非洲黑人性xxxx精品又粗又长| 在线观看日韩欧美| 最新在线观看一区二区三区| 国产成人av教育| 国产精华一区二区三区| 少妇裸体淫交视频免费看高清| av在线天堂中文字幕| 超碰av人人做人人爽久久 | 国产成人aa在线观看| 亚洲五月婷婷丁香| 日本一二三区视频观看| 久久九九热精品免费| 亚洲av一区综合| а√天堂www在线а√下载| 性色av乱码一区二区三区2| 听说在线观看完整版免费高清| 亚洲性夜色夜夜综合| 久久久久久久精品吃奶| 欧美黑人巨大hd| 亚洲精品久久国产高清桃花| 国产精华一区二区三区| 青草久久国产| 免费人成视频x8x8入口观看| 日韩av在线大香蕉| 夜夜爽天天搞| 久久久久九九精品影院| 欧美乱妇无乱码| 欧美激情在线99| 性色av乱码一区二区三区2| 免费大片18禁| 免费大片18禁| 性色av乱码一区二区三区2| 男女之事视频高清在线观看| 最近视频中文字幕2019在线8| 国产精品女同一区二区软件 | 九色国产91popny在线| av黄色大香蕉| 免费人成在线观看视频色| a级一级毛片免费在线观看| 精品久久久久久久人妻蜜臀av| 人人妻人人看人人澡| 亚洲欧美日韩东京热| 国产精品一区二区三区四区免费观看 | 99久久精品热视频| 久久久成人免费电影| 午夜a级毛片| 精品乱码久久久久久99久播| 国产精品一区二区三区四区久久| 哪里可以看免费的av片| 成年女人永久免费观看视频| 九色成人免费人妻av| 九九热线精品视视频播放| 一个人免费在线观看电影| 无遮挡黄片免费观看| 久久国产乱子伦精品免费另类| 精品人妻一区二区三区麻豆 | 桃红色精品国产亚洲av| 一级毛片女人18水好多| 国产69精品久久久久777片| 亚洲人成伊人成综合网2020| 午夜福利成人在线免费观看| 日本一本二区三区精品| 国产亚洲精品久久久久久毛片| 欧美不卡视频在线免费观看| 成人特级黄色片久久久久久久| 两个人的视频大全免费| 欧美又色又爽又黄视频| 热99在线观看视频| 床上黄色一级片| 欧美性感艳星| 中国美女看黄片| 国产精品久久视频播放| 国产欧美日韩一区二区精品| 无人区码免费观看不卡| 成年免费大片在线观看| 成人18禁在线播放| 人妻丰满熟妇av一区二区三区| 色综合婷婷激情| 嫩草影院入口| 桃红色精品国产亚洲av| 国产麻豆成人av免费视频| 欧美成人免费av一区二区三区| 久久伊人香网站| 国内精品久久久久精免费| 国产伦在线观看视频一区| 欧美日韩精品网址| 日韩欧美 国产精品| 成人三级黄色视频| 91麻豆精品激情在线观看国产| 综合色av麻豆| 欧美xxxx黑人xx丫x性爽| 在线观看免费午夜福利视频| 首页视频小说图片口味搜索| 久久久久精品国产欧美久久久| 国产色婷婷99| 久久精品亚洲精品国产色婷小说| 国产精品一区二区三区四区免费观看 | 成人亚洲精品av一区二区| 高清在线国产一区| 亚洲国产欧美网| 亚洲色图av天堂| 国内精品美女久久久久久| 久久精品人妻少妇| 欧美av亚洲av综合av国产av| 97超视频在线观看视频| 搡老妇女老女人老熟妇| 别揉我奶头~嗯~啊~动态视频| 熟女少妇亚洲综合色aaa.| 国产成人av激情在线播放| 9191精品国产免费久久| 一个人看的www免费观看视频| 国产一级毛片七仙女欲春2| 日本免费一区二区三区高清不卡| 99久久精品热视频| 久久久久久大精品| 俺也久久电影网| 久久久久性生活片| 国产成人av激情在线播放| 免费看光身美女| 老汉色∧v一级毛片| 九九久久精品国产亚洲av麻豆| 床上黄色一级片| 亚洲国产精品成人综合色| 成人性生交大片免费视频hd| 色视频www国产| 中文字幕人妻熟人妻熟丝袜美 | 色吧在线观看| 亚洲av一区综合| 精品欧美国产一区二区三| 99热这里只有精品一区| 国产探花在线观看一区二区| 无限看片的www在线观看| 一级黄色大片毛片| 日本五十路高清| 国产私拍福利视频在线观看| 国产伦人伦偷精品视频| 午夜精品在线福利| 99久久综合精品五月天人人| 麻豆一二三区av精品| 久久香蕉精品热| 热99在线观看视频| 最近视频中文字幕2019在线8| 黄片大片在线免费观看| 成年女人毛片免费观看观看9| 欧美大码av| 特级一级黄色大片| 免费在线观看影片大全网站| av欧美777| 久久久久久大精品| 久久久久九九精品影院| 成熟少妇高潮喷水视频| 亚洲无线观看免费| 两个人的视频大全免费| 国产精品免费一区二区三区在线| 欧美性感艳星| 最近在线观看免费完整版| 老司机午夜福利在线观看视频| 精品久久久久久久久久免费视频| 一本综合久久免费| 国产一区二区三区视频了| 国产精品香港三级国产av潘金莲| 亚洲无线在线观看| 啪啪无遮挡十八禁网站| 可以在线观看的亚洲视频| 中文在线观看免费www的网站| 精品一区二区三区av网在线观看| 欧美精品啪啪一区二区三区| 又粗又爽又猛毛片免费看| 一进一出好大好爽视频| 香蕉丝袜av| 熟妇人妻久久中文字幕3abv| 91在线观看av| 国产精品1区2区在线观看.| 好男人电影高清在线观看| 高清日韩中文字幕在线| 久久精品影院6| 国产精品 国内视频| www国产在线视频色| 欧美一级a爱片免费观看看| 在线观看66精品国产| 午夜影院日韩av| 五月玫瑰六月丁香| 精品久久久久久久末码| 国产视频内射| 久久久久久久亚洲中文字幕 | 免费看光身美女| 91字幕亚洲| 男女床上黄色一级片免费看| 天天躁日日操中文字幕| 色视频www国产| av专区在线播放| 亚洲中文字幕一区二区三区有码在线看| 国产精品影院久久| 深夜精品福利| 一边摸一边抽搐一进一小说| 亚洲精品456在线播放app | 国产午夜福利久久久久久| 亚洲成人中文字幕在线播放| 有码 亚洲区| 毛片女人毛片| www.熟女人妻精品国产| 国产精品亚洲美女久久久| 国产三级在线视频| 国产精品三级大全| 国产精华一区二区三区| 在线播放国产精品三级| 又黄又爽又免费观看的视频| 亚洲国产精品sss在线观看| 国产精品自产拍在线观看55亚洲| 亚洲精品久久国产高清桃花| 级片在线观看| 一级黄色大片毛片| 婷婷六月久久综合丁香| 女人十人毛片免费观看3o分钟| 国产日本99.免费观看| 99热精品在线国产| 不卡一级毛片| 美女黄网站色视频| 国产精品 国内视频| 精品99又大又爽又粗少妇毛片 | 在线观看美女被高潮喷水网站 | 香蕉丝袜av| 首页视频小说图片口味搜索| 中文资源天堂在线| 国产av在哪里看| 丰满人妻一区二区三区视频av | www.熟女人妻精品国产| 国语自产精品视频在线第100页| 五月伊人婷婷丁香| 国内毛片毛片毛片毛片毛片| 久99久视频精品免费| 国产激情欧美一区二区| 身体一侧抽搐| 长腿黑丝高跟| 日韩有码中文字幕| 天堂av国产一区二区熟女人妻| 国产精品综合久久久久久久免费| 黄片小视频在线播放| 精品久久久久久久末码| 内地一区二区视频在线| 免费无遮挡裸体视频| 女人高潮潮喷娇喘18禁视频| 黄色视频,在线免费观看| 亚洲成a人片在线一区二区| 亚洲成av人片在线播放无| 久久欧美精品欧美久久欧美| 国产亚洲欧美在线一区二区| 免费观看的影片在线观看| 精品久久久久久久久久久久久| 国产精品久久久久久久电影 | 国产中年淑女户外野战色| 99久久99久久久精品蜜桃| 两个人的视频大全免费| 九色成人免费人妻av| 欧美日韩一级在线毛片| 国产成人欧美在线观看| 小蜜桃在线观看免费完整版高清| 12—13女人毛片做爰片一| 亚洲精品色激情综合| 国产精品嫩草影院av在线观看 | 日韩中文字幕欧美一区二区| 一a级毛片在线观看| 91在线精品国自产拍蜜月 | 色在线成人网| 叶爱在线成人免费视频播放| 很黄的视频免费| 国产又黄又爽又无遮挡在线| 亚洲电影在线观看av| 免费高清视频大片| 成人性生交大片免费视频hd| 99热精品在线国产| 在线观看av片永久免费下载| 麻豆国产97在线/欧美| 黄色片一级片一级黄色片| 免费在线观看日本一区| 99视频精品全部免费 在线| 色综合亚洲欧美另类图片| 久久6这里有精品| 免费看日本二区| 欧美不卡视频在线免费观看| 草草在线视频免费看| 高清毛片免费观看视频网站| 欧美一级a爱片免费观看看| 亚洲专区中文字幕在线| 黄色视频,在线免费观看| 成年免费大片在线观看| 在线十欧美十亚洲十日本专区| 久久精品国产亚洲av香蕉五月| 久久精品91无色码中文字幕| 观看免费一级毛片| 午夜福利在线观看免费完整高清在 | 亚洲黑人精品在线| 久久伊人香网站| 久久精品91蜜桃| 亚洲av成人av| 别揉我奶头~嗯~啊~动态视频| 看片在线看免费视频| 男女午夜视频在线观看| 午夜亚洲福利在线播放| 国产精品亚洲av一区麻豆| 免费看光身美女| 日韩欧美一区二区三区在线观看| 国产精品 国内视频| 国产成人啪精品午夜网站| 在线观看66精品国产| 色综合亚洲欧美另类图片| 桃红色精品国产亚洲av| 亚洲天堂国产精品一区在线| 午夜两性在线视频| 十八禁网站免费在线| 国产成人影院久久av| 夜夜躁狠狠躁天天躁| 国产免费一级a男人的天堂| 一二三四社区在线视频社区8| 日本黄色片子视频| 757午夜福利合集在线观看| 国产精品一及| 国产成人av教育| 国产伦一二天堂av在线观看| 久久精品国产自在天天线| 免费搜索国产男女视频| 黄色丝袜av网址大全| 999久久久精品免费观看国产| 日本 av在线| 午夜老司机福利剧场| 高潮久久久久久久久久久不卡| 国产爱豆传媒在线观看| 亚洲av一区综合| 99久久成人亚洲精品观看| 欧美三级亚洲精品| 少妇的逼好多水| 国产精品亚洲av一区麻豆| av天堂中文字幕网| 久久精品影院6| 美女高潮的动态| 九九在线视频观看精品| 亚洲熟妇中文字幕五十中出| 非洲黑人性xxxx精品又粗又长| 两性午夜刺激爽爽歪歪视频在线观看| 午夜a级毛片| 亚洲国产色片| 欧美日韩国产亚洲二区| 丰满的人妻完整版| 国产美女午夜福利| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品av视频在线免费观看| 亚洲最大成人手机在线| 亚洲中文日韩欧美视频| 亚洲精品成人久久久久久| 日韩高清综合在线| 中文字幕高清在线视频| 欧美bdsm另类| 男插女下体视频免费在线播放| 看免费av毛片| 美女免费视频网站| 真人一进一出gif抽搐免费| 欧美乱码精品一区二区三区| 亚洲美女视频黄频| 午夜精品在线福利| 国产伦一二天堂av在线观看| 黄片小视频在线播放| 一级黄片播放器| 一区二区三区激情视频| 嫩草影院精品99| 日本在线视频免费播放| 久久欧美精品欧美久久欧美| 亚洲aⅴ乱码一区二区在线播放| 亚洲第一电影网av| 国产精品综合久久久久久久免费| 亚洲av第一区精品v没综合| 真实男女啪啪啪动态图| 午夜福利18| 美女cb高潮喷水在线观看| 国产免费男女视频| 香蕉av资源在线| 首页视频小说图片口味搜索| 黄色成人免费大全| 男女床上黄色一级片免费看| 在线看三级毛片| 亚洲人成网站在线播放欧美日韩| 色综合站精品国产| 国内久久婷婷六月综合欲色啪| 一区二区三区高清视频在线| 丝袜美腿在线中文| 久久人妻av系列| 国产成人系列免费观看| 性色avwww在线观看| 国产亚洲精品一区二区www| 色综合亚洲欧美另类图片| 一级黄片播放器| 国产精品98久久久久久宅男小说| 琪琪午夜伦伦电影理论片6080| 亚洲成a人片在线一区二区| 波野结衣二区三区在线 | 男女视频在线观看网站免费| 国产精品野战在线观看| 亚洲成人中文字幕在线播放| 一本久久中文字幕| 久久久国产成人精品二区| 国产色爽女视频免费观看| 2021天堂中文幕一二区在线观| 禁无遮挡网站| 两性午夜刺激爽爽歪歪视频在线观看| 99热只有精品国产| 狂野欧美激情性xxxx| 全区人妻精品视频| 男女之事视频高清在线观看| 国产成人av激情在线播放| 国产一区二区在线观看日韩 | 窝窝影院91人妻| 日韩欧美一区二区三区在线观看| 久久精品国产99精品国产亚洲性色| 午夜两性在线视频| 一夜夜www| 欧美一级a爱片免费观看看| 搡老岳熟女国产| 精品熟女少妇八av免费久了| 一个人看的www免费观看视频| 偷拍熟女少妇极品色| 亚洲欧美日韩高清在线视频| 97人妻精品一区二区三区麻豆| 亚洲人与动物交配视频| 国语自产精品视频在线第100页| 丰满乱子伦码专区| 日本三级黄在线观看| 亚洲精品一区av在线观看| 一级作爱视频免费观看| 99久久九九国产精品国产免费| 久久精品综合一区二区三区| 日韩大尺度精品在线看网址| 午夜a级毛片| 国内精品久久久久久久电影| 精品一区二区三区视频在线观看免费| 观看美女的网站| 窝窝影院91人妻| 日韩欧美一区二区三区在线观看| 久久精品国产99精品国产亚洲性色| 国产爱豆传媒在线观看| 丝袜美腿在线中文| 国产免费男女视频| 丰满乱子伦码专区| 精品久久久久久久久久免费视频| av欧美777| 老司机午夜福利在线观看视频| 悠悠久久av| 欧美xxxx黑人xx丫x性爽| 搞女人的毛片| 哪里可以看免费的av片| 色尼玛亚洲综合影院| 亚洲欧美日韩无卡精品| 欧美日韩一级在线毛片| 九九在线视频观看精品| 性色av乱码一区二区三区2| 国产精品亚洲美女久久久| 九色成人免费人妻av| 久久人妻av系列| 亚洲国产精品999在线| xxxwww97欧美| 国内少妇人妻偷人精品xxx网站| 一本久久中文字幕| 免费无遮挡裸体视频| 一夜夜www| 9191精品国产免费久久| 亚洲内射少妇av| 欧美大码av| 婷婷六月久久综合丁香| 亚洲欧美精品综合久久99| 国产野战对白在线观看| 欧美日韩一级在线毛片| 好看av亚洲va欧美ⅴa在| 亚洲精品粉嫩美女一区| 99久久成人亚洲精品观看| 91av网一区二区| 国内精品久久久久久久电影| 国产高清视频在线播放一区| 男人的好看免费观看在线视频| 日韩精品中文字幕看吧| 亚洲无线在线观看| 国产欧美日韩精品一区二区| 热99在线观看视频| 欧美中文日本在线观看视频| 精品人妻偷拍中文字幕| 亚洲精品影视一区二区三区av| 国产欧美日韩精品一区二区| 国产精品一及| 亚洲精品乱码久久久v下载方式 | 午夜老司机福利剧场| 久久久久国产精品人妻aⅴ院| 亚洲国产欧美人成| 一级毛片高清免费大全| 国产精品免费一区二区三区在线| 精品99又大又爽又粗少妇毛片 | 亚洲人成网站在线播| 欧美日本亚洲视频在线播放| 两人在一起打扑克的视频| 精品电影一区二区在线| 欧美成人一区二区免费高清观看| 国产毛片a区久久久久| 欧美性猛交╳xxx乱大交人| 看片在线看免费视频| 国产欧美日韩精品亚洲av| 国产 一区 欧美 日韩| avwww免费| 蜜桃亚洲精品一区二区三区| 一级a爱片免费观看的视频| 久久久国产成人免费| 精品一区二区三区视频在线观看免费| 国产精品久久久久久人妻精品电影| 欧美最黄视频在线播放免费| 午夜a级毛片| 国产亚洲av嫩草精品影院| 99热这里只有精品一区| 亚洲无线观看免费| 成年女人看的毛片在线观看| 午夜福利在线在线| 丰满人妻一区二区三区视频av | 久久久国产精品麻豆| 免费看日本二区| 日韩欧美 国产精品| 69av精品久久久久久| 国产欧美日韩精品一区二区| 长腿黑丝高跟| 日韩欧美国产在线观看| 国产蜜桃级精品一区二区三区| 精品99又大又爽又粗少妇毛片 | 首页视频小说图片口味搜索| 村上凉子中文字幕在线| 久久久国产成人免费| 女同久久另类99精品国产91| 天天躁日日操中文字幕| 国语自产精品视频在线第100页| 日韩av在线大香蕉| 少妇高潮的动态图| 久久亚洲真实| 国产精品乱码一区二三区的特点| 19禁男女啪啪无遮挡网站| 国产精品一区二区免费欧美| 丰满人妻一区二区三区视频av | 亚洲最大成人中文| 成人午夜高清在线视频| 日本熟妇午夜| 午夜亚洲福利在线播放| 免费看光身美女| 久久久精品欧美日韩精品| 亚洲人成网站高清观看| 久久精品影院6| 国产激情欧美一区二区| 精品乱码久久久久久99久播| 精品一区二区三区av网在线观看| 88av欧美| 亚洲欧美日韩高清在线视频| 长腿黑丝高跟| 高潮久久久久久久久久久不卡| 亚洲熟妇中文字幕五十中出| 天堂av国产一区二区熟女人妻| 在线免费观看不下载黄p国产 | 一本久久中文字幕| 国产极品精品免费视频能看的| 99在线视频只有这里精品首页| 日本黄色视频三级网站网址| 国产精品久久视频播放| 九九在线视频观看精品| 久久天躁狠狠躁夜夜2o2o| 欧美日韩瑟瑟在线播放| 最近视频中文字幕2019在线8| 中文资源天堂在线| 亚洲av一区综合| h日本视频在线播放| 18禁裸乳无遮挡免费网站照片| 国产视频一区二区在线看| 日韩欧美精品免费久久 | 久久精品91蜜桃| 午夜精品久久久久久毛片777|